大变化密度电子散斑干涉条纹图的方向和密度处理方法与流程

技术特征：

1.一种大变化密度电子散斑干涉条纹图的方向和密度处理方法，其特征是，步骤如下：

步骤1：输入一幅大变化电子散斑干涉测量ESPI条纹图f；

步骤2：用变分图像分解BL-Hilbert模型把一幅大变化ESPI条纹图f分解成两幅密度均匀的ESPI条纹图：一幅是只包含低密度条纹的ESPI条纹图u；另一幅是只包含高密度条纹的ESPI条纹图v；

步骤3：分别计算低密度ESPI条纹图u和高密度ESPI条纹图v的条纹方向和密度，其中条纹方向采用梯度法，条纹密度采用累积平方差法；

步骤4：整幅图的条纹方向和条纹密度，通过结合选取低密度ESPI条纹图u和高密度ESPI条纹图v的条纹方向和密度而获得。

2.如权利要求1所述的大变化密度电子散斑干涉条纹图的方向和密度处理方法，其特征是，用变分图像分解BL-Hilbert模型分解一幅大变化ESPI条纹图I方法包括以下步骤：

步骤4-1：输入一幅大变化ESPI条纹图f；

步骤4-2：给定初值和参数，设置迭代初值：低密度条纹的ESPI条纹图u＝f，高密度条纹的ESPI条纹图v＝0，频域场ξ＝0，结合实际处理图片特点选择参数λ、μ、N、N0，其中λ和μ是正则化参数，N是迭代次数，N0为迭代过程中的更新标志；

步骤4-3：固定u和v，计算ξ：极小化频域场其中运算符为计算对角阵的二阶范数，Ψ表示在局域傅里叶框架{ψp,k}p,k下对v的分解，p和k分别代表局域窗口的位置和该窗口的频域坐标，Γ(ξ)＝diagl＝(p,k)γp,k(ξ)是由关于频率场ξ的加权系数γp,k(ξ)构成的对角阵，加权系数γp,k(ξ)定义为：

其中，G_σ＝exp(-(x/σ)²/2)，σ为一尺度参数反映纹理频谱ξk与ξ(xp)的偏差，xp是频域场中的一点，p代表位置，当点xp周围没有明显方向性纹理时，对所有k令γp,k＝1，以此来避免虚假纹理的产生；

步骤4-4：固定ξ和v，计算u：

步骤4-5：固定ξ和u，计算v：

步骤4-6：更新条纹f＝un+1+vn+1重复步骤4-3到步骤4-5；

步骤4-7：当n>N停止迭代，输出此时的u和v。