Интерактивная установка параметров UTM
Самый легкий способ использовать проекцию UTM - через графический пользовательский интерфейс. Вы можете создать или изменить интересующую область UTM с помощью axesmui проекционная панель управления, и получите дополнительную помощь в utmzoneui панель управления.
Можно Shift + кликнуть в окне осей карты или ввести
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.
The
utmzoneuiпанель закрывается, иzoneустановлено значение выбранного поля. Пределы карты обновляются соответственно, и параметры геоида автоматически устанавливаются на соответствующее определение эллипсоида для этой зоны. Можно переопределить выбор по умолчанию, выбрав другой эллипсоид из списка или введя параметры в поле Geoid.
Нажмите кнопку Apply, чтобы закрыть панель управления проекцией.
Затем проекция будет готова для расчетов проекции или команд отображения карты.
Теперь посмотрим базовую карту choropleth из
usstatehishapefile для области в только что выбранной зоне: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 -90mstruct = gcm; [x,y] = projfwd(mstruct, latlim, lonlim)
x =
1.0e+06 *
-1.5029 -0.7829
y =
1.0e+06 *
3.7403 4.5369