2017-07-25 54 views
0

我在做什麼我使用R.我從NetCDF文件創建一個光柵。這是來自NAM模型的天氣預報產品。什麼是從這個netcdf文件lambert適形投影proj4string

什麼問題我不知道如何在創建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 

而且我在​​網格信息。但是,我不知道這些參數意味着什麼,以及如何將它們匹配到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 

非常感謝你。

回答

0

我想回答我自己的問題,因爲它看起來沒什麼興趣。

我已經創建了a gist on Github來解決問題。