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