利用光流場進(jìn)行物體變形相位測量的方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及一種利用光流場進(jìn)行物體變形相位測量的方法。
【背景技術(shù)】
[0002] 光學(xué)干涉測量方法因其獨(dú)特的優(yōu)勢在物體的三維形貌和變形測量中占有十分重 要的地位，如電子散斑干涉、全息干涉、云紋干涉法等。在這些方法中需要解調(diào)出條紋圖 中每一點(diǎn)的相位信息才可以實(shí)現(xiàn)變形測量。解調(diào)相位的方法主要分為時(shí)間相位測量法 和空間相位測量法兩大類。時(shí)間相位測量法，如時(shí)間相移法（Temporal Phase-shifting Method，TPM)等等，具有位移測量精度高，結(jié)構(gòu)簡單等優(yōu)點(diǎn)，但是該類方法需要的條紋圖較 多，耗時(shí)較長，一般情況下只適用于靜態(tài)測量；空間相位測量法是通過在空間上獲取更多信 息來提取相位的，如空間相移法（Spatial Phase-shifting Method, SPM)和Fourier變換 法（Fourier Transform Method, FTM)等等，可實(shí)現(xiàn)對動(dòng)態(tài)過程相位的測量。但是，前者的 光路較為復(fù)雜，且精度較低；后者的處理過程比較繁瑣，很難實(shí)現(xiàn)條紋處理的自動(dòng)化。同時(shí)， 以上方法對于位移場的測量均沒有考慮到時(shí)間參量，因而對動(dòng)態(tài)場的定量測量就無法獲得 滿意的結(jié)果。針對此問題，有研究工作者提出了時(shí)間序列相位法（Time Sequence Phase Method，TSPM)，將時(shí)間參量引入散斑干涉當(dāng)中，成功地提取出了相位信息，能適應(yīng)大變形、 時(shí)變測量。然而該方法在位移方向性上無法給出確定值，故只能用于單調(diào)的時(shí)變場中。除 此之外，以上方法都需要進(jìn)行相位解包絡(luò)操作，這無疑增加了相位提取的運(yùn)算量和誤差。同 時(shí)，由于在條紋稠密區(qū)域無法判斷每個(gè)條紋的細(xì)節(jié)，故以上方法在條紋密集處解調(diào)的相位 誤差都比較大。
[0003] 20世紀(jì)50年代，Gribson等研究人員基于運(yùn)動(dòng)結(jié)構(gòu)重建原理（Structure From Motion, SFM)第一次提出了光流（optical flow)的概念：在隨時(shí)間變化的二維圖像序列 中，運(yùn)動(dòng)的三維場景由于亮度模式的變化而產(chǎn)生的流動(dòng)即為光流。從此以后，國內(nèi)外學(xué)者陸 續(xù)提出了不同的光流方法，在精度，魯棒性和實(shí)時(shí)性等多方面均取得了重大突破。1997年， 國內(nèi)研究人員提出可以將光流場分析法引入光學(xué)干涉計(jì)量當(dāng)中，為光干涉動(dòng)態(tài)測量提供了 新的思路，但是其并未做深入的探討。2011年，J. Vargas等人通過運(yùn)用光流場理論和SPT 算符成功提取出了單幅圖像相位，但是該方法仍需要相位解包等一系列復(fù)雜操作，且不適 用于實(shí)現(xiàn)動(dòng)態(tài)形變測量。

【發(fā)明內(nèi)容】

[0004] 本發(fā)明的目的就是為了解決上述問題，提供一種利用光流場進(jìn)行物體變形相位測 量的方法，只需在空間域中用兩幅條紋圖就可以完成變形相位的提取，無需轉(zhuǎn)換到頻域和 相位解包絡(luò)操作，適用于動(dòng)態(tài)測量，且在條紋越密集的區(qū)域提取相位的效果越好，完全克服 了現(xiàn)有相位解調(diào)方法中條紋過密時(shí)解調(diào)誤差較大的缺點(diǎn)。
[0005] 為了實(shí)現(xiàn)上述目的，本發(fā)明采用如下技術(shù)方案：
[0006] 利用光流場進(jìn)行物體變形相位測量的方法，包括以下步驟：
[0007] 步驟一：根據(jù)光流場理論的基本假設(shè)，獲得兩幀連續(xù)圖像之間的運(yùn)動(dòng)矢量的兩個(gè) 分量u和V ;
[0008] 步驟二：由窗口傅里葉變換法提取第一幀圖像條紋的橫向頻率fx。和縱向頻率 fyO；
[0009] 步驟三：利用所述步驟一得到的運(yùn)動(dòng)矢量的兩個(gè)分量U和V及所述步驟二得到條 紋的橫向頻率fx(:和縱向頻率f y(:計(jì)算得到相位的改變量。
[0010] 所述步驟一的具體方法為：運(yùn)用公式：
[0013] 其中，
工^仁和工沒別為圖像I對于x，y，t的偏導(dǎo)數(shù)，I為t時(shí) 刻像素點(diǎn)（x，y)處灰度值，un+1和vn+1為在圖像的每點(diǎn)第n+1次迭代后所得到的運(yùn)動(dòng)矢量的 兩個(gè)分量u和V的值，即兩幅圖像之間的運(yùn)動(dòng)矢量場；ulP V n為在圖像的每點(diǎn)第η次迭代 后所得到的運(yùn)動(dòng)矢量的兩個(gè)分量u和V的值；α為平滑參數(shù)。
[0014] 所述步驟二的結(jié)果由公式 求得，（fxD，fyD)為使|si(x，y，ξ，η)|取最大值的U，n)，其中，
為空間一點(diǎn)（x，y)對應(yīng)的頻域坐 標(biāo)，SI ( μ，V ; ξ，η )是I (X，y)經(jīng)過窗口傅里葉變換后對應(yīng)頻域的結(jié)果，I (X，y)為一幅輸 入的二維圖像；g(x, y)為窗口函數(shù)，取高斯函數(shù)
odP σ 高斯分布的標(biāo)準(zhǔn)差，其取值決定著高斯窗口的大小。
[0015] 所述步驟三的具體方法為：將運(yùn)動(dòng)矢量的兩個(gè)分量u和V，條紋的橫向頻率fx。和 縱向頻率4。帶入公式
，為點(diǎn)（X。，y。）在Δ t時(shí)間內(nèi) 相位的改變量，求得兩幀連續(xù)圖像的相位該變量。
[0016] 本發(fā)明的有益效果：
[0017] 本發(fā)明只需在空間域中用兩幅條紋圖就可以完成變形相位的提取，無需轉(zhuǎn)換到頻 域和相位解包絡(luò)操作，適用于動(dòng)態(tài)測量，且在條紋越密集的區(qū)域提取相位的效果越好，完全 克服了現(xiàn)有相位解調(diào)方法中條紋過密時(shí)解調(diào)誤差較大的缺點(diǎn)，過程簡單、方便，同時(shí)該方法 也為光干涉動(dòng)態(tài)測量提供了新的途徑。
【附圖說明】
[0018] 圖1為原始干涉條紋模擬圖；
[0019] 圖2 (a)為附加相位為麵>i = Irad后的條紋模擬圖，圖2 (b)為附加相位 Δ灼=0.01 w,,后的條紋模擬圖；
[0020] 圖3為附加相位為Δ約時(shí)所估算的在y = 127處X方向光流場橫截面圖；
[0021 ] 圖4為對附加相位麵\提取的歸一化模擬結(jié)果；
[0022] 圖5 (a)為對附加相位Δ%提取的歸一化理論結(jié)果，圖5 (b)為對附加相位提取 的歸一化模擬結(jié)果；
[0023] 圖6 (a)為實(shí)驗(yàn)獲取的初始條紋圖，圖6 (b)為實(shí)驗(yàn)獲取的附加 π /7相位后的條紋 圖；
[0024] 圖7為實(shí)驗(yàn)結(jié)果圖。
【具體實(shí)施方式】
[0025] 下面結(jié)合附圖與實(shí)施例對本發(fā)明作進(jìn)一步說明。
[0026] 利用光流場進(jìn)行物體變形相位測量的方法，包括以下步驟：
[0027] 步驟一：根據(jù)光流場理論的基本假設(shè)，獲得兩幀連續(xù)圖像之間的運(yùn)動(dòng)矢量的兩個(gè) 分量u和V ;
[0028] 步驟二：由窗口傅里葉變換法提取第一幀圖像條紋的橫向頻率fx。和縱向頻率 fyO；
[0029] 步驟三：利用所述步驟一得到的運(yùn)動(dòng)矢量的兩個(gè)分量U和V及所述步驟二得到條 紋的橫向頻率fx(:和縱向頻率f y(:計(jì)算得到相位的改變量。
[0030] 1運(yùn)動(dòng)矢量的計(jì)算
[0031] 設(shè)在t時(shí)刻像素點(diǎn)（X，y)處灰度值為I (X，y, t)，在（t+Δ t)時(shí)刻該點(diǎn)運(yùn)動(dòng)到新的 位置（χ+ Λ X，y+ Δ y)，灰度值記為I (Χ+ Δ X，y+ Δ y)。根據(jù)圖像一致性假設(shè)，即圖像沿著運(yùn)動(dòng) 軌跡亮度保持不變，滿足

[0032] I (X，y，t) = I (X+Δ X，y+Δ y，t+Δ t) (I)
[0033] 設(shè)U和v分別為該點(diǎn)光流矢量沿x和y方向的兩個(gè)分量，且 I (X+ Δ X，y+ Δ y, t+ Δ t)用泰勒公式展開，得到：
[0040]
Ix、仁和I沒別為圖像I對于X，y，t的偏導(dǎo)數(shù)，它們的值可 以用圖像序列中相鄰圖像目標(biāo)像素點(diǎn)的一階差分估計(jì)得到。公式（5)即為光流基本等式， 由于只有一個(gè)方程，只能求得光流沿梯度方向上的值Ix、1￥和11，而u，V的解是非唯一的，故 需要加入約束條件才能求解速度矢量。根據(jù)Horn-Schunck算法計(jì)算光流運(yùn)動(dòng)場。
[0041 ] Horn-Schunck算法的基本思想是在求解光流時(shí)要求光流本身盡可能的平滑。所謂 平滑，就是在給定的鄰域內(nèi)V2?+V:2v盡可能的小，即
[0043] 綜合式（5)和式（6)，Horn-Schunck算法將光流（u, V)的計(jì)算歸結(jié)為式（7)的變 分問題。
[0045] 可以得到相應(yīng)的歐拉-拉格朗日方程，并利用高斯-賽德爾方法進(jìn)行求解，得到圖 像上每點(diǎn)第（Π +1)次迭代估計(jì)的（Un+1，vn+1)為：
[0048] 根據(jù)公式（8)和（9)迭代所得到的u和V即為兩幀圖像之間的運(yùn)動(dòng)矢量場的兩個(gè) 分量。
[0049] 2條紋頻率的計(jì)算
[0050] 為提取條紋頻率，引入以下窗口傅里葉運(yùn)算：
[0052] 式（10)中，（ξ，η )為空間一點(diǎn)（X，y)對應(yīng)的頻域坐標(biāo)，SI ( μ，V ; ξ，η ) 是I (x，y)經(jīng)過窗口傅里葉變換后對應(yīng)頻域的結(jié)果。g(x，y)為窗口函數(shù)，取高斯函數(shù)
OdP σ y為高斯分布的標(biāo)準(zhǔn)差，其取值決定著高斯窗口的大 小。
[0054] 其中（fx。，fy。）為使I SI (X，y，ξ，n) I取最大值的（ξ，n)，同時(shí)也是干涉條紋圖 中每一點(diǎn)的橫向與縱向條紋頻率。
[0055] 故條紋頻率fx。和fy。的提取方法可以表述為：首先要將條紋圖的條紋頻率 (ξ，n)在其取值范圍內(nèi)均勻離散化，然后把這些值分別帶入|SI(x，y，ξ，11)1，其中能使 si(x，y，ξ，n) I取最大值的