鑷畾涔夊嚱鏁伴兘鍦ㄤ笅闈?/p>
clc
clear all
%
璇诲叆鍥剧墖
pic
1
=
imread
('
lena
1
.
jpg
');
pic
2
=
imread
('
lena
2
.
jpg
');
%
Harris
瑙掔偣妫�娴?/p>
points
1
=
myHarris
(
pic
1
);
points
2
=
myHarris
(
pic
2
);
%
鐢诲嚭
Harris
瑙掔偣
figure
(
1
)
drawHarrisCorner
(
pic
1
,
points
1
,
pic
2
,
points
2
);
%
瑙掔偣鐗瑰緛鎻忚堪
des
1
=
myHarrisCornerDescription
(
pic
1
,
points
1
);
des
2
=
myHarrisCornerDescription
(
pic
2
,
points
2
);
%
瑙掔偣绮楀尮閰?/p>
matchs
=
myMatch
(
des
1
,
des
2
);
%
鑾峰彇鍚勮嚜鍑哄尮閰嶈鐐逛綅缃?/p>
matchedPoints
1
=
points
1
(
matchs
(:,
1
),:);
matchedPoints
2
=
points
2
(
matchs
(:,
2
),:);
%
绮楀尮閰嶈鐐硅繛绾?/p>
figure
(
2
)
drawLinedCorner
(
pic
1
,
matchedPoints
1
,
pic
2
,
matchedPoints
2
);
%
瑙掔偣绮惧尮閰?/p>
[
newLoc
1
,
newLoc
2
]=
pointsSelect
(
matchedPoints
1
,
matchedPoints
2
);
%
绮惧尮閰嶈鐐硅繛绾?/p>
figure
(
3
)
drawLinedCorner
(
pic
1
,
newLoc
1
,
pic
2
,
newLoc
2
);
%
鍥惧儚鎷兼帴
im
=
picMatched
(
pic
1
,
newLoc
1
,
pic
2
,
newLoc
2
);
%
鏄剧ず鎷兼帴鍥惧儚
figure
(
4
)
imshow
(
im
);
set
(
gcf
,'
Color
','
w
');
function points
=
myHarris
(
pic
)
%
鍔熻兘锛氬鎵?/p>
Harris
瑙掔偣
%
杈撳叆锛?/p>
RGB
鍥惧儚鎴?/p>
gray
鍥?/p>
%
杈撳嚭锛氳鐐规墍鍦ㄧ殑琛屻�佺旱鐨?/p>
N
脳
2
鐭╅樀
if length
(
size
(
pic
))==
3
pic
=
rgb
2
gray
(
pic
);
end
pic
=
double
(
pic
);
hx
=[-
1 0 1
];
Ix
=
filter
2
(
hx
,
pic
);