遥感数字图像处理教程matlab,遥感数字图像处理实习报告含Matlab处理代码

遥感数字图像处理实习报告含Matlab处理代码

(24页)

本资源提供全文预览，点击全文预览即可全文预览,如果喜欢文档就下载吧，查找使用更方便哦！

19.9 积分

辽宁工程技术大学《数字图像处理》上机实习报告教学单位 辽宁工程技术大学专 业 实习名称 班 级 测绘研11-3班 学生姓名 路聚峰 孙华牛学 号 471120212 指导教师实习1读取BIP、BILS BSQ文件一、 实验目的用Matlab读取BIP、BIL、BSQ文件，并将结果显示出來。遥感图像包括多个波段，有多种存储格式，但基本的通用格式有3种，即 BSQ、BIL和BIP格式。通过这三种格式，遥感图像处理系统可以对不同传感器 获取的图像数据进行转换。BSQ是像素按波段顺序依次排列的数据格式。BIL 格式中，像素先以行为单位块，在每个块内，按照波段顺序排列像素。BIP格式 屮，以像素为核心，像素的各个波段数据保存在一起，打破了像素空间位置的连 续性。用Matlab读取各个格式的遥感数据，是图像处理的前提条件，只有将图像 读入Matlab I作空间，才能进行后续的图像处理工作。二、 算法描述1. 调用fopen函数用指定的方式打开文件。2. 在for循环中调用fread函数，用指定的格式读取各个像素。3. 用reshape函数，重置图像的行数列数。4. 用imadjust函数调整像素的范围，使其有一定对比度。5. 用imshow显示读取的图像。三、 Matlab源代码1 .读取BSQ的源代码：clear allclclines=400;samples=640;N=6;img=fopen(1D:\sample_BSQ1,1rb1);for i=l:Nbi=fread(imgzlines*samplesz 1uint81);band_cov=reshape(biA samples,lines);band_cov2=band_cov1; band_uint8=uint8(band_cov2);tif=imadjust(band_uint8);mkdir('D:\MATLAB*,'tifbandsl')name=['D:\MATLAB\tifbandsl\tif', int2str(i) A .tif'];imwrite(tifA nameA 1tif1);tilt=['波段',int2str (i)];subplot(3,2,i),imshow(tif);title (tilt);endfclose(img);2 •读取BIP源代码clear allclclines=400;samples=640;N=6;for i=l:Nimg=fopen('D:\MATLAB\sample_BIP','rb');bO=fread(imgA i-lA 1uint81);b=fread(img,lines*samplesz 1uint81, (N-l));band_cov=reshape(b,samples,lines);band_cov2=band_cov1;% x adAband_uint8=uint8(band_cov2);tif=imadjust(band_uint8);mkdir(1E:\MATLAB1z 1tifbands1)name=[* E:\MATLAB\tifbands\tif *,int2str (i), *.tif *];imwrite(tif,name,1tif1); %imwrite(A,filename,-省略部分-r (2fcenter (3,center (4,center (5,center (6,)=crop./n_crop;)=forest./n_forest;) = water./n_water;)=soill./n_soill;) = soil2•/n_soil2;)=soil3./n soil3;max_d = max(max(abs(center-ini_center)));ini center = center;img2 = reshape(classes,y,xz 3);RGB_classes = cat(3fimg2(:,:,1) 1,img2(:f :f 2) 9fimg2(:f :f 3) 9); imshow(uint8(RGB_classes))A title(1kmeans-OAa^a1u1);四、运行结果图1:最小距离分类实习7实习8图2：最大似然分类图3： K-means分类大气校正.反射率、地表亮温反演Habib教授摄影测量课程总结主要讲述了摄影测量里的共线方程。由共线方程求坐标的正算与反算。即由 内方位元素、外方位元素、像点坐标去求解地而点坐标，及由内方位元素、地而 点坐标以及相应的像点坐标去求外方位元素。并考虑到影像扭曲变形的影响，采 用一定的算法消除。兀=兀『『(X-X°)+»d)+ r3「(Z-Z。) (1)• -T 巾・(X -XJ +广23 ・d)+ D3 YZ-z。)=_c "2 YX-XJ + f m +缶竝-乙丿 (2)P 加•(X-XJ +厂23心-匕)"33迄-乙)图1:共线方程示意图(Xo,Y(),Z())投影中心在大地坐标系中的坐标，与地面点同一坐标系，为旋转矩阵，即外方位元素。a(x,y)为像点坐标。Habib教授还主要讲述了两个主要的Matlab程序。第一个是由内外方位元素以 及像点坐标去求地面点坐标，其中考虑到了影像扭曲变形的影响，并采用一定的 算法消除。第二个是由内方位元素，地面控制点坐标，以及近似外方位元素，通 过最小二乘迭代法去求外方位元素。求地方点坐标的代码主要实现步骤：1. 首先，由外方位元素求解旋转矩阵，定义旋转矩阵函数function |R]=Rotation(W, P, K)。2. 由扭曲形变改正算法，定义一个形变函数，用来求解形变量，function [dist.x, dist_y] = Estimate_Distortion(x, y, xp, yp, kl, k2, k3, pl, p2, p3, Al, A2)。3. 将内方位元素，形变改正系数，地面点坐标的文件，加载到工作空间。4. 通过旋转矩阵，将经过形变改正的左右影像的像点坐标转换到地面坐标系 下，使其与地面点坐标在同一坐标系下。得到VL、VR两个向量。5. 由左右影像的外方位元素(XO_L,YO_L,ZO_L)和(XO_R,YO_R,ZO_R)作 差求出基线向量(BX,BY,BZ) oVL(1)6.「x(i)[皿⑵-V/?(2)X⑵.—皿⑶-V/?(3)BX~BY ,通过解方程，求的左右向量的比例系数BZX(l)X⑵"这样便得地面点坐标：XG_L = XO_L + X(1)* VL(1); YG_L = YO_L + X(1)* VL(2);ZG_L = ZO_L + X( 1 严 VL(3); 关 键 词： matlab 图像处理 实习 报告 数字 处理 代码 遥感

 天天文库所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。

关于本文