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)

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

警告消息：在cbind（m [i，]，vals）中：结果的行数不是向量长度（arg 2）的倍数。

但是我不知道为什么。

谢谢！

Answer 1

发生警告是因为文件在坐标参考系统的描述中有一个小错误（在ncdf中为“ grid_mapping”）。它具有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）

1 个答案: