发明名称 GNSS卫星反射信号用于土壤介电常数的精确反演方法
摘要 本发明公开了一种GNSS卫星反射信号用于土壤介电常数的精确反演方法,属于微波遥感技术领域。本发明方法利用GNSS双基雷达接收GNSS卫星的直达信号和反射信号,获得直达信号功率和反射信号功率,进而获得土壤反射率R;获得地表反射系数与相对土壤介电常数的关系;然后根据接收的不同的反射信号,确定对应的反射系数和土壤反射率,设置中间参数,将土壤反射率的式子转化为一元四次方程,求解方程的解,确定土壤介电常数。本发明方法应用范围广,反演精度高,能精确反演土壤介电常数,进而得到准确的土壤湿度。
申请公布号 CN103336104A 申请公布日期 2013.10.02
申请号 CN201310227524.X 申请日期 2013.06.08
申请人 北京航空航天大学 发明人 孙国良;白永星;刘嘉俊
分类号 G01N33/24(2006.01)I 主分类号 G01N33/24(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 祗志洁
主权项 1.一种GNSS卫星反射信号用于土壤介电常数的精确反演方法,其特征在于,包括如下步骤:(1)利用GNSS双基雷达接收GNSS卫星的直达信号和反射信号,获得直达信号功率P<sub>d</sub>和反射信号功率P<sub>r</sub>,然后根据土壤反射率R的定义得到:<img file="FDA00003323881700011.GIF" wi="198" he="144" />(2)得到地表反射系数与相对土壤介电常数ε的关系:<maths num="0001"><![CDATA[<math><mrow><msub><mi>R</mi><mi>v</mi></msub><mo>=</mo><mfrac><mrow><mi>&epsiv;</mi><mi>cos</mi><mi>&theta;</mi><mo>-</mo><msqrt><mi>&epsiv;</mi><mo>-</mo><msup><mi>sin</mi><mn>2</mn></msup><mi>&theta;</mi></msqrt></mrow><mrow><mi>&epsiv;</mi><mi>cos</mi><mi>&theta;</mi><mo>+</mo><msqrt><mi>&epsiv;</mi><mo>-</mo><msup><mi>sin</mi><mn>2</mn></msup><mi>&theta;</mi></msqrt></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>R</mi><mi>h</mi></msub><mo>=</mo><mfrac><mrow><mi></mi><mi>cos</mi><mi>&theta;</mi><mo>-</mo><msqrt><mi>&epsiv;</mi><mo>-</mo><msup><mi>sin</mi><mn>2</mn></msup><mi>&theta;</mi></msqrt></mrow><mrow><mi></mi><mi>cos</mi><mi>&theta;</mi><mo>+</mo><msqrt><mi>&epsiv;</mi><mo>-</mo><msup><mi>sin</mi><mn>2</mn></msup><mi>&theta;</mi></msqrt></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>R<sub>v</sub>是光滑表面的垂直反射系数,R<sub>h</sub>是光滑表面的水平反射系数,θ是GNSS卫星的入射角;(3)根据接收的反射信号,确定对应的反射系数和土壤反射率,土壤反射率为用R<sub>v</sub>和R<sub>h</sub>表达的式子;然后设置中间参数c和x:c=tan<sup>2</sup>θ,<img file="FDA00003323881700014.GIF" wi="377" he="89" />将土壤反射率的式子用中间参数表示,并转化为一元四次方程;当只接受反射信号的左旋圆极化分量、或者接收反射信号的左旋圆极化分量和右旋圆极化分量、或者只考虑反射信号的左旋圆极化分量中的垂直极化分量时,求解一元四次方程的第一个解;当只接受反射信号的右旋圆极化分量时,求解一元四次方程的第二个解;最后,确定对应的土壤介电常数Permittivity:<maths num="0003"><![CDATA[<math><mrow><mi>Permittivity</mi><mo>=</mo><mfrac><mrow><msup><mi>X</mi><mn>2</mn></msup><mo>+</mo><mi>c</mi></mrow><mrow><mn>1</mn><mo>+</mo><mi>c</mi></mrow></mfrac></mrow></math>]]></maths>其中,X表示所求解的第一个解或第二个解。
地址 100191 北京市海淀区学院路37号