matlab在遥感上应用,matlab在遥感数据的应用上机实验

matlab在遥感数据的应用上机实验

《MATLAB 在遥感数据处 理中的应用》 学年学期: 2012—2013 学年 2 学期 院(系、部): 实验名称: MATLAB 在遥感数据 处理中的应用 学 时: 学 号: 姓 名: 班 级: 实验成绩: 指导教师: 辽宁工程技术大学测绘与地理科学学院 2013 年 3 月一、 实验目的及要求 实验目的：熟悉及掌握 MATLAB 中数字图像处理的基本知识，包括如何 读取图像、获取图像相关信息、图像间转化等；了解图像的算术运算在数字图 像处理中的初步应用，体会图像算术运算处理的过程和处理前后图像的变化； 能够对图像进行空间域和频率域的增强、图像的傅里叶变换处理，巩固所学理 论知识；能够使用 MATLAB 软件对图像进行彩色处理，体会图像彩色处理技 术以及对图像处理的效果。能够使用 MATLAB 软件进行图像的分割，通过实 验体会一些主要的分割算子对图像处理的效果，以及各种因素对分割效果的影 响。能够使用数学形态学运算基本方法，了解腐蚀、膨胀、开运算、闭运算取 得的效果，培养处理实际图像的能力。 实验要求：按《MATLAB 在遥感数据处理中的应用上机指导书 》要求，完 成实验一至实验十，并写出操作过程，以图和表的形式给出操作结果。 二、 实验内容 实验一 MATLAB 数字图像基 本处理 (1)实验语句 I=imread( F:\蓝花瓶.jpg ); whos I imfinfo F: \蓝花瓶.jpg imwrite(I, F: \ 蓝花瓶.bmp ); I1=imread( F: \蓝花瓶.bmp ); f=im2bw(I); subplot(231),imshow(I);title( 蓝花瓶. jpg ); subplot(232),imshow(I1);title( 蓝花瓶. bmp ); subplot(233),imshow(f);title( 蓝花瓶二 值图 ); (2)实验成果 实验二 图像的代数运算 1．图像的加法运算 (1)实验语句 imshow(K) I=imread( E:\山水.jpg ); imshow(I) J=imread( E:\树叶.jpg ); figure,imshow(J) K = imadd(I,J); figure,imshow(K); imshow(I) RGB = imread(‘E:\变亮.jpg’) ；%变亮程序 RGB2 = imadd(RGB,50); RGB3 = imadd(RGB,30); subplot(2,2,1);imshow(RGB); subplot(2,2,2);imshow(RGB2); subplot(2,2,3);imshow(RGB3); (2)实验成果 叠加成果：变亮成果： 2．图像的减法运算 (1)实验语句 A=imread( F:\MATLAB 上机\ 山水.jpg ); background = imopen(A,strel( disk ,15)); Ip = imsubtract(A,background); subplot(121),imshow(A);title( 原图像 ); subplot(122),imshow(Ip);title( 减去后图 像 ); (2)实验成果 3. 图像的乘法运算 (1)实验语句 I=imread( E:\山水.jpg ); J = immultiply(I,1.2); subplot(1,2,1);imshow(I); subplot(1,2,2);imshow(J); (2)实验成果 4．图像的除法运算 (1)实验语句 A=imread( F:\实验成果照片\2.jpg ); I = double(A); J= I * 0.43 + 90; A2 = uint8(J); Ip = imdivide(A, A); Imshow(Ip, []); (2)实验成果 实验三 图像增强—灰度变换 (1)实验语句 f=imread( F:\实验成果照片\蓝花瓶灰度.jpg ); g=imhist(f,256); %显示其直方图 g1=imadjust(f,[0 1],[1 0]); figure,imshow(g1) g2=imadjust(f,[0.5 0.75],[0 1]); figure,imshow(g2) g=imread( F:\ 实验成果照片\ 女孩.jpg ); h=log(1+double(g)); %对输入图像对数映射 变换 h=mat2gray(h); %将矩阵 h 转换为灰度 图片 h=im2uint8(h); %将灰度图转换为 8 位图 figure,imshow(h) (2)实验成果 负片图像： 扩展后图像： 原图像： 实验最终图像： 实验四 图像增强—直方图变换 (1)实验语句 I=imread( F:\实验成果照片\蓝花瓶灰度.jpg ); J=histeq(I); imshow(I); title( 原图像 ); figure,imshow(J); title( 直方图均衡化后的图像 ); figure; subplot(1,2,1) ; imhist(I,64); title( 原图像直方图 ) ; subplot(1,2,2); imhist(J,64) ; title( 均衡变换后的直方图 ) ; (2)实验成果 实验五 图像增强—空域滤波(1)实验语句 I=imread( F:\实验成果照片\ 蓝花瓶灰度.jpg ); J = imnoise(I, salt ave1=fspecial( average ,3); ave2=fspecial( average ,5); K = filter2(ave1,J)/255; L = filter2(ave2,J)/255; M = medfilt2(J,[3 3]); N = medfilt2(J,[4 4]); subplot(2,3,1),imshow(I),title( 原图像蓝花瓶 灰度 ); subplot(2,3,2),imshow(J),title( 添加椒盐噪声 J ); subplot(2,3,2),imshow(J),title( 添加椒盐噪声 I ); subplot(2,3,1),imshow(I),title( 原图像蓝花瓶 灰度 I ); subplot(2,3,2),imshow(J),title( 添加椒盐噪声 J ); subplot(2,3,3),imshow(K),title( 均值滤波 3×3 K ); subplot(2,3,4),imshow(L),title( 均值滤波 5×5 K ); subplot(2,3,5),imshow(M),title( 中值滤波 3×3 模板 M ); (2)实验成果 实验六 图像的傅立叶变换 (1)实验语句 I=imread(‘F:\实验成果照片\蓝花瓶灰度.jpg’);imshow(I);