Преобразуйте долготу широты в строки и столбцы сетки данных
[row, col] = setpostn(Z, R, lat, lon)
indx = setpostn(...)
[row, col, indxPointOutsideGrid] = setpostn(...)
[row, col] = setpostn(Z, R, lat, lon) возвращает индексы строки и столбца обычной сетки данных Z для точек, заданных векторами lat и lonR может быть географический объект растровой привязки, вектор ссылки или матрица привязки.
Если 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 матрица привязки, она должна задать (невращательный, нескошенный) отношение, в котором каждый столбец сетки данных падает вдоль меридиана и каждой строки падения вдоль параллели. Точки, выходящие за пределы сетки, проигнорированы в row и col. Все входные углы в градусах.indx = setpostn(...) возвращает индексы Z соответствие точкам в lat и lon. Точки, выходящие за пределы сетки, проигнорированы в indx.
[row, col, indxPointOutsideGrid] = setpostn(...) возвращает индексы lat и lon соответствие точкам вне сетки. Эти точки проигнорированы в row и col.
Что является матричными координатами в topo из Денвера, Колорадо, в (39.7ºN, 105ºW)?
load topo [row,col] = setpostn(topo,topolegend,39.7,105) row = 130 col = 105