遥感图像处理试验x

1、哈尔滨工业大学遥感图像处理及遥感系统仿真实验报告项目名称：遥感图像处理及遥感系统仿真创新姓名：蒋国韬学号： 1140540124院系：电子与信息工程学院专业：遥感科学与技术指导教师：胡悦时间： 2017年 7月实验一：遥感数字图像的增强、实验目的:利用一幅城市多光谱遥感图像，分析其直方图，并利用对比度增强和去相关拉伸方法对 遥感图像进行增强。二、实验过程：.用multibandread语句读取一幅多光谱遥感图像rs_paris.lan （7波段，512x512图像）的可见1, 2, 3波段（分别对应 R, G, B层）；.显示真彩色图像；3.通过研究直方图（imhist ）,分析直接显示的真彩

2、色图像效果差的原因;4.利用对比度增强方法对真彩色图像进行增强（ imadjust,stretchlim）；5.画出对比度增强后的图像红色波段的直方图；）对图像进行增强;File Edit View Insert Took Desktop Wrriciow Help sd- - 1- 一Sfa- I1 I600040002000口 回6.利用 Decorrelation 去相关拉伸方法（decorrstretch7.显示两种图像增强方法的结果图像。D Figure 5| 口回W三、实验分析：(1)高光谱影像由于含有近百个波段，用 matlab自带的图像读写函数imread和imwrite往往不

3、能直接操作，利用matlab函数库中的multibandred函数，可以读取多波段二进制图像。512X512为像素点，7位波段数，bil为图像数组的保存格式，uint8=uint8为转换到matlab的格式，3 2 1的波段分别对应 RGB三种颜色。(2)直接观察真彩复合图像发现，图像的对比度非常低，色彩不均匀。通过观察红绿蓝三色的波段直方图，可以观察到数据集中到很小的一段可用动态范围内，这是真彩色复合图像显得阴暗的原因之一。另外，根据三种颜色的三维散点图，如下File Edit View Insert Tools Desktop Window Help差- I I a a|血夏/电| 口国|

4、 可知红、绿、蓝三维散点的明显线性趋势显示出可见波段数据的高度相关性，于是未增强的真彩色图像显示的像单色图像。(3)图像经过对比度增强后，进行直方图检测发现数据被扩展到更大范围内的可用动态范 围。地表特征变得较为容易识别。(4)图像经过去相关拉伸后，由下图的三色散点图容易看出，此时波段数据已经没有了高度的相关性。经过线性对比度扩展和去相关，图像的地表特征可识别度得到了很大提高，画面中不同波段的不同被夸大了。File Edit View Insert Tools Desktop Window Help差I- I 口号。研 | Q|A一幻夏/口国 | 口四、实验程序：clear all;clc;P

5、aris_img=multibandread(paris.lan,512,512,7,uint8=uint8,128, bil,ieee-le,Band,Direct,3 2 1);figure(1);imshow(Paris_img);r=Paris_img(:,:,3);g=Paris_img(:,:,2);b=Paris_img(:,:,1);figure(2);subplot(311);imhist(r);title(红色波段直方图);subplot(312);imhist(g);title(绿色波段直方图);subplot(313);imhist(b);title(蓝色波段直方图);

6、figure(3);plot3(r(:),g(:),b(:),.1);Stretch_img=imadjust(Paris_img,stretchlim(Paris_img);figure(4);subplot(211);imshow(Stretch_img);subplot(212);imhist(Paris_img(:,:,3);Decorrstretch_img=decorrstretch(Paris_img,Tol,0.01);figure(5);imshow(Decorrstretch_img);rd=Decorrstretch img(:,:,3);gd=Decorrstretch

7、 img(:,:,2);bd=Decorrstretch_img(:,:,1);figure(6);subplot(311);imhist(rd); subplot(312);imhist(gd); subplot(313);imhist(bd);figure(7);plot3(rd(:),gd(:),bd(:),.);实验二：遥感图像的配准、实验目的:利用一幅未配准的机场图像和一幅同样场景正摄影获得的图像进行配准。二、实验过程:1.读取机场.png和正摄影.png两幅未配准图像（imread ）;Figurt 1要/口因T回ditInsert 工0口 k desktop dow Help2.

8、加载配准点.mat文件，其中包含四组控制配准点;3.用fitgeotrans语句找到配准参数;田 moving Pointsdaiible匚J fixedPoirts v4x2 double4.利用找到的配准参数将未配准的图像进行几何校正( imref2d,imwarp );5.显示正摄影图像为基础的配准后的图像( imshowpair )。164.56351,7O&6S.97722319,222278,92024352.07292 gLi44 5113.2990353,5325130.079B1414046284.59353127983829L5312311,9810353.5325Q Fi

9、gure 2File Edit View Insert Tools D&5 ktop Window Help三、实验分析：(1)寻找配准参数的函数 fitgeotrans ,由配准点得到 movingPoints 和 fixedPoints , movingPoints 是图像上想要移动的点的坐标，至少是两个 double型2维点，fixedPoints为目标点，与 movingPoints规模相同，transformationType为变换类型，包括仿射变换、投影变化、相似变 换等，这个函数主要描述了将movingPoints (设大小为：N*2,N=2 )通过某种变换变化到fixedPoi

10、nts来，最后输出了变换矩阵。得到的 fun_para是一个结构体类型，里面包含了变换 矩阵。(2)要求参加配准的图像满足一定灰度相似性，所以原图像对噪声较敏感，在配准后避免 了直接对灰度图像的操作，所以抗噪性较强。四、实验程序：clear all;airport=imread( 机场.png);figure(1);imshow(airport);right=imread(正摄影.png);figure(2);imshow(right);load配准点.mat;fun_para=fitgeotrans(movingPoints,fixedPoints,projective);fixed_siz

11、e=imref2d(size(right);pic_reg=imwarp(airport,fun_para,fixed_size);figure(3);imshowpair(airport,pic_reg,blend);实验三：遥感图像的分割一、实验目的：利用一幅城市多光谱遥感图像，研究其NDVI值的特点，并利用NDVI阈值处理对图像进行分割，获得植被图像。二、实验过程：1.用multibandread 语句读取一幅多光谱遥感图像 rs_paris.lan (7波段，512x512图像)的可见4, 3, 2波段(分别对应 NIR, R, G层);2.显示假彩色图像;3.显示近红外和红光波段的图

12、像;4.5.6.、NIR -R计算NDVI值：；NIR R选取合适阈值分割植被图像;显示分割后的植被图像。W1 Figure 4File Edit Vieiv Insert Tools Desktop Window Hep才oja j iT %猾奴口 tips、实验分析:（1）近红外（NIR）光谱波段（波段 4）和可见光谱的相关性相对较小。和可见光谱一样，卫星地图得到的不可见光谱的信息也可以通过构建和增强RGB复合图像以后直接观察。近红外NIR波段（波段4）很重要，因为该波段在光谱中对叶绿素反射率较高，利用 NIR波 段可以很好地区分出图像中的植被。当和可见光和绿光（分别是波段3和波段2）进行

13、复合形成一个彩色红外复合图像（CIR）时，这一特点比较有用。彩色红外复合图像通常用来识 别植被或评定植被的生长情况。（2）第一次阈值选取时，没有将原始像素的数据格式转为double,由于近红外和红色波段的相关性存在，导致出现NIR与R某些像素相等的情况，NDVI=（近红外波段辐射率-红外辐射率）/ （近红外辐射率+红外波段辐射率）求出来为 0或者1,得到的分割图像如下:这样分割的效果并不明显，于是将近红外和红色波段的数据格式改为了double,让NDVI求出的数值有小数而不是 0或1,再选取合适的阈值(这里选择了 0.01),得到最后分割后 的植被图像。四、实验程序：clear all;clc

14、;Paris_img=multibandread(paris.lan,512,512,7,uint8=uint8,128,bil,ieee-le,Band,Direct,4 3 2);figure(1);imshow(Paris_img);NIR=multibandread(paris.lan,512,512,7,uint8=uint8,128,bil,ieee-le,Band,Direct,4);R=Paris_img(:,:,3);G=Paris_img(:,:,2);figure(2);subplot(311);imhist(NIR);title(近红外波段直方图);subplot(31

15、2);imhist(R);title(红色波段直方图);subplot(313);imhist(G);title(绿色波段直方图);figure(3);subplot(211);imshow(NIR);title(,近红外波段图像,);subplot(212);imshow(R);title(红色色波段图像);NIR_double=im2double(NIR);R_double=im2double(R);NDVI=(NIR_double-R_double)./(NIR_double+R_double);N=size(NDVI);fori=1:Nfor j=1:Nif NDVI(i,j)0.01

16、;t(i,j)=1;endend end figure(4);imshow(t);实验四：遥感图像特征提取一、实验目的：利用一幅道路遥感图像，提取其颜色特征，并利用颜色特征提取方法对图像进行分割，获得道路提取图像。二、实验过程：1.读取道路图像road.jpg;2.利用makecform和applycfor语句获得图像的lab颜色模型;D 2|。| 鬲.分析道路像素特点；.找到合适阈值并对图像进行分割;File Edit View Insert Tools Desktop Window Help21 d J fe C X * S 口囱口.利用形态学方法进一步提取道路信息;6.得到道路提取图像，

17、并进行显示三、实验分析：（1） LAB模式的特点在于灰度信息和色彩信息的分别保存。Lab模式是由国际照明委员会（CIE）于1976年公布的一种色彩模式。RGB模式是一种发光屏幕的加色模式，CMYK模式是一种颜色反光的印刷减色模式。而Lab模式既不依赖光线，也不依赖于颜料，它是CIE组织确定的一个理论上包括了人眼可以看见的所有色彩的色彩模式。Lab模式弥补了 RGB和CMYK两种色彩模式的不足。Lab模式由三个通道组成，但不是 R、G、B通道。它的一个通道是亮度，即 L。另外两个是色彩通道，用 A和B来表示。A通道包括的颜色是从深 绿色（底亮度值）到灰色（中亮度值）再到亮粉红色（高亮度值）；B通

18、道则是从亮蓝色（底 亮度值）到灰色（中亮度值）再到黄色（高亮度值）。因此，这种色彩混合后将产生明亮的色彩。Lab模式所定义的色彩最多， 且与光线及设备无关并且处理速度与RGB模式同样快，比CMYK模式快很多。因此，可以放心大胆的在图象编辑中使用Lab模式。而且，Lab模式在转换成CMYK模式时色彩没有丢失或被替换。因此，最佳避免色彩损失的方法是：应 用Lab模式编辑图象，再转换为 CMYK模式打印输出。当你将 RGB模式转换成CMYK模式时，Photoshop将自动将RGB模式转换为Lab模式，再转换为 CMYK模式。在表达色彩 范围上，处于第一位的是Lab模式，第二位的是 RGB模式，第三位

19、是 CMYK模式。要注意LAB模式在PHOTOSHOP中很多功能都不能用Lab模式的好处在于它弥补了前面两种色彩模式的不足。RGB在蓝色与绿色之间的过渡色太多，绿色与红色之间的过渡色又太少，CMYK模式在编辑处理图片的过程中损失的色彩则更多，而Lab模式在这些方面都有所补偿。Lab模式与RGB模式相似，色彩的混合将产生更亮的色彩。只有亮度通道的值才影响 色彩的明暗变化。可以将Lab模式看作是两个通道的RGB模式加一个亮度通道的模式。Lab模式是与设备无关的，可以用这一模式编辑处理任何一个图片(包括灰图图片)，并且与RGB模式同样快，比CMYK模式则快女?几倍。Lab模式可以保证在进行色彩模式转

20、换时CMYK范围内的色彩没有损失。如果将RGB模式图片转换成 CMYK模式时，在操作步骤上应加上一个中间步骤，即先转换成Lab模式。在非彩色报纸的排版过程中，应用 Lab模式将图片转换成灰度图是经常用到的。(2)将LAB模式转换为可以进行形态学处理的二值图像需要选取合适的阈值，在此次实验中我选择grayshresh函数求出图像的阈值，graythresh (image)函数输入是一幅图像，在我做的图像的变化检测方向就是输入差异图(两幅图像对应位置做差或者做比，或者做对数比)输出就是阈值。在这个函数中，是使用最大类间方差法(即最优全局阈值法OTSU)找到图片的一个合适的阈值(threshold) o它是按图像的灰度特性，将图像分成背景和目标 2部分。 背景和目标之间的类间方差越大，说明构成图像的2部分的差别越大，当部分目标错分为背景或部分背景错分为目标都会导致2部分差别变小。因此，使类间方差最大的分割意味着错分概率最小。再利用im2bw (将灰度图像转换为二值图像)函数，将找到的阈值输入，就可以 把原图变为一个二值图像。(3) bwareaopen函数的作用：删除二值图像BW中面积小于P的对象，默认情况下使用 8邻域，用来做第一步去噪声。之后分别进行了一次膨胀操作和闭操作。膨胀操作和腐蚀操作是形态学两种基本处理方法之一，腐蚀可以分割(isolate