Закодируйте точки данных в обычную сетку данных
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