R和QGIS之间栅格值的差异答案

【问题标题】：Difference in raster values between R and QGISR和QGIS之间栅格值的差异
【发布时间】：2018-04-26 07:43:37
【问题描述】：

我想问一下关于 R 和 QGIS 导入栅格文件的方式的一个基本问题。

我有一个单波段光栅。当我将它导入 R 时，使用 raster 包的“raster”功能，我得到这个像素值范围：

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)

当我将此层堆叠在光栅砖中时，我会得到这些最小值-最大值：

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)

但是，如果我在 QGIS 中加载相同的栅格，则最小值为 228，最大值为 907（我使用“Extent：Full”和“Accuracy: Actual (slower)”选项计算了这些值。

那么，这些差异从何而来？我不明白 R 和 QGIS 在做什么......

【问题讨论】：

标签： pixel raster qgis

【解决方案1】：

对于第一个对象，最小值和最大值是未知的，因为文件没有提供它们（或不正确）。使用RasterLayer r 你可以做到

 r <- setMinMax(r)

看看它们是什么。如果它们与您展示的第二层不同，那么您可能混淆了某些东西。在澄清了这些事情之后，与 QGIS 进行比较可能会有用。为此，您可能需要提供一个示例文件。

【讨论】：

【解决方案2】：

最后，我发现有什么不同！当询问 R 时，我得到了真正的最小值/最大值。相反，QGIS 使用累积计数计算最小值/最大值。当我将“加载最小/最大值”（在“光栅属性”窗口中）设置为“最小/最大值”时，我得到了与 R 显示的相同的值。

【讨论】：