<div dir="ltr"><div><div><br><br>On Tue, Feb 6, 2018 at 10:25 AM, Veronica Andreo <<a href="mailto:veroandreo@gmail.com">veroandreo@gmail.com</a>> wrote:<br>><br>> Sorry, sorry... that was a different band... I realize now. So, for the record, no issues,<br>> gdalinfo -mm T21KZP_20171212T134159_B02.jp2<br>> reports the same range as the one I see after importing in grass.<br>><br>> Computed Min/Max=15.000,17547.000<br>><br>> Again, sorry for the noise and thanks for the -mm hint, MM :)<br><br></div>Good to hear that gdalinfo -mm and GRASS stats are indeed identical.<br><br></div>A hint: use <a href="http://r.info">r.info</a> -s instead of r.univar -g for simple stats (number of non-null cells, min, max, mean, stddev)  on the full raster map, it's faster.<br><div><div><br></div><div>The reason why gdalinfo with and without -mm might report different values is that some software creating raster datasets for some obscure reason estimates min, max, mean, stddev and for some even more obscure reason writes these estimates to metadata. At least regarding min, max one would expect actual min, max, not some estimate with estimated min > actual min and estimated max < actual max.<br></div><div><br></div><div>Markus M<br></div><div>><br>> best,<br>> Vero<br>><br>> 2018-02-05 20:54 GMT+01:00 Veronica Andreo <<a href="mailto:veroandreo@gmail.com">veroandreo@gmail.com</a>>:<br>>><br>>> Forgot the list, sorry...<br>>><br>>> -------<br>>><br>>> Hi Markus,<br>>><br>>> here the output of gdalinfo -mm for the same band reported earlier:<br>>><br>>> gdalinfo -mm T21KZP_20171212T134159_B01.jp2<br>>> Driver: JP2OpenJPEG/JPEG-2000 driver based on OpenJPEG library<br>>> Files: T21KZP_20171212T134159_B01.jp2<br>>> [...]<br>>><br>>> Band 1 Block=192x192 Type=UInt16, ColorInterp=Gray<br>>>   Min=908.000 Max=1511.000   Computed Min/Max=934.000,6934.000<br>>>   Minimum=908.000, Maximum=1511.000, Mean=1138.806, StdDev=74.933<br>>>   Overviews: 915x915, 457x457, 228x228, 114x114<br>>>   Overviews: arbitrary<br>>>   Metadata:<br>>>     STATISTICS_MAXIMUM=1511<br>>>     STATISTICS_MEAN=1138.8064<br>>>     STATISTICS_MINIMUM=908<br>>>     STATISTICS_STDDEV=74.933459275813<br>>>   Image Structure Metadata:<br>>>     COMPRESSION=JPEG2000<br>>>     NBITS=15<br>>><br>>> still differs a lot from what I get after importing into GRASS.<br>>><br>>> System Info                                                                    <br>>> GRASS version: 7.5.svn                                                         <br>>> GRASS SVN revision: r72203                                                     <br>>> Build date: 2018-01-01                                                         <br>>> Build platform: x86_64-pc-linux-gnu                                            <br>>> GDAL: 2.1.3                                                                    <br>>> PROJ.4: 4.9.3                                                                  <br>>> GEOS: 3.6.1                                                                    <br>>> SQLite: 3.20.1                                                                 <br>>> Python: 2.7.14                                                                 <br>>> wxPython: 3.0.2.0                                                              <br>>> Platform: Linux-4.14.14-200.fc26.x86_64-x86_64-with-fedora-26-Twenty_Six<br>>><br>>> best,<br>>> Vero<br>>><br>>><br>>> 2018-02-05 19:49 GMT+01:00 Markus Metz <<a href="mailto:markus.metz.giswork@gmail.com">markus.metz.giswork@gmail.com</a>>:<br>>>><br>>>><br>>>><br>>>> On Mon, Feb 5, 2018 at 1:14 PM, Veronica Andreo <<a href="mailto:veroandreo@gmail.com">veroandreo@gmail.com</a>> wrote:<br>>>> ><br>>>> > I tested also with 74 release branch... same result<br>>>> ><br>>>> > 2018-02-05 13:05 GMT+01:00 Veronica Andreo <<a href="mailto:veroandreo@gmail.com">veroandreo@gmail.com</a>>:<br>>>> >><br>>>> >> Hi devs<br>>>> >><br>>>> >> I'm testing Martin's new add-on, r.sentinel.import. I created an UTM21S location for my data and imported them with:<br>>>> >><br>>>> >> r.sentinel.import input=Downloads/S2B_MSIL1C_20171212T134159_N0206_R124_T21KZP_20171212T201017.SAFE/GRANULE/L1C_T21KZP_A004011_20171212T134158/IMG_DATA/<br>>>> >><br>>>> >> The problem arose when I tried to use i.color.enhance to display an RGB: I got blank screens, so I thought there must be something with the range of values.<br>>>> >><br>>>> >> r.univar map=T21KZP_20171212T134159_B02@PERMANENT                              <br>>>> >> total null and non-null cells: 120560400<br>>>> >> total null cells: 0<br>>>> >> Of the non-null cells:<br>>>> >> ----------------------<br>>>> >> n: 120560400<br>>>> >> minimum: 15<br>>>> >> maximum: 17547<br>>>> >> range: 17532<br>>>> >> mean: 913.778<br>>>> >> mean of absolute values: 913.778<br>>>> >> standard deviation: 134.322<br>>>> >> variance: 18042.4<br>>>> >> variation coefficient: 14.6996 %<br>>>> >> sum: 110165471238<br>>>> >><br>>>> >> that already seemed odd...<br>>>> >> I checked the original band2 jp2 file... Values are different...<br>>>> >><br>>>> >> gdalinfo T21KZP_20171212T134159_B02.jp2<br>>>> >> Driver: JP2OpenJPEG/JPEG-2000 driver based on OpenJPEG library<br>>>> >> Files: T21KZP_20171212T134159_B02.jp2<br>>>> >> ...<br>>>> >> Band 1 Block=1024x1024 Type=UInt16, ColorInterp=Gray<br>>>> >>   Min=602.000 Max=1576.000<br>>>> >>   Minimum=602.000, Maximum=1576.000, Mean=916.799, StdDev=123.509<br>>>> >>   Overviews: 5490x5490, 2745x2745, 1372x1372, 686x686<br>>>> >>   Overviews: arbitrary<br>>>> >>   Metadata:<br>>>> >>     STATISTICS_MAXIMUM=1576<br>>>> >>     STATISTICS_MEAN=916.7992<br>>>> >>     STATISTICS_MINIMUM=602<br>>>> >>     STATISTICS_STDDEV=123.50943477872<br>>>> >>   Image Structure Metadata:<br>>>> >>     COMPRESSION=JPEG2000<br>>>> >>     NBITS=15<br>>>><br>>>> Can you try again with gdalinfo -mm? This forces computation of the actual min/max values, which can differ from what is stored in the metadata (an approximation).<br>>>><br>>>> There can also be something wrong with i.color.enhance<br>>>><br>>>> Markus M<br>>>><br>>>> >><br>>>> >> Just to check, I imported with r.in.gdal and I'm getting exactly the same as with r.sentinel.import.<br>>>> >><br>>>> >> What am I doing wrong here? Why would r.in.gdal change values like that? What else needs to be set? I don't remember to have this problem before and it happens in all bands as far as I've checked.<br>>>> >>  <br>>>> >> Ah, I'm using trunk r72203<br>>>> >><br>>>> >> Thanks in advance for any help :)<br>>>> >><br>>>> >> Cheers,<br>>>> >> Vero<br>>>> >><br>>>> >><br>>>> >><br>>>> ><br>>>> ><br>>>> > _______________________________________________<br>>>> > grass-dev mailing list<br>>>> > <a href="mailto:grass-dev@lists.osgeo.org">grass-dev@lists.osgeo.org</a><br>>>> > <a href="https://lists.osgeo.org/mailman/listinfo/grass-dev">https://lists.osgeo.org/mailman/listinfo/grass-dev</a><br>>>><br>>><br>>><br>><br>><br>> _______________________________________________<br>> grass-dev mailing list<br>> <a href="mailto:grass-dev@lists.osgeo.org">grass-dev@lists.osgeo.org</a><br>> <a href="https://lists.osgeo.org/mailman/listinfo/grass-dev">https://lists.osgeo.org/mailman/listinfo/grass-dev</a><br><br></div></div></div>