数字地图处理实验步骤与成果.docVIP

下面是基于信阳市固始县.jpg遥感图像提取道路处理综合实验的思路和步骤。 由图可知该地区有道路，湖泊，居民地等复杂地物，因此给实验带了很大困难。本次实验只是用到所学的知识提取出道路的大概面貌。 2要在MatLab软件上进行遥感实验数据处理，JPG格式受到了很大的限制，因此下一步必须进行数据格式的转换，由真彩图转换为灰度图像。代码如下： A1 rgb2gray Picture ; figure,imshow A1 ; 得到的图像如下： 图像的颜色有所变化，重要的是数据的编码有很大不同。 3查看该灰度图像的直方图，以便下面进行灰度变换，代码如下： imhist A1 ;得到的直方图如下： 下面进行直方图均衡化处理，代码如下： A2 1*A1+150/255;%增加图像的对比度，便于提取地物； imsubtract I,50 %图像相减imsubtract（一张求减可以调节灰度值、两张相减可以消除一些地物） 对原图像进行直方图均衡化，代码如下： A31 imadjust adapthisteq A2,NumTiles,[10 10] ; A32 imadjust A2,stretchlim A2 ,[] ; 4下面把灰度图像转化成二值图像进行处理，根据产生的图像一个一个值手动进行试探，直到最后产生了最满意的图像，代码如下： A41 im2bw A31,180/255 ; figure,imshow A41 ;%转换前的二值图像 A42 im2bw A32,100/255 ; figure,imshow A42 ;%变换后的二值图像 subplot 1,2,1 ,subimage A41 ,title 变换前的二值图像 ; subplot 1,2,2 ,subimage A42 ,title 变换后的二值图像 ; 5 下面对二值图像进行滤波处理，消除一些噪声，其代码如下： A51 medfilt2 A41,[5,5] ;%5*5中通滤波平滑二值图像 figure,imshow A51 ; H fspecial average,3 ;%3*3均值滤波平滑二值图像 A52 imfilter A41,H,corr,replicate ; figure,imshow A52 ; subplot 1,2,1 ,imshow A52 ,title 5*5中通滤波平滑二值图像 ; subplot 1,2,2 ,imshow A51 ,title 3*3均值滤波平滑二值图像 ; 由上面的实验可以看出效果最好的是5*5中通滤波二值图像，所以一下试验的数据都是采用A52。 6 下面对图像进行樱帽变换，消除背景中那些亮度不一致的背景 7下面对图像进行形态学操作，其代码如下： se9 strel line,4,90 ; B1 imdilate A6,se9 ; figure,imshow B1 ; se8 strel line,3,60 A7 imerode B1,se8 ; figure,imshow A7 ; se7 strel square,2 ; A8 imerode A6,se7 ; figure,imshow A8 ; figure,imshow A9 ; subplot 2,2,1 ,imshow B1 ;title 进行线性90°膨胀 ; subplot 2,2,2 ,imshow A7 ,title 线性60°腐蚀 ; subplot 2,2,3 ,imshow A8 ,title 方阵为2腐蚀 ; subplot 2,2,4 ,imshow A7 ,title 线性0°腐蚀 ; se6 strel line,2,0 ; A9 imerode A7,se6 ; se1 strel diamond,1 ; A10 imerode A8,se1 ; %figure,imshow A10 ; se2 strel rectangle,[2,2] ; A11 imerode A10,se2 ; %figure,imshow A11 ; se3 strel square,3 ; A12 imerode A10,se3 ; %figure,imshow A12 ; subplot 2,2,1 ,imshow A9 ;title 线性0°腐蚀 ; subplot 2,2,2 ,imshow A10 ,title diamond 腐蚀 ; subplot 2,2,3 ,imshow A11 ,title 矩阵腐蚀 ; subplot 2,2,4 ,imshow A12 ,title 方阵3腐蚀 ; se4 strel disk,3 ; A13 imerode A9,se4 ; figure,imshow A13 ;%综合提取后的图像 se5 str