mapoutline
Вычислите схему изображения, на которое геоссылаются, или сетки данных
Синтаксис
[x,y] = mapoutline(R,height,width)
[x,y] = mapoutline(R, sizea)
[x,y] = mapoutline(info)
[x,y] = mapoutline(...,'close')
[lon,lat] = mapoutline(R,...)
outline = mapoutline(...)
Описание
[x,y] = mapoutline(R,height,width) вычисляет схему изображения, на которое геоссылаются, или регулярного набора данных с координатной сеткой в координатах карты. R или 3 2 матрица привязки, задающая 2-мерное аффинное преобразование от внутренних пиксельных координат, чтобы сопоставить координаты или MapCellsReference или MapPostingsReference объект. height и width размеры изображения. x и y 4 1 вектор-столбцы, содержащие координаты карты внешних углов угловых пикселей, в следующем порядке:
(1,1), (height,1), (height, width), (1, width).
[x,y] = mapoutline(R, sizea) принимает sizea = [height, width, ...] вместо height и width.
[x,y] = mapoutline(info) принимает скалярный массив структур с полями
| 3 2 матрица привязки |
| Скалярный номер |
| Скалярный номер |
[x,y] = mapoutline(...,'close') возвращает x и y как векторы 5 на 1, добавляя координаты первого из этих четырех углов в конец.
[lon,lat] = mapoutline(R,...), где R пиксели геоссылок к долготе и широте, а не координатам карты, возвращает схему в географических координатах. Долгота должна предшествовать широте в списке выходных аргументов.
outline = mapoutline(...) возвращает угловые координаты в 4 2 или 5 2 массив.
Примеры
Чертите красную схему, формирующую рисунок Бостонского изображения GeoTIFF, на которое ссылаются на Материк Массачусетса система Координаты на плоскости состояния с модулями ног обзора.
figure
info = georasterinfo('boston.tif');
R = info.RasterReference;
[x,y] = mapoutline(R,R.RasterSize,'close');
hold on
plot(x,y,'r')
xlabel('MA Mainland State Plane easting, survey feet')
ylabel('MA Mainland State Plane northing, survey feet')info = imfinfo('concord_ortho_w.tif');
R = worldfileread('concord_ortho_w.tfw','planar', ...
[info.Height info.Width]);
[x,y] = mapoutline(R, info.Height, info.Width, 'close');
x = x * unitsratio('sf','meter');
y = y * unitsratio('sf','meter');
plot(x,y,'k')