Самый легкий способ использовать проекцию 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', '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 [x,y] = mfwdtran(latlim, lonlim) x = -1.5029e+006 -7.8288e+005 y = 3.7403e+006 4.5369e+006