光栅警告消息：在 cbind(m[i, ], vals) 中：结果的行数不是向量长度的倍数 (arg 2)

Question

我已经下载了这个文件。 'https://www.esrl.noaa.gov/psd/repository/entry/get/hpbl.mon.mean.nc?entryid=synth%3Ae570c8f9-ec09-4e89-93b4-babd5651e7a9%3AL05BUlIvTW9udGhsaWVzL21vbm9sZXZlbC9ocGJsLm1vbi5tZWFuLm5j'

现在我想打开它。

library(raster)

file <- "~/Desktop/input/hpbl/hpbl.mon.mean.nc"

data <- raster(file)

我收到以下警告：警告消息：

Warning message:In cbind(m[i, ], vals) : number of rows of result is not a multiple of vector length (arg 2).

但我想不通为什么。

谢谢！

Answer 1

出现警告是因为文件在坐标参考系统的描述中有一个小错误（"grid_mapping" in ncdf speak）。它有 standard_parallel = 50 50 而不是 standard_parallel = 50。您可以忽略该警告。在未来的版本中它将消失。在这种情况下，您可以通过这样设置来修复 CRS

crs(data) <- "+proj=lcc +x_0=5632642.22547 +y_0=4612545.65137 +lat_0=50 +lon_0=-107 +lat_1=50"

光栅警告消息：在 cbind(m[i, ], vals) 中：结果的行数不是向量长度的倍数 (arg 2)

Raster warning message: In cbind(m[i, ], vals) : number of rows of result is not a multiple of vector length (arg 2)

warnings

r

raster