Драпируйте высоты геоида по топографии

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

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

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

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

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

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

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';