Форумы GIS-Lab.info

nickleb писал(а):Здравствуйте.
Кто-нибудь решал в R Cran следующую задачу:
регулярный по Latitude/Longitude NetCDF-файл
перевести в растр и по этому растру
получить значения в узлах нерегулярной сетки. (?)
Поделитесь опытом, пожалуйста.

grid.XX<_grid.X0+grid.nX*grid.dX
grid.YY<_grid.Y0+grid.nY*grid.dY
grid.XY<-expand.grid(grid.XX,grid.YY)

col<-floor((target$X-grid.X0)/grid.dX)
row<-floor((target$Y-grid.Y0)/grid.dY)

rX<-target$X-(grid.X0+grid.dX*col)
rY<-target$Y-(grid.Y0+grid.dY*col)

Z10<-(row+1-1)*grid.nX+col    Z11<-(row+1-1)*grid.nX+col+1

Z00<-(row-1)*grid.nX+col        Z01<-(row-1)*grid.nX+col+1

gamm писал(а):

nickleb писал(а): 2) Интерполируем в узлы прямоугольной регулярной сетки (это растр, но в табличном виде). Границы берем с запасом, чтобы целевая область уместилось. Интерполируем сплайнами из пакета MBA
Сетку строим так

Код: Выделить всё

grid.XX<_grid.X0+grid.nX*grid.dX
grid.YY<_grid.Y0+grid.nY*grid.dY
grid.XY<-expand.grid(grid.XX,grid.YY)

[/code]

gamm писал(а):Что-то мне подсказывает, что в surfer вы никакой вариограммы не строили, а использовали линейную вариограмму по умолчанию. Т.е. строили сплайн. Поэтому не мудрите, и используте MBA, это тоже сплайны, очень быстрые.

[ Сообщение с мобильного устройства ]

gamm писал(а):

nickleb писал(а):

Интерполируем сплайнами из пакета MBA

nickleb писал(а):а через MBА-package получить в ASCII-файл узлы проинтерполированного 3D-положение изолинии по исходному многоуровенному NetCDF возможно?

gamm писал(а):

nickleb писал(а):а через MBА-package получить в ASCII-файл узлы проинтерполированного 3D-положение изолинии по исходному многоуровенному NetCDF возможно?

ничего не понял - ни что есть на входе, ни что нужно получить на выходе ... NetCDF - сетка, изолиний там нет; и непонятно, к чему относится "многоуровневый"; и что на выходе нужно - сетку? нужны структуры данных, а не форма представления (ASCII)

nickleb писал(а):на входе NetCDF с Temp[Lat/Lon/Depth/] в географии... надо получить глубину изотермы ноль градусов по этому параллелепипеду в задаваемом наборе точек с отличными от исходного набора Lat/Lon-координатами...