我在加载所需的包后导入了栅格堆栈。在调用head()和summary()时,我发现R在导入时以某种方式将栅格读取为NA。当我运行hist()时,我能够看到像素值的分布,因此,数据肯定存在。我已经检查了其他支持链接:NAs returned when extracting values from a raster, Import 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的新手,所以也许我错过了一些东西。