《遥感院matlab报告》由会员分享，可在线阅读，更多相关《遥感院matlab报告（27页珍藏版）》请在金锄头文库上搜索。
1、为了适应公司新战略的发展，保障停车场安保新项目的正常、顺利开展，特制定安保从业人员的业务技能及个人素质的培训计划遥感院matlab报告辽宁工程技术大学数字图像处理上机实习报告教学单位辽宁工程技术大学专业摄影测量与遥感实习名称遥感数字图像处理班级测绘研11-3班学生姓名路聚峰学号指导教师孙华生实习1读取BIP、BIL、BSQ文件一、实验目的用Matlab读取BIP、BIL、BSQ文件，并将结果显示出来。遥感图像包括多个波段，有多种存储格式，但基本的通用格式有3种，即BSQ、BIL和BIP格式。通过这三种格式，遥感图像处理系统可以对不同传感器获取的图像数据进行转换。BSQ是像素按波段顺序依次排列的
2、数据格式。BIL格式中，像素先以行为单位块，在每个块内，按照波段顺序排列像素。BIP格式中，以像素为核心，像素的各个波段数据保存在一起，打破了像素空间位置的连续性。用Matlab读取各个格式的遥感数据，是图像处理的前提条件，只有将图像读入Matlab工作空间，才能进行后续的图像处理工作。二、算法描述1.调用fopen函数用指定的方式打开文件。2.在for循环中调用fread函数，用指定的格式读取各个像素。3.用reshape函数，重置图像的行数列数。4.用imadjust函数调整像素的范围，使其有一定对比度。5.用imshow显示读取的图像。三、Matlab源代码1.读取BSQ的源代码：cle
3、arallclclines=400;samples=640;N=6;img=fopen(D:sample_BSQ,rb);fori=1:Nbi=fread(img,lines*samples,uint8);band_cov=reshape(bi,samples,lines);band_cov2=band_cov;band_uint8=uint8(band_cov2);tif=imadjust(band_uint8);mkdir(D:MATLAB,tifbands1)name=D:MATLABtifbands1tif,int2str(i),.tif;imwrite(tif,name,tif);t
4、ilt=波段,int2str(i);subplot(3,2,i),imshow(tif);title(tilt);endfclose(img);2.读取BIP源代码clearallclclines=400;samples=640;N=6;fori=1:Nimg=fopen(D:MATLABsample_BIP,rb);b0=fread(img,i-1,uint8);b=fread(img,lines*samples,uint8,(N-1);band_cov=reshape(b,samples,lines);band_cov2=band_cov;%a?band_uint8=uint8(band_
5、cov2);tif=imadjust(band_uint8);mkdir(E:MATLAB,tifbands)name=E:MATLABtifbandstif,int2str(i),.tif;imwrite(tif,name,tif);%imwrite(A,filename,fmt)tilt=波段,int2str(i);subplot(3,2,i),imshow(tif);title(tilt);fclose(img);end3.读取BIL的源代码clearallclclines=400;samples=640;N=6;fori=1:Nbi=zeros(lines,samples);forj=
6、1:samplesimg=fopen(D:MATLABsample_BIL,rb);bb=fread(img,(i-1)*640,uint8);b0=fread(img,1*(j-1),uint8);bandi_linej=fread(img,lines,uint8,1*(N*samples-1);fclose(img);bi(:,j)=bandi_linej;endband_uint8=uint8(bi);tif=imadjust(band_uint8);mkdir(D:MATLAB,tifbands)name=D:MATLABtifbandstif,int2str(i),.tif;imwr
7、ite(tif,name,tif);tilt=2?,int2str(i);subplot(3,2,i),imshow(tif);title(tilt);end四、运行结果图1：读取文件的六个波段图实习2均值/中值滤波、边缘信息提取一、实验目的与原理各种图像滤波算子可以实现图像的增强，去噪，边缘提取等。图像增强的目的在于：1.采用一系列技术改善图像的视觉效果，提高图像的清晰度，2.将图像转换成一种更适合于人或机器进行分析处理的形式。它不是以图像保真度为原则，而是通过处理，设法有选择地突出便于人或机器分析某些感兴趣的信息，抑制一些无用的信息，以提高图像的使用价值。图像增强方法从增强的作用域出发，可
8、分为空间域增强和频率域增强。空间域增强就是直接对图像像素灰度进行操作；频率域增强是对图像经傅里叶变换后的频谱成分进行操作，然后经傅里叶逆变换获得所需结果。图像滤波可分为空间域滤波和频率域滤波，前者通过窗口或卷积核进行，它参照相邻像素改变单个像素的灰度值。后者对图像进行傅立叶变换，然后对频谱进行滤波。空间域图像滤波称为平滑和锐化，强调像素与其周围相邻像素的关系。去噪滤波为平滑滤波包括均值滤波和中值滤波。锐化滤波包括罗伯特梯度、索伯尔梯度、拉普拉斯算法、定向检测，用以提取线状地物和边缘。此实验用Matlab采用各种滤波对图像进行了处理，处理结果如下：二、算法描述1.用imread读取图像文件，并用
9、size获取图像的大小。2.设计各种滤波算子。3利用卷积公式对图像的没一个像素进行处理，得到滤波后的图像。4.用imshow显示滤波后的图像。三、Matlab源代码1.均值滤波源码：clearallclcimg=imread();row,column,band=size(img);img0=double(img);f11=1/9;f12=1/9;f13=1/9;f21=1/9;f22=1/9;f23=1/9;f31=1/9;f32=1/9;f33=1/9;img1=img0(:,1,:),img0(:,:,:),img0(:,column,:);img2=img1(1,:,:);img1(:,
10、:,:);img1(row,:,:);filtered=zeros(row,column,band);forii=1:rowforjj=1:columnfiltered(ii,jj,:)=f11*img2(ii,jj,:)+f12*img2(ii,jj+1,:)+f13*img2(ii,jj+2,:)+.f21*img2(ii+1,jj,:)+f22*img2(ii+1,jj+1,:)+f23*img2(ii+1,jj+2,:)+.f31*img2(ii+2,jj,:)+f32*img2(ii+2,jj+1,:)+f33*img2(ii+2,jj+2,:);endendfiltered1=ui
11、nt8(filtered);subplot(1,2,1),imshow(img);title(图1原始RGB图像);subplot(1,2,2),imshow(filtered1);title(图2均值滤波后的图像);imwrite(filtered1,flower_filtered_);2.边缘提取滤波源代码clearallimg=imread();row,column,band=size(img);img0=double(img);f11=1;f12=0;f13=-1;f21=1;f22=0;f23=-1;f31=1;f32=0;f33=-1;img1=img0(:,1,:),img0(:
12、,:,:),img0(:,column,:);img2=img1(1,:,:);img1(:,:,:);img1(row,:,:);filtered=zeros(row,column,band);forii=1:rowforjj=1:columnfiltered(ii,jj,:)=f11*img2(ii,jj,:)+f12*img2(ii,jj+1,:)+f13*img2(ii,jj+2,:)+.第一讲：Matlab与遥感图像处理什么是Matlab？特点是什么？1）MATLAB是MATrixLABoratory的缩写，是一款由美国MathWorks公司推出的商业数学软件。是一种用于算法开发、数
13、据可视化、数据分析以及数值计算的高级语言和交互式环境。2）除了矩阵运算、绘制函数/数据、图像等常用功能外，还可以创建GUI以及对C、C+、Java和FORTRAN等语言的支持。3）主要包括MATLAB主程序、Simulink、Toolbox三大部分。特点：1）友好的编程环境2）自然的程序语言3）强大的科学计算能力4）出色的可视化功能5）丰富的工具箱6）全方位的帮助系统操作界面：命令窗口、工作空间窗口、当前目录窗口、历史命令窗口、一个或多个图形窗口、浏览器打开Help：View菜单中的Help选项、启动平台的树状列表中选择Help项、在Help菜单中选择MatlabHelp菜单项、打开help路
14、径查看html形式的help、在命令界面输入doc图像的数字化：将连续色调的模拟图像经采样量化后转换成数字图像的过程。从模拟图像采样量化编码到数字图像模拟图像：对现实世界的连续表达。数字化意味着数字图像是现实世界的近似表达。数字图像：以有限的二维数值点形式表达的图像。基本单元为像元。*通常，每个像元由一个灰度级来表达，其含义可以是颜色、高度、透明度等。*通常的图像格式包括：每点1样、每点3样、每点4样数字图像处理：【DigitalImageProcessing】又称计算机图像处理，为达到某种目的，利用计算机对数字图像进行的处理过程。DIP的两个主要目的：目视解译、机器感知数字图像处理的三个层次：低级(狭义的图像处理)、中级(图像处理)、高级(图像理解)位深：是指图片中一个像素点占多少位，常见的有8位、24位、32位，就是说分别用这么大的空间来表示一个点的颜色。位深度越大则图
