CN111239209B - 机会信号反射遥感的全极化单次反射仿真方法及系统 - Google Patents

机会信号反射遥感的全极化单次反射仿真方法及系统 Download PDF

Info

Publication number
CN111239209B
CN111239209B CN202010095663.1A CN202010095663A CN111239209B CN 111239209 B CN111239209 B CN 111239209B CN 202010095663 A CN202010095663 A CN 202010095663A CN 111239209 B CN111239209 B CN 111239209B
Authority
CN
China
Prior art keywords
layer
reflection
polarization
dielectric constant
soil
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010095663.1A
Other languages
English (en)
Other versions
CN111239209A (zh
Inventor
吴学睿
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Astronomical Observatory of CAS
Original Assignee
Shanghai Astronomical Observatory of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Astronomical Observatory of CAS filed Critical Shanghai Astronomical Observatory of CAS
Priority to CN202010095663.1A priority Critical patent/CN111239209B/zh
Publication of CN111239209A publication Critical patent/CN111239209A/zh
Application granted granted Critical
Publication of CN111239209B publication Critical patent/CN111239209B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/22—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance
    • G01N27/221—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance by investigating the dielectric properties
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/22—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance
    • G01N27/223—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance for determining moisture content, e.g. humidity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24—Earth materials
    • G01N33/246—Earth materials for water content
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing

Landscapes

  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Electrochemistry (AREA)
  • Electromagnetism (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明提供了一种机会信号反射遥感的全极化单次反射仿真方法，包括：将土壤剖面划分为多个分层，并根据输入参数建立介电常数模型，得到各个分层的介电常数；利用介电常数来获取分层的透过率、反射率和衰减；分别构建光滑和粗糙表面反射率矩阵，利用极化合成的方法计算各极化的反射特性；根据透过率、反射率、衰减及反射特性来获得总接收机波形。本发明还提供了相应的仿真装置。本发明的全极化单次反射仿真方法通过计算土壤各分层的介电常数，并采用极化合成方法，以应用在机会信号反射遥感中，对根区土壤水分建立相应的仿真，计算各种极化的波形信息，大大提高对目标各种信息的获取能力，且通过计算镜像散射和非镜像的漫散射，使计算结果更准确。

Description

机会信号反射遥感的全极化单次反射仿真方法及系统
技术领域
本发明属于陆地和海洋参数测量领域，具体涉及一种机会信号反射遥感的全极化单次反射仿真方法及系统。
背景技术
根区土壤水分处于土壤-植物-大气连续体的核心位置，对干旱半干旱地区作物生长、植被恢复和水土流失过程均具有重要影响。根区土壤水分的获取对于农业以及粮食生产尤为重要，同时土壤水分剖面信息的获取能与当前的陆面/水文模式相契合，更好地支撑数据同化研究、天气与气候预报、水资源管理等应用需求。
微波遥感技术是获取土壤水分的有效监测手段。目前合成孔径雷达和微波辐射计是两种较为常见和成熟的微波遥感方式。土壤水分和海洋盐度观测卫星(SMOS)和土壤水分探测卫星(SMAP)都是工作在L波段的微波遥感卫星。但在L波段，只能穿透浅层地表，因此适用于土壤表层(约5cm)的土壤水分获取，而且由于穿透强度的限制，L波段只能穿透中等厚度的植被冠层。
对于根区土壤水分的获取，则需要使传感器工作在更低的频率范围，如频率范围是240-380MHz的P波段。其中，P波段可以穿透深度约40cm,同时可以很好的规避茂密植被的影响。2012年美国在开展了机载P波段实验，即AirMOSS(Airborne Microwave Observatoryof Sub-canopy and Subsurface)试验。但是P波段的星载观测则面临着天线尺寸过大、RFI(radio-frequency interference)影响和频段分配竞争等挑战。
机会信号反射遥感为P波段根区土壤水分获取提供了新的机遇。机会信号反射遥感是一种现在存在的遥感技术手段，它利用现有的通讯卫星或者导航卫星作为信号源，采用专门的反射信号接收机来接收机会卫星信号经地表反射的信号，对地物参数进行遥感监测。由于使用数字通信卫星作为信号源，不需要研制专门的发射机，遥感过程中，只需要研制反射信号接收机接收地表的反射信号，因此，这种遥感手段相比于不采用机会信号反射遥感的P波段的雷达，可以克服天线尺寸过大、RFI和频段分配竞争等问题，存在着成本低、功耗小、价格便宜、时空分辨率高等诸多优点。
机会信号反射遥感作为一种新型的根区土壤水分遥感手段，其前期的机理仿真分析模型为分析解释卫星观测数据、卫星数据的仿真模拟、卫星数据同化、发展地表参数的定量反演等提供关键的机理工具。然而，目前，现有的仿真分析系统多是针对辐射计或者单站雷达或者GNSS-R遥感的，对于机会反射信号遥感，针对根区土壤水分，尚未发现有相应的仿真分析系统。
发明内容
本发明的目的在于提供一种机会信号反射遥感的全极化单次反射仿真方法及系统，从而通过机会反射信号遥感技术在根区土壤水分发生变化时预测出相应的接收机波形的变化。
为了实现上述目的，本发明提供了一种机会信号反射遥感的全极化单次反射仿真方法，其工作于P波段，包括：
S1：将一待测土壤剖面划分为多个分层，并根据待测土壤剖面的输入参数建立待测土壤剖面的介电常数模型，得到待测土壤剖面的各个分层的介电常数；
S2：利用所述步骤S1中的待测土壤剖面的各个分层的介电常数来获取每个分层的透过率、反射率和每两个相邻分层之间的衰减；
S3：根据各个分层的粗糙程度，分别构建光滑表面反射率矩阵和粗糙表面反射率矩阵；
S4：根据所述步骤S3中的光滑表面反射率矩阵和粗糙表面反射率矩阵，利用极化合成的方法计算各个分层在各个极化状态上的反射特性，该反射特性为各个分层的双站雷达散射截面；
S5：根据所述步骤S2中的每个分层的透过率、反射率和每两个相邻分层之间的衰减以及所述步骤S4中的反射特性来获得单次反射的总接收机波形。
所述全极化单次反射仿真系统的输入参数包括：土壤质地、土壤温度、土壤块密度、地表粗糙度和土壤含水量。
所述待测土壤剖面的介电常数模型是基于混合介质介电常数模型修正得到的，所述介电常数的实部采用线性校正来校正，所述介电常数的虚部ξ"通过修正有效电导率来修正。
所述介电常数的实部的校正公式为：
ξ'＝1.15ξ"-0.68
式中，ξ'和ξ"分别是介电常数的实部和虚部；
有效电导率σeff的修正公式为：
σeff＝0.0467+0.2204ρb-0.4111*S+0.6614*C，
式中，σeff是有效电导率，ρb是土壤块密度，S和C分别是沙土含量和粘土含量。
所述每个分层的透过率、反射率和每两个相邻分层之间的衰减为：
Figure GDA0003551218600000031
Figure GDA0003551218600000032
Figure GDA0003551218600000033
Figure GDA0003551218600000034
Figure GDA0003551218600000035
其中，下标i和j分别是彼此相邻的分层中的第i层和第j层，Tv和Th分别为第i层的垂直和水平极化的透过率，Rv和Rh分别是第i层的垂直和水平极化的反射率；θ为入射角，εi为第i层的介电常数；λ为波长，α为第i层和第j层之间的衰减系数。
在所述步骤S3中，对于分层为光滑表面的情况，则构建光滑表面反射率矩阵，所述光滑表面反射率矩阵是通过假定每层只发生镜像反射，基于菲涅尔公式来构建的；针对分层为粗糙表面的情况，则构建粗糙表面反射率矩阵，所述粗糙表面反射率矩阵是通过建立高级积分方程模型计算粗糙表面的反射率而得到的。
在所述步骤S3中，极化合成的计算公式为：
Figure GDA0003551218600000041
其中，σqp为各个分层的双站雷达散射截面，下标q和p分别为接受和发射的极化状态；ψ和τ分别为椭倾角和椭率角；Yq、Yp为修改的stokes矢量；Q为旋转矩阵；M为光滑表面反射率矩阵或者粗糙表面反射率矩阵。
在所述步骤S5中，所述总接收机波形为：
Figure GDA0003551218600000042
式中，τ为时间延迟，上标dir表示接收机能量中的直射部分，i表示分层中的第i层，N为层数，ρi为第i层的延迟，LC/A为GPS的C/A码片长度，j为复数单位，
Figure GDA0003551218600000043
Ui为第i层的振幅，φi为第i层的相位；
其中，第i层的振幅Ui、相位φi和延迟ρi根据所述每个分层的透过率、反射率、每两个相邻分层之间的衰减以及双站雷达散射截面计算得到。
另一方面，本发明提供一种机会信号反射遥感的全极化单次反射仿真系统，其工作于P波段，包括：
土壤剖面介电常数计算模块，设置为将一待测土壤剖面划分为多个分层，并根据待测土壤剖面的输入参数建立待测土壤剖面的介电常数模型，得到待测土壤剖面的各个分层的介电常数；
层内反射率计算模块、层内透过率计算模块和层间衰减率计算模块，分别设置为所述各个分层的介电常数来获取每个分层的透过率、反射率和每两个相邻分层之间的衰减；
全极化反射率计算模块，根据所述层内反射率计算模块得到的反射率得到全极化的反射特性，其包括：镜像反射模块、非镜像反射模块以及与这两者均相连的极化合成模块，镜像反射模块和非镜像反射模块分别设置为根据分层的粗糙程度构建光滑表面反射率矩阵和粗糙表面反射率矩阵；极化合成模块设置为根据所述光滑表面反射率矩阵和粗糙表面反射率矩阵，利用极化合成的方法计算各个分层在各个极化状态上的反射特性，该反射特性为各个分层的双站雷达散射截面；
单次反射能量计算模块设置为根据所述每个分层的透过率、反射率和每两个相邻分层之间的衰减以及所述各个分层在各个极化状态上的反射特性来获得总接收机波形。
所述全极化单次反射仿真系统的输入参数包括：土壤质地、土壤温度、土壤块密度、地表粗糙度和土壤含水量。
本发明的机会信号反射遥感的全极化单次反射仿真方法通过利用混合介质介电常数模型计算土壤的各个分层的介电常数，对根区土壤水分建立相应的仿真，从而在当根区土壤水分发生变化时，预测出相应的接收机波形的变化，为分析解释卫星观测数据、卫星数据的仿真模拟、卫星数据同化、发展地表参数的定量反演等提供关键的机理工具；此外，本发明的全极化单次反射仿真方法采用极化合成的方法，使其可以应用在机会信号反射遥感中，并且可以计算各种极化的反射率和波形信息，即机会信号卫星发射右旋圆极化信号，接收到圆极化(右旋圆极化和左旋圆极化)和线极化(水平极化和垂直极化)的总接收机波形，由于电磁波的极化对目标的介电常数、物理特性、几何形状和取向等比较敏感，因而极化测量可以大大提高对目标各种信息的获取能力；再者，现有的模型只可以计算每个小分层内的镜像反射，而本发明通过计算每个分层内的入射过来的机会卫星信号的镜像散射和非镜像的漫散射，从而避免了对各个方向的漫散射的忽略，使计算结果更准确。
附图说明
图1为根据本发明的一个实施例的P波段的机会信号反射遥感的全极化单次反射仿真系统的结构示意图。
图2为具有多个分层的土壤剖面的单次反射信号的追踪示意图。
具体实施方式
以下结合具体实施例，对本发明做进一步说明。应理解，以下实施例仅用于说明本发明而非用于限制本发明的范围。
如图1所示，根据本发明的一个实施例的机会信号反射遥感的全极化单次反射仿真方法工作于P波段，其包括以下步骤：
步骤S1：采用一土壤剖面介电常数计算模块1将一待测土壤剖面划分为多个分层，并根据待测土壤剖面的输入参数建立待测土壤剖面的介电常数模型，得到待测土壤剖面的各个分层的介电常数；
由此，模型计算的时候将土壤剖面划分为不同的多个分层，每层内发生单次散射，具有多个分层的土壤剖面的具体示意图如图2所示。所述划分方式为包括等间隔或者任意间隔的划分。
所述全极化单次反射仿真系统的输入参数包括：土壤质地(沙土含量、粘土含量)、土壤温度、土壤块密度、地表粗糙度(相关长度、均方根高度)、土壤含水量。这些参数都可以通过实际测量获得，或者利用现有文献的输入参数，或者使用者根据经验得到的数据。
待测土壤剖面的介电常数模型是基于Dobson等人的混合介质介电常数模型修正得到的。0.3-1.3GHz情况下，需要针对现有的混合介质介电常数模型进行修改，以建立待测土壤剖面的介电常数模型。
所述介电常数的实部通过采用线性校正来校正。
所述介电常数的实部的校正公式为：
ξ'＝1.15ξ"-0.68 (1)
式中，ξ'和ξ"分别是介电常数的实部和虚部。
所述介电常数的虚部ξ"，则通过修正有效电导率来修正。
其中，有效电导率σeff的修正采用以下公式：
σeff＝0.0467+0.2204ρb-0.4111*S+0.6614*C (2)
式中，σeff是有效电导率，ρb是土壤块密度，S和C分别是沙土含量和粘土含量。其中，土壤块密度ρb、沙土含量S和粘土含量C可以是使用者根据相应的实验测得，或者根据经验值给定。
介电常数详细的计算方法参照参考文献【Hallikainen M T.MicrowaveDielectric Behavior of wet soil-Part I:Empirical Models and experimentalobservations[J].IEEE Trans.Geosci.Remote Sensing,1985,23.】。
步骤S2：分别采用层内反射率计算模块2、层内透过率计算模块3和层间衰减率计算模块4，利用所述步骤S1中的待测土壤剖面的各个分层的介电常数来获取每个分层的透过率、反射率和每两个相邻分层之间的衰减。
其中，待测土壤的每个分层的透过率、反射率和每两个相邻分层之间的衰减为：
Figure GDA0003551218600000071
Figure GDA0003551218600000072
Figure GDA0003551218600000073
Figure GDA0003551218600000074
Figure GDA0003551218600000075
其中，下标i和j分别是彼此相邻的分层中的第i层和第j层，Tv和Th分别为第i层的垂直和水平极化的透过率，Rv和Rh分别是第i层的垂直和水平极化的反射率；θ为入射角，εi为第i层的介电常数；λ为波长，α为第i层和第j层之间的衰减系数，其受介电常数的虚部影响。
所述步骤S2求得的待测土壤的每个分层的透过率、反射率和每两个相邻分层之间的衰减可以用来计算下文的步骤S5中的总接收机波形的振幅、延迟。
步骤S3：根据各个分层的粗糙程度，采用镜像反射模块51和非镜像反射模块52分别构建光滑表面反射率矩阵和粗糙表面反射率矩阵。由此，可以为后续步骤S4中的极化合成，也就是使得反射率可以计算各种极化的反射特性做准备。
其中，对于分层为光滑表面的情况，则构建光滑表面反射率矩阵，该分层所对应的光滑表面反射率矩阵是通过假定每层只发生镜像反射，基于菲涅尔公式来构建的。
所述光滑表面反射率矩阵为：
Figure GDA0003551218600000081
其中，下标v和h分别为垂直极化和水平极化，Rv和Rh分别是各个分层的垂直和水平极化的反射率，其可以通过上文中的公式(5)和公式(6)得到；Re是取实部，Im是虚部，*表示共轭。
其中，针对分层为粗糙表面的情况，则构建粗糙表面反射率矩阵，该分层所对应的粗糙表面反射率矩阵是通过建立高级积分方程模型计算粗糙地表的反射率而得到的。
其中，粗糙表面的反射率为：
Figure GDA0003551218600000082
其中，
Figure GDA0003551218600000083
为粗糙表面的反射率的基尔霍夫项，
Figure GDA0003551218600000084
为基尔霍夫项的补偿项，
Figure GDA0003551218600000085
为基尔霍夫项和基尔霍夫项的补偿项这二者的交叉项，下标q和p分别为接受和发射的极化状态，极化状态具体包括圆极化(右旋圆极化和左旋圆极化)和线极化(水平极化和垂直极化)。
基尔霍夫项
Figure GDA0003551218600000086
基尔霍夫项的补偿项
Figure GDA0003551218600000087
和基尔霍夫项和交叉项
Figure GDA0003551218600000088
这三个参数均是采用上文的各个分层的垂直和水平极化的反射率Rv和Rh计算得到的，具体的计算方法在参考文献【Chen,K.S.；Wu,T.D.；Tsang,L.；Li,Q.；Shi,J.；Fung,A.K.Emission of roughsurfaces calculated by the integral equation method with comparison to three-dimensional moment method simulations.IEEE Trans.Geosci.Remote Sens.2003,41,90–101】中给出。
所述粗糙表面的反射率在各个极化状态上的计算结果分别为所述粗糙表面反射率矩阵的元素。
步骤S4：极化合成。采用极化合成模块53，根据所述步骤S3中的光滑表面反射率矩阵和粗糙表面反射率矩阵，利用极化合成的方法计算各个分层在各个极化状态上的反射特性，该反射特性为各个分层的双站雷达散射截面(BRCS)
Figure GDA0003551218600000091
其中，极化合成的计算公式如下：
Figure GDA0003551218600000092
其中，σqp为各个分层的双站雷达散射截面，下标q和p分别为接受和发射的极化状态，极化状态具体包括圆极化(右旋圆极化和左旋圆极化)和线极化(水平极化和垂直极化)；ψ和τ分别为椭倾角和椭率角，不同的极化对应不同的椭倾角和椭率角；Yq、Yp为修改的stokes矢量；Q为旋转矩阵；M为散射矩阵，即为上文所述的光滑表面反射率矩阵或者粗糙表面反射率矩阵，当该分层为光滑表面时，M为光滑表面反射率矩阵，该分层为粗糙表面时，M为粗糙表面反射率矩阵。
由此，将光滑表面反射率矩阵和粗糙表面反射率矩阵作为公式中的M输入，既可以得到光滑表面或者粗糙表面的双站雷达散射截面
Figure GDA0003551218600000093
公式中的
Figure GDA0003551218600000094
就是最终需要的计算结果。
极化合成的具体计算方法可以参考相关文献【Ulaby F T,Elachi C.RadarPolarimetry for Geoscience Applications[J].Geocarto International,1990,5(3):38-38.】。
由此，本发明通过构造4*4的实数矩阵，通过极化合成的方法，使得模型可以计算各种极化的反射率，并通过计算每个分层内的入射过来的机会卫星信号的镜像散射和非镜像的漫散射，从而避免了对各个方向的漫散射的忽略，使计算结果更准确。
步骤S5：根据所述步骤S2中的每个分层的透过率、反射率和每两个相邻分层之间的衰减以及所述步骤S4中的双站雷达散射截面来获得单次反射的总接收机波形。
其中，所述总接收机波形为：
Figure GDA0003551218600000095
式中，τ为时间延迟，上标dir表示接收机能量中的直射部分，i表示分层中的第i层，N为层数，ρi为第i层的延迟，LC/A为GPS的C/A码片长度，j为复数单位，
Figure GDA0003551218600000101
Ui为第i层的振幅，φi为第i层的相位。
其中，第i层的振幅Ui、相位φi和延迟ρi根据所述每个分层的透过率、反射率、每两个相邻分层之间的衰减以及双站雷达散射截面计算得到，具体的表达式参见文献【Boniface,K.,et al.,Comparison of Snow Data Assimilation System with GPSreflectometry snow depth in the Western United States.Hydrological Processes,2015.29(10):p.2425-2437.】以及【Cardellach,E.,et al.,Characterization of dry-snow sub-structure using GNSS reflected signals.Remote Sensing ofEnvironment,2012.124(0):p.122-134.】
所得到的总接收机波形包括各种极化(RR、RL、RV、RH)波形。
由此，利用上述仿真分析方法，可以给出当根区土壤水分发生变化时，预测出相应的接收机波形的变化，进而可以得到针对根区土壤水分的全极化的单次反射仿真模型，从而其前期的机理仿真分析模型为分析解释卫星观测数据、卫星数据的仿真模拟、卫星数据同化、发展地表参数的定量反演等提供关键的机理工具。
如图1所示，本发明提供一种机会信号反射遥感的全极化单次反射仿真系统，其包括土壤剖面介电常数计算模块1、层内反射率计算模块2、层内透过率计算模块3、层间衰减率计算模块4、全极化反射率计算模块5和单次反射能量计算模块6。
其中，土壤剖面介电常数计算模块1设置为将一待测土壤剖面划分为多个分层，并根据待测土壤剖面的输入参数建立待测土壤剖面的介电常数模型，得到待测土壤剖面的各个分层的介电常数；所述全极化单次反射仿真系统的输入参数包括：土壤质地(沙土含量、粘土含量)、土壤温度、土壤块密度、地表粗糙度(相关长度、均方根高度)、土壤含水量。这些参数都可以通过实际测量获得，或者利用现有文献的输入参数，或者使用者根据经验得到的数据。
层内反射率计算模块2、层内透过率计算模块3和层间衰减率计算模块4分别设置为根据土壤剖面介电常数计算模块1得到的各个分层的介电常数来获取每个分层的透过率、反射率和每两个相邻分层之间的衰减。
全极化反射率计算模块5设置为根据所述层内反射率计算模块2得到的反射率得到全极化的反射特性，其包括镜像反射模块51、非镜像反射模块52以及与所述镜像反射模块51和非镜像反射模块52均相连的极化合成模块53，镜像反射模块51和非镜像反射模块52分别设置为根据分层的粗糙程度构建光滑表面反射率矩阵和粗糙表面反射率矩阵，具体地，当分层为光滑表面时，镜像反射模块51设置为构建光滑表面反射率矩阵，当分层为粗糙表面时，非镜像反射模块52设置为构建粗糙表面反射率矩阵；极化合成模块53设置为根据所述光滑表面反射率矩阵和粗糙表面反射率矩阵，利用极化合成的方法计算各个分层在各个极化状态上的反射特性(即双站雷达散射截面BRCS)。
单次反射能量计算模块6设置为根据所述层内反射率计算模块2、层内透过率计算模块3和层间衰减率计算模块4得到的每个分层的透过率、反射率和每两个相邻分层之间的衰减以及所述全极化反射率计算模块5得到的各个分层在各个极化状态上的反射特性来获得总接收机波形；
由此，本发明针对根区土壤水分首先构造了一个土壤剖面，利用混合介质介电常数模型计算这个分层土壤剖面的各个分层的介电常数；然后利用菲涅尔反射系数和高级积分方程模型来分别计算每个分层的镜像散射和非镜像的漫散射；利用极化合成的方法使得模型可以计算各种极化状态上的反射系数；最后将各种极化状态上的反射系数作为模型的一部分放在总接收机波形的公式中进行计算，可以给出当根区土壤水分发生变化时，预测出相应的接收机波形的变化，进而可以得到针对根区土壤水分的全极化的单次反射仿真模型。
以上所述的，仅为本发明的较佳实施例，并非用以限定本发明的范围，本发明的上述实施例还可以做出各种变化。凡是依据本发明申请的权利要求书及说明书内容所作的简单、等效变化与修饰，皆落入本发明专利的权利要求保护范围。本发明未详尽描述的均为常规技术内容。

Claims (10)

1.一种机会信号反射遥感的全极化单次反射仿真方法，其特征在于，其工作于P波段，包括：
步骤S1：将一待测土壤剖面划分为多个分层，并根据待测土壤剖面的输入参数建立待测土壤剖面的介电常数模型，得到待测土壤剖面的各个分层的介电常数；
步骤S2：利用所述步骤S1中的待测土壤剖面的各个分层的介电常数来获取每个分层的透过率、反射率和每两个相邻分层之间的衰减；
步骤S3：根据各个分层的粗糙程度，分别构建光滑表面反射率矩阵和粗糙表面反射率矩阵；
步骤S4：根据所述步骤S3中的光滑表面反射率矩阵和粗糙表面反射率矩阵，利用极化合成的方法计算各个分层在各个极化状态上的反射特性，该反射特性为各个分层的双站雷达散射截面；
步骤S5：根据所述步骤S2中的每个分层的透过率、反射率和每两个相邻分层之间的衰减以及所述步骤S4中的反射特性来获得单次反射的总接收机波形。
2.根据权利要求1所述的机会信号反射遥感的全极化单次反射仿真方法，其特征在于，所述全极化单次反射仿真系统的输入参数包括：土壤质地、土壤温度、土壤块密度、地表粗糙度和土壤含水量。
3.根据权利要求1所述的机会信号反射遥感的全极化单次反射仿真方法，其特征在于，所述待测土壤剖面的介电常数模型是基于混合介质介电常数模型修正得到的，所述介电常数的实部采用线性校正来校正，所述介电常数的虚部ξ"通过修正有效电导率来修正。
4.根据权利要求3所述的机会信号反射遥感的全极化单次反射仿真方法，其特征在于，所述介电常数的实部的校正公式为：
ξ'＝1.15ξ"-0.68
式中，ξ'和ξ"分别是介电常数的实部和虚部；
有效电导率σeff的修正公式为：
σeff＝0.0467+0.2204ρb-0.4111*S+0.6614*C，
式中，σeff是有效电导率，ρb是土壤块密度，S和C分别是沙土含量和粘土含量。
5.根据权利要求1所述的机会信号反射遥感的全极化单次反射仿真方法，其特征在于，所述每个分层的透过率、反射率和每两个相邻分层之间的衰减为：
Figure FDA0003551218590000021
Figure FDA0003551218590000022
Figure FDA0003551218590000023
Figure FDA0003551218590000024
Figure FDA0003551218590000025
其中，下标i和j分别是彼此相邻的分层中的第i层和第j层，Tv和Th分别为第i层的垂直和水平极化的透过率，Rv和Rh分别是第i层的垂直和水平极化的反射率；θ为入射角，εi为第i层的介电常数；λ为波长，α为第i层和第j层之间的衰减系数。
6.根据权利要求1所述的机会信号反射遥感的全极化单次反射仿真方法，其特征在于，在所述步骤S3中，对于分层为光滑表面的情况，则构建光滑表面反射率矩阵，所述光滑表面反射率矩阵是通过假定每层只发生镜像反射，基于菲涅尔公式来构建的；针对分层为粗糙表面的情况，则构建粗糙表面反射率矩阵，所述粗糙表面反射率矩阵是通过建立高级积分方程模型计算粗糙表面的反射率而得到的。
7.根据权利要求1所述的机会信号反射遥感的全极化单次反射仿真方法，其特征在于，在所述步骤S3中，极化合成的计算公式为：
Figure FDA0003551218590000031
其中，σqp为各个分层的双站雷达散射截面，下标q和p分别为接受和发射的极化状态；ψ和τ分别为椭倾角和椭率角；Yq、Yp为修改的stokes矢量；Q为旋转矩阵；M为光滑表面反射率矩阵或者粗糙表面反射率矩阵。
8.根据权利要求1所述的机会信号反射遥感的全极化单次反射仿真方法，其特征在于，在所述步骤S5中，所述总接收机波形为：
Figure FDA0003551218590000032
式中，τ为时间延迟，上标dir表示接收机能量中的直射部分，i表示分层中的第i层，N为层数，ρi为第i层的延迟，LC/A为GPS的C/A码片长度，j为复数单位，
Figure FDA0003551218590000033
Ui为第i层的振幅，φi为第i层的相位；
其中，第i层的振幅Ui、相位φi和延迟ρi根据所述每个分层的透过率、反射率、每两个相邻分层之间的衰减以及双站雷达散射截面计算得到。
9.一种机会信号反射遥感的全极化单次反射仿真系统，其特征在于，其工作于P波段，包括：
土壤剖面介电常数计算模块(1)，设置为将一待测土壤剖面划分为多个分层，并根据待测土壤剖面的输入参数建立待测土壤剖面的介电常数模型，得到待测土壤剖面的各个分层的介电常数；
层内反射率计算模块(2)、层内透过率计算模块(3)和层间衰减率计算模块(4)，分别设置为所述各个分层的介电常数来获取每个分层的透过率、反射率和每两个相邻分层之间的衰减；
全极化反射率计算模块(5)，根据所述层内反射率计算模块(2)得到的反射率得到全极化的反射特性，其包括：镜像反射模块(51)、非镜像反射模块(52)以及与这两者均相连的极化合成模块(53)，镜像反射模块(51)和非镜像反射模块(52)分别设置为根据分层的粗糙程度构建光滑表面反射率矩阵和粗糙表面反射率矩阵；极化合成模块(53)设置为根据所述光滑表面反射率矩阵和粗糙表面反射率矩阵，利用极化合成的方法计算各个分层在各个极化状态上的反射特性，该反射特性为各个分层的双站雷达散射截面；
单次反射能量计算模块(6)设置为根据所述每个分层的透过率、反射率和每两个相邻分层之间的衰减以及所述各个分层在各个极化状态上的反射特性来获得总接收机波形。
10.根据权利要求9所述的机会信号反射遥感的全极化单次反射仿真系统，其特征在于，所述全极化单次反射仿真系统的输入参数包括：土壤质地、土壤温度、土壤块密度、地表粗糙度和土壤含水量。
CN202010095663.1A 2020-02-17 2020-02-17 机会信号反射遥感的全极化单次反射仿真方法及系统 Active CN111239209B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010095663.1A CN111239209B (zh) 2020-02-17 2020-02-17 机会信号反射遥感的全极化单次反射仿真方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010095663.1A CN111239209B (zh) 2020-02-17 2020-02-17 机会信号反射遥感的全极化单次反射仿真方法及系统

Publications (2)

Publication Number Publication Date
CN111239209A CN111239209A (zh) 2020-06-05
CN111239209B true CN111239209B (zh) 2022-07-19

Family

ID=70866285

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010095663.1A Active CN111239209B (zh) 2020-02-17 2020-02-17 机会信号反射遥感的全极化单次反射仿真方法及系统

Country Status (1)

Country Link
CN (1) CN111239209B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4050334A1 (en) * 2021-02-26 2022-08-31 Tata Consultancy Services Limited System and method for root zone soil moisture estimation for vegetation cover using remote sensing
CN114252480B (zh) * 2021-12-21 2024-09-24 中南大学 基于单极化sar数据邻域像素反演裸土区土壤湿度方法
CN116699666B (zh) * 2023-08-03 2023-10-31 北京航空航天大学 星下点观测的星载gnss-r海面风场反演方法及系统

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101738620B (zh) * 2008-11-19 2012-02-15 中国农业科学院农业资源与农业区划研究所 从被动微波遥感数据amsr-e反演地表温度的方法
CN101614818B (zh) * 2009-07-09 2012-01-04 中国科学院遥感应用研究所 一种土壤盐碱化的雷达遥感监测方法
CN103336104B (zh) * 2013-06-08 2015-12-09 北京航空航天大学 Gnss卫星反射信号用于土壤介电常数的精确反演方法
CN103969268B (zh) * 2014-04-04 2016-06-29 中国科学院华南植物园 地表土壤物理参数的被动微波遥感反演方法
CN106569209A (zh) * 2016-10-28 2017-04-19 核工业北京地质研究院 一种利用全极化雷达数据提取土壤含水信息的方法
CN106874621B (zh) * 2017-03-10 2018-02-27 长沙理工大学 一种针叶植被冠层反射率计算方法
CN108761034B (zh) * 2018-05-16 2019-04-30 河海大学 一种多源卫星土壤湿度反演方法
CN109344536A (zh) * 2018-10-30 2019-02-15 中国科学院、水利部成都山地灾害与环境研究所 一种耦合多源数据的多层次被动微波土壤水分验证方法
CN109765247B (zh) * 2018-11-15 2021-06-08 中国农业科学院农业资源与农业区划研究所 一种不同生育阶段麦类作物覆盖区域土壤水分反演方法
CN109977574B (zh) * 2019-04-02 2020-10-20 中国科学院遥感与数字地球研究所 一种基于改进Freeman-Durden极化分解模型的土壤水分反演方法
CN110795895A (zh) * 2020-01-06 2020-02-14 南京邮电大学 利用地表反射信号和随机森林回归算法预测土壤水分的方法

Also Published As

Publication number Publication date
CN111239209A (zh) 2020-06-05

Similar Documents

Publication Publication Date Title
CN111239209B (zh) 机会信号反射遥感的全极化单次反射仿真方法及系统
Leinss et al. Anisotropy of seasonal snow measured by polarimetric phase differences in radar time series
CN108303044B (zh) 一种叶面积指数获取方法及系统
Moon et al. RADARSAT-2 and coastal applications: Surface wind, waterline, and intertidal flat roughness
CN105606631B (zh) 一种盐度卫星双波段亮温数据联合反演土壤水分的方法
Bringi et al. Rainfall estimation from C-band polarimetric radar in Okinawa, Japan: Comparisons with 2D-video disdrometer and 400 MHz wind profiler
CN103336104B (zh) Gnss卫星反射信号用于土壤介电常数的精确反演方法
Al-Ashwal et al. Analysis of bistatic sea clutter-Part I: Average reflectivity
CN101900692A (zh) 大面积土壤湿度测量方法
CN103344653B (zh) 一种基于双gps接收机的土壤湿度实时测量系统及方法
CN111190151B (zh) 扫描模式下多模态小卫星sar的系统参数设计及发射功率优化方法
Wang et al. Assessment of different vegetation parameters for parameterizing the coupled water cloud model and advanced integral equation model for soil moisture retrieval using time series Sentinel-1A data
Reid et al. Mesospheric radar wind comparisons at high and middle southern latitudes
Li et al. A microphysics-based simulator for advanced airborne weather radar development
Srinivasulu et al. Simplified active array L-band radar for atmospheric wind profiling: initial results
Li et al. Extraction of wind direction spreading factor from broad-beam high-frequency surface wave radar data
Asai et al. Validation of X-band multiparameter phased-array weather radar by comparing data from Doppler weather radar with a parabolic dish antenna
Garrison et al. Remote sensing of root-zone soil moisture using I-and P-band signals of opportunity: Instrument validation studies
CN113189611B (zh) 基于太赫兹后向散射信号的土壤温度和湿度检测方法
CN113325418A (zh) 基于合成孔径雷达的降雨反演方法
Green et al. An inversion method for extraction of wind speed from high-frequency ground-wave radar oceanic backscatter
Wu et al. A full-polarization GNSS-R Delay-Doppler-Map (DDM) simulator for bare soil freeze/thaw process detection
Wu et al. Polarization GNSS-reflectometry: potential and possibility
Privette III et al. Utilizing space-based GPS technology to determine hydrological properties of soils
Barnum et al. High-resolution mapping of oceanic wind fields with skywave radar

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant