首页 > 解决方案 > 栅格值都是 NA

问题描述

我有一个正弦坐标系的叶面积指数(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

标签: rraster

解决方案


从您所展示的内容来看,没有任何迹象表明您的所有价值观都是NA. 而且,鉴于可以使用该对象制作正确的绘图,显然情况并非如此。您所显示的是第一个和最后一个单元格值为 NA。这很常见(极点处的 NA 值)。你真正的问题是什么?


推荐阅读