关于碎纸片自动拼接的数学模型大学生数学建模竞赛优秀论文

《关于碎纸片自动拼接的数学模型大学生数学建模竞赛优秀论文》由会员分享，可在线阅读，更多相关《关于碎纸片自动拼接的数学模型大学生数学建模竞赛优秀论文（29页珍藏版）》请在装配图网上搜索。

1、关于碎纸片自动拼接的数学模型摘要本文针对生活中破碎文件的拼接难度大，效率低等现象，从题目所给的情形出发，利用计算机软件把碎纸片图像转化为数字图像，综合运用matlab软件中的数字图像处理方法，建立了以图与图之间的相似程度为基准的数学模型。这个模型的评价标准很简单，就是相似度函数的值。通过比较图像与图像之间的相似度函数的值的大小，就可以得出碎纸片的具体拼接序列。对于问题（1），首先，用matlab软件的imread函数对图像的进行读取，得到数据矩阵为。其次，根据模型的假设（1），找到最右端的碎纸片，并记为。然后，以数据矩阵为基础，引入相似度函数，并求 出相似度函数值。最后，用matlab工具箱中

2、的sort函数把所得到的相似度函数值进行排序，所得到的相似度函数值最小的图像即为与最右端的碎纸片匹配的图像。如此重复18次，即可得附件1的中文图像的排列序号，结果如表1所示。同理可得附件2的英文图像排列序号，结果如表2所示。复原结果图片见论文附件的图1和图2。对于问题（2），同样先找到最右端的11张图像和最上方的19张图像，根据图像的页边距特性确定原图像右上角的第1张图像。利用问题（1）的算法可得最右端的11张图像和最上方的19张图像的排列序号。然后，在问题（1）的算法的基础上，利用图像中的文字的固定间距去改进算法，缩小搜索范围，并在拼接完一行后显示一次结果，由于近似距离计算公式与人主观视觉差

3、异，所以需要人机交互调整结果。如此重复18次，即可得附件3的中文图像的排列序号，结果如表3所示。同理可得附件4的英文图像排列序号，结果如表3所示。对于问题（3），与问题（2）相似，只是碎纸片由单面变为双面。因此在匹配图像时，引入两重相似度函数，以确保正反两面能同时匹配。同时每匹配5张图像显示一次结果，以增加人工干预次数。如此重复若干次，即可得最终的复原图像。关键字 相似度函数 matlab软件 数字图像处理 29一、 问题的重述破碎文件的拼接在司法物证复原、历史文献修复以及军事情报获取等领域都有着重要的应用。传统上，拼接复原工作需由人工完成，准确率较高，但效率很低。特别是当碎片数量巨大，人工拼

4、接很难在短时间内完成任务。随着计算机技术的发展，人们试图开发碎纸片的自动拼接技术，以提高拼接复原效率。请讨论以下问题：1. 对于给定的来自同一页印刷文字文件的碎纸机破碎纸片（仅纵切），建立碎纸片拼接复原模型和算法，并针对附件1、附件2给出的中、英文各一页文件的碎片数据进行拼接复原。如果复原过程需要人工干预，请写出干预方式及干预的时间节点。复原结果以图片形式及表格形式表达（见【结果表达格式说明】）。2. 对于碎纸机既纵切又横切的情形，请设计碎纸片拼接复原模型和算法，并针对附件3、附件4给出的中、英文各一页文件的碎片数据进行拼接复原。如果复原过程需要人工干预，请写出干预方式及干预的时间节点。复原结

5、果表达要求同上。3. 上述所给碎片数据均为单面打印文件，从现实情形出发，还可能有双面打印文件的碎纸片拼接复原问题需要解决。附件5给出的是一页英文印刷文字双面打印文件的碎片数据。请尝试设计相应的碎纸片拼接复原模型与算法，并就附件5的碎片数据给出拼接复原结果，结果表达要求同上。【数据文件说明】(1) 每一附件为同一页纸的碎片数据。(2) 附件1、附件2为纵切碎片数据，每页纸被切为19条碎片。(3) 附件3、附件4为纵横切碎片数据，每页纸被切为1119个碎片。(4) 附件5为纵横切碎片数据，每页纸被切为1119个碎片，每个碎片有正反两面。该附件中每一碎片对应两个文件，共有21119个文件，例如，第一

6、个碎片的两面分别对应文件000a、000b。【结果表达格式说明】复原图片放入附录中，表格表达格式如下：(1) 附件1、附件2的结果：将碎片序号按复原后顺序填入119的表格；(2) 附件3、附件4的结果：将碎片序号按复原后顺序填入1119的表格；(3) 附件5的结果：将碎片序号按复原后顺序填入两个1119的表格；(4) 不能确定复原位置的碎片，可不填入上述表格，单独列表。二、 问题的分析碎纸，即一张纸在外力的作用下被分开的几个小分块。而在实际生活中，往往需要我们把这几个碎纸块还原成一张纸，这就需要用到碎纸拼接技术。随着科学技术的迅速发展，我们可以把碎纸的一些特征用平面扫描仪、数码相机、摄像机等设

7、备记录下来。而如何把一张张图像的特征转化成数字特征，并根据这些数据特征去建立相关的数学模型或计算机算法，从而借用计算机来帮助我们拼接图像，提高效率，就是问题的关键点了。传统的图像碎片自动拼接算法有蚁群优化算法、遗传算法等，我们通过分析数据文件发现这些算法并不适合于本题。原因是以上的算法是基于碎片有不规则边缘的基础上的，而本题中碎纸图像的切痕是规则的，且无文字识别能力。第1步，把所给的图像的数字信息(即像素)用matlab软件读取出来，得到了图像的数据矩阵，通过分析数据后发现，相邻两张图片的边缘像素具有较大的相似度。第2步，利用数字图像处理的方法，结合题目所给数据文件的说明，得到了图与图之间的灰

8、度相关关系的相似度函数，即取出数据矩阵的边缘列与其他图像的数据矩阵边缘列进行最短距离运算。第3步，建立以图与图之间的相似程度为基准的数学模型。这个模型的评价标准很简单，就是相似度函数的值。相似度函数的值越小，就认为两张图像越靠近，即匹配的概率就越大。第4步，通过求相似度函数值，对所得到的函数值进行排序，从中寻找到相似度函数值最小的图像，就得到了最佳匹配图。第5步，每匹配若干张图，显示一次结果，若发现有文字不连续或意思不通的，则进行人工干预。 对于问题（1），由于题中给出的碎纸片是由碎纸机纵向切割而得到的，于是碎片边缘的尖点特征尖角特征，面积特征等几何特点几乎一样，导致无法运用以碎纸片有不规则边

9、缘为基础的传统计算机算法进行拼接。因为边缘相似的碎纸片的拼接，理想的计算机拼接不仅要考虑边缘的匹配还要满足字迹断线或碎片内的内容的相符。然而这种理想方法很难实现，于是利用数字图像处理的方法，建立了以图与图之间的相似程度为基准的数学模型。首先，用matlab软件的imread函数进行图像的读取，得到数据矩阵为。其次，根据问题的假设（1）可以知道当灰度图中后n列灰度值恒等于255时即可认为这一碎纸为最右端的碎纸片，并记为。然后，以数据矩阵为基础，引入相似度函数。最后，用matlab工具箱中的sort函数把所得到的相似度函数值进行排序，所得到的相似度函数值最小的图像即为与最右端的碎纸片匹配的图像。如

10、此重复18次，就可以得到附件1的所有图像的排列序号。附件2的图像拼接算法与附件1 的图像拼接算法一致。由此问题（1）的算法模型确定完毕。 对于问题（2），附件3和附件4的碎片图像是碎纸机采用横切与纵切所得到的，与问题（1）相似的碎纸边缘的特点是一致的，无法采用几何特性拼接。于是可以利用问题（1）的模型，设计出以图像灰度系数相关的算法进行图像的拼接。同样采用matlab软件的imread函数进行图像的读取，转化为灰度图即可得到一个的数据矩阵。由于所被碎纸机纵横切割的是标准纸张，其具有页边距的特性。可以首先在209张碎纸片中寻找到应该位于页面最右端或者最左端的11张碎纸条，因为有页边距的原因，此时

11、它们的灰度图像第后n列灰度值全为255，并且具有相同的边距（即从后n+1行开始均不全为255）。采用同样的方法我们也可以得到位于页面上方的19张碎片。 根据所得到的图像我们可以人工寻找到位于最右端的那一列的第一张碎纸片序号。根据问题（1）的模型的图像拼接的灰度相关方法，我们将位于最右端的那一列的第一张碎纸片与剩余的10张碎片独立进行纵向拼接。得到被碎纸机纵向切割的最右端的一整列碎纸条。但是，在运用灰度度相关关系进行图片的匹配时，将会出现多个与第一张图像匹配的图像。于是我们采用人工干预，在符合条件的若干张（不超过5张）图像中找到能把文字信息完整拼接上的图像。可知最右端的一张碎纸片，以此碎纸条为基

12、准链，从右往左寻找可匹配横向（即第一行）方向图像，除了要根据灰度值相关关系，还要人工干预判断行距是否一致。在此过程中除了问题（1）的模型基础上运用相关算法外，人工干预也起到关键性作用。附件4的图像拼接算法与附件3的图像拼接算法一致。由此问题（2）算法模型确定完毕。问题（3）中所给出的碎纸是由碎纸机切割一页英文印刷文字双面打印文件而得到的。一张碎片有正反两面，并且所给数据中并不能把正反面分开，导致在寻找匹配图像时难度增加。但是，依然可以运用印刷纸张的页边距特性寻找位于最右边碎纸。通过计算机搜索可以找到22张位于页面最右端和最左端的碎片图像。再次运用标准纸的性质，计算机搜索位于页面上方的38张碎纸

13、图像。我们可以通过以找出的碎片图像中找到两幅即位于页面上方与最右边的图像，于是可以得到位于正面右上端的正面的一幅图像，用矩阵表示为，另一幅图像为反面的左上端的图像用矩阵表示为。可以进行人工干预，找到了位于原文件最左端的第1个碎纸图像和最右端的第一个图像。运用模型一的算法寻找与之匹配的图像：第一步，将所得碎片图像分别从上往下纵向寻找匹配图像，因为正反面图像往下搜寻匹配的图像是同一张的，通过两个相似度函数寻找到匹配图像。通过拼接即可得到碎纸机纵向切割的最右的一张碎纸。第二步，以最右端的拼接而成的碎纸图像为基准链从右往左寻找匹配图像。按照解决问题（2）的算法寻找匹配图。在此过程中，会有人工干预从满足

14、灰度值要求的情况下找到行距与基准链相同的匹配图像。通过matlab软件编程可得原文件。三、 模型的假设与符号的说明3.1 模型的假设（1）假设完整的图像是一张有边界的标准纸张纸，即有明显的边界特性。（2）假设完整的图像在切割时和切割后边界整齐，没有不规则的边缘。（3）假设完整的图像的像素点录入过程没有噪声干扰。（4）假设完整的图像中的文字是规则的，即大小一致。3.2 符号的说明四、 模型的建立与求解4.1 问题（1）的模型建立与求解运用matlab软件imread函数进行图像的读取，先读取任意的一个的数据矩阵为。通过将碎纸图像转化为矩阵，运用计算机搜索法寻找最左端碎纸条，因为碎纸是由a4纸切割

15、而成的。根据A4纸页边距的性质，于是可得到当灰度图中后n列灰度值横等于255时，矩阵为：由以上通过图像处理的方法和计算机搜索方法，求得最左端的碎纸。接下来将要通过寻找与此碎纸片相似的右边碎纸片，引入了相似度函数。假设左边的图像灰度图矩阵为右边图像的灰度图矩阵为于是相似度函数为：将最左端的一张碎片与另外碎纸进行匹配，求出其相似度函数值（i，j=1,2,3,4,18）。继而对得出的函数值进行排序，运用matlab工具箱中的sort函数寻找最小函数值。即此时可以得到最佳的匹配图.由问题（1）中附件1所给出的图像可以人工观察得到最右端对应碎纸图像编号为（06），接下来是寻找与第19幅图像匹配的下一张图

16、像。通过求第19幅图像与另外18张图像的相似度函数值，并找出相似度函数值最小的匹配图像，则该图像能与第19幅图匹配。接着往下寻找与第18幅图匹配的第17幅图。第一张图像的程序算法框图为：开始输入图像输入m=1，即取D=图像最后一列-255DIS=D2m=m+1m19?把DIS排序，取最小值输出DIS结束图1 确定问题（1）第一张图像的程序算法框图其他18张图像的程序算法框图为：开始输入第一张图像输入m=1，即取B=DIS图像第一行C=m对应图像最后一行D=B-CDIS=D2m=m+1m19?把DIS排序，取最小值输出DIS结束图2 问题（1）其余18张图像的程序算法框图问题（1）的算法框图的具

17、体代码见附录1。4.2 问题（2）模型的建立与求解运用matlab在图像处理的运用，将209张碎纸图像转化为灰度图，得到第i张图像的数据矩阵为，在第i张图像上点的灰度值为。所得到的矩阵为：通过计算机编程搜索所有图像的数据矩阵，找到前n行的灰度值全为255的图像，所得的数据矩阵为：（）通过人工干预剔除边距不一致的碎纸片，从而筛选出具有相同边距的图像，而所选出的19张等边距图像即为位于最上方的图像碎片。接着运用等边距的特性可同样寻找到位于最右边的碎纸片图像（共11张），其数据矩阵满足：通过以上的步骤筛选，可以得到顶部与最右边的碎纸片图像。通过人工干预结合以上所得的最右边的碎纸片图像（共11张）可得

18、到位于原图像右上端的第1幅图像。接着下一个步骤是寻找与之匹配的纵向或横向图像。由于图像的边缘几何关系一致，在寻找最佳匹配图时，同样可采用问题（1）模型的算法。首先寻找纵向的匹配图，假设最右端的第一个图像的矩阵为。相似度函数为： 每个图像与右上端的第1幅图像都有一个相似度函数关系，再进行排序求出最小函数值从而寻找最佳匹图。然而运用计算机寻找最佳匹配图时，得到相似度函数值的最小值有多个，即有多张可以匹配的图像，此时应采取人工干预读取拼接图像是否符合文字逻辑与行距，从而选取相符合的图像。最上方与最右端的碎纸片图像（分别为19张与11张）的程序算法框图如下具体的程序代码见附录2问题（2）的程序开始输入

19、图像输入m=1，即取D=图像最后20列-255DIS=D2m=m+1m209?把DIS排序，取最小值输出DIS结束 开始输入图像输入m=1，即取D=图像第20行-255DIS=D2m=m+1m209?把DIS排序，取最小值输出DIS结束图3 最上方与最右端的碎纸片图像的程序算法框图于是我们通过寻找得到最右边的图像为基准链，寻找在横向方向上的匹配图像。寻找的方法与以上相同。再次运用相似度函数：继而将所得到的函数值进行排序，选取图像之间相似度函数值最低的图像进行拼接。运用相似度函数为：问题（2）中，附件2的最右端的11张图像的序列程序框图算法如下开始输入DIS图像输入m=1，即取Z=读入DIS图像

20、中的第m图像B=读出的第一张图最后一行Y=取Z的第一行m11?人工确认最后一列第一张图X=Y-BW=SUM(X2)m=m+1输出WW=MIN(W)所在图像QB=读出W最后一行Y=取W的第1行X=Y-BW=SUM(X2)m11?把所有图像组合在一起输出图像结束图4 问题（2）中，附件2的最右端的11张图像的序列程序框图算法框图的具体程序见附录问题（2）的代码，如果把上面的程序框图中数据矩阵B的命令改为取W的第1列，数据矩阵Y的命令改为取W的最后1列，即可得到最上方的11张图像的序列。问题（3）模型的建立与求解 由于所给的的图像数据是来自两面打印文件的碎纸片，所给图像中并不能判断正反面。这样在寻找

21、匹配图像时，很容易就找到满足图像灰度值相关关系一致的多张图像。于是应先找到一条基准链，以双重相似函数值为判断标准，适时进行人工的干预。首先，寻找位于页面右端的碎纸图像，图像矩阵克表示为为。在寻找位于页面上方的图像用矩阵表示。根据问题（2）模型中的算法可得到位于原文件左端和右端的第一张碎纸片。接着，我们引进双重相似度函数，寻找与右端第一个碎纸片图像匹配的纵向方向上的图像。双重相似度函数为：最后，用已得到的位于源文件右端的碎纸图像进行拼接，可得碎纸机纵向切割位于最右端的碎纸片。以此为基准链寻找右边与之匹配的图像，此过程与问题二解法相同，采用matlab与适时的干预。于是可运用附件数据的出结果。五、

22、 模型的结果分析与评价5.1 模型的结果分析5.1.1 问题（1）的结果分析根据问题（1）的程序框图，利用计算机编程可以得到附件1的匹配的图像顺序为表 1 所求中文碎片图像序号的排列方式47381619121621014119131518175根据问题（1）的附件1的数据可得完整图片见附件1中图1根据问题（1）的程序框图，利用计算机编程可以得到附件2的匹配的图像顺序为表 2 所求英文碎片图像序号的排列方式91513164113172561014191281817所拼接的完整图像见附件1中的图2由拼接所得图像可知，符合相关文字特征与内容逻辑关系。于是可以断定通过模型一的算法进行图像的拼接是可行的

23、。采用图像处理方法，以拼接无边缘特性的图像。5.1.2 问题（2）的结果分析通过计算机编程即可恢复原图像所得的结果用表格的形式表示如下：问题（2）附件3的图像排列序号表 3 所求中文碎片图像序号的排列方式049054065143186002057192178118061019078067069099162096131079168100076062142030041023147191038148046161024035081189122103071156083132200017080033202198014128003159082199135012073160094034084183090047

24、121042124144125013182109197016184110187066029064111201005092180048037075007208138158126068175045174000089146102154114040151207155140 19009501102212902809118814106311616307200617702005203605017912008619502600108701813019308816702500800910507401513317020508515216502706020316913403903105110711517607711

25、2149097136164127058043106150021173157181204139145055044206010104098172171059137053056093153070166032196185108117004101113194119123根据以上序号所拼接的完整图像见附件1中的图3同理可得，问题（2）附件4的图像排列顺序根据以上序号所拼接的完整图像见附件1中的图4表格 4 所求中文碎片图像序号的排列方式19107501115419018400210418006420114817019619809411316407810308605110702904015818609802

26、411701919409314108812112610515511415913900112906313815305303812302004110811613607303620713501520802100704906111903314216806207008406001406817413719500804913218109506916716316618811114401704206620501015707414508313408107712820013105212514019308710600414903220406503906714709108010102610000601702814615

27、0005059058092030037046127176182151022057202071165082120175085050160187097203031076043199045173079161179143169054192133118189162197112172156096023099122090185109206003130034013110025027178055018056035016009183152044089048072012177124000102115问题（2）所给的碎纸图像是采用横切与纵切所得，盲目进行搜索匹配图像不但计算量大而且很难用计算机编程实现。于是通过寻找图

28、像特征，找到只用纵切的最左边的碎片为基准链。再次寻找与基准链图像匹配的图像，运用模型一的算法，进行图像的拼接。从拼接的结果可知，符合图像内容与文字特征。于是此算法可行。5.1.3 问题（3）的结果与分析 问题（3）在以双面碎纸片为条件，进行图像的拼接。采用标准纸具有页面边距的特征，即图像的特征。构架新的相似度函数为匹配标准，结合问题（2）的算法得到理想的结果。5.2 模型的评价模型优点：（1）问题（1）的模型从问题的本质出发，摒弃以边缘几何方式的图像匹配方法，而是运用图像处理方式寻找图与图之间在计算机读取的图像的灰度值相关关系，建立了以图像相似函数为标准的图像匹配方式。简化了问题的处理过程，使

29、问题解决过程更加简单清晰。（2）模型充分利用图像的特征信息，能够有的放矢，不是在盲目的搜索，从而简化图像拼接的过程。模型所用算法比传统的灰度相关匹配算法速度要快,鲁棒性也很好。（3）模型解法相对简单，容易采用计算机编程实现。为了防止求解结果具有偶然性错误，多个过程加入了人工干预，误匹配发生的概率小。模型缺点:为了简化问题的需要，忽略了一些因素会导致图像无法准确找到最佳匹配图的因素，如噪声的干扰。计算的代价高，计算量大；该算法需要采用计算机搜索。以及求解相似函数值时计算量庞大。模型改进方向：本文在传统灰度相关匹配算法的基础上,建立了一以图像临界灰度值差异程度为判断标准的图像匹配数学模型。然而本文

30、所提出的模型具有一定的局限性。用碎纸机进行碎纸时得到的的碎纸片边缘特性是一致的，无法运用几何关系对其进行图像的匹配。这样大量的碎纸片运用灰度相关关系寻找相互匹配图像时，可能会出现两张相似度函数值是一样的。此时就比须采取人工干预，这样会浪费很多人力与物力。为了减少人工干预，可以采用一些新技术对碎片进行处理从而获取匹配关系。或者必须增加更多的约束条件，如文字的文字特征，色彩等特点，运用现代优化算法如模拟退火算法，遗传算法等从而寻找最佳匹配图形。六、 参考文献1赵小川.繆远诚.matlab数字图像处理实践，北京：机械工业出版社，2013.2张倩.占君.详解matlab图像函数及其应用，北京：电子工业

31、出版社2011.3赵静.但琦.数学建模与数学实验,北京：高等教育出版社，2008.4田伟丰.陈贝.基于灰度相关和区域特性的图像拼接算法，电子设计工程学报，2011,2（3）；184-1865胡社教.涂桂林.江萍.基于灰度相关图像拼接的改进算法，合肥工业大学学报.2008，6（6）；863-86七、 附录附录1问题（1）的完整程序%把图像存放到matlab工作目录下clearclose allclcfor m=1:19 n=m-1; if n10 filename=strcat(00,int2str(n),.bmp); else filename=strcat(0,int2str(n),.bmp

32、); end im=imread(filename); im=double(im);endfor m=1:19 im=im; d=im(:,72)-255; dis(m)=sum(d.2);endy,ind=sort(dis);l=sort(ind(1)f1=il;b=f1(:,1);b=double(b);for m=1:19 im=im; d=im(:,72)-b; dis(m)=sum(d.2);endy1,ind1=sort(dis);l1=sort(ind1(1)f2=il1;d=f2(:,1);d=double(d);for m=1:19 im=im; e=im(:,72)-d;

33、dis(m)=sum(e.2);endy2,ind2=sort(dis);l2=sort(ind2(1)f3=il2;g=f3(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y3,ind3=sort(dis);l3=sort(ind3(1)f4=il3;g=f4(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y4,ind4=sort(dis);l4=sort(ind4(1)f5=il4;g=f5(:,1);g=d

34、ouble(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y5,ind5=sort(dis);l5=sort(ind5(1)f6=il5;g=f6(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y6,ind6=sort(dis);l6=sort(ind6(1)f7=il6;g=f7(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y7,ind7=sor

35、t(dis);l7=sort(ind7(1)f8=il7;g=f8(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y8,ind8=sort(dis);l8=sort(ind8(1)f9=il8;g=f9(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y9,ind9=sort(dis);l9=sort(ind9(1)f10=il9;g=f10(:,1);g=double(g);for m=1:19 im=im; e

36、=im(:,72)-g; dis(m)=sum(e.2);end;y10,ind10=sort(dis);l10=sort(ind10(1)f11=il10;g=f11(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y11,ind11=sort(dis);l11=sort(ind11(1)f12=il11;g=f12(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y12,ind12=sort(dis);l12=s

37、ort(ind12(1)f13=il12;g=f13(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y13,ind13=sort(dis);l13=sort(ind13(1)f14=il13;g=f14(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y14,ind14=sort(dis);l14=sort(ind14(1)f15=il14;g=f15(:,1);g=double(g);for m=1:19 im=

38、im; e=im(:,72)-g; dis(m)=sum(e.2);end;y15,ind15=sort(dis);l15=sort(ind15(1)f16=il15;g=f16(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y16,ind16=sort(dis);l16=sort(ind16(1)f17=il16;g=f17(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y17,ind17=sort(dis);

39、l17=sort(ind17(1)f18=il17;g=f18(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;y18,ind18=sort(dis);l18=sort(ind18(1)f19=il18;g=f19(:,1);g=double(g);for m=1:19 im=im; e=im(:,72)-g; dis(m)=sum(e.2);end;f=f19,f18,f17,f16,f15,f14,f13,f12,f11,f10,f9,f8,f7,f6,f5,f4,f3,f2,f1;figure;im

40、show(f,)附录2问题（2）程序 部分clearclose allclcfor m=1:209 n=m-1; if n10 filename=strcat(00,int2str(n),.bmp); else if n100 filename=strcat(0,int2str(n),.bmp); else filename=strcat(,int2str(n),.bmp); end end im=imread(filename); im=double(im);endfor m=1:209 im=im; d=im(:,72)-255; k=im(:,71)-255; v=im(:,70)-25

41、5; q=im(:,60)-255; d=double(d); k=double(k); v=double(v); q=double(q); w=d+k+v+q; dis(m)=sum(w.2);endy,ind=sort(dis);l=sort(ind(1:11)%z=imread(123.bmp)a1=imread(141.bmp);%c=imread(036.bmp);b=a1(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2);endy1,ind1=sort(di);l1=sort(ind1(1)

42、;c1=l(l1)a2=ic1;b=a2(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2);endy2,ind2=sort(di);l2=sort(ind2(1);c2=l(l2)a3=ic2;b=a3(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2);endy3,ind3=sort(di);l3=sort(ind3(2);c3=l(l3)a4=ic3;b=a4(180,:);b=double(b);for m=

43、1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2);endy4,ind4=sort(di);l4=sort(ind4(2);c4=l(l4)a5=ic4;b=a5(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2);endy5,ind5=sort(di)a=a1;a2;a3;a4;a5;figure;imshow(a,)/问题（2）的附件3最右端那一列的程序clearclose allclcfor m=1:209 n=m-1; if n10 filename=str

44、cat(00,int2str(n),.bmp); else if n100 filename=strcat(0,int2str(n),.bmp); else filename=strcat(,int2str(n),.bmp); end end im=imread(filename); im=double(im);endfor m=1:209 im=im; d=im(:,72)-255; k=im(:,71)-255; v=im(:,70)-255; q=im(:,60)-255; d=double(d); k=double(k); v=double(v); q=double(q); w=d+k

45、+v+q; dis(m)=sum(w.2);endy,ind=sort(dis);l=sort(ind(1:11)%z=imread(123.bmp)a1=imread(141.bmp);%c=imread(036.bmp);b=a1(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2);endy1,ind1=sort(di);l1=sort(ind1(1);c1=l(l1)a2=ic1;b=a2(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b;

46、 di(m)=sum(d.2);endy2,ind2=sort(di);l2=sort(ind2(1);c2=l(l2)a3=ic2;b=a3(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2);endy3,ind3=sort(di);l3=sort(ind3(2);c3=l(l3)a4=ic3;b=a4(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2);endy4,ind4=sort(di);l4=sort(in

47、d4(2);c4=l(l4)a5=ic4;b=a5(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2);endy5,ind5=sort(di)%l2=sort(ind2(1);%a3=i;%a=a1;a2;a3;a4;a5;%e=a;c;%figure;imshow(a,) l5=sort(ind5(6);c5=l(l5)a6=ic5;b=a6(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2);endy6,ind6

48、=sort(di)l6=sort(ind6(1)c6=l(l6)a7=ic6;b=a7(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2); a=a1;a2;a3;a4;a5;a6;a7; endy7,ind7=sort(di)l7=sort(ind7(5)c7=l(l7)a8=ic7;b=a8(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2); endy8,ind8=sort(di)l8=sort(ind8(1)

49、c8=l(l8)a9=ic8;b=a9(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2); endy9,ind9=sort(di)l9=sort(ind9(1)c9=l(l9)a10=ic9;b=a10(180,:);b=double(b);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2); endy10,ind10=sort(di)l10=sort(ind10(3)c10=l(l10)a11=ic10;b=a11(180,:);b=double(b

50、);for m=1:11 k=l(m); g=ik; d=g(1,:)-b; di(m)=sum(d.2); end a=a1;a2;a3;a4;a5;a6;a7;a8;a9;a10;a11;figure;imshow(a,);y4,ind4=sort(di);第一个编程clearclose allclcfor m=1:209 n=m-1; if n10 filename=strcat(00,int2str(n),.bmp); else if n100 filename=strcat(0,int2str(n),.bmp); else filename=strcat(,int2str(n),.b

51、mp); end end im=imread(filename); im=double(im);enda1=i050;a2=i062;a3=i169;a4=i039;a5=i072;a6=i015;a7=i095;a8=i126;a9=i030;a10=i008;a11=i090;a=a1;a2;a3;a4;a5;a6;a7;a8;a9;a10;a11;a1=i055;a2=i020;a3=i101;a4=i149;a5=i157;a6=i129;a7=i035;a8=i014;a9=i065;a10=i209;a11=i147;b=a1;a2;a3;a4;a5;a6;a7;a8;a9;a10

52、;a11;a1=i066;a2=i079;a3=i077;a4=i047;a5=i084;a6=i004;a7=i085;a8=i183;a9=i112;a10=i139;a11=i103;c=a1;a2;a3;a4;a5;a6;a7;a8;a9;a10;a11;a1=i144;a2=i068;a3=i063;a4=i162;a5=i133;a6=i160;a7=i184;a8=i110;a9=i202;a10=i159;a11=i155;d=a1;a2;a3;a4;a5;a6;a7;a8;a9;a10;a11;a1=i187;a2=i070;a3=i143;a4=i025;a5=i201;a

53、6=i083;a7=i091;a8=i198;a9=i006;a10=i127;a11=i115;e=a1;a2;a3;a4;a5;a6;a7;a8;a9;a10;a11;a1=i003;a2=i100;a3=i031;a4=i036;a5=i018;a6=i200;a7=i048;a8=i017;a9=i093;a10=i069;a11=i041;f=a1;a2;a3;a4;a5;a6;a7;a8;a9;a10;a11;a1=i058;a2=i163;a3=i042;a4=i082;a5=i081;a6=i136;a7=i122;a8=i185;a9=i181;a10=i176;a11=i152;g=a1;a2;a3;a4;a5;a6;a7;a8;a9;a10;a11;a1=i193;a2=i097;a3=i024;a4=i190;a5=i034;a6=i013;a7=i043;a8=i111;a9=i049;a10=i046;a11=i208;h=a1;a2;a3;a4;a5;a6;a7;a8;a9;a10;a11;a1=i179;a2=i132;a3=i148;a4=i123;a5=i203;a6=i074;a7=i125;a8=i188;a9=i038;a10=i175;a11=i156