Unterschied in den Rasterwerten zwischen R und QGIS

Question

Ich möchte fragen, was wahrscheinlich eine grundlegende Frage über die Art und Weise ist, in der R und QGIS Rasterdateien importieren.

Ich habe ein einzelnes Band-Raster. Wenn ich es in R zu importieren, das „Raster“ -Funktion des raster-Paket verwenden, ich diesen Bereich von Pixelwerten erhalten:

class  : RasterLayer 
dimensions : 10980, 10980, 120560400 (nrow, ncol, ncell) 
resolution : 10, 10 (x, y) 
extent  : 6e+05, 709800, 5590200, 5700000 (xmin, xmax, ymin, ymax) 
coord. ref. : +proj=utm +zone=31 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : /data/MTDA/CGS_S2_RADIOMETRY/2017/10/15/S2B_20171015T104525Z_31UFS_TOC_V100/S2B_20171015T104525Z_31UFS_TOC-B02_10M_V100.tif 
names  : S2B_20171015T104525Z_31UFS_TOC.B02_10M_V100 
values  : -32768, 32767 (min, max) 

Wenn ich diese Schicht in einem Raster-Ziegel stapeln, erhalten ich diese min-max-Werte wenn ich das gleiche Raster in QGIS laden, ist der Minimalwert 228 und der Maximalwert ist 907 (I diese Werte mit den Optionen „Umfang: Voll“ berechnet

class  : RasterLayer 
band  : 2 (of 11 bands) 
dimensions : 10980, 10980, 120560400 (nrow, ncol, ncell) 
resolution : 10, 10 (x, y) 
extent  : 6e+05, 709800, 5590200, 5700000 (xmin, xmax, ymin, ymax) 
coord. ref. : +proj=utm +zone=31 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : /tmp/Rtmp882dZS/raster/r_tmp_2017-11-10_172819_11532_86514.grd 

names  : S2B_20171015T104525Z_31UFS_TOC.B02_10M_V100 
values  : -1129, 9994 (min, max) 

jedoch und „Genauigkeit: ist (langsamer:) "

Also, woher kommen diese Unterschiede? Ich verstehe nicht genau, was R und QGIS tun ...

Answer 1

Am Ende fand ich, was ist der Unterschied! Wenn ich nach R frage, bekomme ich die echten Min/Max-Werte. QGIS berechnet stattdessen die Min/Max-Werte mit einer kumulativen Anzahl. Wenn ich "Load min/max values" (im Fenster Raster Properties) auf "Min/Max" setze, habe ich dieselben Werte wie R angezeigt.

Answer 2

Für das erste Objekt sind die Min- und Max-Werte nicht bekannt, da die Datei sie nicht (oder nicht korrekt) bereitstellt. Mit RasterLayerr können Sie tun

r <- setMinMax(r) 

Um zu sehen, was sie sind. Wenn sie nicht wie bei der zweiten Ebene, die du zeigst, werden, hast du wahrscheinlich etwas vermischt. Nachdem diese Dinge geklärt wurden, könnte es nützlich sein, sie mit QGIS zu vergleichen. Dafür müssten Sie wahrscheinlich eine Beispieldatei bereitstellen.