Передайте проекцию карты с помощью 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'});
Получите структуру проекции для ортофотографии и спроектируйте векторы государственной границы к нему (система Координаты на плоскости состояния Массачусетса, американские Ноги Обзора).
proj = geotiffinfo('boston.tif');
lat = [S.Lat];
lon = [S.Lon];
[x, y] = projfwd(proj, lat, lon); Считайте и отобразите 'boston.tif' ортофото изображение.
[RGB, R, bbox] = geotiffread('boston.tif'); figure mapshow(RGB, R) xlabel('MA Mainland State Plane easting, survey feet') ylabel('MA Mainland State Plane northing, survey feet')
Наложите государственную границу и установите пределы карты, чтобы показать немного больше детали.
hold on mapshow(gca, x, y,'Color','black','LineWidth',2.0) set(gca,'XLim', [ 645000, 895000], ... 'YLIm', [2865000, 3040000]);
boston.tif отображают copyright © GeoEye™, все права защищены.
geotiffinfo | mfwdtran | minvtran | projinv | projlist