.clc clear all% 读入图片pic1=imread( lena1. jpg); pic2=imread( lena2. jpg);% Harris角点检测points1=myHarris( pic1);points2=myHarris( pic2);% 画出Harris 角点figure( 1)drawHarrisCorner( pic1, points1, pic2, points2);% 角点特征描述des1=myHarrisCornerDescription ( pic1, points1); des2=myHarrisCornerDescription ( pic2, points2);% 角点粗匹配matchs=myMatch ( des1, des2);% 获取各自出匹配角点位置matchedPoints1=points1( matchs(:,1),:);matchedPoints2=points2( matchs(:,2),:);% 粗匹配角点连线figure( 2)drawLinedC
