CN110515049B - 一种旋转采样综合孔径辐射计相关偏置校正方法 - Google Patents
一种旋转采样综合孔径辐射计相关偏置校正方法 Download PDFInfo
- Publication number
- CN110515049B CN110515049B CN201910725392.0A CN201910725392A CN110515049B CN 110515049 B CN110515049 B CN 110515049B CN 201910725392 A CN201910725392 A CN 201910725392A CN 110515049 B CN110515049 B CN 110515049B
- Authority
- CN
- China
- Prior art keywords
- visibility function
- value
- standard deviation
- offset
- bias
- 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
Links
Images
Classifications
-
- 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/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radiation Pyrometers (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明公开了一种旋转采样综合孔径辐射计相关偏置校正方法,所述方法包括:利用待校正的综合孔径辐射计对一个亮温中心对称的观测目标进行旋转采样,获得各相关基线的可见度函数测量值;任取一个相关基线,对其可见度函数测量值迭代执行校正步骤:计算可见度函数幅度测量值的标准差并确定复相关偏置搜索范围,在该范围内获得使可见度函数幅度标准差最小的复相关偏置最优估计值;对可见度函数测量值进行校正;直至校正后的可见度函数测量值幅度标准差小于预定阈值。本发明的校正方法不依赖额外的内部定标装置,可以有效降低旋转采样综合孔径辐射计的系统复杂度,尤其适用于单元数较多的大规模旋转采样综合孔径辐射计系统。
Description
技术领域
本发明涉及旋转采样综合孔径辐射计领域,尤其涉及一种旋转采样综合孔径辐射计相关偏置校正方法。
背景技术
综合孔径辐射计是被动微波遥感领域中实现微波辐射测量的主要技术手段之一。综合孔径辐射计通过对观测场景亮温分布的空间频率域进行干涉测量,获得可见度函数。然后利用观测场景亮温分布与可见度函数之间的傅里叶变换关系,重建观测场景的亮温图像。相较于真实孔径辐射计,综合孔径技术将大口径实孔径天线稀疏为小口径天线阵列,在实现同等空间分辨率的情况下,可有效降低实际天线物理尺寸,避免大口径实孔径天线制造加工以及空间运载、折叠展开等技术难题。
在实际的综合孔径辐射计系统中,由于不可避免地存在数字相关器通道间互耦和公共本振信号热噪声泄露等非理想性因素,可见度函数测量结果会产生相关偏置误差。对相关偏置误差进行校正是综合孔径辐射计系统误差校正的一个重要环节。目前综合孔径辐射计相关偏置校正方法是采用非相关噪声注入法直接测量各相关基线的偏置误差(参考文献1:I.Corbella et al.,"MIRAS end-to-end calibration:application to SMOSL1processor,"in IEEE Transactions on Geoscience and Remote Sensing,2005,vol.43,no.5,pp.1126-1134.)。
该方法的缺点在于内部定标系统结构复杂,需要内置额外的非相关噪声源,并且配置复杂的噪声功分网络。对于单元数较多的高分辨率综合孔径辐射计系统,该方法难以获得实际应用。
发明内容
本发明的目的在于,克服现有综合孔径辐射计相关偏置校正方法严重依赖内部定标系统的缺点,直接利用观测数据对可见度函数的相关偏置进行校正,无需额外配置内部定标系统,有效降低接收机硬件系统复杂度。
为了实现上述目的,本发明提出了一种旋转采样综合孔径辐射计相关偏置校正方法,所述方法包括:
利用待校正的综合孔径辐射计对一个亮温中心对称的观测目标进行旋转采样,获得各相关基线的可见度函数测量值;
任取一个相关基线,对其可见度函数测量值迭代执行校正步骤:计算可见度函数幅度测量值的标准差并确定复相关偏置搜索范围,在该范围内获得使可见度函数幅度标准差最小的复相关偏置最优估计值;对可见度函数测量值进行校正;直至校正后的可见度函数测量值幅度标准差小于预定阈值;
遍历所有的相关基线,对其可见度函数测量值迭代执行校正步骤,由此校正所有相关基线可见度函数测量值的相关偏置误差。
作为上述方法的一种改进,所述任取一个相关基线,对其可见度函数测量值迭代执行校正步骤,具体包括:
步骤1)对于单元天线i和j组成的相关基线,计算其可见度函数测量值在360°采样范围内的幅度标准差;在正负幅度标准差范围内,以选定的步长遍历该范围,分别搜索该相关基线的可见度函数实部和虚部相关偏置;获得使可见度函数测量值幅度标准差最小的复相关偏置最优估计值;
步骤2)使用复相关偏置最优估计值校正该相关基线的可见度函数测量值;
步骤3)判断校正后的可见度函数幅度标准差的变化量是否小于预定阈值,若判断结果是肯定的,则该相关基线校正结束;否则,将可见度函数测量值更新为校正后的可见度函数测量值,转入步骤1)。
作为上述方法的一种改进,所述步骤1)具体包括:
其中,abs(·)表示取模运算,mean(·)表示取平均值运算,N为该基线可见度函数测量值在360°范围内的采样点数;k为迭代次数,初始值为1;
步骤103)分别采用M+1个可见度函数实部偏置估计值和可见度函数虚部偏置估计值对可见度函数测量值进行校正:
步骤104)分别计算实部偏置校正后和虚部偏置校正后的可见度函数幅度标准差,由此得到使实部偏置校正后的可见度函数幅度标准差最小的实部偏置最优估计值使虚部偏置校正后的可见度函数幅度标准差最小的虚部偏置最优估计值
所述复相关偏置最优估计值为:
其中,1i为虚数符号。
作为上述方法的一种改进,所述步骤2)具体为:
作为上述方法的一种改进,所述步骤3)具体为:
本发明的优点在于:
1、本发明提供的相关偏置校正方法,克服了现有综合孔径辐射计相关偏置校正方法对内部定标系统的依赖,根据亮温中心对称观测目标可见度函数特性,直接利用观测数据对可见度函数的相关偏置进行校正;
2、本发明提供的相关偏置校正方法,外部定标源易于获取,无论是太阳目标还是噪声点源,均适用于地面测试环境和星载观测环境。
附图说明
图1为本发明提出的旋转采样综合孔径辐射计相关偏置校正方法流程图;
图2为实例中亮温中心对称观测目标(太阳)的亮温分布;
图3为实例中的各相关基线可见度函数相关偏置误差分布;
图4为实例中未校正相关偏置误差的可见度函数亮温重建结果;
图5为实例中利用本发明方法进行相关偏置误差校正后的偏置残差分布;
图6为实例中利用本发明方法进行相关偏置误差校正后的可见度函数亮温重建结果。
具体实施方式
下面结合附图和实例对本发明进行详细说明。
如图1所示,本发明提供了一种旋转采样综合孔径辐射计相关偏置校正方法,所述方法包括:
步骤101):利用待校正的综合孔径辐射计对某一亮温中心对称的观测目标(如太阳或噪声点源)进行旋转采样,获得各相关基线可见度函数测量值;
步骤201):计算该基线待校正可见度函数幅度在360°采样范围内的标准差,其计算公式为:
其中,abs(·)表示取模运算,mean(·)表示取平均值运算,N为可见度函数在360°范围内的采样点数,下标ij代表组成相关基线的单元天线编号,一对ij可唯一确定一组相关基线,上标k代表迭代次数,初始值为1。
步骤203):利用式(1)分别计算实部偏置校正后和虚部偏置校正后的可见度函数幅度标准差;
步骤204):分别选取使得偏置校正后的可见度函数幅度标准差在本次迭代中取最小值的实部偏置估计值和虚部偏置估计值,作为本次迭代中的最优复相关偏置估计值;
步骤301):利用步骤204)中得到的最优复相关偏置估计值,对可见度函数进行复偏置误差校正;
步骤401):计算本次迭代中复偏置误差校正后的可见度函数幅度标准差与上次迭代的差异,其计算公式为:
下面进一步通过具体的仿真实例验证说明本发明的方法。
本实例的观测场景为典型的太阳观测场景,符合本发明中要求的亮温中心对称观测目标的应用前提。观测场景亮温分布如图2所示,太阳圆心位于空间域极坐标角度240°处,与坐标原点的距离为0.05,太阳半径为0.0049,太阳亮温为6000K。
可见度函数相关基线长度设置采用等间距分布,最短相关基线长度为3λ,最大相关基线长度为198λ。一个采样周期内每个相关基线在空间频率域复平面上完整覆盖一个采样圆环,以极坐标相角0.5°为采样间隔,每个相关基线共有720个采样点。利用综合孔径辐射计观测亮温与可见度函数的傅里叶变换关系计算各相关基线理想可见度函数,计算公式为:
其中,VFideal(u,v)为连续的理想可见度函数,TB为观测亮温分布,ξ和η为球坐标系方向余弦,u和v为可见度函数采样点在空间频率域中的位置坐标。
对各相关基线理想可见度函数进行离散化,计算公式为:
uij(n)=Re[BLij·exp(1i·θs(n))],n=1,2,…,N
vij(n)=Im[BLij·exp(1i·θs(n))],n=1,2,…,N
其中,是理想可见度函数的离散采样,下标ij为组成该相关基线的天线单元编号,N为该相关基线可见度函数的采样点数,BLij为该相关基线波长归一化基线长度,θs为相关基线在空间频率域复平面上的采样角度,1i代表虚数符号。
设置各相关基线实部和虚部偏置误差分别与该相关基线可见度实部幅度和虚部幅度在一个采样周期内的平均值呈正比,各相关基线偏置误差计算公式为:
其中,和分别为实部和虚部的偏置误差,pbias为误差偏置比例,设置为20%,signreal和signimag为偏置误差的符号,随机选取-1或1中的任意一个值;本实例中各相关基线偏置误差分布如图3所示。
通过对理想可见度函数施加相关偏置误差和随机相位误差,生成误差可见度函数,其误差模型为:
直接对误差可见度函数进行相位定标,其亮温重建结果如图4所示。由于未对相关偏置误差进行校正,可以观察到重建结果中存在以坐标原点为圆心的环状误差。
利用本发明的方法对误差可见度函数进行相关偏置校正,遍历误差点数设为10,迭代停止预定阈值设为1e-5。计算各基线可见度函数残留偏置误差分布如图5所示,可见本发明的方法可以将各相关基线残留偏置误差控制在1e-5内,残留误差占原误差的比例小于3%。对偏置误差校正后的可见度函数进行相位定标,其亮温重建结果如图6所示。对比图4和图6可以看到,重建结果中的环状误差得到了明显的抑制。需要明确的是,偏置误差校正后的重建结果与原始观测亮温分布具有一定差异,这是由于可见度函数在空间频率域有限采样产生的频域截断误差,不涉及本实例中的相关偏置误差造成的影响。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.一种旋转采样综合孔径辐射计相关偏置校正方法,所述方法包括:
利用待校正的综合孔径辐射计对一个亮温中心对称的观测目标进行旋转采样,获得各相关基线的可见度函数测量值;
任取一个相关基线,对其可见度函数测量值迭代执行校正步骤:计算可见度函数幅度测量值的标准差并确定复相关偏置搜索范围,在该范围内获得使可见度函数幅度标准差最小的复相关偏置最优估计值;对可见度函数测量值进行校正;直至校正后的可见度函数测量值幅度标准差小于预定阈值;
遍历所有的相关基线,对其可见度函数测量值迭代执行校正步骤,由此校正所有相关基线可见度函数测量值的相关偏置误差。
2.根据权利要求1所述的综合孔径辐射计相关偏置校正方法,其特征在于,所述任取一个相关基线,对其可见度函数测量值迭代执行校正步骤,具体包括:
步骤1)对于单元天线i和j组成的相关基线,计算其可见度函数测量值在360°采样范围内的幅度标准差;在正负幅度标准差范围内,以选定的步长遍历该范围,分别搜索该相关基线的可见度函数实部和虚部相关偏置;获得使可见度函数测量值幅度标准差最小的复相关偏置最优估计值;
步骤2)使用复相关偏置最优估计值校正该相关基线的可见度函数测量值;
步骤3)判断校正后的可见度函数幅度标准差的变化量是否小于预定阈值,若判断结果是肯定的,则该相关基线校正结束;否则,将可见度函数测量值更新为校正后的可见度函数测量值,转入步骤1)。
3.根据权利要求2所述的旋转采样综合孔径辐射计相关偏置校正方法,其特征在于,所述步骤1)具体包括:
其中,abs(·)表示取模运算,mean(·)表示取平均值运算,N为该基线可见度函数测量值在360°范围内的采样点数;k为迭代次数,初始值为1;
步骤103)分别采用M+1个可见度函数实部偏置估计值和可见度函数虚部偏置估计值对可见度函数测量值进行校正:
步骤104)分别计算实部偏置校正后和虚部偏置校正后的可见度函数幅度标准差,由此得到使实部偏置校正后的可见度函数幅度标准差最小的实部偏置最优估计值使虚部偏置校正后的可见度函数幅度标准差最小的虚部偏置最优估计值
所述复相关偏置最优估计值为:
其中,1i为虚数符号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910725392.0A CN110515049B (zh) | 2019-08-07 | 2019-08-07 | 一种旋转采样综合孔径辐射计相关偏置校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910725392.0A CN110515049B (zh) | 2019-08-07 | 2019-08-07 | 一种旋转采样综合孔径辐射计相关偏置校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110515049A CN110515049A (zh) | 2019-11-29 |
CN110515049B true CN110515049B (zh) | 2021-03-09 |
Family
ID=68624651
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910725392.0A Active CN110515049B (zh) | 2019-08-07 | 2019-08-07 | 一种旋转采样综合孔径辐射计相关偏置校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110515049B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113702804B (zh) * | 2021-07-26 | 2024-02-09 | 西安空间无线电技术研究所 | 一种多通道数字相关器中相关误差校正方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3965501B2 (ja) * | 2004-12-02 | 2007-08-29 | 三菱電機株式会社 | マイクロ波放射計 |
CN102353944B (zh) * | 2011-05-31 | 2013-01-23 | 西安空间无线电技术研究所 | 一种直接相关型全极化微波辐射计接收机定标方法 |
CN105548973B (zh) * | 2015-12-04 | 2017-12-12 | 中国科学院国家空间科学中心 | 一种基于旋转分时采样的综合孔径辐射计相位自定标方法 |
CN106450703A (zh) * | 2016-11-24 | 2017-02-22 | 中国科学院国家空间科学中心 | 一种基于循环子阵的干涉式微波辐射计圆环天线阵列 |
CN109541559B (zh) * | 2018-11-12 | 2020-06-19 | 北京航空航天大学 | 一种基于psk技术的宽带模拟复相关器实时校准方法 |
CN109725317B (zh) * | 2018-12-18 | 2021-06-01 | 中国人民解放军国防科技大学 | 一种基于一维综合孔径微波辐射计的海面亮温成像仿真方法 |
-
2019
- 2019-08-07 CN CN201910725392.0A patent/CN110515049B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110515049A (zh) | 2019-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101598798B (zh) | 一种对高光谱干涉数据进行光谱重建的系统和方法 | |
Pearson et al. | Image formation by self-calibration in radio astronomy | |
Intema et al. | Ionospheric calibration of low frequency radio interferometric observations using the peeling scheme-I. Method description and first results | |
Shirzaei | A wavelet-based multitemporal DInSAR algorithm for monitoring ground surface motion | |
CN109974665B (zh) | 一种针对缺少短波红外数据的气溶胶遥感反演方法及系统 | |
Stober et al. | Retrieving horizontally resolved wind fields using multi-static meteor radar observations | |
US20120188119A9 (en) | Identification and Analysis of Persistent Scatterers In Series of SAR Images | |
Wang et al. | A new likelihood function for consistent phase series estimation in distributed scatterer interferometry | |
Soldo et al. | Mitigation of RFIs for SMOS: A distributed approach | |
CN109471095A (zh) | 基于快速迭代插值的fmcw雷达距离估计算法 | |
CN110109100B (zh) | 一种基于质量图加权的多基线最小二乘相位解缠方法 | |
CN115060208A (zh) | 基于多源卫星融合的输变电线路地质灾害监测方法及系统 | |
CN110515049B (zh) | 一种旋转采样综合孔径辐射计相关偏置校正方法 | |
Muth et al. | Robust separation of background and target signals in radar cross section measurements | |
CN111505637B (zh) | 一种基于二单元扫描干涉仪的自定标近场成像方法及其系统 | |
Xu et al. | Physically unrolling network under contraction integral equation for limited-aperture inverse scattering problem | |
CN109061645B (zh) | 综合孔径辐射计虚拟相关稀疏成像方法及系统 | |
Yin et al. | Super-resolution compressive spherical beamforming based on off-grid sparse Bayesian inference | |
CN110515048B (zh) | 一种旋转采样综合孔径辐射计相位定标方法 | |
Dokara et al. | Metrewave Galactic Plane with the uGMRT (MeGaPluG) Survey: Lessons from the pilot study | |
Esposito et al. | A simple solution for the phase offset estimation of airborne SAR interferograms without using corner reflectors | |
Jin et al. | An on-board external calibration method for aperture synthesis radiometer by rotation | |
CN107607946B (zh) | 三维均匀采样综合孔径辐射计亮温反演方法 | |
CN111121980A (zh) | 一种基于延时相关的二维综合孔径辐射计成像方法及系统 | |
Das et al. | Statistical isotropy violation in WMAP CMB maps due to non-circular beams |
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 |