search

EnMAP-Box / enmap-box

EnMAP-Box source code repository. See https://enmap-box.readthedocs.io for documentation
GNU General Public License v3.0
39 stars 16 forks source link

ENSOMAP - water mask visualisation #328

Closed chabrillat closed 1 year ago

chabrillat commented 1 year ago

Running ENSOMAP In EnMAP-Box 3.11 with QGIS 3.28.1 on my windows laptop, the following bug is appearing: it says that *_water_mask.dat contains only nan values

When we open the same file with ENVI, the file contains 0s and 1s. The water mask file is correct but somehow ENMAP-Box says it is 'nan' although visualisation seems correct too (the file is mostly 0, and only a few or one 1)

In attachment, i have put the water mask file and also the screenshot with the nan.

Please can you investigate this story, image

thank you, Sabine.

janzandr commented 1 year ago

Hi Sabine, I can't see the attached _water_mask.dat file. Please upload or send via email.

chabrillat commented 1 year ago

Hi Andreas,

here it is (i had an issue with the upload, then via email),

thanks

On 16/01/2023 15:40, Andreas Janz wrote:

Hi Sabine, I can't see the attached _water_mask.dat file. Please upload or send via email.

— Reply to this email directly, view it on GitHub https://github.com/EnMAP-Box/enmap-box/issues/328#issuecomment-1384162792, or unsubscribe https://github.com/notifications/unsubscribe-auth/A5I36GC5WIZ4KST6H3FN7ODWSVMXVANCNFSM6AAAAAAT4Z465M. You are receiving this because you authored the thread.Message ID: @.***>

ENVI description = File generated by Hysoma samples = 829 lines = 828 bands = 1 header offset = 0 file type = ENVI Standard data type = 2 interleave = bsq sensor type = Unknown byte order = 0 map info = {UTM, 1.000, 1.000, 580297.000, 4082398.000, 5.0000000000e+000, 5.0000000000e+000, 30, North, units=Meters} x start = 172 y start = 925 coordinate system string = {PROJCS["UTM_Zone_30N", GEOGCS["GCS_North_American_1927", DATUM["D_North_American_1927", SPHEROID["Clarke_1866", 6378206.4, 294.9786982]], PRIMEM["Greenwich", 0.0], UNIT["Degree", 0.0174532925199433]], PROJECTION["Transverse_Mercator"], PARAMETER["False_Easting", 500000.0], PARAMETER["False_Northing", 0.0], PARAMETER["Central_Meridian", -3.0], PARAMETER["Scale_Factor", 0.9996], PARAMETER["Latitude_Of_Origin", 0.0], UNIT["Meter", 1.0]]}

janzandr commented 1 year ago

Thanks @chabrillat for sharing the file. I checked and everything is fine. You just have to use the "Whole raster" and "Actual (slower)" settings, when calculating the data stretch: image

chabrillat commented 1 year ago

great, many thanks! i wonder what was plotted before?! i usually never had to change these options