Передайте проекцию карты с помощью PROJ.4 сопоставьте библиотеку проекции
[ возвращает x,y] =
projfwd(proj,lat,lon)x и y сопоставьте координаты от прямой трансформации проекции. proj структура, задающая проекцию карты. proj может быть mstruct или GeoTIFF info структура. lat и lon массивы координат широты и долготы.
Для полного списка GeoTIFF info и сопоставьте структуры проекции, которые можно использовать с projinv, смотрите страницу с описанием для projlist.
Отобразите государственную границу Массачусетса по ортофотографии Бостона.
Во-первых, считайте географические векторные данные для государственной границы Массачусетса.
S = shaperead('usastatehi','UseGeoCoords',true, ... 'Selector',{@(name) strcmpi(name,'Massachusetts'),'Name'});
Получите информацию проекции об ортофотографии с помощью geotiffinfo. Затем спроектируйте векторы государственной границы.
proj = geotiffinfo('boston.tif');
[x,y] = projfwd(proj,S.Lat,S.Lon);Считайте и отобразите ортофото изображение Бостона.
[A,R] = readgeoraster('boston.tif'); mapshow(A,R) xlabel('MA Mainland State Plane easting, survey feet') ylabel('MA Mainland State Plane northing, survey feet')
Отобразите государственную границу. Покажите больше детали путем установления пределов карты.
mapshow(x,y,'Color','k','LineWidth',2) xlim([ 645000 895000]); ylim([2865000 3040000]);
geotiffinfo | mfwdtran | minvtran | projinv | projlist