MATLAB匹配按列碎开的纸片

纵向切开的纸片的拼接

%% paper_cutting.m

clear, close all

clc

%% import pictures, and save into images{img_num}

[images, img_num, vertical_border] = import_data('.\img_data\pic1\');    % 批量读入图片, 在上上篇中写到.'

%% 求出任意一纸片的左侧和另一任意纸片的右侧的距离roads

roads = zeros(img_num);

header = 0;

tailer = 0;

for i = 1 : img_num

    for j = 1 : img_num

        if i == j

            roads(i, j) = 1000;

            continue

        end

        roads(i, j) = sum(abs(images{i}(:, end) - images{j}(:, 1)));

        if roads(i, j) == 0

            header = j;

            tailer = i;

        end

    end

end

fraction_idx = zeros(img_num, 1);

fraction_idx(1) = header;

fraction_idx(end) = tailer;

for i = 1 : img_num

    if count_white_chunk(images, i, 50)

        fraction_idx(2) = i;

    end

end

%% xiongyali method:

c = roads;

[h, w] = size(c);

c = c'; c = c(:);

Aeq = zeros(h+w, h*w);

for i = 1 : h

    Aeq(i, (i-1) * w + 1 : w*i) = 1;

end

for i=1 : w

    Aeq(h + i, i : w : h*w) = 1;

end

beq = ones(h+w, 1);

A = []; b = [];

[x,fval] = bintprog(c, A, b, Aeq, beq);

X = reshape(x, w, h)';

%% Get the order.

for i = 2 : img_num-1

    fraction_idx(i) = find(X(fraction_idx(i-1), :) ~= 0);

end

%% connect the fractions and show it.

for i = 1 : img_num

    if i == 1

        connected_img = images{fraction_idx(i)};

    else

        connected_img = [connected_img, images{fraction_idx(i)}];

    end

end

figure(1);

imshow(connected_img);

imwrite(connected_img, 'saved1.png', 'png');

title('Connected fractions.');