Импортируйте Бостонские дороги из файла форм и получите структуру проекции из ортофотографии 'boston.tif'.

roads = shaperead('boston_roads.shp');
proj = geotiffinfo('boston.tif');

Преобразуйте дорожные координаты в единицу длины проекции. В этом примере координаты в файле форм дорог исчисляются в метрах. Единицы длины в спроектированной системе координат являются Ногами Обзора США, зарегистрированными в значение поля UOMLength структуры проекции.

x = [roads.X] * unitsratio('survey feet','meter');
y = [roads.Y] * unitsratio('survey feet','meter');

[roadsLat, roadsLon] = projinv(proj, x, y);

Считайте изображение boston_ovr.jpg и worldfile.

RGB = imread('boston_ovr.jpg');
R = worldfileread(getworldfilename('boston_ovr.jpg'));

Векторы контура состояния чтения для Массачусетса из файла форм usastatehi с помощью селектора, чтобы устранить другие состояния.

S = shaperead('usastatehi', 'UseGeoCoords', true, ...
    'Selector',{@(name) strcmpi(name,'Massachusetts'), 'Name'});

Откройте фигуру с Меркаторской проекцией и отобразите изображение, государственную границу и дороги.

figure
axesm('mercator')

geoshow(RGB, R) 
geoshow(S.Lat, S.Lon, 'Color','red') 
geoshow(roadsLat, roadsLon, 'Color', 'green')

Установите контур карты на северные, западные, и южные пределы изображения и восточный предел государственной границы в ограничительной рамке широты изображений.

[lon, lat] = mapoutline(R, size(RGB(:,:,1)));
ltvals = find((S.Lat>=min(lat(:))) & (S.Lat<=max(lat(:))));
setm(gca,'maplonlimit',[min(lon(:)) max(S.Lon(ltvals))], ...
         'maplatlimit',[min(lat(:)) max(lat(:))])
tightmap

boston_ovr.jpg отображает copyright © GeoEye, все права защищены.