Создайте оси карты для Соединенных Штатов Америки
usamap stateusamap(state)usamap 'conus'usamap('conus')usamapusamap(latlim,lonlim)usamap(Z,R)h = usamap(___)h = usamap('all')usamap 'conus' и
usamap('conus') создайте пустую карту оси для совпадающих 48 состояний (то есть, всех состояний, исключая Аляску и Гавайи).
usamap без подарков аргументов меню, из которого можно выбрать одно состояние, округ Колумбия, совпадающие 48 состояний или все состояния.
h = usamap(___)
h = usamap('all')h.
usamap('allequal') совпадает с usamap('all'), но использование 'allequal' будет удалено в будущем релизе.
Сделайте карту Алабамы только.
figure usamap('Alabama') alabamahi = shaperead('usastatehi', 'UseGeoCoords', true,... 'Selector',{@(name) strcmpi(name,'Alabama'), 'Name'}); geoshow(alabamahi, 'FaceColor', [0.3 1.0, 0.675])
Добавьте текст, чтобы маркировать состояние.
textm(alabamahi.LabelLat, alabamahi.LabelLon, alabamahi.Name,... 'HorizontalAlignment', 'center')
Сделайте карту непрерывного континентального массива, который содержит Калифорнию и Монтану.
figure
ax = usamap({'CA','MT'});
set(ax, 'Visible', 'off')
latlim = getm(ax, 'MapLatLimit');
lonlim = getm(ax, 'MapLonLimit');
states = shaperead('usastatehi',...
'UseGeoCoords', true, 'BoundingBox', [lonlim', latlim']);
geoshow(ax, states, 'FaceColor', [0.5 0.5 1])Добавьте метки в каждое состояние.
lat = [states.LabelLat];
lon = [states.LabelLon];
tf = ingeoquad(lat, lon, latlim, lonlim);
textm(lat(tf), lon(tf), {states(tf).Name}, ...
'HorizontalAlignment', 'center')
Сопоставьте совпадающие Соединенные Штаты с различным цветом заливки для каждого состояния.
figure ax = usamap('conus'); states = shaperead('usastatelo', 'UseGeoCoords', true,... 'Selector',... {@(name) ~any(strcmp(name,{'Alaska','Hawaii'})), 'Name'}); faceColors = makesymbolspec('Polygon',... {'INDEX', [1 numel(states)], 'FaceColor', ... polcmap(numel(states))}); %NOTE - colors are random geoshow(ax, states, 'DisplayType', 'polygon', ... 'SymbolSpec', faceColors)
Установите дополнительные настройки отображения.
framem off; gridm off; mlabel off; plabel off;
Сопоставьте США с отдельными осями для Аляски и Гавайев.
figure ax = usamap('all'); set(ax, 'Visible', 'off') states = shaperead('usastatelo', 'UseGeoCoords', true); names = {states.Name}; indexHawaii = strcmp('Hawaii',names); indexAlaska = strcmp('Alaska',names); indexConus = 1:numel(states); indexConus(indexHawaii|indexAlaska) = []; stateColor = [0.5 1 0.5];
Отобразите эти три области.
geoshow(ax(1), states(indexConus), 'FaceColor', stateColor) geoshow(ax(2), states(indexAlaska), 'FaceColor', stateColor) geoshow(ax(3), states(indexHawaii), 'FaceColor', stateColor)
Скройте кадр.
for k = 1:3 setm(ax(k), 'Frame', 'off', 'Grid', 'off',... 'ParallelLabel', 'off', 'MeridianLabel', 'off') end
Все оси, созданные с usamap, инициализируются со сферической моделью Earth, имеющей радиус 6 371 000 метров.
В некоторых случаях usamap использует tightmap, чтобы настроить пределы по осям, трудные вокруг карты. Если вы изменяете проекцию, или только хотите больше пробела вокруг кадра карты, используйте tightmap снова или axis auto.
axes(h(n)), где n = 1, 2 или 3, делает желаемые оси текущими.
set(h,'Visible','on') делает оси видимыми.
axesscale(h(1)) изменяет размер осей, содержащих Аляску и Гавайи к той же шкале как совпадающие состояния.