微波辐射计监测土壤水分的方法技术

本发明专利技术用于微波遥感技术应用领域，特别是微波辐射计监测土壤水分的方法。其步骤为：第一步是准备运行环境和辐射计数据；第二步是反演得到观测点地表介电常数并将其栅格化；第三步是根据介电常数计算得到土壤含水量；第四步是将土壤含水量分布图成图并保存为常见的图像格式。

【技术实现步骤摘要】

本专利技术用于微波遥感技术应用领域领域，特别是利用微波辐射计数据反演土壤水分的方法。
技术介绍
现有利用微波辐射计反演地表土壤水分的方法中，计算地表植被层不透明度时，需要关于地面植被类型分类、植被层含水量的辅助数据，计算公式为τc＝b*ω/cosθ其中b为与植被冠层结构、频率有关的参数，ω为植被层含水量，θ为观测角度。b和ω都需要通过辅助数据得到的。但是在实际应用中，这些辅助数据往往不易获取或者不可靠，即使获得了这部分数据在应用时也比较复杂(如图2)。现有的反演流程中，根据地表植被分类辅助数据计算植被影响的流程为(如图2)。现有技术在计算植被层不透明度时存在的最大的问题是由于辅助数据难以及时获取和更新，只能以经验数值代入运算，无法保证辅助数据与微波辐射计观测数据的时间同步性。例如，参数b是根据地面植被分类数据确定的，但是实际中这一参数会随着植被生长变化而改变，以不同的频率观测这个数值也存在差异；ω作为植被层含水量更是会随着季节、降水的变化而变化。经验计算公式中的b和ω数值都具有时间、气候变化性，用固定的数值来描述辐射计观测时植被层不透明度不能完全反映辐射计观测时刻植被对地表辐射的影响状况。
技术实现思路
本专利技术的目的在于提供。本专利技术提出了在缺少同时相地面植被辅助数据信息的情况下，成功实现利用微波辐射计数据反演得到地表土壤水分的技术流程路线。本专利技术利用多通道之间微波辐射差异实时的对观测点植被覆盖茂密状况做出判断，当存在植被影响时，利用微波辐射亮度值来计算得到实时的地表覆盖的植被的影响，从而实现中等覆盖程度、稀疏覆盖程度的植被地区和裸露地表地区的地表土壤水分的反演。新的流程改进了原有反演流程当中对大量的地表植被辅助数据信息的依赖，并计算观测时刻的地表植被影响。本专利技术可提供一种为在全球范围内低空间分辨率的土壤水分变化的研究与监测工作服务。附图说明图1是本专利技术微波辐射计数据反演土壤水分的方法流程图。图2是计算植被层不透明度的现有技术流程图。图3是计算植被层不透明度的本专利技术技术流程图。具体实施例方式图1的微波辐射计数据反演土壤水分的方法，其过程如下反演工作以微型计算机为物理平台，以微软视窗系统为软件系统环境。操作流程可分为四个步骤第一步是准备运行环境和辐射计数据；第二步是根据观测点地理坐标数据建立投影变换索引表，计算观测点处的植被层不透明度参数，反演运算出地表土壤介电常数；第三步是根据介电常数反演计算得到土壤含水量；第四步是土壤含水量分布图成图并保存为常见的图像文件格式。第一步准备运行环境和辐射计数据。首先接通电源，启动计算机硬件平台进入视窗操作系统软件环境。如果辐射计数据以光盘或者磁带为载体保存，计算机硬件平台需要配备有光盘驱动器或者磁带机，作为辐射计数据导入装置。然后将存储在光盘介质(CD/DVD)(或者磁带)上的微波辐射计测量数据整理到计算机本地硬盘上来。以光盘驱动器为例，首先按键弹出光盘托盘，放入数据光盘，按键收回光盘托盘。在视窗操作系统的光盘驱动器所对应的盘符中即可找到微波辐射计的数据，从中选择进行反演运算的辐射计数据并拷贝到计算机本地磁盘上。这部分是进行反演操作流程的的准备过程。第二步是根据观测点地理坐标数据建立投影变换索引表，并且计算观测点处的植被层不透明度参数，反演运算出地表土壤介电常数。首先从微波辐射计数据当中读取得到观测点地理坐标数据(经度坐标数据和纬度坐标数据)和相同频率下两种极化方式的亮度温度数据(水平极化TBh和垂直极化亮温TBv)，这部分工作可以利用为数据源专门开发的一些工具来完成，以AMSR/E数据为例NASA提供了hdp工具软件，也可以利用ENVI中File-＞Open Image File模块完成。然后分别根据观测点地理坐标数据建立投影变换索引表和根据亮度温度数据计算观测点处植被层不透明度、反演地表土壤介电常数。这一步工作可以利用ENVI的扩展功能模块和IDL环境编程实现。1.根据经度坐标数据和纬度坐标数据，建立观测点地理坐标数据投影变换到栅格的索引表文件。在ENVI中，利用File-＞Open Image File模块把所有观测点的经度数据和纬度数据分别作为两个通道的数据读入，建立索引表时，以经度通道数据作为x轴数据输入，以纬度通道数据作为y轴数据输入。生成的栅格主要参数椭球体选择WGS-84椭球体，坐标系统选择为大地坐标系统，输出像素的大小为0.083度，输出的翻转角度为0度。生成索引表文件的工作可以利用ENVI的Map-＞Groreferencefrom Input Geometry-＞Build GLT模块完成，生成的索引表文件以。GLT为扩展名保存在本地硬盘。2.根据连个极化方式通道的亮度温度数据计算MPDI和植被层不透明度、反演计算地表土壤介电常数，这部分工作可以利用ENVI的波段运算和IDL环境编程实现。1)读取水平极化和垂直极化亮温数据。利用利用File-＞Open Image File模块把所有观测点的水平极化方式和垂直极化方式的亮度温度数据作为TBh和TBv两个通道数据读入。2)计算每个观测点的微波极化差异指数(MPDI)。MPDI的定义为MPDI=TBv-TBhTBv+TBh]]>计算微波极化差异指数(MPDI)根据TBv和TBh两个通道数据，利用ENVI的Basic Tools-＞Band Math模块完成。输入运算表达式时，输入表达式为(b1-b2)/(b1+b2)。选择运算波段的时候，TBv数据通道和MPDI定义式中的b1相对应，TBh数据通道与MPDI定义式中的b2相对应，计算出所有观测点处的微波极化差异指数点数据。3)根据MPDI值确定观测点植被层不透明度。遍历所有观测点，根据每个观测点上的MPDI数值判断观测点处植被覆盖状况，对于浓密植被覆盖地区或者裸露地表、稀疏植被覆盖地区的植被不透明度赋予合理数值，对于中等植被覆盖地区需要根据MPDI值计算植被层的不透明度a.当MPDI＜0.017时，判断观测点为浓密植被覆盖地区，这些地区的植被衰减影响过而大无法消除，观测得到的亮温数据对于地表含水量已经不再敏感，植被层不透明度τc取为-1，用负数表示植被覆盖浓密，无法反演得到地表土壤介电常数；b.当观测点MPDI＞0.026的时候，判断观测点为裸露地表或者稀疏植被覆盖地区，此时不需要考虑植被的衰减影响，植被层不透明度τc取为0值，表示在这些点上不存在或者可以忽略植被的衰减影响；c.当0.017＜MPDI＜0.026时，判断观测点为中等植被覆盖度地区，根据MPDI计算植被层不透明度τcτc=-0.223ln(MPDI/k2)-1.239MPDI-0.0010.12+0.00085k2-0.0547k+0.05411]]>此时植被的影响不能忽略，以MPDI计算植被不透明度，作为接下来的反演土壤介电常数运算的参数；4)代入植被层不透明度，反演地表土壤介电常数。在每个观测点上根据微波辐射传输方程Tbp=Te]]>正演计算观测到的水平和垂直极化亮度温度数据。正演得到的水平和垂直极化亮温与观测值之间的方差最小的时候，对应的地表土壤介电常数就是反演得到的地表土壤介电常数。当以C波段观测的时候，大气的衰减作用可以忽略，不再计算大本文档来自技高网...

【技术保护点】
一种微波辐射计监测土壤水分的方法，其步骤包括：第一步是准备运行环境和辐射计数据；　第二步是建立投影变换索引表，计算观测点处的植被层不透明度参数，反演运算得到地表土壤介电常数；第三步是根据介电常数反演计算得到土壤含水量 ；第四步是土壤含水量分布图成图并保存为常见的图像文件格式。

【技术特征摘要】

【专利技术属性】
技术研发人员：李震，王磊，
申请(专利权)人：中国科学院遥感应用研究所，
类型：发明
国别省市：11[中国|北京]