Полигон, обрисовывающий в общих чертах географический четырехугольник
[lat,lon] = outlinegeoquad(latlim,lonlim,dlat,dlon)
[lat,lon] = outlinegeoquad(latlim,lonlim,dlat,dlon) создает полигон, который прослеживает схему географического четырехугольника, заданного latlim и lonlim. Такой полигон может быть полезен для отображения четырехугольника графически, особенно на проекции, где меридианы и/или параллели не проектируют к прямым линиям. latlim является двухэлементным вектором формы: [southern-limit northern-limit] и lonlim являются двухэлементным вектором формы: [western-limit eastern-limit]. dlat является положительной скалярной величиной, которая задает минимальную вершину, располагающую с интервалами в градусах, чтобы быть примененной вдоль меридианов, которые связали восточные и западные ребра четырехугольника. Аналогично, dlon является положительной скалярной величиной, которая задает минимальную вершину, располагающую с интервалами в градусах долготы, которая будет применена вдоль параллелей, которые связали северные и южные ребра четырехугольника. Выходные параметры lat и lon содержат вершины простого закрытого полигона с по часовой стрелке упорядоченным расположением вершины.
Отобразите основы трех географических четырехугольников, имеющих совсем другие качества сверху простой основной карты:
figure('Color','white')
axesm('ortho','Origin',[-45 110],'frame','on','grid','on')
axis off
load coastlines
geoshow(coastlat, coastlon)
% Quadrangle covering Australia and vicinity
[lat, lon] = outlinegeoquad([-45 5],[110 175],5,5);
geoshow(lat,lon,'DisplayType','polygon','FaceAlpha',0.5);
% Quadrangle covering Antarctic region
antarcticCircleLat = dms2degrees([-66 33 39]);
[lat, lon] = outlinegeoquad([-90 antarcticCircleLat], ...
[-180 180],5,5);
geoshow(lat,lon,'DisplayType','polygon', ...
'FaceColor','cyan','FaceAlpha',0.5);
% Quadrangle covering nominal time zone 9 hours ahead of UTC
[lat, lon] = outlinegeoquad([-90 90], 135 + [-7.5 7.5], 5, 5);
geoshow(lat,lon,'DisplayType','polygon', ...
'FaceColor','green','FaceAlpha',0.5);
Все углы ввода и вывода находятся в модулях степеней. Выберите довольно маленькое значение для dlat (несколько градусов, возможно) при использовании проекции с кривыми меридианами или изогнутыми параллелями.
Чтобы постараться не интерполировать дополнительные вершины вдоль меридианов или параллелей, установите dlat или dlon к значению Inf.
Вставка дополнительных вершин подавлена в полюсах (то есть, если latlim(1) == -90 или latlim(2) == 90. Если lonlim соответствует ширине четырехугольника точно 360 градусов (lonlim == [-180 180], например), то это покрывает всю широтную зону и включает две отдельных, NaN-разделенных части, если также
latlim(1) == -90 или latlim(2) == 90, так, чтобы только одна часть была необходима — полигон, который следует за параллелью по часовой стрелке вокруг одного из полюсов.
latlim(1) == -90 и latlim(2) == 90, так, чтобы четырехугольник охватил целую планету. В этом случае четырехугольник не может быть представлен полигоном долготы широты, и ошибка заканчивается.