发明名称 多频连续波系统的聚焦成像校准方法
摘要 本发明提供了一种多频连续波系统聚焦成像校准方法,从多频连续波系统的二维扫描数据中,抽取出各频率数据,先进行单频聚焦成像,获得最佳聚焦图像。在此图像上抽取大信噪比像素数据,计算出平均相位偏差和平均幅度。再利用计算结果构成校准补偿因子,对多频数据进行相位和幅度校准。校准后的数据进行多频成像处理,能够获得良好的聚焦成像效果,目标轮廓更清晰,背景杂波更少,信噪比有效提高。由于所利用的数据多,信息丰富,校准精度高。本发明是数据驱动的校准方法,无需校准硬件附件,也不需要进行复杂的校准测量。
申请公布号 CN104569936B 申请公布日期 2016.09.07
申请号 CN201510037428.8 申请日期 2015.01.26
申请人 中国人民解放军国防科学技术大学 发明人 黄春琳;陆珉;粟毅;李禹;唐涛
分类号 G01S7/40(2006.01)I 主分类号 G01S7/40(2006.01)I
代理机构 北京中济纬天专利代理有限公司 11429 代理人 胡伟华
主权项 一种多频连续波系统的聚焦成像校准方法,其特征在于包括单频聚焦成像、大信噪比数据抽取、平均相位偏差和平均幅度计算、多频数据的相位和幅度校准四个步骤,第一步,单频聚焦成像从多频连续波系统的二维扫描数据中,抽取出各频率的数据,对每个频率的数据单独进行聚焦成像处理,用下式表示:Y(f<sub>i</sub>)=ψ{R(f<sub>i</sub>)}                 (1)式中R(f<sub>i</sub>)为频率f<sub>i</sub>的单频数据矩阵(f<sub>i</sub>=f<sub>1</sub>,f<sub>2</sub>,…,f<sub>I</sub>),来自x轴和y轴平面二维扫描的多频数据矩阵R′<sub>fxy</sub>;ψ{·}表示聚焦成像处理的运算符,在目标距离d和传播介质介电常数ε两参数已知情况下,直接输入参数值获得最佳聚焦成像结果;若参数未知,采用自聚焦方法,获得最佳聚焦度的成像结果;矩阵Y(f<sub>i</sub>)为聚焦成像后的图像数据矩阵;第二步,大信噪比数据抽取从单频聚焦成像结果Y(f<sub>i</sub>)中抽取出大信噪比像素,构成数组D<sub>k</sub>(f<sub>i</sub>):D<sub>k</sub>(f<sub>i</sub>)=Y<sub>mn</sub>(f<sub>i</sub>),if|Y<sub>mn</sub>(f<sub>i</sub>)|≥λ k=1,2,…,K   (2)式中λ为幅度判断门限,由所允许的最低信噪比决定;m为x轴的像素序号,其中m=1,2,…,M;n为y轴的像素序号,其中n=1,2,…,N;第三步,平均相位偏差和平均幅度计算由抽取的数组D<sub>k</sub>(f<sub>i</sub>),通过期望值计算,得到平均相位偏差和平均幅度:<img file="FDA0001044958820000011.GIF" wi="1149" he="142" /><maths num="0001"><math><![CDATA[<mrow><mover><mi>&alpha;</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>K</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mo>|</mo><msub><mi>D</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001044958820000012.GIF" wi="1086" he="143" /></maths>式中ang{·}为相位计算函数;第四步,多频数据的相位和幅度校准利用第三步计算出的平均相位偏差<img file="FDA0001044958820000021.GIF" wi="159" he="69" />和平均幅度<img file="FDA0001044958820000022.GIF" wi="121" he="62" />构成校准补偿因子:<img file="FDA0001044958820000023.GIF" wi="1118" he="70" />再利用校准补偿因子对多频数据矩阵R′<sub>fxy</sub>进行相位和幅度校准:<maths num="0002"><math><![CDATA[<mrow><msubsup><mi>R</mi><mrow><msub><mi>f</mi><mi>i</mi></msub><mi>mn</mi></mrow><mo>&Prime;</mo></msubsup><mo>=</mo><msubsup><mi>R</mi><mrow><msub><mi>f</mi><mi>i</mi></msub><mi>mn</mi></mrow><mo>&prime;</mo></msubsup><mo>&CenterDot;</mo><mi>&gamma;</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001044958820000024.GIF" wi="942" he="79" /></maths>公式(6)是对多频数据矩阵R′<sub>fxy</sub>的每个单元数据做校准补偿,校准后的数据R″<sub>fxy</sub>能使得目标成像结果的相位均值为零,幅度均值为1,以去除各频率的相位偏差,补偿幅度的不一致性;完成校准后,利用校准后的数据R″<sub>fxy</sub>进行多频成像处理,能够获得无相位偏差和无幅度偏差的成像结果。
地址 410073 湖南省长沙市开福区德雅路109号