本發(fā)明屬于物體應(yīng)變率場測量
技術(shù)領(lǐng)域：
，具體涉及一種基于速度線性擬合的平面內(nèi)應(yīng)變率場的光學(xué)測量方法。
背景技術(shù)：
：物體變形測量是實(shí)驗(yàn)力學(xué)的重要內(nèi)容；與其他測量方法相比，光學(xué)測量方法由于其具有非接觸、全場、高空間分辨率和高測量精度等優(yōu)點(diǎn)，在實(shí)驗(yàn)力學(xué)中扮演著不可替代的重要角色；在實(shí)際應(yīng)用中，物體變形測量常常是指物體平面內(nèi)變形的測量；數(shù)字圖像相關(guān)方法是光學(xué)測量方法中的一種重要方法，是對變形前后采集的物體表面的兩幅圖像(散斑場)進(jìn)行相關(guān)處理，以實(shí)現(xiàn)物體變形測量；該方法不僅可進(jìn)行物體平面內(nèi)變形測量，也能進(jìn)行離面變形測量；試件或結(jié)構(gòu)在受載條件下，經(jīng)過一定階段之后，應(yīng)變率的分布并不均勻；在應(yīng)變率較高區(qū)域，應(yīng)變增加較快，該區(qū)域往往是未來裂紋發(fā)生的區(qū)域；通過探測應(yīng)變率較高區(qū)域的時(shí)空分布規(guī)律，可提前獲知未來破壞區(qū)域的具體位置，這對于材料或結(jié)構(gòu)破壞機(jī)理分析和災(zāi)害預(yù)防研究大有裨益；采用數(shù)字圖像相關(guān)方法，物體平面內(nèi)變形測量方法主要包括：1)基于牛頓—拉菲遜(N-R)方法的數(shù)字圖像相關(guān)方法，同時(shí)獲得物體平面內(nèi)的位移和應(yīng)變，但應(yīng)變計(jì)算結(jié)果的誤差較大；2)基于中心差分方法的數(shù)字圖像相關(guān)方法，通過對位移場進(jìn)行中心差分來獲得應(yīng)變場，但位移場中包含的噪聲會使應(yīng)變計(jì)算結(jié)果的可信度降低；3)基于有限元平滑方法的數(shù)字圖像相關(guān)方法，利用有限元方法對位移場進(jìn)行平滑來獲得應(yīng)變場，但數(shù)學(xué)形式和編程執(zhí)行都較復(fù)雜，不利于推廣和應(yīng)用。技術(shù)實(shí)現(xiàn)要素：針對現(xiàn)有技術(shù)的不足，本發(fā)明提出一種基于速度線性擬合的平面內(nèi)應(yīng)變率場的光學(xué)測量方法，以達(dá)到提高準(zhǔn)確率、可信度和計(jì)算速度目的。一種基于速度線性擬合的平面內(nèi)應(yīng)變率場的光學(xué)測量方法，包括以下步驟：步驟1、采用拍攝設(shè)備采集被測物體某一平面在受載過程中的多張圖像；步驟2、在采集的多張圖像中選定測量圖像，并確定每張測量圖像的拍攝時(shí)刻，進(jìn)而獲得相鄰的兩個(gè)時(shí)刻的測量圖像的間隔時(shí)間；步驟3、將首張測量圖像所在平面作為xoy面，建立三維直角坐標(biāo)系oxyz，并將首張測量圖像中的被測物體所在區(qū)域離散化成無重疊且無間隙的三角形像素塊；步驟4、獲得三角形像素塊三個(gè)頂點(diǎn)的坐標(biāo)，在三角形像素塊的每個(gè)頂點(diǎn)位置布置1個(gè)測點(diǎn)，并在xoy面中每個(gè)三角形像素塊的非頂點(diǎn)位置布置N個(gè)測點(diǎn)，N大于等于0，獲取該測點(diǎn)在x軸和y軸上的坐標(biāo)；步驟5、采用數(shù)字圖像相關(guān)方法，獲得被測物體在受載過程中平面內(nèi)的測點(diǎn)的位移場，根據(jù)相鄰的兩個(gè)時(shí)刻的測量圖像的間隔時(shí)間，獲得平面內(nèi)的測點(diǎn)的速度場；步驟6、根據(jù)每個(gè)三角形像素塊中每個(gè)測點(diǎn)的x坐標(biāo)、y坐標(biāo)及速度，獲得三維直角坐標(biāo)系中每個(gè)三角形像素塊的兩組擬合數(shù)據(jù)點(diǎn)；步驟7、選擇誤差種類，根據(jù)每個(gè)三角形像素塊的每組擬合數(shù)據(jù)點(diǎn)，獲得每個(gè)三角形像素塊的擬合平面的系數(shù)；步驟8、根據(jù)獲得的系數(shù)與速度偏導(dǎo)數(shù)的關(guān)系，獲得每個(gè)三角形像素塊的應(yīng)變率，即獲得被測物體平面內(nèi)的應(yīng)變率場。步驟6所述的根據(jù)每個(gè)三角形像素塊中每個(gè)測點(diǎn)的x坐標(biāo)、y坐標(biāo)及速度，獲得三維直角坐標(biāo)系中每個(gè)三角形像素塊的兩組擬合數(shù)據(jù)點(diǎn)；具體包括以下步驟：步驟6.1、在xoy面上的每個(gè)三角形像素塊中，將每個(gè)測點(diǎn)的x坐標(biāo)作為擬合數(shù)據(jù)點(diǎn)的x坐標(biāo)，將每個(gè)測點(diǎn)的y坐標(biāo)作為擬合數(shù)據(jù)點(diǎn)的y坐標(biāo)，將每個(gè)測點(diǎn)的速度在x軸方向的分量作為擬合數(shù)據(jù)點(diǎn)的z坐標(biāo)，在三維直角坐標(biāo)系oxyz中構(gòu)造每個(gè)三角形像素塊的第一組擬合數(shù)據(jù)點(diǎn)；步驟6.2、在xoy面上的每個(gè)三角形像素塊中，將每個(gè)測點(diǎn)的x坐標(biāo)作為擬合數(shù)據(jù)點(diǎn)的x坐標(biāo)、將每個(gè)測點(diǎn)的y坐標(biāo)作為擬合數(shù)據(jù)點(diǎn)的y坐標(biāo)，將每個(gè)測點(diǎn)的速度在y軸方向的分量作為擬合數(shù)據(jù)點(diǎn)的z坐標(biāo)，在三維直角坐標(biāo)系oxyz中構(gòu)造每個(gè)三角形像素塊的第二組擬合數(shù)據(jù)點(diǎn)。步驟7所述的選擇誤差種類，根據(jù)每個(gè)三角形像素塊的每組擬合數(shù)據(jù)點(diǎn)，獲得每個(gè)三角形像素塊的擬合平面的系數(shù)；平面如下：{u·=A1x+B1y+C1v·=A2x+B2y+C2---(1)]]>其中，表示每個(gè)三角形像素塊的第一個(gè)擬合平面上的每個(gè)點(diǎn)在z方向的坐標(biāo)，表示每個(gè)三角形像素塊的第二個(gè)擬合平面上的每個(gè)點(diǎn)在z方向的坐標(biāo)，x表示在x方向的坐標(biāo)，y表示在y方向的坐標(biāo)，A1、B1和C1表示每個(gè)三角形像素塊的第一個(gè)擬合平面的系數(shù)，A2、B2和C2表示每個(gè)三角形像素塊的第二個(gè)擬合平面的系數(shù)。步驟8所述的獲得的系數(shù)與速度偏導(dǎo)數(shù)的關(guān)系，具體如下：∂u·∂x=A1,∂u·∂y=B1∂v·∂x=A2,∂v·∂y=B2---(2)]]>其中，表示每個(gè)三角形像素塊的第一個(gè)擬合平面上的每個(gè)點(diǎn)在z方向的坐標(biāo)，表示每個(gè)三角形像素塊的第二個(gè)擬合平面上的每個(gè)點(diǎn)在z方向的坐標(biāo)，x表示在x方向的坐標(biāo)，y表示在y方向的坐標(biāo)，A1和B1表示每個(gè)三角形像素塊的第一個(gè)擬合平面的系數(shù)，A2和B2表示每個(gè)三角形像素塊的第二個(gè)擬合平面的系數(shù)；所述的每個(gè)三角形像素塊的應(yīng)變率，公式如下：e·ij=A112(B1+A2)12(B1+A2)B2---(3)]]>其中，表示每個(gè)三角形像素塊的應(yīng)變率張量；i表示坐標(biāo)軸i，為x軸或y軸；j表示坐標(biāo)軸j，為x軸或y軸。本發(fā)明優(yōu)點(diǎn)：本發(fā)明提出一種基于速度線性擬合的平面內(nèi)應(yīng)變率場的光學(xué)測量方法，基于數(shù)據(jù)擬合的思想，優(yōu)于現(xiàn)有的基于中心差分方法的數(shù)字圖像相關(guān)方法，其過程簡單、步驟便捷且易于理解，僅需求解擬合平面方程的系數(shù)便可獲得結(jié)果，提高了測量速度和編程便捷性；本發(fā)明考慮了離散像素塊內(nèi)部測點(diǎn)的信息，提高了最終結(jié)果的準(zhǔn)確率和可信度，并增加了測點(diǎn)選擇的靈活性；本發(fā)明在固體實(shí)驗(yàn)力學(xué)領(lǐng)域有廣泛的應(yīng)用前景。附圖說明圖1為本發(fā)明一種實(shí)施例的物體平面內(nèi)的應(yīng)變率場的光學(xué)測量方法流程圖；圖2為本發(fā)明一種實(shí)施例的選定的單軸壓縮過程中試件表面的測量圖像示意圖，其中，圖(a)為第一張測量圖像，圖(b)為第二張測量圖像，圖(c)為第三張測量圖像；圖3為本發(fā)明一種實(shí)施例的首張測量圖像中的被測物體所在區(qū)域離散化后的三角形像素塊分布圖；圖4為本發(fā)明一種實(shí)施例的每個(gè)測點(diǎn)在首張測量圖像上的位置示意圖；圖5為本發(fā)明一種實(shí)施例的由測點(diǎn)的速度和坐標(biāo)獲得三角形像素塊擬合平面的幾何示意圖；圖6為本發(fā)明一種實(shí)施例的采用本發(fā)明方法獲得的由第二至第三張圖像的應(yīng)變率場填充圖，其中，圖(a)為x方向的線應(yīng)變率場填充圖，圖(b)為y方向的線應(yīng)變率場填充圖，圖(c)為面內(nèi)剪切應(yīng)變率場填充圖；圖7為采用基于N-R方法的數(shù)字圖像相關(guān)方法獲得的由第二至第三張圖像的應(yīng)變率場填充圖，其中，圖(a)為x方向的線應(yīng)變率場填充圖，圖(b)為y方向的線應(yīng)變率場填充圖，圖(c)為面內(nèi)剪切應(yīng)變率場填充圖；圖8為采用基于中心差分方法的數(shù)字圖像相關(guān)方法獲得的由第二至第三張圖像的應(yīng)變率場填充圖，其中，圖(a)為x方向的線應(yīng)變率場填充圖，圖(b)為y方向的線應(yīng)變率場填充圖，圖(c)為面內(nèi)剪切應(yīng)變率場填充圖；圖中，1為三角形像素塊，2為三角形像素塊的頂點(diǎn)，3為測點(diǎn)，4為三角形像素塊，5為測點(diǎn)，6為擬合數(shù)據(jù)點(diǎn)，7為三角形像素塊，8為擬合平面；a為三角形像素塊的第一頂點(diǎn)，b為三角形像素塊的第二頂點(diǎn)，c為三角形像素塊的第三頂點(diǎn)。具體實(shí)施方式下面結(jié)合附圖對本發(fā)明一種實(shí)施例做進(jìn)一步說明。一種基于速度線性擬合的平面內(nèi)應(yīng)變率場的光學(xué)測量方法，方法流程圖如圖1所示，包括以下步驟：步驟1、采用拍攝設(shè)備采集被測物體某一平面在受載過程中的多張圖像；本發(fā)明實(shí)施例中，被測物體為長方體試件，高度為90mm，寬度為50mm，將試件置于平臺上，在試件的一個(gè)表面上采用涂料制作人工散斑，在上端面進(jìn)行位移控制加載，采集該表面的散斑圖；步驟2、在采集的多張圖像中選定測量圖像，并確定每張測量圖像的拍攝時(shí)刻，進(jìn)而獲得相鄰的兩個(gè)時(shí)刻的測量圖像的間隔時(shí)間；本發(fā)明實(shí)施例中，在采集的多張圖像中選定3張測量圖像，選定的測量圖像如圖2中圖(a)至圖(c)所示，獲得每張測量圖像的拍攝時(shí)刻，進(jìn)而獲得首張測量圖像與第二張測量圖像的拍攝間隔時(shí)間為198s，獲得第二張測量圖像與第三張測量圖像的拍攝間隔時(shí)間為2s；步驟3、將首張測量圖像所在平面作為xoy面，建立三維直角坐標(biāo)系oxyz，并將首張測量圖像中的被測物體所在區(qū)域離散化成無重疊且無間隙的三角形像素塊；本發(fā)明實(shí)施例中，如圖3所示(圖中，1表示三角形像素塊，2表示三角形像素塊的頂點(diǎn)，a表示三角形像素塊的第一頂點(diǎn)，b表示三角形像素塊的第二頂點(diǎn)，c表示三角形像素塊的第三頂點(diǎn))，將首張測量圖像中的被測物體所在區(qū)域離散化成5762個(gè)緊密排列的、無重疊且無間隙的三角形像素塊；步驟4、獲得三角形像素塊三個(gè)頂點(diǎn)的坐標(biāo)，在三角形像素塊的每個(gè)頂點(diǎn)位置布置1個(gè)測點(diǎn)，并在xoy面中每個(gè)三角形像素塊的非頂點(diǎn)位置布置N個(gè)測點(diǎn)，N大于等于0，獲取該測點(diǎn)在x軸和y軸上的坐標(biāo)；本發(fā)明實(shí)施例中，如圖4所示(圖中3表示測點(diǎn)，4表示三角形像素塊)，在xoy面中每個(gè)三角形像素塊的非頂點(diǎn)位置布置1個(gè)測點(diǎn)；步驟5、采用數(shù)字圖像相關(guān)方法，獲得被測物體在受載過程中平面內(nèi)的測點(diǎn)的位移場，根據(jù)相鄰的兩個(gè)時(shí)刻的測量圖像的間隔時(shí)間，獲得平面內(nèi)的測點(diǎn)的速度場；具體如下：步驟5.1、根據(jù)每個(gè)測點(diǎn)在首張測量圖像上的位置，采用數(shù)字圖像相關(guān)方法獲得每個(gè)測點(diǎn)在第二張和第三張測量圖像上的位置；步驟5.2、根據(jù)每個(gè)測點(diǎn)在首張測量圖像上的位置和每個(gè)測點(diǎn)在第二張測量圖像上的位置，獲得每個(gè)測點(diǎn)在第二張測量圖像上的位移，采用相同方法獲得每個(gè)測點(diǎn)在第三張測量圖像上的位移；步驟5.3、獲得第二張測量圖像和第三張圖像上的每個(gè)測點(diǎn)的位移差，并根據(jù)第二張測量圖像與第三張測量圖像的拍攝間隔時(shí)間，獲得每個(gè)測點(diǎn)的速度，即獲得平面內(nèi)的測點(diǎn)的速度場；步驟6、根據(jù)每個(gè)三角形像素塊中每個(gè)測點(diǎn)的x坐標(biāo)、y坐標(biāo)及速度，獲得三維直角坐標(biāo)系中每個(gè)三角形像素塊的兩組擬合數(shù)據(jù)點(diǎn)；本發(fā)明實(shí)施例中，如圖5中圖(a)所示，圖中5表示測點(diǎn)，6表示擬合數(shù)據(jù)點(diǎn)，7表示三角形像素塊，獲得三維直角坐標(biāo)系中每個(gè)三角形像素塊的兩組擬合數(shù)據(jù)點(diǎn)，具體包括以下步驟：步驟6.1、在xoy面上的每個(gè)三角形像素塊中，將每個(gè)測點(diǎn)的x坐標(biāo)作為擬合數(shù)據(jù)點(diǎn)的x坐標(biāo)，將每個(gè)測點(diǎn)的y坐標(biāo)作為擬合數(shù)據(jù)點(diǎn)的y坐標(biāo)，將每個(gè)測點(diǎn)的速度在x軸方向的分量作為擬合數(shù)據(jù)點(diǎn)的z坐標(biāo)，在三維直角坐標(biāo)系oxyz中構(gòu)造每個(gè)三角形像素塊的第一組擬合數(shù)據(jù)點(diǎn)；步驟6.2、在xoy面上的每個(gè)三角形像素塊中，將每個(gè)測點(diǎn)的x坐標(biāo)作為擬合數(shù)據(jù)點(diǎn)的x坐標(biāo)、將每個(gè)測點(diǎn)的y坐標(biāo)作為擬合數(shù)據(jù)點(diǎn)的y坐標(biāo)，將每個(gè)測點(diǎn)的速度在y軸方向的分量作為擬合數(shù)據(jù)點(diǎn)的z坐標(biāo)，在三維直角坐標(biāo)系oxyz中構(gòu)造每個(gè)三角形像素塊的第二組擬合數(shù)據(jù)點(diǎn)；步驟7、選擇誤差種類，根據(jù)每個(gè)三角形像素塊的每組擬合數(shù)據(jù)點(diǎn)，獲得每個(gè)三角形像素塊的擬合平面的系數(shù)；所述的誤差種類包括最大誤差、平均誤差和均方根誤差等，屬于本領(lǐng)域的公知常識，可根據(jù)實(shí)際需求自行選擇；步驟7.1、確定每個(gè)三角形像素塊的兩個(gè)擬合平面方程；本發(fā)明實(shí)施例中，擬合出的平面如圖5中圖(b)所示，其中8表示擬合平面；平面方程如下：{u·=A1x+B1y+C1v·=A2x+B2y+C2---(1)]]>其中，表示每個(gè)三角形像素塊的第一個(gè)擬合平面上的每個(gè)點(diǎn)在z方向的坐標(biāo)，表示每個(gè)三角形像素塊的第二個(gè)擬合平面上的每個(gè)點(diǎn)在z方向的坐標(biāo)，x表示在x方向的坐標(biāo)，y表示在y方向的坐標(biāo)，A1、B1和C1表示每個(gè)三角形像素塊的第一個(gè)擬合平面的系數(shù)，A2、B2和C2表示每個(gè)三角形像素塊的第二個(gè)擬合平面的系數(shù)；步驟7.2、選擇均方根誤差，根據(jù)每個(gè)三角形像素塊的兩組擬合數(shù)據(jù)點(diǎn)，獲得每個(gè)三角形像素塊的兩個(gè)擬合平面的均方根誤差，所述均方根誤差最小為擬合數(shù)據(jù)點(diǎn)到擬合平面的垂直距離(即擬合數(shù)據(jù)點(diǎn)沿與z軸平行的方向到擬合平面的距離)的平方和，公式如下：E1=Σk=14(A1xk+B1yk+C1-u·k)2E2=Σk=14(A2xk+B2yk+C2-v·k)2---(4)]]>其中，E1表示每個(gè)三角形像素塊的第一個(gè)擬合平面的均方根誤差；E2表示每個(gè)三角形像素塊的第二個(gè)擬合平面的均方根誤差；xk表示每個(gè)三角形像素塊的擬合數(shù)據(jù)點(diǎn)在x方向的坐標(biāo)；yk表示每個(gè)三角形像素塊的擬合數(shù)據(jù)點(diǎn)在y方向的坐標(biāo)；表示每個(gè)三角形像素塊的第一組擬合數(shù)據(jù)點(diǎn)在z方向的坐標(biāo)；表示每個(gè)三角形像素塊的第二組擬合數(shù)據(jù)點(diǎn)在z方向的坐標(biāo)；k表示每組擬合數(shù)據(jù)點(diǎn)中的第k個(gè)擬合數(shù)據(jù)點(diǎn)，取值范圍為[1，2，3，4]；步驟7.3、通過保證均方根誤差最小，獲得每個(gè)三角形像素塊的兩組最小二乘擬合平面的法方程，具體為：所述保證均方根誤差最小為使偏導(dǎo)數(shù)和和為零，公式如下：(Σk=14xk2)A1+(Σk=14xkyk)B1+(Σk=14xk)C1=Σk=14xku·k(Σk=14xkyk)A1+(Σk=14yk2)B1+(Σk=14yk)C1=Σk=14yku·k(Σk=14xk)A1+(Σk=14yk)B1+4C1=(Σk=14u·k)---(5)]]>(Σk=14xk2)A2+(Σk=14xkyk)B2+(Σk=14xk)C2=Σk=14xkv·k(Σk=14xkyk)A2+(Σk=14yk2)B2+(Σk=14yk)C2=Σk=14ykv·k(Σk=14xk)A2+(Σk=14yk)B2+4C2=(Σk=14v·k)---(6)]]>步驟7.4、求解公式(5)和公式(6)，獲得最小二乘擬合平面的系數(shù)A1、B1、C1、A2、B2和C2；步驟8、根據(jù)獲得的系數(shù)與速度偏導(dǎo)數(shù)的關(guān)系，獲得每個(gè)三角形像素塊的應(yīng)變率，即獲得被測物體平面內(nèi)的應(yīng)變率場；具體包括以下步驟：步驟8.1、將公式(1)中的兩個(gè)最小二乘擬合平面的方程分別求偏導(dǎo)，獲得最小二乘擬合平面的系數(shù)與速度偏導(dǎo)數(shù)的關(guān)系；公式如下：∂u·∂x=A1,∂u·∂y=B1∂v·∂x=A2,∂v·∂y=B2---(2)]]>其中，表示每個(gè)三角形像素塊的第一個(gè)擬合平面上的每個(gè)點(diǎn)在z方向的坐標(biāo)，表示每個(gè)三角形像素塊的第二個(gè)擬合平面上的每個(gè)點(diǎn)在z方向的坐標(biāo)，x表示在x方向的坐標(biāo)，y表示在y方向的坐標(biāo)，A1和B1表示每個(gè)三角形像素塊的第一個(gè)擬合平面的系數(shù)，A2和B2表示每個(gè)三角形像素塊的第二個(gè)擬合平面的系數(shù)；步驟8.2、利用幾何方程獲得速度偏導(dǎo)數(shù)與應(yīng)變率之間的關(guān)系；公式如下：e·ij=12(∂u·i∂Xj+∂u·j∂Xi)---(7)]]>其中，表示每個(gè)三角形像素塊的應(yīng)變率張量；表示在i方向的速度分量；Xj表示在j方向的坐標(biāo)；表示在j方向的速度分量；Xi表示在i方向的坐標(biāo)；i表示坐標(biāo)軸i，為x軸或y軸；j表示坐標(biāo)軸j，為x軸或y軸；步驟8.3、將應(yīng)變率公式(7)轉(zhuǎn)換為工程上的常用形式；具體為：每個(gè)三角形像素塊的應(yīng)變率是一個(gè)張量，對于二維問題，表示為：e·ij=∂u·∂Xx12(∂u·∂Xy+∂v·∂Xx)12(∂u·∂Xy+∂v·∂Xx)∂v·∂Xy---(8)]]>其中，Xx表示在x方向的坐標(biāo)，即為x；Xy表示在y方向的坐標(biāo)，即為y；令則可以表示為：e·ij=ϵ·x12γ·xy12γ·xyϵ·y---(9)]]>其中，表示x方向的線應(yīng)變率；表示面內(nèi)剪切應(yīng)變率；表示y方向的線應(yīng)變率；步驟8.4、將公式(2)與公式(7)進(jìn)行合并，獲得每個(gè)三角形像素塊的應(yīng)變率，即獲得被測物體平面內(nèi)的應(yīng)變率場，公式如下：e·ij=A112(B1+A2)12(B1+A2)B2---(3)]]>本發(fā)明實(shí)施例中，圖6為采用本發(fā)明方法獲得的由第二至第三張圖像的應(yīng)變率場填充圖(包括圖(a)、圖(b)和圖(c))；為了進(jìn)行對比，圖7給出了采用基于N-R方法的數(shù)字圖像相關(guān)方法獲得的由第二至第三張圖像的應(yīng)變率場填充圖(包括圖(a)、圖(b)和圖(c))；圖8給出了采用基于中心差分方法的數(shù)字圖像相關(guān)方法獲得的由第二至第三張圖像的應(yīng)變率場填充圖(包括圖(a)、圖(b)和圖(c))；由圖6～圖8可以發(fā)現(xiàn)，各種應(yīng)變率場呈不均勻分布，尤其是和可以觀察到一些呈帶狀的應(yīng)變率高值區(qū)或低值區(qū)，這些位置的應(yīng)變變化較快，是未來宏觀裂紋發(fā)生的區(qū)域；本發(fā)明實(shí)施例中，利用基于中心差分方法的數(shù)字圖像相關(guān)方法和本發(fā)明方法獲得的和的測量范圍有一定差別，都優(yōu)于利用基于N-R方法的數(shù)字圖像相關(guān)方法獲得的測量結(jié)果；具體而言，利用基于中心差分方法的數(shù)字圖像相關(guān)方法獲得的的測量范圍為(-7～19)×10-3s-1，的測量范圍為(-16～6)×10-3s-1，的測量范圍為(-17～8)×10-3s-1，而采用本發(fā)明方法獲得的的測量范圍為(-15～33)×10-3s-1，的測量范圍為(-21～9)×10-3s-1，的測量范圍為(-29～17)×10-3s-1；相比之下，采用本發(fā)明方法獲得的測量范圍大于利用基于中心差分方法的數(shù)字圖像相關(guān)方法獲得的測量范圍，這使主要的應(yīng)變率高值區(qū)更加突顯，同時(shí)，高值區(qū)的寬度較窄，這意味著本發(fā)明方法能對應(yīng)變率場刻畫得更加細(xì)致；從采用本發(fā)明方法獲得的測量結(jié)果中，可以觀察到較為清晰的呈帶狀的應(yīng)變率高值區(qū)和低值區(qū)，而利用基于中心差分方法的數(shù)字圖像相關(guān)方法獲得的測量結(jié)果較為模糊，應(yīng)變率高值區(qū)顯得較為斑駁；綜上所述，采用本發(fā)明方法獲得的測量結(jié)果優(yōu)于其他兩種方法。當(dāng)前第1頁1 2 3