Для отображения областей, простирающихся на несколько зон UTM, используйте проекцию Меркатора в поперечном аспекте. UTM представляет собой зональную систему координат и предназначен для использования подобно серии карт, выбираемой на соответствующем листе. Хотя можно расширить координаты одной зоны на территорию соседней зоны, это обычно не делается. В этом примере показана поперечная проекция Меркатора, подходящая для Чили. В примере обратите внимание на то, как линия нулевого искажения проекции выровнена с преимущественно осью север-юг страны. Конечно, вы не получаете координаты в метрах, которые бы соответствовали координатам проекции UTM, но результаты будут почти такими же точными. Чтобы разместить нулевую линию искажения точно на средней линии страны, используйте лучшие оценки центрального меридиана и угла ориентации вектора ориентации.
Настройте оси карты с поперечным видом и отобразите карту Чили. Вычислите искажение карты.
figure; latlim = [-60 -15]; centralMeridian = -70; width = 20; axesm('mercator',... 'Origin',[0 centralMeridian -90],... 'Flatlimit',[-width/2 width/2],... 'Flonlimit',sort(-latlim),... 'Aspect','transverse'); land = shaperead('landareas.shp', 'UseGeoCoords', true); geoshow([land.Lat], [land.Lon]); framem gridm; setm(gca,'plinefill',1000) tightmap mdistort scale
Вы можете получать предупреждения о баллах от landareas.shp выпадение за пределы допустимой проекционной области. Такие предупреждения можно игнорировать.