Question

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

我有一个带状光栅。当我将它导入R时，使用＆＃34; raster＆＃34;光栅包的功能，我得到这个像素值范围：

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（我使用选项＆＃34计算这些值;范围：完全＆＃34;和＆＃34;准确度：实际（较慢）＆＃34;。

那么，这些差异来自哪里？我完全不明白R和QGIS在做什么......

Answer 1

对于第一个对象，由于文件未提供（或未正确），因此不知道最小值和最大值。使用RasterLayer r即可

 r <- setMinMax(r)

看看它们是什么。如果它们与你展示的第二层不一样，那么你可能会混淆一些东西。澄清这些内容后，与QGIS进行比较可能会有所帮助。为此，您可能需要提供示例文件。

Answer 2

最后，我发现有什么区别！当询问R时，我得到真实的最小/最大值。相反，QGIS使用累计计数来计算最小/最大值。当我设置＆＃34;加载最小/最大值＆＃34; （在光栅属性窗口中）到＆＃34;最小/最大＆＃34;我得到了R显示的相同值。

R和QGIS之间栅格值的差异

2 个答案: