Преобразуйте долготу широты в строки и столбцы сетки данных
[row, col] = setpostn(Z, R, lat, lon)
indx = setpostn(...)
[row, col, indxPointOutsideGrid] = setpostn(...)
[row, col] = setpostn(Z, R, lat, lon) возвращает индексы строки и столбца обычной сетки данных Z для точек, заданных векторами lat и lon. 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 является матрицей ссылки, он должен задать (невращательный, нескошенный) отношение, в котором каждый столбец сетки данных падает вдоль меридиана и каждой строки падения вдоль параллели. Точки, выходящие за пределы сетки, проигнорированы в 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