Отображение геоидных данных, нанесенных поверх рельефа топографии. В этом примере данные геоида отображаются как атрибут цвета вместо 3-D поверхности.
Загрузите топографические растровые данные и объект ссылки на географические ячейки. Получить геоидные высоты, позвонив egm96geoid и задание ссылочного объекта.
load topo60c
[N,R] = egm96geoid(topo60cR);Создайте объект оси карты, используя стереографическую цилиндрическую проекцию Галла (перспективную проекцию). Использовать 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';