Question

求2013年数学建模国赛B题附件3拼接答案。题

如题。急。多多益善。

Answer 1

去交流群看看，我找到代码了，可是执行不了。
clear;

I=imread('附件1');

if(isgray(I)==0)
disp('请输入灰度图像,本程序用来处理128 *128的灰度图像!');
else
if (size(I)~=[128,128])
disp('图像的大小不合程序要求!');
else
H.color=[1 1 1]; %设置白的画布

figure(H);imshow(I);
title('原图像');
zeroImage=repmat(uint8(0),[128 128]);
figure(H); %为分裂合并后显示的图设置画布
meansImageHandle=imshow(zeroImage);
title('块均值图像');
%%%%%设置分裂后图像的大小由于本图采用了128像素的图
blockSize=[128 64 32 16 8 4 2];
%%设置一个S稀疏矩阵用于四叉树分解后存诸数据
S=uint8(128);
S(128,128)=0;
threshold=input('请输入分裂的阈值(0--1):');%阈值
threshold=round(255*threshold);
M=128;dim=128;

%%%%%%%%%%%%%%%%% 分裂主程序%%%%%%%%%%%

while (dim>1)
[M,N] = size(I);

Sind = find(S == dim);
numBlocks = length(Sind);
if (numBlocks == 0)
%已完成
break;
end
rows = (0:dim-1)';
cols = 0:M:(dim-1)*M;
rows = rows(:,ones(1,dim));
cols = cols(ones(dim,1),:);
ind = rows + cols;
ind = ind(:);

tmp = repmat(Sind', length(ind), 1);
ind = ind(:, ones(1,numBlocks));
ind = ind + tmp;

blockValues= I(ind);
blockValues = reshape(blockValues, [dim dim numBlocks]);

if(isempty(Sind))
%已完成
break;
end
[i,j]=find(S);
set(meansImageHandle,'CData',ComputeMeans(I,S));
maxValues=max(max(blockValues,[],1),[],2);
minValues=min(min(blockValues,[],1),[],2);
doSplit=(double(maxValues)-double(minValues))>threshold;
dim=dim/2;
Sind=Sind(doSplit);
Sind=[Sind;Sind+dim;(Sind+M*dim);(Sind+(M+1)*dim)];
S(Sind)=dim;
end
[i,j]=find(S); % 用来寻找四叉机分解结果中大小为S的块的位置
set(meansImageHandle,'CData',ComputeMeans(I,S)); % 显示分解结果块均值图像
Numberofbloks=length(i); %计算块数

%sizev=size(v);
end
end

程序2

其实很简单的，看懂下面的语句：
PicName1='a.jpg'; % 要合并的图片1

PicName2='b.jpg'; % 要合并的图片2

PicOut='c.jpg'; %合并的结果
IV1=imread(PicName1); % 读入图片1

IV2=imread(PicName2); % 读入图片2

PicData=[IV1; IV2]; % 如果是纵向拼接用这个
PicData=[IV1 IV2]; % 如果是横向拼接用这个
imwrite(PicData, PicOut, 'Quality', 75); % 输出图形\附件1\');

if(isgray(I)==0)
disp('请输入灰度图像,本程序用来处理128 *128的灰度图像!');
else
if (size(I)~=[128,128])
disp('图像的大小不合程序要求!');
else
H.color=[1 1 1]; %设置白的画布

figure(H);imshow(I);
title('原图像');
zeroImage=repmat(uint8(0),[128 128]);
figure(H); %为分裂合并后显示的图设置画布
meansImageHandle=imshow(zeroImage);
title('块均值图像');
%%%%%设置分裂后图像的大小由于本图采用了128像素的图
blockSize=[128 64 32 16 8 4 2];
%%设置一个S稀疏矩阵用于四叉树分解后存诸数据
S=uint8(128);
S(128,128)=0;
threshold=input('请输入分裂的阈值(0--1):');%阈值
threshold=round(255*threshold);
M=128;dim=128;

%%%%%%%%%%%%%%%%% 分裂主程序%%%%%%%%%%%

while (dim>1)
[M,N] = size(I);

Sind = find(S == dim);
numBlocks = length(Sind);
if (numBlocks == 0)
%已完成
break;
end
rows = (0:dim-1)';
cols = 0:M:(dim-1)*M;
rows = rows(:,ones(1,dim));
cols = cols(ones(dim,1),:);
ind = rows + cols;
ind = ind(:);

tmp = repmat(Sind', length(ind), 1);
ind = ind(:, ones(1,numBlocks));
ind = ind + tmp;

blockValues= I(ind);
blockValues = reshape(blockValues, [dim dim numBlocks]);

if(isempty(Sind))
%已完成
break;
end
[i,j]=find(S);
set(meansImageHandle,'CData',ComputeMeans(I,S));
maxValues=max(max(blockValues,[],1),[],2);
minValues=min(min(blockValues,[],1),[],2);
doSplit=(double(maxValues)-double(minValues))>threshold;
dim=dim/2;
Sind=Sind(doSplit);
Sind=[Sind;Sind+dim;(Sind+M*dim);(Sind+(M+1)*dim)];
S(Sind)=dim;
end
[i,j]=find(S); % 用来寻找四叉机分解结果中大小为S的块的位置
set(meansImageHandle,'CData',ComputeMeans(I,S)); % 显示分解结果块均值图像
Numberofbloks=length(i); %计算块数

%sizev=size(v);
end
end

程序2

其实很简单的，看懂下面的语句：
PicName1='a.jpg'; % 要合并的图片1

PicName2='b.jpg'; % 要合并的图片2

PicOut='c.jpg'; %合并的结果
IV1=imread(PicName1); % 读入图片1

IV2=imread(PicName2); % 读入图片2

PicData=[IV1; IV2]; % 如果是纵向拼接用这个
PicData=[IV1 IV2]; % 如果是横向拼接用这个
imwrite(PicData, PicOut, 'Quality', 75); % 输出图形
【希望采纳】

Answer 2

同学你好，我们已经记录了你的ip地址以及你的参赛学校，并向相关单位提交了查询。
根据《全国大学生数学建模竞赛参赛规则》第五条：“抄袭是严重违反竞赛规则的行为，有抄袭行为的参赛队在全国和赛区评阅时视为严重违反竞赛纪律；竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人，包括指导教师，研究及讨论与赛题有关的问题，否则也视为严重违反竞赛纪律。”

Answer 3

请下载附件参考吧，里面有8个文件，针对B题碎片拼接的matlab代码和文献

 

我能帮的只有这么多了，满意还望采纳

Answer 4

这个比较有用。

Answer 5

仅供参考，自己可以再斟酌下。