3d-интерполяция
Добавлено: 23 июл 2018, 09:39
Здравствуйте.
Есть : NetCDF-файл с результатами чужих модельных расчетов некой скалярной характеристики (температура) на неком наборе уровней (глубин). (По сути: "'этажерка" растров на каждом уровне в EPSG:4326 с регулярным шагом в 0.25 градусов по Lon/Lat).
Надо: эту температуру проинтерполировать в некую сетку (нерегулярную по Lon/Lat, но регулярную по азимутальной Ламберта - да пусть даже не в таковую сетку, а в любой иррегулярный набор точек с известной парой Lon/Lat для каждой точки). Горизонты (глубины) целевых slice-вертикалей также отличны от набора уровней в исходном наборе модельных расчётов. Ну, эту вертикальную интерполяцию можно уже сделать и после горизонтальной.
Как, из R не выходя, это сделать?
gdal_wrap'ить не перепроецируя - невежественное кощунство... но, желательно,-то, чтобы в этих целевых slice-точках проинтерполированные значения оказались "родными" для этой исходной модели...
Когда-то, в несколько др. аспекте, это уважаемый gamm освещал, - но, что-то этот пост не нашёл я с ходу.
Да, ещё: линейные размеры [метры] каждой ячейки модели вдоль трёх осей - также есть в специальном файле, прилагаемом к модели.
MBA, akimа - инструментарий из этой "оперы"?...
Есть : NetCDF-файл с результатами чужих модельных расчетов некой скалярной характеристики (температура) на неком наборе уровней (глубин). (По сути: "'этажерка" растров на каждом уровне в EPSG:4326 с регулярным шагом в 0.25 градусов по Lon/Lat).
Надо: эту температуру проинтерполировать в некую сетку (нерегулярную по Lon/Lat, но регулярную по азимутальной Ламберта - да пусть даже не в таковую сетку, а в любой иррегулярный набор точек с известной парой Lon/Lat для каждой точки). Горизонты (глубины) целевых slice-вертикалей также отличны от набора уровней в исходном наборе модельных расчётов. Ну, эту вертикальную интерполяцию можно уже сделать и после горизонтальной.
Как, из R не выходя, это сделать?
gdal_wrap'ить не перепроецируя - невежественное кощунство... но, желательно,-то, чтобы в этих целевых slice-точках проинтерполированные значения оказались "родными" для этой исходной модели...
Когда-то, в несколько др. аспекте, это уважаемый gamm освещал, - но, что-то этот пост не нашёл я с ходу.
Да, ещё: линейные размеры [метры] каждой ячейки модели вдоль трёх осей - также есть в специальном файле, прилагаемом к модели.
MBA, akimа - инструментарий из этой "оперы"?...