【问题标题】:raster values are all NAs栅格值都是 NA
【发布时间】:2019-05-16 00:35:40
【问题描述】:

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

【问题讨论】:

    标签: r raster


    【解决方案1】:

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

    【讨论】:

      猜你喜欢
      • 2013-04-22
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      • 2018-12-01
      • 2015-09-27
      • 2021-09-29
      • 1970-01-01
      • 2019-11-15
      相关资源
      最近更新 更多