CN107607946A - 三维均匀采样综合孔径辐射计亮温反演方法 - Google Patents

三维均匀采样综合孔径辐射计亮温反演方法 Download PDF

Info

Publication number
CN107607946A
CN107607946A CN201710607006.9A CN201710607006A CN107607946A CN 107607946 A CN107607946 A CN 107607946A CN 201710607006 A CN201710607006 A CN 201710607006A CN 107607946 A CN107607946 A CN 107607946A
Authority
CN
China
Prior art keywords
mrow
dimensional
uniform sampling
msub
synthetic aperture
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.)
Granted
Application number
CN201710607006.9A
Other languages
English (en)
Other versions
CN107607946B (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.)
Hubei University of Technology
Original Assignee
Hubei University of Technology
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 Hubei University of Technology filed Critical Hubei University of Technology
Priority to CN201710607006.9A priority Critical patent/CN107607946B/zh
Publication of CN107607946A publication Critical patent/CN107607946A/zh
Application granted granted Critical
Publication of CN107607946B publication Critical patent/CN107607946B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Radiation Pyrometers (AREA)

Abstract

本发明公开了一种三维均匀采样综合孔径辐射计亮温反演方法,其特征在于,所述方法包括如下步骤:S1:对三维均匀采样综合孔径辐射计输出的可见度进行三维傅里叶反变换,得到三维图像体;S2:将反演出的三维图像体映射为二维图像。本发明的三维均匀采样综合孔径辐射计亮温反演方法,可以解决三维均匀采样综合孔径辐射计反演精度差,视场小的难题,有效提升三维均匀采样综合孔径辐射计的反演精度和观测视场,从而推动三维均匀采样综合孔径辐射计的实际应用。

Description

三维均匀采样综合孔径辐射计亮温反演方法
技术领域
本发明涉及微波遥感、目标探测以及深空探测技术领域,具体涉及一种三维均匀采样综合孔径辐射计亮温反演方法。
背景技术
综合孔径辐射计利用多个离散的小天线合成等效的大天线孔径,采用稀疏阵列排布,减少天线的质量和体积,可提高被动微波遥感和探测设备的空间分辨率。为保证反演图像的质量,目前综合孔径辐射计的天线阵列只能排列成一维线性阵列和二维平面阵列。而在目标探测和深空探测等应用中,将天线阵列排列三维阵列是不可避免的。
目前尚未见国内外专门针对三维均匀采样综合孔径辐射计的反演方法的报道和研究。反演精度较差、视场较小是三维均匀采样综合孔径辐射计实际应用的难题。
发明内容
针对三维均匀采样综合孔径辐射计反演精度较差和视场较小的缺陷以及改进需求,本发明提供了一种三维均匀采样综合孔径辐射计亮温反演方法,其目的在于在较大的视场中精确地反演三维均匀采样综合孔径辐射计的场景亮温分布图像。
为了达到上述目的,本发明所采用的技术方案是:一种三维均匀采样综合孔径辐射计亮温反演方法,其特征在于,所述方法包括如下步骤:
S1:对三维均匀采样综合孔径辐射计输出的可见度进行三维傅里叶反变换,得到三维图像体;
S2:将反演出的三维图像体映射为二维图像。
进一步地,步骤S1中,具体实现过程可表示为:
三维天线阵列中,根据公式un=(xi-xj)/λ、vn=(yi-yj)/λ和wn=(zi-zj)/λ可计算出三维天线阵列综合孔径辐射计的第n个基线的三维空间频率(un,vn,zn);(xi,yi,zi)表示第i个单元天线的空间位置坐标,(xj,yj,zj)表示第j个单元天线的空间位置坐标,λ为辐射计接收信号的波长;
假设三维天线阵列中有M个天线单元,则对应的基线数目为个,由于三维均匀采样的天线阵列会产生冗余基线,所以,剔除冗余基线后n=1,2,3,…,N
将三维均匀采样综合孔径辐射计输出的可见度采样值进行三维傅里叶反变换后得到的图像体可表示为:
式(1)中(ξ,η,κ)表示方向余弦,方向余弦的取值范围大于-1小于1,即ξ,η,κ∈[-1,1];
式(1)中Δs表示最小采样体积,可表示为:
Δs=Δu×Δv×Δw (2)
其中,Δu,Δv,Δw分别表示天线阵列在三维空间频率域中均匀采样的最小间隔;
式(1)中V(un,vn,wn)表示三维均匀采样综合孔径辐射计输出的可见度采样值;
式(1)中表示反演出的三维图像体。
进一步地,步骤S2中,具体实现过程可表示为:
反演出的三维图像体的数学公式可表示为:
其中,δ(·)表示单位冲激函数,令并代入式(3),即可将三维图像体映射为二维图像Tm(ξ,η),其中Tm(ξ,η)即为需要反演的二维图像。
与现有技术相比,本发明的有益效果是:
本发明提出一种三维均匀采样综合孔径辐射计亮温反演方法,解决三维均匀采样综合孔径辐射计反演精度差,视场小的难题,有效提升三维均匀采样综合孔径辐射计的反演精度和观测视场,从而推动三维均匀采样综合孔径辐射计的实际应用。
在现有技术路线下,三维的综合孔径天线阵列通过误差校正的方法,将三维可见度校正为二维可见度,然后再通过较为成熟的二维综合孔径辐射计反演算法得到场景亮温分布。但是这种方法只适用于天线阵列尺寸较小的阵列,并且视场范围较小。本发明提出的方法与现有技术路线完全不同,不需要将三维可见度校正为二维可见度。通过三维可见度能直接反演得到场景亮温分布,能适用于天线阵列尺寸较大的阵列,同时视场范围可达全视场。
附图说明
图1为三维均匀采样综合孔径辐射计的结构示意图。
图2为本发明实施例提供的一种三维空间频率域采样模式。
图3为本发明实施例提供的一种模拟观测场景。
图4为本发明实施例提供的一种三维天线阵列综合孔径辐射计反演出的三维图像体。
图5为本发明实施例提供的从三维图像体映射出的二维图像。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明,下面结合实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
图1示出了三维天线阵列综合孔径辐射计的结构,为了便于说明,仅示出了与本发明实施例相关的部分,详述如下:
三维天线阵列综合孔径辐射计包括依次连接的天线阵列1、接收通道阵列2、相关器3;其中天线阵列1为三维均匀采样阵列,即阵列在三维空间频率域中能够均匀采样。接收通道阵列2包括多个接收通道,一个接收通道对应一个天线单元,接收通道将天线单元接收到的信号进行下变频、滤波和放大。相关器3将接收通道输出的信号两两进行复相关。复相关后的输出为可见度函数。
本发明实施例提供的上述三维均匀采样综合孔径辐射计亮温反演方法具体包括下述步骤:
(1)对三维均匀采样综合孔径辐射计输出的可见度采样进行三维傅里叶反变换。具体计算公式为:
式中,un=(xi-xj)/λ、vn=(yi-yj)/λ和wn=(zi-zj)/λ;其中,(xi,yi,zi)表示第i个单元天线的空间位置坐标,(xj,yj,zj)表示第j个单元天线的空间位置坐标,λ为辐射计接收信号的波长。假设三维天线阵列中有M个天线单元,则对应的基线数目为个。由于三维均匀采样的天线阵列会产生冗余基线。所以,剔除冗余基线后n=1,2,3,…,N式中(ξ,η,κ)表示方向余弦,方向余弦的取值范围大于-1小于1,即ξ,η,κ∈[-1,1]。式中Δs表示最小采样体积,可表示为:
Δs=Δu×Δv×Δw
其中,Δu,Δv,Δw分别表示天线阵列在三维空间频率域中均匀采样的最小间隔。式中V(un,vn,wn)表示三维均匀采样综合孔径辐射计输出的可见度采样值。式中表示反演出的三维图像体。
(2)将反演出的三维图像体映射为二维图像。
具体数学过程可表示为:反演出的三维图像体的数学公式可表示为:
式中,δ(·)表示单位冲激函数。令并代入上式,即可将三维图像体映射为二维图像Tm(ξ,η)。其中Tm(ξ,η)即为需要反演的二维图像。
实施例:三维均匀采样综合孔径辐射计系统;该实施例中,天线阵在三维空间频率域中均匀采样。如图2所示。图中的每个小×一个空间频率域采样点。
图像反演方法具体步骤如下:
(1)模拟观察场景如图3所示。并将图3作为一个三维均匀采样综合孔径辐射计的观测场景,此辐射计输出的可见度为V(un,vn,wn)。将其输出的可见度进行三维傅里叶反变换。具体计算公式为:
式中,un=(xi-xj)/λ、vn=(yi-yj)/λ和wn=(zi-zj)/λ;其中,(xi,yi,zi)表示第i个单元天线的空间位置坐标,(xj,yj,zj)表示第j个单元天线的空间位置坐标,λ为辐射计接收信号的波长。在此实例中基线数目为2205,即N=2205。在此实例中,方向余弦的取值范围大于-0.5小于0.5,即ξ,η,κ∈[-0.5,0.5]。式中Δs=5。反演得到的三维图像如图4所示。
(2)令即取三维图像体中单位球面上的反演值。实现从三维图像体到二维图像的映射。如图5所示。
应当理解的是,本说明书未详细阐述的部分均属于现有技术。
应当理解的是,上述针对较佳实施例的描述较为详细,并不能因此而认为是对本发明专利保护范围的限制,本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的范围情况下,还可以做出替换或变形,均落入本发明的保护范围之内,本发明的请求保护范围应以所附权利要求为准。

Claims (3)

1.一种三维均匀采样综合孔径辐射计亮温反演方法,其特征在于,所述方法包括如下步骤:
S1:对三维均匀采样综合孔径辐射计输出的可见度进行三维傅里叶反变换,得到三维图像体;
S2:将反演出的三维图像体映射为二维图像。
2.如权利要求1所述的三维均匀采样综合孔径辐射计亮温反演方法,其特征在于,步骤S1中,具体实现过程可表示为:
三维天线阵列中,根据公式un=(xi-xj)/λ、vn=(yi-yj)/λ和wn=(zi-zj)/λ可计算出三维天线阵列综合孔径辐射计的第n个基线的三维空间频率(un,vn,zn);(xi,yi,zi)表示第i个单元天线的空间位置坐标,(xj,yj,zj)表示第j个单元天线的空间位置坐标,λ为辐射计接收信号的波长;
假设三维天线阵列中有M个天线单元,则对应的基线数目为个,由于三维均匀采样的天线阵列会产生冗余基线,所以,剔除冗余基线后n=1,2,3,…,N
将三维均匀采样综合孔径辐射计输出的可见度采样值进行三维傅里叶反变换后得到的图像体可表示为:
<mrow> <mover> <mi>T</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <mi>&amp;xi;</mi> <mo>,</mo> <mi>&amp;eta;</mi> <mo>,</mo> <mi>&amp;kappa;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>&amp;Delta;</mi> <mi>s</mi> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mi>V</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>n</mi> </msub> <mo>,</mo> <msub> <mi>v</mi> <mi>n</mi> </msub> <mo>,</mo> <msub> <mi>w</mi> <mi>n</mi> </msub> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>n</mi> </msub> <mi>&amp;xi;</mi> <mo>+</mo> <msub> <mi>v</mi> <mi>n</mi> </msub> <mi>&amp;eta;</mi> <mo>+</mo> <msub> <mi>w</mi> <mi>n</mi> </msub> <mi>&amp;kappa;</mi> <mo>)</mo> </mrow> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
式(1)中(ξ,η,κ)表示方向余弦,方向余弦的取值范围大于-1小于1,即ξ,η,κ∈[-1,1];
式(1)中Δs表示最小采样体积,可表示为:
Δs=Δu×Δv×Δw (2)
其中,Δu,Δv,Δw分别表示天线阵列在三维空间频率域中均匀采样的最小间隔;
式(1)中V(un,vn,wn)表示三维均匀采样综合孔径辐射计输出的可见度采样值;
式(1)中表示反演出的三维图像体。
3.如权利要求2所述的三维均匀采样综合孔径辐射计亮温反演方法,其特征在于,步骤S2中,具体实现过程可表示为:
反演出的三维图像体的数学公式可表示为:
<mrow> <mover> <mi>T</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <mi>&amp;xi;</mi> <mo>,</mo> <mi>&amp;eta;</mi> <mo>,</mo> <mi>&amp;kappa;</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>T</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;xi;</mi> <mo>,</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <mi>&amp;delta;</mi> <mrow> <mo>(</mo> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>&amp;xi;</mi> <mn>2</mn> </msup> <mo>-</mo> <msup> <mi>&amp;eta;</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <mi>&amp;kappa;</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
其中,δ(·)表示单位冲激函数,令并代入式(3),即可将三维图像体映射为二维图像Tm(ξ,η),其中Tm(ξ,η)即为需要反演的二维图像。
CN201710607006.9A 2017-07-24 2017-07-24 三维均匀采样综合孔径辐射计亮温反演方法 Expired - Fee Related CN107607946B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710607006.9A CN107607946B (zh) 2017-07-24 2017-07-24 三维均匀采样综合孔径辐射计亮温反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710607006.9A CN107607946B (zh) 2017-07-24 2017-07-24 三维均匀采样综合孔径辐射计亮温反演方法

Publications (2)

Publication Number Publication Date
CN107607946A true CN107607946A (zh) 2018-01-19
CN107607946B CN107607946B (zh) 2020-03-17

Family

ID=61060042

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710607006.9A Expired - Fee Related CN107607946B (zh) 2017-07-24 2017-07-24 三维均匀采样综合孔径辐射计亮温反演方法

Country Status (1)

Country Link
CN (1) CN107607946B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111538001A (zh) * 2020-03-31 2020-08-14 西南电子技术研究所(中国电子科技集团公司第十研究所) 三维天线阵综合孔径辐射计快速亮温反演方法
CN114812849A (zh) * 2022-03-14 2022-07-29 湖北工业大学 一种基于微波辐射测量的电缆接头内部温度测量方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008049948A2 (es) * 2006-10-25 2008-05-02 Universitat Politécnica de Catalunya Sistema híbrido receptor de señales gnss -reflejadas/radiómetro diferencial de pseudo-correlación para la observación pasiva del océano
CN105242304A (zh) * 2014-07-04 2016-01-13 中国石油化工股份有限公司 基于频率空间域波场延拓的地震资料缺口补偿方法
CN106842203A (zh) * 2017-03-07 2017-06-13 湖北工业大学 一种三维天线阵列综合孔径辐射计的图像反演方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008049948A2 (es) * 2006-10-25 2008-05-02 Universitat Politécnica de Catalunya Sistema híbrido receptor de señales gnss -reflejadas/radiómetro diferencial de pseudo-correlación para la observación pasiva del océano
CN105242304A (zh) * 2014-07-04 2016-01-13 中国石油化工股份有限公司 基于频率空间域波场延拓的地震资料缺口补偿方法
CN106842203A (zh) * 2017-03-07 2017-06-13 湖北工业大学 一种三维天线阵列综合孔径辐射计的图像反演方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YUFANG LI等: "Relationship Between Mirrored Aperture Synthesis Radiometers and Aperture Synthesis Radiometers", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 *
李育芳等: "旋转镜像综合孔径空间频率均匀采样方法", 《微波学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111538001A (zh) * 2020-03-31 2020-08-14 西南电子技术研究所(中国电子科技集团公司第十研究所) 三维天线阵综合孔径辐射计快速亮温反演方法
CN111538001B (zh) * 2020-03-31 2023-07-21 西南电子技术研究所(中国电子科技集团公司第十研究所) 三维天线阵综合孔径辐射计快速亮温反演方法
CN114812849A (zh) * 2022-03-14 2022-07-29 湖北工业大学 一种基于微波辐射测量的电缆接头内部温度测量方法

Also Published As

Publication number Publication date
CN107607946B (zh) 2020-03-17

Similar Documents

Publication Publication Date Title
Srivastav et al. A highly digital multiantenna ground-penetrating radar (GPR) system
Akiyama et al. First M87 event horizon telescope results. III. Data processing and calibration
Chang et al. Beam calibration of radio telescopes with drones
Stappers et al. Observing pulsars and fast transients with LOFAR
CN100392426C (zh) 单通道相控阵接收信号重构及处理方法
CN101980048A (zh) 基于天线组阵技术的空间碎片地基雷达系统
CN101685159B (zh) 一种星载sar信号高精度保相成像处理平台的构建方法
CN103323845B (zh) 一种非均匀采样综合孔径辐射计的图像反演方法
Domnin et al. Kharkiv incoherent scatter facility
CN104535998B (zh) 旋转镜像综合孔径辐射计及测量方法
Yeary et al. A brief overview of weather radar technologies and instrumentation
CN102520408B (zh) 一种圆柱阵面三维成像系统的三维成像方法
Wang et al. Effects of anisotropic ionospheric irregularities on space-borne SAR imaging
CN103969628A (zh) 一种基于压缩感知的合成孔径雷达pfa成像算法
CN105158735A (zh) 基于压缩采样阵列的空频二维谱估计方法
Reiner et al. Multipoint observations of solar type III radio bursts from STEREO and Wind
CN107607946A (zh) 三维均匀采样综合孔径辐射计亮温反演方法
Chen et al. MST radars of Chinese meridian project: System description and atmospheric wind measurement
CN107144274A (zh) 在轨x射线脉冲星计时模型构建方法
Zhang et al. Development of China’s spaceborne SAR satellite, processing strategy, and application: take Gaofen-3 series as an example
CN112162158B (zh) 在轨太赫兹探测仪天线指向失配的评估方法及装置
Chen et al. Source positions of an interplanetary type III radio burst and anisotropic radio-wave scattering
CN103512673A (zh) 一种微波探测仪观测月球亮温度数据处理方法
CN115586580A (zh) 基于双站信号反演低电离层d层探测方法
Ramos-Perez et al. PAU-SA: A synthetic aperture interferometric radiometer test bed for potential improvements in future missions

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200317

Termination date: 20210724

CF01 Termination of patent right due to non-payment of annual fee