只有正纬度和经度时导入netcdf在R从PaleoView可能是错误的投影

Question

我使用paleoview软件从过去的气候中下载了一些变量，包括15000BP-10000BP的平均温度（我可以根据要求上传文件，但至少是GB）。

主要问题是，当我读取光栅时，它只包含正纬度和正经度。我从文章中知道，这有一个2.5*2.5的分辨率。

我加载了光栅和ncdf4库，以使用光栅读取它

library(raster)
library(ncdf4)

当我使用以下代码阅读时

r <- raster("mean_temperature-15000BP-10000BP.nc",  varname = "14000BP-13000BP/13300BP")

我得到以下信息

r
class       : RasterLayer 
band        : 1  (of  12  bands)
dimensions  : 72, 144, 10368  (nrow, ncol, ncell)
resolution  : 1, 1  (x, y)
extent      : 0.5, 144.5, 0.5, 72.5  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
names       : Mean.Temperature 
z-value     : 1 
zvar        : 14000BP-13000BP/13300BP 

正如你所看到的，范围只是正的，但我知道它有整个世界的数据，当我绘制地图时，我得到了以下图像：

plot(r)

这清楚地显示了极圈温度极低，南极洲比北极大的情况下的预期模式。

很明显，这里的proj = long lat datum = WGS 84 ellps = WGS 84 tow GS 84 = 0，0，0是错误的，我想如果我算出实际坐标是多少。参考文献。我可以重新投影，得到正确格式的栅格

我尝试使用ncdf4包来尝试从图层中获取更多信息，这就是我所做的：

nc <- nc_open("mean_temperature-15000BP-10000BP.nc")

在使用nc_open函数读取层之后，我看到了变量的名称，以试图了解更多层，这里我显示了前10个

names(nc$var)[1:10]

[1] "window"                  "width"                   "decimals"                "months"                  "latitudes"               "longitudes"              "15100BP-15000BP/15100BP" "15100BP-15000BP/15099BP"
[9] "15100BP-15000BP/15098BP" "15100BP-15000BP/15097BP"

所以如果我继续寻找并转到我得到的纬度和经度名称：

ncatt_get(nc, attributes(nc$var)$names[5])
$units
[1] "degrees north"

和

ncatt_get(nc, attributes(nc$var)$names[6])
$units
[1] "degrees east"

关于如何重新投影此栅格以获得从 -90、90 和经度从 -180 到 180 的任何想法，您在 proj=longlat ellps=WGS84 基准面=WGS84 no_defs towgs84=0，0，0 投影中

Answer 1

要解决这个问题，您可以设置范围:

library(raster)
r <- raster("mean_temperature-15000BP-10000BP.nc",  varname = "14000BP-13000BP/13300BP")
extent(r) <- c(-180,180,-90,90)

但请注意，有时气候数据中的经度从0到360（我想不是这里）：

extent(r) <- c(0,360,-90,90)
r <- rotate(r)

要了解ncdf的更多内容，您可以

print(r)

当然，您只想这样做，因为您非常确定程度，并且您可能希望联系数据提供商报告此问题。