在matlab中进行遥感影像地理坐标的相互转换

在matlab中进行图像处理，一般使用的都是图像本地坐标，以左上角（1，1）开始。处理完成后，如果要将结果在带地理坐标的遥感影像中显示，或者需要输出成shp文件，就需要涉及到本地坐标和地理坐标的转换，我们可以使用pix2map函数来实现这个功能。

clc;clear;close all;
[pic, R] = geotiffread('boston.tif'); %读取tif图像
figure(1),imshow(pic);
% pic即为本地图像，下面对本地图像进行操作
gray=rgb2gray(pic);
bw=im2bw(gray,0.1); %阈值分割
bw=1-bw;
bw1 = bwareaopen(bw,500000); %删除面积小于500000的区域
bw2=imfill(bw1,'holes');  %填充孔洞
figure(2),imshow(bw2);
bw3=bwboundaries(bw2);  %得到边缘轮廓
data=bw3{1,1}; %得到轮廓线的坐标，一个N*2的矩阵,此坐标为本地图像坐标
[lonX,latY]=pix2map(R,data(:,1),data(:,2)); %将本地图像坐标转换为地理坐标
figure(3),axis off;
mapshow('boston.tif');  %显示带地理坐标的遥感影像
mapshow(lonX,latY,'Color','r'); %在遥感影像上显示轮廓