Question

我正在尝试在夏威夷上覆盖包含Aerosol Height数据的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包为lattice对象实现了Raster*方法；所以你要做的就是

library(rasterVis)
levelplot(r)

有关更多使用levelplot的方法，请参见?rasterVis

ssplot也是基于levelplot构建的

spplot(r)

Answer 2

在其他文件上尝试此操作后，对我来说很清楚，我隔离的特定区域没有可用的数据，这导致了我的困惑。

为什么我的栅格图像不显示在levelplot上？

2 个答案: