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

[nickleb]

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

[gamm]

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

[nickleb]

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

gamm, спасибо! с MBA начну.