圖像拼接基于harris檢測(cè)matlab代碼-文庫吧
2025-07-08 21:52 本頁面
【正文】 Ixy=filter2(h,Ixy)。 [heigth,width]=size(pic)。 alpha=。 R=zeros(heigth,width)。 for i=1:heigth for j=1:width M=[Ix2(i,j) Ixy(i,j)。Ixy(i,j) Iy2(i,j)]。 R(i,j)=det(M)alpha*(trace(M)^2)。 end end Rmax=max(max(R))。 pMap=zeros(heigth,width)。 for i=2:heigth1 for j=2:width1 if R(i,j)*Rmax tm=R(i1:i+1,j1:j+1)。 tm(2,2)=0。 if R(i,j)tm pMap(i,j)=1。 end end end end [row,col]=find(pMap==1)。 points=[row,col]。 function drawHarrisCorner(pic1,points1,pic2,points2) % 功能：畫出Harris角點(diǎn)的連接% 輸入： % picpic2：待拼接的圖像 % pointspoints2：Harris角點(diǎn)位置 X1=points1(:,2)。 Y1=points1(:,1)。 X2=points2(:,2)。 Y2=points2(:,1)。 dif=size(pic1,2)。 imshowpair(pic1,pic2,39。montage39。)。 hold on plot(X1,Y1,39。b*39。)。 plot(X2+dif,Y2,39。b*39。)。 set(gcf,39。Color39。,39。w39。)。 function des=myHarrisCornerDescription(pic,points) % 功能：Harris角點(diǎn)特征描述 % 輸入： % pic：原圖像 % points：角點(diǎn)位置 % 輸出： % des：8N的角點(diǎn)特征描述矩陣 if length(size(pic))==3 pic=rgb2gray(pic)。 end len=length(points)。 des=zeros(8,len)。 for k=1:len p=points(k,:)。 pc=pic(p(1),p(2))。 des(1,k)=pic(p(1)1,p(2)1)pc。 des(2,k)=pic(p(1),p(2)1)pc。 des(3,k)=pic
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1