一种利用全极化雷达数据提取土壤含水信息的方法技术

本发明专利技术属于雷达数据遥感分析领域，具体公开一种利用全极化雷达数据提取土壤含水信息的方法，该方法具体包括以下步骤：步骤(1)对全极化雷达数据不同极化方式的数据分别进行预处理；步骤(2)对上述步骤(1)处理后的雷达图像进行RGB彩色合成；步骤(3)上述步骤(2)中得到的RGB彩色图像进行运算，获得土壤相对含水信息图。该方法对全极化雷达数据的不同极化方式下的回波数据进行运算和处理，获取土壤的相对含水量，提取土壤中富水带信息，从而更好地服务于农业生产及地质找水等方向。

【技术实现步骤摘要】

本专利技术属于雷达数据遥感分析领域，具体涉及一种利用全极化雷达数据提取土壤含水信息的方法。
技术介绍
土壤水分是指土壤非饱和层(也称其为渗流层)的水分含量，其能够反映不同地区浅地表的干旱程度，快速地进行大区域的土壤含水信息，不仅可以指导农业生产和水利灌溉，也可以为地质灾害监测和预报提供参考，同时，也可以对部分矿产(如砂岩型铀矿、膏盐)的找矿作用提供帮助。基于卫星遥感技术对土壤水分的时空分布信息进行反演，一直是遥感研究的热点之一。目前的卫星遥感反演土壤含量技术根据遥感测量手段的不同可分为：光学遥感、微波遥感2类。光学遥感主要利用土壤表面光谱反射特性、土壤表面发射率及表面温度来估算土壤水分，其空间分辨率高，可供选择的卫星传感器多，并可提供高光谱数据，但光学方法多建立各种指数与土壤水分之间的经验关系，基于蒸散模型的方法比较复杂，需要输入的参数较多，获取困难，且这些方法属于间接反演土壤水分的方法，其应用有一定的局限性。微波遥感主要建立雷达后向散射系数与土壤介电常数的关系，进而反演出土壤水分，其空间分辨率高，能穿透植被，到达地表5cm左右深度，且能全天候测量，低频波段雷达后向散射系数对土壤水分高度敏感，而在高频波段雷达后向散射系数对植被敏感，L波段反演土壤水分能取得较好的效果。然而，目前的大部分算法都是基于单极化雷达数据，无论是经验模型，半经验模型，还是物理模型，其都基于雷达后向散射系数与土壤含水量的关系，但由于雷达后向散射对地表的几何结构高度敏感，地面的粗糙度和地形起伏都会对反演的土壤含水信息造成干扰。另一方面，相同含水性的土壤在不同极化方式下将产生不同的雷达回波特征，通常情况含水的土壤将在HH和VV极化方式上呈现较强的回波强度，而在VH极化方式下却呈现较弱的回波强度，这表明，可以利用雷达数据的不同极化方式提取土壤含水信息，尤其是在干旱地区寻找浅层地下富水带。如何通过对全极化雷达数据中不同极化方式信息的处理和合成，反演土壤含水信息，提取土壤富水带，减少地形起伏和地表粗糙度的干扰，是亟待解决的问题。
技术实现思路
本专利技术的目的在于提供一种利用全极化雷达数据提取土壤含水信息的方法，该方法对全极化雷达数据的不同极化方式下的回波数据进行运算和处理，获取土壤的相对含水量，提取土壤中富水带信息，从而更好地服务于农业生产及地质找水等方向。实现本专利技术目的的技术方案：一种利用全极化雷达数据提取土壤含水信息的方法，该方法具体包括以下步骤：步骤(1)对全极化雷达数据不同极化方式的数据分别进行预处理；步骤(2)对上述步骤(1)处理后的雷达图像进行RGB彩色合成；步骤(3)上述步骤(2)中得到的RGB彩色图像进行运算，获得土壤相对含水信息图。所述的步骤(1)具体包括如下步骤：(1.1)对雷达系统的RAW数据进行聚焦处理，获得雷达SLC图像；(1.2)对上述步骤(1.1)中得到的SLC图像进行多视处理，获得雷达强度图像；(1.3)对上述步骤(1.2)中得到的雷达强度图像进行几何校正和辐射校正，获取地理编码准确，后向散射系数准确的雷达强度图像，将不同极化方式的雷达图像以不同的文件名命名。所述的步骤(2)具体包括如下步骤：对不同极化的雷达强度图像进行RGB(红绿蓝)彩色合成，其中R＝HH极化，G＝VH极化方式，B＝VV极化方式，以文件名RGB存储图像。所述的步骤(3)具体包括如下步骤：(3.1)对RGB彩色图像进行运算，运算表达式为255*(float(b1)+float(b3))/(float(b1)+(float(b2)+float(b3))；(3.2)对上述步骤(3.1)中得到的运算结果进行密度分割，获得土壤相对含水量的信息。所述的步骤(3.1)中的b1＝R(红色)，b2＝G(绿色)，b3＝B(蓝色)。本专利技术的有益技术效果：(1)本专利技术因为利用的是含水土壤在不同极化方式下的回波特征，因而减小了地表几何特征对反演结果的干扰(2)本专利技术借助于全极化雷达数据，信息丰富，算法却较为简单。本专利技术可应用于雷达遥感领域，服务于农业生产和地质找矿等方向。附图说明图1为本专利技术所提供的一种利用全极化雷达数据提取土壤含水信息的方法的流程图。具体实施方式下面结合附图和实施例对本专利技术作进一步详细说明。如图1所示，本专利技术所提供的一种利用全极化雷达数据提取土壤含水信息的方法，该方法具体包括以下步骤：步骤(1)对全极化雷达数据不同极化方式的数据分别进行预处理。步骤(1.1)对雷达系统的原始(RAW格式)数据进行聚焦处理，获得雷达SLC图像利用雷达数据处理软件(如SARscape)打开数据文件，对雷达系统的RAW格式原始数据()进行聚焦处理，获得雷达单视复型数据(SLC数据)。以新疆乌伦古河地区全极化Alos-2雷达数据为实例，上述打开的数据文件为日本2014年发射升空的Alos-2卫星搭载传感器拍摄新疆乌伦古河地区HH、VH和VV三种极化方式的雷达数据。HH、VH、HV、VV是雷达遥感系统常用四种极化方式。雷达发射的能量脉冲的电场矢量，可以在垂直或水平面内被偏振。无论哪个波长，雷达信号可以传送水平(H)或者垂直(V)电场矢量，接收水平(H)或者垂直(V)或者两者的返回信号。HH极化方式即当雷达以水平电场矢量发射信号同时接收水平的返回信号，VH即以垂直电场矢量发射信号同时接收水平返回的信号，HV即以水平电场矢量发射的信号同时接收垂直返回的信号，VV即以垂直电场矢量发射信号同时接收垂直返回的信号。步骤(1.2)对上述步骤(1.1)中得到的SLC数据进行多视处理，获得雷达强度图像；利用SARscape软件对雷达SLC图像进行多视处理，抑制雷达数据的斑点噪声，获得雷达强度图像。步骤(1.3)对上述步骤(1.2)中得到的雷达强度图像进行几何校正和辐射校正，获取地理编码准确，后向散射系数准确的雷达强度图像，将不同极化方式的雷达图像以不同的文件名命名；HH极化方式记作HH，VH极化方式记作VH，VV极化方式记作VV。例如，下载新疆乌伦古河地区的DEM(数字高程模型)数据，利用SARscape软件对雷达强度图像进行几何校正和辐射校正，获取地理编码准确，后向散射系数准确的三景雷达强度图像；三景图像文件分别以img格式储存，三景图像文件名分别记作Wulunhe-Alos-HH.img，Wulunhe-Alos-VH.img和Wulunhe-Alos-VV.img。步骤(2)对上述步骤(1.3)处理后的雷达图像进行RGB彩色合成对不同极化的雷达强度图像进行RGB彩色合成，其中R＝HH极化，G＝VH极化方式，B＝VV极化方式，以文件名RGB存储图像。例如，利用ENVI软件打开上述新疆乌伦古河地区的三景图像，对这三景图像进行RGB彩色合成，其中R＝Wulunhe-Alos-HH.img，G＝Wulunhe-Alos-VH.img，B＝Wulunhe-Alos-VV.img，获得RGB彩色合成图像，文件名记作Wulunhe-Alos-RGB.img。步骤(3)上述步骤(2)中得到的RGB彩色图像进行运算，获得土壤相对含水信息图步骤(3.1)对RGB彩色图像进行运算，运算表达式为255*(float(b1)+float(b3))/(float(b1)+(float(b2)+float本文档来自技高网...

【技术保护点】
一种利用全极化雷达数据提取土壤含水信息的方法，其特征在于，该方法具体包括以下步骤：步骤(1)对全极化雷达数据不同极化方式的数据分别进行预处理；步骤(2)对上述步骤(1)处理后的雷达图像进行RGB彩色合成；步骤(3)上述步骤(2)中得到的RGB彩色图像进行运算，获得土壤相对含水信息图。

【技术特征摘要】
1.一种利用全极化雷达数据提取土壤含水信息的方法，其特征在于，该方法具体包括以下步骤：步骤(1)对全极化雷达数据不同极化方式的数据分别进行预处理；步骤(2)对上述步骤(1)处理后的雷达图像进行RGB彩色合成；步骤(3)上述步骤(2)中得到的RGB彩色图像进行运算，获得土壤相对含水信息图。2.根据权利要求1所述的一种利用全极化雷达数据提取土壤含水信息的方法，其特征在于：所述的步骤(1)具体包括如下步骤：(1.1)对雷达系统的RAW数据进行聚焦处理，获得雷达SLC图像；(1.2)对上述步骤(1.1)中得到的SLC图像进行多视处理，获得雷达强度图像；(1.3)对上述步骤(1.2)中得到的雷达强度图像进行几何校正和辐射校正，获取地理编码准确，后向散射系数准确的雷达强度图像，将不同极化方式的雷达图像以不同的文件名命名。3.根据权利...

【专利技术属性】
技术研发人员：孟树，叶发旺，张川，刘洪成，郭帮杰，
申请(专利权)人：核工业北京地质研究院，
类型：发明
国别省市：北京;11