Отобразите геоидные высоты из EGM96 модели геоида на 3-D земном шаре. Сначала получите геоидные высоты и географические проводки ссылки объекта. Загрузка данных о широте и долготе береговой линии.

[N,R] = egm96geoid;
load coastlines

Создайте систему координат для отображения 3-D глобуса с помощью axesm. Установите линию зрения для земного шара используя view. Отключите фон осей с помощью axis off. Затем отобразите геоидные высоты и данные о береговой линии.

axesm('globe','Grid','on')
view(60,60)
axis off
meshm(N,R)
plotm(coastlat,coastlon)

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

Создайте образец многоугольника, который содержит отверстие и опирается на поверхность земного шара. Для этого сгенерируйте вершины его внешних и внутренних контуров с помощью outlinegeoquad функция. В качестве первых двух аргументов укажите географические пределы, а в качестве следующих двух аргументов - интервалы между вершинами в степенях. Инвертируйте порядок внутренних граничных вершин, используя flip функция, поэтому они находятся в порядке против часовой стрелки.

[latE,lonE] = outlinegeoquad([-35 35],[-30 30],0.25,0.25);
[latI,lonI] = outlinegeoquad([-15 15],[-15 15],0.25,0.25);
latI = flip(latI);
lonI = flip(lonI);

Объедините вершины в один список путем разделения контуров с NaN значения.

lat = [latE NaN latI];
lon = [lonE NaN lonI];

Векторы lat и lon представление контуров многоугольника, содержащего отверстие. Отобразите контуры земного шара как заполненный многоугольник путем преобразования многоугольника в сетку данных.

Для этого создайте географические камеры ссылки объекта для земного шара и сетку таковых. Замените элементы сетки на многоугольник данных используя vec2mtx функция. Новая сетка содержит 0 с для указания внутренней области многоугольника, 1с для указания контуров и 2с для указания внешней области многоугольника.

R = georefcells([-90 90],[-180 180],0.25,0.25);
V = ones(R.RasterSize);
[V,R] = vec2mtx(lat,lon,V,R,'filled');

Создайте глобус с помощью axesm функция. Отобразите сетку данных как изображение с помощью geoshow функция. Отрегулируйте палитру так, чтобы внутренняя область многоугольника была фиолетовой, а внешняя - белой. Измените линию зрения камеры с помощью view функция, поэтому многоугольник отображается на ближней стороне земного шара.

axesm('globe','Grid','on')
geoshow(V,R,'DisplayType','texturemap')
colormap([0.5 0.5 0.8; 0 0 0; 1 1 1])
axis off
view(100,20)

Внешний вид многоугольников на земном шаре зависит от линии зрения камеры и прозрачности земного шара. Например, сделать глобус слегка прозрачным, используя alpha функция.

alpha(0.6)

При просмотре многоугольника с ближней стороны земного шара внешние граничные вершины появляются в порядке часовой стрелки. При просмотре многоугольника с обратной стороны земного шара внешние граничные вершины появляются в порядке против часовой стрелки. Когда вы поворачиваете глобус, так что многоугольник появляется как с ближней, так и с обратной стороны, тогда многоугольник кажется пересекающимся.