遥感图像及其增强

遥感图像及其增强张超信息与电气工程学院中国农业大学09一月2024主要内容遥感图像特点遥感图像格式遥感图像增强点增强领域增强〔空间域〕低通滤波高通滤波遥感图像与普通图像区别？遥感图像与其它图像区别1、对象不同2、拍摄方式不同〔遥感影像有详细的图像获取参数〕3、采用波段及其范围不同4、遥感图像有3+1种分辨率5、遥感图像经过几何校正后，具有地理坐标系6、遥感影像一般采用非压缩或无损压缩方式存储，特别是定量遥感的图像图像的特征包括图像的自然特征和人工特征自然特征包括：光谱特征：不同类型的地物在各个波段的数字成像，就构成了数字图像的光谱特征。几何特征：包括图像的几何分辨率、图像的纹理结构及图像的几何形状时相特征：在不同时间获得的同一地区各图像之间存在的差异。图像的人工特征即以自然特征为根底，人为定义的某些特性或参数1统计特征〔图像某一波段的最大值、最小值、均值、标准差等〕2直方图特征〔某一波段的直方图分布，见下张幻灯片〕3灰度边缘特征〔反映了图像中目标对象所占的面积大小和形状〕4线和角点特征〔道路和河流等线状地物，角点代表的局部结构关系不因视角而改变，图像几何校正/配准中有用〕5纹理描述特征，图像局部不规那么而在宏观上有规律的特性称为纹理特性以统计等方法计算的各种纹理测度图像统计直方图图像直方图的定义〔1〕一个灰度级别在范围[0，L-1]的数字图像的直方图是一个离散函数p(rk)=nk/nn是数字图像的像素总数nk是数字图像中第k个灰度级的像素总数rk是第k个灰度级，k=0,1,2,…,L-1图像及其直方图例如直方图定义〔2〕一个灰度级别在范围[0，L-1]的数字图像的直方图是一个离散函数p(rk)=nk由于rk的增量是1，直方数可表示为：p(k)=nk即，图像中不同灰度级像素出现的次数两种图像直方图定义的比较p(rk)=nk/n〔第一种定义〕p(rk)=nk〔第二种定义〕第一种定义优点–使函数值正那么化到[0,1]区间，成为实数函数–函数值的范围与数素的总数无关–给出灰度级rk在图像中出现的概率密度的统计不同图像直方图比照直方图用途：评价图像质量图像直方图谷做为灰度阈值可以得到合理的物体的边界直方图用途：图像分割图像边界阈值选择遥感数字图像的存储格式根本图像格式:BMPTIFFJPEG遥感图像存储方式:BSQBIPBIL数字图像存储图像描述信息一般以结构或类来描述主要包括：宽、高等图像数据数组顺序存放的连续数据BMPheadData图像的表示:计算机表示headData头和数据分别在两个文件中ENVI图像格式typedefstructtagBITMAPFILEHEADER{//bmfhWORDbfType; //图像类型BMPDWORDbfSize; //图像大小。

WORDbfReserved1;WORDbfReserved2;DWORDbfOffBits; //数据偏移量。}BITMAPFILEHEADER;BMP图像的描述信息BMP图像文件包括：文件头；信息头；图像数据三局部文件头图像的表示:计算机表示BMP图像的描述信息-信息头typedefstructtagBITMAPINFOHEADER{//bmihDWORDbiSize; //数据结构的尺寸

LONGbiWidth;//位图以象素为单位的宽。

LONGbiHeight;//位图以扫描行为单位的高

WORDbiPlanes;//目标设备的平面数，必须为1WORDbiBitCount//每一个色平面象素的位数1、4DWORDbiCompression;//压缩类型。BI_RGB,BI_REL4DWORDbiSizeImage;//图像以字节为单位的尺寸

LONGbiXPelsPerMeter;//x方向的分辨率dpmLONGbiYPelsPerMeter;//y方向的分辨率dpmDWORDbiClrUsed;//实际使用色彩数

DWORDbiClrImportant;//重要色彩数

}BITMAPINFOHEADERBMP文件中的图像数据图像阵列记录了图像的每一个像素值，在生成图像文件时，Windows从位图的左下角开始，从左到右从下到上，将图像的像素值一一记录下来。图像阵列有压缩和非压缩两种存储格式。根本图像格式-TIFFTIFF是TagImageFileFormat的缩写，是由Aldus公司和Microsoft公司共同开发设计的图像文件格式。TIFF图像文件主要由三局部组成：文件头、标识信息区和图像数据区由于应用了标识的功能，TIFF图像文件才能够实现多幅图像〔多波段〕的存储TIFF文件的特点善于应用指针的功能，可以存储多幅图像文件内数据区没有固定的排列顺序，只规定文件头必须在文件前端，对于标识信息区和图像数据区在文件中可以随意存放可以制定私人用的标识信息除了RGB模式外，TIFF图像文件还能够接受CMYK等多种不同的图像模式可以存储多份调色板数据，调色板的数据类型和排列顺序较为特殊能提供多种不同的压缩数据的方法，便于使用者选择图像数据可以分割成几个局部分别存档JEPGJPEG〔JointPhotographicExpertsGROUP〕是由国际标准组织〔ISO：InternationalStandardizationOrganization〕和国际电报咨询委员会〔CCITTConsultationCommitteeoftheInternationalTelephoneandTelegraph〕为静态图像所建立的第一个国际图像压缩标准MPEG〔连续帧〕JPEG压缩流程符号解码器DCT逆向变换DCT正向变换量化器符号编码器构造8x8的子图合成8x8的子图输入图像NxN压缩图像压缩的图像解压图像颜色空间转换零偏置转换颜色空间转换零偏置转换遥感图像存储格式1、BSQ〔BandSequential〕格式按图像的波段顺序存储，每个波段以像元的行列号排列2、BIL〔BandInterleavedbyline〕格式按波段顺序相间的排列各行像元数据。3、BIP〔BandInterleavedbyPixel〕格式按波段顺序相间的排列各列像元数据B1B2B3B1B2B3B1B2B3BSQBILBIP遥感图像的增强处理遥感图像增强处理的目的改善图像的质量，提高图像目视效果，突出所需要的信息，为进一步遥感目视判读做预处理工作例如：图像比照度不够，图像有噪声、图像边界模糊等具体讲：改善图像的灰度等级，提高图像的比照度，平滑滤波消除噪声，突出地物对象边缘等遥感图像增强分类空间增强点运算线性变换；非线性变换；直方图调整邻域运算：图像平滑钝化；图像锐化频率域增强图像平滑图像锐化

图像增强:点运算增强点运算——灰度级变换增强灰度级变换特点获取变换函数的方法固定函数、交互样点插值灰度级变换的应用亮度调整、比照度拉伸、密度分割灰度级变换的分析点运算〔pointoperation〕特点对于一幅输入图像，经过点运算后，产生的输出图像的灰度值仅由相应输入像素点的灰度值决定，与周围像素不发生直接联系点运算又称为比照度增强、比照度拉伸或灰度变换包括：线性变换，非线性变换〔包括直方图调整〕线性变换线性变换；分段线性变换y=ax+b255021623255023150输入输出输入输出50150255非线性变换指数变换y=beax+c对数变换y=blog(ax+1)+c

图像增强:点运算增强获取变换函数的方法〔1〕固定函数：指数函数、正弦函数、分段线性方程、对数函数，如显示傅立叶的s=clog(1+|r|) 2550255142021623

图像增强:点运算增强获取变换函数的方法〔2〕交互样点插值用过点的三次样条插值曲线，获得变换函数Photoshop软件调整图像

图像增强:点运算增强灰度级变换的应用〔1〕亮度调整——加亮、减暗图像25512825521825512825532加亮减暗

图像增强:点运算增强灰度级变换的应用〔2〕比照度拉伸——提高、降低比照度2554825502551282551422180降低局部提高

图像增强:点运算增强灰度级变换的应用〔2〕提高比照度通常通过直方图得到两个拐点的位置降低比照度降低比照度一般用于输出设备的灰度级小于输入图像的灰度级的情况，如在8位的计算机屏幕上显示11位的卫星影像时。

图像增强:点运算增强灰度级变换的应用〔2〕局部提高、局部降低比照度大于45度提高；小于45度降低255482550255128255142196021623线性非线性

图像增强:点运算增强灰度级变换的应用〔3〕灰度级切片〔密度分割，彩色编码〔伪彩色〕〕255482550134176密度分割densityslice

所谓密度分割将图像的象元按灰度进行归类，归类函数可用以下函数曲线表示：i2550密度分割后伪彩色处理分割后各类图斑用不同颜色表示出来，以示区别。这种处理又称作伪彩色〔颜色查找表〕处理。密度分割的关键是确定像元灰度分割的阈值。直方图调整直方图均衡化直方图规定化〔直方图匹配〕直方图均衡化一种调节图像比照度质量的算法使用的方法是灰度级变换：s=T(r)将原图像的直方图通过变换函数变换为均匀的直方图；然后按均匀直方图修改原图像；最终获得一幅灰度均匀分布的新图像。直方图均衡化示意直方图均衡化步骤1、统计原图像每一灰度级的像元数和累积像元数2、根据如下变换函数式计算每一灰度级xa均衡化后对应的新值，并四舍五入取整，得到新的灰度级xbL-1为变换后灰度级N为像素总个数3、以新值替代原灰度值，形成均衡化后的新图像4、根据原图像像元统计值对应找到新图像像元统计值，做出新直方图直方图规定化(匹配)直方图规定化是指使一幅图像的直方图变成规定形状的直方图，而对图像进行变换的增强方法规定的直方图可以是一幅图像的直方图，通过变换，使两辐图像的亮度变化规律尽可能的接近规定的直方图也可以是指定函数形式的直方图，从而使变换后图像的亮度变化尽可能的服从这种函数分布在遥感图像镶嵌以及融合中经常用到，又称为直方图匹配直方图规定化直方图均衡化是直方图规定化的桥梁

图像增强:点运算增强灰度级变换的分析灰度级变换对层次的影响由于变换是在有限个灰度级上进行的，因此有时候会造成层次的丧失 改变图像原有的真实的光谱响应关系，使图像不能反映不同地类地物反射特性改变图像视觉效果，对于遥感图像计算机自动分类或者遥感定量分析一般不做空域邻域增强空域邻域卷积运算图像平滑去噪图像锐化卷积频域增强与空域模板增强的关系卷积的离散表达式，根本上可以理解为模板运算的数学表达方式 M-1N-1g(x,y)=f*h=f(m,n)h(x–m,y–n)m=0n=0因此，卷积的冲击响应h(x,y)，被称为空域卷积模板，在模板相对原点是对称时成立

图像增强:空域过滤器空域过滤器空域过滤器的根本概念钝化〔平滑〕过滤器锐化过滤器123456789

图像增强:空域过滤空域过滤处理的根本概念空域过滤及过滤器的定义使用空域模板进行的图像处理，被称为空域过滤。模板本身被称为空域过滤器

图像增强:空域过滤空域过滤分类处理效果分类空域过滤器非线性过滤器线性过滤器带通低通高通中值最小值最大值锐化过滤器平滑过滤器

图像增强:空域过滤线性过滤器的定义线性过滤器是线性系统和频域过滤概念在空域的自然延伸。其特征是结果像素值的计算由以下公式定义： R=w1z1+w2z2+…+wnzn其中：wii=1,2,…,n是模板的系数 zii=1,2,…,n是被计算像素及其邻域像素的值123456789

图像增强:空域过滤低通滤波器主要用途：平滑图像、去除噪音高通滤波器主要用途：边缘增强、边缘提取带通滤波器主要用途：删除特定频率、增强中很少用

图像增强:空域过滤非线性过滤器的定义使用模板进行计算，结果值直接取决于像素邻域的值，而不使用如下乘积和的计算:

R=w1z1+w2z2+…+wnzn中值滤波主要用途：平滑图像、去除噪音计算公式：R=mid{zk|k=1,2,…,9}

图像增强:空域过滤最大值滤波，去除胡椒噪声〔暗〕主要用途：寻找最亮点计算公式：R=max{zk|k=1,2,…,9}最小值滤波：去除盐噪声〔亮〕主要用途：寻找最暗点计算公式：R=min{zk|k=1,2,…,9}123456789一般用中值滤波，最大值，最小值很少用

图像增强:空域过滤平滑过滤器平滑过滤器的主要用途低通空域滤波中值滤波

图像增强:空域过滤平滑过滤器的主要用途对图像处理前，删去无用的细小细节连接中断的线段和曲线降低噪音钝化处理，恢复过分锐化的图像图像创艺〔阴影、软边、朦胧效果〕

图像增强:空域过滤过滤器模板系数的设计根据空域中低通冲激响应函数的图形来设计模板的系数例如，选择高斯函数作为冲激函数0

图像增强:空域过滤器设计模板系数的原那么1〕大于02〕模板的和为1，所以下面的模板要除以251111111111111111111111111模板大小一般为奇数3×3；5×5；7×70

空域滤波数字卷积实验以以下图为例，试看卷积实验过程。1011111210101111121111111210102001112111111121211111212121010131112浅灰字卷积模板元素值暗红数字图像灰度值2142121卷积模板〔低通滤波器〕101011121110111112111235231310345834122335231210101313121010131112

滤波结果图

图像增强:空域过滤模板尺寸对过滤器效果的影响模板尺寸越大，图像越模糊，图像细节丧失越多低通空域滤波的缺点和问题图像处理的目的是去除噪音，但是，低通滤波在去除噪音的同时也钝化了边和锋利的细节

图像增强:空域过滤中值滤波中值滤波的原理用模板区域内像素的中值，作为结果值R=mid{zk|k=1,2,…,9}强迫突出的亮点〔暗点〕更象它周围的值，以消除孤立的亮点〔暗点〕中值过滤器是非线性的

图像增强:空域过滤中值滤波算法的实现将模板区域内的像素排序，求出中值。例如：3x3的模板，第5个是中值， 5x5的模板，第13个是中值， 7x7的模板，第25个是中值， 9x9的模板，第41个是中值。对于同值像素，连续排列。如〔10，15，20，20，20，20，20，25，100〕最终值20

图像增强:空域过滤高效的排序算法很重要，如：冒泡法排序中值滤波算法的特点在去除噪音的同时，可以比较好地保存边的锐度和图像的细节

图像增强:空域过滤锐化过滤器锐化过滤器的主要用途高通滤波微分过滤器

图像增强:空域过滤锐化过滤器的主要用途图像识别中，分割前的边缘提取锐化处理模糊的图像图像出图创艺〔突出边界轮廓的图像〕图像的目标识别、提取基于特征的遥感图像融合中，融合前特征增强

图像增强:空域过滤根本高通滤波过滤器模板系数的设计过滤器效果的分析高通空域滤波的缺点和问题算法实现和提高效率

图像增强:空域过滤过滤器模板系数的设计根据空域中高通冲激响应函数的形状来设计模板的系数0Laplace

图像增强:空域过滤器设计模板系数的原那么1〕中心系数为正值，外围为负值2〕系数之和为01-118-11-111-11-111-1-1-1-1-1-1-1-1-1-1-1-1-18-1-1-1-1-1-1-1-1400-100-1

图像增强:空域过滤高通过滤器效果的分析常数或变化平缓的区域，结果为0或很小，图像整体很暗，亮度被降低了在暗的背景上边缘被增强了图像的整体比照度降低了计算时会出现负值，要进行相应处理

图像增强:空域过滤高通空域滤波的缺点和问题高通滤波在增强了边的同时，丧失了图像的层次和亮度〔图像变暗〕

图像增强:空域过滤微分过滤器微分过滤器的原理过滤器扩大因子及模板系数的设计微分过滤器效果的分析

图像增强:空域过滤微分过滤器的原理均值产生钝化的效果，而均值与积分相似，由此而联想到，微分同样能产生相反的效果，即锐化的效果在图像处理中应用微分最常用的方法是计算梯度。函数f(x,y)在(x,y)处的梯度为一个向量：

f=[f/x,f/y]

图像增强:空域过滤微分过滤器的原理计算这个向量的大小为： f=mag(f)=[(f/x)2+(f/y)2]1/2考虑一个3x3的图像区域，z代表灰度级，上式在点z5的值可用数字方式近似。最简单的是用x方向的差〔z5–z6〕和y方向的差〔z5–z8〕，组合为：f[(z5-z6)2+(z5-z8)2]1/2z2z8z5z3z9z6z1z7z4

图像增强:空域过滤微分过滤器的原理用绝对值替换平方和平方根有：

f|z5-z6|+