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

下面是基于信阳市固始县.jpg遥感图像提取道路处理综合实验的思路和步骤。实验步骤：1从GoogleEarth 软件上截下我的家乡信阳市固始县的遥感图像，并保存为jpg格式（为了便于实验处理）。代码如下：clear,close allPicture=imread(信阳市固始县.jpg);figure,imshow(Picture);截的图如下：由图可知该地区有道路，湖泊，居民地等复杂地物，因此给实验带了很大困难。本次实验只是用到所学的知识提取出道路的大概面貌。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=strel(line,2,60);B2=imdilate(A13,se5);B3 = bwmorph(A13,skel,Inf);figure,imshow(B3);E=Picture;E(A10)=255;figure,imshow(E);subplot(2,2,1),imshow(A13);title(圆形半径3腐蚀);subplot(2,2,2),imshow(B2),title(线性60腐蚀);subplot(2,2,3),imshow(B3),t