将栅格导入R值显示为NA？

Question

为什么R为导入的栅格的值显示NA？

#load libraries
library(raster)

#set working directory
setwd("C:/Users/Documents/myData")

#import raster data 
r <- raster("myRasterdata.tif")

返回：     class：RasterLayer

dimensions: 34073, 19663, 669977399  (nrow, ncol, ncell)

resolution  : 30, 30  (x, y)

extent : -969193.5, -379303.5, 277627.9, 1299818  (xmin, xmax, ymin, ymax)

coord. ref. : +proj=lcc +lat_1=49 +lat_2=77 +lat_0=49 +lon_0=-95 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0 

data source : C:/Users/Documents/myData/myRasterdata.tif 

names       : dem

values      : 0, 1857  (min, max)

正如您所看到的，它可以很好地导入文件，并显示存在0到1857之间的值，但是当我使用函数head（r）时，我会返回所有NA值。

同样，当我导入许多栅格并构建堆栈时，同样当我尝试读取GTAL函数时。

因此，当我尝试使用多边形图层执行“提取”功能时，我会返回所有NA值。

Answer 1

您可以使用trim功能删除NA值。

r<-trim(r)