四章遥感图像处理10GIS210618课件

第四章遥感图像处理4.1光与色4.2图像校正辐射校正几何校正4.3图像增强4.4多源信息复合多源遥感数据融合遥感数据与非遥感数据的融合第四章遥感图像处理4.1光与色14.1光与色光:在0.38-0.76μm范围内,能被眼睛感觉的、产生视觉现象的辐射。0.7μm为红色0.58μm为黄色0.51μm为绿色0.47μm为蓝色4.1光与色光:在0.38-0.76μm范围内,能2颜色视觉

在光亮条件下，人眼能分辨各种颜色；当光谱亮度降低到一定程度，人眼的感觉便是无彩色的，光谱变成不同明暗的灰带。亮度对比：视场中对象相对于背景的明亮程度。适当改变对比度，可以提高视觉效果。

颜色对比：视场中相邻区域不同颜色的相互影响。颜色的对比会产生不同的视觉效果。L背景C=|L对象-L背景|L对象C=△对象颜色视觉L背景C=|L对象-L背景|L对象C=3颜色对比：两种颜色相互影响的结果，使每种颜色会向其影响色的补色变化。在两种颜色的边界，对比现象更为明显。例如：在一块品红的背景上放一小块白纸或灰纸，用眼睛注视白纸中心几分钟，白纸会表现出绿色。如果背景是黄色，白纸会出现蓝色。颜色对比：两种颜色相互影响的结果，使每种颜色会向其影响色的补4颜色的性质色调：颜色的种类明度：颜色的明亮程度饱和度：表示颜色纯洁程度的物理量所有颜色是对某段波长有选择地反射而对其他波长吸收的结果。

如：如果物体反射可见光中0.6m以上的波长，则看起来是红色，反射从0.55m其反射率偏低便成成了棕红色。颜色的性质由色调、明度和饱和度来描述。颜色的性质色调：颜色的种类所有颜色是对某段波5四章遥感图像处理10GIS210618课件6四章遥感图像处理10GIS210618课件7四章遥感图像处理10GIS210618课件8颜色立体橙色调蓝绿白饱和度黄黑紫蓝绿红明度为了形象地描述颜色三个特性之间的关系，通常用颜色立体来表现一种理想化的示意图。颜色立体橙色调蓝绿白饱和度黄黑紫蓝绿红明度为9孟赛尔颜色立体：中央轴代表无彩色的明度等级，顶部白为10，底部黑为0，从0至10共分为10个明度级；水平剖面是色调，沿顺时针方向分别为红、红黄、黄、黄绿、绿、绿蓝、蓝、蓝紫、姿、紫红10中色调；颜色离开中央轴的水平距离代表饱和度的变化。任何颜色在孟塞尔系统中都可以用三个坐标值：色调、明度和彩度（饱和度）表示。孟赛尔颜色立体：10加色法加法三原色：红、绿、蓝加色法：将三原色中的两种或三色等量相加而产生另一种颜色的方法。加色法的应用：

彩色显示器的显象管++++加色法加法三原色：++++11减色法减法三原色：

黄、品红、青

即蓝、绿、红的补色减色法：由白色光减去加法三原色中的一种或两种色光而产生另外一种颜色的方法。减色法的应用：遥感影像的色彩摄影、彩色印刷、电影胶片的印染。---减色法减法三原色：---12彩色合成加色法彩色合成减色法彩色合成光学增强处理光学增强处理彩色合成光学增强处理13加色法彩色合成合成仪法：是将不同波段的黑白透明片分别放入有红、绿、蓝滤光片的光学投影通道中精确配准和重叠，生成彩色影像的过程。

分层曝光法：指利用彩色胶片具有的三层乳剂，使每一层乳剂依次曝光的方法。加色法彩色合成合成仪法：是将不同波段的黑白透明片分别放入有红14减色法彩色合成染印法：是一种使用特别浮雕片、接受纸和冲显染印药制作彩色合成影像的方法。印刷法：利用普通胶印设备，直接使用不同波段的遥感底片和黄、品红、青三种油墨，经分色、加网、制版，套印成彩色合成图像。重氮法：利用重氮盐的化学反应处理彩色单波段影像透明片的方法。利用减色法原理，使白光经过多种（层）乳剂或染料或滤色片等而反射或投射出来的合成彩色。减色法彩色合成染印法：是一种使用特别浮雕片、接受纸和冲显染印15光学增强处理相关掩膜处理方法：把几何位置完全配准的原片，制成不同密度、不同反差的正片或负片，通过它们的各种不同叠加方案改变原有影像的显示效果，达到信息增强目的的方法。改变对比度：使用两张同波段同地区的负片（或正片）进行合成，一张反差适中，另一张反差较小，合成后反差一般加大，从而提高了对比度；显示动态变化：不同时相同一地区的正负片影像叠合掩膜，当被叠合影像反差相同时，凡密度发生变化的部分就是动态变化的位置。边缘突出：先将两张反差相同的正片和负片叠合，叠合配准后，再沿希望突出的线形特征的垂直方向错位。目的在于突出线形特征。光学增强处理相关掩膜处理方法：把几何位置完全配准的原片，制成164.2数字图像的校正遥感数字图像处理：利用计算机对遥感图像及其资料进行的各种技术处理。数字图像处理的优点：快捷、准确、客观地提取遥感信息适应地理信息系统的发展4.2数字图像的校正遥感数字图像处理：17一、数字图像及其直方图数字图像：遥感数据有光学图像和数字图像之分。数字图像是能被计算机存储、处理和使用的用数字表示的图像。数字化：将连续的图像变化，作等间距的抽样和量化。通常是以像元的亮度值表示。数字量和模拟量的本质区别：连续变量，离散变量。数字图像的表示：矩阵函数一、数字图像及其直方图数字图像：遥感数据有光学图像和数字图像18四章遥感图像处理10GIS210618课件19令：x∈[0，xmax]，y∈[0，ymax]从而有：f(x,y)=f(0,0)f(1,0)f(1,N-1)f(M-1,0)f(0,1)f(0,N-1)f(1,1)f(M-1,1)……….….f(M-1,N-1)距阵中的每一元素代表图像中的一个像元。f(x,y)为像元的亮度值。令：x∈[0，xmax]，y∈[0，ymax]从而有：f(20灰度直方图：以每个像元为单位，表示图像中各亮度值或亮度值区间像元出现的频率的分布图。直方图的作用：直观地了解图像的亮度值分布范围、峰值的位置、均值以及亮度值分布的离散程度。直方图的曲线可以反映图像的质量差异。正态分布：反差适中，亮度分布均匀，层次丰富，图像质量高。偏态分布：图像偏亮或偏暗，层次少，质量较差。注意：图像直方图是描述图像质量的可视化图表。在图像处理中，可以通过调整图像直方图的形态，改善图像显示的质量，以达到图像增强的目的。灰度直方图：以每个像元为单位，表示图像中各亮度值或亮度值区间21直方图直方图22二、辐射校正（Radiometriccorrection）辐射畸变：

地物目标光谱反射率的差异在实际测量时，受到传感器本身、大气辐射等其他因素的影响而发生改变。这种改变称为辐射畸变。影响辐射畸变的因素传感器本身的影响：导致图像不均匀，产生条纹和噪音。大气对辐射的影响图像的灰度←辐射强度←地物波谱反射率←地物不同。二、辐射校正（Radiometriccorrection23大气影响的定量分析：大气的主要影响是减少了图像的对比度，使原始信号和背景信号都增加了因子，图像质量下降。大气影响的粗略校正：通过简单的方法去掉程辐射度（散射光直接进入传感器的那部分），从而改善图像质量。直方图最小值去除法回归分析法:校正的方法是将波段b中每个像元的亮度值减去a，来改善图像，去掉程辐射。

大气影响的定量分析：大气的主要影响是减少了图像的对比度，使24直方图最小值去除法基本思路:每幅图像上都有辐射亮度或反射亮度应为0的地区，而事实上并不等于0,说明亮度最小值必定是这一地区大气影响的程辐射度增值。校正方法:将每一波段中每个像元的亮度值都减去本波段的最小值。使图像亮度动态范围得到改善，对比度增强，从而提高了图像质量。直方图最小值去除法基本思路:每幅图像上都有辐射亮度或反射亮度25四章遥感图像处理10GIS210618课件26四章遥感图像处理10GIS210618课件27三、几何校正几何畸变：遥感图像的几何位置上发生变化，产生诸如行列不均匀，像元大小与地面大小对应不准确，地物形状不规则变化等变形。几何畸变是平移、缩放、旋转、偏扭、弯曲等作用的结果。三、几何校正几何畸变：遥感图像的几何位置上发生变化，产生诸如281、遥感影像变形的原因遥感平台位置和运动状态变化的影响：

航高、航速、俯仰、翻滚、偏航。地形起伏的影响：产生像点位移。地球表面曲率的影响：一是像点位置的移动；二是像元对应于地面宽度不等，距星下点愈远畸变愈大，对应地面长度越长。大气折射的影响：产生像点位移。地球自转的影响：产生影像偏离。1、遥感影像变形的原因遥感平台位置和运动状态变化的影响：航29遥感平台运动状态变化引起的误差遥感平台运动状态变化引起的误差30地形起伏引起的误差地形起伏引起的误差31地球表面曲率引起的误差1、像点位置的移动2、像元对应地面宽度不同地球表面曲率引起的误差1、像点位置的移动32大气折射引起的误差像点位置的移动大气折射引起的误差像点位置的移动33地球自传引起的误差影像偏离地球自传引起的误差影像偏离342、几何畸变校正基本思路：把存在几何畸变的图像，纠正成符合某种地图投影的图像，且要找到新图像中每一像元的亮度值。具体步骤

1）计算校正后每一点所对应原图中的位置；2）计算每一点的亮度值。2、几何畸变校正基本思路：把存在几何畸变的图像，纠正成符合某35几何校正流程图开始显示图像文件启动几何校正模块选择计算转换模型采集地面控制点图像重采样检验地面控制点结束几何校正流程图开显示图像文件启动几何校正模块选择计算转换模型36控制点的选取数目的确定:齐次多项式的项数N与其阶数n的关系N=（n+1）（n+2）/2地面控制点(GCP)个数L与N之间的关系L≥N选择的原则易分辨、易定位的特征点：道路的交叉口，河流汇流处，湖心岛，飞机场，桥，水库大坝等。特征变化大的地区应多选些。尽可能满幅均匀选取。控制点的选取37四章遥感图像处理10GIS210618课件38四章遥感图像处理10GIS210618课件39四章遥感图像处理10GIS210618课件40几种采样方法的优缺点：几种采样方法的优缺点：414.3图像增强辐射增强空间滤波彩色增强图像运算多光谱变换增强目视效果、提高图像质量和突出所需信息，有利于分析判读或作进一步的处理。4.3图像增强辐射增强增强目视效果、提高图42辐射增强（对比度变换）：通过改善图像像元的辐射亮度值以提高图像质量的方法一、辐射增强（对比度变换）直方图是图像的灰度——像素数统计图，即对于每个灰度值，统计在图像中具有该灰度值的像素个数，并绘制成图形。直方图是图像的灰度——像素数统计图，即对于每个灰度值，统计在图像中具有该灰度值的像素个数，并绘制成图形。辐射增强（对比度变换）：通过改善图像像元的辐射亮度值以提高图43一幅适合用统计方法分析的图像其直方图应是正态分布的，因此，对比度变换就是将亮度值不符合正态分布的图像尽量使其亮度呈正态分布。一幅适合用统计方法分析的图像其直方图应是正态44四章遥感图像处理10GIS210618课件45线性变换变换后亮度值xb1015103020OO'变换前亮度值xaxb=b2-b1a2-a1(xa-a1)+b1例：一图像亮度值为0-15，现在把它拉伸到0-30。线性方程：Xb=2Xa线性变换变换后亮度值xb1015103020OO'变换前亮度46变换后图像变换前后直方图对比变换后图像变换前后直方图对比47四章遥感图像处理10GIS210618课件48非线性变换xb=（axa+1）+c对数变换xb=beaxa+c指数变换

在亮度值较高的部分拉伸，即扩大亮度间隔。在亮度值较低的部分压缩，即缩小亮度间隔。在亮度值较高的部分压缩，即缩小亮度间隔。在亮度值较低的部分拉伸，即扩大亮度间隔。非线性变换xb=（axa+1）+c对数变换xb=bea49二、空域邻域增强空间邻域增强，即空域邻域卷积运算图像平滑图像锐化二、空域邻域增强空间邻域增强，即空域邻域卷积运算50图像卷积运算：在空间域上对图像作局部检测选定一卷积函数——模板，M×N图像移动模板，对每个像元进行卷积运算

R=w1z1+w2z2+…+wnzn其中：wii=1,2,…,n是模板的系数

zii=1,2,…,n是被计算像素及其邻域像素的值123456789图像卷积运算：在空间域上对图像作局部检测12345678951以下图为例，试看卷积实验过程。1011111210101111121111111210102001112111111121211111212121010131112浅灰字卷积模板元素值暗红数字图像灰度值2142121

卷积模板（低通滤波器）101011121010111112101235231210345834122335231211101313121010131112

滤波结果图以下图为例，试看卷积实验过程。101152中值滤波（有时也用均值滤波）主要用途：平滑图像、去除噪音计算公式：R=mid{zk|k=1,2,…,9}最大值滤波，去除胡椒噪声（暗）主要用途：寻找最亮点计算公式：R=max{zk|k=1,2,…,9}最小值滤波：去除椒盐噪声（亮）主要用途：寻找最暗点计算公式：R=min{zk|k=1,2,…,9}123456789一般用中值滤波，最大值，最小值很少用平滑：图像中出现某些亮度变化过大的区域，或出现不该有的亮点时，可以减少变化，使亮度平缓。中值滤波（有时也用均值滤波）123456789一般用中值滤波53中值滤波算法的实现将模板区域内的像素排序，求出中值。例如：3x3的模板，第5个是中值，

5x5的模板，第13个是中值，

7x7的模板，第25个是中值，

9x9的模板，第41个是中值。对于同值像素，连续排列。例如：（10，15，20，20，20，20，20，25，100）最终值为：20中值滤波算法的实现54锐化突出图像中的边缘、线状目标、亮度突变区域。锐化后的图像失去原图像的特征。具体方法很多：罗伯特梯度，索伯尔梯度，拉普拉斯算法，定向检测等。相关算法详见：课本p118锐化突出图像中的边缘、线状目标、亮度突变区域。55四章遥感图像处理10GIS210618课件56把数字图像组合转换成彩色图形，或者把各种增强或分类图像组合叠加，以彩色图像显示出来。（彩色的视觉分辨能力比黑白高）方法：假彩色密度分割；彩色合成。三、彩色变换三、彩色变换57概念：单波段黑白遥感图像可按亮度分层，对每层赋予不同的色彩，使之成为一幅彩色图像。这种方法又叫密度分割。分层方案的确定：分层方案与地物光谱差异对应合适，可以较好地区分地物类别。处理过程效果分析1、单波段彩色变换（假彩色密度分割）概念：单波段黑白遥感图像可按亮度分层，对每层赋予不同的色彩，58处理过程输入图像显示直方图确定分割的等级数，并计算分割的间距像元亮度值转换为像元新值赋色处理过程输入图像显示直方图确定分割的等级数，并计算分割的间距59四章遥感图像处理10GIS210618课件60以不同的色彩表示图像的色调变化，增强了图像的显示能力同一地物或现象可能被分割成两种不同密度并以不同的颜色显示出来，或同一色彩却表示两种以上不同的地物，造成判读错误。效果分析以不同的色彩表示图像的色调变化，增强了图像的显示能力效果分析61单波段彩色变换（伪彩色增强）红色转换器绿色转换器蓝色转换器

IR(x,y)

IG(x,y)

IB(x,y)输入三个电子枪，可得到其颜色内容由3个变换函数调制的混合图像单波段彩色变换（伪彩色增强）红色转换器绿色转换器蓝色转换器62a)红色变换特性b)绿色变换特性IRL0

L/2 L

f(x,y)

L/2 L

f(x,y)IGL0L/2 L

f(x,y)IBL0L/2 Lf(x,y)ILc)蓝色变换特性d)合成特性a)红色变换特性b)绿色变换特性I63四章遥感图像处理10GIS210618课件642真彩色增强从彩色到彩色的一种转换目的1）引起人们的特别关注2）人眼对不同颜色的灵敏度不同2真彩色增强从彩色到彩色的一种转换65多波段色彩变换概念：利用计算机将同一地区不同波段的图像存放在不同通道的存储器中，并依照彩色合成原理，分别对各通道的图像进行单基色变换，在彩色屏幕上进行叠置，从而构成彩色合成图像。合成方案：彩色合成图像分为真彩色图像和假彩色图像。多波段色彩变换概念：利用计算机将同一地区不同波段的图像存放在66TM1TM2TM3TM4TM5TM6TM7LandsatTM5sub-sceneshowingtheregionaroundtheAlpinforschungszentrumRudolfshütteTM1TM2TM3TM4TM5TM6TM7LandsatT67TM7,4,1TM5,7,2TM5,4,3TM4,3,2TMTMTMTM68概念：两幅或多幅单波段影像，完成空间配准后，通过一系列运算，可以实现图像增强，达到提取某些信息或去掉某些不必要信息的目的。

原理：地物不同波段的光谱差异。

差值运算比值运算加、减、乘、除四则运算等四、图像运算概念：两幅或多幅单波段影像，完成空间配准后，通过一系列运算，69差值运算f(i,j)=f1(i,j)-f2(i,j)四、图像运算差值运算f(i,j)=f1(i,j)-f2(i,j)四、图70差值运算目的是突出两图像亮度差大的部分例如：用红外波段减红波段，植被的反射率在这两个波段的差异很大，而土壤和水在这两个波段的反射率差值很小，因此相减后的图像可以把植被信息突出出来。

植被光谱曲线差值运算目的是突出两图像亮度差大的部分植被光谱曲线71四章遥感图像处理10GIS210618课件72四章遥感图像处理10GIS210618课件73``74差值运算另外应用于同一地区不同时相的动态监测：

森林火灾后的过火面积洪水淹没区域、受灾面积及损失城市的扩展差值运算另外应用于同一地区不同时相的动态监测：75四章遥感图像处理10GIS210618课件76四章遥感图像处理10GIS210618课件77在进行卫星图象处理时,遇到一山体阴影,其数据组合如下:请选择恰当的处理方法去处阴影.采用比值方法Rmn（j，K）=Fm（j，K）/Fn（j，K）则：光照区：Rmn（j，K）=28/42=0.66

阴影区:Rmn（j，K）=22/34=0.65两区的亮度值近似,从而消除阴影.

亮度值地类4波段5波段光照区2842阴影区2234比值运算去除地形影响植被指数两幅同样行、列数的图像，对应像元的亮度值相除（除数不为0）就是比值运算。该运算常用于突出遥感影像中的植被特征、提取植被类型或估算植被生物量，这种算法的结果称为植被指数。在进行卫星图象处理时,遇到一山体阴影,其数据组合如下:请选择78Minnaert校正图像地形去除Minnaert校正图像地形去除79四章遥感图像处理10GIS210618课件80T=45.399-0.004*E-0.816*L图像运算模型：T=45.399-0.004*E-0.816*LT=45.399-0.004*E-0.816*L图像运算模型81五、多光谱变换多光谱变换：通过选择不同波段遥感数据之间的变换函数，经过分析、计算和变换，达到数据压缩和增强（提取）地物信息的目的。实质：对遥感图像实行线性变换，使光谱空间的坐标按一定规律进行旋转。X=X2XiX1Xn=[x1，x2，xi，xn]T多光谱空间：只表示各波段光谱之间的关系，没有图像空间的意义。…五、多光谱变换多光谱变换：通过选择不同波段遥感数据之间的变换82K—L变换

K—L变换又称主成分变换，是利用变换矩阵对多光谱图像进行线性组合，最终产生一组新的多光谱图像，其目的是数据压缩和图像增强。Y=AXy1y2y3yn=Φ11Φ21Φ12Φ22Φ1nΦ2nΦn1Φn2ΦijΦmn·x1x2xixnX为变换前的多光谱空间像元矢量，Y为变换后的主成分空间像元矢量，A为变换矩阵。使输出图像Y的各个分量之间具有较小的相关性K—L变换K—L变换又称主成分变换，是利83ΔK—L变换的特点(简单应用举例因子反演)①从几何意义来看，变换后的主分量空间坐标系与变换前的空间坐标系相比旋转了一个角度，而且新坐标系一定指向数据信息量较大的方向(主分量方向)。②就变换后的新波段主分量而言，K—L变换后的各个新波段主分量包括的信息量不同，呈逐渐减少趋势。其中，第一主分量集中了最大的信息量，常常占80％以上，第二、第三主分量的信息量依次快速递减，到第n分量信息几乎为0。

③由于K—L变换对不相关的噪声没有影响，所以信息减少时，便突出了噪声，最后的分量几乎全是噪声。所以这种变换又可分离出噪声。ΔK—L变换的特点84K—T变换(缨帽变换)(1)K—T变换是Kauth—Thomas变换的简称，这种变换也是一种线性组合变换，其变换公式为：Y=BX这里X为变换前的多光谱空间的像元矢量，y为变换后的新坐标空间的像元矢量，B为变换矩阵。一种坐标空间发生旋转的线性变换，但旋转后的坐标轴不是指向主成分方向，而是指向了与地面景物有密切关系的方向。主要针对TM的1至5和第7波段，低分辨率的热红外(第6波段)波段不予考虑。

B与矢量相乘后得到新的6个分量，

K—T变换(缨帽变换)(1)K—T变换是Kauth—Thom85经过研究，新分量中的前三个分量与地面景物的关系明确：yl为亮度，实际是TM的6个波段的加权和，反映了总体的反射值。y2为绿度，从变换矩阵B的第二行系数看，对波长较长的红外波段5和7即有很明显的抵削，剩下4与1，2，3，刚好是近红外与可见光部分的差值，反映了绿色生物量的特征。y3为湿度，这个分量反映了可见光至近红外波段(1至4)与波长较长的红外(5，7)波段的差值，而5，7两波段对土壤湿度和植被湿度最为敏感，便于反映出湿度特征。y4,y5,y6这三个分量与地物没有明确的对应关系，因此K—T变换后只取前三个分量。这样也实现了数据压缩。

K—T变换的研究主要针对TM数据和以前使用过的MSS数据，目前应用范围较窄，但它抓住了地面景物，特别是植键和土壤在多光谱空间中的特征。对于扩大陆地卫星TM影像数据分析在农业方面应用有重要意义。经过研究，新分量中的前三个分量与地面景物的关系明确：y86K—T变换Y=BXB为变换距阵、X为变换前的多光谱空间像元矢量，Y为变换后的新坐标空间像元矢量。K—T变换Y=BXB为变换距阵、X为变换前的多光谱空间像元矢87数据融合HLS变换流程图

亮度（Lightness）、色调（Hue）、饱和度(Saturation）数据融合HLS变换流程图亮度（Lightness）、色调（88数据融合的目的：提高图像空间分辨率，增强特征显示能力，

即：在提高融合图像分辨率的同时，尽可能保持光谱信息，以提高图像分析精度。（1）模型的选取目前常用的彩色模型分为两类，一类是面向硬设备的RGB模型，一类是面向彩色处理的HLS模型.（2）变换的步骤：

步骤一：将多光谱图像的三个波段进行彩色合成RGB；

步骤二：将产生的RGB图像进行HLS变换，得到亮度（I）、色调（H）、饱和度（S）三幅图像；

步骤三：高分辨率色图像与HLS变换后的亮度图像进行直方图匹配，生成与亮度分量图像有相似直方图分布特性的图像；

步骤四：用第三步生成的全色图像代替亮度（I）图像，并同H、S图像进行HIS反变换，得到融合图像。

数据融合的目的：提高图像空间分辨率，增强特征显示能力，89RGB显示HLS显示在HLS坐标中进行有目的的处理，然后反变换到RGB坐标中进行显示，使图像彩色增强获得更佳的效果，并可使不同分辨率的图像进行复合获得最佳的综合显示效果。RGB显示HLS显示在HLS坐标中进行有目的的处理，然后反变90多光谱2.44m全色0.61m主成分变换法QuickBird影像融合效果多光谱2.44m全色0.61m主成分变换法QuickBird91几种典型正变换的表达式彩色变换正变换备注单六角锥变换I=DS=(D-d)/DH=(5+(R-B)/(R-G))/6当R=D.G=dH=(1-(R-G)/(R-B))/6当R=D.B=dH=(1+(G-R)/(G-B))/6当G=D.B=dH=(3-(G-B)/(G-R))/6当G=D.R=dH=(3+(B-G)/(B-R))/6当B=D.R=dH=(5-(B-R)/(B-G))/6当B=D.G=dD=max[R,G,B]D=min[R,G,B]三角形变换I=(R+G+B)/3H=(G-B)/(3(I-B))S=1-B/I,当B=dH=(B-R)/(3(I-R))S=1-R/I.当R=dH=(R-G)/(3(I-G))S=1-G/I当G=d

D=min[R,G,B]球体变换I=(M+m)/2S=(M-m)/(2-M-m)H=60(2+b-g).当r=MH=60(4+r-b).当g=MH=60(6+g-r).当b=MD=max[R,G,B]D=min[R,G,B]r=(D-R)/(D-d)g=(D-G)/(D-d)b=(D-B)/(D-d)M=max[r,g,b]m=min[r,g,b]圆柱体变换I=(R+G+B)H=tg-1（）＋CS=几种典型正变换的表达式彩色变换正变换备注单六角锥变换I=D92融合后的遥感影像（盐渍化）中巴资源卫星的4、3、2波段（20m×20m）与ETM的8波段（15m×15m）融合而成的图像融合后的遥感影像（盐渍化）中巴资源卫星的4、3、2波段（2