Самый легкий способ использовать проекцию UTM через графический интерфейс пользователя. Можно создать или изменить сферу интересов UTM с axesmui
панель управления проекции, и добирается, дополнительная помощь формируют utmzoneui
панель управления.
Вы можете Shift +click в окне осей карты или типе axesmui
отобразить панель управления проекции. Здесь вы запускаете с нуля:
figure
axesm utm
axesmui
Поле Map Projection установлено в cyln: Universal Transverse Mercator (UTM)
.
Примечание
Для UTM и карт UPS, поле Aspect установлено в normal
и не может быть изменен. При попытке задать transverse
, ошибка заканчивается.
Нажмите кнопку Zone, чтобы открыть utmzoneui
панель. Кликните по карте около своей сферы интересов, чтобы выбрать зону:
Обратите внимание на то, что, в то время как можно открыть utmzoneui
панель управления из командной строки, затем необходимо вручную обновить фигуру с именем зоны, которое это возвращает с setm
команда:
setm(gca,'zone',ans)
Нажмите кнопку Accept.
utmzoneui
завершения панели и zone
поле установлено в то, которое вы выбрали. Пределы карты обновляются соответственно, и параметры геоида автоматически устанавливаются на соответствующее определение эллипсоида для той зоны. Можно заменить выбор по умолчанию путем выбора другого эллипсоида из списка или путем ввода параметров в поле Geoid.
Нажмите Apply, чтобы закрыть панель управления проекции.
Проекция затем готова к вычислениям проекции или командам отображения карты.
Представление Now choropleth основывает карту от usstatehi
файл форм для области в зоне, которую вы только выбрали:
states = shaperead('usastatehi.shp', 'UseGeoCoords', true); framem faceColors = makesymbolspec('Polygon',... {'INDEX', [1 numel(states)],... 'FaceColor', polcmap(numel(states))}); geoshow(states,'DisplayType', 'polygon',... 'SymbolSpec', faceColors)
То, что вы видите, зависит от зоны, которую вы выбрали. Предыдущее отображение для зоны 18T, который содержит фрагменты Новой Англии и средние Атлантические состояния.
Можно также вычислить спроектированные декартовы координаты UTM от широт и долгот:
[latlim, lonlim] = utmzone('15S')
latlim = 32 40 lonlim = -96 -90
mstruct = gcm; [x,y] = projfwd(mstruct, latlim, lonlim)
x = 1.0e+06 * -1.5029 -0.7829 y = 1.0e+06 * 3.7403 4.5369