Векторные данные для географических объектов, которые охватывают полюс, должны столкнуться или пересечь линию перемены даты. В то время как тулбокс правильно отображает такие многоугольники, они могут вызвать проблемы для функций как пересечение многоугольника и Логические операции, которые работают с Декартовыми координатами. Когда эти многоугольники обработаны как Декартовы координаты, результаты пересечения линии перемены даты в побочном линейном сегменте и многоугольник, отображенный, когда закрашенной фигуре не заполнили внутреннюю часть правильно.

antarctica = shaperead('landareas', 'UseGeoCoords', true,...
    'Selector', {@(name) strcmp(name,'Antarctica'), 'Name'});
figure
plot(antarctica.Lon, antarctica.Lat)
ylim([-100 -60])

[latf, lonf] = flatearthpoly(antarctica.Lat', antarctica.Lon');
figure; mapshow(lonf, latf, 'DisplayType', 'polygon')
ylim([-100 -60])
xlim([-200 200])
axis square 

Разрывы координаты долготы в Антимеридиане могут перепутать операции присвоения на многоугольниках. Подготовить географические данные к использованию с polybool или для рендеринга закрашенной фигуры, сокращения многоугольники в Антимеридиане с flatearthpoly функция. flatearthpoly функция возвращает многоугольник с точками, вставленными, чтобы следовать за Антимеридианом до полюса, пересечь долготы в полюсе и возвратиться к Антимеридиану вдоль другого ребра Антимеридиана.

Создайте ортогональное представление Земли и постройте береговые линии на ней.

axesm ortho
setm(gca,'Origin', [60 170]); framem on; gridm on
load coastlines
plotm(coastlat,coastlon)

Сгенерируйте маленький круг, который охватывает Северный полюс, и окрасьте его в желтый.

[latc,lonc] = scircle1(75,45,30);
patchm(latc,lonc,'y')

Сгладьте маленький круг с помощью flatearthpoly функция.

[latf,lonf] = flatearthpoly(latc,lonc);

Постройте круг сокращения, который вы только сгенерировали как пурпурная линия.

plotm(latf,lonf,'m')

Сгенерируйте второй маленький круг, который не включает полюс.

[latc1, lonc1] = scircle1(20, 170, 30);

Сгладьте круг и постройте его как красную линию. Обратите внимание на то, что второй маленький круг, который не покрывает полюс, отсекается в две части вдоль Антимеридиана. Многоугольник для первого маленького круга построен в координатах на плоскости, чтобы проиллюстрировать его сглаженную форму. flatearthpoly функция принимает, что внутренняя часть сглаженного многоугольника находится в полушарии, которое содержит большинство его точек ребра. Таким образом многоугольник производится flatearthpoly не покрывает больше, чем полушарие.

[latf1,lonf1] = flatearthpoly(latc1,lonc1);
plotm(latf1,lonf1,'r')