碎纸片的拼接复原数学模型

1、碎纸片的拼接复原数学模型摘要首先对图片进行灰度化处理，然后转化为0-1二值矩阵，利用矩阵行（列）偏差函数，建立了基于最小二乘法的碎纸片拼接数学模型，并利用模型对图片进行拼接复原。本文主要研究了规则碎纸片的拼接复原问题。首先利用二值法、环形像素点匹配等算法建立基于像素点数值匹配模型，然后利用MATLAB软件对碎纸片像素点进行数字化处理，得到各碎纸片的像素点数值矩阵，再利用MATLAB软件编程进行矩阵特征优化匹配得到复原图。针对问题一，当两个数字矩阵列向量的偏差函数最小时，对应两张图片可以左右拼接。经计算，得到附件1的拼接结果为：9,15,13,16,4,11,3,17,2,5,6,10,14,1

2、9,12,8,18,1,7附件2的拼接结果为:4,7,3,8,16,19,12,1,6,2,10,14,11,9,13,15,18,17,5针对问题二：本文采用聚类分析方法建立优化模型。将所有图片的像素转化为矩阵，其次采用人工干预的方式找出第一行第一张碎纸片，然后将所有图片灰度二制化，将灰色像素值转化为黑色像素值，然后由上至下找出黑白像素的分界线，根据分界线位置的差异值最小聚合出第一行所有的图片，同理聚合出第一列所有的图片，将第一行所有图片利用问题一的拼接算法求解出第一行图片拼接顺序。以第一行第一张碎纸片下侧像素点与其他碎纸片上侧像素点像素匹配差异值最小为目标函数，依次穷举第一列所有的碎纸片，

3、从而找出第二行第一张碎纸片。根据第一行排好顺序的碎纸片以及第二行第一张碎纸片，采用上下侧边缘像素匹配以及左右侧边缘匹配，依次找出并排列好第2行，其他行依理重复步骤（2）、（3）即可。完成匹配后，进行适当的人工干预即可确定碎纸片的复原图。（见附录3、4图11，图12）。最后，本文对所建模型的优点和缺点进行了客观的评价，认为本文研究的结果在图像拼接问题中有一定的参考价值。关键词：灰度处理，图像二值化，最小二乘法，聚类分析，碎纸片拼接一、问题重述碎纸片的拼接复原技术在司法鉴定、历史文献修复与研究、军事情报获取以及故障分析等领域都有着广泛的应用。近年来，随着德国“斯塔西”文件的恢复工程的公布，碎纸文件

4、复原技术的研究引起了人们的广泛关注。传统上，拼接复原工作需由人工完成，准确率较高，但效率很低。特别是当碎片数量巨大，人工拼接很难在短时间内完成任务。随着计算机技术的发展，人们试图开发碎纸片的自动拼接技术，以提高拼接复原效率。对于一页印刷文档，针对不同的破碎方法，讨论下列三个问题：（1）将给定的一页印刷文字文件纵切，建立碎纸片拼接复原模型和算法，并针对附件1、附件2给出的中、英文各一页文件的碎片数据进行拼接复原。（2）对于碎纸机既纵切又横切的情形，设计碎纸片拼接复原模型和算法，并针对附件3、附件4给出的中、英文各一页文件的碎片数据进行拼接复原。（3）对于双面打印文档，研究如何进行碎纸片的拼接复原

5、问题。附件5给出的是一页英文印刷文字双面打印文件的碎片数据。要求尝试设计相应的碎纸片拼接复原模型与算法，并就附件5的碎片数据给出拼接复原结果。二、模型的基本假设1、待拼接的碎纸片来自同一页印刷文字文件。2、待拼接复原的碎纸片是规整的矩形。3、模型中的碎纸片长度、宽度和面积都相等。4、附件中照片都是同标准拍摄。三、符号说明表1 符号说明符号符号说明灰度值红色绿色蓝色矩阵裁截距裁截文字长度行间距裁截空白距离字体高度四、问题分析将不规则的文档碎纸片进行拼接，一般是利用碎纸片的边缘曲线，尖点、尖角、面积等几何特征，搜索与之匹配的相邻碎纸片。但对于边缘形状相似的碎纸片，这种基于边界几何特征的拼接方法失效

6、，拼接时不但要考虑待拼接碎纸片边缘是否匹配，还要判断碎片内的字迹断线或碎片内的文字内容是否匹配。本问题给定的碎纸片有以下几个特点：1、每一张碎纸片都是规整的矩形；2、所有的碎纸片的长度、宽度都相等，形状是完全一样的；3、每一张碎纸片里都包含着文字（汉字、英文），不存在空白的碎纸片；4、不同的碎纸片之间没有重叠部分。由于碎纸片的形状相同，因而不能针对碎纸片的几何特征建立数学模型；碎纸片间无重叠，也不能利用图像融合技术进行图像配准。根据上述分析，我们考虑将图片进行数字化处理，根据每张碎纸片上的边缘文字特征进行匹配，也就是利用图片边缘文字的像素进行最优化匹配。五、模型的建立与求解5.1问题一的建模与

7、算法 由于碎纸片本身不具有体现其拼接特性的数字特征，我们需要将其数字化、矩阵化，将问题转化为矩阵之间的相关性。5.1.1图片的灰度处理为了对碎纸片进行数字化，我们将图像进行灰度处理，取出图像中每一个像素点的灰度值。通过MATLAB的imread函数读入灰度图， dither函数转化为二值图，进一步得到图像数据矩阵，每张图片对应一个的灰度矩阵。5.1.2图片的二值化处理将图片进行灰度处理以后，每个像素的灰度值介于0255之间。灰度值不能直接用于文字图片的拼接，还须进行二值化处理。将图片放入直角坐标系，规定：若点的像素灰度值大于或等于T，该点用数值表示，并将其设定为白色；若点的像素灰度值小于T，该

8、点用数值表示，并将其设定为黑色。由此得到像素点的二值化函数：其中，为预先设定的全局灰度阈值。于是，每张图片的灰度矩阵转化为下列的数字矩阵：,其中5.1.3最小二乘法1、图片左右拼接的数学模型 设分别表示左右放置的两张图片对应的数字矩阵，定义前一个矩阵的最后一列与后一个矩阵的第一列之间的偏差函数为：其中，分别表示矩阵第列和第列的元素。对于给定的矩阵，若存在矩阵，使得与之间的偏差函数达到最小，则称与可以匹配，此时与对应的图片可以左右拼接。2、图片上下拼接的数学模型 类似地，设分别表示上下放置的两张图片对应的数字矩阵，定义上面矩阵的最后一行与下面矩阵的第一行之间的偏差函数为：其中，分别表示矩阵第行和

9、第行的元素。 对于给定的矩阵，若存在矩阵，使得与之间的偏差函数达到最小，则称与可以匹配，此时与对应的图片可以上下拼接。我们称上述基于数字矩阵之间列（或行）距离的图片拼接模型为最小二乘法拼接复原模型。5.1.4算法与求解（一）算法思想第一步，对附件中的19幅图片分别进行灰度处理，然后取灰度阈值T=125，进行二值化，得到19个数字矩阵，即图片的数字化。第二步，对上述19个数字矩阵进行检测，若存在一个矩阵的最左侧一列元素全是1，根据破碎图片的特点，则该图片即为从左边起第一张碎纸片，记为。第三步，计算与其余18张图片对应矩阵的列偏差值。若存在，使得达到最小，则即位第二张图片。重复上述的步骤，依次得到

10、所有碎纸片的排列，即可拼接成完整图片。（二）附件1、2的拼接复原结果附件1和附件2的拼接顺序如下表：（附件1的算法程序见附录一，复原图片见附录二；附件2的算法程序见附录三，复原图片见附录四）表2 附件1拼接顺序91513164113172561014191281817表3 附件2拼接顺序4738161912162101411913151817552问题二模型的建立与求解5.2.1问题二分析该问题的目标是建立最合适的碎纸片拼接复原模型和算法来复原既纵切又横切的图片。问题二与问题一的区别在于碎纸片的切割由纵切变成了既有纵切又有横切的，所以问题一采用的人工干预找出第一张图片及简单地通过前后两张碎纸片

11、的偏差值进行排列会产生相当大的误差，因此我们拟采用聚类分析方法建立优化模型来拼接附加3和附加4中的碎纸片。（1）我们人工干预选出一组上侧留有空白的图片作为备选图片，以及选取一张左侧跟上侧都留有空白的图片作为拼接的第一张图片，其次将选取的第一张图片为基准采用相应的拼合算法不断的向左向下匹配下一张图片。对于碎纸片上的文字处理可采取以下方式： 图1：附录3碎纸片文字像素处理图首先用Matlab中的imread（）函数读取图片，将图片像素转化为相应的矩阵，其次将矩阵中处于0255之间的数值所在的行的全部数值转化为0（即将所有灰色像素点转换为黑色像素点，并将其所在的行全部转化为黑色像素点），然后由上至下

12、找出黑白像素的分界线，并依次将第一张图片的分界线与其他图片的分界线进行差异匹配，从而分类聚合出第一行所有的图片。（2）利用问题一中的拼接算法，以第一张图片右侧边缘像素点与第一行所有图片左侧边缘像素点像素匹配差异值最小为目标函数，建立优化模型，从而排列出第一行所有图片的顺序。（3）同理，将剩余图片像素矩阵中处于0255之间的数值所在的列的全部数值转化为0，然后由左至右找出黑白像素的分界线，并依次将第一张图片的分界线与其他图片的分界线进行差异匹配，从而分类聚合出第一列所有的图片。（4）改进问题一中的算法，以第一张图片下侧边缘像素点与第一列所有图片上侧边缘像素点像素匹配差异值最小为目标函数，建立优化

13、模型，从而确定出第二行第一张碎纸片，然后将第二行第一张碎纸片的右侧边缘像素点与待排图片左侧边缘像素点以及与上一行中对应图片下侧边缘像素点与待排图片上侧边缘像素点的像素匹配差异值最小为目标函数，建立优化模型，依次穷举出所有未排序的图片，找出第二行第一张碎纸片。其他行碎纸片依照此理类推即可排出顺序。（5）所有图片匹配完成后，对于拼接错误的图片进行相应的人工干预即可求得最后的图像复原图。5.2.2问题二模型准备（1）图片像素的量化和图片贴近度考虑到要根据图片之间的像素贴近度进行匹配，需要用Matlab软件将碎纸片像素量化，将像素分布转化为相应的数据。为此，我们采用了imread（）函数，将各图片转化

14、成由像素构成的矩阵，方便了matlab对图片匹配的处理操作。图片贴近度用以描述两张图片的相近程度，主要由两张图片的像素决定。两张图片像素匹配差异值越小，图片贴近度越高；反之，图片贴近度越低。（2）图片聚类分析方法数学建模中，研究对样品或指标分类的一种多元统计方法，是依据研究对象的个体的特征进行分类的方法。聚类分析的基本思想是认为我们研究的样本或指标之间存在着程度不同的相似性。于是根据一批样本的观察指标，具体找出想死程度较大的样本聚为一类，关系密切的聚到一个小的分类单位，关系疏远的聚到一个大的分类单位，直到所有样本聚合完毕。这样就形成了一个个类，大大减少了数据处理的难度。在图片拼接过程中，将碎纸

15、片图像量化得到的矩阵中数值处于0255（白色的像素值为255，黑色的像素值为0）之间的数所在的行或列中所有的数值全部转化为0，找出黑白像素的分界线，根据分界线的匹配度进行聚合分类。5.2.3问题二模型的建立图2：问题二模型求解流程图（1）碎纸片预处理 预处理的目的是将待处理物体表示为适合于计算机处理的形式。对于碎纸片来说，就是将其像素化、数字化，为此，我们利用Matlab软件中的imread（）函数对碎纸片进行量化，由于题目中所给的碎纸片图像分辨率都是，所以我们将每张碎纸片都转换成了的矩阵。（2）碎纸片聚类分析及匹配 在碎纸片聚类过程中，首先将所有图片得到的像素矩阵灰度二制化，将灰色像素值转化

16、为黑色像素值，由上至下找出黑白像素分界线，根据上侧跟左侧都留有白的特性人工干预找到第一张图片，然后以与第一张图片黑白像素分界线与其他图片黑白像素分界线差异值最小为目标函数建立优化模型，目标函数如公式（1）所示。.目标函数如下:我们根据分界线位置的差异值最小，对碎片预处理得到的像素矩阵，聚类分析得到第一行所有的图片，同样运用聚合分析的思想得到第一列所有的图片，根据问题一的匹配方法求出第一行图片的拼接顺序。拼接效果如图3所示。 图3：聚类出的第一行的拼接效果图然后以第一行第一张碎纸片下侧像素点与其他碎纸片上侧像素点像素匹配差异值最小为目标函数，目标函数如公式（2）所示。目标函数： 依次穷举第一列所

17、有的碎纸片，找出第二行第一张碎纸片。其次，根据第一行排好顺序的碎纸片以及第二行第一张碎纸片，采用上下侧边缘像素匹配以及左右侧边缘匹配，目标函数如公式（3）所示。目标函数：依次找出并排列好第2行，拼接效果图如图4所示：图4：第二行的拼接效果图重复上述两个步骤即可得到其余各行的排列顺序。（3）人工干预及绘图观察通过Matlab的imshow()对聚类分析及匹配得到的图片顺序的图像，调整不吻合的图片，逐步使得得到的图像准确。5.2.4问题二模型的求解利用聚类分析求解出的附件3跟附件4第一行图片排列顺序见下表：表4：附录3中采用聚合分析得出的第一行图片顺序49546514318625719217811

18、81909511221292891188141 表5：附录4中采用聚合分析得出的第一行图片顺序1917511154190184210418064106414932204653967147利用聚类分析求解出的附件3跟附件4第一列图片排列顺序见下表： 表6：附录3中采用聚合分析得出的第一列图片顺序49611683814941252977189表7：附录4中采用聚合分析得出的第一列图片顺序1912018619159202087013217181附件3图片最终拼接顺序见表8，拼接图见附录六。表8 附件3拼接顺序4954651431862571921781181909511221292891188141

19、611978676999162961317963116163726177205236168100766214230412314719150179120861952608718381484616124358118912210313019388167258910574141283159821991351273160203169134393151107115176943484183904712142124144771121499713616412758431251318210919716184110187661061502117315718120413914529641112015921804837

20、755544206101049817217159720813815812668175451741137535693153701663219671156831322001780332021981513317020585152165276089146102154114401512071551401851081174101113194119123附件4图片最终拼接顺序见表9，拼接图见附录八。表9 附件4拼接顺序191751115419018421041806410641493220465396714720114817019619894113164781039180101261006172814686

21、511072940158186982411715055958923037461271919493141881211261051551141761821512257202711658215913911296313815353381231201758550160187972033120411081161367336207135157643199451737916117914320821749611193314216862169541921331181891621971127084601468174137195847172156962399122901851091321819569167163166

22、1881111442063130341311025271781714266205101577414583134551856351691831524481771282001315212514019387894872121771240102115六、模型分析与评价6.1模型的优点（1）在进行拼接4张图片之前都进行了简单地人工干预来确定第一张图片，由此可以大大提高后来图片拼接模型的运行效率。（2）问题一建立的通过像素匹配差异值最小为目标函数，通过循环来穷举剩余图片，比较相邻图片贴近度来建立模型求解，简单明了，易于理解，并且可操作性较强。（3）问题二采取的聚类分析方法层次分明，逻辑清晰，并且能有效减少

23、运算量，提高精确度，因此提高的拼接碎纸片的效率。（4）总的来说，我们采用的模型理解较为容易，也便于操作，对碎纸片拼接的问题有较高的参考价值。6.2模型的缺点（1）我们在拼接碎纸片过程当中自始至终都运用了位图像素点的概念，要理解建立的模型，需要对计算机的图像处理知识有一定的了解。（2）我们在各个问题碎纸片拼接之前进行了人工干预选出第一张，如果选择的第一张不够准确，可能造成建立的整个模型完全错误，因此有一定的难度来选择第一张图片，这给使用此模型的人增加了对图片观察力和分析力的要求。（3）此模型中需要适当的人工干预来调整初步确定的图片模型，因此需要使用此模型的人了解模型的细节以便进行适时的人工干预。

24、6.3模型的改进（1）由于问题一要求排列的图片较少，因此可以依此选择各个图片作为第一张，自动的进行拼接，最后从每张图片作为第一张的得到图片序列中，根据总的贴近度最大原则，选出最优的排列。（2）我们问题二中仅采用从上到下、从左到右的方法进行搜索排序，如果再进行从从上到下，从右到左进行搜索排序，应该可以通过人工干预对比来更容易的实现排序。参考文献1姜启源等，数学模型（第四版），北京：高等教育出版社，2011 2罗智中，基于线段扫描的碎纸片边界检测算法研究J仪器仪表学报，2011,32（2）：289-294.2白宗文，基于HALCON与图像拼接的文物修复系统设计与实现J，电子设计工程，2013,21

25、（9）：24-26.3Huei-Yung等，Reconstruction of shredded document based on image feature matching，Expert Systems with Applications，第39卷，第3期，3324-3332， 2012。4黄添强，陈智文，苏立超等利用内容连续性的数字视频篡改检测J，南京大学学报（自然科学版），2011,47（5）：493-503.5罗智中，基于线段扫描的碎纸片边界检测算法研究J仪器仪表学报，2011,32（2）：289-294.6李利军，李云伟，基于图像灰度的拼接技术研究J，计算机与数字工程，2007,

26、35（9）：128-130.7贾海燕，朱良家，周宗潭等，一种碎纸自动拼接中的形状匹配方法J，计算机仿真，2006,23（11）：180-183.附录附录一：clear,clc;A=zeros(19,19); for j=1:19str=C:UsersAdministratorDesktop第六轮题目第六轮题目B附件1-副本; I=imread(str,num2str(j),.bmp); i2=dither(I); a=i2(:,72,1); str=C:UsersAdministratorDesktop第六轮题目第六轮题目B附件1-副本; for i=1:19 I=imread(str,num

27、2str(i),.bmp); i2=dither(I);mi=i2(:,1,1); ni=a-mi; A(j,i)=sqrt(dot(ni,ni); endenda=0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0; for i=1:19 k=A(i,1); a(i)=1; for j=1:19 if(kA(i,j) k=A(i,j); a(i)=j; end endendfirst=1;for i=1:19 for j=1:19 if(A(i,j)=0) first=j; end endend str=C:UsersAdministratorDesktop第六轮题目

28、第六轮题目B附件1-副本;z= imread(str,num2str(first),.bmp);for m=1:18 m=1;str=C:UsersAdministratorDesktop第六轮题目第六轮题目B附件1-副本;x = imread(str,num2str(first),.bmp); y=imread(str,num2str(a(first),.bmp); x=z; z = x,y; first=a(first);endimshow(z)附录二：附录三：clear,clc;A=zeros(19,19); for j=1:19str=C:UsersAdministratorDeskt

29、op第六轮题目第六轮题目B附件2; I=imread(str,num2str(j),.bmp); i2=dither(I); a=i2(:,72,1); str=C:UsersAdministratorDesktop第六轮题目第六轮题目B附件2; for i=1:19 I=imread(str,num2str(i),.bmp); i2=dither(I);mi=i2(:,1,1); ni=a-mi; A(j,i)=sqrt(dot(ni,ni); endenda=0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0; for i=1:19 k=A(i,1); a(i)=1

30、; for j=1:19 if(kA(i,j) k=A(i,j); a(i)=j; end endendfirst=1;for i=1:19 for j=1:19 if(A(i,j)=0) first=j; end endend str=C:UsersAdministratorDesktop第六轮题目第六轮题目B附件2;z= imread(str,num2str(first),.bmp);for m=1:18 m=1;str=C:UsersAdministratorDesktop第六轮题目第六轮题目B附件2;x = imread(str,num2str(first),.bmp); y=imre

31、ad(str,num2str(a(first),.bmp); x=z; z = x,y; first=a(first);endimshow(z)附录四：附录五：clear,clc;for i=1:209 if(i=10) imgname=strcat(D:附件附件300,num2str(i-1),.bmp); elseif(i=100) imgname=strcat(D:附件附件30,num2str(i-1),.bmp); else imgname=strcat(D:附件附件3,num2str(i-1),.bmp); end I(:,:,i)=imread(imgname); left(:,i

32、)=I(:,1,i); right(:,i)=I(:,72,i); up(i,:)=I(1,:,i); down(i,:)=I(180,:,i);endleft=int32(left);right=int32(right);up=int32(up);down=int32(down);mAns=zeros(1,19);mAns(1)=15;for i=2:19 xsd=zeros(1,209); for j=1:209 if j=mAns(i-1) continue end xsd(j)=max(sum(abs(right(:,mAns(i-1)-left(:,j); for tmp_i=1:2

33、09 if sum(up(tmp_i,:)-255)=0 xsd(tmp_i)=; end end for k=1:i-1 xsd(mAns(k)=; end end format long g xsd; ma,maid=min(xsd); mAns(i)=maid; %人工干预 if mAns(i-1)=200 mAns(i)=136; end if mAns(i-1)=204 mAns(i)=170; end if mAns(i-1)=170 mAns(i)=135;end if mAns(i-1)=161 mAns(i)=204; end if mAns(i-1)=40 mAns(i)=

34、32;end if mAns(i-1)=32 mAns(i)=52;end if mAns(i-1)=116 mAns(i)=177; endendmyans=zeros(11,19);myans(1,:)=mAns;for hang=2:11 xsd=zeros(1,209); for j=1:209 xsd(j)=max(sum(abs(down(myans(hang-1,1),:)-up(j,:); for tmp_i=1:209 if sum(left(:,tmp_i)-255)=0 xsd(tmp_i)=; end end for tmp_i=1:hang-1 for tmp_j=1

35、:19 xsd(myans(tmp_i,tmp_j)=; end end end format long g xsd; ma,maid=min(xsd); maid; myans(hang,1)=maid; if myans(hang-1,1)=95 myans(hang,1)=126;end if myans(hang-1,1)=8 myans(hang,1)=72; end if myans(hang-1,1)=72 myans(hang,1)=90; end if myans(hang-1,1)=62 myans(hang,1)=169; end for i=2:19 for j=1:2

36、09 if j=myans(hang,i-1) continue endxsd(j)=max(sum(abs(right(:,myans(hang,i-1)-left(:,j)+sum(abs(down(myans(hang-1,i-1),:)-up(j,:); for tmp_i=1:hang-1 for tmp_j=1:19 xsd(myans(tmp_i,tmp_j)=; end end for tmp_j=1:i-1 xsd(myans(hang,tmp_j)=; end end ma,maid=min(xsd); maid; myans(hang,i)=maid; if myans(

37、hang,i-1)=95 myans(hang,i)=35; end if myans(hang,i-1)=85 myans(hang,i)=184;end if myans(hang,i-1)=91 myans(hang,i)=48;end if myans(hang,i-1)=48 myans(hang,i)=122;end if myans(hang,i-1)=125 myans(hang,i)=145; end if myans(hang,i-1)=78 myans(hang,i)=113;end if myans(hang,i-1)=113 myans(hang,i)=150;end

38、 if myans(hang,i-1)=128 myans(hang,i)=59; end if myans(hang,i-1)=59 myans(hang,i)=44;end if myans(hang-1,i)=95 myans(hang,i)=126;end if myans(hang,i-1)=14 myans(hang,i)=183;end if myans(hang,i-1)=17 myans(hang,i)=185; end if myans(hang,i-1)=205 myans(hang,i)=140; end if myans(hang,i-1)=65 myans(hang

39、,i)=112;end if myans(hang,i-1)=112 myans(hang,i)=202;end if myans(hang,i-1)=93 myans(hang,i)=181;end if myans(hang,i-1)=76 myans(hang,i)=56;end if myans(hang,i-1)=207 myans(hang,i)=11;end if myans(hang,i-1)=11 myans(hang,i)=105; end if myans(hang,i-1)=173 myans(hang,i)=172; end if myans(hang,i-1)=20

40、9 myans(hang,i)=139;end if myans(hang,i-1)=159 myans(hang,i)=127;end if myans(hang,i-1)=127 myans(hang,i)=69;end if myans(hang,i-1)=57 myans(hang,i)=94;end if myans(hang,i-1)=94 myans(hang,i)=154; end if myans(hang,i-1)=154 myans(hang,i)=71; end if myans(hang,i-1)=71 myans(hang,i)=167; end if myans(

41、hang,i-1)=167 myans(hang,i)=33;end if myans(hang,i-1)=157 myans(hang,i)=84;end if myans(hang,i-1)=84 myans(hang,i)=133;end if myans(hang,i-1)=18 myans(hang,i)=81;end if myans(hang,i-1)=147 myans(hang,i)=103;end if myans(hang,i-1)=152 myans(hang,i)=208;end if myans(hang,i-1)=109 myans(hang,i)=118;end

42、 if myans(hang,i-1)=118 myans(hang,i)=5;end if myans(hang,i-1)=66 myans(hang,i)=144;end if myans(hang,i-1)=96 myans(hang,i)=12;end if myans(hang,i-1)=23 myans(hang,i)=130;end if myans(hang,i-1)=79 myans(hang,i)=68;end if myans(hang,i-1)=100 myans(hang,i)=163;end if myans(hang,i-1)=80 myans(hang,i)=6

43、4;end if myans(hang,i-1)=164 myans(hang,i)=73;end if myans(hang,i-1)=73 myans(hang,i)=7;end if myans(hang,i-1)=7 myans(hang,i)=178;end if myans(hang,i-1)=53 myans(hang,i)=37;end if myans(hang,i-1)=2 myans(hang,i)=88;end if myans(hang,i-1)=88 myans(hang,i)=19;end if myans(hang,i-1)=39 myans(hang,i)=1

44、49;end if myans(hang,i-1)=149 myans(hang,i)=47;end if myans(hang,i-1)=47 myans(hang,i)=162;end if myans(hang,i-1)=36 myans(hang,i)=82;end if myans(hang,i-1)=89 myans(hang,i)=168;end endendmyans1=myans(8:11,1:19);myans(1:7,1:19);myans=myans1; myans-1;% A=reshape(ans,1,209);% for i=1:209% B(i)=I(A(i)+

45、1);% % i=i+1;% end% % f=horzcat(I(A(1:209);% imshow(B)imshow(I(:,:,50) I(:,:,55) I(:,:,66) I(:,:,144) I(:,:,187) I(:,:,3) I(:,:,58) I(:,:,193) I(:,:,179) I(:,:,119) I(:,:,191) I(:,:,96) I(:,:,12) I(:,:,23) I(:,:,130) I(:,:,29) I(:,:,92) I(:,:,189) I(:,:,142);I(:,:,62) I(:,:,20) I(:,:,79) I(:,:,68) I

46、(:,:,70) I(:,:,100) I(:,:,163) I(:,:,97) I(:,:,132) I(:,:,80) I(:,:,64) I(:,:,117) I(:,:,164) I(:,:,73) I(:,:,7) I(:,:,178) I(:,:,21) I(:,:,53) I(:,:,37); I(:,:,169) I(:,:,101) I(:,:,77) I(:,:,63) I(:,:,143) I(:,:,31) I(:,:,42) I(:,:,24) I(:,:,148) I(:,:,192) I(:,:,51) I(:,:,180) I(:,:,121) I(:,:,87

47、) I(:,:,196) I(:,:,27) I(:,:,2) I(:,:,88) I(:,:,19);I(:,:,39) I(:,:,149) I(:,:,47) I(:,:,162) I(:,:,25) I(:,:,36) I(:,:,82) I(:,:,190) I(:,:,123) I(:,:,104) I(:,:,131) I(:,:,194) I(:,:,89) I(:,:,168) I(:,:,26) I(:,:,9) I(:,:,10) I(:,:,106) I(:,:,75);I(:,:,72) I(:,:,157) I(:,:,84) I(:,:,133) I(:,:,201) I(:,:,18) I(:,:,81) I(:,:,34) I(:,:,203) I(:,:,199) I(:,:,16) I(:,:,134) I(:,:,171) I(:,:,206) I(:,:,86) I(: