Что такое proj4string для конформной проекции Ламберта из этого файла netcdf

Что я делаю, я использую R. Я создаю растр из файла NetCDF. Это прогноз погоды по модели NAM.

В чем проблема, я не знаю, как указать параметр proj4string при создании SpatialPoints. Например, я ищу что-то вроде этого:

+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

Большое спасибо.

1 ответ

Решение

Я хотел бы ответить на свой вопрос, так как, похоже, это не заинтриговало меня.

Я создал суть на Github для решения проблемы.

Другие вопросы по тегам