CN101788665A - 一种被动毫米波焦平面成像装置 - Google Patents
一种被动毫米波焦平面成像装置 Download PDFInfo
- Publication number
- CN101788665A CN101788665A CN 201010124342 CN201010124342A CN101788665A CN 101788665 A CN101788665 A CN 101788665A CN 201010124342 CN201010124342 CN 201010124342 CN 201010124342 A CN201010124342 A CN 201010124342A CN 101788665 A CN101788665 A CN 101788665A
- Authority
- CN
- China
- Prior art keywords
- cij
- feed
- millimeter wave
- focal plane
- image formation
- 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
Links
Images
Landscapes
- Radiation Pyrometers (AREA)
Abstract
一种被动毫米波焦平面成像装置,属于被动毫米波成像领域,解决现有成像系统存在的像素数有限、成像分辨率较低的问题。本发明包括N个馈源和N路毫米波辐射计通道,N个馈源以线阵形式排列在聚焦天线的焦平面上,每路毫米波辐射计通道由射频放大器、射频滤波器、混频器和模数转换器依次串联构成,N路中频数字信号送入计算机进行自相关和互相关处理,计算得到天线温度TAi和干涉测量亮温Tcij,计算机再对天线温度TAi处理得到像素值Ti,将干涉测量亮温Tcij插入像素值Ti和Tj之间,共同构成焦平面图像。本发明将相干成像和非相干成像进行有效结合,增加了图像的像素,进而提高了成像的空间分辨率。
Description
技术领域
本发明属于被动毫米波成像领域,具体涉及一种被动毫米波焦平面成像装置。
背景技术
被动毫米波(Passive Millimeter-wave,PMMW)成像系统具有全天时、全天候的工作能力,并且因为不发射电磁波,具有很高的隐蔽性和抗干扰性,对人体安全,在反恐斗争、场景监控、安全检查、军事侦察等领域具有广阔的应用前景。早期的被动毫米波成像系统由于毫米波器件不成熟,都是采用单通道机械扫描结构,随着毫米波器件的不断进步以及被动毫米波技术应用范围的日益扩大,这类系统的一个自然发展是采用多通道并行扫描以提高数据率及成像速度,最终的目标是构建大型焦平面阵列凝视场景成像,即被动毫米波焦平面成像系统。
被动毫米波焦平面成像与红外焦平面成像原理上类似,属于非相干的直接成像,其原理是将多个接收单元以一定的阵列排列在聚焦天线的焦平面上,利用各馈源的偏焦不同产生不同指向的波束覆盖视域,每一波束对应场景中的一个像点,相应接收单元的输出信号即对应着该像点的亮度温度。
目前,被动毫米波焦平面成像系统在国外已走向样机阶段,如美国TRW空间与电子设备集团已研制成功W波段焦平面PMMW实时成像仪,其系统频率为89GHz,空间分辨率为0.5度,视场角为15°×10°,接收机阵列数量达到1040,图像像素数为26×40,接收机采用单片集成MMIC的直接检波式接收机;Lockheed Martin公司的PMMW成像仪样机接收天线采用卡塞格伦天线,水平方向线形排列34个接收机,垂直方向机械扫描,角分辨率为0.5°,最终可获得34×44个像素,帧频为10Hz;Millivision手持式PMMW扫描成像仪采用可旋转镜片实现扫描成像,频率94GHz,视域26°×26°,空间分辨率5cm(1.6m远处),像素数为26×40,帧频为10Hz。这些样机都采用基于非相干测量的常规被动毫米波焦平面成像,即一个馈源输出对应一个像素,所形成的像素数有限,成像的分辨率还比较低。
发明内容
本发明提供一种被动毫米波焦平面成像装置,解决现有被动毫米波焦平面成像系统存在的像素数有限、成像分辨率较低的问题,将相干成像和非相干成像进行有效结合,获得比基于非相干测量技术的常规焦平面成像方法更多的像素,进而提高成像的分辨率。
本发明的一种被动毫米波焦平面成像装置,包括N个馈源和N路毫米波辐射计通道,N个馈源以线阵形式排列在聚焦天线的焦平面上,每路毫米波辐射计通道由射频放大器、射频滤波器、混频器和模数转换器依次串联构成,N路毫米波辐射计通道输出的N路中频数字信号Xi(n)送入计算机进行自相关处理,得到自相关输出值Vi:Q为每路中频数字信号Xi(n)的采样个数,i=1~N,N=2~2000,其特征在于:
计算机对所述N路中频数字信号Xi(n)中每相邻两路进行互相关处理,得到互相关输出值Гij:
然后,计算机对自相关输出值Vi和互相关输出值Γij进行处理,得到天线温度TAi和干涉测量亮温Tcij:
最后,计算机根据干涉测量亮温Tcij对天线温度TAi进行处理,得到像素值Ti:
其中干涉测量亮温Tcij=Tcji;Ωi为第i号馈源波束的立体角,Ωcij为第i号和j号馈源波束重叠部分的立体角,Ωcij=Ωcji;
将干涉测量亮温Tcij插入像素值Ti和Tj之间,以N个像素值Ti和N-1个干涉测量亮温Tcij共同构成焦平面图像。
所述的被动毫米波焦平面成像装置,所述每路毫米波辐射计通道中,混频器和模数转换器之间,可以依次串联中频滤波器和中频放大器。
如图1所示,馈源5排列在焦平面4上,其中相邻排列的馈源i和馈源j通过透镜3聚焦,存在波束交叉区域2,馈源i和馈源j的波束分别指向物平面1上的AC区以及BC区,则物平面1上C区的辐射会同时进入两个馈源,两馈源的输出进行相关后得到的结果将只包括C区的亮温信息。由于馈源波束足够窄,可将相邻波束重叠部分对应的场景亮温当成一定值。
根据辐射测量原理,自相关输出值Vi与馈源接收到的来自场景的毫米波辐射功率成正比,即Vi=GPi,G为通道增益(射频放大器增益或者射频放大器和中频放大器增益之和),功率Pi可表示为Pi=kTAiΔf,其中k为波尔兹曼常量,Δf为系统带宽,TAi为天线温度,它是视在亮温分布按天线功率方向图加权得到的,可见自相关输出值Vi与天线温度TAi成正比,根据自相关输出值Vi可以得到天线温度TAi
由于不同物点的毫米波辐射互不相关,对相邻两个馈源的输出信号进行互相关处理,互相关输出值Гij只包含来自两馈源波束重叠或交叉部分的辐射亮温信息。
根据干涉测量理论,对来自相邻馈源的信号进行互相关处理,其结果可表示为:
其中Tcij为待测的波束重叠部分的亮温,Tcij=Tcji;fni(θ,φ)与fnj(θ,φ)分别为馈源i和馈源j的归一化场强方向图;Di与Dj分别为馈源i和馈源j的最大方向性;ΔS(θ,φ)为同一物点的毫米波辐射进入两馈源的路程差;λ为接收机工作频率所对应的波长。
上式的积分项为一定值,设为参数S,即Гij=STcij,则互相关输出和待测的波束重叠部分的亮温成线性关系,即可由互相关输出Γij得到干涉测量亮温Tcij。
本发明的工作过程为:
(1)N个馈源以线阵形式排列在聚焦天线的焦平面上,利用各馈源的偏焦不同产生不同指向的波束覆盖视域,接收场景的毫米波辐射;第i个馈源接收毫米波辐射输出的信号为Ui(t),
其中i=1,2,...,N,N为大于1的整数,t为时间变量;
(2)将每路馈源输出的信号分别经放大滤波以及混频处理后,得到中频信号Yi(t),i=1,2,...,N,对中频模拟信号Yi(t)进行采样得到中频数字信号Xi(n),i=1,2,...,N,输入计算机进行存储;
(3)计算机将中频数字信号Xi(n)进行自相关处理,得到一路自相关输出值Vi,
同时,计算机对每两路来自相邻馈源的中频数字信号进行互相关处理得到一组互相关输出Гij,其中i、j表示两相邻馈源的编号,Гij表示来自第i号和j号馈源的信号相关的结果,Гij=Гji;互相关输出Гij与波束重叠部分的亮温成线性关系,
(4)自相关输出值Vi通过定标得到一路天线温度TAi,i=1,2,...,N;以天线温度作为一个像素的度量值,共有N个像素度量值;
互相关输出Гij通过定标得到干涉测量亮温Tcij;
(5)用干涉测量亮温Tcij,修正步骤(4)得到的N个像素度量值,作为像素值Ti:
其中干涉测量亮温Tcij=Tcji;Ωi为第i号馈源波束的立体角,Ωcij为第i号和j号馈源波束重叠部分的立体角,Ωcij=Ωcji;
(6)以干涉测量亮温Tcij作为新生成像素值,将其插入到第i号及第j号馈源对应的像素值之间;
以步骤(5)得到的N个像素值和步骤(4)的N-1个新生成像素值共同构成焦平面图像。
若需扩大成像视场范围,则以0.5倍半功率波束宽度为步长转动波束指向进行扫描,扫描的方向根据馈源排列的形式确定。
本发明利用排列在聚焦天线焦平面上的馈源接收场景的被动毫米波辐射,通过对相邻馈源波束存在的重叠区域,采用互相关处理产生新的像素信息,将其融合到传统的焦平面图像中,增加了图像的像素;并且利用获得的相邻波束重叠部分的亮温值对常规焦平面成像中的天线温度进行修正,修正后的天线温度对应的场景测量区域小于原来非相干测量的区域,提高了成像的空间分辨率。
附图说明
图1为被动毫米波焦平面成像示意图;
图2为本发明实施例的结构示意图;
图3为本发明实施例中的计算机处理过程示意图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
如图2所示,本发明包括2个馈源和2路毫米波辐射计通道,2个馈源以线阵形式排列在聚焦天线的焦平面上,在垂直于阵列排列方向上结合扫描进行成像,每路毫米波辐射计通道由射频放大器、射频滤波器、混频器和模数转换器依次串联构成,2路毫米波辐射计通道输出的2路中频数字信号送入计算机进行自相关和互相关处理,最终得到毫米波图像。
所述每路毫米波辐射计通道中,混频器和模数转换器之间,可以依次串联中频滤波器和中频放大器。
聚焦天线反射面采用抛物柱面,馈源选用矩形喇叭,其半功率波束宽度为18度,两馈源间隔5.2cm;
毫米波辐射计通道中,射频放大器中心频率为36.4GHz、增益为55dB,射频滤波器带宽>100MHz,混频器中的本振频率为36.52GHz,中频滤波器带宽>100MHz,中频放大器增益为47dB;模数转换器采用双通道8bit高速AD卡,采样频率为250MHz,最大存储深度为8Mbytes。
图3为本发明实施例中的计算机处理过程示意图。
(1)聚焦天线反射面采用抛物柱面,馈源选用矩形喇叭,其半功率波束宽度为18度,两馈源间隔5.2cm,接收场景的毫米波辐射;馈源接收毫米波辐射输出的信号分别为U1(t),U2(t);
(2)高频前端模块中的射频放大器的中心频率为36.4GHz,增益为55dB,系统带宽100MHz,混频器中的本振频率为36.52GHz,中频放大器增益为47dB,高频前端的净增益为102dB,则信号U1(t),U2(t)分别经高频前端后,得到中频信号Y1(t),Y2(t),其中心频率为120MHz;
(3)用8bit高速AD卡对中频信号Y1(t),Y2(t)进行采样存储,AD卡的采样频率为250MHz,最大存储深度为8Mbytes,采样后得到的数字信号为X1(n),X2(n),n=1,2,...,Q;
(4)分别对数字信号X1(n),X2(n)平方(相当于平方律检波器),然后对Q个平方的数据求平均得到输出值V1,V2,
又V1=GkTA1Δf,V2=GkTA2Δf,G为通道增益,在这里为102dB,k为波尔兹曼常量,Δf为系统带宽,在这里为100MHz,由此可以得到天线温度TA1,TA2,以此2个天线温度作为2个像素的度量值。
(5)对数字信号X1(n),X2(n)进行互相关处理得到输出Г12,其中Г12表示来自第1号和第2号馈源的信号相关的结果,Гc21=Γc12;
又Г12=STc12,在这里参数S=1.34,由此通过互相关输出Г12得到干涉测量亮温Tc12;
(6)用干涉测量亮温Tc12,修正步骤(4)得到的像素度量值,作为像素值T1,T2:
其中干涉测量亮温Tc12=Tc21;Ω1为第1号馈源波束的立体角,Ω2为第2号馈源波束的立体角,均为18度;Ωc12为第1号和第2号馈源波束重叠部分的立体角,为9度;Ωc12=Ωc21;
(7)用像素值T1,T2和干涉测量亮温值Tc12构成新的毫米波图像;
若需扩大成像视场范围,则以0.5倍半功率波束宽度为步长转动波束指向进行扫描,扫描的方向根据馈源排列的形式确定。
Claims (2)
1.一种被动毫米波焦平面成像装置,包括N个馈源和N路毫米波辐射计通道,N个馈源以线阵形式排列在聚焦天线的焦平面上,每路毫米波辐射计通道由射频放大器、射频滤波器、混频器和模数转换器依次串联构成,N路毫米波辐射计通道输出的N路中频数字信号Xi(n)送入计算机进行自相关处理,得到自相关输出值Vi:Q为每路中频数字信号Xi(n)的采样个数,i=1~N,N=2~2000,其特征在于:
计算机对所述N路中频数字信号Xi(n)中每相邻两路进行互相关处理,得到互相关输出值Гij:
然后,计算机对自相关输出值Vi和互相关输出值Гij进行处理,得到天线温度TAi和干涉测量亮温Tcij:
最后,计算机根据干涉测量亮温Tcij对天线温度TAi进行处理,得到像素值Ti:
其中干涉测量亮温Tcij=Tcji;Ωi为第i号馈源波束的立体角,Ωcij为第i号和j号馈源波束重叠部分的立体角,Ωcij=Ωcji;
将干涉测量亮温Tcij插入像素值Ti和Tj之间,以N个像素值Ti和N-1个干涉测量亮温Tcij共同构成焦平面图像。
2.如权利要求1所述的被动毫米波焦平面成像装置,其特征在于:
所述每路毫米波辐射计通道中,混频器和模数转换器之间,依次串联中频滤波器和中频放大器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101243426A CN101788665B (zh) | 2010-03-16 | 2010-03-16 | 一种被动毫米波焦平面成像装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101243426A CN101788665B (zh) | 2010-03-16 | 2010-03-16 | 一种被动毫米波焦平面成像装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101788665A true CN101788665A (zh) | 2010-07-28 |
CN101788665B CN101788665B (zh) | 2011-12-14 |
Family
ID=42531946
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101243426A Active CN101788665B (zh) | 2010-03-16 | 2010-03-16 | 一种被动毫米波焦平面成像装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101788665B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102003997A (zh) * | 2010-09-09 | 2011-04-06 | 北京航空航天大学 | 光学滤波的毫米波频率扫描成像仪 |
CN102288959A (zh) * | 2011-06-20 | 2011-12-21 | 北京理工大学 | 一种毫米波焦平面被动成像系统及方法 |
CN104698461A (zh) * | 2015-03-18 | 2015-06-10 | 中国科学院空间科学与应用研究中心 | 一种毫米波被动安检成像仪及人体安检成像系统 |
CN106154345A (zh) * | 2015-03-27 | 2016-11-23 | 常州吉赫射频电子技术有限公司 | 椭球面被动毫米波成像系统 |
CN106533494A (zh) * | 2016-12-22 | 2017-03-22 | 华讯方舟科技(湖北)有限公司 | 一种微波变频器及其微波变频电路 |
CN107037495A (zh) * | 2017-03-17 | 2017-08-11 | 芜湖博高光电科技股份有限公司 | 一种超宽带双模毫米波主被动复合探测器 |
CN108761553A (zh) * | 2018-03-16 | 2018-11-06 | 电子科技大学 | 用于安检的被动毫米波双通道同步成像系统及其成像方法 |
CN109029367A (zh) * | 2018-08-30 | 2018-12-18 | 中国科学院长春光学精密机械与物理研究所 | 一种基于拓展目标点的凝视成像方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644770A (zh) * | 2009-09-07 | 2010-02-10 | 哈尔滨工业大学 | 被动式毫米波成像系统 |
-
2010
- 2010-03-16 CN CN2010101243426A patent/CN101788665B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644770A (zh) * | 2009-09-07 | 2010-02-10 | 哈尔滨工业大学 | 被动式毫米波成像系统 |
Non-Patent Citations (6)
Title |
---|
《2007年全国微波毫米波会议论文集》 20071231 张光锋 等 3mm波段被动毫米波成像研究 1354-1357 1-2 , 2 * |
《制导与引信》 20031231 钱嵩松 等 被动毫米波成像综述 29-36 1-2 第24卷, 第4期 2 * |
《探测与控制学报》 20040930 钱嵩松 等 用于隐匿武器检测的PMMW焦平面成像系统设计 58-61 1-2 第26卷, 第3期 2 * |
《电视技术》 20081031 万金 等 3毫米十六通道焦平面辐射计数据采集与测试系统设计 1-4 1-2 第48卷, 第10期 2 * |
《系统仿真学报》 20071231 熊祖彪 等 综合孔径微波辐射探测系统仿真平台设计 5663-5666 1-2 第19卷, 第24期 2 * |
《系统仿真学报》 20090430 胡飞 等 微波辐射探测仿真系统并行处理研究 2057-2060 1-2 第21卷, 第7期 2 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102003997A (zh) * | 2010-09-09 | 2011-04-06 | 北京航空航天大学 | 光学滤波的毫米波频率扫描成像仪 |
CN102003997B (zh) * | 2010-09-09 | 2011-10-26 | 北京航空航天大学 | 光学滤波的毫米波频率扫描成像仪 |
CN102288959A (zh) * | 2011-06-20 | 2011-12-21 | 北京理工大学 | 一种毫米波焦平面被动成像系统及方法 |
CN104698461A (zh) * | 2015-03-18 | 2015-06-10 | 中国科学院空间科学与应用研究中心 | 一种毫米波被动安检成像仪及人体安检成像系统 |
CN106154345A (zh) * | 2015-03-27 | 2016-11-23 | 常州吉赫射频电子技术有限公司 | 椭球面被动毫米波成像系统 |
CN106533494A (zh) * | 2016-12-22 | 2017-03-22 | 华讯方舟科技(湖北)有限公司 | 一种微波变频器及其微波变频电路 |
CN106533494B (zh) * | 2016-12-22 | 2019-05-03 | 华讯方舟科技(湖北)有限公司 | 一种微波变频器及其微波变频电路 |
CN107037495A (zh) * | 2017-03-17 | 2017-08-11 | 芜湖博高光电科技股份有限公司 | 一种超宽带双模毫米波主被动复合探测器 |
CN107037495B (zh) * | 2017-03-17 | 2019-04-23 | 芜湖博高光电科技股份有限公司 | 一种超宽带双模毫米波主被动复合探测器 |
CN108761553A (zh) * | 2018-03-16 | 2018-11-06 | 电子科技大学 | 用于安检的被动毫米波双通道同步成像系统及其成像方法 |
CN109029367A (zh) * | 2018-08-30 | 2018-12-18 | 中国科学院长春光学精密机械与物理研究所 | 一种基于拓展目标点的凝视成像方法 |
CN109029367B (zh) * | 2018-08-30 | 2020-07-07 | 中国科学院长春光学精密机械与物理研究所 | 一种基于拓展目标点的凝视成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101788665B (zh) | 2011-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101788665B (zh) | 一种被动毫米波焦平面成像装置 | |
Ruze | Antenna tolerance theory—A review | |
CN102052967B (zh) | 多像元超导探测器接收系统及太赫兹信号检测方法 | |
CN102288959A (zh) | 一种毫米波焦平面被动成像系统及方法 | |
CN106848546A (zh) | 一种宽带双极化天线阵列装置及高分辨测向方法 | |
CN107957574A (zh) | 基于ifft和混合匹配追踪的时分地基mimo滑坡雷达成像方法 | |
CN111458051B (zh) | 基于像素级光谱分光探测器的三维温度场测量系统和方法 | |
US4499470A (en) | Method of measuring sea surface water temperature with a satellite including wideband passive synthetic-aperture multichannel receiver | |
CN103399315A (zh) | 实孔径相控阵雷达高分辨探测成像方法 | |
CN108535725A (zh) | 星载多频段一维综合孔径一维实孔径的微波辐射探测方法 | |
CN110764059B (zh) | 一种收发垂直波束三坐标相控阵雷达方法 | |
CN104122273A (zh) | 基于多通道频带合成的辐射计 | |
US7402794B2 (en) | Radiometer imaging system and method thereof | |
CN113253305A (zh) | 阵列天线获取卫星入射信号导向矢量的方法 | |
CN106644099A (zh) | 一种多光谱温度场测量装置 | |
CN106450703A (zh) | 一种基于循环子阵的干涉式微波辐射计圆环天线阵列 | |
CN107466063A (zh) | 一种通信卫星多波束无线测试方法 | |
JP2004336390A (ja) | アダプティブアレーおよび測位装置 | |
RU2008112092A (ru) | Способ повышения разрешающей способности радиотепловых изображений | |
CN109188381B (zh) | 一种基于外部噪声注入的毫米波干涉式综合孔径成像系统快速校准方法 | |
CN114200448A (zh) | 综合孔径辐射计波数域近场成像方法及设备 | |
Huang et al. | A novel millimeter wave synthetic aperture radiometer passive imaging system | |
Schreiber et al. | VESAS: a novel concept for fully-electronic passive MW imaging | |
Guo et al. | Millimeter-wave MIMO imaging radar with an ultra-high angular resolution of 0.6 degree | |
Tran et al. | UHF/C-band testing of AOA estimation using MUSIC algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |