Отображение геоидных данных, выделенных по топографическому рельефу. В данном примере отобразите данные геоида как атрибут цвета вместо поверхности 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';