arcsi CLOUDS product segfaulting for landsat images

[Dan Slayback]

I've been trying to generate the CLOUDS layer using arcsi.py, for Landsat7 imagery, and keep getting segmentation faults when it tries to add a histogram/color table.  Here's an example call:

arcsi.py -i LE70340322013188EDC00_MTL.txt -s ls7 -f GTiff --outbasename arcsi_cld -o . -p TOA CLOUDS -t /tmp

And the last part of the output:
...
Generate Cloud Mask
Opening: ./arcsi_cld_rad_toa.tif
Opening: ./arcsi_cld_thermal.tif
Opening: ./arcsi_cld_sat.tif
Apply first pass FMask to classifiy initial clear sky regions...
Started .0..10..20..30..40..50..60..70..80..90..100. Complete.
Adding colour table and histogram to image
Get Image Min and Max.
Started .0..10..20..30..40..50..60..70..80..90..100. Complete.
Get Image Histogram.
Started .0..10..20..30..40..50..60..70..80..90..100. Complete.
Adding Histogram and Colour Table to image file
Segmentation fault

Is this a known issue or should this be working? I've tried on several L7 and L8 scenes, and get the same crash at the same place. Also, same issue when calling rsgislib.imagecalibration.applyLandsatTMCloudFMask directly from an interactive python session (3.4.1, everything fully updated with conda).

thanks!

[Daniel Clewley]

Hi Dan,

I think the problem is with support for histograms and colour tables in Geotiffs, have you tried using ERDAS imagine or KEA format as an output? Most of the testing has been carried out using the KEA format (.kea).

Thanks,

Dan

--
You received this message because you are subscribed to the Google Groups "RSGISLib Support" group.
To unsubscribe from this group and stop receiving emails from it, send an email to rsgislib-suppo...@googlegroups.com.
To post to this group, send email to rsgislib...@googlegroups.com.
Visit this group at http://groups.google.com/group/rsgislib-support.
For more options, visit https://groups.google.com/d/optout.

[Dan Slayback]

Aha - yes, that did it! I switched to KEA output and no problems.

thanks very much for the help.

Dan