遥感图像处理-基本原理-厉青课件

二、遥感图像处理基本原理第一章概述第二章遥感图像处理基本原理第三章遥感图像预处理第四章图像增强第五章

图像分类二、遥感图像处理基本原理第一章概述第一章概述一、遥感图像处理技术的概况和发展二、遥感数字图像处理的主要内容、特点及应用三、数字图像处理设备和系统第一章概述一、遥感图像处理技术的概况和发展第一章概述

一、遥感图像处理技术的概况和发展图像处理技术是随着计算机技术的发展而开拓出来的一个新的应用领域，就是把图像转换成一个数据矩阵存放进计算机对之进行处理。具体说，利用数字计算机或其它数字硬件，对从图像信息转换而来的电信号进行某些数学运算，以期提高图像的实用性，从而达到人们所要求的某些预期的结果。第一章概述

数字图像处理起源于20世纪20年代，首先用于改变数字化报纸的图片。新闻照片经过数字化，由海底电缆开始在纽约和伦敦之间传送，为了用电缆传输图片，首先要进行编码，然后在接收端由专用打印设备进行图片重建。此时的技术能力，仅把图像编码为五个亮度级，到1929年才提高到15个等级。但那时图像的处理还不是用计算机进行的。从六十年代起，随着电子计算机技术的进步，数字图像处理获得了飞跃的发展。1964年美国喷气推进实验室（JPL）开始用计算机处理航天照片，以消除畸变，提高质量，取得了惊人的效果。从此以后，遥感数字图像处理技术和设备得到了迅速的发展。

硬件：主处理机：①60年代到70年代上半期，遥感图像处理的实现一般都是在已有的大、中型计算机上配备必要的外围设备，利用图像处理软件来完成图像处理与分析功能。②70年代下半期到80年代初期，专用图像处理机的出现，将许多图像处理算法固化，一部分图像处理与分析功能可以用硬件快速实现，从而使遥感图像处理速度大幅度地提高。③80年代初到90年代初期，超小型机在图像处理中占有统治地位，图像运算速度也有了很大的提高。④近年来，在遥感图像处理系统的构成方面，硬件发展已进入了第四个阶段，即在工作站上建立大规模专业化的遥感图像处理系统，如建立在工作站上的ERMAPPER。而一般的应用用户，又可以建立以PC机为核心遥感图像处理系统，专用图像处理机在遥感图像处理系统中地位已经开始下降。软件：软件功能反映了遥感技术应用水平①遥感图像处理系统与GIS的融合②多种遥感资料的处理。③非遥感资料的图像处理。硬件：二、遥感数字图像处理的主要内容、特点及应用

1、数字图像处理的内容

a、图像数字化及编码压缩Imagedigitization.Coding.Compressionb、图像恢复或校正Imagerestorationorcorrectionc、

图像增强Imageenhancementd、图像分析（或图像分割）ImageanalysisorImagesegmentation

二、遥感数字图像处理的主要内容、特点及应用2、遥感数字图像处理的优点

a、图像信息损失低，处理的精度高b、抽象性强，再现性好ｃ、通用性广，灵活性高ｄ、能充分利用遥感图像信息，对图像灰度可分解为很多级，可检测出微小的图像细节，并能对图像信息作定量分析。而目视判读对黑白图像仅能分辨出十余级灰度。ｅ、处理周期短，速度快，适于进行动态分析。而且计算机数字图像处理是以数学理论为基础的，方法广泛，可以充分开发利用不同形式的遥感图像信息。计算机技术的发展和普及是当今世界的大趋势，计算机数字处理正是顺应这个大趋势的，其图像处理的技术、方法及质量可不断发展提高。

2、遥感数字图像处理的优点3、遥感数字图像处理的应用

遥感数字图像处理，是遥感技术的主要组成部分，它在遥感技术应用中起着重要的作用．1．利用遥感数字图像处理技术获得满足一定精度要求的各种图件2．利用遥感数字图像处理技术，快速准确地提取所需信息3．利用遥感数字图像处理技术，为图像的解译奠定基础3、遥感数字图像处理的应用

三、数字图像处理设备和系统

遥感图像处理系统

通用计算机型

专用处理机型

运算处理设备

输入设备

输出设备

存储器（外存）

图像软件

系统软件应用软件

三、数字图像处理设备和系统

遥感图像处理软件介绍

1

、ERDASIMAGINGERDASIMAGING是美国ERDAS公司开发的遥感图像处理系统，它以其先进的图像处理技术，友好、灵活的用户界面和操作方式，面向广阔应用领域的产品模块，服务于不同层次用户的模型开发工具以及高度的RS/GIS集成功能，为遥感及相关应用领域的用户提供了内容丰富而功能强大的图像处理工具，代表了图像处理系统未来发展趋势。

ERDASIMAGING软件构成：核心模块IMAGNEEssentials：基本地图生产系统，提供了众多工具生产以图像可视化与分析（增强、分类、专题图制作、三维景观显示等），信息集成（影像库管理，矢量数据查询检索等）和简单几何纠正为主的高质量专题地图，另外还拥有易于使用应用界面扩展工具EML宏语言。IMAGNEAdvantage：一个易于使用的集成环境，提供图像镶嵌、地表DEM内插、高级图像解译（增强、滤波、高光谱图像处理、傅立叶变换）和正射纠正工具，对图像数据进行综合分析。IMAGNEprofessional：快速图像检测与分类，特别功能包括无需编程就能扩展图形空间数据建模、最大似然、模糊决策等，提供雷达图像基本分析和完整而专业易用的分类工具。ERDASIMAGING软件构成：2、ENVI

ENVI是美国RSI开发的一套功能齐全的遥感图像处理系统，是处理、分析并显示多光谱数据、高光谱数据和雷达数据的高级工具，能够进行地物目标识别，直观且功能强大，完全由IDL开发，方便灵活，可扩展性强，获2000年美国权威机构NIMA遥感软件测评第一。它具备如下特点：接受更多的传感数据格式

ENVI对于要处理的图像波段数据没有限制，可以处理最先进的卫星格式，如Landsat-7、IKONOS、SPOT、RADARSAT、NASA、NOAA、ERSO和TERRA，并准备接受未来所有传感器的信息。2、ENVI强大的多光谱影像处理功能

ENVI能够提取图像信息，具备进行图像预处理、图像增强、掩膜、图像计算、图像分类及后处理、图像变换和滤波、图像镶嵌、融合等功能的全套完整的遥感影像处理工具。此外，ENVI提供丰富完备的投影软件包，可支持各种投影类型。同时，ENVI还创造性地将一些高光谱数据处理方法用于多光谱影像处理。更便捷地集成栅格和矢量数据

ENVI包含所有基本的遥感影像处理功能，如校正、定标、波段运算、分类对比增强、滤波、变换、边缘检测及制图输出功能，并加注汉字。ENVI具有对遥感影像进行配准和正射校正的功能可以给影像添加地图投影，并与各种GIS数据套合。ENVI的矢量工具可以进行屏幕数字化、栅格和矢量叠合，建立新的矢量层、编辑点、线、多边形数据，缓冲区分析，创建并编辑属性并进行相关矢量层的属性查询。强大的多光谱影像处理功能高光谱影像处理

ENVI拥有先进的高光谱和多光谱分析工具，可识别出图像中纯度最高的像元，通过与已知波谱库的比较确定未知波谱的组成。使用ENVI自带的波谱库，也可自定义波谱库，甚至可以组合使用线形波谱分离和匹配滤波技术进行亚像元分解，以消除匹配误差，得到精确的结果快速处理雷达数据ENVI可以快速处理SAR数据，提取CEOS信息并浏览RADARSAT和ERS－1数据。用天线阵列校正、斜距校正、自适应滤波等功能提高数据的利用。纹理分析还可分段分析SAR数据。如果是极化雷达数据，可从SIR－C和AIRSAR压缩数据中选择自己的极化和工作频率。此外，还可以浏览和比较感兴趣区的极化信号，并创建幅度图像和相位图像。高光谱影像处理可扩展性强，拥有第四代可视化语言IDL

ENVI具有强大的底层二次开发工龄－IDL（InteractivedataLanguage）交互式数据语言开发。IDL语言是面向矩阵、语法简单的第四代可视化语言，致力于科学数据的可视化和分析，是跨平台应用开发的最佳选择。IDL集可视、交互分析、大型商业开发为一体，为用户提供了完善、灵活、有效的开发环境。ENVI的完全开发版本具有IDL底层环境，可随时与IDL进行衔接。用户可以用第四代科学计算评议为其编写扩展功能。

可扩展性强，拥有第四代可视化语言IDL3、PCI

PCI软件是加拿大PCI公司开发的用于图像处理、几何制图、GIS、雷达数据分析以及资源管理和环境监测的多功能软件系统。它以其丰富的软件模块、灵活的编程能力和便利的数据可操作性代表了图像处理软件系统的发展趋势和技术先导，支持所有的数据格式、适用于各种硬件平台。PCI不仅可用于卫星和航空遥感图像的处理，还可应用于石油勘探、矿产资源开发、林业、农业、测绘、土地资源调查评估与管理、自然灾害动态监测、环境保护、城市规划、铁路交通、沙漠治理、工程建设、气象预报等非常广泛的领域。3、PCI

PCI不断地将最新、最前沿的科学研究成果转化为软件产品之中。它是支持多传感器正射投影校正和从立体像对中提取DEM图像的软件；最早提供神经网络和模糊逻辑的精确分类方法；率先能够处理多达1024个通道的多光谱数据及推出处理极化数据和模拟图像的RADAR软件模块；PCI软件能够运行于多种平台且拥有最齐全的功能模块：经过十几年的发展，PCI公司的产品结构发生了极大的变化，已经从一个单纯的遥感图像处理软件产品，发展到包括GIS，WEBGIS，专业制图、数字摄影测量等系列产品的地理资讯公司。从图像处理的产品来讲，PCI总体已形成三个系列的产品。即专业遥感图像处理产品、专业雷达信号处理及分析产品、数字摄影测量产品。其中，PCI专业遥感图像处理产品按用户不同的需求，包括从普通应用到专业分析共三个商业软件包及五个专业扩展模块。PCI不断地将最新、最前沿的科学研究成第一软件包IMAGEW0RKS该软件包主要由三部分组成：影像处理和多光谱分类工具(IMAGE—WORK／MULTISPECTRALCLASSIFICATION)、几何校正工具(GCPWORKS)和直接数据读取(GEOGATEWAY)。第二软件包EASI／PACEIMAGEPROCESSINGKIT／VISUALMODE—ILER该软件包在包含了IMAGEWORKS的全部功能的基础上又增加了VISUALMODUDLE、XPACEKENNAL、IMAGEPROCESING、GEOMETRICCORRECTION、MULTILAYER、VECTORUTILTIES、ACK专业制图、地形分析、航片立体像对DEM提取和磁带输入输出的功能。此外，本软件包可以追加PCI专业扩展模块

第三软件包EASI／PACEIMACEPB0EASI／PACEIMAGEPRO除包含EASI／PACEIMAGEPROCESSINGKIT／VISUALMODELLER软件包全部的功能外，还具有多光谱分析、雷达分析、AVHRR轨道领航者、大气校正、高光谱分析、神经元网络分类器、目标图像库管理者、三维可视化飞行模拟。并且，本软件包还可以追加PCI专业扩展模块。第一软件包IMAGEW0RKSPCI遥感图像处理软件的五个专业扩展模块该软件包括5个专业扩展模块，分别为大气校正(AtmosphericCorrection)、图像锁数据融合(Imagelockdatafusion)、极化雷达(PolarimetricRadarAnalysis)、PCI作者(PCIAuthor)、软件工具箱目标库(SoftwareToolboxObjectLibraries)。这五个模块的任意一个，均可以单独与其他软件包配合使用。PCI遥感图像处理软件的五个专业扩展模块第二章遥感图像处理基本原理第一节遥感成像过程第二节遥感图像模式及图像函数第三节遥感图像的数据结构第四节遥感图像的基本统计分析第五节遥感图像的信息特征第二章遥感图像处理基本原理第一节遥感成像过程第二章遥感图像图像的基本概念：图像（image）是对客观对象的一种相似性的描述或写真，它包含了被描述或写真对象的信息，是最主要的信息源。遥感图像是通过远距离探测而记录的地球表面、大气层、以及其它星球表面等物体在不同的电磁波波段所反射或发射的能量的分布和时空变化的产物。

第二章遥感图像根据波段的多少

单波段图像：图像上每个点只有一个亮度值多波段图像：图像上每个点不只一个特性超波段图像：图像上每个点具有几十甚至几百个特性根据波段的多少根据图像的明暗程度和空间坐标的连续性划分：

数字图像（digitalimage）

：被计算机存储、处理和使用的图像，是一种空间坐标和灰度均不连续、用离散数学表示的图像，它属于不可见图像

模拟图像(analog)：空间坐标和明暗程度都连续变化的、计算机无法直接处理的图像，它属于可见图像模拟图像数字图像模数转换，记A/D转换数字图像模拟图像数模转换，记D/A转换遥感数字图像（digitalimage）是指以数字形式表述的遥感影像。根据图像的明暗程度和空间坐标的连续性划分：

第一节遥感成像过程

对遥感图像进行数字处理其实质是通过各种运算及变化来提取和增强图像中的专题信息，以达到识别目标（地物、专题信息）的目的。遥感图像数字处理方法与方案的选择是依据遥感专题的目的和遥感图像本身的特征来确定的。而遥感图像的本身特征又在很大程度上取决于图像的形成过程，因而图像的形成过程对于遥感图像的数字处理有重要的意义。遥感图像处理-基本原理-厉青课件遥感图像处理-基本原理-厉青课件

我们所取得和使用的图像数据或像片是地物的辐射信息经过整个系统的记录、传送、加工及转化后的产品。从把运载工具发射到预定轨道后传感器开始工作起，到由专用的处理系统生产出遥感数字产品（计算机兼容磁带、磁盘等）和模拟产品（各种比例尺的遥感相片）的全过程和全部设备构成一个遥感系统。从把运载工具发射到预定轨道后传感器开始工作起，

决定遥感图像的信息内容和质量的最主要的因素是所用电磁波及相应的传感器的性能。轨道高度和瞬时视场（或地面分辨率）的大小也直接影响图像内容的详细程度和精度。遥感图像的获取的方式：摄影成像扫描成像雷达成像决定遥感图像的信息内容和质量的最主要的因素是第二节遥感图像模式及图像函数

一、遥感图像模式

遥感图像记录的是遥感传感器对地观测的结果，遥感图像的灰度值反映了地物的反射和发射电磁波能力，遥感图像的灰度值与地物的成分、结构等以及遥感传感器的性质之间存在着某种内在联系，这种内在联系可以用一函数关系来表达，也就是说客观存在遥感图像的模式。第二节遥感图像模式及图像函数作为一个普遍性的理论模式，遥感图像一般可以用这样一个公式来描述（一幅遥感图像上任一点（x，y）所记录的波谱辐射量L）：L（x,y;λ,t,p）=ρ(x，y;λ,t,p)I（x，y；λ,t）+[1-ρ(x，y;λ，t,p)]e(λ)反射的波谱辐射量发射的波谱辐射量式中L－(x,y)点的波谱辐射量；x,y－图像的空间位置坐标；t－图像的成像时间；λ－传感器电磁波工作波段的波长；p－电磁波段极化性质；ρ－(x,y)点的波谱反射率；I－(x,y)点的入射辐射量；e(λ)－黑体的波谱发射本领

即：遥感图像是关于观测区空间位置变化（x,y坐标）、传感器电磁波工作波段的波长、成像的时间以及电磁波极化性质的函数。

遥感图像处理-基本原理-厉青课件

讨论：①在可见光及近红外波段，物体自身发射的辐射量可以忽略不计，则L（x,y;λ,t,p）=ρ(x，y;λ,t,p)I（x，y；λ，t）此时，入射量I取决于太阳光照条件和传感器的几何特征，而ρ则反映物体的性质。②中热红外波段，中红外，反射和发射都要考虑，热红外，考虑发射。摄像时间是一个重要因素，夜晚摄取的主要是地物的热发射;白天的反射部分处于不同的波段，一般是通过不同的探测器分别记录为不同的图像。③微波遥感有主动雷达系统及被动的微波扫描系统，前者是利用人工辐射源产生的脉冲波经地面反射后的回波，而后者则是接收地物发射的微波热辐射。雷达系统接收的反射信息与一般的反射辐射具有不同的性质和特点。

二、图像函数固定成像时间t，如果不考虑L（x,y;λ,t,p）的物理意义，地物的辐射量度是一个关于空间坐标位置的函数，也就是说图像就是一个二维的空间分布，可以把图像函数L用关于（x,y）的实函数f（x,y）来表示。图像函数f（x,y）所代表的实际上是二维空间内物体所反射或发射的辐射能量的分布，是真实情况的直接反映，还不是传感器记录下来的数据或像片的图像本身。因为传感器所接收的电磁波辐射能量至少还要受到辐射电磁波与大气层的相互作用的影响，因此，传感器所记录下来的物体辐射特征及几何特征与实际的地物辐射特征之间还有差别，它们之间还存在着一个对应变换关系。即：g(x,Y)=Tf(x,y)

式中：g(x,Y)－描述传感器接受的电磁波辐射特征的函数，即图像函数f(x,y)－反映实际地物电磁波辐射特征的函数，即物面函数；T－f(x,y)与g(x,Y)之间的对应变换关系。即成像系统：二、图像函数景物图像成像系统

f(x,y)g(x,Y)

g(x,y)就是遥感图像处理中真正要遇到的图像函数，它是一个关于空间坐标（x,y）点地物的某种物理量的一种度量，它的大小称为图像点的亮度、灰度、密度等，它既反映了图像灰度的大小，也反映了图像灰度的分布。遥感图像函数g(x,y)有以下三个基本特点：①函数值物理意义的明确性。②函数定义域的限定性。③函数值的限定性。景物遥感图像函数g(x,y)有以下三个基本特点：①函数值物理意义的明确性遥感图像的灰度值即图像函数的值代表了地物电磁波辐射的一种度量，也就是说图像函数值主要反映的是地物的光谱特征。②函数定义域的限定性由于每一种遥感传感器都具有一定的视域，因而它所获得的图像大小也是有限的，所以图像函数只在图像范围内有效。即它在空间上是限定的，包含在一个矩形区域内，一般坐标原点在左上角。③函数值的限定性图像函数g(x,y)的取值都在一定的值域内，即0≤g(x,y)≤RmaxRmax—地物的最大辐射量，g(x,y)没有负值，而且不大于最大的辐射量。[0，Rmax]称为灰度空间，通常g(x,y)=0定为黑色，g(x,y)=Rmax定为白色所有中间值都是由黑连续地变为白时的灰度等级遥感图像函数g(x,y)有以下

在考虑图像模式L（x,y;λ,t,p）的简化过程中，图像函数g(x,y)只考虑了空间坐标位置的变化，在实际的遥感图像中，其它三个变量往往通过别的途径来表示：①电磁波波长λ的变化。②成像时间t的变化。③电磁波极化性质ρ的变化。

在考虑图像模式L（x,y;λ,t,p）的简化过程中，图三、图像函数的数字化图像函数g(x,y)的数字化就是把连续函数图像转换为数字图像，图像函数的数字化，通常包括以下两个方面的内容：三、图像函数的数字化图像函数g(x,y)的1、抽样：

进行像幅空间坐标的数字化，即沿像幅X轴和Y轴等距离地分割为M＊N个像元（即网格化），并分别量测这些像元的平均灰度值，这一过程称为抽样。2、量化：对量测的灰度值进行数字化。即将灰度值转换成二进制字码代表的某一灰度级。经过量化后，数字图像的灰度值不是像元的平均辐射值，而是像元平均辐射值所在的编码区间的级数。灰度级的级数i一般用2的指数m，即i=2m（m=1，2，…，8）当m=1时，灰度有黑白两级（二值化）当m=8时，则灰度有256级，0—255，0为黑，255为白。如MSS的量化级26TM、SPOT的量化级281、抽样：

进行像幅空间坐标的数字化，即沿像幅X遥感图像处理-基本原理-厉青课件灰度、坐标、物理图象和数字图象灰度、坐标、物理图象和数字图象讨论：

①采样间隔越小，数字图像与模拟图像越接近，图像数字化后图像的失真也就越小，但是采样间隔越小，图像数据量也就越大，这必然要影响遥感图像数字处理的速度。②灰度等级越多，对地物辐射光谱描述的越精确，但灰度等级的无限增加又会影响到图像的数据量，进而影响遥感图像数据的传输与处理，使遥感图像数字处理复杂化。讨论：3、数字图像

由离散化的坐标和灰度值组成的M×N数字矩阵：g（0，0），g(0,1),……,g(0,N-1)

g(x,y)＝g（1，0），g(1,1),……,g(1,N-1)

g（M-1，0），g(M-1,1),……,g(M-1,N-1)

即为数字化图像，其中每一个格网称为一像元，它在M×N数字阵中，用行、列号和灰度值表示。

3、数字图像

遥感图像处理-基本原理-厉青课件第三节遥感图像的数据结构

CCT磁带记录图像格式有三种：1、BIP（BandIntervalbyPixel）记录格式按象元波段交叉式记录2、BIL（BandIntervalbyLine）记录格式按行波段交叉式记录3、BSQ（BandSequential）记录格式按波段顺序式记录总之：BIP2：有利于作子区处理BSQ：便于单波段处理和提取BIL：适用于多波段运算和分类

第三节遥感图像的数据结构遥感图像处理-基本原理-厉青课件遥感图像处理-基本原理-厉青课件遥感图像处理-基本原理-厉青课件第四节遥感图像的基本统计分析

一、图像的统计量基本统计量是指一组数值或一幅单波段图像数据的中心趋势统计量和变化程度统计量。1、

均值：是一幅图像中所有像元亮度值的算术平均值2、

中值：是指图像中所有不同亮度值的中间值3、

众数：是图像中出现最多的一个亮度值

第四节遥感图像的基本统计分析4、

数值域：亮度值的动态变化范围，即最大和最小亮度值之间的差值gange=gmx-gmin5、

方差和标准差：

是由每个像元值Xi与均值的差异所累积形成的总的离散程度4、

数值域：亮度值的动态变化范围，即最大和最小亮度值之

方差和标准差描述了像元值与图像平均值的离散程度,是图像信息量大小的标志。方差是衡量一个波段中像元值差别的大小或分散的程度，一定程度上也是信息量多少的量度。方差小，亮度值比较集中，图像反差小；方差大，亮度值比较分散，图像反差大。6、反差：C1＝DNmax/DNminC2＝DNmax-DNminC3＝Sb（标准差）标准差大，反差大

方差和标准差描述了像元值与图像平均值的离散程度,是图二、

图像的直方图1、图像直方图的基本概念：直方图是图像中的每个波段亮度值的分布曲线。图像直方图的横坐标表示图像的灰度级变化，直方图的纵坐标表示图像中某个灰度级像元数目占整个图像象元数目的百分比或累计百分比。直方图是图像灰度分布的直观描述，它能够反映图像的信息量及分布特征。2、图像直方图的基本类型①频数直方图：纵坐标值是某个灰度级的像元在图像中出现的百分数，频数直方图相当于概率密度曲线。二、

图像的直方图遥感图像处理-基本原理-厉青课件遥感图像处理-基本原理-厉青课件遥感图像处理-基本原理-厉青课件②累积直方图：纵坐标是小于或等于特定灰度级像元在图像中的百分数，相当于累积概率密度曲线。②累积直方图：纵坐标是小于或等于特定灰度级像元在图像中的百3、正态分布由于图像数据具有随机性质，在图像像元数目相当大，而地物类型差异不很大的情况下，遥感数据像自然界的其它现象一样，接近正态分布。正态分布曲线：μ－均值σ－标准差3、正态分布4、直方图的偏斜与量度

实际图像的直方图总是和正态曲线之间有着不同程度的差异，这种差异称为直方图的偏斜。直方图发生偏斜时，遥感图像的均值、中值以及众数将明显地不一致。偏斜度公式SK可用下式计算：若SK为正值，称为正偏斜或向右偏斜；若SK为负值，称为负偏斜或向左偏斜。

4、直方图的偏斜与量度

三、多波段数据的统计特征多波段图像数据，不仅每个波段有其本身的统计特征，而且波段之间也存在着统计关系，其中最常采用的是协方差矩阵和相关矩阵。1、协方差与协方差矩阵：

三、多波段数据的统计特征2、相关系数与相关矩阵相关系数是描述多波段图像的相关程度的统计量。对遥感图像来说，二个波段之间彼此包含的信息很多，称相关程度高，若彼此包含不同信息，则说相关程度低。相关系数表示两个波段图像所包含的信息内容的重叠程度，是多波段图像彩色合成的重要依据。2、相关系数与相关矩阵

3、MSS图像的相关特点①两个波段之间的相关系数,与它们所在波谱中的距离有关。②r65对植被的反映最敏感，其次是r75。③相距最远的MSS4与MSS7的相关系数，说明两波段的信息内容差别最大。④多数情况下r76＞r65⑤r65、r75、r64、r74值的大小与植被发育程度成反比。⑥地形是影响相关程度的一个重要因素。地形起伏较大的地区，波段之间的相关系数普遍增大。⑦图像中含有较大面积水体时，r65、r75、r64显著降低，⑧时相的变化，也常常引起相关系数的变化。

3、MSS图像的相关特点

4、TM图像的相关特点（1）TM数据相关分组性①TM1、2、3相关系数较大，属于一组②TM5、7均为红外波段，相关系数较大，成为一组③TM4与其它波段的相关性都很小，单独成为一组④TM6因其几何分辨率较其它波段都低，无法与其它波段进行相关分析，也单独为一组。（2）波段相关系数大小受自然地理条件、气候、地质环境变化等因素的影响。

图像的增强处理在很大程度上是为了压缩或排除，不同波段之间的相关性和信息的重复。

4、TM图像的相关特点

第五节遥感图像的信息特征一、遥感图像的波谱特征信息与测度

波谱信息是地物识别的重要依据，遥感图像的波谱信息主要是通过灰度等级数目来度量的，灰度级数目越大，反映波谱信息越丰富。二、遥感图像的空间形态信息与测度

遥感图像的空间结构信息指的是亮度在空间上的变化特征，图像中的点、线、边缘、纹理都属于空间形态信息。影响空间信息的重要因素是遥感图像的几何分辨率三、遥感图像的时间特征信息与测度遥感图像的时间特征信息指的是不同时相遥感图像的光谱特征和空间结构的差异。遥感图像中的时间特征信息在目视解译时，是通过不同时相图像的对比分析来辨识的。卫星对地重复观测的周期越短、频率越高，则对地物光谱特征随时间变化的特征了解得就越多，因此，时间分辨率是衡量时间特征信息丰富程度的一种测度。

第五节遥感图像的信息特征

遥感图像特征：几何特征、物理特征和时间特征其表现参数为：空间分辨率、光谱分辨率、辐射分辨率和时间分辨。

图像的空间分辨率：指像元所代表的地面范围的大小，即扫描仪的瞬时视场，或地面物体能分辨的最小单元。

图像的波谱分辨率：指传感器在接收目标辐射的波谱时能分辨的最小波长间隔。间隔愈小，分辨率愈高。成像光谱仪在可见光至红外波长范围内，被分割成几百个窄波段，具有很高的光谱分辨率。

图像的辐射分辨率：指传感器接收波谱信号时，能分辨的最小的辐射度差。在遥感图像上表现为每一像元的辐射量化级。（由位数决定的记录辐射值的数值范围-，如8bit(0~255)。）

图像的时间分辨率：指对同一地点进行遥感采样的时间间隔，即采样的时间频率，也称重访周期。遥感图像特征：几何特征、物理特征和时间特征分辨率空间分辨率（Spatialresolution）:也称地面分辨率波谱分辨率（Spectralresolution）时间分辨率（Temporalresolution）辐射分辨率（Radiometricresolution）：由位数决定的记录辐射值的数值范围-，如8bit(0~255)。分辨率空间分辨率（Spatialresolution）:也图像的空间分辨率是指像素所代表的地面范围的大小，即扫描仪的瞬时视场，或地面物体能分辨的最小单元。图像的空间分辨率是指像素所代表的地面范围的大小，即扫描仪的瞬波谱分辨率是指传感器在接收目标辐射的波谱时能分辨的最小波长间隔。波谱分辨率是指传感器在接收目标辐射的波谱时能分辨的最小波长间Radiometric

Resolution8-bit(0-255)9-bit(0-511)10-bit(0-1023)000辐射分辨率是指传感器接收波谱信号时，能分辨的最小辐射度差。在遥感图像上表现为每一像元的辐射量化等级RadiometricResolution8-bit9-bTemporal

ResolutionJune1,2001June17,2002July3,2003RemoteSensorDataAcquisition16days时间分辨率指对同一地点进行遥感采样的时间间隔，即采样的时间频率，也称重访周期TemporalResolutionJune1,200第三章预处理第一节辐射校正、大气较正第二节几何校正的原理与方法第三节遥感图像数字镶嵌第三章预处理第一节辐射校正、大气校正

由于传感器响应特性和大气的吸收、散射以及其他随机因素影响，导致图像模糊失真，造成图像的分辨率和对比度相对下降，这些都需要通过辐射校正复原。

为什么要进行辐射纠正？1）传感器本身的特性，2）大气对于电磁辐射的衰减；（散射和吸收）3）地形因子的影响——阴影4)其它生态环境因子图象不能全部真实地反映不同地物地特征，影响了数字图象的质量。第一节辐射校正、大气校正一、系统辐射校正

（一）光学摄影机内部辐射误差校正光学摄影机内部辐射误差主要是由镜头中心和边缘的透射光的强度不一致造成的，它使得在图像上不同位置的同一类地物有不同的灰度值，如航空像片时常边缘较暗就是因此形成的。设原始图像灰度值为g，校正后的图像灰度值为g’则有g’＝g／cosθ式中，θ为像点成像时光线与主光轴夹角。遥感图像处理-基本原理-厉青课件

(二)光电扫描仪内部辐射误差的校正

光电扫描仪的内部辐射误差主要有二类：一类是光电转换误差；另一类是探测器增益变化引起的误差。通常采用楔校准处理方法加以消除。现在以陆地卫星可见光波段为例，校准模型：(二)光电扫描仪内部辐射误差的校正二、大气校正

所谓大气校正就是指消除主要由大气散射、吸收引起的辐射误差的处理过程。前提假设：大气散射的影响主要在短波波段，红外波段中清洁的水体几乎不受影响，反射率值应当为0。由于散射影响，而使得水体的反射率不等于0，推定是由于受到了天空辐射项的影响。二、大气校正(一)回归分析法由于大气散射主要影响短波部分，波长较长的波段几乎不受影响，因此可用其校正其他波段数据。作法如下：在不受大气影响的波段(如TM5)和待校正的某一波段图像中，选择由最亮至最暗的一系列目标，将每一目标的两个待比较的波段灰度值提取出来进行回归分析。例如，用TM1波段与TM5波段比较，作出回归直线图，(一)回归分析法回归方程为(3.9)(3.10)(3.11)(3.12)(3.9)(3.10)(3.11)(3.12)以上各式中T1、T5表示TM1与TM5波段灰度值，T’l为TM1波段校正后的灰度值，对于任一波段其校正量ai是第i波段回归分析的截距。实际校正时只要依据以上公式通过TM1、TM5波段灰度值计算出a1、b1值，再根据式(3.12)减去改正量即可求出每个像元TM1波段校正后灰度值。(二)直方图校正法另一种简单方法是通过灰度直方图对比找出校正量。如图所示，由于大气散射影响只作用于短波段，对可见光以外的红外波段几乎没有影响(如TM5波段)，因此如果影像范围内存在灰度值为零的地物(例如深水体，高山背阴处等)，其灰度直方图往往从原点开始，而其他波段灰度直方图离原点有一段距离(例如图中线段a1)，这段距离即为大气散射引起的灰度直方图漂移值。依此为改正量进行校正，相当于从每个像元灰度值中减去这个数值。其他波段的校正量可以由此类推求得。遥感图像处理-基本原理-厉青课件遥感图像处理-基本原理-厉青课件

(三)基于大气辐射传输模型方法(三)基于大气辐射传输模型方法三、噪声消除遥感图像的噪声源主要有大气传输信道的噪声及传感器内所产生的噪声。大气传输信道中由于大气的湍流扰动影响产生随机噪声；传感器的噪声源包含有转换和滤波过程中产生的噪声，例如，光电检测系统的电流不稳定性所表现的散粒噪声，即电流的无规则起伏现象，低光亮时为泊松分布，高光亮时为高斯分布。传感器的另一噪声源为滤波器电路中的电阻、电容产生的热噪声，这种噪声也具有高斯分布、零均值的随机特性；热红外波段传感器因各部分的温度变化也产生噪声效应；此外，摄影胶片记录信息时存在的胶片颗粒噪声，即由胶片上溴化银颗粒受等值曝光时因银颗粒的大小、形状不完全一致，而且分布不均匀引起的。以上这些噪音均属于对图像的高频干扰，因此可以采用滤波方法子以消除。三、噪声消除1、条纹消除遥感影像中条纹误差判定和消弱的常用方法有：平均值法、直方图法及在垂直扫描线方向上采用最近邻点法或三次褶积法。下图是采用三次褶积法进行条纹误差的效果图。1、条纹消除2、斑点消除校正后的斑点亮度值取其邻域像元亮度值的平均值或用三次褶积法进行修正。通常边缘附近的斑点不进行校正，影像四周的像元不进行斑点校正。2、斑点消除

第二节几何校正的原理与方法

遥感图像在获取过程中由于多种原因导致景物中目标物相对位置的坐标关系在图像中发生变化，这种变化称为几何畸变。纠正各种原因引起的畸变是几何校正的基本内容。一、遥感图像几何畸变来源遥感图像的变形误差总的可分为内部误差和外部误差两类。内部误差主要是由于传感器自身的性能、结构等因素造成；外部误差指的传感器以外的各因素所造成的，例如地球曲率、地形起伏、地球旋转等因素所引起的变形误差等。第二节几何校正的原

(一)传感器成像几何形态影响传感器一般的成像几何形态有中心投影、全景投影、斜距投影以及平行投影等几种不同类型。其中，全景和斜距投影产生图像变形规律可以通过与正射投影的图像相比较获得。1．全景投影变形红外机械扫描仪以及采用CCD直线阵列作为检测器的推帚式传感器的每一条扫描线都相当于中心投影，其成像面相当于全景缝隙摄影机的投影面，是一个圆柱面，称为全景面。因此，扫描视场角越大，边缘变形越大，图b示地面方格网成像后的变形状况，这种畸变又称全景畸变。2．斜距投影变形

侧视雷达属斜距投影，其成像变形规律如图c所示。

(一)传感器成像几何形态影响遥感图像处理-基本原理-厉青课件遥感图像处理-基本原理-厉青课件(二)传感器外方位元素变化畸变传感器外方位元素变化是指决定遥感平台姿态的6个自由度：三轴方向(X，Y，Z)及姿态角(φ,ω,κ)，其中任何一个发生变化，都会给遥感图像带来不同变形，这种畸变是成像瞬间的综合影响。对于不同类型的传感器，外方位元素变化带来的畸变可能不尽相同。下图以多光谱扫描图像为例说明这6个自由度变化对地面一个方格网图像成像后带来畸变的表现形式。

(二)传感器外方位元素变化畸变遥感图像处理-基本原理-厉青课件(三)地球自转的影响地球自转对于瞬时光学成像遥感方式没有影响，对于扫描成像则造成图像平行错动。如图所示：△ye＝tevφ其中，△ye为图像错动量；

te为扫描整景图像时间(te＝L／Rω)；

vφ是纬度为φ时该点地球自转线速度；L为像幅地面长度；R为地球平均半径6378km；

ω为卫星运行平均角速度。

遥感图像处理-基本原理-厉青课件(四)地球曲率影响地球曲率影响可以通过图说明，设OA0为成像基准面，A为地表一点。在考虑地球曲率影响情况下，A与OA0存在着由地球曲率引起的高差h，A在OA0代表的平面上投影点为A0，由于高差h的存在使得A点在像平面Fa0上产生像点位移。一般来说，在星下点视场角比较小，扫描范围又比较小时地球曲率影响可以忽略，此时可以看成近垂直投影。

(四)地球曲率影响二、几何校正原理遥感图像的数字纠正是通过计算机对图像每个像元逐个地解析纠正处理完成的，所以能够精确地改正线性和非线变形误差。它一般包括两方面：一是像元坐标变换，二是像元灰度值重新计算（重采样）。二、几何校正原理遥感图像的数字纠正是通过计算

1.地形图的准备

原则上要求所用地形图的比例尺应大于遥感影像制图的比例尺。对分辨率小于5m的影像制图，应采用1∶5万的地形图纠正；对于分辨率大于5m的影像制图，应采用1∶1万的地形图纠正。（一）、准备工作

1.地形图的准备（一）、准备工作

2.校正图像的准备根据影像数据分析与预处理的结果，首先需确定是否为多景数据处理。多景数据处理的原则为：时间相近的图像，可先镶嵌后再进行几何处理；获取时间差别较大的图像，应分别进行几何处理再镶嵌。

其次生成供选取控制点的图像。可以对图像进行增强以改善目视效果，有利于地物点的确定。也可以选择某一时相的TM彩色合成（743、543、741等）图像，作为供选取控制点的影像。2.校正图像的准备（二）、纠正变换函数的建立

用以建立影像坐标和地面坐标（或地图）间的数学关系，即输入图像与输出图像间的坐标变换关系。这种坐标变换关系，通常有两种互逆的表达式法：（二）、纠正变换函数的建立用以建立影像坐标

1.直接法方案

从原始图像阵列出发，按行列的顺序依次对每个原始像素点位求其在地面坐标系（也是输出图像坐标系）中的正确位置：

X＝Fx(x,y)Y＝Fy(x,y)式中Fx、Fy为直接纠正变换函数。按照原始图像的阵列，依次对每个象元(x，y)进行变换纠正，求得图像的位置（X，Y），同时把原图像(x，y)的灰度值送到新图像（X，Y）的位置上。

1.直接法方案

2.间接法方案

从空白的输出图像阵列出发，亦按行列的顺序依次对每个输出象元点位反求其在原始图像坐标的位置。x＝Gx(X,Y)y＝Gy(X,Y)式中Gx、Gy为间接纠正变换函数。同时把上式所算得的原始图像点位上的亮度值取出填回到空白输出图像点阵中相应的象元点位上去.由于计算的（x，y）不一定刚好位于原始图像的某个象素中心上，必须经过灰度内插确定（x，y）的灰度值。

2.间接法方案（三）、纠正后数字图像的边界范围

经过纠正后的图像仍为数字图像，它与原始图像的形状和方向都不一致，所以纠正变换前，必须为计算机输出图像预留一定的存储空间和该空间边界的地图坐标定义值，即必须预先确定纠正后数字图像的边界范围。其方法如下图：（三）、纠正后数字图像的边界范围遥感图像处理-基本原理-厉青课件根据X＝Fx(x,y)

Y＝Fy(x,y)

求出原始图像四个角点（a、b、c、d）纠正后的对应点（a′、b′、c′、d′）的大地坐标（Xa′,Ya′）,（Xb′,Yb′）,（Xc′,Yc′）,（Xd′,Yd′），然后求其最大和最小值，即：X1＝min(Xa′，Xb′,Xc′,Xd′)X2＝max(Xa′，Xb′,Xc′,Xd′)(3)Y1＝min(Ya′,Yb′,Yc′,Yd′)Y2＝max(Ya′,Yb′,Yc′,Yd′)由（3）式确定的大地坐标值即为纠正后图像的边界范围。根据X＝Fx(x,y)

为了把该边界范围转换为计算机中纠正图像的贮存数组空间，必须在其中划分出格网，每个网代表一个输出象元。为此应根据精度要求定义输出象元的地面尺寸△X和△Y，并以边界范围左上角A点为输出图像的坐标原点，以x′轴（AC边）表示图像行号，以y′轴（AB边）表示图像列号。图像总的行列数M和N由下式确定：

为了把该边界范围转换为计算机中纠正图像的贮存数输出图像在A－x′y′坐标系中，每个象元将以其所在的行列号来确定其位置。行列号的范围为：

x′＝1，2，……，My′＝1，2，……，N

任意一个输出图像坐标（行列号）与坐标的关系为：

输出图像在A－x′y′坐标系中，每个象元将以其所在的行列号或式中：Xp，Yp－纠正后象元中的地面坐标，xp′，yp′－纠正后象元中的图像坐标（行列号）或式中：Xp，Yp－纠正后象元中的地面坐标，（四）、纠正后数字图像灰度值的重采样

以间接法纠正方案为例，假如输出图像阵列中的任意一个象元在原始图像中的投影点的坐标值为整数时，便可简单地将整数点位上原始图像上的灰度值直接取出，填入输出图像中。但当投影点位的坐标不为整数时，则投影点的灰度值需根据周围阵列象元的灰度确定，这种方法就称为灰度值重采样。所谓重采样，是相对于遥感信息获取时已进行过一次采样而言的。

（四）、纠正后数字图像灰度值的重采样

1.邻近法：以离投影点最近的一个象元的灰度值作为输出象元的灰度值。这种方法保持了原来的亮度值不变，即光谱信息不变，但却搬动了空间位置，相比之下，几何精度差。但方法简单，计算速度快，灰度失真大。1.邻近法：2.双线性内插法：以投影点周围邻近四个象元的灰度值，确定输出象元的灰度值，其简单式为：2.双线性内插法：

gx′y′－输出象元的灰度值gi－邻近点I的灰度值pi－邻近点对投影点的权值（pi＝1/di，di为邻近点至投影点的距离）这种方法的计算量比邻近法要大，但几何上比较准确，即输出象元灰度值的保真度较好。

3.双三次卷积法：以投影点周围邻近16个象元的灰度值，确定输出象元的灰度值，其简单式为：

这种方法精度最高，但计算量较大这种方法精度最高，但计算量较大（五）、多项式纠正法：用多项式近似地描述纠正前后相应点的坐标关系，并利用控制点的图像坐标和参考坐标系中的理论坐标，按最小二乘法原理求解出多项式中的系数，然后以此多项式对图像进行几何纠正

（五）、多项式纠正法：其一般简化式为：Xi＝Fx(xi,yi)＝a0+a1xi+a2yi+a3xiyi+a4xi2+a5yi2Yi＝Fy(xi,yi)＝b0+b1xi+b2yi+b3xiyi+b4xi2+b5yi2或xi＝Gx(Xi,Yi)＝c0+c1xi+c2yi+c3xiyi+c4xi2+c5yi2yi＝Gy(Xi,Yi)＝d0+d1xi+d2yi+d3xiyi+d4xi2+d5yi2其中：X，Y：为控制点在参考坐标系中的理论坐标；x,y：为同名控制点对应的原始图像坐标（行列号）ai、bi、ci、di：为多项式系数。

利用有限的控制点的已知坐标，解求多项式的系数，确定变换函数。然后将各个像元带入多项式进行计算，得到纠正后的坐标。其一般简化式为：利用有限的控制点的已知坐标实际工作中多项式最高次数一般取2就足以满足精度要求了。可以看出，只要求解出ai、bi、ci、di这四组系数就可以直接利用上式进行纠正了。实际工作中常利用待定系数法求解ai、bi、ci、di。选择待纠正图像和地形图（或其它参考地图、图像）上均可清晰辨认像点作为地面控制点（groundcontrolpoint,GCP）读取控制点在图像和地形图（或其它参考地图、图像）上坐标，代替上式中的x,y,X,Y,得出一个关于ai和bi、ci和di方程组。可以看出，控制点个数不得少于方程中ai、bi、ci、di的数目i。这样通过求解ai和bi、ci和di方程组，可确定上式系数。多项式系数求出后，根据上式可以求解原始图像任一像元的坐标，并对图像灰度进行内插，获取某种投影的纠正图像。

为此必须选择控制点

实际工作中多项式最高次数一般取2就足以满足精度要求了。可选择控制点时，应遵循以下原则：

①均匀分布：一般先在图像的四角和对角线交点处选择控制点，然后逐渐加密，保证均匀分布。②特征明显：尽可能选在固定的地物交叉点上，无精确定位的标志情况下，利用半固定的地形地物交叉点（山顶、河流交叉处）。如道路交叉点，标志物，水域的边界，山顶，小岛中心，机场等。③足够数量：控制点数量每景宜在25～35个左右，山区或丘陵区适当增加。数量应当超过多项式系数的个数（(n+1)*(n+2)/2）。当控制点的个数超过多项式的系数个数时，采用最小2乘法进行系数的确定，使得到的系数最佳。

ENVI:>(degree+1)^2选择控制点时，应遵循以下原则：ENVI:>(deg

采用多项式纠正法应注意以下几点：第一，多项式纠正法的精度与地面控制点(GCP)的精度、分布、数量及纠正范围有关。GCP的位置精度越高，则几何纠正精度越高。对于一般齐次多项式，GCP的个数至少不得低于多项式的系数个数。适当增加GCP的数量，可以提高几何纠正的精度，但过多地增加GCP的数量，不仅不会显著提高纠正精度，而且会增大选择GCP的工作量，有时甚至难以选出大量的GCP。一幅TM影像要纠正到一个像元以内的精度，一般需要30个左右均匀分布的GCP；一幅SPOT卫星影像采用三次多项式拟合，用13—17个GCP，其纠正精度可达两个像元左右。·第二，GCP应尽可能在整幅图像内均匀分布，否则会在GCP密集区几何纠正精度较高，而在GCP分布稀疏区将出现较大的拟合误差。

采用多项式纠正法应注意以下几点：

第三节遥感图像数字镶嵌

实际工作中时常需要将多幅卫片或相邻航片镶嵌在一起，早年通过裁切光学像片拼接在一起，误差自然很难避免。现在通过数字镶嵌方法可以较好地实现“无缝”拼接。数字镶嵌在理论和方法上与几何校正类似，有几点需要注意：第一，镶嵌要有足够宽的重叠区，最好不少于图像的l／5如果过于狭窄会影响镶嵌精度，特别是图像边缘会出现扭曲。第二，镶嵌时相邻的图像往往色调(对彩色图像)或灰度值不一致(对灰度图像)，此时需要其中一幅图像进行诸如“直方图匹配”等处理。第三，如果镶嵌后需要某种地图投影变换，最好依据该地图投影方式先分幅校正，然后再镶嵌，这样精度会更高些。航空像片扫描进计算机后也可以进行数字镶嵌，但由于中心投影方式(特别是大比例尺航片)带来的边缘变形导致镶嵌很难满足精度要求，此时可考虑先进行正射投影变换再进行镶嵌第三节遥感图像数字镶嵌准备工作：1）

要参加拚接的图像必须具有统一的坐标系，即首先进行图像的几何纠正。2）

图像灰度的调整对于彩色图像，需要从红绿蓝三个波段分别进行灰度的调整；对于多波段的图像文件，进行一一对应的多个波段的灰度调整。灰度调整的方法：进行交互式的图像拉伸，进行图像直方图的规定化(匹配)，或者进行更加复杂的类似变化（如线性回归）。准备工作：遥感图像处理-基本原理-厉青课件确定输出图像的大小；确定参加图像拚接的方法；有无羽化（Feathering）。

参加拚接的图像如没有羽化，则用上面的图像灰度值直接替代重复部分的灰度值；如果有羽化，则重复部分的灰度值由两个图像的灰度值确定。羽化方法：一是边缘羽化，二是用户自定义羽化。

图像拚接的操作：确定输出图像的大小；图像拚接的操作：1）边缘羽化用户需要键入羽化的距离。羽化部分的灰度值采取两个图像的线性加权值。边缘处Image#20%,Image#1100%,AveragingDistance处，Image#2占100%,Image＃1占0％；中间部位线性组合。1）边缘羽化羽化部分的灰度值采取两个图像的线性加权值。2)自定义羽化首先必须在图像上使用ANN定义Cutline.灰度值的计算方法：Cutline处，Image#2占100%,Image#1占0%,AveragingDistance处，Image#2占0%,Image＃1占100％；中间地带，线性插值。2)自定义羽化首先必须在图像上使用ANN定义Cutline第四章图像增强Imageenhancement第四章图像增强图像增强Imageenhancement目的：突出图像中的有用信息，扩大不同影像特征（亮度值）之间的差别，以提高对图像的解译和分析能力，使之更适合实际应用。类型划分：按照主要增强的信息内容可分为：

波谱特征增强：突出不同地物之间波谱特征的差别空间特征增强：突出空间形态特征、边缘、线条、纹理结构特点时间信息增强：提取多时相图像中波谱与空间特征随时间变化的信息

按照参加处理运算的波段影像数目可分为：

单波段图像增强：对单一波段数据处理，如反差扩展多波段图像增强：对多个波段数据处理，如彩色合成图像增强Imageenhancement

方法：从图像处理的数学形式看：点处理（点运算）邻域处理（局部运算）点处理：点运算公式：g’(x,y)=I［g(x,y)］I为与(x,y)位置无关的线性或非线性变换。点运算就是将原图像中每一个像元的灰度值，通过变换转换成输出图像中一个新的灰度值，其图像中各个像元的位置并不改变。

方法：从图像处理的数学形式看：

邻域处理：局部运算形式：g’(x,y)=g(x,y)*h(x,y)h(x,y)为某种增强处理算子，符号＊为卷积运算。

局部运算是针对一个像元点周围的一个小领域的所有像元而进行，输出值的大小以像元点在原图像中的灰度值大小有关外，还决定于它邻近像元点的灰度值大小。

邻域处理：局部运算形式：g’(x,y)=g(x,y)*第五章

图像分类

第五章

图像分类

第一节图像分类的基本原理

一、图像分类的概念图像分类是将图像的所有像元按其性质分为若干个类别的技术过程。分类的实质是把多维特征空间划分为若干区域（子空间），每个区域相当于一类（依据各类地物的光谱特征分析来选择特征参数），然后将影像内各个像元划分到各个子空间中区，从而实现分类。以区别图象中所含的多个目标物为目的，对每个像元或者比较均匀的像元组给出对应其特征的名称。

第一节图像分类的基本原理一、图像分类的概念特征名称一般称类别，例如土地利用/覆盖中的居民地、水域、园地、林地、耕地等。参与分类的多个特征量（向量）所定义的空间，称为特征空间，例如2个波段定义的2维特征空间，3个波段定义的3维特征空间等。特征名称一般称类别，例如土地利用/覆盖中的居民地、水域、园

目视判读：以图像的空间特征（地物的几何特征和光谱特征的空间反映）为判读依据。

计算机判读：以图像像元的灰度（地物光谱特征的直接反映）为判读依据，即通过对图像灰度级的统计、运算、对比和归纳来实现地物的判读和分类。目视判读：以图像的空间特征（地物的几何特征和光谱特计算机判读是建立在以下思想基础上实现的：①同类地物具有相同（似）的光谱特征，不同的地物其光谱特征具有明显的差别。由于影响地物光谱特征因素很多，所以影像的判读都是建立在统计分析的基础上的。②同类地物的图像灰度概率在单波段（称一维空间）符合正态分布规律。③多维图像（多波段）中的一个像元值（灰度）向量，在几何上相当于多维空间中的一个点，而同类地物的像元值，既不集中于一点，也绝非是杂乱无章的分布，而是相对地密集在一起形成一个点群（一个点群就是地物的一种类别）。计算机判读是建立在以下思想基础上实现的：二、地物点群的分布1、理想情况不同类别地物的集群在某一个子空间（即某一子波段图像中）的投影（即亮度范围或一维直方图分布范围）是完全可以区分开的。这时只要在单波段图像中利用影像亮度分割的措施便可以达到地物分类的目的。2、典型情况不同类别地物的点群在任一子空间中的投影都是相互重叠的，但在总的特征空间中都是完全可以区分的。这就意味着利用单波段图像不能实现完善的分类，只有利用多光谱图像，才有可能精确地分类。二、地物点群的分布3、一般情况无论在总空间或子空间中，不同类别点群都有重叠现象。这时重叠部分的影像在分类时，会出现一定程度的误差，这是遥感图像分类中最常见的情况。

3、一般情况无论哪种情况，分类问题的核心都是确定地物类别点群在特征空间的中心位置、边界范围及类与类之间的边界。

位置（一个点群的中心）－计算图像灰度的均值向量（即数学期望）

范围－计算图像灰度的标准差向量（或协方差矩阵)。即点群的离散程度。

边界－应用边界函数，鉴别图像像元的类别归属。无论哪种情况，分类问题的核心都是确定地物类别点三、监督分类和非监督分类1、监督分类

先在图像中选择有代表性的训练区（应是波谱特征比较均一的地区），由训练组数据得出各个类别的统计数据，然后根据这些统计数据对整个图像进行分类，即把图像中各个像元点归化到给定类中。2、非监督分类根据图像数据的本身统计特征及自然点群的分布情况划分类别，各类别的统计特征是未知的，一般是提供少数阈值对分类过程给以部分控制。

监督分类与非监督分类的最大区别在于，监督分类首先给定类别，而非监督分类由图像数据的统计特征来决定。

三、监督分类和非监督分类

监督分类和非监督分类相同之处：它们都是以图像的灰度为基础，通过统计、计算灰度的均值、方差、协方差(标准差)(或协方差矩阵)等特征参数进行分类的。为了提高分类精度，通常还采用一些“辅助信息”参入计算机的自动分类识别.例如，由于遥感图像时间特征的影响，考虑采用多时相的图像数据。由于地形的影响考虑采用数字地面模型的数据。在具有地理信息系统的情况下，采用地理信息系统支持的分类等。

监督分类和非监督分类相同之处：它们都是以图像的灰度为基谢谢大家！谢谢大家！二、遥感图像处理基本原理第一章概述第二章遥感图像处理基本原理第三章遥感图像预处理第四章图像增强第五章

图像分类二、遥感图像处理基本原理第一章概述第一章概述一、遥感图像处理技术的概况和发展二、遥感数字图像处理的主要内容、特点及应用三、数字图像处理设备和系统第一章概述一、遥感图像处理技术的概况和发展第一章概述

一、遥感图像处理技术的概况和发展图像处理技术是随着计算机技术的发展而开拓出来的一个新的应用领域，就是把图像转换成一个数据矩阵存放进计算机对之进行处理。具体说，利用数字计算机或其它数字硬件，对从图像信息转换而来的电信号进行某些数学运算，以期提高图像的实用性，从而达到人们所要求的某些预期的结果。第一章概述

数字图像处理起源于20世纪20年代，首先用于改变数字化报纸的图片。新闻照片经过数字化，由海底电缆开始在纽约和伦敦之间传送，为了用电缆传输图片，首先要进行编码，然后在接收端由专用打印设备进行图片重建。此时的技术能力，仅把图像编码为五个亮度级，到1929年才提高到15个等级。但那时图像的处理还不是用计算机进行的。从六十年代起，随着电子计算机技术的进步，数字图像处理获得了飞跃的发展。1964年美国喷气推进实验室（JPL）开始用计算机处理航天照片，以消除畸变，提高质量，取得了惊人的效果。从此以后，遥感数字图像处理技术和设备得到了迅速的发展。

硬件：主处理机：①60年代到70年代上半期，遥感图像处理的实现一般都是在已有的大、中型计算机上配备必要的外围设备，利用图像处理软件来完成图像处理与分析功能。②70年代下半期到80年代初期，专用图像处理机的出现，将许多图像处理算法固化，一部分图像处理与分析功能可以用硬件快速实现，从而使遥感图像处理速度大幅度地提高。③80年代初到90年代初期，超小型机在图像处理中占有统治地位，图像运算速度也有了很大的提高。④近年来，在遥感图像处理系统的构成方面，硬件发展已进入了第四个阶段，即在工作站上建立大规模专业化的遥感图像处理系统，如建立在工作站上的ERMAPPER。而一般的应用用户，又可以建立以PC机为核心遥感图像处理系统，专用图像处理机在遥感图像处理系统中地位已经开始下降。软件：软件功能反映了遥感技术应用水平①遥感图像处理系统与GIS的融合②多种遥感资料的处理。③非遥感资料的图像处理。硬件：二、遥感数字图像处理的主要内容、特点及应用

1、数字图像处理的内容

a、图像数字化及编码压缩Imagedigitization.Coding.Compressionb、图像恢复或校正Imagerestorationorcorrectionc、

图像增强