数学模型在遥感图像处理中的应用。
1 绪论。数学模型是近些年发展起来的新学科,是数学理论与实际问题相结合的一门科学。它是关于部分现实世界和为一种特殊目的而作的一个抽象的、简化的结构。
具体来说,数学模型就是为了某种目的,用字母、数学及其它数学符号建立起来的等式或不等式以及图表、图像、框图等描述客观事物的特征及其内在联系的数学结构表达式。 数学模型是针对参照某种事物系统的特征或数量依存关系,采用数学语言,概括地或近似地表述出的一种数学结构,这种数学结构是借助于数学符号刻划出来的某种系统的纯关系结构。
遥感是以航空摄影技术为基础,在20世纪60年代初发展起来的一门新兴技术。它是利用遥感器从空中来探测地面物体性质的,根据不同物体对波谱产生不同响应的原理,识别地面上各类地物,具有遥远感知事物的意思。从广义上说是泛指从远处探测、感知物体或事物的技术。
即不直接接触物体本身,从远处通过仪器(传感器)探测和接收来自目标物体的信息(如电场、磁场、电磁波、**波等信息),经过信息的传输及其处理分析,识别物体的属性及其分布等特征的技术。也就是利用地面上空的飞机、飞船、卫星等飞行物上的遥感器收集地面数据资料,并从中获取信息,经记录、传送、分析和判读来识别地物。当前遥感形成了一个从地面到空中,乃至空间,从信息数据收集、处理到判读分析和应用,对全球进行探测和监测的多层次、多视角、多领域的观测体系,成为获取地球资源与环境信息的重要手段。
遥感信息应用是遥感的最终目的。在遥感信息应用中,最重要的是对遥感图像进行数据处理。遥感影像数字图像处理的内容主要有:
1、图像恢复:即校正在成像、记录、传输或回放过程中引入的数据错误、噪声与畸变。包括辐射校正、几何校正等;
2、数据压缩:以改进传输、存储和处理数据效率;
3、影像增强:突出数据的某些特征,以提高影像目视质量。包括彩色增强、反差增强、边缘增强、密度分割、比值运算、去模糊等;
4、信息提取:从经过增强处理的影像中提取有用的遥感信息。包括采用各种统计分析、集群分析、频谱分析等自动识别与分类。
通常利用专用数字图像处理系统来实现,且依据目的不同采用不同算法和技术。
2 遥感图像处理中相关的数学模型。
在遥感图像处理中,几乎在每个部分都有关于数学模型的应用。其中在图像恢复、影像增强以及信息提取中的应用最为广泛。
2.1 图像恢复中的数学模型。
由于遥感图像在接收过程中受到传感器响应特性、太阳高度、地形倾斜以及大气吸收与散射的影响,这就需要对在应用遥感影像之前对影像进行预处理。预处理过程主要是对图像进行辐射校正。根据影响因子的不同,校正过程分为系统辐射校正与大气辐射校正,用到了不同的数学模型。
2.1.1系统辐射校正中的数学模型。
传感器响应特性与太阳辐射两方面引起系统辐射畸变。因此,系统辐射校正从两面进行校正。
1)传感器响应特性引起的畸变校正模型。
由于光电变换系统的灵敏度特性通常有很高的重复性,所以可以定期地在地面测量器特性,根据测量值可对其进行辐射畸变校正。如对landsat 卫星的mss图像和tm图像可以按公式(2-1)的数学模型对传感器的输出()进行校正。
式中—已校正过的亮度值;
—传感器输出的辐射亮度值;
和—探测器能够输出的最大和最小辐射亮度。
—最大辐射分辨率。
探测器增益变化引起的辐射误差通常采用楔校准处理方法加以消除,现以陆地卫星可见光波段为例,校准模型为:
式中—校准后的输出亮度值;
—未校准的输入亮度值;
—滤波增益,取决于传感器响应因素;
—滤波偏移值,取决于传感器系统大气干扰;
—太阳角校正系数,在此作为常数。
2)太阳辐射引起的畸变校正模型。
太阳高度角引起的畸变校正是将太阳光线清晰恩照射时获取的图像校正为太阳光线垂直照射时获取的图像。太阳高度角可根据成像时间、季节和地理位置来确定,即:
式中—太阳赤纬(成像时太阳直射点的地理纬度);
—时角(地区经度与成像时太阳直射点地区经度的经差).
2.1.2 大气辐射校正中的数学模型。
大气辐射校正主要是消去程辐射,在消去程辐射过程中,主要采用图像分析模型与回归分析模型。由于程辐射主要发生在短波波长,把近红外波段作为无散射影响的标准图像,通过对不同波段图像的对比分析来计算大气影响。设红外波段为a,现需求其他波段相应的亮度最小值,这些波段设为b.
分别以a,b波段的像元亮度值为坐标,做二维光谱空间,两个波段中对应像元在坐标系内用一个点表示。结合图像分析模型,得到的回归模型如下:
式中—直线在轴上的截距;
—斜率,且。
其中,,分别为a,b 波段亮度的均值。
式中—波段a中的亮度为零处在波段b中所具有的亮度,可认为就是波段b的程辐射度。
校正的方法就是将波段b中的所有像元值都减去这个截距值,来改善图像,去掉程辐射。
2.2 遥感影像增强中的数学模型。
对遥感影像预处理之后,由于图像中含有噪声影响了用户的分析识别,这就需要图像增强。图像增强的主要目的是改变图像的灰度等级,提高图像对比度;消除边缘和噪声,平滑图像,突出边缘或线状地物,锐化图像;合成彩色图像;压缩图像数据量,突出主要信息等。图像空间域增强分为点增强与邻域增强。
在对不同范围的增强中用到了许多不同的数学模型,比如线性变换模型与非线性变换模型、卷积函数模型、不同算子函数模型等。
2.2.1 图像点增强中的数学模型。
对图像中的点增强根据运算中的变换函数的不同分为线性变换模型与非线性变换模型。线性变换模型为简单的线性关系式:
式中—原始图像的灰度图像值变量;
—扩展后的灰度值变量;
—斜率(扩展系数);
常数。经过线性变换模型后起到对图像灰度值拉伸的效果。
非线性变换的函数很多,如对数变换、指数变换、平方根变换、三角函数变换等,常用的有指数变换和对数变换。
指数变换函数的意义是在灰度值较高的部分扩大灰度间隔,属于拉伸,而在灰度值较低的部分缩小灰度值间隔,属于压缩,其数学模型为:
式中—变换前图像每个像元的灰度值;
—变换后图像每个像元的灰度值;
—可调参数,可以改变指数函数曲线的形态,从而实现不同的拉伸比例。
对数变换函数与指数变换相反,它的意义是在灰度值较低的部分拉伸,而在灰度值较高的部分压缩,其数学模型为:
式中符号意义同前。
2.2.2 图像邻域增强中的数学模型。
在对遥感图像利用之前需要空间滤波。它是以重点突出图像上的某些特征为目的,如突出边缘或线性地物等,也可以有目的地去除某些特征,如抑制图像上获取和传输过程中产生的各种噪声。在进行增强运算时,多采用空间卷积技术(又称掩膜技术),即在原图像上移动“活动窗口”,逐块进行局部运算,以实现平滑和锐化的目的。
卷积运算是在空间域上对图像进行邻域检测的运算。卷积运算的模板数学模型为:
在利用遥感影像信息提取时,需要先进行锐化处理。常见的锐化算子有roberts梯度算子、sobel梯度算子和laplace算子。这些算子都用到了梯度法,使用的是微分数学模型。
图像函数在像元点处的梯度定义为一个矢量,即:
梯度的模的数学模型为:
设,,则:
这种梯度模型又称为水平垂直差分模型,另外有一种罗伯特梯度模型,它是一种交叉差分计算模型,它的梯度模型表达式为:
与前述不同,拉普拉斯算子属于二阶导数算子,即:
对于离散的数字图像,二阶导数可以用二阶差分近似计算,推导出laplace算子的数学模型表达式:
在遥感图像处理中,利用以上数学模型可以对图像中的每一个像元计算梯度值,最终产生一个梯度图像,达到突出边缘——锐化的目的。
2.3 遥感信息提取中的数学模型。
在对遥感图像进行以上处理之后,可以提取遥感影像中的信息进行实际应用。在对遥感信息提取中,都需要对遥感影像中的地物进行分类。鉴于常规的目视判读技术难以发挥卫星影像多波段和多时相的优势和克服其较低的空间分辨率的缺点,也不能满足实时处理大量信息的要求,使得应用数学模型进行数字遥感图像的计算机分类识别,具有越来越重要的意义,成为遥感图像处理和分析领域中最活跃的分支。
20世纪70年代的计算机分类识别的数学模型,根据不同的地物类别具有不同的光谱反射特征和遥感数据固有的随机性特点,选择不同的时相和波段组合,应用统计模式识别方法,对图像像元进行逐点分类。这些数学模型中常用的有最大似然分类模型、最小距离分类模型、平行六面体模型等。这一类数学模型简单易行,但是识别精度不稳定,这就要求有进一步的改善。
20世纪80年代的计算机分类识别模型,就是在更全面准确地模拟判读员这些经验的基础上诞生的。这些数学模型大致分为两类:
l)第一类数学模型在图像识别中除考虑各地类像元光谱特征外,着重于像元空间特征的描述。空间特征可以分为两种类型:纹理特征和上下文特征。
纹理特征指的是图像色调的变化,它反映了相邻像元之间光谱值变化的统计关系;上下文特征指的是一个(或一组)像元与图像内其它像元的空间关系,这种空间关系通过像元之间的距离、方向、连通性和内含性得到反映。随着遥感卫星空间分辨率大幅度提高,像元的空间信息越来越丰富,建立纹理和上下文信息提取和分析的数学模型对于改善分类精度具有越来越重要的意义。
数学模型作业
问题 chapter 4 2 解 将大学生数量为34,29,42,21,56,18,71的区分别标号为1,2,3,4,5,6,7区,划出区与区之间的如下相邻关系图 记为第i区大学生人数,用0 1变量表示 i,j 区的大学生由一个销售 点 图书 i 建立lindo模型 max 63x12 76x13 ...
数学模型作业
2007011207 无78王萌。7.要从雨中从一处沿直线跑道另一处,若雨速为常数且方向不变,试建立数学模型讨论是否跑得越快,淋雨量越烧。将人体简化为一个长方体,高a 1.5m 颈部以下 宽b 0.5m,厚c 0.2m。设跑步距离d 1000m,跑步最大速度vm 5m s,雨速u 4m s,降雨量w...
数学模型作业
问题分析若要为该电力公司制定本月和下月的生产经营计划,则必须在不损坏发电站的发电设施的情况下使发电站能得到的利润最大化,而能得到的利润又与水库a中的水以及水库b中的水有直接联系,其中发电站a每千度电需水25m,发电站b每千度电需水20m,而由水电站的地理位置,可知水库a的水在a电站发完电后还要流入b...