为什么我的光栅图像不能用 levelplot 显示？

Question

我正在尝试叠加包含夏威夷上空气溶胶高度数据的 netcdf4 栅格。示例文件可用 here. 我感兴趣的变量是纬度、经度、时间和气溶胶高度。这是一些可重现的数据。

s1 <- data.frame(as.vector(lon), as.vector(lat), as.vector(ah))
s1
#    as.vector.lon. as.vector.lat. as.vector.ah.
#1       -127.45199      -79.15431            NA
#2       -126.99632      -79.16919            NA
#3       -126.54577      -79.18321            NA
#4       -126.10027      -79.19641            NA
#5       -125.65974      -79.20880            NA
#6       -125.22412      -79.22042            NA
#7       -124.79333      -79.23129            NA

crsLatLon <- "+proj=longlat +datum=WGS84"
ex <- extent(c(-180,180,-90,90))

#empty raster with 0.1 degree resolution
pmraster <- raster(ncol=360*10, nrow=180*10, crs=crsLatLon,ext=ex)

#fills the empty raster with values from dataframe, s1
pmraster <- rasterize(s1[,1:2], pmraster, s1[,3], fun=mean, na.rm=T)

show(pmraster)
#class      : RasterLayer 
#dimensions : 1800, 3600, 6480000  (nrow, ncol, ncell)
#resolution : 0.1, 0.1  (x, y)
#extent     : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
#crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
#source     : r_tmp_2020-06-26_114048_16840_75885.grd 
#names      : layer 
#values     : 0.1314196, 9424.118  (min, max)

#specifies region over Hawaii
exHI <- extent(c(-180,-140,10,30))

levelplot(crop(pmraster,exHI))

#Error: $ operator is invalid for atomic vectors
#In addition: Warning messages:
#1: In min(x) : no non-missing arguments to min; returning Inf
#2: In max(x) : no non-missing arguments to max; returning -Inf
#3: In min(x) : no non-missing arguments to min; returning Inf
#4: In max(x) : no non-missing arguments to max; returning -Inf

任何人都可以帮助解释为什么我收到此错误消息以及我如何继续生成所需的光栅图像？提前致谢！

Answer 1

这是一个最小且可重现的示例：

library(raster)
library(lattice)
f <- system.file("external/test.grd", package="raster")
r <- raster(f)
levelplot(r)
#Error in UseMethod("levelplot") : 
#  no applicable method for 'levelplot' applied to an object of class "c('RasterLayer', #'Raster', 'BasicRaster')"

lattice 包方法不知道 Raster* 是什么。因此你不能像这样使用 levelplot 。好消息是 rasterVis 包为 Raster* 对象实现了 lattice 方法；所以你需要做的就是

library(rasterVis)
levelplot(r)

并参阅 ?rasterVis 了解更多使用 levelplot 的方法

ssplot 也是基于 levelplot

spplot(r)

Answer 2

在对不同的文件进行尝试后，我清楚地知道在我隔离的特定区域没有可用数据，这让我感到困惑。

为什么我的光栅图像不能用 levelplot 显示？

Why won't my raster image display with levelplot?

raster

geospatial

satellite

rasterizing