资源简介
最小二乘法匹配,matlab代码,相关系数最大,单点最小二乘匹配
![](http://www.nz998.com/pic/75080.jpg)
代码片段和文件信息
%---------------------------------------------%
% %
% 工作室提供代做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
相关资源
- Matlab路面裂缝识别69319
- 高灵敏度GPS接收机MATLAB仿真,附捕获
- 基于MATLAB的质点弹道计算与外弹道优
- 阵列天线的matlab仿真
- MATLAB 经典程序源代码大全
- MATLAB小波软阈值去噪代码33473
- 天线阵的波束形成在MATLAB仿真程序及
- 非线性SVM算法-matlab实现
- 《MATLAB 智能算法超级学习手册》-程序
- 组合导航matlab程序
- 读取txt文件内容matlab代码实现
- Matlab实现基于相关的模板匹配程序
- matlab优化工具箱讲解
- 基于MATLAB的快速傅里叶变换
- 光纤传输中的分布傅立叶算法matlab实
- 基于matlab的图像处理源程序
- matlab 椭圆拟合程序
- 算术编码解码matlab源代码
- optical_flow 光流法 matlab 实现程序
- 引导图像滤波器 Matlab实现
- 分形几何中一些经典图形的Matlab画法
- OFDM系统MATLAB仿真代码
- SVM工具箱(matlab中运行)
- 图像小波变换MatLab源代码
- LU分解的MATLAB实现
- 冈萨雷斯数字图像处理matlab版(第三
- 替代数据法的matlab程序
- 用matlab实现的多站定位系统性能仿真
- 通过不同方法进行粗糙集属性约简m
- k近邻算法matlab实现
评论
共有 条评论