Контурная сетка в локальной системе с результатами широты-долготы
[ возвращает линии и многоугольники, содержащие контурные линии и полигоны заливки контуров, соответственно. Эта функция является неграфической. Можно построить график значений возврата с помощью contourLines,contourPolygons]
= geocontourxy(X,Y,Z,lat0,lon0,h0)geoshow, при желании.
[___] = geocontourxy(___, задает пары "имя-значение", которые управляют аспектами операции. Имена параметров могут быть сокращены и не чувствительны к регистру.Name,Value)
Задайте координаты множества x и Y и создайте горизонтальные линии и горизонтальные полигоны.
X = -150000:10000:150000; Y = 0:10000:300000; [xmesh, ymesh] = meshgrid(X/50000, (Y - 150000)/50000); Z = 8 + peaks(xmesh, ymesh); lat0 = dm2degrees([ 21 18]); lon0 = dm2degrees([-157 49]); h0 = 300; levels = 0:2:18; [contourLines, contourPolygons] = geocontourxy(X,Y,Z,lat0,lon0,h0, ... 'LevelList',levels,'XYRotation',120)
contourLines =
8x1 geoshape vector with properties:
Collection properties:
Geometry: 'line'
Metadata: [1x1 struct]
Vertex properties:
(8 features concatenated with 7 delimiters)
Latitude: [1x329 double]
Longitude: [1x329 double]
Height: [1x329 double]
Feature properties:
ContourLevel: [2 4 6 8 10 12 14 16]
contourPolygons =
9x1 geoshape vector with properties:
Collection properties:
Geometry: 'polygon'
Metadata: [1x1 struct]
Vertex properties:
(9 features concatenated with 8 delimiters)
Latitude: [1x651 double]
Longitude: [1x651 double]
Height: [1x651 double]
Feature properties:
LowerContourLevel: [0 2 4 6 8 10 12 14 16]
UpperContourLevel: [2 4 6 8 10 12 14 16 18]
Отобразите Гавайи на карте, добавьте маркер, а затем отобразите многоугольники, возвращенные geocontourxy на карте.
figure usamap([18.5 22.5],[-161 -154]) hawaii = shaperead('usastatehi', 'UseGeoCoords', true,... 'Selector',{@(name) strcmpi(name,'Hawaii'), 'Name'}); geoshow(hawaii) geoshow(lat0,lon0,'DisplayType','point','Marker','o',... 'MarkerEdgeColor','k','MarkerFaceColor','k','MarkerSize',10) cmap = parula(1 + length(levels)); for k = 1:length(contourPolygons) lat = contourPolygons(k).Latitude; lon = contourPolygons(k).Longitude; geoshow(lat,lon,'Display','polygon', ... 'FaceColor',cmap(k,:),'FaceAlpha',0.5,'EdgeColor','none') end geoshow(contourLines.Latitude,contourLines.Longitude,'Color','black')
