一种基于地形和遥感影像融合技术提取冲沟方法
【专利摘要】本发明提供一种基于地形和遥感影像融合技术提取冲沟方法，该方法是基于遥感RGB（红、绿、蓝）色彩空间同HIS（色相、色彩亮度、饱和度）色彩空间之间的变换，以线性标准化地表粗糙度作为太阳西北方向地形阴影图SRM（Shaded?Relief?Model）的权重，以作为色彩亮度分量权重，求和构成新的色彩亮度分量，并基于新的色彩亮度图像，进行HIS色彩空间向RGB色彩空间变换，融合地形信息和遥感图像信息，使得遥感图像中冲沟表现为凹陷，山峰表现为凸起。基于融合后的遥感图像，结合基于DEM（数字高程模型）计算的冲沟沟缘线处坡度阈值数据，进行冲沟沟缘线的解译。同传统的基于遥感图像冲沟解译方法相比较，本发明使遥感二维图像具有符合人类视觉习惯的地形信息、沟谷的遥感图像特征清晰、解译沟缘线精度高的优点。
【专利说明】一种基于地形和遥感影像融合技术提取冲沟方法
[0001]【技术领域】:
本发明涉及遥感【技术领域】，具体地讲是一种将地形信息和遥感影像融合技术应用于冲沟沟缘线的提取。该发明可用于基于光学遥感图像，进行冲沟解译、沟蚀研究及制图场合。
[0002]【背景技术】:
冲沟侵蚀是一种重要的土壤侵蚀方式，它不但是导致土地退化的一个重要过程，还是江河泥沙的一个主要来源。目前国内外学者主要利用地面手工测量法(Casali et al.，2006; Vandekerckhove et al., 2001 ；)、侵蚀针(桩)监测法(1nita et al., 2006 ；Martinez et al.2003)、RTK_GPS 测量法(何福红 et al., 2005; Yougqiuffu et al.，2005)、航片和数字高程模型(Martinez et al.,2003 ；Harley et al.，1999)等方法来定量研究沟蚀的发生演变规律。自20世纪60年代初航空摄影技术出现以来，遥感技术已经被广泛应用于土壤侵蚀、冲沟沟蚀研究等领域。利用航天遥感图像解译土壤侵蚀的研究可以追溯到1940年(Smith et al.，1943)。随后的数十年里，诸多学者利用航空照片及摄影测量技术对冲沟侵蚀进行研究(Nachtergaele and Poesen, 1999; Betts and DeRose, 1999 ；Martinez-Casasnovas et al., 2003)。研究冲沟的方式是对遥感图像，目视解译冲沟沟缘线，由于沟缘线是一条重要的坡度分界线，其上部为坡度较为平缓的沟间地类型，其下部是坡度较为陡峻的沟谷地类型，而遥感二维图像提供了的三维地形信息量不足以准确解译出冲沟沟缘线，因而精度较低。另外，对北半球来说，由于资源卫星大多为太阳同步极轨卫星，其成像时间为地方时十点半。成像时，阳光从东南向射入，山脊两侧的南向坡形成光照面，北向坡为阴影面，而传统的上北下南的影像构图方式，使得阴影面位于光照面上方，这就形成了视觉上的反立体现象(Saraf et al., 1996; Rudnicki, 2000; Patterson, 2004,Boffu et al, 2012)。反立体现象,使得遥感图像中的冲沟在视觉上表现为凸起山脊,而山脊则表现为凹陷的冲沟，增加了准确判读冲沟的难度。因而单纯基于遥感图像解译冲沟信息效果有限。
[0003]数字高程模型(DEM)是地表三维地形的数字表达方式，是用一组有序数值阵列形式表示地面高程的一种实体地面模型能，是计算诸如坡度、坡向、坡度变化率等因子的数学模型。但利用DEM数据难以精确解译冲沟沟缘线。目前，基于DEM数据计算机自动识别冲沟沟缘线的方法，但由于数据源、地貌类型的多样性及沟缘线的多级性和层次性均导致自动识别沟缘线算法精度和通用性差。
[0004]因此，将符合人类视觉习惯的三维地形信息融合到遥感图像中，将增加遥感二维图像的三维地形信息，且使得图像空间几何关系及纹理特征符合人类视觉习惯，有助于精确识别冲沟；并辅以由DEM数据计算的冲沟沟缘线坡度阈值数据，将有助于提高解译冲沟沟缘线的精度。目前三维地形信息同遥感图像相融合的方法主要有图像南北向旋转法、像元值逆转法、SRM (Shaded Relief Model)参与的HIS融合3类。像元值逆转法和SRM参与的HIS融合法通过修改图像像素灰度值来实现正立体校正，但均存在很明显的光谱信息损失。图像南北向旋转法尽管不存在光谱信息损失问题，但构图方式采用上南下北构图方式，导致较难推广使用。[0005]
【发明内容】
:
本发明的目的是克服上述已有技术的不足，而提供一种基于地形和遥感影像融合技术提取冲沟方法，主要解决SRM参与的HIS融合法方法有明显的遥感图像光谱信息损失的问题。
[0006]本发明的技术方案是:一种基于地形和遥感影像融合技术提取冲沟的方法，其特殊之处在于，包括以下步骤:
(I )将图幅相同的光学遥感红、绿、蓝三波段(RGB三波段)图像和DEM (数字高程模型)数据的坐标系统一成高斯-克吕格平面直角坐标系，并进行重采样，使两者具有相同的空间分辨率；
(2 )将光学遥感红、绿、蓝三波段图像(RGB色彩空间)，变换至彩色亮度(Intensity)、色相(Hue)和饱和度(Saturation)空间(HIS色彩空间)，获得彩色亮度i色相If和饱和
度S三分量数据；
(3 )根据遥感图像成像时刻(年月日时分秒)及图像中心点经纬度坐标，计算遥感图像成像时刻的太阳方位角< 和太阳高度角a;
(4 )基于DEM数据，计算地表粗糙度，获得地表粗糙度数据SJ，并之进行线性标准化处理(式2)，(高斯-克吕格平面直角系)获得线性标准化地表粗糙度数据(式3或式4); (5)基于DEM数据，以遥感图像成像时刻太阳方位角名加,τ (180° )和遥感图像成
像时刻太阳高度角β为参数计算研究区太阳西北方向地形阴影图.SM/(Shaded ReliefModel)；
(6 )太阳西北方向地形阴影图SRV和彩色亮度J分量数据分别以线性标准化地表粗糙度Cf和I—fT作为权重，合并为新的彩色亮度分量数据(式I)；r ^SRM - {1-?Ι') ⑴；
(7 )将(2 )获得的色相I和饱和度S.以及(6 )获得的新亮度J.三分量数据进行HIS->RGB变换，反变换至RGB色彩空间；该步骤实现了将地形信息同遥感图像相融合，形成地形信息同遥感信图像融合图像；
(8 )基于DEM数据(高斯-克吕格平面直角系)，计算每个像素的坡度，获得地形坡度数
据 ;
(9 )野外对冲沟沟缘线坡度值进行实测和统计，确定冲沟沟缘线坡度最小阈值数据S..!TSfl ，
(I O)根据野外实测冲沟沟缘线坡度最小阈值S.，对地形坡度数据Sfcif进行阈值分
害I]，将冲沟沟缘线分割出来，获得冲沟沟缘线图斑数据；
(I I)坡度阈值切割获得的冲沟沟缘线图斑数据，结合地形信息同遥感信图像相融合的图像，基于目视解译方法，精确解译冲沟沟缘线。
[0007]进一步的，所述的遥感数据(图像)为航天、航空遥感图像中的可见光波段(红光、绿光和蓝光)。[0008]进一步的，所述的地形坡度数据和地表粗糙度数据是基于DEM数据来计算的；从地形学角度出发，将地面凹凸不平的程度定义为粗糙度，一般定义为地表单元曲面面积与投影面积之比(式2);地表粗糙度反映了地表地形复杂程度，取值大于等于I ;地表坡度为0，即平坦地区的地表粗糙度值为1，而地表坡度大，即区域(山岭)地表粗糙度值高；线性标准化地表粗糙度Cf是由(式3)或(式4)确定，Cf取值(0，I ]。:
【权利要求】
1.一种基于地形和遥感影像融合技术提取冲沟的方法，其特征在于，包括以下步骤: (I )将图幅相同的光学遥感红、绿、蓝三波段(RGB三波段)图像和DEM (数字高程模型)数据的坐标系统一成高斯-克吕格平面直角坐标系，并进行重采样，使两者具有相同的空间分辨率； (2)将光学遥感红、绿、蓝三波段图像(RGB色彩空间)，变换至彩色亮度(Intensity)、色相(Hue)和饱和度(Saturation)空间(HIS色彩空间)，获得彩色亮度J、色相孖和饱和度S三分量数据； (3 )根据遥感图像成像时刻(年月日时分秒)及图像中心点经纬度坐标，计算遥感图像成像时刻的太阳方位角#和太阳高度角α; (4 )基于DEM数据，计算地表粗糙度，获得地表粗糙度数据20,并之进行线性标准化处理，获得线性标准化地表粗糙度数据范.; (5)基于DEM数据，以遥感图像成像时刻太阳方位角I加;r (180° )和遥感图像成像时刻太阳高度角β为参数计算研究区太阳西北方向地形阴影图5ΚΛ/(Shaded ReliefModel)； (6 )太阳西北方向地形阴影图5ΚΛ/和彩色亮度分量数据分别以线性标准化地表粗糙度和1-2T作为权重,合并为新的彩色亮度分量数据(式I)；
I=ω - a1 M SRM - (1- Err) * I (I); (7 )将(2 )获得的色相Jf和饱和度S以及(6 )获得的新亮度厂三分量数据进行HIS->RGB变换，反变换至RGB色彩空间；该步骤实现了将地形信息同遥感图像相融合，形成地形信息同遥感信图像融合图像； (8 )基于DEM数据，计算每个像素的坡度，获得地形坡度数据Shm ； (9 )野外对冲沟沟缘线坡度值进行实测和统计，确定冲沟沟缘线坡度最小阈值数据S ^.nun ， (1O)根据野外实测冲沟沟缘线坡度最小阈值Sraut,对地形坡度数据进行阈值分害I]，将冲沟沟缘线分割出来，获得冲沟沟缘线图斑数据； (11)坡度阈值切割获得的冲沟沟缘线图斑数据，结合地形信息同遥感信图像相融合的图像，基于目视解译方法，精确解译冲沟沟缘线。
2.根据权利要求1所述的一种基于地形和遥感影像融合技术提取冲沟方法，其特征在于，所述的遥感数据(图像)为航天、航空遥感图像中的可见光波段(红光、绿光和蓝光)。
3.根据权利要求1所述的一种基于地形和遥感影像融合技术提取冲沟方法，其特征在于，所述的地形坡度数据和地表粗糙度数据是基于DEM数据来计算的(式2): m - l/cos(Slope_rad)(2) 其中，Slqje_rad为坡度数据，以弧度为单位； 线性标准化地表粗糙度是由式(3)确定:CTi = tix/maxicr)(3) 其中，iimx(m.)为ar最大管j i取值(O，I]； 式(3)亦可以采用式(4)来完成线性标准化:
ετ'* = (err -miii(Ern/max(07)(4)。
4.根据权利要求1所述的一种基于地形和遥感影像融合技术提取冲沟方法，其特征在于，所述的HIS色彩空间 向RGB色彩空间变换过程中，新的彩色亮度/ *分量数据取值是太阳西北方向地形阴影图5ΚΛ/和原始彩色亮度/分量数据分别以粗糙度和1-+CH乍为权重进行取值；即由式(I)来确定。
5.根据权利要求1所述的一种基于地表粗糙度和阴影模型的遥感图像反立体校正方法，其特征在于，所述HIS色彩空间包括HSV色彩空间、HLS色彩空间；所述RGB色彩空间向HIS色彩空间变换，在有的遥感图像处理软件系统中又被称为HSV变换、HLS变换、USGSMunsell HSV变换；在HIS色彩空间中彩色亮度I分量同HLS色彩空间中的颜色明亮度(Lightness)分量、HSV色彩空间中的色彩明度(Value)分量、USGS Munsell HSV色彩空间中的色彩明度(Value)分量相对应，均表征色彩强弱(明亮)，彩色亮度I(Intensity)分量及公式(3)中的J,分别等价为HLS色彩空间中的颜色明亮度(Lightness)分量、HSV色彩空间中的色彩明度(Value )分量。
6.根据权利要求1所述的一种基于地形和遥感影像融合技术提取冲沟方法，其特征在于，所述的冲沟沟缘线图斑数据是基于野外实测沟缘线处坡度确定冲沟沟缘线坡度阈值Skei ,通过对地形坡度数据进行阈值分割而获得，其形态特征零散、破碎，不便于目视解译冲沟。
【文档编号】G01C11/00GK103940407SQ201410049422
【公开日】2014年7月23日 申请日期:2014年2月13日 优先权日:2014年2月13日
【发明者】王涛, 何福红, 蒋卫国, 顾丽娟, 张振华 申请人:鲁东大学