碎纸片的拼接复原分析模型

1、碎纸片拼接复原的设计与实现摘要 破碎文件的拼接在司法物证复原、历史文献修复以及军事情报获取等领域都有着重要的应用。由于人工拼接效率较低，我们利用MATLAB软件编写程序，实现碎纸片拼接技术的计算机化，实现批量拼接，以节省人力和时间。（一）为寻找最吻合拼接方案，利用MATLAB软件中的imread函数，实现了碎纸片与矩阵之间的形式转化以便于碎纸图片的拼接比对。数值0至255表示图中某一像素点由黑到白的变化程度。再根据纸张的边界留白较多，通过计算每个矩阵第一列中各向量的元素和，可将所得和数值最大的列向量所在矩阵对应的碎纸片确定为左边界。经计算知：008图为整体图片的左边界。根据使吻合参数最小的原则

2、，可计算出下一张图片。重复此步骤，以此类推，每次都挑选出剩余图片中与前一幅图片吻合参数最小的作为与之相连接的碎纸图片。最后可得到附件1的答案矩阵为：008014012015003010002016001004005009013018011007017000006附件2与附件1的模型求解方法相同，最后得到：附件2的答案矩阵为：003006002007015018011000005001009013010008012014017016004（二）由于每片独立的拼接无法达到最佳效果。故我们以缩小比对范围的形式来进行优化。附件3类比于问题一，用相同的方法找到所有碎纸片的左边界，共11张，分别04906

3、1168038071014094125029007089根据行特征筛选出每行的碎纸图片以保证全部图片均得到分类。此时选取数量最少一行作为入手点，利用MATLAB软件进行图片拼接。但由于碎纸图片的行特征值有误差，故图片大块拼接正确，但与实际情况有细微差别。通过简单的人工检测得到准确行的排列顺序，由于行信息充足，借用第一问最终达到正确拼接效果。结果矩阵如 表12 所示。第二种情况的关键点是找出准确的行特征，由于汉字与英文书写格式不同，汉字均为方块字，易于定位。而英文由于特殊的书写方式，上下位置不同，不易于找到行特征值。如：英文字母最为密集行中点位置。需将附件4图片对应的矩阵转化为0-1列向量（空白

4、行为0，反之则为1）。利用此方法筛选出位于同一行的碎纸片并进行纸片拼接。后续拼接方法同附件3。结果矩阵如 表13 所示。（三）结合双面信息处理边界，得到边界特点。根据元音字母中心位置得出行高，从而更准确筛选同行图片。双面信息同时校准，更容易得到拼接顺序，人工干预少。结果矩阵如 表5 所示。模型一简单易处理，适用于含大量信息的碎纸片拼接且准确度高；模型二针对文字内容的中英文差异分别利用吻合参数和行高作为标准来筛选图片；模型三深入生活实际，考虑日常生活中反正面印刷情况并结合英文印刷特点，实用性高，双面信息同时校准，人工干预少。关键字：碎纸拼接 MATLAB 吻合参数 灰度一、问题重述破碎文件的拼接

5、在司法物证复原、历史文献修复以及军事情报获取等领域都有着重要的应用。传统上，拼接复原工作需由人工完成，准确率较高，但效率很低。特别是当碎片数量巨大，人工拼接很难在短时间内完成任务。随着计算机技术的发展，人们试图开发碎纸片的自动拼接技术，以提高拼接复原效率。请讨论以下问题：1. 对于给定的来自同一页印刷文字文件的碎纸机破碎纸片（仅纵切），建立碎纸片拼接复原模型和算法，并针对附件1、附件2给出的中、英文各一页文件的碎片数据进行拼接复原。如果复原过程需要人工干预，请写出干预方式及干预的时间节点。复原结果以图片形式及表格形式表达（见【结果表达格式说明】）。2. 对于碎纸机既纵切又横切的情形，请设计碎纸

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

7、19个碎片。(4) 附件5为纵横切碎片数据，每页纸被切为1119个碎片，每个碎片有正反两面。该附件中每一碎片对应两个文件，共有21119个文件，例如，第一个碎片的两面分别对应文件000a、000b。【结果表达格式说明】复原图片放入附录中，表格表达格式如下：(1) 附件1、附件2的结果：将碎片序号按复原后顺序填入119的表格；(2) 附件3、附件4的结果：将碎片序号按复原后顺序填入1119的表格；(3) 附件5的结果：将碎片序号按复原后顺序填入两个1119的表格；(4) 不能确定复原位置的碎片，可不填入上述表格，单独列表。二、模型假设1、 每张碎纸片的形状都是完全相同的长方形且每张碎纸片形状规则

8、；2、 文字打印清晰，无缺墨断墨情况；3、 纸片边缘整齐，无重叠，无损耗；4、 扫描过程中每张碎纸片的位置都是完全平行的，不会出现倾斜的情况；5、 假设恰好能完全拼接，即碎片无缺失，也没有其他碎片混杂；6、 纸片无倒转；7、 碎片文字均为相同字号，字号大小适中；8、 文字印刷体行高、行间距相同；9、 页边距非0，但较小。三、字符说明1、表示第张图片所转化得到的矩阵；2、表示第个图所形成的矩阵的第行列对应灰度值；3、表示第个图所形成的矩阵的第列；4、表示第个图和第个图之间的吻合参数（其含义在模型分析中说明），其中两张图的吻合参数越低表示这两张图越吻合；5、表示图的0-1列向量四、模型分析本文针对

9、三种碎片拼接类型分别建立数学模型。模型一简单易处理，适用于单片含大量信息的碎纸片拼接且准确度高；模型二针对文字内容的中英文差异分别利用吻合参数和行高作为标准来筛选图片，且以缩小比对范围的形式应用于庞大数据量的整理进而进行碎片拼接复原；模型三深入生活实际，考虑日常生活中反正面印刷情况并结合英文印刷特点，实用性高，准确度大。问题一： 利用MATLAB程序中的imread函数将附件1的19张碎纸片分别转化成矩阵形式，每一张碎纸片都可转化为一个的矩阵。for循环的使用可以很大程度上节省时间、人力，再根据边界准则通过计算首列数字变化来确定其边界。最后根据吻合参数的大小确定最后的碎纸片拼接顺序。求解步骤如

10、下：确定最后整体图片的左边界将碎纸片转换为数字矩阵求出其他图片与左边界图片的吻合系数比较所有吻合系数，系数最小所对应的图片为第二张图片重复步骤3，求出其他矩阵与第二张图片的吻合系数，吻合系数最小的为第三张图片依照上述方法，求得最后所有图片的拼接顺序完成拼接效果图（图1：问题一求解思路图）问题二：同样利用第一问中求得图片左边界的方法，附件3、附件4中两幅图片的左边界。找出与左边界11张图片同行的其余碎纸图片。编程筛选出每行的图片，利用问题一中的方法进行计算机拼接，但由于拼接存在一定的误差，所以部分图片拼接不完整正确，在此情况之下进行人工干预，观察图片内容，根据内容手动移动碎纸图片进行拼接。求解步

11、骤如下：找到与每幅边界图片同行的剩余图片找出左边界中的图片（共11张）选取所有行中图片数量最少的一行进行图片拼接，若拼接不准确则引入人工干预确定第行的图片拼接顺序，并在其他行删除在第行中出现的图片在剩余行中重复步骤3和步骤4，直至所有行图片拼接结束将每行看做整体，进而完成完成拼接效果图将碎纸片转化为数字矩阵（图2：问题二求解思路图）问题三：对于有双面打印的碎纸片，根据边缘留白情况，可筛选出22个左边界，并根据字母中点所在行的行高分为11类。并根据行高对其他图片进行筛选分类。引入人工干预，完成效果拼接。求解步骤如下：找出左边界中的图片（共22张）将边界图片进行分类，共分为11类由所在行的两张边界

12、图片筛选出所在行的其他碎纸图片将每行碎纸图片利用程序进行计算机拼接，拼接不完全，引入大量人工干预完成拼接效果图将碎纸片转化为数字矩阵同时考虑正反两面。加大拼接的正确性（图3：问题三求解思路图）五、模型建立与求解5.1附件1中碎纸片拼接图片是否能拼接，重点要看边缘的吻合程度，于是我们首先给出定义：表示第个图所形成的矩阵的第行列对应灰度值，则第个图与第个图所分别形成矩阵的第行列的差的绝对值，即记为第个图与第个图在边缘点的吻合参数，其中两张图在该点的吻合参数越低表示这两张图该点越吻合。对于两图边缘列，则有，为边缘的吻合参数，其中两张图边缘的吻合参数越低表示这两张图越吻合；故，我们用作为全局的吻合参数

13、，其中的限定条件为具体求解过程如下：1）为将19张碎纸片拼凑成一幅完整图片，需对由每一张碎纸片转化得到的数字矩阵进行处理。其中000图转化为矩阵,001图转化为矩阵， ，018图转化为矩阵。（见数据表1）2）第个图所形成的矩阵的第1列就是第个图的左边界，第72列表示第个图的右边界，列中数字的变化表示图边缘文字的灰度变化。灰度是指黑白图像中点的颜色深度，范围一般从0到255。白色为255，黑色为0。将每个矩阵首列向量的每个元素相加求和，因为完整图片的左右边界均为空白，故求和数值越大代表空白越多，通过计算找出19个矩阵中，首列数字无变化的（即首列向量和最大的）矩阵为由008图所形成矩阵。则可确定0

14、08图为最后完整图的左边界。3）再利用for循环语句将矩阵的最后一个列向量与另外每个矩阵的第一个列向量元素对应相减，会得到由每个元素差值组成的新向量。再将得到的新向量中的每个元素的绝对值相加求和得到一个数值。不妨将此求和所得数值称为两向量的吻合参数，记为。4）吻合参数最小的矩阵所对应的碎纸片即为第二张。从而我们得到数学模型：第二个图为：达到该最小值时所对应的矩阵为第个矩阵，我们用表示为所选出的第个矩阵的标码，以便后续求解；5）依据第二张碎纸片所对应矩阵的最后一列与其他矩阵第一列对应元素相减，再绝对值求和进行吻合参数的比较，从而第三个图为：达到该最小值时所对应的矩阵为第个矩阵，我们用表示为所选出

15、的第个矩阵的标码，以便后续求解；以此类推:排序后的第个图为：，达到该最小值时所对应的矩阵为第个矩阵，我们用表示为所选出的第个矩阵的标码；6）以此类推，第18个图为：，达到该最小值时所对应的矩阵为第个矩阵，对应第个图为排序后第18个图，剩余的一个自然为最后一图即图片右边界。通过实际拼接实践，我们可以知道，局部拼接复原结果即全局最优拼接复原结果。最后我们得到附件1的拼接顺序：表1：008014012015003010002016001004005009013018011007017000006（复原图片见附录1.1）5.2附件2中碎纸片拼接附件2的处理办法与附件1相同，利用所有矩阵的第一列先找到所

16、有碎纸片中的左边界，再根据吻合参数确定碎纸片的拼接顺序。最后我们得到附件2的拼接顺序：表2：003006002007015018011000005001009013010008012014017016004（复原图片见附录1.2）由于附件1和附件2中碎纸片的行高较高，所以在吻合参数的计算过程中相对误差很小，利用在MATLAB软件中编写的程序即可求出附件1和附件2所有碎纸片的正确拼接顺序，故在此问中不涉及人工干预。5.3附件3中碎纸片拼接1）利用第一问在MATLAB中编写程序求出位于左边界的图片的思想，运行程序可以找到附件3中位于第一列中的图片共有11张，分别为049图，061图，168图，03

17、8图，071图，014图，094图，125图，029图，007图，089图。恰好符合矩阵的格式要求。即说明这11张图片为整体图片第一列的所有组成部分。2）用MATLAB编写循环语句，以已经确定的第一列的11张图片为标准分为11类，分别找出与第一列11个图空白行位置相近似的图片归为一类（此时可认为同一类中图片的文字行是对齐的，即同类图片理论上应处在相同水平位置）与图007（标号为8）空白位置相近的图片有22个，编号为 8115334654576971729094127138139154159167175176197209与图014（标号为15）空白位置相近的图片有57个，编号为 15459101

18、3262832334041526171747583869094102103106108109114115116118120124126129135136141147152153154155156160161166167170171177186195197200204206208与图029（标号为30）空白位置相近的图片有35个，编号为 30611172238454956606567769399105107110111112140146151158172173174181182185188198202205207与图038（标号为39）空白位置相近的图片有22个，编号为3991015252636

19、477275828990104106123131149162168190194与图049（标号为50）空白位置相近的图片有35个，编号为 5031217222329555866679296107110111119130140142144146151158174179182185187188189191193198205与图061（标号为62）空白位置相近的图片有22个，编号为 6271520213753646870727379809097100117132163164178与图071（标号为72）空白位置相近的图片编号有47个7291016171822262833345761677175818

20、48694106107110111133134140146151153154157158166167171174182185188197198199201203205206与图089（标号为90）空白位置相近的图片编号有69个9059101416171822262833344120857616771727581848694102103106107109110111114115118120124126133134140141146147151152153154155156157158166167171174182183185186188195197198199201203205206与图094（

21、标号为95）空白位置相近的图片编号有22个 95153543444859727885909198113122125128137145150165184与图125（标号为126）空白位置相近的图片编号有19个 126141517226790107110111146158174182183185188198205与图168（标号为169）空白位置相近的图片编号有37个 16921517192224273142516367778788101103107110111121140143146148151158174180182185188192196198205由于图片的切割都不完整导致空行较多，故通过

22、程序运行出的结果每一行的符合要求的图片都张。3）有一组恰为19块图，但通过第一问的模型对这19块进行横向拼接时，发现有三块明显不位于同一行，将这三块拿出，分别为,图014，图016，图089。再找相对较少（22块）的几行，发现其中重复的恰为首列中的某3个图，删去后，得到这几列的19个图分别进行拼接，其中首图为007图的拼接结果为：图4（人工干预）我们发现所选图片恰为该行图片，但拼接结果由于碎片边缘笔划变化导致相连部分误差较大，从而出现计算机拼接错误，人工干预将错误部分进行更改后，该行顺序为：表3：0072081381581260681750451740001370530560931530701

23、66032196同理，可得到其余行向找出为22块行的拼接顺序，分别为表4：038148046161024035081189122103130193088167025008009105074061019078067069099162096131079063116163072006177020052036再处理行中找出35或37个图片的行，首先通过程序将之前已使用过的图片排除，则首图为图029的行，剩余图片为表5：029005010037044048055059064076092098104111139151171172180201206用问题一中模型对该21块图进行拼接，发现图139和图150不

24、属于该行，将两图删除后，再次进行拼接，得到上一步中情况（有个别拼接不准确），（人工干预）将其调换位置，得到该行顺序为表6：029064111201005092180048037075055044206010104098172171059同理，得到首图为图049，图168的两行排列为表7：049054065143186002057192178118190095011022129028091188141168100076062142030041023147191050179120086195026001087018再用以上结果筛选47个图的行，剩余图片为表8：0710080090150170250

25、27032033060074080083085105132133152156165170198200202205用问题一中模型对该25块图进行拼接，发现图008、图009、图025、图032、图074、图105不属于该行，将两图删除后，再次进行拼接，得到上一步中情况（有个别拼接不准确），（人工干预）将其调换位置，得到该行顺序为:表9：071156083132200017080033202198015133170205085152165027060而对于一行中出现69块和67块的两行，重复较多，于是我们选用另一种方式筛选。当得到每块图所对应矩阵后，我们将有全空白行记为0，存在黑色部分行记为1，则

26、用mathematica得到该图片的一个0-1列向量，例如图089的对应0-1列向量为：其中，我们挑选出几列进行部分对比，左侧为其对应上部几行对比，右侧为其对应中部几行对比。与左边对应中间行位置00000111100000011110000001111000000111100000011110000001111000000111100000011111000001111100000111110000011111000001111100011111110011111110001111111000111111100011111110001111111000111111101011111110101

27、1111111101111111110111111111011111111101111111110111111111011111111101111111110111111111011111111101111111110111100111011110000101111000010111100001011110000观察所选出对应图片的0-1列向量对比，不难发现，其1-1对应，0-0对应较准确。若两个图片属于同一行，则其空白位置高度应对应相等，字体高度对应相等，针对这一特性，我们对两个列向量做内积后，同行的内积不变，建立模型得：与图089同行的图片应满足：尽量小，我们将结果由小到大排列，前10个即

28、为与图089同行的图片同样的与图014同行的图片应满足：尽量小，我们将结果由小到大排列，前10个即为与图014同行的图片从而对首图为图089和图014的两行又做了筛选，得到表10：089146102154114040151207155140185108117004101113194119123014128003159082199135012073160203169134039031051107115176则剩余为首图为图125一行中图片，用问题一中模型排序为表11：1250131821091970161841101870661061500211731571812041391453）完成行图片的

29、拼接，后续只相当于拼接矩阵，此时所有图片仅有横切情况，类似于第一问中仅纵切情况，利用MATLAB变成完成最后拼接效果。附件3的拼接顺序： 表12：04905406514318600205719217811819009501102212902809118814106101907806706909916209613107906311616307200617702005203616810007606214203004102314719105017912008619502600108701803814804616102403508118912210313019308816702500800910507

30、407115608313220001708003320219801513317020508515216502706001412800315908219913501207316020316913403903105110711517609403408418309004712104212414407711214909713616412705804312501318210919701618411018706610615002117315718120413914502906411120100509218004803707505504420601010409817217105900720813815812

31、6068175045174000137053056093153070166032196089146102154114040151207155140185108117004101113194119123该模型虽能准确拼出复原后图片，但在每一行挑选过程中不能分别准确的找出对应的19个图片，若对的每个列矩阵进行调整，相信可以更准确地快速地拼出图片。5.4附件4中碎纸片拼接 我们首先希望同样通过附件3的方法，即通过碎纸片左边界留白的情况，筛选出11张图片作为第一列。但是由于英语文章的特殊性，每个单词之间都有一定的间隔，加之碎纸片的高度较小，以至于不止11张纸条的左边界是空白，于是我们引入一种新的方法。

32、首先，当得到每块图所对应矩阵后，如果某行全为空白，则将该行记为0，否则记为1，则用mathematica将该图片转化为一个0-1列向量。从而筛选1元素较多的即英文字母密集的图片，由这些图片（人工干预）逐次筛选后，所留英文字母密集且行位置不同的，作为我们拼接的起点。再由其确定原因字母的中心行，进而确定出行高，由密集图片行高的不同位置，两边延伸，进行拼接，由于某些行为段落结尾处，其后空白较多，需人工干预，同样可拼接出11行英文字母，最后同附件3的处理方式，对11个行转化成矩阵进行拼接。附件4的拼接顺序表13：191075011154190184002104180064106004149032204

33、065039067147201148170196198094113164078103091080101026100006017028146086051107029040158186098024117150005059058092030037046127019194093141088121126105155114176182151022057202071165082159139001129063138153053038123120175085050160187097203031020041108116136073036207135015076043199045173079161179143208

34、0210070490611190331421680621690541921331181891621971120700840600140681741371950080471721560960230991220901851091321810950691671631661881111442060031300340131100250271781710421662050101570741450831340550180560350160091831520440810771282001310521251401930870890480720121771240001021155.5 附件5的拼接复原对于双面打印

35、的英文印刷图片 1）将碎纸图片转化为的数字矩阵2）确定附件5中所有碎纸图片中的边界图片，依据为纸张的边界留白更多这一准则，计算每个矩阵对应列向量和，形成行向量，可将所得行向量元素一次达到最大值的向量所在矩阵对应的碎纸片确定为左边界。3）又由于所有碎纸图片分为正反两面，每页纸被切为个碎片，故经过筛选之后共得到边界图片22张。其中有11张边界图片为反面的右边界。则这22张碎纸图片可分为11类。每一类的分类标准为两张碎纸图片中字母所在行灰度吻合度最高。依此标准得到的这11类边界分布。4）每一类的边界有正反两个边界，依据这两张边界图片中字母所占行灰度，在其他剩余碎纸图片中寻找与之灰度吻合参数最小的作为

36、同一行。由于英文字母本身占格位置的特殊性，误差较汉字来说要大。在进行同行碎纸图片拼接时准确率不高。此时，需引入大量的人工干预来进行拼接工作。5）经过大量的人工干预，最后我们可以得到：附件5的拼接顺序为表14：136b047a020a164b081b189b029a018b108a066a110a174b183b150a155a140a125a111b078b005a152a147a60b059a014a079a144a120b022a124b192a025b044a178a076b036a010b089a143b200b086b187b131b056b138a045a137b061b094b0

37、98a121a038a030a042b084b153a186b083a039b097a175a072b093a132b087a198b181b034a156a206b173b194b169b161a011b199b090a203b162b002a139b070b041a170b151b001b166b115b065b191a037b180a149b107a088b013a024a057a142a208a064b102b017b012a028b154b197a158a058a207a116b179b184b114a035a159a073b193b163a130a021b202a053b177b0

38、16b019b092b190b050a201a031a171b146a172a122a182b040a127a188a068b008b117b167a075b063b067a046a168a157a128a195a165b105a204b141a135b027a080b000b185a176a126b074b032a069a004a077a148b085b007b003b009b145a082b205a015b101a118b129b062a052a071b033b119a160b095a051b048a133a023b054b196b112a103a055b100b106b091a049b0

39、26b113a134a104a006a123a109a096b043a099a复原图片见附录图3.1和图3.2六、模型评价与改进优点：1问题一模型，简单、易操作。适用于单张碎片所含文字信息量大的图片拼接复原；2问题二模型，对于碎片上信息量较少，但数量庞大的情况，以减少比对范围来达到理想拼接复原效果，且效率较高；3问题三中模型，考虑了现实生活中的所有考虑因素，并运用MATLAB软件编程得出了结果，应此结果准确性高，误差小，切合生活实际；4实现碎纸片拼接技术的计算机化，实现批量拼接，以节省人力和时间。5假设合理，模型具有坚实可靠的数学基础。缺点：1.假设条件很多，这就直接导致与实际情况仍有一定的差

40、距，同时许多实际因素尚未纳入考虑的范围之内，没有考虑周全，因此导致一系列计算数据误差较大。2.因为缺乏某些现实统计数据的支持，导致人为主观因素在模型的建立计算过程中的影响举足轻重，因此在选择处理方法时也受到条件上的众多限制。故此模型并不完全贴合实际，不利于实际的应用。3.由于时间关系，一些改进思路来不及实现。模型改进：1.对于分类不够精确这一问题，我们应用mathematica软件对英文字母进行识别（程序见附录末），已识别出这些在英文书写中出现频率较高且仅占用中间行的字母，以此易于定位，可以提高分类的精确度，从而解决了这一问题。七、参考文献1.蔡旭辉刘卫国蔡立燕，MATLAB基础与应用教程，人

41、民邮电出版社，2009.8;2.尚明姝，基于特征点约束关系的图像拼接算法J，电视技术，第36卷第13期:2012;3.姜丽凤巩建学,基于特征点的图像拼接算法研究与实现,电脑开发与应用,第24卷 第6期： 2011年; 4.姜启源，数学模型（第四版），北京：高等教育出版社，2003.8;5.罗智中，基于文字特征的文档碎纸片半自动拼接，计算机工程与应用，2012;6.方静，基于特征点的图像拼接算法优缺点比较，电脑知识与技术，第9卷，第07期：2013.3。附录一：（针对问题一）附件1图片生成程序：clcclearp=cell(1,18);for i=1:10 imageName=strcat(I:

42、题目附件100,num2str(i-1),.BMP); pi = imread(imageName);endfor i=11:19 imageName=strcat(I:题目附件10,num2str(i-1),.BMP); pi = imread(imageName);endsump=;Answer=;for i=1:19 Answer(i)=i;endfor i=1:19 sump(i)=(sum(pi(:,1)+sum(pi(:,2);endA=find(sump=max(sump);C=Answer(A);Answer(A)=1;Answer(1)=C;B=pA;pA=p1;p1=B;f

43、or i=2:19 sump=; for j=i:19 a=pi-1(:,72); b=pj(:,1); s=abs(a-b); sump(j)=sum(s); for k=1:(i-1) sump(k)=9999999; end endA=find(sump=min(sump);C=Answer(A);Answer(A)=Answer(i);Answer(i)=C;B=pA;pA=pi;pi=B;endAnswer-1q=p1;for i=2:19q=q,pi;end imshow(q)附件1复原结果图片：图1.1附件2图片生成程序：clcclearp=cell(1,18);for i=1:10 imageName=strcat(I:题目附件200,num2str(i-1),.BMP); pi = imread(imageName);endfor i=11:19 imageName=strcat(I:题目附件20,num2str(i-1),.BMP); pi = imread(imageName);endsump=;Answer=;for i=1:19 Answer(i)=i;endfor i=1:19 sump(i)=(sum(pi(:,1)+sum(pi(:,2);endA=find(sump=max(sump);C=Answer