Закодируйте точки данных в обычную сетку данных
Z = imbedm(lat, lon, value, Z, R)
Z = imbedm(lat, lon, value, Z, R, units)
[Z, indxPointOutsideGrid] = imbedm(...)
Z = imbedm(lat, lon, value, Z, R) сбрасывает определенные записи обычной сетки данных, ZR может быть географический объект растровой привязки, вектор ссылки или матрица привязки.
Если 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. Все углы находятся в модулях степеней.
Z = imbedm(lat, lon, value, Z, R, units) задает модули векторов lat и lon, где units любой допустимый угловой модульный вектор символов ('degrees' по умолчанию).
[Z, indxPointOutsideGrid] = imbedm(...) возвращает индексы lat и lon соответствие точкам вне сетки в переменной indxPointOutsideGrid.
Создайте простую карту сетки и встройте новые значения в нее:
Z = ones(3,6)
Z =
1 1 1 1 1 1
1 1 1 1 1 1
1 1 1 1 1 1
refvec = [1/60 90 -180]
refvec =
0.0167 90.0000 -180.0000
newgrid = imbedm([23 -23], [45 -45],[5 5],Z,refvec)
newgrid =
1 1 1 1 1 1
1 1 5 5 1 1
1 1 1 1 1 1