最小二乘法匹配，matlab代码，相关系数最大，单点最小二乘匹配.zip

最小二乘法匹配，matlab代码，相关系数最大，单点最小二乘匹配

%---------------------------------------------%
%					      %
%          工作室提供代做matlab仿真	      %
%					      %
%  详情请访问：http://cn.mikecrm.com/5k6v1DP  %
%					      %
%---------------------------------------------%

clc;
clear;

%% 数据读取
imglef = imread（‘LOR50.bmp‘）; imgrig = imread（‘LOR49.bmp‘）;
[rowlef collef] = size（imglef）;
[rowrig colrig] = size（imgrig）;

% 选点
figure（1）;
imshow（imglef）;
title（‘image on the left‘）;
hold on
[xpick ypick] = ginput（1）;
xpick = fix（xpick）;
ypick = fix（ypick）;
hold on
plot（xpick ypick ‘y+‘）;

% 设置目标区大小，计算目标区步长 
n = 9; step_n = （n-1）/2;
% 构建选点目标区
winobjlef = xpick - step_n; winobjrig = xpick + step_n;
winobjbot = ypick - step_n; winobjtop = ypick + step_n;
% 忽略边缘点
if winobjlef>1 && winobjrig1 && winobjtop    winobj = imglef（winobjbot : winobjtop winobjlef : winobjrig）;
end

% 计算选点最邻近角点
corpt = corner（winobj ‘Harris‘）;
corpt（: 1） = corpt（: 1） - （step_n+1） + xpick;
corpt（: 2） = corpt（: 2） - （step_n+1） + ypick;
plot（corpt（: 1） corpt（: 2） ‘r+‘）;
[numcorpt ~] = size（corpt）;
if numcorpt == 1
    corptnear = [corpt（: 1） corpt（: 2）];
elseif numcorpt > 1
   distsquare = （corpt（: 1）-xpick）.*（corpt（: 1）-xpick）...
              + （corpt（: 2）-ypick）.*（corpt（: 2）-ypick）;
       rowmin = find（distsquare == min（distsquare））;
    corptnear = [corpt（rowmin 1） corpt（rowmin 2）];
end

%% 相关系数法进行影像粗匹配
% 构建角点目标区
winobjlef = corptnear（1） - step_n; winobjrig = corptnear（1） + step_n;
winobjbot = corptnear（2） - step_n; winobjtop = corptnear（2） + step_n;
% 忽略边缘点
if winobjlef>1 && winobjrig1 && winobjtop    winobj = imglef（winobjbot : winobjtop winobjlef : winobjrig）;
end

% 根据左右片特点构建搜索区
rho = zeros（size（winobj））;
for x = step_n+1 : fix（0.7*colrig）-step_n
    for y = step_n+1 : rowrig-step_n
        % 在搜索区中构建同名窗口
        winrig = imgrig（y-step_n : y+step_n x-step_n : x+step_n）;
        % 计算同名窗口与目标窗口的相关系数
        rho（xy） = corr2（winobj winrig）;
    end
end

rhomax = max（max（rho））;
[xmax ymax] = find（rho == rhomax）;
winrig = imgrig（ymax-step_n : ymax+step_n xmax-step_n : xmax+step_n）;

%% 最小二乘法进行影像精匹配
% 赋参数初值
h0 = 0; h1 = 1;
a0 = 0; a1 = 1; a2 = 0;
b0 = 0; b1 = 0; b2 = 1;

% 中心化
g1 = double（winobj） - mean2（winobj）;
g2 = double（winrig） - mean2（winrig）;
g1in = g1（2 : n-1 2 : n-1）; 
g2in = g2（2 : n-1 2 : n-1）;
rho = abs（corr2（g1 g2））;

iter = 0; flag = 1;
while （flag）
    % 构建误差方程
    countpixel = 1; C = zeros（（n-2）*（n-2） 8）; L = zeros（（n-2）*（n-2） 1）;
    for i = 2 : n-1
        for j = 2 : n-1
            % 计算影像灰度梯度值
            gradx = 0.5 * （g2（i   j+1） - g2（i   j-1））;
            grady = 0.5 * （g2（i+1 j  ） - g2（i-1 j  ））;
            % 计算系数阵与常数阵
            c1 = 1;   c2 = g2（i j）;
            c3 = double（gradx）; c4 = i * double（gradx）; c5 = j * double（gradx）;
            c6 = double（grady）; c7 = i * double（grady）; c8 = j * double（grady）;
         

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-05-09 16:37  最小二乘法匹配，matlab代码，相关系数最大，单点最小二乘匹配\
     文件        5701  2019-05-09 16:39  最小二乘法匹配，matlab代码，相关系数最大，单点最小二乘匹配\m.m