Географические координаты растровых элементов
[, где lat,lon] = geographicGrid(R,gridOption)gridOption 'gridvectors', возвращает lat и lon как векторы. Координаты растрового элемента (i,j) (lat(i),lon(j)). Значение по умолчанию для gridOption 'fullgrid', который возвращает lat и lon как 2D массивы.
Импортируйте данные о вертикальном изменении [1] для области вокруг Южного Пика Валуна в Колорадо как массив и географический объект ссылки регистраций. Найдите координаты каждого элемента в массиве.
[Z,R] = readgeoraster('n39_w106_3arc_v2.dt1');
[lat,lon] = geographicGrid(R);Создайте карту с пределами широты и долготы, которые совпадают с пределами данных. Отобразите данные с помощью соответствующей палитры.
usamap(R.LatitudeLimits,R.LongitudeLimits) surfm(lat,lon,Z) demcmap(Z)
[1] Данные о вертикальном изменении, используемые в этом примере, являются любезностью американской Геологической службы.
Создайте географический объект ссылки ячеек для растра 3 на 4 со значениями широты в области значений [0, 30] степени и значения долготы в области значений [-20, 20] степени. Получите координаты растровых элементов и возвратите их как векторы-строки.
R = georefcells([0 30],[-20 20],[3 4]);
[lat,lon] = geographicGrid(R,'gridvectors')lat = 1×3
5 15 25
lon = 1×4
-15 -5 5 15
Если вы не задаете второй аргумент как 'gridvectors', затем geographicGrid функция возвращает 2D массивы по умолчанию.
[latFull,lonFull] = geographicGrid(R)
latFull = 3×4
5 5 5 5
15 15 15 15
25 25 25 25
lonFull = 3×4
-15 -5 5 15
-15 -5 5 15
-15 -5 5 15