頻域稀疏采樣和激光成像方法*

崔岸婧 李道京 吳疆 周凱 高敬涵

1) (中國科學(xué)院空天信息創(chuàng)新研究院,微波成像技術(shù)國家級重點實驗室,北京 100190)

2) (中國科學(xué)院大學(xué),北京 100049)

激光的單色性和自然圖像頻譜稀疏且集中在低頻區(qū)間的特點,使圖像頻譜稀疏采樣成像成為可能.基于小規(guī)模激光探測器,引入?yún)⒖技す?本文提出了頻域稀疏采樣激光成像方法.介紹了頻域稀疏采樣激光成像的原理和成像系統(tǒng)結(jié)構(gòu),推導(dǎo)了激光回波重構(gòu)復(fù)頻譜的表達(dá)式,給出了重構(gòu)頻譜和復(fù)圖像的仿真結(jié)果并分析了信號參數(shù)對重構(gòu)效果的影響,同時采用相干系數(shù)、均方誤差和結(jié)構(gòu)相似度來評價其重構(gòu)效果.規(guī)模為256 × 256的激光回波復(fù)圖像仿真表明,5 個拼接1/4 × 1/4 規(guī)模頻域探測器組成的近似十字型稀疏采樣結(jié)構(gòu),在約31.25% (5/16)的頻域稀疏采樣條件下,仍可獲得較好的重構(gòu)頻譜和重構(gòu)復(fù)圖像.

1 引言

目前大規(guī)模可見光探測器相機(jī)已達(dá)到30k ×30k 量級,高分辨率寬幅成像的需求在推動探測器規(guī)模不斷擴(kuò)大的同時,也產(chǎn)生了數(shù)據(jù)量大、傳輸和處理困難等問題.

前端先大數(shù)據(jù)量采集,然后后端再數(shù)據(jù)壓縮,這似乎是一個前后矛盾的過程.激光具有窄帶性和單色性,由傅里葉光學(xué)成像4f實驗[1,2]和文獻(xiàn)[3,4]可知,利用激光探測獲取的目標(biāo)信號在頻域稀疏且頻譜集中在低頻區(qū)間.在此基礎(chǔ)上,利用傅里葉透鏡將激光圖像變換到頻域,在其低頻區(qū)間實施稀疏采樣,進(jìn)一步通過反演實現(xiàn)激光成像具有可行性.

將小規(guī)模面陣探測器設(shè)置在激光回波頻譜低頻區(qū)間,用頻域稀疏采樣的方式可等效對激光回波復(fù)圖像的二維低頻濾波,在僅丟失圖像部分高頻信息的情況下,以小幅犧牲分辨率為代價減少圖像數(shù)據(jù)量,或可大幅緩解探測器規(guī)模和高分寬幅成像的矛盾.近年快速發(fā)展的計算成像技術(shù)[5]為頻域稀疏采樣激光成像思路提供了一定程度的理論和實踐支持,如文獻(xiàn)[6,7]將傅里葉疊層成像技術(shù)用于遠(yuǎn)距離遙感成像問題.

文獻(xiàn)[8]探討了圖像頻域稀疏激光成像問題,給出了一些初步仿真結(jié)果,基于激光本振相干探測技術(shù)體制,其技術(shù)實現(xiàn)原理清楚.2020 年美國Point Cloud 公司基于硅光芯片的FMCW 激光雷達(dá)相干陣列探測器(陣元規(guī)模為512(32 × 16))[9],使頻域稀疏采樣激光成像成為可能.考慮到目前大量應(yīng)用的激光探測器仍然采用沒有激光本振的直接探測體制,主要實現(xiàn)信號平方律檢波功能,為便于區(qū)分,本文將其簡稱為直接探測器,本文在文獻(xiàn)[8]和時域相干探測[10]的基礎(chǔ)上研究了基于直接探測器的頻域稀疏采樣激光成像問題.

2 頻域稀疏采樣激光成像原理

基于直接探測器的頻域稀疏采樣成像系統(tǒng)結(jié)構(gòu)如圖1 所示,其中F 為經(jīng)接收望遠(yuǎn)鏡聚焦后的激光回波復(fù)圖像,下文將激光回波復(fù)圖像簡稱為激光回波,f表示傅里葉透鏡的焦距.激光器發(fā)射相同的兩束激光,其中一束用于照射目標(biāo)以獲取激光回波,另一束作為參考激光,經(jīng)空間相位調(diào)制后用于激光回波頻譜的重構(gòu).激光回波經(jīng)過傅里葉透鏡處理,其頻譜和參考激光的相干光由直接探測器接收.空間光調(diào)制器采用時分方式,對參考激光實現(xiàn)0°和90°相移處理,以獲取激光回波頻譜的實部和虛部.為減少數(shù)據(jù)量,采用小規(guī)模面陣直接探測器,并在低頻區(qū)由5 個1/4×1/4 規(guī)模頻域探測器構(gòu)成十字型.

圖1 基于直接探測器的頻域稀疏采樣成像系統(tǒng)結(jié)構(gòu)Fig.1.Structure of the frequency-domain sparse sampling imaging system based on direct detector.

在對基于頻域直接探測器的復(fù)圖像重構(gòu)過程進(jìn)行分析前,做出以下假設(shè):

1)面陣探測器中各直接探測像元位于x-y平面上;

2)頻域采樣過程中,激光器、目標(biāo)和傅里葉透鏡之間的距離固定,激光由激光器傳播至傅里葉透鏡所經(jīng)過的光程為定值;

3)在空間光調(diào)制器對參考激光做0°和90°相移過程中,激光回波復(fù)圖像保持不變;

4)參考激光的幅度遠(yuǎn)大于激光回波頻譜的幅度.

激光回波為復(fù)信號,經(jīng)過傅里葉透鏡處理后所得頻譜為復(fù)數(shù),因此可設(shè)置激光回波頻譜的表達(dá)式為

雙路參考激光的表達(dá)式為

參考激光經(jīng)空間光調(diào)制器相移90°后的表達(dá)式為

式中,a1(x,y)和a2(x,y) 分別為回波頻譜和參考激光在探測器平面上的幅度,θ1(x,y,t)=2πf0t+φ1(x,y)和θ2(x,y,t)=2πf0t+φ2(x,y) 分別表示回波頻譜和參考激光在探測器平面上的相位,f0為激 光 載 頻,t為 快 時 間,φ1(x,y)和φ2(x,y) 分 別 為回波頻譜和參考激光在探測器平面上坐標(biāo) (x,y) 位置上的初始相位.

文獻(xiàn)[10]表明,同源同頻(同波長)參考激光和激光回波具有時域相干性,傅里葉透鏡對每個快時間時刻的激光回波復(fù)圖像在二維空間域做傅里葉變換,并未改變激光回波的時域相干性,因此激光回波頻譜與參考激光在時域也具有相干性.

激光回波頻譜和參考激光在探測器平面上相干后的光強(qiáng)[11,12]為

激光回波頻譜和相移90°后的參考激光在探測器平面上相干后的光強(qiáng)為

式中,*表示信號的共軛,θ1(x,y,t)-θ2(x,y,t)=φ1(x,y)-φ2(x,y)消除了激光載波的影響.

由假設(shè)3)可得a2(x,y)?a1(x,y),因此(4)式和(5)式中的(x,y) 可忽略.用功率計可獲取參考激光的光強(qiáng)(x,y),經(jīng)換算可將(4)式和(5)式轉(zhuǎn)化為激光回波頻譜實部和虛部的表達(dá)式

由(6)式和(7)式組合可得復(fù)頻譜

因θ1(x,y,t)-θ2(x,y,t) 僅為x-y平面上的變量,該復(fù)頻譜在任意時刻固定.(8)式所示復(fù)頻譜與激光回波頻譜a1(x,y)ejφ1(x,y) 存在相位差,當(dāng)參考激光的初始相位φ2(x,y)=0 時,激光回波頻譜恢復(fù)效果較好;當(dāng)參考激光的初始相位φ2(x,y)=π/2時,激光回波頻譜的幅度仍能較好恢復(fù),但重構(gòu)頻譜與激光回波頻譜正交.

當(dāng)設(shè)置參考激光的初始相位為0 時,(8)式重構(gòu)頻譜經(jīng)傅里葉逆變換可得激光回波,即實現(xiàn)圖像的重構(gòu).

3 頻譜與圖像重構(gòu)的仿真分析

目前涉及激光雷達(dá)回波復(fù)圖像研究的文獻(xiàn)不多,與微波成像雷達(dá)[13]類似,激光成像雷達(dá)回波信號[14]的初始相位也應(yīng)由目標(biāo)斜距決定,因此激光回波復(fù)圖像的相位與其幅度應(yīng)具有強(qiáng)相關(guān)性,部分傅里葉疊層成像文獻(xiàn)[15,16]對復(fù)圖像的幅度和相位不建立相關(guān)性,會造成復(fù)圖像物理意義不明確.

受實驗條件所限并使問題簡化,本文仿真中假定激光回波復(fù)圖像的相位與其灰度圖相同,傅里葉疊層成像[4,6,7]研究工作中也常做這種假定.

圖2 給出本文仿真分析所用的像素規(guī)模為256 × 256 的復(fù)圖像幅度圖,設(shè)置該復(fù)圖像的相位圖與幅度圖一致,且相位的變化范圍設(shè)置為0— 2π.根據(jù)第2 節(jié)中推導(dǎo)仿真激光回波頻譜和圖像的重構(gòu).本文在仿真中設(shè)置參考激光和激光回波頻譜的幅度關(guān)系為a2(x,y)=30μ[a1(x,y)],a2(x,y) 為不隨x和y變化的固定值,參考激光的初始相位φ2(x,y)=0,其中μ(·) 表示取均值,實際工作中可在沒有參考激光的條件下用頻域直接探測器大致估計a1(x,y) 的均值.

圖2 激光回波復(fù)圖像幅度圖Fig.2.Laser echo complex image amplitude diagram.

3.1 面陣探測器重構(gòu)頻譜和圖像

在面陣中滿布256 × 256 個直接探測器的情況下,仿真激光回波頻譜和圖像的重構(gòu),仿真結(jié)果如圖3 至圖5 所示.由仿真結(jié)果可見,面陣探測器重構(gòu)圖像與激光回波圖像基本一致.

圖3 面陣探測器頻譜實部重構(gòu) (a)頻譜幅度;(b)頻譜相位;(c)圖像Fig.3.Reconstruction of the real part of the plane array detector signals’ spectrum:(a) Amplitude of the spectrum;(b) phase of the spectrum;(c) image.

3.2 稀疏面陣探測器重構(gòu)頻譜和圖像

用5 個1/4×1/4 規(guī)模頻域直接探測器拼接構(gòu)成十字型的情況下,仿真激光回波頻譜和圖像的重構(gòu),探測器范圍如圖6 所示,仿真結(jié)果如圖7 至圖9所示.圖9 與圖5 對比可見,面陣探測器重構(gòu)復(fù)圖像效果優(yōu)于十字型探測器重構(gòu)復(fù)圖像,頻域稀疏導(dǎo)致激光回波頻譜部分信息損失.

圖4 面陣探測器頻譜虛部重構(gòu) (a)頻譜幅度;(b)頻譜相位;(c)圖像Fig.4.Reconstruction of the imaginary part of the plane array detector signals’ spectrum:(a) Amplitude of the spectrum;(b) phase of the spectrum;(c) image.

圖5 面陣探測器頻譜重構(gòu) (a)頻譜幅度;(b)頻譜相位;(c)圖像Fig.5.Reconstruction of the plane array detector signals’ spectrum:(a) Amplitude of the spectrum;(b) phase of the spectrum;(c) image.

圖6 5個1/4 × 1/4 規(guī)模頻域探測器拼接構(gòu)成十字型探測范圍Fig.6.Cross detection range constituted by five 1/4 × 1/4 scale frequency domain detectors.

圖7 稀疏面陣探測器頻譜實部重構(gòu) (a)頻譜幅度;(b)頻譜相位;(c) 圖像Fig.7.Reconstruction of the real part of the sparse plane array detector signals’ spectrum:(a) Amplitude of the spectrum;(b) phase of the spectrum;(c) image.

圖8 稀疏面陣探測器頻譜虛部重構(gòu) (a)頻譜幅度;(b)頻譜相位;(c)圖像Fig.8.Reconstruction of the imaginary part of the sparse plane array detector signals’ spectrum:(a) Amplitude of the spectrum;(b) phase of the spectrum;(c) image.

圖9 稀疏面陣探測器頻譜重構(gòu) (a)頻譜幅度;(b)頻譜相位;(c)圖像Fig.9.Reconstruction of the sparse plane array detector signals’ spectrum:(a) Amplitude of the spectrum;(b) phase of the spectrum;(c) image.

3.3 重構(gòu)效果評價

本文采用相干系數(shù)[17]、均方誤差和結(jié)構(gòu)相似度[18,19]評價激光回波復(fù)圖像和頻譜的重構(gòu)效果.以下將在面陣探測器和稀疏面陣探測器的條件下,根據(jù)(8)式重構(gòu)所得頻譜及其經(jīng)傅里葉逆變換所得復(fù)圖像分別定義為面陣探測器重構(gòu)頻譜、面陣探測器重構(gòu)復(fù)圖像、稀疏面陣探測器重構(gòu)頻譜和稀疏面陣探測器重構(gòu)復(fù)圖像.

相干系數(shù)表示重構(gòu)復(fù)圖像(頻譜)與激光回波復(fù)圖像(頻譜)之間的線性關(guān)系,其范圍為 [0,1],當(dāng)相干系數(shù)接近1 時,復(fù)圖像(頻譜)的重構(gòu)效果越好,其定義式為

式中,E(·) 表示計算數(shù)學(xué)期望,也可用平均值替代,S為激光回波復(fù)圖像或頻譜,為重構(gòu)復(fù)圖像或頻譜.

均方誤差表征重構(gòu)復(fù)圖像(頻譜)與激光回波復(fù)圖像(頻譜)的差異程度,其表達(dá)式為

式中N表示復(fù)圖像(頻譜)的規(guī)模,本文中N=256×256 ;和si分別表示重構(gòu)復(fù)圖像(頻譜)和激光回波復(fù)圖像(頻譜)中第i個單元的值.均方誤差越小,表示重構(gòu)效果越好.

結(jié)構(gòu)相似度從對比度、亮度和結(jié)構(gòu)等角度評價復(fù)圖像(頻譜)的重構(gòu)效果,其表達(dá)式為

式中,D(·)表示計算方差,cov(·) 表示計算協(xié)方差,C1和C2均為常數(shù)矩陣.結(jié)構(gòu)相似度的范圍為 [0,1],當(dāng) SSIM=1 時,兩幅復(fù)圖像(頻譜)完全相同;當(dāng)SSIM=0時,兩幅復(fù)圖像(頻譜)不相關(guān).

面陣探測器和稀疏面陣探測器重構(gòu)頻譜和復(fù)圖像的相干系數(shù)、均方誤差和結(jié)構(gòu)相似度如表1所示.由相干系數(shù)和均方誤差可見稀疏面陣探測器的頻譜重構(gòu)效果較好,但由于利用傅里葉逆變換重構(gòu)復(fù)圖像對頻譜的精確度較高,因此稀疏面陣探測器重構(gòu)復(fù)圖像結(jié)構(gòu)相似度明顯降低,均方誤差增大.

表1 面陣探測器和稀疏面陣探測器頻譜與復(fù)圖像重構(gòu)效果Table 1.Spectrum and complex image reconstruction effect of the plane array detectors and the sparse plane array detectors.

3.4 參數(shù)變化的影響分析

上述仿真和分析均基于參考激光和激光回波頻譜幅度滿足a2(x,y)=30μ[a1(x,y)],參考激光的初始相位φ2(x,y)=0 的條件,以下基于相干系數(shù)分析幅度比值和φ2(x,y) 對重構(gòu)效果的影響.

圖10 為φ2(x,y)=0條件下,不同R對應(yīng)的稀疏面陣探測器重構(gòu)復(fù)圖像的相干系數(shù)曲線.

圖10 R 與稀疏面陣探測器重構(gòu)復(fù)圖像相干系數(shù)的變化曲線Fig.10.Change curves of the relationship between R and correlation coefficient of the reconstructed complex image of the sparse plane array detectors.

由仿真結(jié)果可知,參考激光幅度和激光回波頻譜幅度均值的比值R的增大可提升復(fù)圖像的重構(gòu)性能.當(dāng)R →30時,相干系數(shù)的增大趨于平穩(wěn).R=30條件雖在一定程度上犧牲了探測器的動態(tài)范圍,但能夠保證獲得較好的重構(gòu)結(jié)果.

圖11 為R=30條件 下,不 同φ2(x,y) 對應(yīng)稀疏面陣探測器重構(gòu)復(fù)圖像的相干系數(shù)曲線.結(jié)果表明,參考激光的初始相位對激光回波復(fù)圖像的重構(gòu)效果沒有明顯的影響.

圖11 φ2(x,y) 與稀疏面陣探測器重構(gòu)復(fù)圖像相干系數(shù)的變化曲線Fig.11.Change curves of the relationship betweenφ2(x,y)and correlation coefficient of the reconstructed complex image of the sparse plane array detectors.

4 結(jié)論

本文研究了基于直接探測器的稀疏頻域采樣成像方法并給出了仿真結(jié)果,分析結(jié)果表明該方法具有一定可行性,后續(xù)將開展驗證實驗.

為基于直接探測器重構(gòu)復(fù)頻譜,本文方法采用時分方式對參考激光實現(xiàn)的空間相位調(diào)制,這對高數(shù)據(jù)率成像造成了一定影響,因此本文方法適用于場景變化緩慢、需求對時效性要求不高的情景.與此同時,該方法要求參考激光的幅度較大,這在一定程序上降低了探測器的動態(tài)范圍.

本文在復(fù)頻譜重構(gòu)過程中引入?yún)⒖技す?其思路和數(shù)字全息成像[20,21]有相近之處.目前的數(shù)字全息成像采樣在圖像域完成,其重構(gòu)工作主要是為了恢復(fù)圖像的幅度(灰度),對其相位的重構(gòu)效果較差,將本文方法用于數(shù)字全息成像,不以稀疏采樣減少數(shù)據(jù)量為目的時,有可能獲得比現(xiàn)有數(shù)字全息方法更好的復(fù)圖像重構(gòu)效果,持續(xù)相關(guān)研究工作具有重要意義.