来自此netcdf文件的lambert conformal投影的proj4string是什么

时间:2017-07-25 22:24:12

标签: r raster projection netcdf

我正在做什么我正在使用R.我正在从NetCDF文件创建一个栅格。这是NAM模型的天气预报产品。

问题是什么我在创建SpatialPoints时不知道如何指定proj4string参数。例如,我正在寻找类似的东西:

+proj=aea +lat_1=20 +lat_2=60 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs

我在grid description处有网格信息。但是,我不知道这些参数的含义以及如何将它们与lat_0等匹配。

因此,要为从NetCDF文件创建的栅格创建正确的投影信息,如何指定NetCDF的投影?

print(nc_open('NAM_forecast.nc'))的部分输出如下:

float gridlat_0[xgrid_0,ygrid_0]   
            La1: 12.1899995803833
            Lo1: 226.541000366211
            Lov: 265
            Dx: 12.1909999847412
            Dy: 12.1909999847412
            Latin1: 25
            Latin2: 25
            units: degrees_north
            grid_type: Lambert Conformal (secant, tangent, conical or bipolar)
            long_name: latitude
            corners: 12.1899995803833
             corners: 14.3420877456665
             corners: 57.3284339904785
             corners: 54.5653381347656
float gridlon_0[xgrid_0,ygrid_0]   
            La1: 12.1899995803833
            Lo1: 226.541000366211
            Lov: 265
            Dx: 12.1909999847412
            Dy: 12.1909999847412
            Latin1: 25
            Latin2: 25
            units: degrees_east
            grid_type: Lambert Conformal (secant, tangent, conical or bipolar)
            long_name: longitude
            corners: -133.458999633789
             corners: -65.1255569458008
             corners: -49.4160041809082
             corners: -152.878616333008

非常感谢。

1 个答案:

答案 0 :(得分:0)

我想回答我自己的问题,因为它看起来并没有引起人们的兴趣。

我已经创建a gist on Github来解决问题。