此 netcdf 文件中用于兰伯特共形投影的 proj4string 是什么
What's the proj4string for lambert conformal projection from this netcdf file
我在做什么 我正在使用 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
非常感谢。
我想回答我自己的问题,因为它看起来并没有引起太多兴趣。
我已经创建了 a gist on Github 来解决这个问题。
我在做什么 我正在使用 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
非常感谢。
我想回答我自己的问题,因为它看起来并没有引起太多兴趣。
我已经创建了 a gist on Github 来解决这个问题。