Question

我有一个具有正弦坐标系的叶面积指数（MODIS产品）的栅格。我使用raster（）函数打开它，它显示的值在0到254之间，这正是我在ArcMap中打开它时所看到的。绘制栅格也可以正确显示它。但是，当我调用其中的values（）时，它将所有单元格都转换为NA。我还尝试了head（），tail（），getValues（），而我所得到的只是NA。这是我的代码：

A2000057<- raster ('MODIS/LAI/clip/A2000057.Lai_500m.tif')

A2000057

class       : RasterLayer 
dimensions  : 3038, 6292, 19115096  (nrow, ncol, ncell)
resolution  : 463.3127, 463.3127  (x, y)
extent      : 12449213, 15364376, -4359773, -2952229  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs 
data source : C:GIS work/MODIS/LAI/clip/A2000057.Lai_500m.tif 
names       : A2000057.Lai_500m 
values      : 0, 254  (min, max)

values(A2000057)

   [1] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
  [35] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
  [69] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
.
.
.
.
 [953] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
 [987] NA NA NA NA NA NA NA NA NA NA NA NA NA NA

 [ reached getOption("max.print") -- omitted 19114096 entries ]

head(value)

[1] NA NA NA NA NA NA

tail(value)

[1] NA NA NA NA NA NA

NAvalue(A2000057)

[1] -Inf

Answer 1

从显示的内容来看，没有什么表明您的所有值都是NA。并且，鉴于可以使用该对象绘制正确的图，显然情况并非如此。您显示的是第一个和最后一个单元格值为NA。这很常见（极点处的NA值）。您的真实问题是什么？

栅格值均为NA

1 个答案: