91av视频/亚洲h视频/操亚洲美女/外国一级黄色毛片 - 国产三级三级三级三级

  • 大小: 3KB
    文件類型: .zip
    金幣: 2
    下載: 0 次
    發布日期: 2021-05-16
  • 語言: Matlab
  • 標簽: matlab??

資源簡介

最小二乘法匹配,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

評論

共有 條評論