资源描述：
《数字地图处理实验步骤与成果》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、下面是基于信阳市固始县.jpg遥感图像提取道路处理综合实验的思路和步骤。实验步骤：1从GoogleEarth软件上截下我的家乡信阳市固始县的遥感图像，并保存为jpg格式（为了便于实验处理）。代码如下：clear,closeallPicture=imread('信阳市固始县.jpg');figure,imshow(Picture);截的图如下：由图可知该地区有道路，湖泊，居民地等复杂地物，因此给实验带了很大困难。本次实验只是用到所学的知识提取出道路的大概面貌。2要在MatLab软件上进行遥感实验数据处理，JPG格式受到了很
2、大的限制，因此下一步必须进行数据格式的转换，由真彩图转换为灰度图像。代码如下：A1=rgb2gray(Picture);figure,imshow(A1);得到的图像如下：图像的颜色有所变化，重要的是数据的编码有很大不同。3查看该灰度图像的直方图，以便下面进行灰度变换，代码如下：imhist(A1);得到的直方图如下：下面进行直方图均衡化处理，代码如下：A2=1*A1+150/255;%增加图像的对比度，便于提取地物；imsubtract(I,50)%图像相减imsubtract（一张求减可以调节灰度值、两张相减可以消除
3、一些地物）对原图像进行直方图均衡化，代码如下：A31=imadjust(adapthisteq(A2,'NumTiles',[1010]));A32=imadjust(A2,stretchlim(A2),[]);4下面把灰度图像转化成二值图像进行处理，根据产生的图像一个一个值手动进行试探，直到最后产生了最满意的图像，代码如下：A41=im2bw(A31,180/255);figure,imshow(A41);%转换前的二值图像A42=im2bw(A32,100/255);figure,imshow(A42);%变换后的二
4、值图像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');figu
5、re,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);figur
6、e,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),
7、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=imero
8、de(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('矩阵腐蚀'
