将 .adf 文件转换为 NetCDF

Question

我有一组包含空间分布数据的 .adf 文件。文件如下所示：

dblbnd.adf
hdr.adf
metadata.xml
prj.adf
sta.adf
w001001.adf
w001001x.adf

我可以将大文件 (w001001.adf) 作为光栅读入 R。这是栅格特征：

class      : RasterLayer 
dimensions : 1692, 3611, 6109812  (nrow, ncol, ncell)
resolution : 6, 6  (x, y)
extent     : -10833.68, 10832.32, -4713.677, 5438.323  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
source     : C:/Users/ahsa361/Documents/CEDS_Data/gridding-preprocessing/input/Corbett/w001001.adf 
names      : w001001 
values     : 0, 5637.337  (min, max)

我想将其转换为 180 度纬度的 NetCDF。 x 360 度。分辨率为 0.5 度。我该怎么做？

Answer 1

您的范围与您的坐标参考系统不匹配

#extent     : -10833.68, 10832.32, -4713.677, 5438.323  (xmin, xmax, ymin, ymax)
#crs        : +proj=longlat +datum=WGS84 +no_defs 

如果没有更多信息，它真的不可能为您提供太多帮助。 crs 很可能是错误的，您需要用正确的替换它。您是否将其设置为错误的值？ （这是此问题的最常见原因）。如果你知道它应该是什么，你可以这样设置

 crs(x) <- "correct crs definition"

范围也有可能是错误的，您也可以将其设置为正确的值。但根据你所说的，你不知道，所以你需要弄清楚这些事情，然后我们才能帮助如何解决你的问题。

要获得全局范围和接近 0.1 度的分辨率，您可以这样做

 extent(x) <- c(-180,180,-90,90)

鉴于您拥有的行数和列数，您可以：

library(raster)
raster(nrow=1692, ncol=3611)
#class      : RasterLayer 
#dimensions : 1692, 3611, 6109812  (nrow, ncol, ncell)
#resolution : 0.09969538, 0.106383  (x, y)
#extent     : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
#crs        : +proj=longlat +datum=WGS84 +no_defs