Импортируйте данные о вертикальном изменении, найдите координаты импортированных данных, затем отобразите данные как поверхность.

Сначала импортируйте данные о вертикальном изменении [1] для области вокруг Горы Вашингтон как массив и объект ссылки ячеек карты. Подготовьте данные к графическому выводу, заменив пропускающий значения данных с 0. Затем найдите координаты каждого элемента в массиве.

[Z,R] = readgeoraster('MtWashington-ft.grd','OutputType','double');
info = georasterinfo('MtWashington-ft.grd');
Z = standardizeMissing(Z,info.MissingDataIndicator);

[X,Y] = worldGrid(R);

Отобразите данные о вертикальном изменении как поверхность с помощью соответствующей палитры. Установите EdgeColor свойство к 'none' удалить mesh из поверхности.

surf(X,Y,Z,'EdgeColor','none')
demcmap(Z)
xlabel('x (meters)')
ylabel('y (meters)')
zlabel('Elevation (feet)')

[1] Данные о вертикальном изменении, используемые в этом примере, являются любезностью американской Геологической службы.

Создайте объект ссылки ячеек карты для растра 3 на 4 с x значениями в области значений [7000, 7400] метры и y значения в области значений [2700, 3300] метры. Получите координаты растровых элементов и возвратите их как векторы-строки.

R = maprefcells([7000 7400],[2700 3300],[3 4]);
[X,Y] = worldGrid(R,'gridvectors')

X = 1×4

        7050        7150        7250        7350

Y = 1×3

        2800        3000        3200

Если вы не задаете второй аргумент как 'gridvectors', затем worldGrid функция возвращает 2D массивы по умолчанию.

[xFull,yFull] = worldGrid(R)

xFull = 3×4

        7050        7150        7250        7350
        7050        7150        7250        7350
        7050        7150        7250        7350

yFull = 3×4

        2800        2800        2800        2800
        3000        3000        3000        3000
        3200        3200        3200        3200