Многоугольник, обрисовывающий в общих чертах географический четырехугольник
[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
, так, чтобы четырехугольник охватил целую планету. В этом случае четырехугольник не может быть представлен многоугольником долготы широты, и ошибка заканчивается.