3d-интерполяция

Вопросы по статистическому пакету R. Не обязательно гео.
Ответить
nickleb
Гуру
Сообщения: 870
Зарегистрирован: 22 май 2010, 20:20
Репутация: 102

3d-интерполяция

Сообщение nickleb » 23 июл 2018, 09:39

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

gamm
Гуру
Сообщения: 2987
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 696
Ваше звание: программист
Откуда: Казань

Re: 3d-интерполяция

Сообщение gamm » 23 июл 2018, 12:47

если нужно послойно, то делаете из исходных слоев XYZ, и интерполируете на целевую сетку в пакете MBA. Акима - это линейная интерполяция в треугольниках, даже не билинейная - получите артефакты.

nickleb
Гуру
Сообщения: 870
Зарегистрирован: 22 май 2010, 20:20
Репутация: 102

Re: 3d-интерполяция

Сообщение nickleb » 23 июл 2018, 13:35

gamm писал(а):
23 июл 2018, 12:47
если нужно послойно, то делаете из исходных слоев XYZ, и интерполируете на целевую сетку в пакете MBA. Акима - это линейная интерполяция в треугольниках, даже не билинейная - получите артефакты.
gamm, спасибо! с MBA начну.

Ответить

Вернуться в «R»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость