遥感数字图像处理重点整理.pdf

《遥感数字图像处理重点整理.pdf》由会员分享，可在线阅读，更多相关《遥感数字图像处理重点整理.pdf（23页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、遥感数字图像处理遥感数字图像处理第一章概论第一章概论1、图像图像定义：IMAGE，指通过镜头等设备得到的视觉形象或以某一技术手段再现于二维画面上的视觉信息，是二维数据阵列的光学模拟。分类：按人眼的视觉可视性：可见图像可见图像照片、素描、油画不可见图像不可见图像不可见光成像如紫外线、红外线、不可见测量值如温度、人口密度等的分布图按图像的敏感程度和空间坐标的连续性：数字图像数字图像指用电脑存储和处理的图像，是一种空间坐标和灰度不连续、以离散数字原理表达的图像，不见图不见图像像模拟图像模拟图像又称光学图像光学图像，指空间坐标的明暗程度连续变化的、电脑无法直接处理的图像，可见图像可见图像模拟图像数字图

2、像：模/数转换A/D 转换数字图像模拟图像：数/模转换D/A转换2、像素像素定义：是 A/D 转换的取样点，是电脑图像处理的最小单元，每个像素具有特定的空间位置空间位置和属性特征属性特征。3、遥感数字图像遥感数字图像定义：数字形式的遥感图像，不同的地物能够反射或辐射不同波长的电磁波，利用这种特性，遥感系统可以产生不同的遥感数字图像。像素值：称为亮度值或灰度值、DN 值，量化的整数灰度就是数字量值。亮度值的高低由遥感传感器所探测到的地物电磁波的辐射强度决定。具有相对应的意义，仅在图像内才能进行相互比较。遥感数字图像与照片的差异遥感数字图像与照片的差异照片来自于模拟方式通过摄影系统产生没有像素没有

3、行列结构没有扫描行0 表示没有数据任何点都没有编号摄影受电磁光谱的成像范围限制一旦获取了照片，颜色就是确定的具有红、绿、蓝 3 个通道遥感数字图像来自于数字方式通过扫面和数码照相机产生基本构成单位是像素具有行和列可能会观察到扫描行0 是数值，不表示没有数据每个点都有确定的数字编号可以是电磁光谱的任意范围颜色没有特定的规则，在处理过程中可以根据需要通过合成产生多个波段3-80004、遥感数字图像处理遥感数字图像处理定义：是利用电脑图像处理系统对遥感图像中的像素进行系统操作的过程。传统的模拟图像受媒介大小的限制无法完全表述这些信息，也很难进行信息的进一步处理，只有经数字化后才能有效地进行信息分析和

4、处理，数字图像处理极大地提高了图像处理的精度和信息提取的效率。主要内容：图像增强图像增强：压抑、去除噪声，增强整体图像或突出图像中的特定地物的信息，使图像更容易理解、解释和判读。灰度拉升、平滑、锐化、彩色合成、主成分K-L变换、K-T 变换、代数运算、图像融合等增强过程本身不会增加数据中原有的信息内容，仅仅是突出了特定的图像特征，使得图像更易于可视化的解释和理解。图像校正：图像校正：也称为图像恢复、图像复原，主要是对传感器或环境造成的退化图像进行模糊消除、噪声消除、几何失真或非线性校正。校正方法：辐射校正和几何纠正。信息提取：信息提取：根据地物光谱特征和几何特征，确定不同地物信息的提取规则。主

5、要包括图像分割、分类等方法，处理结果为分类专题图。5、遥感数字图像处理系统遥感数字图像处理系统硬件系统硬件系统：电脑、数字化设备、大容量存储器、显示器、输出设备、操作台软件系统软件系统：ERDAS IMAGINEC 语言、ENVIIDL 语言、PCI Geomatica、ER Mapper6、数字图像处理的发展和两个观点两个观点离散方法离散方法：一幅图像的存储和表示均为数字形式，因为数字是离散的空间域空间域图像处理以图像平面本身为参考，直接对图像中的像素进行处理。连续方法连续方法：图像通常源自物理世界，服从可用连续数学描述的规律，因此具有连续性。频率域频率域基于傅里叶变换，频率域图像处理时对傅

6、里叶变换后产生的反映频率信息的图像进行处理。第二章遥感数字图像的获取和存储第二章遥感数字图像的获取和存储1、遥感系统遥感系统：是一个从地面到空中乃至整个空间，从信息收集、存储、传输、处理到分析、判读、应用的技术体系，主要包括遥感试验遥感试验、信息获取信息获取传感器、遥感平台、信息传输信息传输、信息处理信息处理、信息应用信息应用等 5个部分。2、传感器传感器信息获取的核心定义：收集和记录电磁辐射能量信息的装置，是信息获取的核心部件，如航空摄影机航空摄影机、多光谱扫描仪多光谱扫描仪、成像成像等分类：按工作方式是否具有人工辐射源：被动方式被动方式被动遥感：以太阳辐射和地物自然辐射为辐射源，不需人工辐

7、射源，如各种摄像机、扫描仪、辐射计等，其工作波段集中在可见光和红外区。主动方式主动方式主动遥感：具有人工辐射源，主动向目标发射强大的电磁波，然后传感器接受目标反射的回波，如各种形式的雷达，其工作波段集中在微波区。按数据的记录方式：成像方式成像方式：成像传感器接受的目标电磁辐射信号转换成数字或模拟图像，是目前最常见的传感器类型。按成像原理可以分为摄影成像和扫描成像两类。摄影成像摄影成像：在打开快门后的一瞬间几乎同时收集目标上所有的反射光，聚焦到胶片上成为一幅影像，并记录下来如全景摄影机、多光谱摄影机、缝隙摄影机、框幅摄影机等。扫描成像扫描成像：扫描方式的传感器逐点逐行地收集信息，各点的信息按一定

8、顺序先后进入传感器，经过一段时间后才能收集完一幅图像的全部信息。目标面扫描目标面扫描的方式：收集系统直接对目标面扫描，一点一行顺序收集目标面上各单位的信息，然后拼成一幅图像。光学-机械扫描仪多光谱扫描仪MSS 专题制图仪TM成像雷达全景雷达和侧视雷达影像面扫描影像面扫描的方式：收集系统不直接对地面扫描，而是先用光学系统将目标的辐射信息在靶面上聚集形成一幅图像然后利用摄像管中的电子束对靶面扫描来收集其数据，或依靠电荷耦合器件CCD组成的阵列进行电子扫描来获得数据。电视摄影机、固体扫描仪非成像方式非成像方式：记录的是一些物理参数，不能产生图像，如可见光-近红外辐射计、热红外辐射计、微波辐射计、微波

9、高度计、微波散射计3、电磁波电磁波传感器紫外可见光红外微波特点对紫外线光感测用感光胶片或光电探测器作为感测元件，地面分辨率高，但只能在晴朗的白天使用近红外常使用感光胶片，中、远红外常使用光学机械扫描仪；具有昼夜工作能力具有昼夜工作能力，能提供高分辨率图像，并传感器实例紫外摄影机可见光摄影机、TV 摄影机等红外光学扫描仪合成孔径雷达波长范围/nm503803807607601.0*1061.0*1061.0*109与日照、云层遮挡无光，有一定穿透能力4、传感器的分辨率传感器的分辨率定义：指传感器区分自然特征相似或光谱特征相似的相邻地物的能力。高分辨率意味着区分能力强，能够区分小的相邻地物，低分辨

10、率意味着能够获取大范围的平均辐照度。1辐射分辨率辐射分辨率：是传感器区分反射或发射的电磁波辐射强度差异的能力高辐射分辨率意味着可以区分信号强度的微小差异。表示：可见光、近红外波段：用等效反射率表示热红外波段：用噪声等效温差、最小可探测温差和最小可分辨温差表示。遥感图像：图像的量化位数可以看做是辐射分辨率的近似表述。2光谱分辨率光谱分辨率：是传感器记录的电磁光谱中特定波长的范围和数量波长范围越窄，波段数越多，光谱分辨率越高，地物越容易被区分和识别。但并非波段越多越好，会产生“掩盖”地物的辐射特性，不利于快速探测和识别地物。3空间分辨率空间分辨率：是遥感图像上能够详细区分的最小单元的尺寸或大小，即

11、传感器能把两个目标物作为清晰的实体记录下来的两个目标物之间最小的距离。是表征图像分辨地物目标细节能力的指标。也可以用地面分辨率 像素所代表的地面实际尺寸的大小 来表示。表示：像素：是将地物信息离散化而形成的格网单元，在遥感图像中，单位为米，像素为正方形，其大小与遥感空间分辨率高低密切相关，空间分辨率越高，像素越小。解像力解相率：用单位距离能分辨的线宽或间隔相等的平行线的条数来表示，单位为线/毫米或线对/毫米瞬间视角场 IFOV角分辨率：是传感器的瞬间视场，它和遥感平台高度H 共同决定了地面的分辨单元，即像素大小 P，P=2Htan(IFOV/2)分类：高高空间分辨率：小于10M，SPOTQUI

12、CKBIRDIKONOS，重访周期长数天，能够反映地物明确的几何信息，适用于特定地区进行定点监测，多用于数字城市和工程制图。中中空间分辨率：10100M，ASTERTM，重访周期数周，具有较多的光谱信息，便于进行土地利用和土地覆盖、资源、地表景观等方面的研究。低低空间分辨率：大于 100M，NOAAMODIS，重访周期短数小时，适用于进行大范围的环境遥感监测，例如洪水、火灾、云和沙尘暴等。4时间分辨率时间分辨率：对同一目标进行重复探测时，相邻两次探测的时间间隔一类是传感器本身设计的时间，受卫星运行规律影响，不能改变一类是人为设计的时间分辨率，它一定等于或小于卫星传感器本身的时间分辨率在遥感图像

13、应用中意义重大，可以进行动态监测和预报，进行自然历史便器和驱动力分析，可以提高成像率和解像力，提高地物的识别精度。5、数字化包括采样、量化采样、量化采样采样：将空间上连续的图像变换成离散点即像素的操作。采样间隔和采样孔径的大小是重要的参数。采样间隔影响着图像表示地物的真实性，间隔越小，图像越接近于真实，但采用成本及后处理的成本也越高，图像存储所需要的空间越大。采样时，连续的图像空间被划分为网格，并对各个网格内的辐射值进行测量。量化量化：将像素灰度值转换成整数灰度级的过程。图像数据 M 行，N 列，量化位数为 g，存储空间 M*N*g/8 字节。2=M，其中 N 为量化级，M 为量化位数量化灰度

14、级通常得到的都是量化之后的遥感图像，图像中的像素值被称为数字值DN 值、灰度级或亮度值。量化影响着图像细节的再现程度，间隔越大，细节损失越多，图像的棋盘化效应越明显，量化影响着图像细节的可分辨程度，量化位数越高，细节的可分辨率越高，保持图像大小不变，降低量化位数减少了灰度级会导致假的轮廓。6、遥感图像类型遥感图像类型不相干图像不相干图像：光学遥感所产生的图像，包括多光谱图像多波段图像，每个采样位置包括多个波段的值、高光谱图像和高空间分辨率图像属于被动遥感，图像受大气状况影响很大，限制了再多云多雨地区的应用。相干图像相干图像：微波遥感所产生的图像，多取决于传播的条件，微波遥感属于主动遥感，其穿透

15、能力强，不受天气的影响，可以全天候工作。7、遥感数字图像的级别遥感数字图像的级别0 级产品：未经过任何校正的原始图像数据。1 级产品：经过了初步辐射校正的图像数据。2 级产品：经过了系统级的几何校正我们到手的数据3 级产品：经过了几何精校正，保密级别8、元数据元数据定义：是关于图像数据特征的表述，是关于数据的数据。包括图像获取的日期和时间、投影参数、几何校正精度、图像分辨率、辐射校正参数元数据与图像数据同时发布，或者嵌入到图像文件中，或者是单独的文件。也称为头文件，多为二进制格式或随机。9、遥感图像数据格式遥感图像数据格式设图像数据位 N 列，M 行，K 个波段1BSQBSQ 格式格式是像素按

16、波段顺序依次排列的数据格式，先按照波段顺序分块排列，在每个波段块内，再按照行列顺序排列，保证了像素空间位置的连续性。B1(1,1)(1,2)(1,3)(1,N)(2,1)(2,2)(2,3)(2,N)BK(1,1)(1,2)(1,3)(1,N)(2,1)(2,2)(2,3)(2,N)(M,1)(M,2)(M,3)(M,N)2BILBIL 格式格式像素先以行为单位分块，在每个块内，按照波段顺序排列像素，像素的空间位置在列的方向上是连续的。第一行B1(1,1)(1,2)(1,3)(1,N)B2(1,1)(1,2)(1,3)(1,N)BK(1,1)(1,2)(1,3)(1,N)第 M 行B1(M,1

17、)(M,2)(M,3)(M,N)B2(M,1)(M,2)(M,3)(M,N)BK(M,1)(M,2)(M,3)(M,N)3BIPBIP 格式格式以像素为核心，像素的各个波段数据保存在一起，打破了像素空间位置的连续性。保持行的顺序不变，在列的方向上按列分块，每个块内为当前像素不同波段的像素值。B1B2 BKB1B2 BK第一行(1,1)(1,1)(1,1)(1,2)(1,2)(1,2)第二行(2,1)(2,1)(2,1)(2,2)(2,2)(2,2)第 M 行(M,1)(M,1)(M,1)(M,2)(M,2)(M,2)4TIFFTIFF 图像格式图像格式5HDFHDF 数据格式数据格式优势：独立

18、于操作平台的可移植性；超文本；自我描述性；可扩展性6GeoTIFFGeoTIFF 图像格式图像格式独立性和拓展性的特点，条理清晰、结构严谨，容易实现与其他遥感影像格式的转换。7陆地资源卫星陆地资源卫星 L5L5 的数据格式的数据格式10、图像文件图像文件坐标：左上角像素的坐标从0 开始，向右向下按整数递增，经过集合校正后图像文件坐标用地图坐标来表示大小=图像行数*图像列数*每个像素的字节数*波段数*辅助参数一般为1 8 位量化位 1 个字节，16 位为 2 个11、数字图像分辨率数字图像分辨率定义：图像上的点被映射或指定到给定的空间里的数量通常是以英寸、厘米、像素为单位，是图像中最小可分辨距离

19、。第三章遥感数字图像的表示和统计描述第三章遥感数字图像的表示和统计描述1 1、遥感图像模型、遥感图像模型定义：遥感图像是传感器通过探测地物电磁波辐射能量所得到的的图像，反映连续变化的物理场，虽然波段不同，记录的辐射能量、成像的方式以及成像系统等也有差异，但还是可以从理论角度归纳得到一个具有普遍意义的模型，成为遥感图像模型。图像上的能量图像上的能量=目标发射的能量目标发射的能量+目标反射的能量目标反射的能量陆地遥感图像模型：陆地遥感图像模型：L(x,y,t,p)=1 (,p)()+(,p)(,)其中，(,p)为目标物的反射率；()为黑体的电磁波发射能力；(,)为入射的辐射量；p表示极化方向，主要

20、用于微波成像；代表波长；为成像时间。在可见光和近红外波段，白天物体自身发射的辐射量课忽略不计，可以简化为L(x,y,t,p)=(,p)(,)2 2、多源图像、多源图像定义：在同一地区，随时间、波段和极化方向不同而获得的多个图像的组合，叫做多源图像分类：多波段图像多光谱图像多时相图像多极化图像3、图像函数 f(x,y)实际上代表在二维空间内物体反射或发射电磁波辐射能量的分布，不是传感器实际记录的图像数据。设g(x,y)表示二维空间的图像函数，则对应的变换关系可以表示为g(x,y)=Tf(x,y)，其中T 表示某种由地物到图像的变换；g(x,y)是遥感图像处理后产生的图像函数。g(x,y)具有连续

21、性；定义域的限定性；函数值的限定性；函数值物理意义的明确性。4 4、遥感图像的数字表示、遥感图像的数字表示1确定的确定的写出图像函数表达式，对于数字图像，表示成矩阵或向量图像的矩阵表示二值图像：每个像素的取值为 0 或 1 的图像，没有颜色的概念，数值仅包括 0 和 1，是逻辑运算之后的结果，0 用来表示背景假，1 用来表示前景目标真。灰度图像：每个像素由一个量化的灰度值灰度级来描述的图像，单波段图像为灰度图像。对于8 位量化而言，灰度值 0 为黑色，255 为白色。彩色图像：每个像素由红、绿、蓝RGB三原色构成的图像，对于多光谱遥感图像，可通过 RGB 合成产生彩色图像。图像的向量表示按行或

22、列排序，优点是可以直接利用向量分析的有关理论和方法。2统计的统计的用一种平均特征来表示图像的数字特征可作为区分或识别图像中地物的依据，将图像看做具有正态分布的随机变量，可以使用统计学方法对图像进行统计描述。单波段图像的统计特征单波段图像的统计特征基本统计特征反映像素值平均信息的统计参数：均值、中值、众数反映像素值变化信息的统计参数：方差、变差最大值与最小值之差、反差比照度直方图：是灰度级的函数，描述的是图像中各个灰度级像素的个数。对于数字图像来说，直方图实际就是灰度值概率密度函数的离散化图形。H(i)=(=0,1,L 1)性质：反映了图像灰度的分布规律可以通过修改图像的直方图来改变图像的反差。

23、任何一幅特定的图像都有唯一的直方图与之对应，不同的图像可以有相同的直方图。如果一幅图像仅包括两个不相连的区域，并且每个区域的直方图已知，则整幅图像的直方图是这两个区域的直方图之和遥感图像的数据服从或接近正态分布，直方图的形态与正态分布的曲线形态类似应用：可以有目的地改变直方图形态来改善图像的比照度。峰值偏向灰度大，图像偏亮；反之二阶矩表示灰度级的比照度；三阶矩表示直方图的偏斜度；四阶矩表示直方图的峰度多波段图像的统计特征多波段图像的统计特征协方差：协方差矩阵P48值越高表示两个波段图像之间的协变性越强相关系数：相关矩阵P48直方图匹配：可用欧氏距离来衡量5、窗口、邻域和卷积窗口窗口：对于图像中

24、任一像素x,y，以此为中心，按上下左右对称所设定的像素范围，称为窗口邻域邻域：中心像素周围的行列称为该像素的邻域，按照与中心像素相邻的行列总数来命名邻域运算邻域运算：对于中心像素x,y，其值 f(x,y)表示，可按照相邻性规则通过计算产生卷积卷积：是空间域上针对特定窗口进行的运算，是图像平滑、锐化中使用的基本的计算方法卷积公式g(i,j)=或卷积核 是相邻像素对中心像素影响=1=1(,)(,)h(x,y)为窗口模板程度的表述，边缘的处理方法：设为0 值；保留原值；按对称原则从图像中取值6、滤波滤波狭义：滤波是指改变信号中各个频率分量的相对大小，或者分立出来加以抑制，甚至全部滤除某些频率分量的过

25、程广义：把某种信号处理成为另一种信号的过程Y(f)=X(f)H(f)，滤波的概念主要用在频率域中，在空间域中即为卷积运算。7、纹理纹理定义：图像的某些局部性质，或是对局部区域中像素之间关系的一种度量。纹理基元按照某种确定性的规律或只是按照某种统计规律重复排列组成的。三要素：某种局部的序列性，不断重复；序列基本部分非随机排列组成；各部分大致是均匀的统一体，纹理区域内任何地方都有大致相同的结构尺寸分类：人工纹理有规则的，是自然背景上的符号排列组成的，可以使线条、点、字母、数字自然纹理不规则的，是具有重复排列现象的自然景物，如森林、草地之类的纹理在局部区域内呈现不规则形，而在宏观上又表现出某种规律。

26、纹理的性质有均匀性、密度、粗细度、粗糙度、规律性、线性度、定向性、方向性、频率、相位。传统的纹理特征描述方法：统计方法利用图像的特征求出特征值，实现对纹理特征的描述，自相关函数、灰度共生矩阵、灰度级行程长、滤波模板、随机模型、分形模型结构方法用结构方法力图找到纹理基元，一般只适用于规则性较强的人工纹理8、空间自相关函数方法空间自相关函数可用来对纹理的粗糙程度进行描述。一般地，粗纹理的自相关函数随距离的变化比较缓慢，细纹理变化比较快9、共生矩阵是用来描述纹理中灰度基元之间空间联系的基础，反映了纹理中灰度分布的性质。基于共生矩阵，可定义许多基于统计法的纹理特性。案例 P55常用指标有行平均、列平均

27、、行标准差、列标准差。第四章图像显示和拉伸第四章图像显示和拉伸1、颜色颜色非彩色消色：黑、白以及从黑过渡到白的一系列的灰色，它们对光谱上各个波长的反射是没有选择性的，称之为中性色。黑白系列的非彩色只能反映物质的光反射率的变化，其在视觉上的感觉是亮度的变化。彩色：是指除黑白系列以外的各种颜色。有色调、明度、色度色调、明度、色度三个基本特性。色调(hue)：是色彩最重要、最基本的特征。根据色调来称呼某种颜色为红色、黄色、绿色。色调是经物体外表反射或投射后到达神经的色光确定的，对于单色光的色调，可以用其波长来确定，而由混合光组成的色彩的色调，可以用组成混合光的各种波长光量的比例来确定。明度：是颜色的

28、亮度在人们视觉上的反映，是从人的感觉上来说明颜色的性质明度与亮度是有区别的，亮度在光度学的概念中可以用光度计测量，与物体反射的光量有关，与人的视觉无关色度：是水中溶解性的物质或胶状物质呈现的米黄色乃至黄褐色的程度，包含有色调信息，但没有亮度信息。2、颜色空间模型颜色空间模型RGB 模型：是根据人眼锥体接受光线的方法构造成的模型，可生成256*256*256=16M 种颜色中心波长 R-700nm、G-546.1nm、B-435.8nmCMYK 模型：是彩色胶片的染料和印刷油墨所形成的颜色空间HIS 模型、LAB 模型3、颜色模型颜色模型一种颜色模型是用一个三维坐标系统及这个系统中的一个子空间来

29、表示，系统中的每个颜色都由一个单点表示。RGB 模型彩色监视器和彩色摄像机：基于笛卡尔坐标系统，黑色在原点，白色在离原点最远的角上，蓝绿色、紫红色和黄色在另 3 个顶角CMY 模型彩色打印机：公式：P62YIQ 模型彩色电视广播：公式：P62 去掉了亮度Y和颜色信息I 和 Q间的紧密联系HIS 模型图像显示处理：强度成分在图像中与颜色信息无关；色调和饱和度成分与人们获得颜色的方式密切相关4、图像的显示图像的显示1显示过程显示过程CPU 从存储介质中读取数据，并以位的形式保存在图像处理器的缓冲存储器图像存储器内。读书器读取存储器里的位数据，然后将此值赋给彩色查找表CLUT 中的 RGB。模数转换

30、器 D/A 将 RGB 数据值变为适当的模拟信号 连续的电信号，该模拟信号用以调整 RGB 电子枪的强度，控制着每个像素在视频 CRT阴极摄像管屏幕显示的亮度。图 P632全色显示全色显示8 位图形处理器有一个相连的 8 位256 个元素的查找表CLUT，该 CLUT 中 RGB 排列均呈 0-255 渐变顺序，每个元素都有对应的 RGB 值。3彩色显示彩色显示电子显示法软拷贝：用彩色监视器显示彩色硬拷贝设备显示：相减混色原理显示彩色图像，使用CMYK 颜色空间。相加混色基本规律：红+绿=黄 红+蓝=紫蓝+绿=青 红+蓝+绿=白三基色：红蓝绿三补色：青紫黄真彩色：选择波段与 RGB 对应假彩色

31、：输入波长与 RGB 不同图像显示使用颜色索引P65 图仔细研究5、图像的彩色合成彩色增强图像的彩色合成彩色增强伪彩色合成伪彩色合成：把单波段灰度图像中的不同灰度级按特定的函数关系变换成彩色，然后进行彩色图像显示，主要通过密度分割来实现。密度分割法：对单波段遥感图像按灰度分级，对每级赋予不同的色彩，使之成为一幅彩色图像。经过密度分割后，图像的可分辨力得到明显提高，如果分级与地位光谱特性的差异对应较好，可以较准确地区分出地物类别。真彩色合成真彩色合成：合成后图像的颜色更接近于自然色，与人对地物的视觉感觉相适应，更容易对地物进行识别。假彩色合成假彩色合成：假彩色合成使用的数据是多波段图像，假彩色合

32、成选用的波段应该与地物的光谱特征为出发点，不同的波段合成方法，用来突出不同的地物信息。标准假彩色方案：TM(432)、MSS(421)、SPOT(321)，植被红色、水体与黑色或蓝色、城镇为深色模拟真彩色合成模拟真彩色合成：由于蓝光容易受大气中气溶胶的影响，有些传感器舍弃了蓝波段，因此通过彩色合成无法得到真彩图像，可通过某种形式的运算得到模拟的红绿蓝3 个通道，然后通过彩色合成近似地产生真彩色图像。SPOT IMAGE公司：红色用 XS2 表示，绿色用XS1+XS2+XS3/3 的波段运算来实现，蓝色采用 XS1波段代替，绿波段单做蓝波段，红波段不变，绿波段用三个波段的平均值代替。ERDAS

33、IMAGING软件：红色用 XS2 表示，绿色用(XS1*3+XS3)/4 波段，蓝色用 XS1 波段。不确定参数法：引入全色波段(P)，红色(aP+(1-a)*XS3)，绿色 2*P*XS2/(XS1+XS2)，蓝色2*P*XS1/(XS1+XS2)，a 引入为了防止出现过饱和0.1-0.5，假设 XS1 和 XS2为 0，则会出现不合理结果。6、图像拉伸图像拉伸拉伸用来改善图像显示的比照度，以波段为处理对象，通过处理波段中单个像素值来实现增强的效果。1灰度拉伸灰度拉伸线性拉伸：全域线性拉伸公式P772%拉伸a 取灰度级 2%，b 取灰度级 98%分段线性拉伸公式 P78灰度窗口切片：为了将

34、某一区间的灰度级和其他部分背景分开清除背景：不在灰度窗口范围内的像素赋值为最小灰度级，在灰度窗口范围内的像素都赋值为最大灰度级保留背景：不在灰度窗口范围内的像素保留灰度值，在灰度窗口范围内的像素都赋值为最大灰度级非线性拉伸：指数变换：对于图像中亮的部分，指数变换扩大了灰度间隔，突出了细节，对于暗的部分，缩小了灰度间隔，弱化了细节对数变换：拉伸图像中暗的部分，压缩亮的部分多波段拉伸：对各个波段分别进行线性或非线性拉伸，再综合2图像均衡化图像均衡化基本思想：对原始图像的像素灰度做某种映射变换，使变换后图像灰度的概率密度呈均匀分布，即变换后图像的灰度级均匀分布，这意味着图像灰度的动态范围得到了增加，

35、从而提高了图像的比照度。基本步骤：统计原图像每一灰度级的像素数和累积像素数计算每一灰度级 xa 均衡化后对应的新值，并对其四舍五入取整，得到新灰度级xb以新值替代原灰度级，形成均衡化后的新图像根据原图像像素统计值对应找到新图像像素统计值，作出新直方图算例，笔记本上有特点：各灰度级中像素出现的频率近似相等原图像上像素出现频率校的灰度级被合并，实现压缩；像素出现频率高的灰度级被拉伸，突出了细节信息3直方图规定化直方图匹配直方图规定化直方图匹配原理：对两个直方图都作均衡化，变成归一化的均匀直方图，以此直方图做中介，再对参考图像作均衡化的逆运算。具体步骤：P84这种方法经常作为图像镶嵌或应用遥感图像进

36、行动态变化研究的预处理工作，可以部分消除由于太阳高度角或大气影响造成的相邻图像的色调差异，从而可以降低目视解译的错误。第五章图像校正第五章图像校正1 1、辐射校正、辐射校正定义：消除图像中依附在幅亮度中的各种失真的过程成为幅亮度校正，简称辐射校正目的：尽可能消除因传感器自身条件、薄雾等大气条件、太阳位置和角度条件及某些不可防止的噪声引起的传感器的测量值与目标的光谱反射率或光谱幅亮度等物理量之间的差异；尽可能恢复图像的本来面目，为遥感图像的分类、分割、解译等后续工作做好准备内容：传感器端的辐射校正、大气校正、地表辐射校正2 2、辐射传输、辐射传输基本概念P87-88重要：辐射通量单位时间内通过某

37、一外表的辐射能量成为辐射通量，单位W幅照度单位时间内单位面积上接受的辐射能量，单位W/M2幅亮度/幅照度沿辐射方向、单位面积、单位立体角上的辐射通量，单位W/(M2*SR)在介质内部，反射率、吸收率和透射率的和为1反照率界面反射的辐照度与内部反射的辐照度之和/入射的辐照度3 3、电磁波的大气传输、电磁波的大气传输1可见光和红外传输由于空气分子和悬浮颗粒的散射，可见光在大气层传输时会被削弱，传感器接收到大气散射部分的电磁波称为程辐射或路径辐射。分子散射对波长较短的电磁波UV 紫外线和蓝光有影响，对红外线没有影响。瑞利散射悬浮颗粒散射量的大小和角度的变化取决于悬浮颗粒的总容量，粒径分布、介电常数和

38、微粒的形状。悬浮颗粒散射的影响很难校正，因此过度浑浊不能进行大气校正条件是标准大气和洁净空气。2热红外传输CO2、H2O、N2O、O3 是主要的吸收气体获取地表温度和辐射率的另一种方法是别离窗口技术劈窗方法，它利用具有不同透射率的大气窗口进行计算，NOAA 先进的高精度辐射计(AVHRR)图像可以作为地物反演的数据4 4、辐射传输理论、辐射传输理论传感器接受目标物辐射或反射的电磁波，由此形成的遥感原始图像与目标相比存在失真，这是因为在太阳大气目标大气传感器的辐射传输过程中存在着许多干扰因素，使接受的信号不能准确反映地表物理特征光谱反射率、光谱幅亮度等。1因素：A 大气分子及气溶胶的瑞利散射和米

39、氏散射、分子及气溶胶的吸收、散射以及散射吸收的耦合作用。B 外表因素的奉献传感器高度较低那么外表因素的奉献不可忽略C 地形因素的奉献目标的高度和坡向D 太阳辐射光谱的影响太阳本身为一个黑体，其光谱辐射遵循普朗克定律2基本流程 P923大气辐射传输方程公式P93可见光-短红外光谱区0.4-2.5um地球本身的辐射可以忽略，只考虑太阳光的辐射传输传感器所能接受的太阳光：太阳光直射到地表后地表的反射辐射、被大气散射辐射的太阳光在地表的反射辐射、大气上行散射辐射程辐射/路径辐射中红外和热红外区域传感器接收的能量主要来源于地球的热辐射，其能量为：地表热发射辐射、大气下行热发射辐射被地表反射后的辐射、大气

40、上行热发射辐射。4应用大气辐射传输模型进行遥感图像的大气校正需要解决两个问题A 有关大气介质特征参数的获取B 具体使用的大气辐射传输模型5 5、辐射误差、辐射误差产生原因：传感器的响应特性、外界自然环境包括大气云和雾、太阳辐射1传感器的响应特性光学摄影机引起的辐射误差、光电扫描仪引起的辐射误差光电转换误差、探测器增益2大气电磁波在大气传输中，受到大气中分子和微小粒子的作用。这些分子和微小粒子对光波多次作用的结果即散射，它岁电磁波波长和散射体大小的不同而不同。散射：选择性散射：波长越短散射越强瑞利散射由远小于光波长的气体分子所引起，大小与波长的四次方成反比米氏散射由大小与波长相当的颗粒气溶胶：烟

41、、水蒸气、霾引起，也成为气溶胶散射，大小与波长成反比非选择性散射：尘埃、雾、云以及大小超过光波长10 倍的颗粒引起，对各种波长予以同等散射。散射增加了到达卫星传感器的能量，从而降低了遥感图像的反差 最大亮度值/最小亮度值，降低了图像的分辨率，因此必须进行校正。低分辨率的空间范围较大，不认为图像中各处的大气散射是均匀的，往往要分区校正。3太阳辐射太阳位置太阳高度角、方位角、地形起伏4其他误差各检测器的差异、干扰、故障灯等6 6、系统辐射误差校正、系统辐射误差校正1光学镜头的非均匀性引起的边缘减光现象的校正边缘减光：由于透镜光学特性的非均匀性，在成像平面上边缘部分比中间部分暗2条纹条纹主要是由检测

42、器引起的，条纹误差判定和消除的常用方法：平均值法、直方图法及在垂直扫描线方向上采用最近邻点法或三次褶积法。3斑点斑点误差主要是由噪声或磁带的误码率等造成，在图像中往往是分散和孤立的。校正后的斑点亮度取其领域像素亮度值的平均值或用三次褶积法进行修正。边缘附近的斑点不进行消除，图像四周的像素不进行斑点消除。4灰度一致化应该先完成几何精纠正，具有相同的地理坐标等概率变换优点：简单，变换效果较好利用了重叠部分的全部像素灰度，属于非线性变换缺点：按概率分布进行像素灰度变换，因而存在位置配准误差线性灰度变换该方法简单易行，N 足够大时有一定精度7 7、传感器端的辐射校正、传感器端的辐射校正1可见光和近红外

43、波段的辐射定标2红外波段的辐射校正热红外波段，星上传感器入瞳出接收的总辐射：通过大气向上幻术的直接地面辐射、由大气自身向上传输的辐射、大气向下辐射到达地面再经地面反射后通过大气向上传输的辐射。3图像的灰度级和辐亮度图像上的像素值为灰度级，实际的电磁波辐射强度为辐亮度，在图像数字化的时候，电磁波辐亮度被量化为灰度级，进行反演的时候需要将灰度级转换为辐亮度。8 8、大气校正、大气校正消除由大气散射引起的辐射误差的处理过程称为大气校正。1 1统计学方法统计学方法通常将野外实地光谱测试获得的无大气影响的辐射值与卫星传感器同步观测结果进行回归分析，确定校正量主要有：内部平均法：校正后为相对反射率值，考虑

44、乘性奉献公式P103优点：可以大大消除地形阴影和其他整体亮度的差异缺点：该方法假设地面变化时充分异构的，光谱反射特性的空间变化会相互抵消。假设不成立时，得到的光谱则具有虚假性平场域法：校正后为相对反射率值，考虑乘性奉献公式P103要求图像具有一个光谱反射率曲线变化相对平坦的、比较均一的区域，该区域的平均光谱受大气辐射、大气散射和吸收影响的共同控制。两个重要假设：区域的平均光谱没有明显的吸收特征；区域辐射光谱主要反映当时大气条件下的太阳光谱。缺点：不适合大量多条带高光谱数据的处理人工查找方法有一定的随意性经验线性法：校正后为绝对反射率值，考虑加性奉献公式P103需要两个以上光谱均一、有一定面积大

45、小的目标分别为按目标和量目标，作为定标点定标点要求：选择尽可能各向同性的均一地物，且面积足够大；地物在光谱上要跨越尽可能宽的反射光谱段，明暗目标之间有足够的差异尽可能与研究区域保持同一海拔高度仅适用于包含地面实况数据的图像。实测光谱回归方法：2 2辐射传递方程计算法辐射传递方程计算法辐射传输模型：一个模拟的太阳辐射光谱开始，计算太阳高度的辐射影响根据成像时间以及大气散射和吸收。6S 模型利用辐射传递方程通常只能得到近似解，改良的方法是在获取图像的同时，利用搭载在同一平台上的其他传感器获取气溶胶密度和水蒸气浓度数据，然后利用这些数据进行大气校正。3 3波段比照法波段比照法理论依据：大气散射的选择

46、性，即大气散射对短波影响大，对长波影响小回归分析法：在不受大气影响的波段和待校正的某一波段图像中，选择最黑区域通常为高山阴影区中的一系列目标，将每一目标的两个待比照的波段亮度值提取出来进行回归分析暗像素法。如果不是在各个波段全黑，回归分析得到的拟合曲线就存在散射的影响。直方图法：如果图像内包括暗色地物或地形阴影，可从各个波段中减去其最小的亮度值或一个阴暗地区的平均亮度值进行校正。值得注意的是假设是找不到全黑的区域，则两种方法都不能使用。9 9、地面辐射校正、地面辐射校正1太阳辐射校正：主要校正太阳高度角导致的辐射误差公式法P107和波段比值法2地形辐射校正：校正坡度和坡向带来的影响1010、几

47、何校正、几何校正1几何误差的来源静态误差：成像过程中，传感器相对地球外表呈静止状态时所具有的各种误差内部误差：由于传感器自身的性能、技术指标偏离标准数值造成外部误差：传感器本身处于正常工作的条件下，由传感器意外的各因素造成动态误差：由于成像过程中地球的旋转所造成的图像误差2 2几何精纠正几何配准几何精纠正几何配准是把不同传感器具有几何精度的图像、地图或数据集中的相同地物元素精确地彼此匹配、叠加在一起的过程。相对纠正：基础数据集是图像绝对纠正：基础数据集是标准的地图几何纠正必须在遥感图像信息提取之前进行纠正内容：系统几何纠正、投影变形纠正、几何精纠正重要性：只有在进行纠正后，才能对图像信息进行各

48、种分析，制作满足测量和定位要求的各类遥感专题图在同一区域，应用不同传感器、不同光谱范围及不同成像时间的各类图像数据进行电脑自动分类、地物特征的变化检测或其他应用处理时，必须进行图像间的空间配准，保证不同图像间的几何一致性。利用遥感图像进行地形图测图或更新时，要求图像具有较高的地理坐标精度3基本原理回避成像的空间几何过程，直接利用地面控制点数据对遥感图像的几何畸变本身进行数学模拟，并且认为遥感图像的总体畸变可以看做是挤压、扭曲、缩放、偏移以及更高层次的基本变形的综合作用的结果。基本技术：同名坐标变换方法，即通过在基础数据和图像中分别寻找地面控制点的同名坐标，并借此建立变换关系来进行几何精纠正4操

49、作步骤准备工作输入原始数据图像确定工作范围选择地面控制点选择地图投影匹配地面控制点和像素位置评估纠正精度坐标变换重采样输出纠正后图像5 5地面控制点地面控制点GCPGCP控制点数目和分布：最小值按位置系数多少来确定，2030 个；分布均匀图像中控制点确实定：容易分辨、相对稳定、特征明显的位置如道路交叉点、河流弯曲或分叉处、海岸线弯曲处、湖泊边缘、飞机场跑道地面控制点坐标确实定：地形图或现场实测获得控制点的数目不如控制点的分布对纠正结果的影响大6 6纠正方法纠正方法多项式纠正方程和共线方程多项式纠正方程原理直观、计算简单，对地面相对平坦的图像具有足够高的纠正精度，较为常用注意问题：多项式纠正的精

50、度与地面控制点的精度、分布和数量及纠正的范围有关精度越高、分布越均匀、数量越多越好采用多项式纠正时，在GCP 处拟合较好，但在其他点的误差可能会较大，平均误差小不能保证各点误差都小多项式阶数确实定，取决于图像中几何变形程度的认识7重采样像素位置的变换、像素值的变换像素位置的变换：按选定的纠正方程把原始图像中的各个像素变换到输出图像相应的位置上去直接成图法正解法间接成图法反解法：根据纠正公式计算规则网的地理坐标 XY在原始图像中对应的位置xy，根据 xy与周围像素之间的关系内插产生新的像素值，然后把像素值写到XY中。内插计算像素值的过程成为数字图像的重采样重采样方法：最近邻重采样：优点：算法简单