Question

如何在netcdf文件中进行逆坐标转换？我的网格有75个经度值和36个纬度值：

nc<-create.n("filename.nc")
#Dimentions
dim.def.nc(nc,"lon",75)
dim.def.nc(nc,"lat",36)
dim.def.nc(nc,"time",365)
#Vars
var.def.nc(nc,"Observation","NC_FLOAT", c(1,0,2))
var.def.nc(nc,"lon","NC_FLOAT", c(0))
var.def.nc(nc,"lat","NC_FLOAT", c(1))
var.def.nc(nc,"time","NC_FLOAT", c(2))
(...)

根据unidata中的文档，应该可以让netcdf进行从（lat，lon）到（x，y）的逆变换，但我不知道我怎么能执行此操作。我想将我的lat long网格转换为Lambert Conformal网格。

Answer 1

这是我继续重新投影netCDF文件的方式：基本上我抓住了我想要的经度，纬度和数据并创建了一个shapefile，然后我可以使用包rgdal，{{1 }和sp。（此处的示例使用NOAA下载的数据here）

maptools

netcdf文件反投影

1 个答案: