CN107480728A - 一种基于傅里叶残差值的打印文件的鉴别方法 - Google Patents
一种基于傅里叶残差值的打印文件的鉴别方法 Download PDFInfo
- Publication number
- CN107480728A CN107480728A CN201710748676.2A CN201710748676A CN107480728A CN 107480728 A CN107480728 A CN 107480728A CN 201710748676 A CN201710748676 A CN 201710748676A CN 107480728 A CN107480728 A CN 107480728A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msubsup
- mfrac
- msup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000012850 discrimination method Methods 0.000 title claims abstract description 16
- 238000000034 method Methods 0.000 claims abstract description 13
- 238000000605 extraction Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000012549 training Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000002596 correlated effect Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 5
- 239000000284 extract Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种基于傅里叶残差值的打印文档鉴别方法,解决打印文档之间没有相同字符或者相同字符较少的情况下的鉴别问题。具体方法是:获取一张打印文档的文档图像;提取文档图像的文本行,找到最大的打印区域;提取文档图像的傅里叶残差图,在残差图上提取灰度,熵的统计特征,对这些特征利用分类器进行分类。基于傅里叶残差值的打印文档鉴别方法,可以很好的区分不同打印机之间的打印特性,同时,也不需要打印文档之间有太多相同字符,提高了打印文档鉴别方法的适用范围和准确率。
Description
技术领域
本发明涉及打印文件的源机认定技术的领域,尤其涉及基于傅里叶残差值的打印文件的鉴别方法。
背景技术
由于打印机在生产时的参数配置差异以及日常使用中的损耗情况的不同,因此,打印出的文档也不相同,这是打印文档机源认定技术的基础,所谓打印文档的鉴别,就是判断两份打印文档是否来自同一台打印机。
随着计算机技术的发展和打印机的普及,打印文档取代手写文件作为承载信息的媒介已经成为一种趋势。与此同时,对于打印文件的修改、伪造等犯罪也正在兴起。有越来越多的工作投入到对打印文档的鉴别上。
目前打印文档识别技术中,有的运用小波分析的方法提取单个字符的频域特征做训练样本,使用分类器进行单字符识别。有的基于纹理信息,利用Gabor滤波器提取字符特征进行识别。其识别率很高,但它的维数高,计算量大。但是这些方法都限制在对打印的字符的识别,无法达到打印文档内容无关的识别,因此提出一种打印文档内容无关的鉴别方法是很有必要的。
打印文档字符的傅里叶残差信息承载了其打印机的特性,可利用这些残差信息鉴别两份打印文件是否来自同一打印文件。同时通过对这些残差信息统计建模,以解决打印文档内容无关的鉴别问题。
发明内容
发明目的:本发明所要解决的技术问题是针对现有技术的不足,提供一种更简单泛化能力更好的用于打印文档识别的方法,从而达到对内容无关的打印文档的鉴别。
为了解决上述技术问题,本发明采用的技术方案是:
一种基于傅里叶残差值的打印文档鉴别方法,其特征在于,包括以下步骤:
步骤1,获取一张打印文档的文档图像;
步骤2,提取文档图像的文本行,找到最大的打印区域;
步骤3,提取文档图像的傅里叶残差图,具体方法是:
步骤31,将文档图像转化为灰度图,记为Igray(x,y),其中(x,y)代表文档图像的像素点;
步骤32,对灰度图做傅里叶变换转换到频率域,记为F(u,v),其中(u,v)代表频率域空间坐标:
式中,j为复数虚部单位;
步骤33,对傅里叶图在频率域做拉普拉斯滤波:
G(u,v)=-(u2+v2)F(u,v),
其中,G(u,v)为傅里叶-拉普拉斯滤波增强后的图像;
步骤34,重建灰度图像IRG(x,y)和拉普拉斯滤波后的图像IRFL(x,y),将两个图像取绝对差,得到残差图像Ires(x,y):
Ires(x,y)=IRG(x,y)-IRFL(x,y)
步骤4,对承载打印机特性的残差图提取灰度,熵的统计特征:
步骤41,提取残差图像的灰度特征,设残差图一共有M行N列,将每一行,每一列的均值分别记为和
步骤42,将每一行,每一列的均值分别与残差图中的每一列,每一行计算自相关系数,记为prow,pcol,:
式中,Ires(i,.)为残差图第i行的像素值,Ires(.,j)为残差图第j列的像素值;
步骤43,在均值和自相关系数的基础上提取统计特征,一共14个统计特征,分别为 的方差,偏度和峰度,记为prow、pcol的均值方差,偏度和峰度,记为σprow,σpcol,γ1prow,γ1pcol,γ2prow,γ2pcol:
其中,分别为prow、pcol的均值;
步骤44,提取残差图的信息熵的相关特征,将每一行,每一列的熵,分别记为
步骤45,计算残差图的每一列的值和每一列熵的相关系数,每一行的值和每一行熵的相关系数,分别记为qcol和qrow:
式中,Ires(i,.)为残差图第i行的像素值,Ires(.,j)为残差图第j列的像素值;
步骤46,对残差图的局部信息熵和全局信息熵提取特征,计算相关系数,均值,方差,偏度,峰度的14个统计特征;
其中,分别为qrow、qcol的均值;
步骤46,提取最后一个特征是关于prow,pcol的相对值,记为f29:
步骤5,通过分类器分类,判别文档图像所属的打印机。
所述步骤2还包括对文档图像的倾斜做简单处理。
所述步骤2中,利用Tesseract OCR软件提取文档的文本行。
所述步骤2中,利用扫描仪获取打印文档的文档图像。
步骤1获取的文档图像像素为400dpi。
本发明步骤5中所述通过分类器分类,判别打印文档所属的打印机包括如下步骤:
在提取所有文档的特征的基础上,利用SVM对训练集进行训练,找出每种打印机的判别边界,然后对测试集中的文档进行判断,以判断最有可能属于的打印机类别,以达到鉴别的目的。
与现有技术相比,本发明具有如下优点:
本发明是通过重建文档的灰度图和锐化增强后的图像获取残差信息,对残差信息提取统计特征,来达到一种更简单泛化能力更好的打印机文档鉴别方法。
(1)本发明提出了傅里叶残差,通过在频率域结合了傅里叶的拉普拉斯滤波方法,增强了代表不同打印机噪声模式的系数,同时在重建灰度信息和锐化后的信息的基础上,剔除了打印文档中的字符图像特征,只留下与打印机相关的残差信息。这个残差信息可以很好的用来描述不同打印机的打印特征,即使是通过提取简单的统计特征,结合简单的线性分类器就可以获得不错的效果。
(2)本发明将灰度统计特征和信息熵特征结合起来,利用结合后的特征对打印文档进行鉴别,只通过残差图的灰度统计特征进行判断的效果不是很好,通过挖掘残差图每一行的信息熵相关信息,同时构建残差图整体的信息熵特征,可以对打印机的不同打印模式进行更好的建模,因此,这种做法较以往只利用灰度统计特征的方法更加精确。不需要打印文档之间有太多相同字符,提高了打印文档鉴别方法的适用范围和准确率。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述和或其他方面的优点将会变得更加清楚。
图1为本发明流程图。
图2为输入扫描的文档实例。
图3为输入实例文档的提取文本行。
图4为包含噪声信息的残差图。
具体实施方式:
本发明一种基于傅里叶残差值的打印文档鉴别方法的基本出发点是根据对不同打印机的打印特性,产生的打印噪声差异不同,通过傅里叶残差图提取出这些噪声,并进行统计建模,根据统计特征的差异性,进行分类,从而达到对打印文档的鉴别。
下面结合附图对本发明做更加详细的解释:
1.扫描文档
对一张打印的文档,通过扫描仪扫描,获取像素为400dpi的图像,如图2所示。
2.找出文本行
利用Tesseract OCR软件提取文档的文本行,找到最大的打印区域,同时对文档的倾斜做简单处理。如图3所示。
3.构造残差图
首先将文档图像利用傅里叶变换转换到频域得到F图,然后在频域内做拉普拉斯滤波操作,得到G图,接着对F图做逆傅里叶变换,转换到空域,得到R-G图,同时对G图也做逆傅里叶变换,得到R-FL图,将R-G图减去R-FL图,就能得到包含噪声信息的傅里叶残差图,如图4所示。
4.特征提取
残差图的特征提取主要考虑到要做到打印内容的无关性,结合了灰度和熵的一些统计特征,来挖掘出噪声信息的每一列,每一行的一些相同特征。
对包含噪声信息的残差图提取自相关,熵的统计特征。前14个特征分别是每一行,每一列均值的方差,偏度,峰度,每一行,每一列的均值分别与残差图中的每一列,每一行计算自相关系数的均值,方差,偏度,峰度。另外14个特征是每一行,每一列熵的方差,偏度,峰度,每一行,每一列的熵分别与残差图中的每一列,每一行计算自相关系数的均值,方差,偏度,峰度。最后一个特征是相关系数的一个相对值。
5.分类器
利用SVM对训练集进行训练,找出每种打印机的判别边界,然后对测试集中的文档进行判断,以判断最有可能属于的打印机类别,以达到鉴别的目的。
实施例:
本实施例包括以下部分:
1.输入文本文档
本实例的数据集来自于DPKI打印文档数据集,该数据集主要包括回执单,合同以及论文三类,每类60张,分别用20种不同的打印机打印,共计1200个文档,记为DPKI-20。图2就是其中的一个文档。另一个数据集是我们自己搜集了5台打印机,同时打印DPKI数据集中的文档,然后扫描保存为图像,记为DPKI-5。
2.实验
将数据集中的文档中80%做样本,20%做测试。构造出残差图后提取自相关和熵的统计特征,最后用SVM分类。实验结果用20类打印机的平均精确率,平均召回率,平均F1值三个标准进行评价,如表1所示。
表1分类鉴别结果
数据集 | 平均精确率 | 平均召回率 | 平均F1值 |
DPKI-20 | 0.90 | 0.90 | 0.90 |
DPKI-5 | 0.91 | 0.91 | 0.90 |
从表1可见,本方法在打印机种类较多的情况下,依然可以达到很高的分类准确率,有很好的鲁棒性。
Claims (6)
1.一种基于傅里叶残差值的打印文档鉴别方法,其特征在于,包括以下步骤:
步骤1,获取一张打印文档的文档图像;
步骤2,提取文档图像的文本行,找到最大的打印区域;
步骤3,提取文档图像的傅里叶残差图,具体方法是:
步骤31,将文档图像转化为灰度图,记为Igray(x,y),其中(x,y)代表文档图像的像素点;
步骤32,对灰度图做傅里叶变换转换到频率域,记为F(u,v),其中(u,v)代表频率域空间坐标:
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>v</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>M</mi>
<mi>N</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>I</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>a</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mi>u</mi>
<mi>x</mi>
</mrow>
<mi>M</mi>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mi>y</mi>
</mrow>
<mi>N</mi>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mrow>
式中,j为复数的虚部单位;
步骤33,对傅里叶图在频率域做拉普拉斯滤波:
G(u,v)=-(u2+v2)F(u,v),
其中,G(u,v)为傅里叶-拉普拉斯滤波增强后的图像;
步骤34,重建灰度图像IRG(x,y)和拉普拉斯滤波后的图像IRFL(x,y),将两个图像取绝对差,得到残差图像Ires(x,y):
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>R</mi>
<mi>G</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>M</mi>
<mi>N</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>u</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>v</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>v</mi>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mi>u</mi>
<mi>x</mi>
</mrow>
<mi>M</mi>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mi>y</mi>
</mrow>
<mi>N</mi>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>R</mi>
<mi>F</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>M</mi>
<mi>N</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>u</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>v</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>v</mi>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mi>u</mi>
<mi>x</mi>
</mrow>
<mi>M</mi>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mi>y</mi>
</mrow>
<mi>N</mi>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mrow>
Ires(x,y)=IRG(x,y)-IRFL(x,y)
步骤4,对承载打印机特性的残差图提取灰度,熵的统计特征:
步骤41,提取残差图像的灰度特征,设残差图一共有M行N列,将每一行,每一列的均值分别记为和
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>;</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>j</mi>
<mo>&le;</mo>
<mi>N</mi>
<mo>,</mo>
</mrow>
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<mi>M</mi>
<mo>,</mo>
</mrow>
步骤42,将每一行,每一列的均值分别与残差图中的每一列,每一行计算自相关系数,记为prow,pcol,:
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>C</mi>
<mi>o</mi>
<mi>r</mi>
<mi>r</mi>
<mi>c</mi>
<mi>o</mi>
<mi>e</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>,</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mo>.</mo>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>C</mi>
<mi>o</mi>
<mi>r</mi>
<mi>r</mi>
<mi>c</mi>
<mi>o</mi>
<mi>e</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mo>,</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mo>.</mo>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
式中,Ires(i,.)为残差图第i行的像素值,Ires(.,j)为残差图第j列的像素值;
步骤43,在均值和自相关系数的基础上提取统计特征,一共14个统计特征,分别为 的方差,偏度和峰度,记为prow、pcol的均值方差,偏度和峰度,记为σprow,σpcol,γ1prow,γ1pcol,γ2prow,γ2pcol:
<mrow>
<msubsup>
<mi>&sigma;I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msubsup>
<mi>&sigma;I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&sigma;p</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>3</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msub>
<mi>&sigma;p</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>4</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>3</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>3</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>4</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>4</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>3</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>3</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>4</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>4</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
其中,分别为prow、pcol的均值;
步骤44,提取残差图的信息熵的相关特征,将每一行,每一列的熵,分别记为
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mo>)</mo>
<mo>*</mo>
<msub>
<mi>log</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>;</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<mi>M</mi>
</mrow>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msub>
<mi>log</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>;</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>j</mi>
<mo>&le;</mo>
<mi>N</mi>
</mrow>
步骤45,计算残差图的每一列的值和每一列熵的相关系数,每一行的值和每一行熵的相关系数,分别记为qcol和qrow:
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>C</mi>
<mi>o</mi>
<mi>r</mi>
<mi>r</mi>
<mi>c</mi>
<mi>o</mi>
<mi>e</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>,</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mo>.</mo>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>C</mi>
<mi>o</mi>
<mi>r</mi>
<mi>r</mi>
<mi>c</mi>
<mi>o</mi>
<mi>e</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mo>,</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mo>.</mo>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
式中,Ires(i,.)为残差图第i行的像素值,Ires(.,j)为残差图第j列的像素值;
步骤46,对残差图的局部信息熵和全局信息熵提取特征,计算相关系数,均值,方差,偏度,峰度的14个统计特征;
<mrow>
<msubsup>
<mi>&sigma;E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>5</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msubsup>
<mi>&sigma;E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>6</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>5</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>5</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>6</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>6</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>5</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>5</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>6</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>6</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&sigma;q</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>7</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msub>
<mi>&sigma;q</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>8</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>q</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>7</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>7</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>q</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>8</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>8</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>q</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>7</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>7</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>q</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>8</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mn>8</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
其中,分别为qrow、qcol的均值;
步骤46,提取最后一个特征是关于prow,pcol的相对值,记为f29:
<mrow>
<msub>
<mi>f</mi>
<mn>29</mn>
</msub>
<mo>=</mo>
<mn>100</mn>
<mo>*</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msub>
<mi>p</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msub>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mi>o</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
步骤5,通过分类器分类,判别文档图像所属的打印机。
2.根据权利要求1所述的一种基于傅里叶残差值的打印文档鉴别方法,其特征在于,所述步骤2还包括对文档图像的倾斜做简单处理。
3.根据权利要求1所述的一种基于傅里叶残差值的打印文档鉴别方法,其特征在于,步骤5中所述通过分类器分类,判别打印文档所属的打印机包括如下步骤:
在提取所有文档图像的特征的基础上,利用SVM对训练集进行训练,找出每种打印机的判别边界,然后对测试集中的文档图像进行判断。
4.根据权利要求1所述的一种基于傅里叶残差值的打印文档鉴别方法,其特征在于,所述步骤2中,利用Tesseract OCR软件提取文档的文本行。
5.根据权利要求1所述的一种基于傅里叶残差值的打印文档鉴别方法,其特征在于,所述步骤2中,利用扫描仪获取打印文档的文档图像。
6.根据权利要求1所述的一种基于傅里叶残差值的打印文档鉴别方法,其特征在于,步骤1获取的文档图像像素为400dpi。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710748676.2A CN107480728B (zh) | 2017-08-28 | 2017-08-28 | 一种基于傅里叶残差值的打印文件的鉴别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710748676.2A CN107480728B (zh) | 2017-08-28 | 2017-08-28 | 一种基于傅里叶残差值的打印文件的鉴别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107480728A true CN107480728A (zh) | 2017-12-15 |
CN107480728B CN107480728B (zh) | 2019-02-26 |
Family
ID=60602766
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710748676.2A Active CN107480728B (zh) | 2017-08-28 | 2017-08-28 | 一种基于傅里叶残差值的打印文件的鉴别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107480728B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021042509A1 (zh) * | 2019-09-06 | 2021-03-11 | 平安科技(深圳)有限公司 | 文本图像角度纠偏方法、装置及计算机可读存储介质 |
CN113012018A (zh) * | 2021-04-12 | 2021-06-22 | 中船重工汉光科技股份有限公司 | 一种隐形办公打印文本加密及解密方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101755274A (zh) * | 2007-06-01 | 2010-06-23 | 先进追踪和寻踪公司 | 用于保证文件安全的方法和设备 |
CN102289665A (zh) * | 2011-09-01 | 2011-12-21 | 华中师范大学 | 基于墨粉堆积纹理分析的打印文件鉴别方法 |
CN104700106A (zh) * | 2015-03-30 | 2015-06-10 | 武汉珞珈博研科技有限责任公司 | 一种基于信息挖掘与信息融合的打印文件鉴别方法 |
CN104732548A (zh) * | 2015-04-07 | 2015-06-24 | 武汉大学 | 一种基于纹理合成的打印文件鉴别方法 |
CN105069458A (zh) * | 2015-07-29 | 2015-11-18 | 武汉理工大学 | 基于切比雪夫傅里叶矩的复杂Logo图案特征提取方法 |
US20160371842A1 (en) * | 2015-06-16 | 2016-12-22 | Thomson Licensing | Method and apparatus for computing an estimate position of a micro-image produced by a micro-lens of an array of micro-lenses of an optical acquisition system |
-
2017
- 2017-08-28 CN CN201710748676.2A patent/CN107480728B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101755274A (zh) * | 2007-06-01 | 2010-06-23 | 先进追踪和寻踪公司 | 用于保证文件安全的方法和设备 |
CN102289665A (zh) * | 2011-09-01 | 2011-12-21 | 华中师范大学 | 基于墨粉堆积纹理分析的打印文件鉴别方法 |
CN104700106A (zh) * | 2015-03-30 | 2015-06-10 | 武汉珞珈博研科技有限责任公司 | 一种基于信息挖掘与信息融合的打印文件鉴别方法 |
CN104732548A (zh) * | 2015-04-07 | 2015-06-24 | 武汉大学 | 一种基于纹理合成的打印文件鉴别方法 |
US20160371842A1 (en) * | 2015-06-16 | 2016-12-22 | Thomson Licensing | Method and apparatus for computing an estimate position of a micro-image produced by a micro-lens of an array of micro-lenses of an optical acquisition system |
CN105069458A (zh) * | 2015-07-29 | 2015-11-18 | 武汉理工大学 | 基于切比雪夫傅里叶矩的复杂Logo图案特征提取方法 |
Non-Patent Citations (3)
Title |
---|
DEBADATTA PATI 等: "Speaker Information from Subband Energies of Linear Prediction Residual", 《2010 NATIONAL CONFERENCE ON COMMUNICATIONS (NCC)》 * |
尚士泽: "面向纸质文件的数字化被动取证研究", 《中国博士学位论文全文数据库 社会科学Ⅰ辑》 * |
陈庆虎等: "基于纹理合成的打印文件鉴别", 《科学技术与工程》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021042509A1 (zh) * | 2019-09-06 | 2021-03-11 | 平安科技(深圳)有限公司 | 文本图像角度纠偏方法、装置及计算机可读存储介质 |
CN113012018A (zh) * | 2021-04-12 | 2021-06-22 | 中船重工汉光科技股份有限公司 | 一种隐形办公打印文本加密及解密方法 |
CN113012018B (zh) * | 2021-04-12 | 2022-02-25 | 中船重工汉光科技股份有限公司 | 一种隐形办公打印文本加密及解密方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107480728B (zh) | 2019-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101491074B (zh) | 使用与快速遮盖归一化互相关性相匹配的模板进行文档认证 | |
AU2008335636B2 (en) | Identification and verification of an unknown document according to an Eigen image process | |
Gebhardt et al. | Document authentication using printing technique features and unsupervised anomaly detection | |
RU2658272C2 (ru) | Определение подлинности защищенных документов и мобильное устройство для определения подлинности | |
Zhou et al. | Hough technique for bar charts detection and recognition in document images | |
Bourlai et al. | Restoring degraded face images: A case study in matching faxed, printed, and scanned photos | |
JP2003132358A (ja) | 画像処理方法、装置およびシステム | |
Brink et al. | Towards robust writer verification by correcting unnatural slant | |
CN107480728B (zh) | 一种基于傅里叶残差值的打印文件的鉴别方法 | |
CN106709915B (zh) | 一种图像重采样操作检测方法 | |
Khanna et al. | Intrinsic signatures for scanned documents forensics: effect of font shape and size | |
CN108921006B (zh) | 手写签名图像真伪鉴别模型建立方法及真伪鉴别方法 | |
Wu et al. | A printer forensics method using halftone dot arrangement model | |
CN107609482B (zh) | 一种基于汉字笔画特征的中文文本图像倒置判别方法 | |
CN104732548A (zh) | 一种基于纹理合成的打印文件鉴别方法 | |
CN109087234A (zh) | 一种文本图像中的水印嵌入方法和装置 | |
Benhamza et al. | Image forgery detection review | |
Wang et al. | Fourier-residual for printer identification | |
US7532756B2 (en) | Grayscale character dictionary generation apparatus | |
CN104463792A (zh) | 一种用于去除二值化图像中噪点的方法及其装置 | |
Chhabra et al. | Detecting fraudulent bank checks | |
Schreyer et al. | Intelligent Printing Technique Recognition and Photocopy Detection for Forensic Document Examination. | |
CN101854461B (zh) | 一种利用半色调信息检测文档真实性的打印文档取证方法 | |
JP5998090B2 (ja) | 画像照合装置、画像照合方法、画像照合プログラム | |
Lin et al. | Locating license plate based on edge features of intensity and saturation subimages |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |