[发明专利]用于电子散斑干涉条纹图的方向偏微分方程滤波方法

权利要求书

1.一种用于电子散斑干涉条纹图的方向偏微分方程滤波方法，包括下列步骤：

步骤1：输入一幅电子散斑干涉条纹图像u；

步骤2：对该电子散斑干涉条纹图像u进行离散化：假设图像u的大小为M×N，u_i,j(1≤i≤M,1≤j≤N)为(i,j)点处的灰度值，时间步长为Δt，在方程的演化过程中，t_n=nΔt时刻的演化图像u(i,j,nΔt)表示为用一个前向差分来计算，即构造方向偏微分方程的离散格式为：

ui,jn+1=ui,jn+Δt·((uxx)i,jncos2θi,j+(uyy)i,jnsin2θi,j+2(uxy)i,jnsinθi,jcosθi,j)]]>

其中，θ_i,j表示条纹方向与x轴方向的夹角，u_x、u_y分别表示图像u沿着x轴、y轴方向的一阶偏导数，u_xx、u_yy分别表示图像u沿着x轴、y轴方向的二阶偏导数；

步骤3：计算条纹方向与x轴方向的夹角θ_i,j；

步骤4：自适应得到时间步长Δt；

步骤5：自适应得到迭代次数N_c；

步骤6：对于每次迭代，按照公式(ux)i,jn=ui+1,jn-ui-1,jn2,(uy)i,jn=ui,j+1n-ui,j-1n2,]]>(uxx)i,jn=ui+1,jn-2ui,jn+ui-1,jn,(uyy)i,jn=ui,j+1n+ui,j-1n-2ui,jn]]>以及

(uxy)i,jn=ui+1,j+1n-ui+1,j-1n-ui-1,j+1n+ui-1,j-1n4,]]>

求出图像u每个像素的一阶偏导数u_x，u_y以及二阶偏导数u_xx，u_xy和u_yy；

步骤7：基于步骤2中的方向偏微分方程的离散格式，求出图像u每个像素的数值解

步骤8：重复步骤6和步骤7，直到达到设置的最大迭代次数N_c停止迭代，此时的数值解即为滤波图像。