MATLAB图像拼接——怎么用MATLAB做拼图？

如下图所示，共有120个子图，如何用MATLAB拼接起来？

1 读取数据

首先读取所有图像数据，存在X中，并用imshow显示。

clear all
%% 读取图片数据并显示
X = {
   };
for i = 1:120
    X{
   i} = imread(['G:\图像处理\matlab\拼接图像\附件2-一般图片\00',num2str(i),'.bmp']);
end
for i=1:120
    subplot(12,10,i)
    imshow(X{
   i})
end

如图所示，这是需要拼接的120个子图

2 预处理

%% 转化为二值图
X_01 = {
   };
for i = 1:120
    X_01{
   i} = im2bw(rgb2gray(X{
   i}));
end
%示意
subplot(131);imshow(X{
   1})
subplot(132);imshow(rgb2gray(X{
   1}))
subplot(133);imshow(im2bw(rgb2gray(X{
   i})))

原图、rgb2gray变成的灰度图、以及二值图对比：

3 计算子图之间的相关性

只考虑边缘与边缘之间的关系。建立一个120x120的矩阵，里面分别存放着子图与子图之间的相关性。
如：如果图A左侧和图B右侧的二值接近，那么图A就在图B右边。

%%  计算分数
%score中保存着四种分数，分别对应两个图块的四种接法
score = zeros(120,120,4);
s=zeros(120,120,4);
for k = 1:4
    for i = 1:120        
        img1 = X_01{
   i};
        for j = 1:120
            if i~=j
                img2 = X_01{
   j};
                %s(i,j,k) = img12img2(img1,img2,k);
                if k==1
                    score(i,j,k) = sum((img1(:,1)-img2(:,end)).^2);
                elseif k==2
                    score(i,j,k) = sum((img1(1,:)-img2(end,:)).^2);
                elseif k==3
                    score(i,j,k) = sum((img1(:,end)-img2(:,1)).^2)