projinv

Обратная проекция карты с помощью PROJ.4 сопоставьте библиотеку проекции

Описание

пример

[lat,lon] = projinv(proj,x,y) возвращает значения широты и долготы от обратной трансформации проекции. proj структура, задающая проекцию карты. proj может быть проекция карты mstruct или GeoTIFF info структура. x и y массивы координаты карты x-y. Для полного списка GeoTIFF info и сопоставьте структуры проекции, которые можно использовать с projinv, смотрите страницу с описанием для projlist.

Примеры

свернуть все

Импортируйте Бостонские дороги из файла форм и получите структуру проекции из '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 изображение и файл привязки.

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, все права защищены.

Входные параметры

свернуть все

Сопоставьте проекцию, заданную как скалярная структура проекции карты (mstruct) или структура информации о GeoTIFF.

Типы данных: struct

x-координаты одной или нескольких точек в спроектированной системе координат, заданной как скалярное значение, вектор, матрица или массив N-D. Размер должен совпадать с размером y входной параметр.

Типы данных: single | double

y-координаты одной или нескольких точек в спроектированной системе координат, заданной как скалярное значение, вектор, матрица или массив N-D. Размер должен совпадать с размером x входной параметр.

Типы данных: single | double

Выходные аргументы

свернуть все

Геодезические широты одной или нескольких точек, возвращенных как скалярное значение, вектор, матрица или массив N-D, в модулях степеней.

Долготы одной или нескольких точек, возвращенных как скалярное значение, вектор, матрица или массив N-D, в модулях степеней.

Смотрите также

| | | |

Представлено до R2006a

Для просмотра документации необходимо авторизоваться на сайте