基于灰度投影的图像配准MATLAB

MATLAB代码，基于灰度在水平和竖直方向进行投影计算累计像素值，得到两条曲线，根据两条曲线可以进行图像配准。

clear; clc;
% result:
% vertical: 52
% horizon: 48

%% read img data
img1 = imread（‘one.jpg‘）;
img2 = imread（‘two.jpg‘）;
img1 = rgb2gray（img1）;
img2 = rgb2gray（img2）;

[row column] = size（img1）;

%% project the src（source） img in horizon and vertical
horizonImg1 = zeros（1row）;
horizonImg2 = zeros（1row）;
verticalImg1 = zeros（1 column）;
verticalImg2 = zeros（1 column）;

for i=1:row
    horizonImg1（1i） = sum（img1（i:））;
    horizonImg2（1i） = sum（img2（i:））;
end

for i=1:column
    verticalImg1（1i） = sum（img1（:i））;
    verticalImg2（1i） = sum（img2（:i））;
end

%% plot the project of src img
figure
% subplot（211）
verticalImg1_projection_src = verticalImg1;
plot（verticalImg1_projection_src‘-or‘‘markersize‘3）
hold on
% title（‘vertical-1‘）
% subplot（212）
verticalImg2_projection_src = verticalImg2;
plot（verticalImg2_projection_src‘-ob‘‘markersize‘3）
% title（‘vertical-2‘）
title（‘vertical - src‘）
xlabel（‘方向‘）
ylabel（‘the sum of pixel value‘）

figure
% subplot（211）
horizonImg1_projection_src = horizonImg1;
plot（horizonImg1_projection_src‘-or‘‘markersize‘3）
hold on;
% title（‘horizon-1‘）
% subplot（212）
horizonImg2_projection_src = horizonImg2;
plot（horizonImg2_projection_src‘-ob‘‘markersize‘3）
% title（‘horizon-2‘）
title（‘horizon - src‘）
xlabel（‘方向‘）
ylabel（‘the sum of pixel value‘）

%% compute the error in horizon the image 1 is referenced
region = 200; % 匹配总范围
stepHorizon = 2; % 步长
start = 2; % 匹配起点
errorHorizon = inf * ones（1 floor（region-start）/stepHorizon+1）;
for position=start:stepHorizon:region
    horizonImg1_ = horizonImg1（11:row - position + 1）;
    horizonImg2_ = horizonImg2（1position:row）;
    errorHorizon（1（position-start）/stepHorizon+1） = sum（abs（horizonImg1_ - horizonImg2_））;
end
minErrorHorizon = min（errorHorizon）;
minErrorIndexHorizon = stepHorizon * （find（errorHorizon == minErrorHorizon） - start） + start; % may be not only one

%% compute the error in vertical the image 1 is referenced too
regionVertical = 200;
stepVertical = 2;
start = 2;
errorVertical = inf * ones（1 floor（regionVertical - start）/stepVertical + 1）;
for position=start:stepVertical:regionVertical
    verticalImg1_ = verticalImg1（11:column-position+1）;
    verticalImg2_ = verticalImg2（1position:column）;
    errorVertical（1 （position - start） / stepVertical + 1） = sum（abs（verticalImg2_ - verticalImg1_））;
end
minErrorVertical = min（errorVertical）;
minErrorIndexVertical = stepVertical * （find（errorVertical == minErrorVertical） - start） + start; % may be not only one

%% change the no-reference image
img2 = imread（‘two.jpg‘）;
img2 = img2（minErrorIndexHorizon（1）:end minErrorIndexVertical（1）:end:）;
figureimshow（img2）
title（‘align - img2‘）

img1 = imread（‘one.jpg‘）;
img1 = img1（1:row-minErrorIndexHorizon 1:column-minErrorIndexVertical :）;
figure imshow（img1）;
title（‘align - img1‘）

[row column] = size（img1）;
horizonImg1 = zeros（1row）;
horizonImg2 =

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-07-31 23:02  灰度投影配准图像-MATLAB\
     文件       22017  2018-07-30 22:08  灰度投影配准图像-MATLAB\1.jpg
     文件       21100  2018-07-30 22:08  灰度投影配准图像-MATLAB\22.jpg
     文件        3824  2018-07-31 22:47  灰度投影配准图像-MATLAB\alignTwoImages.m
     文件       23236  2018-07-31 22:40  灰度投影配准图像-MATLAB\align_h.jpg
     文件       19591  2018-07-31 22:47  灰度投影配准图像-MATLAB\align_img1.jpg
     文件       18965  2018-07-31 22:47  灰度投影配准图像-MATLAB\align_img2.jpg
     文件       32425  2018-07-31 22:39  灰度投影配准图像-MATLAB\align_v.jpg
     文件       48272  2018-07-30 22:06  灰度投影配准图像-MATLAB\one.jpg
     文件       23590  2018-07-31 22:39  灰度投影配准图像-MATLAB\src_h.jpg
     文件       26420  2018-07-31 22:51  灰度投影配准图像-MATLAB\src_v.jpg
     文件       45258  2018-07-30 22:07  灰度投影配准图像-MATLAB\two.jpg