[发明专利]基于能量最小化框架的文档图像倾斜度检测与纠正方法

权利要求书

1.基于能量最小化框架的文档图像倾斜度检测与纠正方法，该方法包括三个过程：

(1)计算前景像素状态信息

在计算前景像素状态信息之前，首先对扫描得到的文档图像进行二值化，用黑色表示前景，白色表示背景像素，给定一个前二值文档图像I之后，前景像素的状态计算过程如下：一个边界框定义为一个文档图像的边界，用P记作整个前景像素的集合，(W，H)记作文档图像I的大小，那么对每一个前景像素p∈P，它的状态信息为s_p＝(x_p，y_p，w_p，h_p)，其中x_p，y_p，w_p，h_p分别为p到图像最左、最上、最右和最下边的距离；

(2)利用直线拟合技术估算初始倾斜度

接下来利用文档图像中最外围的前景像素的状态信息来估算初始倾斜度，一个边界框有四个边：上、下、左和右，对每一边都能得到其最外围的前景像素状态信息子集，用以下方式得到上边最外围前景像素状态信息子集，记为

TS=∪i=1Wsi---(1)]]>

y_i∈s_i∩y_i＝min{y_p|y_p∈s_p∩x_p∈s_p∩x_p＝i}

将TS中每个元素的x_p做为x坐标，y_p作为y坐标画图，图中下面的点拟合成一条直线，边界框剩下的三边都经过该处理，将使用直线拟合的技术来估算文档图像的初始倾斜度，在直线拟合之前先对TP进行采样，用如下方式将TP划分为N个互不重叠的部分STS_i：

TS=∪i=1NSTSi]]>

STSi={sp|xp∈sp∩(i-1)×WN+1≤xp≤i×WN}]]>

N＝32，接下来用如下方式构建一个子集FTS，即计算每个部分STS_i中y_i最小的那个前景像素状态信息：

FTS=∪i=1N(si∩si∈STSi)---(3)]]>

y_i∈s_i∩y_i=min{y_p|y_p∈s_p∩s_p∈STS_i}

进行采样操作完之后，需要通过以下方式进一步消除没用的状态信息得到有效的用来进行直线拟合的状态信息VTS

VTS={sp|sp∈FTS∩yp∈sp∩yp<H3}---(4)]]>

然后采用穷举搜索的方式在VTS中做直线拟合直到找到两个状态信息使得有最多的其他状态信息到由这两个状态信息确定的直线之间的距离小于指定的阈值D；

对边界框的四边都进行直线拟合之后得到四条直线，接下来就是找到拟合最好的那条直线，同时该直线对应的倾斜角就是文档图像的初始倾斜角；用{l_t，l_b，l_i，l_r}记作拟合的四条直线，{LS_t，LS_b，LS_i，LS_r}记作靠近相应直线的状态信息，要是某条直线对应的状态信息的个数小于M，在下面的操作中将不再考虑该直线，对每条直线l_i，计算所有直线状态信息LS_i和直线l_i之间的距离之和，然后用如下方式计算比值R_i：

Ri=SDi(f(LSi))2---(5)]]>

其中f(·)计算一个集合中元素的个数，{R_t，R_b，R_i，R_r}中的最小值对应的直线就是最佳拟合的直线，最终最佳拟合的直线对应的倾斜角就是文档图像的初始倾斜角，其中D＝5，M＝5；

(3)使用能量最小化过程计算最终倾斜度

得到初始倾斜角之后，然后用能量最小化过程计算最终倾斜角，如下式所示：

S^=argminSE(S)---(6)]]>

该过程包括能量函数构建和能量最小化两个步骤，如下式所示：

E(S)＝ωE_B(S)+(1-ω)E_F(S)      (7)

ω＝0.98，其中E_B(S)考虑了全局背景信息，如下式所示：

设置和φ(·)为：

Y_i＝{s_p|s_p∈S∩y_p∈s_p∩y_p＝i}     (9)

X_j＝{s_p|s_p∈S∩x_p∈s_p∩x_p＝j}

其中Sgn(·)是一个符号函数，定义为：

E_F(S)反应了全局的前景信息，如下式所示：

E_F(S)＝δ(S)+λ(S)      (11)

设置δ(·)和λ(·)为：

δ(S)=1MY1f(Y)Σk=1f(Y)(f(Yk)-Y‾)2]]>

Y‾=1f(Y)Σk=1f(Y)f(Yk)]]>      (12)

λ(S)=1MX1f(X)Σk=1f(X)(f(Xk)-X‾)2]]>

X‾=1f(X)Σk=1f(X)f(Xk)]]>

其中M_Y＝max{f(Y_i)|Y_i∈Y}，M_X＝max{f(X_i)|X_i∈X}，且f(·)计算一个集合中元素的个数；

构造完能量函数以后，根据直线拟合得到的初始倾斜角和所有前景像素的状态信息，通过反复地计算能量函数和旋转前景像素的状态信息，找到使得能量函数值最小时所旋转的角度，这个角度就是最终的倾斜角；

状态信息的旋转过程如下：

S′＝rotate(S，θ)     (13)

其中rotate(·)计算每个前景像素的状态信息s_p∈S旋转后的结果s_p′，计算过程如下：

xp′=(xp-W2)cosθ-(yp-H2)sinθ+W2]]>

yp′=(xp-W2)sinθ+(yp-H2)cosθ+H2---(14)]]>

w_p′＝W-x_p′，h_p′＝H-y_p′

。