flatearthpoly
Вставьте точки вдоль линии перемены даты, чтобы орудовать шестами
Синтаксис
[latf,lonf] = flatearthpoly(lat,lon)
[latf,lonf] = flatearthpoly(lat,lon,longitudeOrigin)
Описание
[latf,lonf] = flatearthpoly(lat,lon) обрезки NaN-разделенные полигоны, заданные векторами широты и долготы lat и lon к пределам [-180 180] в долготе и [-90 90] в широте, вставляя прямые сегменты вперед +/-меридианы на 180 градусов и в полюсах. Вводы и выводы в градусах.
[latf,lonf] = flatearthpoly(lat,lon,longitudeOrigin) сосредотачивает пределы долготы на долготе, заданной скалярной долготой longitudeOrigin.
Примеры
Векторные данные для географических объектов, которые охватывают полюс, неизбежно столкнутся или пересекут линию перемены даты. В то время как тулбокс правильно отображает такие полигоны, они могут вызвать проблемы для функций как пересечение полигона и Логические операции, которые работают с Декартовыми координатами. Когда эти полигоны обработаны как Декартовы координаты, результаты пересечения линии перемены даты в побочном линейном сегменте и полигон, отображенный, когда закрашенной фигуре не заполнили внутреннюю часть правильно.
antarctica = shaperead('landareas', 'UseGeoCoords', true,...
'Selector', {@(name) strcmp(name,'Antarctica'), 'Name'});
figure; plot(antarctica.Lon, antarctica.Lat); ylim([-100 -60])
Полигоны могут быть переформатированы более соответственно для Декартовых координат с помощью функции flatearthpoly. Результат напоминает отображение карты на цилиндрической проекции. Полигон соответствует линии перемены даты, раскрывается к полюсу, развертывается через долготы в полюсе и следует за линией перемены даты до другой стороны пересечения линии перемены даты.
[latf, lonf] = flatearthpoly(antarctica.Lat', antarctica.Lon'); figure; mapshow(lonf, latf, 'DisplayType', 'polygon') ylim([-100 -60]) xlim([-200 200]) axis square
Советы
Топология полигона для входных векторов должна быть допустимой. Это означает, что вершины для внешних звонков (основной полигон или “островные” полигоны) должны быть в по часовой стрелке порядке, и любые внутренние звонки (“озера”) должны запуститься в против часовой стрелки порядке для функции, чтобы работать правильно. Можно использовать функцию ispolycw, чтобы проверять, соответствуют ли lat, векторы lon этому критерию, и poly2cw и функциям poly2ccw, чтобы исправить какой-либо что выполнение в неправильном направлении.