CN111538001B - 三维天线阵综合孔径辐射计快速亮温反演方法 - Google Patents
三维天线阵综合孔径辐射计快速亮温反演方法 Download PDFInfo
- Publication number
- CN111538001B CN111538001B CN202010243267.9A CN202010243267A CN111538001B CN 111538001 B CN111538001 B CN 111538001B CN 202010243267 A CN202010243267 A CN 202010243267A CN 111538001 B CN111538001 B CN 111538001B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- antenna array
- fanyan
- wen
- aperture radiometer
- 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
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
- G01S13/904—SAR modes
-
- 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
- G01S13/9004—SAR image acquisition techniques
- G01S13/9011—SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth
Abstract
本发明公开的一种三维天线阵综合孔径辐射计快速亮温反演方法,旨在解决当前三维天线阵反演精度差、视场范围小的问题。本发明通过下述技术方案实现:首先对K元天线阵进行三维空间排布,读取阵列参数;利用三维天线阵的位置坐标计算基线bx,by,bz;其次构造投影函数F(bz);然后根据相关天线的输出结果,获取视场中亮温空间分布视函数S(bx,by,bz),由S(bx,by,bz)和F(bz),计算SF(bx,by);通过对SF(bx,by)进行逆傅里叶变换,获取视场空间范围内亮温分布图像。本发明可以实现亮温的快算反演,改善三维SAIR反演精度差,视场小的问题;同时,具有很强的工程适应性。
Description
技术领域
本发明涉及遥感、探测技术领域,更具体地,涉及一种三维天线阵综合孔径辐射计SAIR (synthetic aperture interferometric radiometer)亮温反演方法。
背景技术
综合孔径辐射计SAIR是一种采用稀疏天线阵列的辐射计系统,它通过小孔径阵列天线等效合成大天线,系统中的小天线阵元输出信号经过两两复相关处理后可以获取视场中亮温空间分布的谱分量(也称视函数),进而通过反演算法,由谱域采样值重构亮温分布图像。综合孔径辐射计可实现对较大视场范围和高空间分辨率的亮温分布成像,同时采用较小的天线单元和分布式稀疏阵列取代大口径窄波束定向天线,降低了天线加工、安装及扫描伺服系统的难度,在高分辨率空间微波遥感、大气遥感、射电天文观测及基于亮温分布的辐射计探测系统中具有广泛的应用。
传统的一维线阵和二维U型、T型、L型、Y型及圆形平面阵虽然在一定程度上可以满足反演图像的质量要求,但是其布阵方式存在较大的局限性。三维综合孔径辐射计突破了二维综合孔径辐射计对天线阵元平面分布的限制,使得天线阵配置更加灵活,可增加对谱域采样密度,同时也消除了平面天线阵在偏离阵面法向时分辨率降低的限制,增大了辐射计视场范围。对于定向天线组阵,由于天线体积大,阵面尺寸很大,阵元分布很难保证在同一平面上,实际上应视为三维阵列。实际的天线布阵不可能是理想的一维线阵或者二维平面阵,天线元在垂直天线阵平面的位移等引起的幅相不一致都可归入接收通道传输函数的幅相误差内,天线元位置的误差将使基线数值产生误差,引起重构的SAIR图像出现畸变。
发明内容
本发明的目的是针对现有技术存在的问题和不足之处,提供一种反演速度快,误差小,反演精度高的三维天线阵综合孔径辐射计快速亮温反演方法,解决当前三维天线阵反演精度差、视场范围小的问题。
本发明的上述目的可以通过以下措施来达到,一种三维天线阵综合孔径辐射计快速亮温反演方法,具有如下技术特征:在综合孔径辐射计SAIR的天线阵中,设置K元天线阵的三维空间排布,建立相应的xoy坐标位置关系;读取阵列参数,根据阵列配置和任意两幅相关天线的位置坐标(x i,y i,z i)、(x j,y j,z j),计算三维基线bx,by,bz,其中,bx= (x i-x j)/λ 0,by= (y i-y j)/λ 0,bz= (z i-z j)/λ 0,λ 0为中心频率对应的波长;基于冲击函数δ构造投影函数表达式:
,其中,e为自然常数,j为虚数单位,δ<<1;根据相关天线的输出结果,获取视场中三维视函数S(bx,by,bz),利用S(bx,by,bz)和投影函数F(bz)计算二维视函数SF(bx,by),然后对二维视函数SF(bx,by)进行逆傅里叶变换IFFT,获得视场空间范围内亮温分布图像,实现综合孔径辐射计亮温反演。
本发明相比于现有技术具有如下有益效果:
反演速度快。本发明采用基于傅里叶换理论及相应的快速算法,对相关天线测量得到的可视度函数进行频域采样,实现综合孔径辐射计亮温分布图像的准确、快速反演计算。传统的三维SAIR反演算法一般会涉及收敛过程,难以实现快速反演;而本发明专利提出的SAIR反演算法并未涉及收敛过程,执行完一次算法就可以得到亮温的反演图像。
反演精度高。本发明的三维综合孔径辐射计采用三维天线阵,每个阵元独立接收空间辐射信号,任意两个相关天线的输出结果可以得到空间亮温分布的一个谱域采样值。当天线阵配置满足一定条件,可以得到近似均匀的谱域样本集,基于冲击函数δ的宽频谱特点,在进行FFT/IFFT变化时不会丢失频域分量,可以改善三维SAIR反演精度差,视场小的问题,从而克服传统基于收敛算法三维SAIR反演在进行FFT/IFFT变换时存在频谱丢失的问题。
适应范围宽。本发明专利提出的SAIR反演算法,具有很强的工程适应性,可以适应各种形式的三维T阵,U阵,十字阵,Y阵,以及规则或者不规则的天线分布。
附图说明
图1为三维天线阵综合孔径辐射计SAIR快速亮温反演方法的天线阵排列示意图。
图2为三维天线阵综合孔径辐射计SAIR快速亮温反演方法的基线模式分布示意图。
图3本发明的三维天线阵综合孔径辐射计快速亮温反演方法的算法流程图。
图4是点目标位置坐标为(50,50)的反演结果示意图。
图5是亮温反演结果演示图,其中左图为原始图像,右图为三维天线阵反演后的图像。
具体实施方式
参阅图1-3。根据本发明,在综合孔径辐射计SAIR的天线阵中,设置K元天线阵的三维空间排布,建立相应的xoy坐标位置关系;读取阵列参数,根据阵列配置和任意两幅相关天线的位置坐标(x i,y i,z i)、(x j,y j,z j),计算三维基线bx,by,bz,其中,bx= (x i-x j)/λ 0,by=(y i-y j)/λ 0,bz= (z i-z j)/λ 0;基于冲击函数δ构造投影函数表达式:
,其中,e为自然常数,j为虚数单位,δ<<1;根据相关天线的输出结果,获取视场中三维视函数S(bx,by,bz),利用S(bx,by,bz)和投影函数F(bz)计算二维视函数SF(bx,by),同时对二维视函数SF(bx,by)进行逆傅里叶变换IFFT,获得视场空间范围内亮温分布图像,实现综合孔径辐射计亮温反演。
在可选的实施例中,以K元T形天线阵为例,在xoy平面投影中,x轴方向上的天线阵元数为M,y轴正方向上的天线阵元数为N,且满足M+N-1 = K。
如附图1所示,三维T阵SAIR的天线阵共有7个天线单元(图中单元位置用黑色球体表示)。SAIR任意两个天线阵元的相对位置矢量与中心频率对应的波长λ 0之比是二维空间频率取样点,称为基线。在可选的实施例中,三维T阵SAIR的基线如附图2所示。
三维视函数S(bx,by,bz)可以通过投影函数F(bz)变换为二维视函数SF(bx,by),并且三维视函数S(bx,by,bz)与投影函数F(bz)的乘积为:。
三维视函数的表达式为
其中,为玻尔兹曼常数,μ,ν,γ均为给定坐标系(x,y,z)下的方向余弦,/>为经过接收机滤波后中心频率f 0附近的亮温,/>为天线单元的功率方向图。
在t时刻,空间坐标(μ,ν,γ)方向辐射经过接收机滤波后中心频率f 0的毫米波信号可以表示为,/>的功率谱密度为/>,且,
其中,μ,ν,γ均为给定坐标系(x,y,z)下的方向余弦,θ为球坐标系下位置矢量与z轴正方向的夹角,φ为球坐标系下位置矢量在xoy平面投影矢量与x轴正方向的夹角,为玻尔兹曼常数,/>为频率f 0附近的亮温。
若t时刻的辐射信号的坐标位置为,/>时刻的坐标位置为/>,则不同方向、不同时刻的辐射信号不相关,即/>,其中/>为/>的共轭信号。
t时刻,天线1(x 1,y 1,z 1)和天线2(x 2,y 2,z 2)收到的空间毫米波辐射信号分别为
其中,和/>分别是天线1、2在μv坐标系下的电压方向图,积分式中的指数相代表信号相对于坐标原点处信号的相移。
将收到的毫米波辐射信号和/>进行复相关运算,可以得到
,令
,
,
,
,
,
其中,为/>的共轭信号,/>为/>的共轭电压方向图,为天线单元的功率方向图。则,,因此,三维视函数/>是的傅里叶变换,或一对天线接收信号的复相关运算可产生的一个谱分量。当对谱域有足够的采样点后,可以通过逆傅里叶变换重构空间亮温分布;略去常数因子/>和/>,上式三维视函数可以表示为
,由于含有γ,为消去γ,进行加权平均运算,两边同时乘以投影函数F(bz),并对基线bz积分,则有二维视函数
为了反演得到中心频率f 0附近的亮温并在不降低分辨率的前提下消除函数/>的影响,构造如下的函数,/>则,
其中,δ<<1,例如可取δ=0.001。
对f(γ)进行逆傅里叶变换,可以得到投影函数F(bz),
在三维稀疏阵的情况下,通常只在/>非零,故
,对/>进行傅里叶逆变换,可以得到中心频率f 0附近的亮温/>。
针对如附图1所示的T阵在三维空间中的一种配置方式,计算附图1所对应的SAIR的基线模式;其次,构造投影函数F(bz)。由相关接收天线得到三维视函数S(bx,by,bz),根据三维视函数S(bx,by,bz)和投影函数F(bz),计算二维视函数,即,,对二维视函数/>进行IFFT2变换,得到反演亮温。除在γ= 0附近外可实现近似良好的亮温分布反演;在γ= 1附近处于视场边沿,天线增益趋于零,从而存在较大误差,去掉这一段反演值即可。至此算法实现了空间亮温反演。
以上为本发明的一种优选实施方式,但本发明的保护范围并不受限于此,任何熟悉本领域的技术人员在本发明创新构思的前提下所做出的变化或替换,均在本发明的保护范围之内。
Claims (10)
1.一种三维天线阵综合孔径辐射计快速亮温反演方法,具有如下技术特征:在综合孔径辐射计SAIR的天线阵中,设置K元天线阵的三维空间排布,建立相应的xoy坐标位置关系;读取阵列参数,根据阵列配置和任意两幅相关天线的位置坐标(x i,y i,z i)、(x j,y j,z j),计算三维基线bx,by,bz,其中,bx = (x i-x j)/λ 0,by = (y i-y j)/λ 0,bz = (z i-z j)/λ 0,λ 0为中心频率对应的波长;基于冲击函数δ构造投影函数表达式:
,其中,e为自然常数,j为虚数单位,δ<<1;根据相关天线的输出结果,获取视场中三维视函数S(bx,by,bz),利用S(bx,by,bz)和投影函数F(bz)计算二维视函数SF(bx,by),然后对二维视函数SF(bx,by)进行逆傅里叶变换IFFT,获得视场空间范围内亮温分布图像,实现综合孔径辐射计亮温反演。
2.如权利要求1所述的三维天线阵综合孔径辐射计快速亮温反演方法,其特征在于,在xoy平面投影中,x轴方向上的天线阵元数为M,y轴方向上的天线阵元数为N,且满足M+N-1 =K。
3.如权利要求1所述的三维天线阵综合孔径辐射计快速亮温反演方法,其特征在于,综合孔径辐射计任意两个天线阵元的相对位置矢量与中心频率f 0对应的波长λ 0之比是二维空间频率取样点,称为基线。
4.如权利要求1所述的三维天线阵综合孔径辐射计快速亮温反演方法,其特征在于,三维视函数S(bx,by,bz)通过投影函数F(bz)变换为二维视函数SF(bx,by),并且三维视函数S(bx,by,bz)与投影函数F(bz)的乘积为:。
5.如权利要求1所述的三维天线阵综合孔径辐射计快速亮温反演方法,其特征在于,三维视函数的表达式为
其中,为玻尔兹曼常数,μ,ν,γ为给定坐标系(x,y,z)下的方向余弦,/>为经过接收机滤波后中心频率f 0附近的亮温,/>为天线单元的功率方向图。
6.如权利要求4所述的三维天线阵综合孔径辐射计快速亮温反演方法,其特征在于,在t时刻,空间坐标(μ,ν,γ)方向辐射经过接收机滤波后中心频率f 0的毫米波信号为,/>的功率谱密度为/>,/>为玻尔兹曼常数,且
其中,θ为球坐标系下位置矢量与z轴正方向的夹角,φ为球坐标系下位置矢量在xoy平面投影矢量与x轴正方向的夹角。
7.如权利要求5所述的三维天线阵综合孔径辐射计快速亮温反演方法,其特征在于,当对谱域有足够的采样点后,通过逆傅里叶变换IFFT重构空间亮温分布,略去常数因子和,三维视函数表示为/>。
8.如权利要求6所述的三维天线阵综合孔径辐射计快速亮温反演方法,其特征在于,为消去γ,进行加权平均运算,两边同时乘以投影函数F(bz),并对基线bz积分,则有二维视函数/>
。
9.如权利要求1所述的三维天线阵综合孔径辐射计快速亮温反演方法,其特征在于,为了反演得到中心频率f 0附近的亮温并在不降低分辨率的前提下消除函数/>的影响,构造如下的函数,/>则,/>,其中,δ<<1。
10.如权利要求1所述的三维天线阵综合孔径辐射计快速亮温反演方法,其特征在于,在三维稀疏阵的情况下,三维视函数在/>非零,二维视函数,将二维视函数/>进行逆傅里叶变换,得到中心频率f 0附近的亮温/>。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010243267.9A CN111538001B (zh) | 2020-03-31 | 2020-03-31 | 三维天线阵综合孔径辐射计快速亮温反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010243267.9A CN111538001B (zh) | 2020-03-31 | 2020-03-31 | 三维天线阵综合孔径辐射计快速亮温反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111538001A CN111538001A (zh) | 2020-08-14 |
CN111538001B true CN111538001B (zh) | 2023-07-21 |
Family
ID=71974829
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010243267.9A Active CN111538001B (zh) | 2020-03-31 | 2020-03-31 | 三维天线阵综合孔径辐射计快速亮温反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111538001B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113533408A (zh) * | 2021-07-21 | 2021-10-22 | 杭州电子科技大学 | 一种改善并行磁共振重建图像质量的变密度数据采样方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102914774A (zh) * | 2012-09-18 | 2013-02-06 | 华中科技大学 | 一种综合孔径辐射计图像反演方法 |
CN103323845A (zh) * | 2013-05-13 | 2013-09-25 | 华中科技大学 | 一种非均匀采样综合孔径辐射计的图像反演方法 |
CN103955602A (zh) * | 2014-04-09 | 2014-07-30 | 华中科技大学 | 一种综合孔径微波辐射计阵列因子成型方法 |
CN105738851A (zh) * | 2016-02-02 | 2016-07-06 | 华中科技大学 | 一种综合孔径微波辐射计联合校正方法 |
CN106842203A (zh) * | 2017-03-07 | 2017-06-13 | 湖北工业大学 | 一种三维天线阵列综合孔径辐射计的图像反演方法 |
CN107167807A (zh) * | 2017-05-23 | 2017-09-15 | 湖北工业大学 | 三维天线阵列综合孔径辐射计分段图像反演方法 |
CN107607946A (zh) * | 2017-07-24 | 2018-01-19 | 湖北工业大学 | 三维均匀采样综合孔径辐射计亮温反演方法 |
CN107797110A (zh) * | 2017-09-13 | 2018-03-13 | 南京邮电大学 | 综合孔径投影辐射的系统成像方法及综合孔径投影辐射计 |
-
2020
- 2020-03-31 CN CN202010243267.9A patent/CN111538001B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102914774A (zh) * | 2012-09-18 | 2013-02-06 | 华中科技大学 | 一种综合孔径辐射计图像反演方法 |
CN103323845A (zh) * | 2013-05-13 | 2013-09-25 | 华中科技大学 | 一种非均匀采样综合孔径辐射计的图像反演方法 |
CN103955602A (zh) * | 2014-04-09 | 2014-07-30 | 华中科技大学 | 一种综合孔径微波辐射计阵列因子成型方法 |
CN105738851A (zh) * | 2016-02-02 | 2016-07-06 | 华中科技大学 | 一种综合孔径微波辐射计联合校正方法 |
CN106842203A (zh) * | 2017-03-07 | 2017-06-13 | 湖北工业大学 | 一种三维天线阵列综合孔径辐射计的图像反演方法 |
CN107167807A (zh) * | 2017-05-23 | 2017-09-15 | 湖北工业大学 | 三维天线阵列综合孔径辐射计分段图像反演方法 |
CN107607946A (zh) * | 2017-07-24 | 2018-01-19 | 湖北工业大学 | 三维均匀采样综合孔径辐射计亮温反演方法 |
CN107797110A (zh) * | 2017-09-13 | 2018-03-13 | 南京邮电大学 | 综合孔径投影辐射的系统成像方法及综合孔径投影辐射计 |
Non-Patent Citations (1)
Title |
---|
黄成丽. 毫米波综合孔径辐射计图像反演研究.《中国优秀硕士学位论文全文数据库 (信息科技辑)》.2016,I140-417. * |
Also Published As
Publication number | Publication date |
---|---|
CN111538001A (zh) | 2020-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6911861B2 (ja) | 物体検知装置および物体検知方法 | |
CN111538000B (zh) | 均匀圆阵列综合孔径辐射计亮温反演成像方法 | |
CN110794471B (zh) | 一种毫米波稀疏阵列远程监视成像方法及系统 | |
Sakamoto et al. | Fast imaging method for security systems using ultrawideband radar | |
CN109884627B (zh) | 任意线阵构型的近程毫米波快速三维成像方法 | |
Yang et al. | A comparative study of compressed sensing approaches for 3-D synthetic aperture radar image reconstruction | |
Molaei et al. | Fourier compatible near-field multiple-input multiple-output terahertz imaging with sparse non-uniform apertures | |
CN111538001B (zh) | 三维天线阵综合孔径辐射计快速亮温反演方法 | |
CN107797110B (zh) | 综合孔径投影辐射的系统成像方法及综合孔径投影辐射计 | |
Wang et al. | Efficient ADMM framework based on functional measurement model for mmW 3-D SAR imaging | |
US7402794B2 (en) | Radiometer imaging system and method thereof | |
CN113933834B (zh) | 圆柱扫描微波成像方法 | |
Chen et al. | An efficient mmW frequency-domain imaging algorithm for near-field scanning 1-D SIMO/MIMO array | |
CN108375767B (zh) | 基于反射板组合的镜像综合孔径辐射计成像方法 | |
Wang et al. | Azimuth improved radar imaging with virtual array in the forward-looking sight | |
CN109581526A (zh) | 一种应用于摄像式毫米波人体安检仪的实时跟踪快速成像方法 | |
CN110231625B (zh) | 一种基于多尺度融合的综合孔径成像方法 | |
CN108375758B (zh) | 基于外部单点源的镜像综合孔径辐射计误差校正方法 | |
CN108375768B (zh) | 基于极化信息的镜像综合孔径辐射计成像方法 | |
CN108445458A (zh) | 一种合成孔径雷达轨道误差消除方法和装置 | |
RU2316786C1 (ru) | Способ наблюдения за поверхностью на базе многоканальной бортовой рлс | |
CN114185047A (zh) | 一种基于最优极坐标变换的双基sar动目标重聚焦方法 | |
Salarkaleji et al. | Metamaterial-based microwave tomography and remote sensing using linear sampling method | |
CN107607946B (zh) | 三维均匀采样综合孔径辐射计亮温反演方法 | |
Almutiry | Wideband tomographic super-resolution radar image |
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 |