Передайте проекцию карты, пользующуюся библиотекой проекции карты PROJ.4
[x,y] =
projfwd(proj,lat,lon)[ возвращает 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