微米顆粒衍射成像過程仿真與特征分析

王文進,方 敬,楊 翔,向 理,田 芃,胡新華

(1. 湖南理工學院 信息光子學與空間光通信湖南省重點實驗室,湖南 岳陽 414006;2. 湖南理工學院 物理與電子科學學院,湖南 岳陽 414006;3. 東卡羅萊納大學 物理系,美國 北卡羅萊納州 27858)

相干散射光強度的空間分布即衍射圖像與散射體的形態和折射率分布高度相關,所以基于散射光成像理論的檢測技術被廣泛應用于生物、大氣、海洋科學等領域[1,2]. 近年來,Yu[3～5]等人開發出一種可以同時獲得散射體在不同偏振方向上的衍射圖的衍射成像式流式細胞儀,如圖1 所示. 入射線偏振光被流體室中的微粒散射后進入數值孔徑為0.55 的顯微物鏡,然后經過鏡筒透鏡成像于CCD 傳感器上. 通過分析衍射圖像特征可以實現包括生物細胞在內的不同結構散射體的區分和結構參數變化分析[6～8]. 但如何利用理論模型對微粒散射光衍射成像過程進行快速、準確分析還沒有很好地解決[9]. 目前研究相干光傳輸成像的主要方法是計算麥克斯韋方程組或求解波動方程,求解復雜且誤差大[10]. 所以必須建立一個快速、精確的理論仿真模型獲得微粒的衍射圖像并對衍射圖像的紋理特征進行提取、分析[11]. 圖像紋理特征提取、分析的方法有很多,其中灰度共生矩陣(GLCM)被證明是一種非常有效的手段[12].

圖1 偏振衍射成像式流式細胞儀

本文基于嚴格的光散射模型獲得了由微球、微圓柱、微橢球組成的1965 個復雜散射體在流體室和顯微物鏡之間的光強分布,并利用高斯混合模型(GMM)對圖像的GLCM參數進行了聚類分析,同時研究了散射體結構參數與GLCM 參數之間的關聯.

1 研究方法

散射體在流體聚焦作用下依次流進流體室. 當散射體經過入射激光的聚焦區域時會產生相干的散射光,散射光經過主介質和玻璃以后進入成像單元,如圖2 所示.

圖2 微粒衍射成像過程仿真

首先,定義散射體三維空間內任一點的折射率

其中nsr為折射率的實部決定散射體的散射系數;nsi為折射率的虛部決定散射體的吸收系數;λ為入射光的波長.

通過這種方式可以實現任意結構、非均勻散射體光學模型的建立. 然后,基于離散偶極子近似理論(ADDA),計算強度為Ii的平行激光入射時空間任意點的光強

其中i1(θ,φ)與i2(θ,φ)分別為與偏振方向平行和垂直散射面的散射強度函數;r為空間任意一點距離散射中心的距離;S ij(θs,φs)為微粒的散射振幅函數即穆勒矩陣,這里i、j分別取值為1、2、3、4;θ、φ為散射光線的散射角和方位角.

通過矩陣元的不同組合運算可以實現任意偏振光散射問題的求解,以入射光偏振方向為45°時S 偏振方向的成像為例,有

最后,利用MATLAB 程序將I(r,θ,φ)投射到入射平面Γin上,Γin位于流體室內是虛擬平面. 入射平面上的光強分布由具有一定尺寸和像素分布的圖像表示,傳播方向由(θs,φs)確定. ZEMAX 基于成像單元模型,利用幾何光學的方法分析散射光經顯微物鏡系統后的光強分布.

圖3 為直徑1～9.5 μm 的微球顆粒的實驗和仿真圖像,通過對比可以看出兩者基本一致. 在文[3]的研究中,通過仿真模型對單球微粒的直徑分布進行了分析,其結果與本實驗結果以及聚乙烯微球供應商提供的數據高度吻合,由此證實了模型的可靠性.

圖3 仿真模型的實驗驗證(A)和圖像特征分析(B)

獲得衍射圖像之后,GLCM 被用于圖像紋理特征的提取和分析. 通過GLCM 可以獲得與圖像相關的16 個參數,每一個參數都具有不同的統計意義,但相互之間也具有一定的相關性. 根據GLCM 參數之間的斯皮爾曼相關系數(大于0.85),可以將16 個參數分為六組,見表1. 結合微粒光散射原理,同一組參數都對應散射體類似的形態特征. 所以,本文分別從六組參數中選取角二階矩(ASM)、熵(ENT)、對比度(CON)、相關性(COR)、協方差(VAR)、逆矩差(IDM)、和熵(SEN)等7 個參數對衍射圖像的紋理特征進行了研究.

表1 GLCM 參數分組結果

其中pi和c分別是散射體和不同分類集合在GLCM 空間內的位置向量和中心.

2 仿真分析

2.1 復雜散射體聚類

散射體形態結構見表2. 表2 中共有1965 個散射體,每個散射體都是由單個或兩個基本單元按照不同排列方式組成. 基本組成單元有球、圓柱、橢球三種不同形狀,每種基本組成單元包括1.558 和1.4 兩種不同折射率. 圖4 為基本組成單元不同的散射體的主要衍射特征. 由于基本組成單元的折射率是均勻的,所以三者的衍射圖像特征類似,都是以條紋狀和網格狀為主. 但是衍射圖像特征之間也有差別,這是由散射體之間結構的差異所引起的.

圖4 不同組成單元散射體的衍射條紋特征

表2 散射體形態結構

利用GLCM 參數可以將衍射圖像紋理反應出來的散射體特征進行量化研究. 表2 中所有散射體的衍射圖像首先歸一化為8 位圖像,然后沿水平、豎直、主對角線、副對角線四個方向提取每一張衍射圖像的GLCM 參數并計算平均值,以用于散射體聚類和結構特征研究. 表2 中的所有散射體、折射率為1.558 的散射體、折射率為1.4 的散射體按照基本組成單元進行無監督聚類,結果見表3.

表3 散射體聚類結果

從表3 中可以看出,當散射體中包括折射率變化時,平均準確率為65.8%. 這主要是因為球和橢球組成的散射體以及圓柱和橢球組成的散射體之間有一部分沒有準確聚類. 當散射體的折射率為1.558 時,聚類準確率可達到95.7%; 當散射體的折射率為1.4 時,聚類準確率為90.5%. 雖然折射率為1.4 的散射體數量不多,其聚類結果不具有統計意義,但是可以看出折射率的變化對聚類結果有明顯的影響. 這說明當散射體折射率均勻時,利用GLCM 參數可以有效分析散射體幾何結構變化對衍射特征的影響.

2.2 散射體結構參數與GLCM 參數定量關聯

研究表明單微球尺寸變化與衍射圖像特征之間具有定量關聯[11]. 本文以雙圓柱模型為例,研究復雜散射體與其衍射特征之間的關系. 雙圓柱的主要參數包括中心間距、中心連線與z軸方向之間的夾角. 圖5所示為雙圓柱中心連線沿z軸方向,中心間距在3.1～6 μm之間變化時衍射圖像特征參數CON的變化. 圖6 所示為雙圓柱中心間距為4.5 μm,中心連線方向與z軸夾角從5～175°變化時衍射圖像特征參數COR 的變化. 圖中圓點為GLCM 參數的實際數值,虛線為擬合的曲線,并給出擬合曲線的數學表達式及擬合優度判定系數R2.

圖5 圓柱中心間距離與CON 的線性關系

圖6 圓柱中心連線方向與COR 的線性關系

從圖5 可以看出,雙圓柱中心間距(C)與衍射圖像特征參數CON 之間具有定量關聯. 當C從3.1 μm 逐漸增加至6 μm 時,CON 也逐漸增加. 此外,CON 隨C變化的散點圖與趨勢直線的擬合優度R2= 0.99,這說明兩者具有良好的線性關系. 圖6 為雙圓柱中心連線方向與z軸夾角(θ)和衍射圖像特征參數(COR)之間的定量關聯. COR 關于θ= 90°對稱,這是因為衍射模型中入射光沿z方向并且雙圓柱散射體具有良好的軸對稱關系. 另外,圖中四次曲線與COR 隨θ變化的散點圖同樣具有良好的擬合關系.

3 結束語

本文基于離散偶極子近似算法,結合光跡追蹤軟件,對復雜結構散射體的衍射過程進行了快速、精確仿真,得到了不同結構散射體的衍射圖像. 并利用衍射圖像量化特征實現了散射體的無監督聚類分析,建立衍射圖像紋理特征參數與散射體三維結構參數之間的定量關聯. 研究表明,衍射圖像特征雖然很難用視覺來量化和解釋,但是非常適用于微米量級散射體的免標記分類與識別以及結構參數的定量分析.