导入栅格堆栈会返回较大的NA值

时间:2020-06-03 14:25:33

标签: r extract raster na

我在加载所需的包后导入了栅格堆栈。在调用head()和summary()时,我发现R在导入时以某种方式将栅格读取为NA。当我运行hist()时,我能够看到像素值的分布,因此,数据肯定存在。我已经检查了其他支持链接:NAs returned when extracting values from a rasterImport raster to R - values appear as NAs?。但我无法理解出了什么问题。

library(raster)
x <- brick("raster.tiff")
x
#class      : RasterBrick 
#dimensions : 1610, 2640, 4440744, 4  (nrow, ncol, ncell, nlayers)
#resolution : 30, 30  (x, y)
#extent     : 694781.4, 773171.4, 1261280, 1311920  (xmin, xmax, ymin, ymax)
#crs        : +proj=utm +zone=43 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +towgs84=0,0,0 
#source     : file-path.grd 
#names      :        x.1,        x.2,          x.3,          x.4 
#min values :     275.2141,     357.1080,     370.5995,     592.5593 
#max values :     2122.801,     5525.976,     6406.745,     6870.159 

x[1]
#                x.1        x.2         x.3            x.4
#[1,]           NA           NA           NA           NA

summary(x)
#               x.1        x.2           x.3         x.4
#Min.        305.2251     425.7844     401.7832     641.0283
#1st Qu.     565.8536     796.9203     887.6340    1941.3184
#Median      652.2755     918.4586    1071.1154    2192.9951
#3rd Qu.     758.0400    1100.3450    1380.2045    2464.9850
#Max.       1287.9128    2654.8047    3272.8123    4073.0056
#NA's    1433504.0000 1433504.0000 1433504.0000 1433504.0000

另外,我是R的新手,甚至是stackoverflow的新手,所以也许我错过了一些东西。

0 个答案:

没有答案
相关问题