Кодирование точек данных в регулярную сетку данных
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