【发布时间】: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
【问题讨论】: