全国大学生数学建模竞赛—B题—碎纸片的拼接复原

1、2013高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了全国大学生数学建模竞赛章程和全国大学生数学建模竞赛参赛规则（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛

2、规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： B 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 西华大学 参赛队员 (打印并签名) ：1. 杨尚安 2. 刘 洋 3. 叶 军 指导教师或指导教师组负责人 (打印并签名)： （论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。） 日期：

3、 2013 年 09 月 15 日赛区评阅编号（由赛区组委会评阅前进行编号）：2013高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：碎纸片的拼接复原摘要本文通过分析题中相关要求及条件，建立数学模型解决了各种规则碎纸片的拼接复原问题。针对问题一， 首先将题中所给图片导入matlab软件，利用imread函数得到每张图片的文字灰度像素矩阵，再取出所有矩阵左、右列，建立像素绝对差拟配模型，得到拟配程度最高的两幅图片，进

4、行拼接，出现不合理拼接情况则进行人工干预，最后重复上述过程，完成全部拼接并导出图像。针对问题二，首先将全部碎片导入matlab软件，经过处理得到每张碎片中符号距离碎片上下端的像素位，再根据分类聚类思想，利用excel表格处理，将所有具有“相同”像素位的图片分为一组，得到11个分组，然后在每一个分组中建立左右连接点数目最匹配模型，再配合人工干预，将所有碎片拼接为一行图像，最后将这11行图像利用问题一中模型拼接为最终图像并打印结果。针对问题三，首先建立一种基于K-Means局部最优性的高效聚类模型，然后根据模型利用matlab，将所给图片全部导入分类，分好类并人工调整补充后再利用matlab在每一

5、组分类中利用问题二模型在人工干预情况下得出原始图像并打印结果。关键词：像素绝对差拟配模型 左右连接点数目最匹配模型 人工干预一、问题重述破碎文件的拼接在司法物证复原、历史文献修复以及军事情报获取等领域都有着重要的应用。传统上，拼接复原工作需由人工完成，准确率较高，但效率很低。特别是当碎片数量巨大，人工拼接很难在短时间内完成任务。随着计算机技术的发展，人们试图开发碎纸片的自动拼接技术，以提高拼接复原效率。请讨论以下问题：1. 对于给定的来自同一页印刷文字文件的碎纸机破碎纸片（仅纵切），建立碎纸片拼接复原模型和算法，并针对附件1、附件2给出的中、英文各一页文件的碎片数据进行拼接复原。如果复原过程需

6、要人工干预，请写出干预方式及干预的时间节点。复原结果以图片形式及表格形式表达（见【结果表达格式说明】）。2. 对于碎纸机既纵切又横切的情形，请设计碎纸片拼接复原模型和算法，并针对附件3、附件4给出的中、英文各一页文件的碎片数据进行拼接复原。如果复原过程需要人工干预，请写出干预方式及干预的时间节点。复原结果表达要求同上。3. 上述所给碎片数据均为单面打印文件，从现实情形出发，还可能有双面打印文件的碎纸片拼接复原问题需要解决。附件5给出的是一页英文印刷文字双面打印文件的碎片数据。请尝试设计相应的碎纸片拼接复原模型与算法，并就附件5的碎片数据给出拼接复原结果，结果表达要求同上。二、模型假设1、假设全

7、部碎纸片边缘光滑2、假设字符色调一致3、假设字符间距相同，没有特殊情况4、假设除字符外，页面没有其他地方具有任何色彩5、假设英文字符书写标准，大小写字号均相同三、符号说明表示灰色像素矩阵表示灰色像素矩阵的列数表示灰色像素矩阵的行数表示第几个碎片表示某个像素点表示某灰度像素点为黑色还是白色表示灰色像素矩阵最右边列表示灰色像素矩阵最左边列表示某个碎片灰色像素矩阵最左列与另一个碎片灰色像素矩阵最右列的差的绝对值的和四、模型建立与求解4.1问题一问题分析整体来看，本问题要求利用数学模型，改原有手动拼接技术为自动或半自动拼接技术，完成题中所给的相应碎纸片的拼接复原工作。具体操作，考虑所给碎纸片内容仅有汉

8、字或英文，而没有颜色、大小、字形之分。因此，只能利用碎纸片中相应的文字特征进行操作，考虑碎纸片扫描进入在计算机后是以图片的形式存在，而图片又是以像素的情况组成。所以，首先可将图片导入matlab中，以其像素为基点，得到每个图片的像素矩阵，每一像素矩阵即可表示该图片的特征。为了利用图片像素矩阵完成图片的拼接，考虑问题一只是将原图分为了19列，每一列具有1980像素，首先可根据左端全为空白，找出原图最左一列碎片，然后利用拼接好的图片最右列像素点去匹配未拼接图片的最左列像素点，使得拼接最为吻合的即为需要拼接的图片，然后拼接，再重复上述过程，直到拼接完成。具体操作流程如下：碎纸片Matlab像素矩阵求

9、解程序取出最左与最右列拼接导出结果判断吻合程度导入吻合不吻合像素绝对差拟配模型图1 问题一解答流程图数据处理将图片导入matlab中，然后编写程序（具体代码见附录1），可得每个碎纸片灰度像素矩阵（碎片000局部像素点如下）。图2 碎片000局部灰度像素点列像素绝对差拟配模型建立令碎片导入matlab编程计算所得的灰色像素矩阵为：由于碎片像素为72*1980，因此矩阵也是72*1980的，矩阵每一列数据即为碎片相应列像素值，其中每个像素点表示此处为黑色或白色，用表示某灰度像素点为黑色还是白色，即：令表示灰色像素矩阵最右边列，那么令表示灰色像素矩阵最左边列，则令表示某个碎片灰色像素矩阵最左列与另一

10、个碎片灰色像素矩阵最右列的差的绝对值的和。那么有根据上述模型即可确定某一碎片灰度像素矩阵最右边列与其余未拼接碎片最左边列的绝对差值，下面讨论因差值不同而产生的匹配问题。1、最左列的确定：当出现某一碎片灰度像素矩阵最左列均为255时，那么说明该碎片为原始图像的最左列。2、假设出现情况，那么首先将对应的碎片与该基准碎片进行拼接，若拼接不合适，这时就需要人工干预，换对应的碎片与基准碎片进行拼接。情况如下：这是不确定的，而进行人工干预选择对应的碎片后，将会出现下面情况：这样就能正确的完成两个碎片的拼接。3、假设出现情况，这与上述情况相同。因此，人工干预方式及时间选择也相同。4.1.4像素绝对差拟配模型

11、求解对于附件一中碎片复原，根据上述模型，利用matlab软件，求解可得008碎片最左端矩阵列与006碎片最右端矩阵列均为：，因此，可知008碎片为复原图最左一个碎片，006碎片为复原图最右端碎片。其余求得所有最小的距离的值，根据的值，可将碎片进行复原。复原结果如下表，复原图像见附录2。008014012015003010002016001004005009013018011007017000006表1 问题一中文复原表格序列对于附录二英文复原，与上求解过程雷同，利用matlab可得复原结果如下表，复原图像见附录3。0030060020070150180110000050010090130100

12、08012014017016004表2 问题一英文复原表格序列4.1.5问题一综合分析 综上所述，对于问题一的求解过程，未使用人工干预。本文除使用对问题所给的碎片进行复原外，同时对具有相同属性的其他图形碎片也进行了复原，效果良好，模型稳定，可推广到所有只进行竖切的文档恢复。4.2问题二中文碎片复原问题分析综合分析。由于考虑问题二在问题一的基础上将碎片分的更加的细小，那么碎片的灰色像素矩阵数据在原有的基础上将会变得少很多，考虑使用问题一方法及模型，那么首先就要构造出与问题一相同的19个竖碎片，因此考虑将所有碎片分为19组，但经过试验分为19组后，由于空白出现太多，在每组中将11个碎片拼接在一起是

13、相当困难的。因此，转变思想，考虑将所给所有碎片分为11个组，在每个分组中将19张碎片拼接在一起，然后在将11个分组拼接在一起完成最后解答。具体操作。要想将11*19张图片分为11组，考虑文字具有行高的性质，分组中所拼接的19张碎片，所有文字具有的行高应该都是相同的。根据这一思想，可将所有碎片导入matlab中，编程计算可得每张碎片符号距离碎片上下端的像素位，并将所有结果导入excel中，然后根据分类与聚类思想，利用excel表格处理，将碎片符号距离碎片上下端的像素位“相同”（不是绝对相等，允许误差前后波动两个像素）的点分为一组，对于出现空白位置误差较大的点可根据单边距离进行分类与聚类，若根据单

14、边无法确定具体分入那组，那么就同时分入可能的分组中。分组完成后那么每个分组中的图片定能拼接为一行图片，那么我们可建立左右连接点数目最匹配模型，结合人工干预，将每个分组中图片拼接在一起。最后利用问题一中模型可将11个分组拼接在一起得到原图。具体流程如下图：人工干预拼接为一行209张碎片导入matlab计算字符距离上下边界的像素位分类与聚类思想得到11个分组取任一分组左右连接点数目最匹配模型像素绝对差拟配模型完成拼接，导出图像图3 问题二解答流程图数据处理将209张碎片导入matlab中，编程得到每张碎片灰色像素矩阵，然后在利用矩阵编写程序得到每张碎片字符距离上下边界的像素位，并将其导入excel

15、中（具体代码见附录4），下表为000至016结果：（其余碎片结果见附录5）碎片编号012345678910111213141516上距离022372512000000392589316下距离150005809159600000030表3 碎片字符距离上下边界的像素位情况表得到像素位上下边缘距离后可根据上下距离“相等”（不是绝对相等，允许误差前后波动两个像素）原则，利用excel表格处理将所给数据分为11组。其中距上边缘距离为0，下边缘距离为21的一个分组为下表：（其余分组见附录6）34424347587784949718390136112127124121144149164表4 某一属于同行碎

16、片的分组情况表在每一分组内，再利用matlab编程计算每张碎片左端与右端具有的可连接点数目（采用四舍五入原则） (具体代码见附录7) ，下表为上一分组数据的左右连接点数目：（其他分组连接点数目见附录8）4(34)73(42)64(43)06(47)35(58)34(77)27(84)00(94)49(97)31(183)99(90)63(136)33(112)75(127)66(124)53(121)35(144)47(149)93(164)5表5 某一同行碎片的左右连接点数目情况表 左右连接点数目最匹配模型本模型属于半自动模型，需人工干预，具体步骤如下：1、选取任一分组左右连接点数目情况表，

17、观察左右连接点数；2、选取左端连接点数目为0的碎片作为最左端碎片，并将该图片作为基准图片；3、观察基准图片右端连接点数目，从未拼接图片左端连接点数目中找寻与该数目最接近的碎片，人工控制，观察是否可连接。若可连接则拼接上，并将新拼接上碎片作为基本图片，若不可连接，则重新找寻符合要求的碎片，观察是否可连接；4、重复3步骤，直到将图片全部连接完成。模型求解以上述模型为标准，考虑数据处理中那行连接过程。首先，寻找19个点钟左端连接数为0的点，找到（94）号碎片，将其作为基准图片，观察其右边连接点数为4，从其余碎片中找寻发现（34）（43）（77）左端连接数均为4，因此，通过人工干预，观察图片字样走势发

18、现只有（34）号碎片符合要求，再将（34）号图片作为基准图片，其右端连接点数为7，从未连接碎片中找寻发现（84）（149）号均为7 ，同理（84）碎片作为基准图片，以此类推即可得到该分组图片排序为：（94）（34）（84）（183）（90）（47）（121）（42）（124）（144）（77）（112）（149）（97）（136）（164）（127）（58）（43）其具体碎片拼接图形如下：然后根据上述模型，以相同的办法结合附录6中分组情况即可将全部11个分组中图片的连接情况找出，然后利用问题一中像素绝对差拟配模型即可拼接处原图，得到原图表格连接情况如下，具体图像见附录9。04905406514

19、318600205719217811819009501102212902809118814106101907806706909916209613107906311616307200617702005203616810007606214203004102314719105017912008619502600108701803814804616102403508118912210313019308816702500800910507407115608313220001708003320219801513317020508515216502706001412800315908219913501207

20、316020316913403903105110711517609403408418309004712104212414407711214909713616412705804312501318210919701618411018706610615002117315718120413914502906411120100509218004803707505504420601010409817217105900720813815812606817504517400013705305609315307016603219608914610215411404015120715514018510811700

21、4101113194119123表6 汉字碎片拼接情况表4.2.5中文碎片恢复综合分析由于解决本问题使用的左右连接点数目最匹配模型，属于半自动模型。因此，对本文的恢复进行了人工干预。恢复此中文文档，本模型一共进行了9次人工干预。干预方式为：终止程序继续运行，将程序拼接过程恢复至上一步（出现碎片拼接不吻合时的前一步），然后将程序用于拼接的碎片导出，再恢复程序继续运行，找到该步拼接吻合碎片并拼接后，再将导出碎片重新导入继续运行程序。干预时间节点:干预时间节点即对每行碎片单独拼接时，出现碎片拼接不吻合情况时的节点。4.3问题二英文碎片复原4.3.1问题分析对于附录四英语碎片恢复，由于英文与汉字写法不

22、同，英语中弧线居多，而汉字中直线居多。因此，可以采用另一种方式对英文碎片进行拼接，依然考虑问题一中的像素绝对差拟配模型，可首先任意选择一张基础碎片，然后利用该模型进行适应性匹配，匹配过程中加以人工干预。具体操作流程如下：全部拼接完成某一英文碎片像素绝对值拟配模型人工干预拼接完成全部拼接成功寻找下一个图4 英文碎片复原流程图模型建立1.1像素绝对值拟配模型令碎片导入matlab编程计算所得的灰色像素矩阵为：由于碎片像素为72*180，因此矩阵也是72*180的，矩阵每一列数据即为碎片相应列像素值，其中每个像素点表示此处为黑色或白色。令表示灰色像素矩阵最右边列，那么令表示灰色像素矩阵最左边列，则令

23、表示某个碎片灰色像素矩阵最左列与另一个碎片灰色像素矩阵最右列的差的绝对值的和。那么有根据上述模型即可确定某一碎片灰度像素矩阵最右边列与其余未拼接碎片最左边列的绝对差值，下面讨论因差值不同而产生的匹配问题。1、最左列的确定：当出现某一碎片灰度像素矩阵最左列均为255时，那么说明该碎片为原始图像的最左列。2、假设出现情况，那么首先将对应的碎片与该基准碎片进行拼接，若拼接不合适，这时就需要人工干预，换对应的碎片与基准碎片进行拼接。情况如下：这是不确定的，而进行人工干预后将会出现下面情况：这样就能正确的完成两个碎片的拼接。3、假设出现情况，这与上述情况相同。因此，人工干预方式及时间选择也相同。1.2人

24、工干预在进行像素绝对值拟配模型计算后，将会得到与基准碎片拼接度最大的几个碎片，然后利用这几个碎片可进行人工干预，具体人工干预模型如下：1、首先将程序计算得到的拟配程度最大的碎片与基准碎片进行拼接；2、人工判断拼接是否合理；3、若拼接合理则进行下一次拟配模型计算，若拼接不合理则找寻第一步中与基准碎片拟配差一点的碎片进行拼接；4、直到找到拼接成功的点才结束本次拼接，并将新拼接上的图片作为基本图片利用模型寻找拟配度最高的碎片，返回第一步。4.3.3模型求解结合上述模型1.1及1.2可计算得到问题二英文碎片复原图表格如下，具体图像见附录101910750111541901840021041800641

25、060041490322040650390671472011481701961980941131640781030910801010261000060170281460860511070290401581860980241171500050590580920300370461270191940931410881211261051551141761821510220572020711650821591390011290631381530530381231201750850501601870972030310200411081161360730362071350150760431990451730

26、79161179143208021007049061119033142168062169054192133118189162197112070084060014068174137195008049172156096023099122090185109132181095069167163166188111144206003130034013110025027178171042066205010157074145083134055018056035016009183152044081077128200131052125140193087089048072012177124000102115表7 问

27、题二英文碎片复原图表格表示4.2.4英文碎片恢复综合分析由于英文碎片相似程度高于中文图片。所以每一次以基准图片找寻最佳匹配图形时很多时候出现多张图片符合匹配，因此，对此英文碎片的恢复进行了人工干预。恢复此英文文档，本模型一共进行了39次人工干预。干预方式为：终止程序继续运行，将程序拼接过程恢复至上一步（出现碎片拼接不吻合时的前一步），然后将该步程序用于拼接的碎片导出，再恢复程序继续运行，直到找到与该基准碎片拼接吻合碎片并拼接完成。干预时间节点:当出现与基准碎片匹配不吻合时。4.3问题三问题分析考虑问题三附录中所给图片具有正反面，却不知每一个序号中a是正面还是b是正面，这也真是问题二英语复原与问

28、题三双面复原的区别。因此，问题二中所用的分类与聚类的方法不能完成分组。为了完成分组，我们可考虑使用一种更加严密，严苛的分类方法，只要分类完成，那么再使用问题二连接图片的办法即可实现图片的复原。4.3.2模型建立许多聚类算法的基本框架是搜索与合并。如在层次方法中需要搜索两个距离最近的类簇然后合并；而基于密度的聚类算法则不断地搜索高密子区域，然后利用连通性将其合并到当前聚类结果中。很明显，搜索过程需要面对整个样本集合，通常会导致算法低效。如DBSCAN需要测试每个对象是否是核心对象，并对每个核心对象搜索其直接密度可达的对象，如果没有空间索引的辅助，DBSCAN算法的复杂度为。实际上，现有的很多聚类

29、算法已经关注到这个问题，如CURE算法利用采样方法减小搜索空间，而Chameleon算法则通过图划分算法将样本对象聚类为大量相对较小的子簇。具体到本文，我们采用了随机采样和K-Means算法高效地聚出大量的高密子簇，后续处理都基于这些构造出的高密子簇进行，无须直接面对所有的样本。我们称这种算法框架为构造与合并。因此，K-MeanSCAN算法的处理流程为：随机采样、预聚类、合并和后处理。聚类算法的本质是密度估计问题，K-MeanSCAN算法的核心思想是，增加K-Means算法中高斯混合模型的高斯分量数目以提升密度估计的精度，并利用K-Means聚类的局部最优和结果敏感性的特征进行高斯分量(即高密

30、子簇)的合并剪枝，克服过拟合问题，最终实现有效的聚类。K-Means 算法实质上是一种将聚类视为密度估计问题的概率方法。在概率方法中，假设样本来自于如下形式的混合模型:式中，是待估计的参数向量；条件概率密度，称为分量密度，表示类别的概率密度形式，且参数向量未知；先验概率称为混合因子。为了简化问题，K-Means算法进一步假设:(1) 每个类别的概率密度形式为球形高斯分布，即且，未知，（2）每个样本唯一地属于一个类别；(3) 假设所有类别的混合因子相等。于是,混合模型简化为该简化模型可以通过最大似然方法求解,对于观测样本X=(x1,xn),相应的对数似然函数为：最大化该对数似然函数等价于最小化上

31、式的欧氏距离平方项,即得到 K-Means 的误差平方和准则:通过迭代优化上述的误差平方和准则,K-Means算法最终可以估计出每个高斯分量的均值向量和协方差矩阵。式中,是类簇的样本数目。要获得更有效的样本密度估计,的值自然是越小越好.但是, 的值不仅取决于样本的分类情况,而且与类别数目有关.当类别数目给定时, 的值由样本的分类情况所决定,且存在一个最小值对应于最优的样本类别划分.如果类别数目和高斯混合模型假设与实际问题相匹配时,最小值必定很小,从而可以很好地近似样本密度;而如果模型假设不合理,则最小值可能依然很大,对样本分布的近似效果较差。对于任意形状的类簇,很明显不能直接要求数据分布满足高

32、斯混合模型的假设,否则会导致最小的误差平方和很大.实际上,高斯混合模型具有很强的表达能力,如果高斯分量密度的数目足够大,则高斯混合模型几乎可以近似任意一种概率分布。换言之,随着类别数目的增加,相应的会减小。简单证明如下: 类别数目的增加,必然导致最终每个类簇的形状缩小,对应于高斯分量的减小；而从公式(3)和公式(4)可以得出,即最小的误差平方和正比于；因此, 的增加会导致最小的误差平方和的减小。 极端情况下, ,则每个样本点都是一个类簇,即,说明此时的经验误差为0,但是此时,模型的推广能力极差。根据统计学习理论,经验误差最小并不等于期望误差最小,经验风险只有在样本数无穷大时才趋近于期望风险.因

33、此,经验误差最小不能保证分类器的推广能力,需要找到经验风险最小和推广能力最大的平衡点。同样,利用足够多的高斯分量组成的混合模型来描述数据会导致过拟合的问题,影响模型的推广能力。因此,在K-MeanSCAN算法中,我们采用过拟合-剪枝的策略进行聚类,即首先使用分量足够多的高斯混合模型来较好地近似样本分布,然后通过合并一些高斯分量的剪枝策略来处理过拟合问题。4.3.3模型求解根据上述K-Means局部最优性的高效聚类模型，利用matlab编程计算后结合人为补充调整可得图片分组情况。其中一种分组情况见下表，其余分组情况见附录11。108a110a125a066a047a020a150a029a140

34、a155a111b136b174b183b164b189b018b081b078b表8 一种碎片分组情况4.3.4最终求解上述模型求解结束后，结合问题二英文求解模型及人工干预时机，可将原图还原，其具体还原图像见附录12，具体还原图像表格如下：12345678910078b111b125a140a155a150a183b174b110a066a089a010b036a076b178a044a025b192a124b022a186b153a084b042b030a038a121a098a094b061b199b011b161a169b194b173b206b156a034a181b088b107a

35、149b180a037b191a065b115b166b001b114a184b179b116b207a058a158a197a154b028b146a171b031a201a050a190b092b019b016b177b165b195a128a157a168a046a067a063b075b167a003b007b085b148b077a004a069a032a074b126b023b133a048a051b095a160b119a033b071b052a099a043a096b109a123a006a104a134a113a026b111213141516171819108a018b02

36、9a189b081b164b020a047a136b120b144a079a014a059a060b147a152a005a137b045a138a056b131b187b086b200b143b198b087a132b093a072b175a097a039b083a151b170b041a070b139b002a162b203b090a012a017b102b064b208a142a057a024a013a053b202a021b130a163a193b073b159a035a117b008b068b188a127a040a182b122a172a176a185a000b080b027a13

37、5b141a204b105a062a129b118b101a015b205a082b145a009b049b091a106b100b055b103a112a196b054b表9 问题3其中一面恢复图表格表示12345678910136a047b020b164a081a189a029b018a108b066b005b152b147b060a059b014b079b144b120a022b143a200a086a187a131a056a138b045b137a061a083b039a097b175b072a093b132a087b198a181a090b203a162a002b139a070a04

38、1b170a151a001a013b024b057b142b208b064a102a017a012b028a035b159b073a193a163b130b021a202b053a177a172b122b182b040b127b188b068a008a117a167b105b204a141b135a027b080a000a185b176b126a009a145b082a205b015a101b118a129a062b052b054a196a112b103b055a100a106a091b049a026a111213141516171819110b174a183a150b155b140b125b

39、111a078a124a192b025a044b178b076a036b010a089b094a098b121b038b030b042a084a153b186a034b156b206a173a194a169a161b011a199a166a115a065a191b037a180b149a107b088a154a197b158b058b207b116a179a184a114b016a019a092a190a050b201b031b171a146b075a063a067b046b168b157b128b195b165a074a032b069b004b077b148a085a007a003b071a

40、033a119b160a095b051a048b133b023a113b134b104b006b123b109b096a043b099b表10问题3另一面恢复图表格表示4.3.5问题三综合分析对于本问题，由于英语碎片出现正反两面，分类时情况复杂。因此，恢复原文时进行了多次人工干预。干预方式：对正反英文碎片进行分类结束后，对未归类的碎片，采用人工归类方法，进行多次干预。在每一类中进行碎片拼接时，出现拼接不匹配时，终止程序继续运行，将程序拼接过程恢复至上一步（出现碎片拼接不吻合时的前一步），然后将程序用于拼接的碎片导出，再恢复程序继续运行，找到该步拼接吻合碎片并拼接后，再将导出碎片重新导入继续运行

41、程序。干预时间节点：分类完成后，对未进行归类的碎片进行人工分类时，还有当进行拼接出现拼接碎片不吻合时。五、模型评价与推广5.1模型的优点使用该模型，可以很好的提高工作效率，将原来使用纯手工拼接过程变得简单合理，而且拼接过程全部在计算机上进行，使得结果更加准确，从碎片出发复原原图像，以点代面的做法，使思路更加简单明了；运用matlab计算及绘出图形，用数形结合的方法来进行分析，模型思路更加清晰，更有说服力。5.2模型的缺点模型求解过程中，很多地方使用了人工干预，使得整个过程不能实现完全的智能化。有时的四舍五入也可能是结果有些偏差。5.3模型的推广该题所建模型求解的是与相似图形的匹配问题，在实际生

42、活中，与图像相似，匹配，识别相关的问题，均可以使用该模型进行运算。对于本题问题三，我们也可从另一方面入手考虑，由于每个碎片为180*72像素的，所含信息量较少，因此，使用问题一的方法直接求解不行。如果我们在保证图片大小不变的情况下，增加图片的像素或者寻求另一种模式使碎片分组更加细，那么直接运用方法一求解规则碎片复原将变得相当的容易。六、参考文献1 2 七、附录附录1，使用软件名称：matlab(daima1)A=imread('附件1000.bmp');%区图片的相对路径附录2，使用软件名称：matlab(daima2)function Pho1clcclear allA=im

43、read('附件1000.bmp');%区图片的相对路径m,n=size(A);pho=ones(m,1)*255;sort=zeros(1,19);x=0;left=zeros(m,1);right=zeros(m,1);local=0;Used=zeros(19,1);B=zeros(1,19);for i=1:19 %首先找到最左两边的碎片 if i<11 str=strcat('附件100',num2str(i-1),'.bmp'); else str=strcat('附件10',num2str(i-1),'

44、.bmp'); end a=imread(str); left=a(:,1);% 取第一列 if length(find(left>250)>1700 %当第一列所有数据为白色时，就定义为这个碎片在图片的最左边 pho=cat(2,pho,a); Used(i)=1; x=x+1; %图片的顺序 sort(x)=i-1; right=a(:,n);% 取最后一列 break; endendfor i=1:18 min=999999999999999999; for k=1:19 if Used(k)=1 if k<11 str=strcat('附件100

45、9;,num2str(k-1),'.bmp');%拼接图片路径 else str=strcat('附件10',num2str(k-1),'.bmp'); end a=imread(str); left=a(:,1);% 取第一列 distance=abs(right-left); temp=sum(sum(distance); if temp<min max=temp; local=k; end end if Used(local)=1 && k=19 if local<11 str1=strcat('附件10

46、0',num2str(local-1),'.bmp'); else str1=strcat('附件10',num2str(local-1),'.bmp'); end goal=imread(str1); pho=cat(2,pho,goal); %将得到的碎片加入到图片中 right=goal(:,n); %取右边 Used(local)=1; %将这个碎片设置为已经使用 x=x+1; %图片的顺序 sort(x)=local-1; end endendpho(:,1)=;imshow(pho)imwrite(pho,'附件1tr

47、ue1.bmp','bmp')sortend附件3，使用软件名称：matlab(daima3)%-附件一拼接-%function Pho2A=imread('附件2000.bmp');%取图片的相对路径m,n=size(A);pho=ones(m,1)*255;sort=zeros(1,19);x=0;left=zeros(m,1);right=zeros(m,1);local=0;Used=zeros(19,1);for i=1:19 %首先找到最左两边的碎片 if i<11 str=strcat('附件200',num2str(i-1),'.bmp'); else str=strcat('附件20',num2str(i-1),'.bmp'); end a=imread(str); left=a(:,1); % 取第一列 if length(find(left>