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

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

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


下一篇:MATLAB中文件的读写和数据的导入导出