规则碎纸片的拼接复原模型文档格式.docx
《规则碎纸片的拼接复原模型文档格式.docx》由会员分享，可在线阅读，更多相关《规则碎纸片的拼接复原模型文档格式.docx（18页珍藏版）》请在冰豆网上搜索。
3.假设纸面上除了文字外无其他杂点黑点。
4.假设切割过程中纸片没有损坏。
5.假设切割后的纸片不是歪斜的。
6.假设原图的左边（起始碎片的左边）与原图的右边（终始碎片的右边）边界线为白色。
7.假设原始图像本身的内容是连续的具有意义的。
8.假设所有的碎片中都存在相互匹配的碎片。
9.假设所有碎片都在一次扫描下完成。
三．符号说明
：
第i个碎片
纵向切割时，碎片的灰度矩
:
横纵切割时，碎片左边的灰度矩
横纵切割时，碎片右边的灰度矩
横纵切割时，碎片上边的灰度矩
横纵切割时，碎片下边的灰度矩
灰度矩
与
的相关系数
四．模型的建立与求解
5.1.1问题一的分析
如果是一个整体分裂成的两个碎片,那么两个碎片裂口对应的部分的颜色（对灰度图像而言,就是灰度信息）会有很强的相关性,这是碎片间除了轮廓外信息外最重要的信息之一。
由于碎片的轮廓是完全一致的，我们只考虑碎片裂口对应的部分颜色对图像匹配的影响。
完成图像碎片匹配过程后,相互匹配的图像碎片就被找到了,同时也得到了图像碎片的匹配段,并在碎片间建立了对应的匹配信息表,并记录了构成匹配对的匹配线段的相关信息如匹配起始点终止点等。
接下来的工作就是将相互匹配的图像碎片拼接起来,恢复图像的原貌。
5.1.2问题一的模型建立
设碎片
的每条边的灰度矩为
，
，i=1,2,…19,j=1或2，
=0或1。
i代表第i个碎片，j=1表示碎片左边，j=2表示碎片右边，
=0代表黑色，
=1代表白色。
由于每个碎片的左边只能与另一个碎片的右边相匹配，此时两向量之间的相关系数可表示为
采用相关系数这一测度可以测量两个描绘子之间的相似性,如果已知用于比较多个描绘子
，就可以计算出各个描绘子之间的相互关系,并找出最大相关系数
。
灰度矩相似性最高的两条边,就可以作为真正的匹配边,同时可以排除掉其他干扰匹配对。
5.1.3问题一的模型求解
首先通过碎纸片左右边界线是否为单一色（白色）找出起始碎片和终止碎片，将拥有与起始碎片右边相似度最高的边的碎片作为起始匹配图像，运用迭代法将剩余碎片匹配拼接。
利用matlab软件求得i从1到19时
相对应的
匹配图像，将复原后匹配碎片的顺序列入下表：
碎片
顺序
8
14
12
15
3
10
2
16
1
4
5
9
13
18
11
7
17
6
表一：
附件1的复原结果
表二：
附件2的复原结果
5.1.4问题一的模型结果的分析与检验
经过灰度矩的相关系数匹配，模型一最终精确的拼接复原了碎纸片。
在拼接过程中并没有进行人工干预，使得计算工作更为简便。
但是对于较为复杂的图像，纯粹地依靠计算机寻找匹配图像时，一般只能找到匹配图像的一个区域，并不能唯一确定具体的结果，从而增大误差。
5.2.1问题二的分析
问题二中不仅有纵向切割还有横向切割，碎片切割的复杂化，使得按照问题一的模型，不能让计算机自动地完成拼接过程。
在此文中，对这类边缘相同的碎纸片的拼接，理想的计算机拼接过程应与人工拼接过程类似，拼接时不需要考虑待拼接碎纸片边缘形状是否匹配，但要判断碎片内的字迹断线或碎片内的文字内容是否匹配，然而由于理论和技术的限制，让计算机具备类似人那种识别碎片边缘的字迹断线、以及理解碎片内文字图像含义的智能几乎不太可能。
但是利用现有的技术，完全可以获取碎片文字所在行的几何特征信息，比如文字行的行高、文字行的间距等信息，拼接碎片时如利用这些信息进行拼接，其拼接效率应该比较好。
由于此文字文档的文字行方向为单一水平方向，如果碎片内的文字行在碎片边缘断裂，那么与它相邻的碎纸片在边缘处一定有相同高度、相同间距的文字行，凭此特征可以很容易地从形状相似的多碎片中挑选出相邻碎片。
因文字行的高度特征、间距特征的识别比字迹断线识别和文字图像的理解实现起来要容易得多，利用碎片内文字行特征拼接形状相似的碎纸片理论上是可行的。
另一方面由于计算机数字分析图像能力的缺陷，让计算机对碎片进行完全意义上的自动化拼接也几乎不太可能，为保证拼接的准确性，需要在拼接过程中加入人工干扰过程。
一般而言拼接碎片时先利用计算机搜索与目标碎片匹配的未拼接碎片，并根据匹配程度按顺序显示待选碎片，我们再根据人脑进一步分析结果舍弃或拼接待选碎片。
这种半自动拼接方法综合利用了计算机高速计算能力以及人的文字图像识别和理解能力，拼接效率比纯人工高，拼接准确性也好于纯计算机拼接法。
5.2.2问题二的模型建立
为提高模型的准确性，假设汉字宽度与高度比值1/3~3。
这意味着每个文字图像与其他文字图像之间有空白点，文字图像宽度与高度的比值在1/3~3之间，如果碎片内有英文单词，应将英文单词图像拆分成类汉字图像，即将英文单词图像分割成宽度与高度近似的类汉字图像。
每个碎片像素高度H=180，宽度w=72。
碎片文字行的方向为水平方向，先将两个碎片按文字行方向（水平方向）线位置对齐,然后计算文字行方向线与碎片边界的交点与处于同一水平位置的另一个碎片交点的距离
如果碎片没有对齐，则距离相等的连续点的个数比对齐位置的距离相等连续点的个数要少，即距离相等的连续点的个数越多，图像越匹配。
由于各碎片边缘都为垂直线（或水平线），因此距离相等等价于相同高度下（或宽度）点的像素相同。
设X,Y,H,L分别为碎纸片的左右上下边的0-1灰度矩阵，函数F（
）表示对应
的对数，目标函数为maxF（
），maxF（
）。
5.2.3问题二的模型求解
（1）根据边界灰度矩阵
=0，
=0，从候选碎片中筛选出起始碎片
；
（2）计算该碎片上（或下）边界与未接碎片下边界（或上边界）可能拼接位置处的距离相等的连续点的个数；
（3）将距离相等的连续交点个数按升序排列，按升序依次显示未拼接碎片，人工选择碎片并拼接到计算机屏幕上；
（4）判断是否完成本列碎片的拼接，如果继续拼接，转
（2），否则转（5）；
=
（5）以第一列为起始碎片，进行11次下面步骤循环；
（6）计算该碎片左（或右）边界与未接碎片右边界（或左边界）可能拼接位置处的距离相等的连续点的个数；
（7）将距离相等的连续交点个数按升序排列，按升序依次显示未拼接碎片，人工选择碎片并拼接到计算机屏幕上；
（8）判断是否完成本行碎片的拼接，如果继续拼接，转（6），否则结束此次循环。
将复原后匹配碎片的顺序列入下表：
（注：
加黑序号为人工筛选序号）
128
159
82
199
135
73
160
203
169
134
39
31
51
107
115
176
94
34
84
183
90
47
121
42
124
144
77
112
149
97
136
164
127
58
43
125
182
109
197
184
110
187
66
106
150
21
173
157
181
204
139
145
29
64
111
201
92
180
48
37
75
55
44
206
104
98
172
171
59
208
138
158
126
68
175
45
174
137
53
56
93
153
70
166
32
198
71
156
83
132
200
80
33
202
133
170
205
85
152
165
27
60
89
146
102
154
114
40
151
207
155
140
185
108
117
101
113
194
119
123
49
54
65
143
186
57
192
178
118
190
95
22
129
28
91
188
141
61
19
78
67
69
99
162
96
131
79
63
116
163
72
177
20
52
36
168
100
76
62
142
30
41
23
147
191
179
120
86
195
26
87
38
148
46
161
24
35
81
189
122
103
130
193
88
167
25
105
74
表三：
附件3的复原结果
196
50
表四：
附件4的复原结果
5.2.4问题二的模型结果的分析与检验
该模型采用了基于边界交点距离的碎片半自动拼接算法，该算法不依赖于碎片几何特征，实现简单，可靠性比较好。
对于碎片文字是汉字的碎片，该模型能更精确地找出匹配碎片的图像，而对英文碎片，搜索的精度略低，人工筛选工作量稍大。
该半自动拼接算法最终求得的复原图像符合实际情况，且拼接计算工作量在允许范围内，效果显著。
5.3.1问题三的分析
较之问题一和问题二的碎片数据是单面打印文件，问题三研究的是在横纵切割的基础上加上正反面碎片的双面打印文件，显得更为复杂。
在碎片数量增加一倍的情况下，若单纯地以匹配条件去寻找匹配对象，则工作量会很大。
所以我们想到引进一种搜索算法，可以使搜索到的匹配对象范围更为精准，从而缩小工作量，尤其是人工筛选的工作量。
经过研究，决定采用广度搜索算法，在此基础上，又考虑到之前提取碎片边缘像素点，将之化为0-1灰度矩，有失图像的真实性，便决定提取碎片原有的像素作出灰度矩，而再用广度搜索算法搜索满足匹配条件的碎片对象，从而更为精准迅速地得出结果。
5.3.2问题三的模型建立
设X,Y,H,L为边界的实际灰度矩，d（
，d（
，目标函数为mind（
），mind（
5.3.3问题三的模型求解
对于一个给定的
，利用广度搜索算法求出minf（
）时的
，从而确定碎片
的最优匹配碎片
，再经由人工最终确定
的匹配图像。
将所求碎片复原序号填入下表：
136a
047b
020b
005b
152b
147b
143a
200a
086a
083b
039a
097b
090b
293a
162a
013b
024b
057b
035b
159b
073a
172b
122b
182a
105b
204a
141b
009a
145b
082a
054a
196a
112b
表五：
附件5碎片正面复原结果
表六：
附件5碎片反面复原结果
5.3.4问题三的模型结果的分析与检验
由于本模型中采用的灰度矩是对各边像素的真实提取，并且运用了广度搜索算法搜索最优匹配项，使得该模型的效率得到了提高，能在可接受的时间范围内求得匹配结果。
尽管已在问题二的基础上对模型进行了改进，仍不能脱离人工干预实现计算机的全自动化，在碎片数量多的情况下会有较大的工作量。
五．模型的改进与推广
本文对算法的计算工作量进行了分析，分析表明，对大多数碎片的实际碎纸片，拼接计算工作量在允许范围内，如果对算法作些改进，拼接计算量可大幅度减小。
本文虽然对二维图像碎片复原方法的关键技术进行了研究,但依然还有很多地方需要进一步深入探讨。
本文主要针对规则碎片的特征提出了基于边缘特征点的复原算法,对特征点的利用可以推广到任意不规则碎片类型,但还需进一步研究。
另外,由于本文只是针对灰度图像进行了实际研究并取得了较好效果,理论上讲,彩色图像和灰度图像适用于同样的方法,但彩色图像包含了更丰富的内容信息,无论是算法研究还是在实证研究方面都还有很大的研究空间,有更强的实用性,也具有更大的挑战性。
本文主要针对二维规则碎片,应用前景受限,如果能把二维的工作扩展到三维空间,将会有更加广泛的应用空间。
将规则碎片推广到针对多边形碎片复原的复原匹配算法，由于多边形碎片的特殊性,综合利用了多边形碎片的形状和灰度特征进行匹配,这是整个算法的关键,也是创新点。
图像碎片复原问题除了上边提到的以外,还可以推广到很多其他领域,比如,破碎的艺术品复原,破碎的纸币复原,公安机关破案中的破碎证物复原等等。
从理论上讲,图像碎片复原技术的核心技术是碎片匹配拼接问题,是模式识别技术的一个新发展,超出了传统模板匹配的局限,对这个问题的深入研究将促进机器识别,计算机视觉等相关领域的发展,开辟模式识别新的应用领域。
因此,不管从理论角度还是从实践角度来看,本课题的相关研究工作都具有积极的意义。
六．模型的优缺点分析与评价
虽然我们进行实验仿真所涉及的纸片数量不多，但是它充分表明了所进行的碎纸拼接复原技术研究的可行性。
虽然文中所用算法在实验中取得了比较好的效果,但还需挖掘潜力,进一步优化。
因此，如何减少错误匹配的数量，提高效率，增加适应性，也是下一步的工作重点。
在实际应用中，由于设备等条件的限制，一次处理往往不能完全完成拼接任务。
因此，必须开发合理的数据库，以充分利用每次拼接得到的结果。
虽然在进行简单拼接复原的时候，计算机可以满足要求，但是在实际应用中，通常会涉及对大量纸片数据进行管理和处理工作。
目前，我们主要还是以实现功能为主，下一步必须优化程序结构，改善用户界面，提高程序的交互能力，真正实现快速有效的计算机辅助碎纸
自动拼接复原系统。
本文对碎纸自动拼接进行了研究，在前人研究的基础上寻找一些新的研究方法进行研究，或者是对一些传统的方法进行应用与改进，可进行深入探讨，将其应用于计算机辅助碎纸自动拼接中，以取得很好的效果。
参考文献
[1]曹弋，MATLAB教程与实训，北京，机械工业出版社，2012年
[2]盛骤，概率论与数理统计第四版，北京，高等教育出版社，2008年
[3]邓薇，MATLAB函数速查手册，北京，人民邮电出版社，2010年
[4]姚泽清，全国大学生数学建模竞赛赛题与优秀论文评析（2005年-2011年B题），北京，国防工业出版社，2012年
[5]贾海燕，碎纸自动拼接关键技术研究，2005年
[6]先毅，图像碎片复原方法的研究，2010年
[7]楼兵军，基于区域的图像匹配算法的关键技术研究，2006年
附件
附录一：
附件1的复原图
附件2的复原图
附录二：
附件3的复原图
附件4的复原图
附录三：
问题一的程序代码
clc;
clear
A=cell（1,19）;
fori=0:
18;
ifi<
=9
NO=strcat（'
00'
num2str（i））;
elseifi<
=99
0'
elseNO=num2str（i）;
end
dir=[NO,'
.bmp'
];
a=imread（dir,'
bmp'
）;
end
fork=1:
19
B{1,k}=im2bw（A{1,k}）;
end
%fork=1:
%imshow（B{1,k},[]）
%end
C=cell（1,19）;
M=zeros（1980,19）;
N=zeros（1980,19）;
fori=1:
19;
m=find（B{1,i}（:
1）==0）;
ifsum（m）==0;
c=i;
M（1:
end,i）=B{1,i}（:
1）;
n=find（B{1,i}（:
72）==0）;
ifsum（n）==0
d=i;
N（1:
end,i）=B
升级会员 