CN107167807A - 三维天线阵列综合孔径辐射计分段图像反演方法 - Google Patents

三维天线阵列综合孔径辐射计分段图像反演方法 Download PDF

Info

Publication number
CN107167807A
CN107167807A CN201710366740.0A CN201710366740A CN107167807A CN 107167807 A CN107167807 A CN 107167807A CN 201710366740 A CN201710366740 A CN 201710366740A CN 107167807 A CN107167807 A CN 107167807A
Authority
CN
China
Prior art keywords
mrow
msub
visual field
subinterval
antenna array
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
CN201710366740.0A
Other languages
English (en)
Other versions
CN107167807B (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 CN201710366740.0A priority Critical patent/CN107167807B/zh
Publication of CN107167807A publication Critical patent/CN107167807A/zh
Application granted granted Critical
Publication of CN107167807B publication Critical patent/CN107167807B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radiation Pyrometers (AREA)

Abstract

本发明公开了三维天线阵列综合孔径辐射计分段图像反演方法。所述方法包括:将观测视场区间划分为N个子区间;在划分的N个视场子区间中分别放入一个外部校正源;依次测量外部校正源在N个视场子区间内三维天线阵列综合孔径辐射计输出的可见度函数,得到可见度函数的相位角;仿真计算三维天线阵列综合孔径辐射计只有两维坐标情况下的可见度函数的相位;依次将N个视场子区间测量的相位角减去仿真得到的相位角,得到第n个视场子区间的相位补偿角;使用三维天线阵列综合孔径辐射计对目标场景进行观测,依次对划分的N个视场子区间进行图像反演;将这N个视场子区间中的反演图像拼接成整个视场内的反演图像。本方法精度较高,具有较好的应用前景。

Description

三维天线阵列综合孔径辐射计分段图像反演方法
技术领域
本发明涉及微波遥感、目标探测以及深空探测技术领域,具体涉及一种三维天线阵列综合孔径辐射计分段图像反演方法。
背景技术
综合孔径辐射计利用多个离散的小天线合成等效的大天线孔径,采用稀疏阵列排布,减少天线的质量和体积,可提高被动微波遥感和探测设备的空间分辨率。为保证反演图像的质量,目前综合孔径辐射计的天线阵列只能排列成一维线性阵列和二维平面阵列。而在目标探测和深空探测等应用中,将天线阵列排列三维阵列是不可避免的。
目前尚未见国内外专门针对三维天线阵列综合孔径辐射计的反演方法的报道和研究。反演精度较差是三维天线阵列综合孔径辐射计实际应用的难题。
发明内容
针对三维天线阵列综合孔径辐射计反演精度较差的缺陷以及改进需求,本发明提供了一种三维天线阵列综合孔径辐射计分段图像反演方法,其目的在于精确地反演三维天线阵列综合孔径辐射计的场景亮温分布图像。
为了达到上述目的,本发明所采用的技术方案是:一种三维天线阵列综合孔径辐射计分段图像反演方法,其特征在于,所述方法包括如下步骤:
S1:根据三维天线阵列中,单元天线的空间位置坐标计算所有基线的三维空间频率;
S2:将观测视场区间划分为N个子区间;
具体数学过程可表示为:观测视场范围为F=[ξmin≤ξ≤ξmaxmin≤η≤ηmax],其中(ξ,η)表示方向余弦,θ表示辐射计观测方向的俯仰角,表示辐射计观测方向的方位角,ξmin和ξmax分别表示ξ在观测视场内的ξ方向上最小值和最大值,ηmin和ηmax分别表示η在观测视场内的η方向上最小值和最大值;
观测视场被划分为N个子区间Fn,n=1,2,3,…,N,且Fn∈F;
S3:在划分的N个视场子区间Fn中分别放入一个外部校正源,其中外校正源的位置为:Pn=(ξnn),n=1,2,…,N,且Pn∈Fn
S4:依次测量外部校正源在N个视场子区间内三维天线阵列综合孔径辐射计输出的可见度函数,从而得到可见度函数的相位角,记为其中,为外部校正源在第n个视场子区间测量可见度,φn为测量可见度的相位角,angle()表示求复数的相位角,或称为复数的辐角;
S5:在划分的N个视场子区间内,外部校正源的位置不变,不考虑系统误差,不考虑三维天线阵列综合孔径辐射计空间频域中的第三维坐标w,仿真计算三维天线阵列综合孔径辐射计只有两维坐标(u,v)情况下的可见度函数的相位,记为Φnn=1,2,…,N;
S6:依次将N个视场子区间测量的相位角减去仿真得到的相位角,即△φn=φnn,△φn称为第n个视场子区间的相位补偿角;
S7:使用三维天线阵列综合孔径辐射计对目标场景进行观测,然后依次对划分的N个视场子区间进行图像反演;
S8:N个视场子区间中的反演图像都可由步骤S7依次得到,最后将这N个视场子区间中的反演图像拼接成整个视场内的反演图像。
进一步地,步骤S1中,具体数学过程可表示为:三维天线阵列中,单元天线的空间位置坐标为(x,y,z),根据公式uk=(xi-xj)/λ、vk=(yi-yj)/λ和wk=(zi-zj)/λ可计算出三维天线阵列综合孔径辐射计的第k个基线的三维空间频率(uk,vk,zk);(xi,yi,zi)表示第i个单元天线的空间位置坐标,(xj,yj,zj)表示第j个单元天线的空间位置坐标,λ为辐射计接收信号的波长;假设三维天线阵列中有M个天线单元,则对应的基线数目为个,所以
进一步地,步骤S2中,视场的子区间可以均匀划分,或者非均匀划分。
进一步地,步骤S5中,具体的计算公式为:
式中表示外部校正源在第n个视场子区间内,计算得到的第k个可见度函数,M表示三维阵列的天线数目,(uk,vk)表示第k个基线的空间频率,T(ξ,η)表示修正亮温,Fi(ξ,η)表示第i个天线的归一化方向图,Torig(ξ,η)表示外部校正源的亮温,因此,第k个可见度函数的相位角为:
进一步地,步骤S7中,具体反演过程如下:
1)将三维天线阵列综合孔径辐射计测量的可见度函数的相位角减去相位补偿角△φn
2)再将第n个视场子区间中观测场景的直流分量从测量可见度中减去,得到交变可见度函数,即
△V(u,v)=V-(u,v)-TDirct·Verr(u,v)
式中△V(u,v)表示交变可见度函数,V-(u,v)表示减去相位补偿角以后的测量可见度函数,TDirct表示观测场景的直流分量,可表示为:
式中V-(0,0)表示测量的零基线,F(ξ,η)表示天线的归一化方向图;
Verr(u,v)表示天线方向图和系统消条纹的误差可见度,可表示为:
式中表示综合孔径辐射计的消条纹函数,f0表示系统的中心频率;
3)对交变可见度进行亮温反演,得到交变反演图像,即,并将△T作为反卷图像Td的初始值,即
4)对交变反演图像△T进行二维傅里叶变换,得到三维天线阵列综合孔径辐射计的仿真可见度函数,再对仿真可见度函数进行亮温反演生成一副中间图像,即
式中Vsim(uk,vk)表示三维天线阵列综合孔径辐射计对观测场景△T仿真得到的可见度函数,Vsim(uk,vk)可表示为:
5)用交变反演图像△T减去γ倍的中间图像T′,得到残差图像Tres=△T-γT′,γ一般为大于0,小于1的常数;并将Tres作为残差图像的初始值,即
6)更新残差图像,即
7)将残差图像加上反卷图像,更新反卷图像,即
8)判断是否满足迭代结束的条件,若满足迭代结束条件则进入9);若不满足迭代结束的条件则跳转至6);迭代结束的条件为:计算残差图像Tres的范数||Tres||,||·||表示计算范数,若||Tres||<stop_flag,则迭代结束,其中stop_flag为迭代停止标识,可根据应用环境自行设定;
9)将γ倍的反卷图像加上观测场景直流分量,得到n个视场子区间的反演图像,即其中Td的值为迭代结束时的反卷图像的值。
进一步地,步骤5)中,γ=0.1。
与现有技术相比,本发明的有益效果是:本发明提出的三维天线阵列综合孔径辐射计分段图像反演方法,解决了三维天线阵列综合孔径辐射计反演精度差的难题,有效提升了三维天线阵列综合孔径辐射计的反演精度,从而推动三维天线阵列综合孔径辐射计的实际应用,使得三维天线阵列综合孔径辐射计在被动微波遥感、目标探测,以及深空探测等领域的实际应用变得可能。本发明提出的算法不仅精度较高,而且可进行并行计算,具有较好的应用前景。
附图说明
图1为三维天线阵列综合孔径辐射计的结构示意图。
图2为本发明实施例提供的一种三维T形阵列的顶视图。
图3为本发明实施例提供的一种三维T形阵列的立体图。
图4为本发明实施例提供的一种三维天线阵列综合孔径辐射计的所有
空间频率分布。
图5为视场子区间的均匀划分。
图6为观测场景。
图7为反演图像。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明,下面结合实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
图1示出了三维天线阵列综合孔径辐射计的结构,为了便于说明,仅示出了与本发明实施例相关的部分,详述如下:
三维天线阵列综合孔径辐射计包括依次连接的天线阵列1、接收通道阵列2、相关器3;其中天线阵列1为三维天线阵列,即阵列中的所有天线单元都位于三维空间中。接收通道阵列2包括多个接收通道,一个接收通道对应一个天线单元,接收通道将天线单元接收到的信号进行下变频、滤波和放大。相关器3将接收通道输出的信号两两进行复相关。复相关后的输出为可见度函数。
本发明实施例提供的基于上述三维天线阵列综合孔径辐射计分段图像反演方法具体包括下述步骤:
(1)根据三维天线阵列中,单元天线的空间位置坐标计算所有基线的三维空间频率。根据计算公式uk=(xi-xj)/λ、vk=(yi-yj)/λ和wk=(zi-zj)/λ可计算出三维天线阵列综合孔径辐射计的第k个基线的三维空间频率(uk,vk,zk);其中,(xi,yi,zi)表示第i个单元天线的空间位置坐标,(xj,yj,zj)表示第j个单元天线的空间位置坐标,λ为辐射计接收信号的波长。如果三维天线阵列中有M个单元天线,则基线数目为个,所以
(2)将观测视场区间划分为N个子区间。
具体数学过程可表示为:观测视场范围为F=[ξmin≤ξ≤ξmaxmin≤η≤ηmax],其中(ξ,η)表示方向余弦,θ表示辐射计观测方向的俯仰角,表示辐射计观测方向的方位角。ξmin和ξmax分别表示ξ在观测视场内的ξ方向上最小值和最大值,ηmin和ηmax分别表示η在观测视场内的η方向上最小值和最大值。
观测视场被划分为N个子区间Fn,n=1,2,3,…,N,且Fn∈F。视场的子区间可均匀划分,也可以非均匀划分。
(3)在划分的N个视场子区间Fn中分别放入一个外部校正源。其中外校正源的位置为:Pn=(ξnn),n=1,2,…,N,且Pn∈Fn
(4)依次测量外部校正源在N个视场子区间内三维天线阵列综合孔径辐射计输出的可见度函数,从而得到可见度函数的相位角,记为其中,为外部校正源在第n个视场子区间测量可见度,φn为测量可见度的相位角,angle()表示求复数的相位角(或称为复数的辐角)。
(5)在划分的N个视场子区间内,外部校正源的位置不变,不考虑系统误差,不考虑三维天线阵列综合孔径辐射计空间频域中的第三维坐标w,仿真计算三维天线阵列综合孔径辐射计只有两维坐标(u,v)情况下的可见度函数的相位,记为Φnn=1,2,…,N。
具体的计算公式为:
式中表示外部校正源在第n个视场子区间内,计算得到的第k个可见度函数,M表示三维阵列的天线数目。(uk,vk)表示第k个基线的空间频率,T(ξ,η)表示修正亮温,Fi(ξ,η)表示第i个天线的归一化方向图,Torig(ξ,η)表示外部校正源的亮温。因此,第k个可见度函数的相位角为:
(6)依次将N个视场子区间测量的相位角减去仿真得到的相位角,即△φn=φnn。△φn称为第n个视场子区间的相位补偿角。
(7)使用三维天线阵列综合孔径辐射计对目标场景进行观测,然后依次对划分的N个视场子区间进行图像反演。具体反演过程如下:
a)将三维天线阵列综合孔径辐射计测量的可见度函数的相位角减去相位补偿角△φn
b)再将第n个视场子区间中观测场景的直流分量从测量可见度中减去,得到交变可见度函数,即
△V(u,v)=V-(u,v)-TDirct·Verr(u,v)
式中△V(u,v)表示交变可见度函数,V-(u,v)表示减去相位补偿角以后的测量可见度函数,TDirct表示观测场景的直流分量,可表示为:
式中V-(0,0)表示测量的零基线,F(ξ,η)表示天线的归一化方向图。
Verr(u,v)表示天线方向图和系统消条纹的误差可见度,可表示为:
式中表示综合孔径辐射计的消条纹函数,f0表示系统的中心频率。
c)对交变可见度进行亮温反演,得到交变反演图像。即,并将△T作为反卷图像Td的初始值,即
d)对交变反演图像△T进行二维傅里叶变换,得到三维天线阵列综合孔径辐射计的仿真可见度函数,再对仿真可见度函数进行亮温反演生成一副中间图像,即
式中Vsim(uk,vk)表示三维天线阵列综合孔径辐射计对观测场景△T仿真得到的可见度函数,Vsim(uk,vk)可表示为:
e)用交变反演图像△T减去γ倍的中间图像T′,得到残差图像Tres=△T-γT′,γ一般为大于0,小于1的常数。通常γ=0.1。并将Tres作为残差图像的初始值,即
f)更新残差图像,即
g)将残差图像加上反卷图像,更新反卷图像,即
h)判断是否满足迭代结束的条件,若满足迭代结束条件则进入i)。若不满足迭代结束的条件则跳转至f)。迭代结束的条件为:计算残差图像Tres的范数||Tres||,||·||表示计算范数。若||Tres||<stop_flag,则迭代结束,其中stop_flag为迭代停止标识,可根据应用环境自行设定。
i)将γ倍的反卷图像加上观测场景直流分量,得到n个视场子区间的反演图像,即其中Td的值为迭代结束时的反卷图像的值。
(8)N个视场子区间中的反演图像都可由步骤(7)依次得到,最后将这N个视场子区间中的反演图像拼接成整个视场内的反演图像。
实施例:三维天线阵列的综合孔径辐射计系统;该实施例中,天线阵为一个三维T形阵列。其所有单元天线分布在三维空间中,图2示出该三维天线阵列的顶视图(即从上往下看的正投影),图3示出该三维天线阵列的三维立体图。这两个图中的每个小圆圈代表一个单元天线。
图像反演方法具体步骤如下:
(1)根据三维天线阵列中,单元天线的空间位置坐标计算所有基线的三维空间频率。根据计算公式uk=(xi-xj)/λ、vk=(yi-yj)/λ和wk=(zi-zj)/λ可计算出三维天线阵列综合孔径辐射计的第k个基线的三维空间频率(uk,vk,zk);其中,(xi,yi,zi)表示第i个单元天线的空间位置坐标,(xj,yj,zj)表示第j个单元天线的空间位置坐标,λ为辐射计接收信号的波长。图4示出了三维天线阵列所有的空间频率分布。
(2)将观测视场区间F划分为9个子区间,即Fn,n=1,2,3,…,9,且Fn∈F,F1∪F2…∪F9∈F。此例中观测视场的子区间为均匀划分。如图5所示。
(3)在划分的9个视场子区间中依次放入一个外部校正源。此例中,放入的外部校正源为一个噪声源。
(4)依次测量噪声源在这9个视场子区间内三维T形阵列输出的可见度函数,从而得到可见度函数的相位角,记为其中,为外部校正源在第n个视场子区间测量可见度,φn为测量可见度的相位角。
(5)在划分的9个视场子区间内,噪声源的位置不变,不考虑系统误差,不考虑三维天线阵列综合孔径辐射计空间频域中的第三维坐标w,仿真计算三维天线阵列综合孔径辐射计只有两维坐标(u,v)情况下的可见度函数的相位,记为Φn n=1,2,…,N。
具体的计算公式为:
式中表示外部校正源在第n个视场子区间内,计算得到的第k个可见度函数,M表示三维阵列的天线数目。(uk,vk)表示第k个基线的空间频率,T(ξ,η)表示修正亮温,Fi(ξ,η)表示第i个天线的归一化方向图,δ(ξ-ξ′,η-η′)表示噪声源的亮温。因此,第k个可见度函数的相位角为:
(6)依次将9个视场子区间测量的相位角减去仿真得到的相位角,即△φn=φnn。△φn称为第n个视场子区间的相位补偿角。
(7)使用此三维T形阵列,对图6所示的目标场景进行观测,然后依次对划分的9个视场子区间进行图像反演。具体反演过程如下:
a)将三维天线阵列综合孔径辐射计测量的可见度函数的相位角减去相位补偿角△φn
b)再将第n个视场子区间中观测场景的直流分量从测量可见度中减去,得到交变可见度函数,即
△V(u,v)=V-(u,v)-TDirct·Verr(u,v)
式中△V(u,v)表示交变可见度函数,V-(u,v)表示减去相位补偿角以后的测量可见度函数,TDirct表示观测场景的直流分量,可表示为:
式中V-(0,0)表示测量的零基线,F(ξ,η)表示天线的归一化方向图。
Verr(u,v)表示天线方向图和系统消条纹的误差可见度,可表示为:
式中表示综合孔径辐射计的消条纹函数,f0表示系统的中心频率,在本例中中心频率为1.4GHz。
c)对交变可见度进行亮温反演,得到交变反演图像。即,并将△T作为反卷图像Td的初始值,即
d)对交变反演图像△T进行二维傅里叶变换,得到三维天线阵列综合孔径辐射计的仿真可见度函数,再对仿真可见度函数进行亮温反演生成一副中间图像,即
式中Vsim(uk,vk)表示三维天线阵列综合孔径辐射计对观测场景△T仿真得到的可见度函数,Vsim(uk,vk)可表示为:
e)用交变反演图像△T减去γ倍的中间图像T′,得到残差图像Tres=△T-γT′,在本例中,γ=0.1。并将Tres作为残差图像的初始值,即
f)更新残差图像,即
g)将残差图像加上反卷图像,更新反卷图像,即
h)判断是否满足迭代结束的条件,若满足迭代结束条件则进入i)。若不满足迭代结束的条件则跳转至f)。迭代结束的条件为:计算残差图像Tres的范数||Tres||,||·||表示计算范数。若||Tres||<stop_flag,则迭代结束,其中stop_flag为迭代停止标识,在本例中stop_flag=0.01。
i)将γ倍的反卷图像加上观测场景直流分量,得到n个视场子区间的反演图像,即其中Td的值为迭代结束时的反卷图像的值。在本例中,总的迭代次数为45次。
(8)9个视场子区间中的反演图像都可由步骤(7)依次得到,最后将这9个视场子区间中的反演图像拼接成整个视场内的反演图像,如图7所示。
应当理解的是,本说明书未详细阐述的部分均属于现有技术。
应当理解的是,上述针对较佳实施例的描述较为详细,并不能因此而认为是对本发明专利保护范围的限制,本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的范围情况下,还可以做出替换或变形,均落入本发明的保护范围之内,本发明的请求保护范围应以所附权利要求为准。

Claims (6)

1.一种三维天线阵列综合孔径辐射计分段图像反演方法,其特征在于,所述方法包括如下步骤:
S1:根据三维天线阵列中,单元天线的空间位置坐标计算所有基线的三维空间频率;
S2:将观测视场区间划分为N个子区间;
具体数学过程可表示为:观测视场范围为F=[ξmin≤ξ≤ξmaxmin≤η≤ηmax],其中(ξ,η)表示方向余弦,θ表示辐射计观测方向的俯仰角,表示辐射计观测方向的方位角,ξmin和ξmax分别表示ξ在观测视场内的ξ方向上最小值和最大值,ηmin和ηmax分别表示η在观测视场内的η方向上最小值和最大值;
观测视场被划分为N个子区间Fn,n=1,2,3,…,N,且Fn∈F;
S3:在划分的N个视场子区间Fn中分别放入一个外部校正源,其中外校正源的位置为:Pn=(ξnn),n=1,2,…,N,且Pn∈Fn
S4:依次测量外部校正源在N个视场子区间内三维天线阵列综合孔径辐射计输出的可见度函数,从而得到可见度函数的相位角,记为其中,为外部校正源在第n个视场子区间测量可见度,φn为测量可见度的相位角,angle()表示求复数的相位角,或称为复数的辐角;
S5:在划分的N个视场子区间内,外部校正源的位置不变,不考虑系统误差,不考虑三维天线阵列综合孔径辐射计空间频域中的第三维坐标w,仿真计算三维天线阵列综合孔径辐射计只有两维坐标(u,v)情况下的可见度函数的相位,记为Φn n=1,2,…,N;
S6:依次将N个视场子区间测量的相位角减去仿真得到的相位角,即△φn=φnn,△φn称为第n个视场子区间的相位补偿角;
S7:使用三维天线阵列综合孔径辐射计对目标场景进行观测,然后依次对划分的N个视场子区间进行图像反演;
S8:N个视场子区间中的反演图像都可由步骤S7依次得到,最后将这N个视场子区间中的反演图像拼接成整个视场内的反演图像。
2.如权利要求1所述的三维天线阵列综合孔径辐射计分段图像反演方法,其特征在于,步骤S1中,具体数学过程可表示为:三维天线阵列中,单元天线的空间位置坐标为(x,y,z),根据公式uk=(xi-xj)/λ、vk=(yi-yj)/λ和wk=(zi-zj)/λ可计算出三维天线阵列综合孔径辐射计的第k个基线的三维空间频率(uk,vk,zk);(xi,yi,zi)表示第i个单元天线的空间位置坐标,(xj,yj,zj)表示第j个单元天线的空间位置坐标,λ为辐射计接收信号的波长;假设三维天线阵列中有M个天线单元,则对应的基线数目为个,所以
3.如权利要求1所述的三维天线阵列综合孔径辐射计分段图像反演方法,其特征在于,步骤S2中,视场的子区间可以均匀划分,或者非均匀划分。
4.如权利要求1所述的三维天线阵列综合孔径辐射计分段图像反演方法,其特征在于,步骤S5中,具体的计算公式为:
<mrow> <msubsup> <mi>V</mi> <mi>k</mi> <mi>n</mi> </msubsup> <mo>=</mo> <mi>&amp;Delta;</mi> <mi>S</mi> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> </munder> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> </munder> <mi>T</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;xi;</mi> <mi>n</mi> </msub> <mo>,</mo> <msub> <mi>&amp;eta;</mi> <mi>m</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>k</mi> </msub> <msub> <mi>&amp;xi;</mi> <mi>n</mi> </msub> <mo>+</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <msub> <mi>&amp;eta;</mi> <mi>m</mi> </msub> <mo>)</mo> </mrow> </mrow> </msup> </mrow>
式中表示外部校正源在第n个视场子区间内,计算得到的第k个可见度函数,M表示三维阵列的天线数目,(uk,vk)表示第k个基线的空间频率,T(ξ,η)表示修正亮温,Fi(ξ,η)表示第i个天线的归一化方向图,Torig(ξ,η)表示外部校正源的亮温,因此,第k个可见度函数的相位角为:
5.如权利要求1所述的三维天线阵列综合孔径辐射计分段图像反演方法,其特征在于,步骤S7中,具体反演过程如下:
1)将三维天线阵列综合孔径辐射计测量的可见度函数的相位角减去相位补偿角△φn
2)再将第n个视场子区间中观测场景的直流分量从测量可见度中减去,得到交变可见度函数,即
△V(u,v)=V-(u,v)-TDirct·Verr(u,v)
式中△V(u,v)表示交变可见度函数,V-(u,v)表示减去相位补偿角以后的测量可见度函数,TDirct表示观测场景的直流分量,可表示为:
<mrow> <msub> <mi>T</mi> <mrow> <mi>D</mi> <mi>i</mi> <mi>r</mi> <mi>c</mi> <mi>t</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msup> <mi>V</mi> <mo>-</mo> </msup> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mn>0</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;Integral;</mo> <msub> <mo>&amp;Integral;</mo> <msub> <mi>F</mi> <mi>n</mi> </msub> </msub> <msup> <mrow> <mo>|</mo> <msub> <mi>F</mi> <mi>n</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;xi;</mi> <mo>,</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mi>d</mi> <mi>&amp;xi;</mi> <mi>d</mi> <mi>&amp;eta;</mi> </mrow> </mfrac> </mrow>
式中V-(0,0)表示测量的零基线,F(ξ,η)表示天线的归一化方向图;
Verr(u,v)表示天线方向图和系统消条纹的误差可见度,可表示为:
<mrow> <msub> <mi>V</mi> <mrow> <mi>e</mi> <mi>r</mi> <mi>r</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>&amp;Integral;</mo> <msub> <mo>&amp;Integral;</mo> <msub> <mi>F</mi> <mi>n</mi> </msub> </msub> <msub> <mi>F</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;xi;</mi> <mo>,</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <msup> <msub> <mi>F</mi> <mi>j</mi> </msub> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mi>&amp;xi;</mi> <mo>,</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <mover> <mi>r</mi> <mo>~</mo> </mover> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mrow> <mi>u</mi> <mi>&amp;xi;</mi> <mo>+</mo> <mi>v</mi> <mi>&amp;eta;</mi> </mrow> <msub> <mi>f</mi> <mn>0</mn> </msub> </mfrac> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <mi>u</mi> <mi>&amp;xi;</mi> <mo>+</mo> <mi>v</mi> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> </mrow> </msup> <mi>d</mi> <mi>&amp;xi;</mi> <mi>d</mi> <mi>&amp;eta;</mi> </mrow>
式中表示综合孔径辐射计的消条纹函数,f0表示系统的中心频率;
3)对交变可见度进行亮温反演,得到交变反演图像,即,并将△T作为反卷图像Td的初始值,即
4)对交变反演图像△T进行二维傅里叶变换,得到三维天线阵列综合孔径辐射计的仿真可见度函数,再对仿真可见度函数进行亮温反演生成一副中间图像,即
<mrow> <msup> <mi>T</mi> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> </munder> <msub> <mi>V</mi> <mrow> <mi>s</mi> <mi>i</mi> <mi>m</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>k</mi> </msub> <mo>,</mo> <msub> <mi>v</mi> <mi>k</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>k</mi> </msub> <mi>&amp;xi;</mi> <mo>+</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> </mrow> </msup> </mrow>
式中Vsim(uk,vk)表示三维天线阵列综合孔径辐射计对观测场景△T仿真得到的可见度函数,Vsim(uk,vk)可表示为:
<mrow> <msub> <mi>V</mi> <mrow> <mi>s</mi> <mi>i</mi> <mi>m</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>k</mi> </msub> <mo>,</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mo>&amp;Integral;</mo> <msub> <mo>&amp;Integral;</mo> <msub> <mi>F</mi> <mi>n</mi> </msub> </msub> <msub> <mi>&amp;Delta;TF</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;xi;</mi> <mo>,</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <msubsup> <mi>F</mi> <mi>j</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mi>&amp;xi;</mi> <mo>,</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <mover> <mi>r</mi> <mo>~</mo> </mover> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mrow> <mi>u</mi> <mi>&amp;xi;</mi> <mo>+</mo> <mi>v</mi> <mi>&amp;eta;</mi> </mrow> <msub> <mi>f</mi> <mn>0</mn> </msub> </mfrac> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>k</mi> </msub> <mi>&amp;xi;</mi> <mo>+</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> </mrow> </msup> <mi>d</mi> <mi>&amp;xi;</mi> <mi>d</mi> <mi>&amp;eta;</mi> </mrow>
5)用交变反演图像△T减去γ倍的中间图像T′,得到残差图像Tres=△T-γT′,γ一般为大于0,小于1的常数;并将Tres作为残差图像的初始值,即
6)更新残差图像,即
<mrow> <msubsup> <mi>T</mi> <mrow> <mi>r</mi> <mi>e</mi> <mi>s</mi> </mrow> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <msubsup> <mi>T</mi> <mrow> <mi>r</mi> <mi>e</mi> <mi>s</mi> </mrow> <mi>k</mi> </msubsup> <mo>-</mo> <mi>&amp;gamma;</mi> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> </munder> <mo>&amp;lsqb;</mo> <mo>&amp;Integral;</mo> <msub> <mo>&amp;Integral;</mo> <msub> <mi>F</mi> <mi>n</mi> </msub> </msub> <msubsup> <mi>T</mi> <mrow> <mi>r</mi> <mi>e</mi> <mi>s</mi> </mrow> <mi>k</mi> </msubsup> <msub> <mi>F</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;xi;</mi> <mo>,</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <msubsup> <mi>F</mi> <mi>j</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mi>&amp;xi;</mi> <mo>,</mo> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> <mover> <mi>r</mi> <mo>~</mo> </mover> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mrow> <mi>u</mi> <mi>&amp;xi;</mi> <mo>+</mo> <mi>v</mi> <mi>&amp;eta;</mi> </mrow> <msub> <mi>f</mi> <mn>0</mn> </msub> </mfrac> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>k</mi> </msub> <mi>&amp;xi;</mi> <mo>+</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> </mrow> </msup> <mi>d</mi> <mi>&amp;xi;</mi> <mi>d</mi> <mi>&amp;eta;</mi> <mo>&amp;rsqb;</mo> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>k</mi> </msub> <mi>&amp;xi;</mi> <mo>+</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mi>&amp;eta;</mi> <mo>)</mo> </mrow> </mrow> </msup> </mrow>
7)将残差图像加上反卷图像,更新反卷图像,即
8)判断是否满足迭代结束的条件,若满足迭代结束条件则进入9);若不满足迭代结束的条件则跳转至6);迭代结束的条件为:计算残差图像Tres的范数||Tres||,||·||表示计算范数,若||Tres||<stop_flag,则迭代结束,其中stop_flag为迭代停止标识,可根据应用环境自行设定;
9)将γ倍的反卷图像加上观测场景直流分量,得到n个视场子区间的反演图像,即其中Td的值为迭代结束时的反卷图像的值。
6.如权利要求5所述的三维天线阵列综合孔径辐射计分段图像反演方法,其特征在于,步骤5)中,γ=0.1。
CN201710366740.0A 2017-05-23 2017-05-23 三维天线阵列综合孔径辐射计分段图像反演方法 Expired - Fee Related CN107167807B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710366740.0A CN107167807B (zh) 2017-05-23 2017-05-23 三维天线阵列综合孔径辐射计分段图像反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710366740.0A CN107167807B (zh) 2017-05-23 2017-05-23 三维天线阵列综合孔径辐射计分段图像反演方法

Publications (2)

Publication Number Publication Date
CN107167807A true CN107167807A (zh) 2017-09-15
CN107167807B CN107167807B (zh) 2019-10-11

Family

ID=59816513

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710366740.0A Expired - Fee Related CN107167807B (zh) 2017-05-23 2017-05-23 三维天线阵列综合孔径辐射计分段图像反演方法

Country Status (1)

Country Link
CN (1) CN107167807B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110515048A (zh) * 2019-08-07 2019-11-29 中国科学院国家空间科学中心 一种旋转采样综合孔径辐射计相位定标方法
CN111538001A (zh) * 2020-03-31 2020-08-14 西南电子技术研究所(中国电子科技集团公司第十研究所) 三维天线阵综合孔径辐射计快速亮温反演方法
CN111597730A (zh) * 2020-05-26 2020-08-28 华中科技大学 一种一维综合孔径辐射计天线阵列排布方法
CN111929651A (zh) * 2020-07-30 2020-11-13 西安空间无线电技术研究所 一种综合孔径辐射计系统相位残差标定方法
CN113687446A (zh) * 2021-07-19 2021-11-23 西安空间无线电技术研究所 一种基于星载综合孔径微波辐射计的近场外定标方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101975947A (zh) * 2010-10-22 2011-02-16 华中科技大学 二维镜像综合孔径辐射成像方法
CN102621532A (zh) * 2012-03-31 2012-08-01 华中科技大学 基于阵列旋转的综合孔径辐射计可见度相位误差校正方法
CN105738851A (zh) * 2016-02-02 2016-07-06 华中科技大学 一种综合孔径微波辐射计联合校正方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101975947A (zh) * 2010-10-22 2011-02-16 华中科技大学 二维镜像综合孔径辐射成像方法
CN102621532A (zh) * 2012-03-31 2012-08-01 华中科技大学 基于阵列旋转的综合孔径辐射计可见度相位误差校正方法
CN105738851A (zh) * 2016-02-02 2016-07-06 华中科技大学 一种综合孔径微波辐射计联合校正方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LI FENG 等: "ARRAY CONFIGURATION OPTIMIZATION FOR ONE-DIMENSIONAL NONUNIFORM APERTURE SYNTHESIS RADIOMETERS", 《2016 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM》 *
李育芳 等: "旋转镜像综合孔径空间频率均匀采样方法", 《微波学报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110515048A (zh) * 2019-08-07 2019-11-29 中国科学院国家空间科学中心 一种旋转采样综合孔径辐射计相位定标方法
CN110515048B (zh) * 2019-08-07 2021-04-09 中国科学院国家空间科学中心 一种旋转采样综合孔径辐射计相位定标方法
CN111538001A (zh) * 2020-03-31 2020-08-14 西南电子技术研究所(中国电子科技集团公司第十研究所) 三维天线阵综合孔径辐射计快速亮温反演方法
CN111538001B (zh) * 2020-03-31 2023-07-21 西南电子技术研究所(中国电子科技集团公司第十研究所) 三维天线阵综合孔径辐射计快速亮温反演方法
CN111597730A (zh) * 2020-05-26 2020-08-28 华中科技大学 一种一维综合孔径辐射计天线阵列排布方法
CN111597730B (zh) * 2020-05-26 2024-03-19 华中科技大学 一种一维综合孔径辐射计天线阵列排布方法
CN111929651A (zh) * 2020-07-30 2020-11-13 西安空间无线电技术研究所 一种综合孔径辐射计系统相位残差标定方法
CN111929651B (zh) * 2020-07-30 2024-03-15 西安空间无线电技术研究所 一种综合孔径辐射计系统相位残差标定方法
CN113687446A (zh) * 2021-07-19 2021-11-23 西安空间无线电技术研究所 一种基于星载综合孔径微波辐射计的近场外定标方法
CN113687446B (zh) * 2021-07-19 2024-02-09 西安空间无线电技术研究所 一种基于星载综合孔径微波辐射计的近场外定标方法

Also Published As

Publication number Publication date
CN107167807B (zh) 2019-10-11

Similar Documents

Publication Publication Date Title
CN107167807A (zh) 三维天线阵列综合孔径辐射计分段图像反演方法
US7119739B1 (en) Near field to far field DF antenna array calibration technique
Nelles et al. A parameterization for the radio emission of air showers as predicted by CoREAS simulations and applied to LOFAR measurements
Reisswig et al. Unambiguous determination of gravitational waveforms from binary black hole mergers
CN106291129A (zh) 相控阵天线远场方向图测量方法
CN102073035B (zh) 星载干涉合成孔径雷达系统仿真性能评估方法
CN103323845B (zh) 一种非均匀采样综合孔径辐射计的图像反演方法
CN106842203B (zh) 一种三维天线阵列综合孔径辐射计的图像反演方法
CN103235301A (zh) 基于复数域平差理论的POLInSAR植被高度反演方法
CN104063544B (zh) 近场效应误差分析方法
CN103091660A (zh) 基于相位跳变修正的虚拟基线测向法
CN104569625B (zh) 一种基于可转动辅助天线的大型天线方向图测量方法
CN105136073B (zh) 一种在边坡形变监测中的气象校正模型
CN104391183A (zh) 近场测量天线远区场特性快速计算方法
CN109375198B (zh) 一种低频声呐阵阻抗特性校准方法
Zilles et al. Radio Morphing: towards a fast computation of the radio signal from air showers
CN105388449A (zh) 一种衡量天线罩对天线阵列测向性能影响的方法
Noorikuhani et al. Wide-angle and relativistic effects in Fourier-space clustering statistics
CN109085556A (zh) 一种基于一二阶峰比值的高频地波雷达浪场形成方法
CN106168662A (zh) 基于极大似然估计的被动传感器的误差配准方法及装置
CN113030847A (zh) 一种用于双通道测向系统的深度学习数据集生成方法
CN108061920A (zh) 探地雷达建模的方法
Fujii et al. High-precision three-dimensional field mapping of a high resolution magnetic spectrometer for hypernuclear spectroscopy at JLab
CN110208741B (zh) 一种基于多圆阵测相的超视距单目标直接定位方法
CN107390296B (zh) 一种用于被动合成孔径成像仪系统的快速标定方法

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: 20191011

Termination date: 20210523

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