Здравствуйте.
Есть : NetCDF-файл с результатами чужих модельных расчетов некой скалярной характеристики (температура) на неком наборе уровней (глубин). (По сути: "'этажерка" растров на каждом уровне в EPSG:4326 с регулярным шагом в 0.25 градусов по Lon/Lat).
Надо: эту температуру проинтерполировать в некую сетку (нерегулярную по Lon/Lat, но регулярную по азимутальной Ламберта - да пусть даже не в таковую сетку, а в любой иррегулярный набор точек с известной парой Lon/Lat для каждой точки). Горизонты (глубины) целевых slice-вертикалей также отличны от набора уровней в исходном наборе модельных расчётов. Ну, эту вертикальную интерполяцию можно уже сделать и после горизонтальной.
Как, из R не выходя, это сделать?
gdal_wrap'ить не перепроецируя - невежественное кощунство... но, желательно,-то, чтобы в этих целевых slice-точках проинтерполированные значения оказались "родными" для этой исходной модели...
Когда-то, в несколько др. аспекте, это уважаемый gamm освещал, - но, что-то этот пост не нашёл я с ходу.
Да, ещё: линейные размеры [метры] каждой ячейки модели вдоль трёх осей - также есть в специальном файле, прилагаемом к модели.
MBA, akimа - инструментарий из этой "оперы"?...
3d-интерполяция
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
-
- Гуру
- Сообщения: 4069
- Зарегистрирован: 15 окт 2010, 08:33
- Репутация: 1064
- Ваше звание: программист
- Откуда: Казань
Re: 3d-интерполяция
если нужно послойно, то делаете из исходных слоев XYZ, и интерполируете на целевую сетку в пакете MBA. Акима - это линейная интерполяция в треугольниках, даже не билинейная - получите артефакты.
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 0 гостей