遥感图像识别与信息提取

关于遥感图像识别与信息提取第1页，共90页，2023年，2月20日，星期六内容提要1.遥感图像地物特征与识别遥感图像地物特征典型地物的反射光谱特征2.遥感图像的目视解译目视解译标志目视解译方法与步骤3.遥感数字图像信息提取遥感数字图像的概念遥感数字图像处理4.遥感数字图像定量信息反演

第2页，共90页，2023年，2月20日，星期六1遥感图像地物特征与识别（1）遥感图像地物特征地物的反射光谱特性地物的发射光谱特性地物的透射光谱特性（2）典型地物的反射光谱特征第3页，共90页，2023年，2月20日，星期六（1）遥感图像地物特征

在可见光与近红外波段，地表物体自身的辐射几乎等于零。地物发出的波谱主要以反射太阳辐射为主。太阳辐射到达地面之后，物体除了反射作用外，还有对电磁辐射的吸收作用。电磁辐射未被吸收和反射的其余部分则是透过的部分，即：到达地面的太阳辐射能量＝反射能量＋吸收能量＋透射能量一般而言，绝大多数物体对可见光都不具备透射能力，而有些物体如水，对一定波长的电磁波透射能力较强，特别是对0.45~0.56μm的蓝绿光波段，一般水体的透射深度可达10~20m，清澈水体可达100m的深度。对于一般不能透过可见光的地面物体，波长5cm的电磁波却有透射能力，如超长波的透射能力就很强，可以透过地面岩石和土壤。第4页，共90页，2023年，2月20日，星期六（1）遥感图像地物特征地物的反射光谱特性反射率：地物的反射能量Pe占总入射能量Po的百分比，称为反射率ρ反射类型镜面反射（Specularreflection):入射波与反射波在同一平面内，入射角与反射角相等时，所形成的反射现象漫反射（Diffusereflection):无论入射波方向如何，其反射波分散到各个方向第5页，共90页，2023年，2月20日，星期六（1）遥感图像地物特征方向反射：朗伯体表面实际上是一个理想化的表面。它被假定为介质是均匀的、各向同性的，并在遥感中多用以作为近似的自然表面。

电磁波反射的三种形式第6页，共90页，2023年，2月20日，星期六（1）遥感图像地物特征2地物的发射光谱特性发射率:地物发射电磁辐射的能力，以黑体辐射作为基准－指单位面积上观测地物发射的某一波长的辐射通量密度；－指与观测地物同温度下黑体的辐射通量密度；第7页，共90页，2023年，2月20日，星期六（1）遥感图像地物特征

通常，根据发射率与波长的关系，将地物分为三种类型

黑体，其发射率ελ=1，即黑体发射率对所有波长都是一个常数，并且等于1。灰体，其发射率ελ=常数＜1。即灰体的发射率始终小于1，发射率ελ不随波长变化。选择性辐射体，其发射率ελ＜1，发射率ελ随波长而变化。地物的发射波谱特性曲线

第8页，共90页，2023年，2月20日，星期六（1）遥感图像地物特征太阳辐射第9页，共90页，2023年，2月20日，星期六（1）遥感图像地物特征3地物的透射光谱特性透射率即地物透射的能量与入射总能量的百分率，称之为透射率

τ=Eτ/Eο×100%3第10页，共90页，2023年，2月20日，星期六（2）典型地物的反射光谱特征1岩石的反射光谱特征

岩石的波谱特征是地质遥感的基础，不同的矿物成分、矿物含量、风化程度、含水状况、颗粒大小、表面的光滑程度、色泽等都会影响到其反射波谱特征。第11页，共90页，2023年，2月20日，星期六（2）典型地物的反射光谱特征2、土壤的反射光谱特征自然状况的土壤表面的反射率没有明显的峰值和谷值，一般来说土质越细，反射率越高，有机质含量越高和含水量越高反射率越低。此外土壤的肥力也会对反射率产生影响。第12页，共90页，2023年，2月20日，星期六（2）典型地物的反射光谱特征不同含水量的土壤反射光谱曲线

三种土壤反射波谱曲线比较

第13页，共90页，2023年，2月20日，星期六（2）典型地物的反射光谱特征3、水体的反射光谱特征水体的反射主要在蓝光波段，其他波段吸收都很强，特别在近红外以后水体便成为一个吸收体。光谱反射特性可能包括来自三方面的贡献：水的表面反射、水体底部物质的反射和水中悬浮物质的反射。光谱吸收和透射特性不仅与水体本身的性质有关，而且还明显地受到水中各种类型和大小的物质——有机物和无机物的影响。第14页，共90页，2023年，2月20日，星期六（2）典型地物的反射光谱特征水体的反射光谱特征

第15页，共90页，2023年，2月20日，星期六（2）典型地物的反射光谱特征4、植被的反射光谱特征

植被的反射波谱曲线规律性明显而独特，主要分为三个波段：可见光波段（0.4~0.76m)有一个小的反射峰，位于0.55m(绿光波段）处，两侧蓝光波段（0.45m

）和红光波段(0.67m)则有两个吸收带；近红外波段（0.7~0.8m）有一反射徒坡，至1.1m附近有一峰值；中红外波段（1.3~2.5m）受绿色植物含水量的影响，吸收率大增，反射率下降，在1.45、1.95和2.7m为中心是水的吸收带，开成低谷。第16页，共90页，2023年，2月20日，星期六（2）典型地物的反射光谱特征矿区红杉林反射曲线的蓝移现象

绿色植物有效光谱响应特征

第17页，共90页，2023年，2月20日，星期六2遥感图像目视解译

目视解译是用肉眼或借助于简单的工具如放大镜、立体镜、投影观察器等，直接由肉眼来识别图像特性，从而提取有用信息，即人把物体与图像联系起来的过程。具备的基本知识：专业知识、地理区域知识、遥感系统知识。

第18页，共90页，2023年，2月20日，星期六（1）目视解译标志1、直接解译标志色调（tone）：色调是识别目标地物的基本依据。由于地物属性的差别，在遥感图像上表现出色调上的差别。一般来说由于人眼的局限性，在图像可分出16个灰度级。颜色（colour）：由于目标地物在不同波段中反射或发射电磁辐射能量的差异性，由此而表现在彩色图像颜色的差异性。颜色也是目视解译的基本标志之一。图型（pattern）：目标地物有规律排列而成的图形结构。第19页，共90页，2023年，2月20日，星期六（1）目视解译标志阴影（shadow）：阴影是遥感图像上由于电磁辐能量被遮挡而产生的辐射能量减弱。由于阴影的存在，可据此地物的性质或高度，应注意的是阴影的形状与大小受到辐射能量入射角的影响。形状（shape）：是指目标地物在遥感影像上呈现的外部轮廓。由于不同地物的顶视平面的差异，可据此判断目标地物的性质。纹理（texture）：是指遥感图像中目标地物内部色调有规则变化造成的影影结构。大小（size）：指遥感图像上目标地物的形状、面积与体积的度量关系。是遥感图像上测量目标地物最重要的数量特征之一。同时根据其大小可以推断地物的属性。第20页，共90页，2023年，2月20日，星期六（1）目视解译标志2、间接解译标志位置（site）：指目标地物分布的地点。任何地物与其周围地理环境（geographicalEnvironment）或生态环境(ecologicalEnvironment)总是存在着一定的空间关系，并受到周围环境的某种程度的制约。相关布局（association）：是指多个目标地物间的空间配置关系。第21页，共90页，2023年，2月20日，星期六2遥感图像解译方法与步骤1、目视解译的认知过程自下向上过程图像信息获取特征提取识别证据选取自上向下过程特征匹配提出假设图像辨识第22页，共90页，2023年，2月20日，星期六2遥感图像解译方法与步骤2、图像解译方法遥感资料的选择及影像处理遥感资料的选择资料类型选择波段选择时间选择比例尺选择遥感图像的处理影像放大影像数字化图像处理第23页，共90页，2023年，2月20日，星期六2遥感图像解译方法与步骤

目视解译的方法

直判法是根据遥感影像目视判读直接标志，直接确定目标地物属性与范围的一种方法。对比分析法对比分析法包括同类地物对比分析法、空间对比分析法和时相动态对比法。不仅是同一遥感影像图进行对比，而且要借助不同时相的遥感影像图之间进行对比。综合推理法综合考虑遥感图像多种解译特征，分析、推理某种目标地物的方法。第24页，共90页，2023年，2月20日，星期六2遥感图像解译方法与步骤信息复合法利用专题图和地形图等信息与遥感影像图重合，根据这些辅助信息识别遥感影像图上目标地物的方法。地理相关分析法根据地理环境中各种地理要素之间的相互依存，相互制约的关系，借助专业知识，分析和推理某种地理要素性质、类型、状况与分布的方法。第25页，共90页，2023年，2月20日，星期六2遥感图像解译方法与步骤

3、解译步骤准备工作

包括资料收集、分析、整理和处理初步解译、建立解译标志包括路线路勘，制订解译对象的专业分类系统和建立解译标志

室内解译野外验证

包括解译结果校核检查，样品采集和调绘补测成果整理

包括编绘成图，资料整理和文字总结第26页，共90页，2023年，2月20日，星期六2遥感图像解译方法与步骤目视解译步骤第27页，共90页，2023年，2月20日，星期六3遥感数字图像信息提取1遥感数字图像的概念图象、数字图象遥感数字图象2遥感数字图像处理遥感数字图像的获取遥感数字图像预处理遥感数字图像的变换、增强和融合遥感数字图像分析第28页，共90页，2023年，2月20日，星期六(1)遥感数字图像的概念图像、数字图像

物理世界中客观对象的一种表示数字图像客体或可见图像的数字表述。它实际上是具有某种数值的一些点按行(横)和列(纵)排成的二维矩阵模拟图像指空间坐标和明暗程度都连续变化的、计算机无法直接处理的图像第29页，共90页，2023年，2月20日，星期六(1)遥感数字图像的概念遥感数字图像指以遥感方式获得的以数字形式表述的遥感影像遥感数字图像最基本的单位是像素像素的属性特征常用灰度值来表示，即该像素位置上亮暗程度的整数值

第30页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理一、遥感数字图像的获取

通过接收、记录目标物的电磁波特征的仪器，即传感器获得的

数字图像获取过程第31页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理二、遥感数字图像预处理

1、辐射校正消除图像数据中依附在辐射亮度里的各种失真的过程称为辐射校正。完整的辐射校正包括遥感器校正、大气校正，以及太阳高度和地形校正。

遥感器纠正：遥感器的设计大气辐射纠正：地形辐射纠正：需要DEM地物反射模型纠正：需要和成像时刻取得同步的地面地物光谱测量数据第32页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理

大气纠正方法以红外波段最低值校正可见光波段回归法经验线性法第33页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理以红外波段最低值校正可见光波段前提假设：大气散射的影响主要在短波波段，红外波段中清洁的水体几乎不受影响，反射率值应当为0。由于散射影响，而使得水体的反射率不等于0，推定是由于受到了天空辐射项的影响。直方图法确定纠正方法：差值法第34页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理回归法

原理：选择可见光和红外波段进行2维散点图，建立线性回归方程。第35页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理经验线性法第36页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理2、几何校正

造成几何位置的畸变有原因：遥感器本身引起的畸变外部因素引起的畸变处理过程中引起的畸变

第37页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理

遥感器本身引起的几何畸变与遥感器的结构、特性和工作方式不同而异。这些因素主要包括：

1）

透镜的辐射方向畸变像差；

2）透镜的切线方向畸变像差；

3）

透镜的焦距误差；

4）

透镜的光轴与投影面不正交；

5）

图像的投影面非平面；

6）

探测元件排列不整齐；

7）

采样速率的变化；

8）

采样时刻的偏差;9)扫描镜的扫描速度变化。

遥感器本身引起的畸变第38页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理影响图像变形的外部因素包括：

1）

地球的曲率

2）

大气密度差引起的折光

3）

地形起伏

4）

地球自传

5）遥感器轨道位置和姿态等

外部因素引起的畸变第39页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理遥感器轨道位置和姿态引起的误差地球自传引起的误差第40页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理地球曲率和地形起伏引起的误差遥感器轨道位置和姿态引起的误差第41页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理

遥感图像再处理过程中产生的误差，主要是由于处理设备产生的噪声引起的。传输、复制光学数字

处理过程中引起的畸变第42页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理遥感图像的几何纠正方法遥感图象的几何粗处理和精处理遥感图像的几何纠正按照处理方式分为光学纠正和数字纠正光学纠正主要用于早期的遥感图像的处理中，现在的应用已经不多。除了对框幅式的航空照片（中心投影）可以进行比较严密的纠正以外，对于大多数动态获得的遥感影像只能进行近似的纠正第43页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理数字图象几何纠正：通过计算机对离散结构的数字图像中的每一个像元逐个进行纠正处理的方法。几何精校正是指利用地面控制点使遥感图像的几何位置符合某种地理系统，与地图配准，并调整亮度值。也就是在遥感图像的像元与地面实际位置之间建立数学关系，将畸变图像空间中的全部像元转换到校正图像空间去。内容一般包括两个方面：一个是图像像元空间位置的变换，另一个是像元灰度值的重采样。

第44页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理几何精校正具体步骤：选取地面控制点（GCP）选择空间变换函数重采样和内插最邻近内插法双线性内插三次卷积内插第45页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理数字图象的纠正过程

纠正的函数可有多种选择：多项式方法、共线方程方法、随机场内插方法等等。其中多项式方法的应用最为普遍。第46页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理三、遥感数字图像的变换、增强和融合

1.遥感数字图像变换概念：图像变换指的是将图像从空间域转换到变换域例如频率域的过程。进行图像变换的目的就是为了使图像的处理过程简化。作用：通过图像变换简单而有效的实现增强处理通过图像变换可以对图像进行特征抽取。第47页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理快速傅立叶变换（FFT）目的进行数据压缩、图像的增强、特征提取方法信号处理中的频率域分析方法步骤选择适当的变换函数进行傅里叶变换分析变换的结果进行傅里叶逆变换第48页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理K－L变换目的减少图像波段之间的相关性，去除多余的信息，减少图像的数据量方法统计学中的正交变换方法性质和特点步骤进行数据统计进行主成分分析进行旋转变换和逆变换第49页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理K－T变换目的分离和消除干扰信息突出研究的专题信息方法几何中的坐标旋转方法性质和特点步骤作两个波段的散点图分析灰度值的变换特点第50页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理2.遥感数字图像增强

1、反差增强图像反差增强又称对比度增强灰度拉伸根据原图像的直方图确定需要做拉伸变换的灰度值区间，然后把这一（或一些）灰度值区间按某种直线或曲线方程关系拉伸或压缩而成为变换后的灰度值区间直方图线性拉伸直方图分段线性拉伸非线性拉伸变换

第51页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理直方图均衡化

变换后的直方图接近均匀分布。即图象中每一灰度级的像元数目大致相同。使得面积最大的地物细节得以增强，而面积小的地物与其灰度接近的地物进行合并，形成综合地物。减少灰度等级换取对比度的增大。原直方图：均衡后直方图：第52页，共90页，2023年，2月20日，星期六原图及其直方图均衡后的图和直方图第53页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理直方图匹配

将原始图象转换为给定直方图的图象各种拉伸以后虽然对于感兴趣的地物提高了分辨能力，但同时也造成了信息损失。第54页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理2、空间域图像增强空间域图像增强处理是应用某种数学模式直接改变图像像元灰度值的变换。这种变换与像元的坐标无关，只改变像元点的灰度深浅程度。即改变了图像的对比度，这样图像中的某些信息被突出（可能也有些信息被压抑了）达到了图像增强的目的。可分为单点处理和邻域处理

第55页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理边缘增强突出目标的轮廓或边缘信息，主要是通过微分法实现的梯度法和拉普拉斯算子法平滑滤波低频增强的空间域滤波技术。它可以滤掉由于孤立的单点噪声而引起的灰度偏差中值滤波定向滤波

第56页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理3、频率域图像增强

高通滤波（锐化）低通滤波（平滑）带通滤波（突出地物）同态滤波（改善图像质量）第57页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理4、代数运算增强

基本代数运算差值运算两个波段的差值比值运算两个波段的比值应用第58页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理差值运算利用不同地物之间光谱的特征有利于突出目标与背景反差小的信息红外波段-红波段，突出植被实例B4-B3第59页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理比值运算两个波段的比值目的减轻地形的干扰如果阴坡的坡度大于太阳高度角，则难于处理扩展地物，特别是植被间的光谱差异抑制大气的影响减少了数据之间的相关性第60页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理植被指数NDVITM4/TM3常用的比值第61页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理二应用实例NDVI第62页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理TM图像的5/2突出水体第63页，共90页，2023年，2月20日，星期六

（2）遥感数字图像处理5、彩色图像增强假彩色密度分割将灰度按照指定的间隔分割为不同的级，对新的密度分级分别赋予不同的颜色第64页，共90页，2023年，2月20日，星期六第65页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理②彩色增强在一个输入像元的灰度级上分别进行三个独立的色变换在显示屏幕上产生一个彩色合成图像。从而达到明显的分辨和识别效果

？第66页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理③假彩色合成目的综合不同波段的特征，突出研究对象的差异原理原则信息量最大相关性最小差异最大方法信息量分析相关分析方差分析其他第67页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理TM图像4-3-2合成实例第68页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理3.遥感数据融合实质在统一地理坐标系中将对同一目标检测的多幅遥感图像数据采用一定的算法，生成一幅新的、更能有效表示该目标的图像信息原理及过程预处理主要包括遥感影像的大气校正、辐射校正及空间配准数据融合根据融合目的和融合层次智能地选择合适的融合算法，将空间配准的遥感影像数据（或提取的图像特征或模式识别的属性说明）进行有机合成，得到目标的更准确表示或估计第69页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理数据融合分类像元级融合像元级融合是一种最低水平的融合特征级融合特征级融合是一种中等水平的融合决策级融合决策级融合是最高水平的融合数据融合方法IHS变换、代数法、Brovey变换、图像回归法、主成分变换（PCT）、小波变换等第70页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理三级融合层次下的融合方法

像元级特征级决策级代数法熵法专家系统IHS变换表决法神经网络小波变换聚类分析Bayes估计Brovey变换Bayes估计模糊聚类法主成分变换神经网络法可靠性理论回归模型法加权平均法基于知识的融合法Kalman滤波法Dempater-shafer推理法Dempater-shafer推理法第71页，共90页，2023年，2月20日，星期六TM741RadarSatTM和RadarSat复合第72页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理四、遥感数字图像分析图像理解从图像中提取有用信息的过程包括图像分类图像分割特征分析输出内容明确的数据或图形图像特点以具体的知识和目的为基础，结果更接近于解释的成果和实际的应用第73页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理图像信息获取信息处理特征提取分类

图像分类第74页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理1、遥感图像分类技术原理同类地物在相同的条件下（光照、地形等）应该具有相同或相似的光谱信息和空间信息特征。不同类的地物之间具有差异同类地物的像元在数目较大时，其特征量的分布类型接近正态分布方法利用像素的灰度值之间的关系，综合考虑矢量数据，考虑属性数据。对像素的进行分类以随机变量的统计分析为基础第75页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理特点多变量的图像分类分类在一定的比例尺上进行原则多变量综合考虑在特征空间的位置可以用均值表示，离散程度可以用方差表示分类的实质是将特征空间分为若干子区域，每个子区域为一类第76页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理指标（特征）波段其他属性地形特征调查数据…统计量均值方差距离度量离差平方和第77页，共90页，2023年，2月20日，星期六（2）遥感数字图像处理2、具体分类方法非监督分类最大似然