《遥感数字图像处理》实验报告

《遥感数字图像处理》实验报告

《遥感技术原理与应用》期末报告

研究生《遥感技术原理与应用》

期末考试报告

题目：利用TM遥感数据进行土地覆盖分类和制图

专业：地图学与地理信息系统

2015．12

一、研究方法

缨帽变换：也称K-T变换，是一种特殊的主成分变换。但与主成分不同，其旋转轴不是指向主成分方向，而是指向与地面景物有密切关系的方向，特别是与植物生长过程和土壤有关。传统的NDVI植被信息提取方法受到影像空间分辨率的限制，对影像上信息量少的植被（如道路两旁的行道树、居民小区中的绿地等）提取效果不佳。缨帽变换对区分不同类型植被类型如树、灌木、草地、农作物等非常有效，此次试验具有较好的应用。

支持向量机分类法：是建立在统计学习理论的VC维理论和结构风险最小原理基础上的，根据有限的样本信息在模型的复杂性（即对特定训练样本的学习精度）和学习能力（即无错误地识别任意样本的能力）之间寻求最佳折中，以求获得最好的推广能力。

最大似然分类法：假设每一个波段的每一类统计都呈正态分布，计算给定像元属于某一训练样本的似然度，像元最终被归并到似然度最大的一类当中。

二、研究内容及数据

对富民县散旦乡TM影像进行信息挖掘后突出植被和水体等地物信息；结合二调数据，选择样本，分别用最大似然和支持向量机（SVM）分类法对散旦乡进行分类，通过对比分类精度，比较两种分类方法的优缺点。

数据：对富民县进行裁剪后得到的散旦乡Landsat TM影像；富民县二类调查小班数据；富民县县行政区数据。

三、研究过程

1.裁剪研究区域

将富民县行政区数据导入ArcGIS软件中，根据属性表查找得到散旦乡数据，导入ENVI，再利用ENVI提供的不规则裁剪工具进行裁剪得到散旦乡TM影像（4，3，2假彩色合成），见图1、2。

图1 散旦乡在富民县的位置

图2 研究区原始影像

2.缨帽变换

在主菜单Transforms→Tassled Cap中使用缨帽变换对研究区影像进行正交变换，变换结果包括亮度“Brightness”、绿度“Greenness”、第三波段“Third”三个波段信息。

图3 缨帽变换结果

3.归一化植被指数——NDVI的提取

NDVI(Normal Differential Vegetation Index)，其表达式为NDVI=(NIR-R)/(NIR+R)，是基于近红外波段与红外波段的归一化比值指数。利用波段运算工具Band Math对散旦乡影像进行NDVI运算（影像中3，4波段分别代表公式中的R和NIR），得到植被指数影像。

图4 NDVI提取结果

4.图像合成

对缨帽变换得到的绿度，NDVI得到的植被指数，以及原散旦乡TM影像的4波段进行合成，得到一张新的散旦乡影像，信息挖掘前后对比见图5。通过对比可以看出，植被及水体地物均得到了增强，为接下来影像的分类工作提供了方便。

图5 原始影像（左）与信息挖掘后（右）对比

5.选择训练样本

在ArcGIS中，依据属性表中的class字段，根据分类要求提取6个类别的图层数据；然后导入ENVI，叠加显示在影像上（图上红色范围），然后在小班范围内勾绘的训练样本（图上蓝色区域），参考提取的小数据勾画训练样本，见图6。

图6 选择训练样本

6.影像分类

为避免背景参与分类，使用主菜单下Basic tool→masking→build masking生成掩膜文件，然后对影像分别按照最大似然和支持向量机两种方法进行分类得到分类结果图，如图7所示。

图7 最大似然（左）与支持向量机（右）分类结果图7.分类后处理

分类结果中，不可避免会产生一些面积很小的图斑，需要对这些小图斑进行处理。在主菜单Classification→Post Classification中，选择Majority/Minority工具将小图斑合并

拉曼光谱

拉曼光谱实验报告 一、实验目的 1. 了解拉曼光谱的基本原理、主要部件的功能； 2. 了解拉曼光谱对所观察与分析样品的要求； 3. 了解拉曼光谱所观察材料的微观组织结构和实际应用； 4. 初步掌握制样技术和观察记录方法 二、实验仪器原理 1928年C.V.拉曼实验发现，当光穿过透明介质被分子散射的光发生频率变化，这一现象称为拉曼散射，同年稍后在苏联和法国也被观察到。在透明介质的散射光谱中，频率与入射光频率υ0相同的成分称为瑞利散射；频率对称分布在υ0两侧的谱线或谱带υ0±υ1即为拉曼光谱，其中频率较小的成分υ0－υ1又称为斯托克斯线，频率较大的成分υ0+υ1又称为反斯托克斯线。靠近瑞利散射线两侧的谱线称为小拉曼光谱；远离瑞利线的两侧出现的谱线称为大拉曼光谱。瑞利散射线的强度只有入射光强度的10-3，拉曼光谱强度大约只有瑞利线的10-3。小拉曼光谱与分子的转动能级有关，大拉曼光谱与分子振动-转动能级有关。拉曼光谱的理论解释是，入射光子与分子发生非弹性散射，分子吸收频率为υ0的光子，发射υ0－υ1的光子（即吸收的能量大于释放的能量），同时分子从低能态跃迁到高能态（斯托克斯线）；分子吸收频率为υ0的光子，发射υ0+υ1的光子（即释放的能量大于吸收的能量），同时分子从高能态跃迁到低能态（反斯托克斯线)。分子能级的跃迁仅涉及转动能级，发射的是小拉曼光谱；涉及到振动-转动能级，发射的是大拉曼光谱。与分子红外光谱不同，极性分子和非极性分子都能产生拉曼光谱。激光器的问世，提供了优质高强度单色光，有力推动了拉曼散射的研究及其应用。拉曼光谱的应用范围遍及化学、物理学、生物学和医学等各个领域，对于纯定性分析、高度定量分析和测定分子结构都有很大价值。 拉曼效应起源于分子振动（和点阵振动）与转动，因此从拉曼光谱中可以得到分子振动能级(点阵振动能级)与转动能级结构的知识。用虚的上能级概念可以说明了拉曼效应： 设散射物分子原来处于基电子态，当受到入射光照射时，激发光与此分子的作用引起的极化可以看作为虚的吸收，表述为电子跃迁到虚态（Virtual state），虚能级上的电子立即跃迁到下能级而发光，即为散射光。设仍回到初始的电子态，则有如图所示的三种情况。因而散射光中既有与入射光频率相同的谱线，也有与入射光频率不同的谱线，前者称为瑞利线，后者称为拉曼线。在拉曼线中，又把频率小于入射光频率的谱线称为斯托克斯线，而把频率大于入射光频率的谱线称为反斯托克斯线。

遥感数字图像处理实习1

（1）以多波段组合方式将GeoTIFF格式的白银市TM原始数据转换为ENVI Standard 格式： 利用Basic Tools/Layer Stacking弹出对话框然后Import File，弹出对话框，导入GeoTIFF格式的TM原始数据，选择波段1、2、3、4、5和7， 点击OK，利用Choose选择输出路径及文件名，同时可以利用Reorder Files对输入的文件根据自己的需要进行调换顺序，点击OK输出ENVI Standard格式的数据。 （2）查询并记录影像文件的基本信息、投影信息，以及各个波段直方图信息，然后编辑头文件： 利用Basic Tools/Resize Data弹出对话框里面选择要查看的影像，左 边会出现其基本信息，如图所示：也有投影信息，既可以用来看单波段的也可以看合成后整个影像的信息。在对话框下，合成影像的名字上右击，选择Quick Statistics弹出对话框，在此对话框中点击Select Plo下拉菜单，选择单波段或者多波段的直方图，相应的对话框中会出现直方图（在结果与分析中记录），还可以右击选择edit修改横、纵坐标的单位。 同样的在合成影像的名字上右击，选择Edit Head，弹出对话框

然后点击Edit Attributes/Band Name弹出对话框，选中波段输入修改 后名字，点击OK即可进行波段名字的修改。点击Edit Attributes/Wavelengths弹出进行相应的波长的修改。 （3）在View视窗中，利用影像缩小、放大、漫游工具识别影像中的土地利用/土地覆盖类型: 可以结合当地的google earth上高分辨率的遥感影像，进行识别，利用Viewer视窗下Tools/SPEAR/Google Earth/Jump to Location可以在google earth上显示View主视窗中相应选中地物对应的位置。 （4）利用Viewer视窗打开影像，分别选取4、3、2和7、4、2波段组合进行假彩色合成，观察实习内容中所要求地物的色调变化： 利用File/Open Image File，选择第1步合成的ENVI Standard 格式的数据，弹出对话框，在其中选择RGB Color，将R、G、B分别设为4、3、2波段，点击Load Band，在Viewer#1中出现了4、3、2波段组合的假彩色图像，再在此窗口中，点击Display/New Display,弹出Viewer#2，选择RGB Color,将R、G、B分别设为7、4、2波段，点击Load Band，在Viewer#2中出现了7、4、2波段组合的假彩色图像，在Viewer窗口中右击选择Link Displays，弹出对话框，点击OK，可以把两个窗口中同一位置进行连接起来， 即其中一个窗口放大、缩小、漫游到某个位置，另外一个也跟着漫游到其相对应的位置。这样可以进行地物色调变化的对比。 （5）提取6种地物在不同波段的数值（Digital Number，DN），做光谱剖面图： 在Viewer视窗中Tools/Profile/Z Profile(Spectrum)弹出对话框,在其 Options下拉菜单中勾选Plot Key,对话框中出现了Viewer视窗中选中的目标地物的X，Y坐标，然后勾选Collect Spectra，鼠标箭头变为十字箭头，在目标地物中取九个点（本来图上就有一个，总共是十个点），然后在选择File/Save Plot As/ASCII弹出对话框 ，点击Select All Items，利用Choose选择输出路径和文件名，点击 OK，将其保存为.txt格式。选六种地物，重复以上操作，提取不同波段的数值（Digital Number，DN）。将.txt格式的文件用excel打开，然后用插入函数中的average函数求出每种地物的平均DN值，然后做出光谱剖面（光谱图如结果与分析中所示）。 （6）使用Excel制作6种地物的样本特征光谱统计表： 在Excel中分别使用插入函数中的AVERAGE、VAR、STDEV、MAX和MIN函数求出各地物样本DN值在各个波段的平均值、方差、标准差、最大值和最小值。然后，在07版Excel 的“Microsoft Office 按钮”，单击“Excel 选项”。“加载项”，然后在“管理”框中，选择“Excel 加载项”，单击“转到”弹出“加载宏”，在弹出来的对话框中选择“分析工具库”，并点击确定。然后从“工具”中找到“数据分析”，从“数据分析”对话框中选择“协方差”，并导入某种地物需求协方差的数据区域并选择“逐行”进行，最后选择数据输出区域并确定，则可得该地物的协方差矩阵。同理，在从“数据分析”对话框中选择“相关系数”，进行相应操作，可求得相关系数矩阵。（在结果与分析中附有个地物的样本特征光谱统计表）（7）制作散点图： 在Excel中，打开6种地物的样本DN数据（5步骤产生的），选择band2和band4做散

数字信号处理实验报告

一、实验名称：基本信号的产生 二、实验目的：I 利用MATLAB 产生连续信号并作图 II 利用MATLAB 产生离散序列并作图 III 利用MATLAB 进行噪声处理 三、 实验内容： I 利用MATLAB 产生下列连续信号并作图 ①X(t)=-2u(t-1),-1=0); plot(t,x); 图形如右： ② X(t)=-(e^-0.1t)*sin(2/3*t),0

-1.5-1 -0.5 0.5 1 1.5 2 II 利用MATLAB 产生下列离散序列并作图 ① X(t)=1,-5<=t<=5 else 0,-15<=t<=15 MATLAB 程序如下: k= -15: 15; x=[zeros(1,10),ones(1,11),zeros(1,10)]; stem(k,x) 图形如下： ② X(t)=0.9^k*(cos(0.25*pi*k)+sin(0.25*pi*p),-20

拉曼光谱实验报告

拉曼光谱实验 姓名学号 何婷21530100 李玉环21530092 宋丹21530111 [实验目的] 1、了解Raman光谱的原理和特点； 2、掌握Raman光谱的定性和定量分析方法； 3、了解Raman光谱的谱带指认。 4、了解显微成像Raman光谱。 [仪器和装置] 1、显微Raman光谱系统一套，拉曼光谱仪的型号为SPL-RAMAN-785 USB2000+的拉曼光谱仪，自带785nm激光； 2、带二维步进电机平移台一台（有控制器一台）； 3、PT纳米线样品； 4、光谱仪软件SpectraSuite； 5、步进电机驱动软件； 6、摄像头（已与显微镜集成在一起）。 [实验内容] 1、使用显微Raman系统及海洋光谱软件对单根或多根纳米线进行显微Raman光谱测量， 对测量的图和标准图进行比较，并通过文献阅读对PT纳米线Raman（测量和标准）的谱峰进行指认。 2、使用显微拉曼扫描系统进行二维样品表面拉曼信号收集，并生成样品表面特定波长处的 拉曼信号强度三维图，模拟样品表面拉曼表征。选择多个拉曼波长对样品形状进行观察。[实验结果及分析]

观察PbTiO3的拉曼散射谱并比对具体的拉曼散射光谱数据进行分析，可以找到以上10个拉曼散射峰，分别位于784.54nm，794.94 nm，798.60 nm，802.90 nm，806.84 nm，811.91 nm，817.10 nm，825.29 nm，832.44 nm，879.69nm附近，对应的Raman Shift分别是-7.46 cm-1 159.28 cm-1 216.94 cm-1 284.00 cm-1 344.82 cm-1 422.21 cm-1 500.44 cm-1 621.90 cm-1 725.97 cm-1 1371.21 cm-1。 （通过Raman Shift=1/λ入射-1/λ散射计算得到） PT纳米线Raman测量的谱峰指认： 分析可知，-7.46 cm-1 159.28 cm-1 216.94 cm-1 284.00 cm-1 344.82 cm-1 422.21 cm-1 500.44 cm-1 621.90 cm-1 725.97 cm-1附近的9个振动模，分别对应于PbTiO3的A1（1TO），E（1LO），E（2TO），B1+E，A1（2TO），E（2LO）+A1（2LO），E（3TO）A1（3TO），A1（3LO）声子模。 位于159.28 cm-1附近的模对应PbTiO3纳米线表面的TiO6八面体相对于Pb的振动；位于500.44 cm-1附近的模分别对应于表面Ti-O或Pb-O键的振动；位于725.97 cm-1附近的模对应于TiO6八面体中Ti-O键的振动。而位于284.00 cm-1的振动模为静模。此外，在725.97 cm-1处PbTiO3还具有额外的Raman振动模，可能与该相中含有大量且复杂的晶胞结构有关。据报道，复杂钙钛矿结构中氧八面体的畸变或八面体内B位离子的移动在某种程度上会破坏平移对称性，引起相邻晶胞不再具有相似的局部电场和极化率。 位于-7.46 cm-1处的拉曼峰强度增强，相比标准PbTiO3纳米线，其余拉曼峰强度均减弱。798nm处样品表面拉曼信号三维强度图：

数字图像处理实验报告

数字图像处理实验报告 实验一数字图像基本操作及灰度调整 一、实验目的 1)掌握读、写图像的基本方法。 2)掌握MATLAB语言中图像数据与信息的读取方法。 3)理解图像灰度变换处理在图像增强的作用。 4)掌握绘制灰度直方图的方法，理解灰度直方图的灰度变换及均衡化的方 法。 二、实验内容与要求 1.熟悉MATLAB语言中对图像数据读取，显示等基本函数 特别需要熟悉下列命令：熟悉imread()函数、imwrite()函数、size()函数、Subplot（）函数、Figure（）函数。 1)将MATLAB目录下work文件夹中的forest.tif图像文件读出.用到imread， imfinfo 等文件，观察一下图像数据，了解一下数字图像在MATLAB中的处理就是处理一个矩阵。将这个图像显示出来（用imshow）。尝试修改map颜色矩阵的值，再将图像显示出来，观察图像颜色的变化。 2)将MATLAB目录下work文件夹中的b747.jpg图像文件读出，用rgb2gray() 将其 转化为灰度图像，记为变量B。 2.图像灰度变换处理在图像增强的作用 读入不同情况的图像，请自己编程和调用Matlab函数用常用灰度变换函数对输入图像进行灰度变换，比较相应的处理效果。 3.绘制图像灰度直方图的方法，对图像进行均衡化处理 请自己编程和调用Matlab函数完成如下实验。 1)显示B的图像及灰度直方图，可以发现其灰度值集中在一段区域，用 imadjust函 数将它的灰度值调整到[0，1]之间，并观察调整后的图像与原图像的差别，调整后的灰

度直方图与原灰度直方图的区别。 2) 对B 进行直方图均衡化处理，试比较与源图的异同。 3) 对B 进行如图所示的分段线形变换处理，试比较与直方图均衡化处理的异同。 图1.1 分段线性变换函数 三、实验原理与算法分析 1. 灰度变换 灰度变换是图像增强的一种重要手段，它常用于改变图象的灰度范围及分布，是图象数字化及图象显示的重要工具。 1) 图像反转 灰度级范围为[0, L-1]的图像反转可由下式获得 r L s --=1 2) 对数运算：有时原图的动态范围太大，超出某些显示设备的允许动态范围， 如直接使用原图，则一部分细节可能丢失。解决的方法是对原图进行灰度压缩，如对数变换： s = c log(1 + r )，c 为常数，r ≥ 0 3) 幂次变换： 0,0,≥≥=γγc cr s 4) 对比拉伸：在实际应用中,为了突出图像中感兴趣的研究对象，常常要求 局部扩展拉伸某一范围的灰度值，或对不同范围的灰度值进行不同的拉伸处理，即分段线性拉伸： 其对应的数学表达式为：

遥感数字图像处理

遥感数字图像处理-要点 1．概论 遥感、遥感过程 遥感图像、遥感数字图像、遥感图像的数据量 遥感图像的数字化、采样和量化 通用遥感数据格式（BSQ、BIL、BIP） 遥感图像的模型：多光谱空间 遥感图像的信息内容： 遥感数字图像处理、遥感数字图像处理的内容 遥感图像的获取方式主要有哪几种? 如何估计一幅遥感图像的存储空间大小? 遥感图像的信息内容包括哪几个方面？ 多光谱空间中，像元点的坐标值的含义是什么？ 与通用图像处理技术比较，遥感数字图像处理有何特点？ 遥感数字图像处理包括那几个环节？各环节的处理目的是什么？ 2.遥感图像的统计特征 2.1图像空间的统计量 灰度直方图：概念、类型、性质、应用 最大值、最小值、均值、方差的意义 2.2多光谱空间的统计特征 均值向量、协方差矩阵、相关系数、相关矩阵的概念及意义波段散点图概念及分析 主要遥感图像的统计特征量的意义 两个重要的图像分析工具：直方图、散点图 3.遥感数字图像增强处理 图像增强：概念、方法 空间域增强、频率域增强 3.1辐射增强：概念、实现原理 直方图修正，线性变换、分段线性变换算法原理 直方图均衡化、直方图匹配的应用 3.2空间增强 邻域、邻域运算、模板、模板运算 空间增强的概念 平滑（均值滤波、中值滤波）原理、特点、应用 锐化、边缘增强概念

方向模板、罗伯特算子、索伯尔算子、拉普拉斯算子的算法和特点? 计算图像经过下列操作后，其中心象元的值： – 3×3中值滤波 –采用3×3平滑图像的减平滑边缘增强 –域值为2的3×1平滑模板 – Sobel边缘检测 – Roberts边缘检测 –模板 3.3频率域处理 高频和低频的意义 图像的傅里叶频谱 频率域增强的一般过程 频率域低通滤波 频率域高通滤波 同态滤波的应用 3.4彩色增强 彩色影像的类型：真彩色、假彩色、伪彩色

信号处理实验报告、

第一题 如何用计算机模拟一个随机事件，并估计随机事件发生的概率以计算圆周率π。 解： （一）蒙特卡洛方法可用于近似计算圆周率：让计算机每次随机生成两个0到1之间的数，看以这两个实数为横纵坐标的点是否在单位圆内。生成一系列随机点，统计单位圆内的点数与总点数，（圆面积和外切正方形面积之比为π：4），当随机点取得越多时，其结果越接近于圆周率。 代码： N=100000000; x=rand(N,1); y=rand(N,1); count=0; for i=1:N if (x(i)^2+y(i)^2<=1) count=count+1; end end PI=vpa(4*count/N,10) PI = 3.1420384

蒙特卡洛法实验结果与试验次数相关，试验次数增加，结果更接近理论值 （二）18世纪，法国数学家布丰和勒可莱尔提出的“投针问题”，记载于布丰1777年出版的著作中：“在平面上画有一组间距为d的平行线，将一根长度为l （l

激光拉曼实验报告

激光拉曼及荧光光谱实验 一、实验目的 1、 了解激光拉曼的基本原理和基本知识以及用激光拉曼的方法鉴别物质成分和分子结构的原理； 2、 掌握LRS – II 激光拉曼/荧光光谱仪的系统结构和操作方法； 3、 研究四氯化碳CCL 4、苯C 6H 6等物质典型的振动—转动光谱谱线特征。 二、实验原理 2.1 基本原理 分子有振动。原子分双子的振动按经典力学的观点可以看成是简谐振子，其能量为 A 是振幅，k 是力常数。按照量子力学，简谐振子的能量是量子化的， t=0，1，2，3，···，是振动量子数，f 是振子的固有振动频率。如果在同一电子态中，有振动能级的跃迁，那么产生的光子能量 hf t t E E h )('12-=-=ν 波数为 CO 在红外部分有4.67微米、2.35微米、1.58微米等光谱带，其倒数之比近似为1： 2：3。当Δt=1时，测得的ν ~反映了分子键的强弱。 分子有转动。双原子分子的转动轴是通过质心而垂直于联接二原子核的直线的。按照经典力学，转动的动能是 式中P 是角动量，Ｉ是转动惯量， 222211r m r m I += 可以证明 I P I E 2212 2= =ω2 2 2 121r r m m m m I μ=+= 2222 1212 1 kA kx mv E =+ = 2 12 1m m m m m += hf t E )2 1(+=m k f π21= ,3,2,)(1 ~12ωωωωλ ν =?=-'=-= =t c f t t hc E E

上式中ｒ１，ｒ２和ｒ分别代表两原子到转轴的距离及两原子之间的距离，μ称为约化质量。按照量子力学，角动量应等于 代入上式得 此式可以从量子力学直接推得，Ｊ称为转动量子数。当Ｊ＝０，１，２，３，···等值时，相应的Ｊ（Ｊ＋１）＝０，２，６，１２，···，所以能级的间隔是I h 228π的２，４，６，８，···倍。 实验和理论都证明纯转动能级的跃迁只能在邻近能级之间，就是ΔＪ＝±１。所得 光谱的波长应该有下式表达的值： 谱线波数（ν ~）的间隔是相等的。ＨCL 分子远红外吸收谱中，曾观察到很多条吸收线，这些线的波数间隔应该是２Ｂ，实验测得：Ｂ＝１０．３４厘米 －１ ，所以由此求得 转动惯量Ｉ，进而求得ＨCL 分子中原子之间的核间距这一重要数据。 多原子分子的转动可以近似地看作刚体的转动，这涉及到多个转轴的不同的转动惯量。其谱线结构较为复杂，只有直线型的分子和对称高的分子转动曾研究出一些结果。在分析化学领域中提供了一些分析样品的标准特征谱线可供实验参照。 光通过透明的物体时，有一部分被散射。如果入射光具有线状谱，散射光的光谱中 除有入射光的谱线外，还另有一些较弱的谱线，这些谱线的波数ν '~等于入射光某一波数0~ν加或减一个数值，即10~~~ννν±='。新出现谱线的波数与入射光的波数之差发现与光源无关，只决定于散射物。如果换一个光源，0~ν不同了，但如果散射物不变换，那么0~~νν-'还是等于原来的1~ν，散射光的波数变动反映了散射物的性质。由于散射光的波数等于入射光的波数与另一数值1 ~ν组合的数值，所以这样的散射称作组合散射。 可以在紫外或可见区观测分子的振动和转动能级，通过选择波长在可见光波段的激 ,2,1,0,2) 1(=+=J h J J P π ) 1(82 2+= J J I h E πIc h B J BJ J J J J Ic h hc E E 2''''2'8, ,3,2,12)]1()1([8~1 ππνλ= ==+-+=-==

图像处理实验报告

重庆交通大学 学生实验报告 实验课程名称数字图像处理 开课实验室数学实验室 学院理学院年级信息与计算科学专业 2 班学生姓名李伟凯学号631122020203 开课时间2014 至2015 学年第 1 学期

实验（一）图像处理基础 ?实验目的 学习Matlab软件的图像处理工具箱，掌握常用的一些图像处理命令；通过编程实现几种简单的图像增强算法，加强对图像增强的理解。 ?实验内容 题目A.打开Matlab软件帮助，学习了解Matlab中图像处理工具箱的基本功能；题目B.掌握以下常见图像处理函数的使用： imread( ) imageinfo( ) imwrite( ) imopen( ) imclose( ) imshow( ) impixel( ) imresize( ) imadjust( ) imnoise( ) imrotate( ) im2bw( ) rgb2gray( ) 题目C.编程实现对图像的线性灰度拉伸y = ax + b，函数形式为：imstrech(I, a, b)； 题目D.编程实现对图像进行直方图均衡化处理，并将实验结果与Matab中imhist 命令结果比较。 三、实验结果 1）．基本图像处理函数的使用： I=imread('rice.png'); se = strel('disk',1); I_opened = imopen(I,se); %对边缘进行平滑 subplot(1,2,1), imshow(I), title('原始图像') subplot(1,2,2), imshow(I_opened), title('平滑图像') 原始图像平滑图像

哈尔滨工程大学 语音信号处理实验报告

实 验 报 告 实验课程名称： 语音信号处理实验 姓名： 班级： 20120811 学号： 指导教师 张磊 实验教室 21B#293 实验时间 2015年4月12日 实验成绩 实验序号 实验名称 实验过程 实验结果 实验成绩 实验一 语音信号的端点检测 实验二 语音信号的特征提取 实验三 语音信号的基频提取

实验一 语音信号的端点检测 一、实验目的 1、掌握短时能量的求解方法 2、掌握短时平均过零率的求解方法 3、掌握利用短时平均过零率和短时能量等特征，对输入的语音信号进行端点检测。 二、实验设备 HP 计算机、Matlab 软件 三、实验原理 1、短时能量 语音信号的短时能量分析给出了反应这些幅度变化的一个合适的描述方法。对于信号)}({n x ，短时能量的定义如下： ∑ ∑∞ -∞ =∞ -∞ =*=-= -= m m n n h n x m n h m x m n w m x E )()()()()]()([222 2、短时平均过零率 短时平均过零率是指每帧内信号通过零值的次数。对于连续语音信号，可以 考察其时域波形通过时间轴的情况。对于离散信号，实质上就是信号采样点符号变化的次数。过零率在一定程度上可以反映出频率的信息。短时平均过零率的公式为： ∑∑-+=∞ -∞=--= ---=1)] 1(sgn[)](sgn[2 1 ) ()]1(sgn[)](sgn[21N n n m w w m n m x m x m n w m x m x Z 其中，sgn[.]是符号函数，即 ? ? ?<-≥=0)(10)(1 )](sgn[n x n x n x

激光拉曼光谱实验报告

激光拉曼光谱实验报告 摘要：本实验研究了用半导体激光器泵浦的3Nd + ：4YVO 晶体并倍频后得到的532nm 激 光作为激发光源照射液体样品的4CCL 分子而得到的拉曼光谱，谱线很好地吻合了理论分析的4CCL 分子4种振动模式，且频率的实验值与标准值比误差低于2%。又利用偏振片及半波片获得与入射光偏振方向垂直及平行的出射光，确定了各振动的退偏度，分别为、、、，和标准值0和比较偏大。 关键词:拉曼散射、分子振动、退偏 一， 引言 1928年，印度物理学家拉曼（）和克利希南（）实验发现，当光穿过液体苯时被分子散射的光发生频率变化，这种现象称为拉曼散射。几乎与此同时，苏联物理学家兰斯别而格（）和曼杰尔斯达姆（）也在晶体石英样品中发现了类似现象。在散射光谱中，频率与入射光频率0υ相同的成分称为瑞利散射，频率对称分布在0υ两侧的谱线或谱带01υυ±即为拉曼光谱，其中频率较小的成分01υυ-又称为斯托克斯线，频率较大的成分01υυ+又称为反斯托克斯线。这种新的散射谱线与散射体中分子的震动和转动，或晶格的振动等有关。 拉曼效应是单色光与分子或晶体物质作用时产生的一种非弹性散射现象。拉曼谱线的数目，位移的大小，谱线的长度直接与试样分子振动或转动能级有关。因此，与红外吸收光谱类似，对拉曼光谱的研究，也可以得到有关分子振动或转动的信息。目前拉曼光谱分析技术已广泛应用于物质的鉴定，分子结构的研究谱线特征。 20世纪60年代激光的问世促进了拉曼光谱学的发展。由于激光极高的单色亮度，它很快被用到拉曼光谱中作为激发光源。而且基于新激光技术在拉曼光谱学中的使用，发展了共振拉曼、受激拉曼散射和番斯托克斯拉曼散射等新的实验技术和手段。 拉曼光谱分析技术是以拉曼效应为基础建立起来的分子结构表征技术,其信号来源于分子的振动和转动。它提供快速、简单、可重复、且更重要的是无损伤的定性定量分析，无需样品准备，样品可直接通过光纤探头或者通过玻璃、石英、和光纤测量。拉曼光谱的分析方向有定性分析、结构分析和定量分析。

matlab图像处理实验报告

图像处理实验报告 姓名：陈琼暖 班级：07计科一班 学号：20070810104

目录： 实验一：灰度图像处理 (3) 实验二：灰度图像增强 (5) 实验三：二值图像处理 (8) 实验四：图像变换 (13) 大实验：车牌检测 (15)

实验一：灰度图像处理题目：直方图与灰度均衡 基本要求： （1） BMP灰度图像读取、显示、保存； （2）编程实现得出灰度图像的直方图； （3）实现灰度均衡算法. 实验过程： 1、BMP灰度图像读取、显示、保存； ?图像的读写与显示操作：用imread( )读取图像。 ?图像显示于屏幕：imshow( ) 。 ?

2、编程实现得出灰度图像的直方图； 3、实现灰度均衡算法； ?直方图均衡化可用histeq( )函数实现。 ?imhist(I) 显示直方图。直方图中bin的数目有图像的类型决定。如果I是个灰度图像，imhist将 使用默认值256个bins。如果I是一个二值图像，imhist使用两bins。 实验总结： Matlab 语言是一种简洁，可读性较强的高效率编程软件，通过运用图像处理工具箱中的有关函数,就可以对原图像进行简单的处理。 通过比较灰度原图和经均衡化后的图形可见图像变得清晰，均衡化后的直方图形状比原直方图的形状更理想。

实验二：灰度图像增强 题目：图像平滑与锐化 基本要求： （1）使用邻域平均法实现平滑运算； （2）使用中值滤波实现平滑运算； （3）使用拉普拉斯算子实现锐化运算. 实验过程： 1、 使用邻域平均法实现平滑运算； 步骤：对图像添加噪声，对带噪声的图像数据进行平滑处理； ? 对图像添加噪声 J = imnoise(I,type,parameters)

《遥感数字图像处理》试卷

东南大学2008—2009学年考试试题 课程名称：遥感数字图像处理 学号姓名成绩 一、单项选择题（2分×20=40分） 1.遥感技术是利用地物具有完全不同的电磁波（）或（）辐射特征来判断地物目标和自然现象。 A.反射发射 B.干涉衍射 C.反射干涉 D.反射衍射 2.TM6所采用的10.4~12.6um属于（）波段。 A.红外 B.紫外 C.热红外 D.微波 3.彩红外影像上（）呈现黑色，而（）呈现红色。 A.植被 B. 水体 C.干土 D.建筑物 4.影响地物光谱反射率的变化的主要原因包括（）。 A. 太阳高度角 B.不同的地理位置 C. 卫星高度 D.成像传感器姿态角 5.红外姿态测量仪可以测定（）。 A. 航偏角 B. 俯仰角 C.太阳高度角 D. 滚动角 6.下面遥感卫星影像光谱分辨率最高的是（）。 A. Landsat-7 ETM+ B.SPOT 5 C.IKONOS-2 D. MODIS 7.下面采用近极地轨道的卫星是（）。 A. Landsat-5 B. SPOT 5 C. 神州7号 D. IKONOS-2 8.下面可获取立体影像的遥感卫星是（）。 A. Landsat-7 B.SPOT 5 C.IKONOS-2 D. MODIS 9.侧视雷达图像的几何特征有（）。 A.山体前倾 B.高差产生投影差 C.比例尺变化 D. 可构成立体像对 10.通过推扫式传感器获得的一景遥感影像，在（）属于中心投影。 A.沿轨方向 B. 横轨方向 C. 平行于地球自转轴方向 D. 任意方向 11. SPOT 1-4 卫星上装载的HRV传感器是一种线阵（）扫描仪。 A. 面阵 B. 推扫式 C. 横扫式 D. 框幅式 12.（）只能处理三波段影像与全色影像的融合。 A.IHS变换 B.KL变换 C. 比值变换 D. 乘积变换 13.（）是遥感图像处理软件系统。 A. AreInfo B.ERDAS C. AUTOCAD D. CorelDRAW 14.一阶哈达玛变换相当于将坐标轴旋转了（）。 A.30° B. 45° C. 60° D.90° 15.遥感影像景物的时间特征在图像上以（）表现出来。 A. 波谱反射特性曲线 B.空间几何形态 C. 光谱特征及空间特征的变化 D.偏振特性 16.遥感传感器的分辨率指标包括有（）。 A.几何分辨率 B.光谱分辨率 C.辐射分辨率 D.时间分辨率 17.遥感图像构像方程是指地物点在图像上的（）和其在地物对应点的大地坐标之间的数学关系。 A.投影差 B. 几何特征 C.图像坐标 D. 光谱特征

工程信号处理实验报告

( 2011-2012 学年 第二学期) 重庆理工大学研究生课程论文 课程论文题目： 《工程信号处理实验报告》 课程名称 工程信号处理实验 课程类别 □学位课 非学位课 任课教师 谢明 所在学院 汽车学院 学科专业 机械设计及理念 姓名 李文中 学 号 50110802313 提交日期 2012年4月12日

工程信号处理实验报告 姓名：李文中学号：50110802313 实验报告一 实验名称：数据信号采集及采样参数选定 1实验目的 1.1了解信号采集系统的组成，初步掌握信号采集系统的使用。 1.2加深对采样定理的理解，掌握采样参数的选择方法 1.3了解信号采集在工程信号处理中的实际应用，及注意事项。 2 实验原理 2.1 模数转换及其控制 对模拟信号进行采集，就是将模拟信号转换为数字信号，即模/数（A/D）转换，然后送入计算机或专用设备进行处理。模数转换包括三个步骤：（1）采样，（2）量化，（3）编码。采样，是对已知的模拟信号按一定的间隔抽出一个样本数据。若间隔为一定时间 T，则称这种采样为等时间间隔采样。除特别注明外，一般都采用等时间间隔采样；量化，是一种用有限字长的数字量逼近模拟量的过程。编码，是将已经量化的数字量变为二进制数码，因为数字处理器只能接受有限长的二进制数。模拟信号经过这三步转换后，变成了时间上离散、幅值上量化的数字信号。A/D转换器是完成这三个步骤的主要器件。 在信号采集系统中，A/D 转换器与计算机联合使用完成模数转换。用计算机的时钟或用软件产生等间隔采样脉冲控制 A/D 转换器采样。A/D 转换器通过内部电路进行量化与编码，输出有限长的二进制代码。信号采集系统中，通常由以 A/D转换器为核心的接口电路及控制软件，进行信号采集控制。 *注这部分是由本实验所用的信号采集器自动完成的，以上也是实验器材-信号采集器的部分工作原理。以后实验中就不再赘述。 2.2 信号采集的参数选择

激光拉曼光谱仪实验报告

实验六 激光拉曼光谱仪 【目的要求】 1.学习和了解拉曼散射的基本原理； 2.学习使用激光拉曼光谱仪测量CCL 4的谱线； 【仪器用具】 LRS-3型激光拉曼光谱仪、CCL 4、计算机、打印机 【原 理】 1. 拉曼散射 当平行光投射于气体、液体或透明晶体的样品上，大部分按原来的方向透射 而过，小部分按照不同的角度散射开来，这种现象称为光的散射。散射是光子与物质分子相互碰撞的结果。由于碰撞方式不同，光子和分子之间会有多种散射形式。 ⑴ 弹性碰撞 弹性碰撞是光子和分子之间没有能量交换，只是改变了光子的运动方向，使得散射光的频率与入射光的频率基本相同，频率变化小于3×105HZ ，在光谱上称为瑞利散射。瑞利散射在光谱上给出了一条与入射光的频率相同的很强的散射谱线，就是瑞利线。 ⑵ 非弹性碰撞 光子和分子之间在碰撞时发生了能量交换，这不仅使光子改变了其运动方向，也改变了其能量，使散射光频率与入射光频率不同，这种散射在光谱上称为拉曼散射，强度很弱，大约只有入射线的10-6。 由于散射线的强度很低，所以为了排除入射光的干扰，拉曼散射一般在入射线的垂直方向检测。散射谱线的排列方式是围绕瑞利线而对称的。在拉曼散射中散射光频率小于入射光频率的散射线被称为斯托克斯线；而散射光频率大于入射光频率的散射线被称为反斯托克斯线。斯托克斯线和反斯托克斯线是如何形成的呢？在非弹性碰撞过程中，光子与分子有能量交换, 光子转移一部分能量给分子, 或者从分子中吸收一部分能量，从而使它的频率改变，它取自或给予散射分子的能量只能是分子两定态之间的差值21E E E -=?。在光子与分子发生非弹性碰撞过程中，光子把一部分能量交给分子时，光子则以较小的频率散射出去，称为频率较低的光(即斯托克斯线)，散射分子接受的能量转变成为分子的振动或转动能

数字图像处理实验报告92184

数字图像处理试验报告 实验二：数字图像的空间滤波和频域滤波 姓名：XX学号：2XXXXXXX 实验日期：2017 年4 月26 日1.实验目的 1. 掌握图像滤波的基本定义及目的。 2. 理解空间域滤波的基本原理及方法。 3. 掌握进行图像的空域滤波的方法。 4. 掌握傅立叶变换及逆变换的基本原理方法。 5. 理解频域滤波的基本原理及方法。 6. 掌握进行图像的频域滤波的方法。 2.实验内容与要求 1. 平滑空间滤波： 1) 读出一幅图像，给这幅图像分别加入椒盐噪声和高斯噪声后并与前一张图显示在同 一图像窗口中。 2) 对加入噪声图像选用不同的平滑（低通）模板做运算，对比不同模板所形成的效果， 要求在同一窗口中显示。 3) 使用函数 imfilter 时，分别采用不同的填充方法（或边界选项，如零填 充、’replicate’、’symmetric’、’circular’）进行低通滤波，显示处理后 的图像。 4) 运用for 循环，将加有椒盐噪声的图像进行10 次，20 次均值滤波，查看其特点, 显示均值处理后的图像（提示:利用fspecial 函数的’average’类型生成均值滤 波器）。 5) 对加入椒盐噪声的图像分别采用均值滤波法，和中值滤波法对有噪声的图像做处理， 要求在同一窗口中显示结果。 6) 自己设计平滑空间滤波器，并将其对噪声图像进行处理，显示处理后的图像。 2. 锐化空间滤波 1) 读出一幅图像，采用3×3 的拉普拉斯算子 w = [ 1, 1, 1; 1 – 8 1; 1, 1, 1] 对其进行滤波。 2) 编写函数w = genlaplacian(n)，自动产生任一奇数尺寸n 的拉普拉斯算子，如5 ×5的拉普拉斯算子 w = [ 1 1 1 1 1 1 1 1 1 1 1 1 -24 1 1 1 1 1 1 1 1 1 1 1 1] 3) 分别采用5×5，9×9，15×15和25×25大小的拉普拉斯算子对blurry_moon.tif

信号处理实验报告

数字信号处理 第四次实验报告 一、 实验目的 1.了解离散系统的零极点与系统因果性能和稳定性的关系 2.观察离散系统零极点对系统冲激响应的影响 3.熟悉MATLAB 中进行离散系统零极点分析的常用子函数 4.加深对离散系统的频率响应特性基本概念的理解 5.了解离散系统的零极点与频响特性之间的关系 6.熟悉MATLAB 中进行离散系统分析频响特性的常用子函数，掌握离散系统幅频响应和相频响应的求解方法。 二、实验过程 9.2已知离散时间系统函数分别为 ) 7.05.0)(7.05.0(3 .0)(1j z j z z z H ++-+-= )1)(1(3 .0)() 8.06.0)(8.06.0(3 .0)(32j z j z z z H j z j z z z H ++-+-= ++-+-= 求这些系统的零极点分布图以及系统的冲击响应，并判断系统因果稳定性。 %---------第一式-----------------------------------------------------------------------------% z1=[0.3,0]';p1=[-0.5+0.7j,-0.5-0.7j]';k=1; %z1零点向量矩阵，p1极点向量矩阵,k 系统增益系数---------------------------% [bl,al]=zp2tf(z1,p1,k); %将零极点增益函数转换为系统传递函数 subplot(3,2,1),zplane(bl,al); %zplane 显示离散系统的零极点分布图 ylabel('极点在单位圆内'); subplot(3,2,2),impz(bl,al,20); %impz 绘制系统的冲激响应图 %---------第二式-----------------------------------------------------------------------------% z2=[0,3,0]';p2=[-0.6+0.8j,-0.6-0.8j]'; %z2零点向量矩阵，p2极点向量矩阵---------------------------------------------------% [b2,a2]=zp2tf(z2,p2,k); %将零极点增益函数转换为系统传递函数 subplot(3,2,3),zplane(b2,a2); %zplane 显示离散系统的零极点分布图 ylabel('极点在单位圆上'); subplot(3,2,4),impz(b2,a2,20); %impz 绘制系统的冲激响应图 %---------第三式-----------------------------------------------------------------------------%

拉曼光谱实验报告

成绩 评定 教师 签名 嘉应学院物理学院近代物理实验 实验报告 实验项目：拉曼光谱 实验地点： 班级： 姓名： 座号： 实验时间：年月日

图2 ν? 0ν ν? 斯托克斯线 瑞利线 反斯托克斯线 一、实验目的： 1、 了解拉曼散射的基本原理 2、 学习使用拉曼光谱仪测量物质的谱线，知道简单的谱线分析方法。 二、实验仪器和用具： RBD 型激光拉曼光谱仪 三、实验原理： 按散射光相对于入射光波数的改变情况，可将散射光分为瑞利散射、布利源散射、拉曼散射；其中瑞利散射最强，拉曼散射最弱。在经典理论中，拉曼散射可以看作入射光的电磁波使原子或分子电极化以后所产生的，因为原子和分子都是可以极化的，因而产生瑞利散射，因为极化率又随着分子内部的运动（转动、振动等）而变化，所以产生拉曼散射。 在量子理论中，把拉曼散射看作光量子与分子相碰撞时产生的非弹性碰撞过程。在弹性碰撞过程中，光量子与分子均没有能量交换，于是它的频率保持恒定，这叫瑞利散射，如图（1a ）；在非弹性碰撞过程中光量子与分子有能量交换，从而使它的频率改变，它取自或给予散射分子的能量只能是分子两定态之间的差值12E E E ?=-，当光量子把一部分能量交给分子时，频率较低的光为斯托克斯线，散射分子接受的能量转变成为分子的振动或转动能量，从而处于激发态1E ，如图（1b ），这时的光量子的频率为0ννν'=-?；光量子从较大的频率散射，称为反斯托克斯线，这时的光量子的频率为0ννν'=+?。 最简单的拉曼光谱如图2所示，中央的是瑞利散射线，频率为0ν，强度最强；低频一侧的是斯托克斯线，强度比瑞利线的强度弱很多；高频的一侧是反斯托克斯线，强度比斯托克斯线的 图（1a ） 0h ν ()0h νν+? 0h ν ()0h νν-? 图（1b ） （上能态是虚能态，实 际不存在。这样的跃迁 过程只是一种模型实 际并没有发生） 0h ν 0h ν 0h ν 0h ν

数字图像处理实验报告实验三

中南大学 数字图像处理实验报告实验三数学形态学及其应用

实验三 数学形态学及其应用 一．实验目的 1.了解二值形态学的基本运算 2.掌握基本形态学运算的实现 3.了解形态操作的应用 二．实验基本原理 腐蚀和膨胀是数学形态学最基本的变换，数学形态学的应用几乎覆盖了图像处理的所有领域，给出利用数学形态学对二值图像处理的一些运算。 膨胀就是把连接成分的边界扩大一层的处理。而收缩则是把连接成分的边界点去掉从而缩小一层的处理。 二值形态学 I(x,y), T(i,j)为 0/1图像Θ 腐蚀：[]),(&),(),)((),(0,j i T j y i x I AND y x T I y x E m j i ++=Θ== 膨胀：[]),(&),(),)((),(0 ,j i T j y i x I OR y x T I y x D m j i ++=⊕== 灰度形态学T(i,j)可取10以外的值 腐蚀： []),(),(min ),)((),(1 ,0j i T j y i x I y x T I y x E m j i -++=Θ=-≤≤ 膨胀： []),(),(max ),)((),(1 ,0j i T j y i x I y x T I y x D m j i +++=⊕=-≤≤ 1.腐蚀Erosion: {}x B x B X x ?=Θ: 1B 删两边 2B 删右上 图5-1 剥去一层（皮） 2.膨胀Dilation: {}X B x B X x ↑⊕:＝ 1B 补两边 2B 补左下 图5-2 添上一层（漆） 3.开运算open ：

B B X ⊕Θ=)(X B 4.闭close ：∨ Θ⊕=B B X X B )( 5.HMT(Hit-Miss Transform:击中——击不中变换) 条件严格的模板匹配 ),(21T T T =模板由两部分组成。1T ：物体，2T ：背景。 {} C x x i X T X T X T X ??=?21, 图5-3 击不中变换示意图 性质： （1）φ=2T 时，1T X T X Θ=? （2）)()()(21T X T X T X C Θ?Θ=? C T X T X )()(21Θ?Θ= )/()(21T X T X ΘΘ= 6.细化/粗化 (1)细化（Thin ） C T X X T X XoT )(/??=?= 去掉满足匹配条件的点。 图5-4 细化示意图 系统细化{}n B oB XoB T Xo ))(((21=， i B 是1-i B 旋转的结果（90?，180?，270?）共8种情况 适于细化的结构元素 1111000d d I = d d d L 10110 0= (2)粗化（Thick ） )(T X X T X ??=? 用(){}0,01=T (){}0,12=T 时，X X X T X =?=? X 21 1 1 2 3 T ? XoT X ? X X ?T X ΘT T ⊕