Закодируйте точки данных в обычную сетку данных
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) сбрасывает определенные записи обычной сетки данных, Z. 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. Все углы находятся в модулях степеней.
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