Получение высоты геоида и объекта привязки географических проводок из модели геоида EGM96. Затем создайте график 3-D контура с 40 уровнями контура, используя проекцию Миллера. Устранение лишнего пробела с помощью tightmap.

[N,R] = egm96geoid;
axesm miller
contour3m(N,R,40,'LineColor','k')
tightmap

Добавьте геоид как поверхность, увеличьте пропорции осей с помощью daspectи просмотрите карту в 3-D.

hold on
geoshow(N,R,'DisplayType','surface')
daspect([1 1 50])
view(3)

В этом примере показано, как очертить на карте оси топографии и батиметрии Южной Азии и северной части Индийского океана с интервалом контура 500 метров.

Загрузите растровые данные отметки и объект ссылки на географические ячейки. Обрезать данные в районе Южной Азии.

load topo60c
latlim = [0 50];
lonlim = [35 115];
[Z,R] = geocrop(topo60c,topo60cR,latlim,lonlim);

Отображение данных в виде контурного графика с использованием стандартной конформной конической проекции Ламберта.

figure
axesm('lambertstd','MapLatLimit',latlim,'MapLonLimit',lonlim)
tightmap
axis off
contour3m(Z,R,'k','LevelStep',500)

Отображение топографии в виде поверхности. Применение карты цветов, соответствующей данным отметок, с помощью demcmap функция. Затем добавьте заголовок.

geoshow(Z,R,'DisplayType','surface')
demcmap(Z)
title({'South Asia Topography and Bathymetry', ...
   'with 500 m Contours'});

Просмотр отображения в 3-D.

set(gca,'DataAspectRatio',[1 1 40000])
view(3)