如下图所示,共有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