Drape Geoid Высоты над топографией

Отображение геоидных данных, выделенных по топографическому рельефу. В данном примере отобразите данные геоида как атрибут цвета вместо поверхности 3-D.

Загрузка топографических растровых данных и географических камер ссылки объекта. Получите геоидные высоты по вызову egm96geoid функция и установка ссылки объекта.

load topo60c
[N,R] = egm96geoid(topo60cR);

Создайте объект осей с помощью стереографической цилиндрической проекции Gall (перспективная проекция). Использование meshm постройте график цветного отображения изменений геоида, но задайте topo60c в качестве конечного аргумента, чтобы задать каждой геоидной сетке камеры высоту (z значение) соответствующей топографической сетки камеры. Низкие геоидные высоты показаны как синие высокие таковые как жёлтые.

axesm gstereo;
meshm(N,R,size(N),topo60c)

Для ссылки постройте морские береговые линии в черном цвете, поднимите их повышение до 1000 метров (достаточно высоко, чтобы очистить поверхность в их окрестностях), и разверните карту, чтобы заполнить систему координат.

load coastlines
plotm(coastlat,coastlon,'k')
zdatam(handlem('allline'),1000)
tightmap

Из-за вертикального вида и отсутствия подсветки топографический рельеф не виден, но является частью данных о поверхности рисунка. Вывести его путем преувеличения рельефа, а затем установить вид с юга-юго-востока.

daspectm('m',200); tightmap
view(20,35)

Удалите ограничивающий прямоугольник, осветите свет на поверхности (используя положение по умолчанию, смещение вправо от точки обзора) и выполните тонирование снова с помощью затенения Гуро.

ax = gca;
ax.Box = 'off';
camlight;
lighting Gouraud

Наконец, установите перспективу, чтобы немного сходиться (перспектива по умолчанию ортогональная). Заметьте, что геоид отражает рельеф крупных горных цепей, таких как Анды, Гималаи и Срединно-Атлантический хребет. Можно также увидеть, что большие площади с высокими или низкими геоидными высотами являются не просто результатом топографии.

ax.Projection = 'perspective';