发明名称 |
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>ϵ</mi><mi>cos</mi><mi>θ</mi><mo>-</mo><msqrt><mi>ϵ</mi><mo>-</mo><msup><mi>sin</mi><mn>2</mn></msup><mi>θ</mi></msqrt></mrow><mrow><mi>ϵ</mi><mi>cos</mi><mi>θ</mi><mo>+</mo><msqrt><mi>ϵ</mi><mo>-</mo><msup><mi>sin</mi><mn>2</mn></msup><mi>θ</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>θ</mi><mo>-</mo><msqrt><mi>ϵ</mi><mo>-</mo><msup><mi>sin</mi><mn>2</mn></msup><mi>θ</mi></msqrt></mrow><mrow><mi></mi><mi>cos</mi><mi>θ</mi><mo>+</mo><msqrt><mi>ϵ</mi><mo>-</mo><msup><mi>sin</mi><mn>2</mn></msup><mi>θ</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号 |