Извлеките значения сетки данных для заданных местоположений
val = ltln2val(Z, R, lat, lon)
val = ltln2val(Z, R, lat, lon, method)
val = ltln2val(Z, R, lat, lon) интерполирует обычную сетку данных Z со ссылкой на векторный R в точках, заданных векторами широты и долготы, lat и lon. R может быть географическим растровым объектом ссылки, вектором ссылки или матрицей ссылки.
Если R является географическим растровым объектом ссылки, его свойство RasterSize должно быть сопоставимо с size(Z).
Если R является вектором ссылки, это должно быть 1 3 с элементами:
[cells/degree northern_latitude_limit western_longitude_limit]
Если R является матрицей ссылки, это должно быть 3 2 и преобразовать растровые индексы строки и столбца к или от географических координат согласно:
[lon lat] = [row col 1] * R
Если R является матрицей ссылки, он должен задать (невращательный, нескошенный) отношение, в котором каждый столбец сетки данных падает вдоль меридиана и каждой строки падения вдоль параллели. Интерполяция ближайшего соседа используется по умолчанию. NaN возвращен для точек вне пределов сетки или для которого lat или lon содержат NaN. Все углы находятся в модулях степеней.
val = ltln2val(Z, R, lat, lon, method) где method задает тип интерполяции: 'bilinear' для линейной интерполяции, 'bicubic' для кубичной интерполяции или 'nearest' для самой близкой соседней интерполяции.
Найдите повышения в topo сопоставленными с тремя европейскими городами — Милан, Берн и Прага (повышения topo исчисляются в метрах):
load topo % The city locations, [Milan Bern Prague] lats = [45.45; 46.95; 50.1]; longs = [9.2; 7.4; 14.45]; elevations = ltln2val(topo,topolegend,lats,longs) elevations = 313 1660 297