2013年全国大学生数学建模竞赛B题碎纸片复原（含word论文和源代码资源）

已于 2024-08-03 14:51:26 修改
阅读量961 收藏 8
点赞数 5
分类专栏： 数学建模 文章标签： 数学建模
于 2024-07-21 15:50:50 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_35500719/article/details/140588152
版权

文章目录

一、部分题目

2013高教社杯全国大学生数学建模竞赛题目

B题 碎纸片的拼接复原

破碎文件的拼接在司法物证复原、历史文献修复以及军事情报获取等领域都有着重要的应用。传统上，拼接复原工作需由人工完成，准确率较高，但效率很低。特别是当碎片数量巨大，人工拼接很难在短时间内完成任务。随着计算机技术的发展，人们试图开发碎纸片的自动拼接技术，以提高拼接复原效率。请讨论以下问题：

问题1 对于给定的来自同一页印刷文字文件的碎纸机破碎纸片（仅纵切），建立碎纸片拼接复原模型和算法，并针对附件1、附件2给出的中、英文各一页文件的碎片数据进行拼接复原。如果复原过程需要人工干预，请写出干预方式及干预的时间节点。复原结果以图片形式及表格形式表达（见【结果表达格式说明】）。

问题2 对于碎纸机既纵切又横切的情形，请设计碎纸片拼接复原模型和算法，并针对附件3、附件4给出的中、英文各一页文件的碎片数据进行拼接复原。如果复原过程需要人工干预，请写出干预方式及干预的时间节点。复原结果表达要求同上。

问题3 上述所给碎片数据均为单面打印文件，从现实情形出发，还可能有双面打印文件的碎纸片拼接复原问题需要解决。附件5给出的是一页英文印刷文字双面打印文件的碎片数据。请尝试设计相应的碎纸片拼接复原模型与算法，并就附件5的碎片数据给出拼接复原结果，结果表达要求同上。

二、部分论文

请添加图片描述

请添加图片描述
请添加图片描述
请添加图片描述
请添加图片描述
请添加图片描述
请添加图片描述
请添加图片描述
请添加图片描述
请添加图片描述
请添加图片描述
请添加图片描述
请添加图片描述

三、部分源代码

完整代码 https://github.com/yan-fanyu/CUMCM-Paper-And-SourceCode
%分类算法，一类为一个裁剪行
num=209;
for i=1:num  %读取所有碎纸片
    if i<11
        [img0(:,i*72-71:i*72),cmap0(:,i*3-2:i*3)]=imread(strcat('00',num2str(i-1)),'bmp');
    elseif i<101
        [img0(:,i*72-71:i*72),cmap0(:,i*3-2:i*3)]=imread(strcat('0',num2str(i-1)),'bmp');
    else
        [img0(:,i*72-71:i*72),cmap0(:,i*3-2:i*3)]=imread(num2str(i-1),'bmp');
    end
end
biaozhun=zeros(1,num);
for i=1:num
    temp=img0(:,72*i-71:72*i);
    temp2=zeros(180,1);
    for j=1:180
        if(sum(temp(j,:))==72*255)
            temp2(j)=255;
        end
    end
    for j=1:180
        if(temp2(j+1)==255 && temp2(j)~=255) %寻找第一个变为255的行数(黑变白)
            b=j;
            break;
        end
    end
    biaozhun(i)=b;
end

biaozhunpaixu=sort(biaozhun);
fenlei=zeros(11,19);
for i=1:11 %分类，19个分一类
    fenlei(i,:)=find(biaozhun<=biaozhunpaixu(i*19) & biaozhun>=biaozhunpaixu(i*19-18));
end


%寻找第一片
diyipian=zeros(1,11);
k=1;
for i=1:209
    if(sum(sum(img0(:,i*72-71:i*72-62)))==255*180*10)
        diyipian(k)=i;
        k=k+1;
    end
end
%把实际第一片换到第一位 
for i=1:11
    [m,n]=find(fenlei==diyipian(i));
    fenlei(m,n)=fenlei(m,1);
    fenlei(m,1)=diyipian(i);
end


gezhangtouying;
num=418;
%寻找第一个完整行的下边界作为分类依据
%计算第一个完整行的高度，分为三类:小字母[23,25] 一头或一尾或大写[34,39] [47,52]
biaozhun=zeros(num,1);
for i=1:num
    for j=1:179
        if(touying(j,i)==0 && touying(j+1,i)~=0)
            for k=1:52
                if(touying(j+k,i)~=0 && touying(j+1+k,i)==0)
                    break;
                end
            end
            %k为该完整行的高度
            if(k>=23 && k<=27)%对于小字母k+j为要记录的底部
                biaozhun(i)=k+j;
            elseif(k>=31 && k<=42)
                biaozhun(i)=k+j;
            elseif(k>=45 && k<=54)
                biaozhun(i)=k+j-14;
            end
            if(k<23)
                continue;
            end
            break;
        end
    end
end
temp=[4	8	11	29	79	102	107	136	142	146	157	159	161	163	168	174	175	199	204	222	240	259	270	272	278	289	295	306	307	312	313	322	330	336	342 347	363	377	378		386	399	400	405	410	413	417];
[t,n]=size(temp);
for i=1:n
    biaozhun(temp(i))=biaozhun(temp(i))-14;
end
biaozhun(22)=84;%逗号
biaozhun(11)=84;%?
biaozhun(275)=(96)
biaozhun(69)=107;%下面半个
biaozhun(127)=86;
biaozhun(235)=86;
biaozhun(271)=33;%18
biaozhun(279)=55;%24
biaozhun(305)=84;%23
biaozhun(55)=32;%45
biaozhun(108)=76;%23
biaozhun(225)=55;%27
biaozhun(255)=85;%25
biaozhun(132)=53;%26
biaozhun(139)=95;%38
biaozhun(213)=55;%41
biaozhun(242)=84;%34
biaozhun(282)=74;%40
biaozhun(304)=53;%19
biaozhun(379)=75;%28
biaozhun(126)=109;%15
biaozhun(135)=85;
biaozhun(285)=63;%32
biaozhun(405)=76;%19?
%biaozhun(74)=biaozhun(74)-12+14;
%biaozhun(93)=biaozhun(93)-11;
%biaozhun(105)=biaozhun(105)-11;
%biaozhun(173)=89;%逗号
%biaozhun(95)=86;%逗号
%[0 2 22 25 33 34 62 40 74 102 110 113 140 151 167 169 173 201 208]-14
%!!73-12 104-11 92-11
%73
%寻找每行开头5列像素的255
diyipian=zeros(1,22);
k=1;
for i=1:num
    if(sum(sum(img0(:,i*72-71:i*72-60)))==255*180*12)
        diyipian(k)=i;
        k=k+1;
    end
end
%减64
for i=1:num
    if biaozhun(i)>66
        biaozhun(i)=biaozhun(i)-64*floor(biaozhun(i)/64)
    end
end

%计算每块第一块黑色区域的下边界
%65为行间距
%for i=1:num
%    if(biaozhun(i)-65*floor(biaozhun(i)/65)>11)
%        biaozhun(i)=biaozhun(i)-65*floor(biaozhun(i)/65);
%    end
%end

biaozhunpaixu=sort(biaozhun);


完整代码 https://github.com/yan-fanyu/CUMCM-Paper-And-SourceCode

%第一块黑色的下边缘
%{
for i=1:num
    for j=1:179
        if(touying(j,i)-touying(j+1,i)>=9 && touying(j+1,i)<14)
            biaozhun(i)=j;
            break;
        end
    end
end
ans=sort(biaozhun);
%}


%像素数直方图
num=209;

for i=1:num  %读取所有碎纸片
    if i<11
        [img0(:,i*72-71:i*72),cmap0(:,i*3-2:i*3)]=imread(strcat('00',num2str(i-1)),'bmp');
    elseif i<101
        [img0(:,i*72-71:i*72),cmap0(:,i*3-2:i*3)]=imread(strcat('0',num2str(i-1)),'bmp');
    else
        [img0(:,i*72-71:i*72),cmap0(:,i*3-2:i*3)]=imread(num2str(i-1),'bmp');
    end
end
touying=zeros(180,209);%touying用每一列记录对应图片黑色像素点数
for i=1:209
    for j=1:180
        touying(j,i)=length(find(img0(j,i*72-71:i*72)~=255));
    end
end


%计算投影
num=418;

for j=1:num  %读取所有碎纸片
    if mod(j,2)==1
        si='a';
    else
        si='b';
    end
    i=ceil(j/2);
    if i<11
        [img0(:,j*72-71:j*72),cmap0(:,j*3-2:j*3)]=imread(['00',num2str(i-1),si],'bmp');
    elseif i<101
        [img0(:,j*72-71:j*72),cmap0(:,j*3-2:j*3)]=imread(['0',num2str(i-1),si],'bmp');
    else
        [img0(:,j*72-71:j*72),cmap0(:,j*3-2:j*3)]=imread([num2str(i-1),si],'bmp');
    end
end
touying=zeros(180,num);%touying用每一列记录对应图片黑色像素点数
for i=1:num
    for j=1:180
        touying(j,i)=length(find(img0(j,i*72-71:i*72)~=255));
    end
end
%image(img0(:,1:72));
%colormap(cmap0(:,1:3));


num=19;
sign=hang(i);
for i=1:num  %读取所有碎纸片
    if sign(i)<11
        [img0(:,i*72-71:i*72),cmap0(:,i*3-2:i*3)]=imread(strcat('00',num2str(sign(i)-1)),'bmp');
    elseif sign(i)<101
        [img0(:,i*72-71:i*72),cmap0(:,i*3-2:i*3)]=imread(strcat('0',num2str(sign(i)-1)),'bmp');
    else
        [img0(:,i*72-71:i*72),cmap0(:,i*3-2:i*3)]=imread(num2str(sign(i)-1),'bmp');
    end
end
img1=im2double(img0)*255;

for i=1:19
    for j=1:19 %deta255(i,j)i的右侧和j的左侧的M距离
        deta255(i,j)=sum(abs(img1(:,i*72)-img1(:,j*72-71)));
    end
end
%人为干预
%deta255(13,7)=inf;deta255(13,6)=inf;deta255(8,7)=inf;
deta255(10,2)=0
now=1;
img(:,1:72)=img0(:,now*72-71:now*72);
shunxu(1)=sign(now);
now_temp(1)=now;
for i=2:num
    deta255(:,now)=inf;
    now=find(deta255(now,:)==min(deta255(now,:)))
    now_temp(i)=now;
    shunxu(i)=sign(now);
    img(:,i*72-71:i*72)=img0(:,now*72-71:now*72);
end
figure;
image(img);
colormap(cmap0(:,1:3));
hold on;
for i=1:num
    text((i-1)*72,150,int2str(now_temp(i)));
end
for i=1:num
    plot([72*i-1,72*i-1],[0,180])
end
%img0(:,1:72:1368)=0;

完整代码 https://github.com/yan-fanyu/CUMCM-Paper-And-SourceCode

img_gehang=zeros(180*11,72*19);
img0=zeros(180*11,72*19);
cmap0=zeros(180*11,72*19);
for i=1:11
for j=1:19  %读取所有碎纸片
    if gehang(i,j)<11
        [img0(180*i-179:180*i,j*72-71:j*72),cmap0(256*i-255:256*i,j*3-2:j*3)]=imread(strcat('00',num2str(gehang(i,j)-1)),'bmp');
    elseif gehang(i,j)<101
        [img0(180*i-179:180*i,j*72-71:j*72),cmap0(256*i-255:256*i,j*3-2:j*3)]=imread(strcat('0',num2str(gehang(i,j)-1)),'bmp');
    else
        [img0(180*i-179:180*i,j*72-71:j*72),cmap0(256*i-255:256*i,j*3-2:j*3)]=imread(num2str(gehang(i,j)-1),'bmp');
    end
end
end
image(img0);
colormap(cmap0(1:256,1:3));
touying=zeros(1980,1);
touying(find(sum(img0,2)==1368*255))=255;
hold on
axis off
for i=1:11
plot([0 1368],[i*180 i*180])
end
for i=1:180
    for j=1:11
        touying1(i,j)=touying((j-1)*180+i);
    end
end
for i=1:19
plot([i*72 i*72],[0 1980])
end


%对于纵向剪裁的纸片的拼接，适用于中、英文
img0=zeros(1980,1368);cmap0=zeros(256,19*3);img=zeros(1980,1368);
for i=1:19%读取所有碎纸片
    [img0(:,i*72-71:i*72),cmap0(:,i*3-2:i*3)]=imread(num2str(i-1),'bmp');
end

deta255=zeros(19);
for i=1:19
    for j=1:19 %deta255(i,j)i的右侧和j的左侧
        deta255(i,j)=abs(sum(abs((img0(:,i*72)-img0(:,j*72-71)))));
    end
end
img(:,1:72)=img0(:,9*72-71:9*72);
now=9;%第一片序号
for i=2:19
    now=find(deta255(now,:)==min(deta255(now,:)));
    img(:,i*72-71:i*72)=img0(:,now*72-71:now*72);
end
image(img);
colormap(cmap0(:,1:3));
hold on
num=19
for i=1:num
    text((i-1)*72,1980,int2str(now_temp(i)-1));
end
for i=1:num
    plot([72*i-1,72*i-1],[0,1980])
end
axis off
完整代码 https://github.com/yan-fanyu/CUMCM-Paper-And-SourceCode

四、完整word论文和源代码

请添加图片描述

2013 数学建模国赛 b题（详细易懂版）
qq_22655703的博客
08-27 2415
本文主要围绕碎纸片拼接问题，建立了边界。
2013数学建模国赛B题思考过程和模型还有代码（适合初学者）
09-17
这个是2013数学建模国赛B题的一个完整的解答过程，包括了对题目的思考，模型的建立，代码的编写，论文的写法。当然还包括了最终的结果-----完整的报告。
2013高教社杯全国大学生数学建模竞赛题目B题（含附件）
09-13
2013高教社杯全国大学生数学建模竞赛题目B题（含附件）
数学建模2013年国赛B题碎纸片复原（纵切和横纵切两问）
06-20
数字图像大作业，效果感觉一般，做个记录。第二问效果不是很好，横纵切英文附件行间归类需重新完善
Python实现图片拼接的代码
09-20
本文通过实例代码给大家介绍了python实现图片拼接的方法，非常不错，具有一定的参考借鉴借鉴价值，需要的朋友参考下吧
2013年全国大学数学建模竞赛B题碎纸片复原（word论文+matlab源代码）
最新发布
06-17
2013年全国大学数学建模竞赛B；题碎纸片复原word论文；matlab源代码；国家一等奖；国家二等奖；可编辑word论文文档；包含数据处理过程；碎纸片的拼接复原方法研究 摘 要 本文研究了不同切割方式下中、英文碎纸片的拼接复原方法。文中充分利用不同图片边缘数据信息的匹配度、行间距及字高等文字排版数据建立了一系列拼接复原模型。 其中对于纵横切碎片复原问题，按照先找各裁剪行的组成碎片，然后实现行内拼接，最后对裁剪行之间的拼接顺序进行求解。 对仅纵切方式下的中、英文碎纸片的复原问题，本文建立了基于边缘列灰度值向量相似性的图像拼接复原模型：提取图像矩阵的左、右边缘灰度列向量，引入Minkowski距离，来衡量一个图像右边缘和另一个左边缘的匹配度，依据最佳匹配进行拼接。将仅纵切方式下中、英文碎片的复原结果列入表1和表2中，将复原结果图列入附录5和附录8中。 在纵向且横向切割方式下的单面图片复原问题中，引入分治的思想，缩小匹配范围，解决了由于碎片长度减小导致的碎片匹配度低的问题。首先对中、英文分别建立了基于像素点投影和灰度计数直方图的分类方式，通过投影和灰度计数直方图将二维的字体转化为一维的像
2013年全国大学生数学建模B题代码(不大)
09-16
针对B题碎片拼接的matlab代码和文献 我能帮的只有这么多了，满意还望采纳
2013数学建模B题
qq_43966202的博客
09-01 3051
参考链接：https://wenku.baidu.com/view/08647684ce2f0066f53322c7.html 思路 针对附件一： 根据汉字的矩形特点，41×41，记录下边界残缺汉字宽度进行匹配，越接近于41匹配度就越高。（首先人工干预找出最左边碎片） 针对附件二： 由于英文不具备汉字矩形特点，根据边缘像素点值的相近程度进行匹配。 针对附件三： 先分行，再对行进行排序。 根据上边缘...
2013年全国大学生数学建模竞赛b题
09-13
2013年全国大学生数学建模竞赛b题的题目，希望能对你有所帮助。
（数学建模）2013年国赛B题-碎纸片复原python代码
热门推荐
weixin_47696370的博客
06-20 1万+
数字图像处理的期末大作业效成绩出来了，感觉一般般，做个记录目录，对应到数学建模2013年国赛B题
2013数学建模B题算法
09-14
数学建模必备 碎纸片拼接 碎纸片复原 2013数学建模
2013年全国大学生数学建模B题matlab代码
09-16
请下载附件参考吧，里面有8个文件，针对B题碎片拼接的matlab代码和文献 我能帮的只有这么多了，满意还望采纳
2013数学建模B题附件
09-13
建模必备，2013年参加国赛的同学做B题一定要用到的材料
2013年数学建模B题问题一代码
09-16
2013年数学建模B题附件1代码
2013数学建模B题答案（附图）
09-15
2. 对于碎纸机既纵切又横切的情形，请设计碎纸片拼接复原模型和算法，并针对附件3、附件4给出的中、英文各一页文件的碎片数据进行拼接复原。如果复原过程需要人工干预，请写出干预方式及干预的时间节点。复原结果表达要求同上。 3. 上述所给碎片数据均为单面打印文件，从现实情形出发，还可能有双面打印文件的碎纸片拼接复原问题需要解决。附件5给出的是一页英文印刷文字双面打印文件的碎片数据。请尝试设计相应的碎纸片拼接复原模型与算法，并就附件5的碎片数据给出拼接复原结果，结果表达要求同上。
数学建模资料整理
shaozheng0503的博客
03-11 2031
（4）学习相关软件。在欣赏优秀论文的过程中，我十分注意模型的适用范围，举个例子来说，对于预测类的题目，比较常用的预测模型有时间序列模型、灰色预测模型、贝叶斯预测模型、神经网络预测模型等，这些模型并不是对所有的数据都是适用的，有些模型需要先对数据进行剔除、平均等处理，这些细节需要特别注意，一旦不注意就会影响整篇论文的质量。因为数学建模本身就是一个学习的过程，在短短3天时间里，将陌生的知识转化成自己的知识是具有挑战的，更何况还要对模型进行改进，但是正是这样，我们才能不断接触新知识，不断培养自己的学习能力。
2013全国大学生数学建模B题matlab代码
09-16
2013全国大学生数学建模B题matlab代码
2013年国赛碎纸片b题
05-15
2013年全国信息学奥林匹克竞赛（以下简称国赛）碎纸片B题是一道考察递推算法的题目。题目主要考查学生对递推算法的掌握及思维能力。 题目大意是：给出一张 $n\times m$ 的纸片，每个位置上初始写有1或0，依次进行 $t$ 次操作，每次操作选择一个初始为1的位置，将这个位置的1撕掉并在其上、下、左、右四个方向相邻位置加上1，如果该位置上原本是0，则不执行操作。最后统计剩余纸片上初始为1的位置数目。 对于这道题目，首先需要建立模型，将题目的要求转化为数学模型，进而进行求解。建立模型时，需要仔细思考，设计合理的递推式，同时需要保证递推式的正确性。在本题中，需要注意到纸片上的点在进行一次撕纸操作时，只会影响到相邻的四个位置。因此我们可以设置一个 $n\times m$ 的数组存储每个位置的值，然后按照题目要求进行 $t$ 次操作。每次选取初始位置后，通过递推式进行处理，同时需要注意判断位置是否越界，以及位置的初始值是否为1。最后统计数组中初始值为1的位置数目即为答案。 在解题过程中，需要分析算法的时间复杂度和空间复杂度。在此题中，时间复杂度为 $O(nmt)$，空间复杂度为 $O(nm)$。因此，需要在保证正确性的前提下，尽可能优化算法，减小时间和空间的占用。 对于本题，需要注意一些细节问题，如初始纸片上元素为0时的统计，纸片尺寸变化对边界的影响等。经过认真分析和细心实现，就能得到正确的答案。 总之，本题难度较大，需要考生具有一定的数学建模和编程能力。通过这种方式锻炼了学生的思维能力和算法水平，为科学研究和工程实践提供了有力的支撑。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

亦是远方

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值