歸一化互相關(guān)

實(shí)訓(xùn)題目用歸一化互相關(guān)注冊(cè)一幅圖像指導(dǎo)老師盛利元專(zhuān)業(yè)班級(jí)電信0901姓名耿寶印學(xué)號(hào) 1404090114handletostart(seeGCBO)reservedhandletostart(seeGCBO)reserved-tobedefinedinafutureversionstructurewithhandlesanduserdata(seeRegisteringanImageUsingNormalizedCross-Correlation這次實(shí)訓(xùn)，是在參考matlab樣題后進(jìn)行的。剛接觸這個(gè)題目時(shí)，自己對(duì)什么是互相關(guān)也不知道，經(jīng)過(guò)查資料，相關(guān)說(shuō)簡(jiǎn)單點(diǎn)，就是相似性。構(gòu)思；由一及多。選擇兩張?zhí)囟ǖ膱D片，大致選定區(qū)域，然后用normxcorr2判定相關(guān)性，找出補(bǔ)償點(diǎn)，對(duì)圖片進(jìn)行相似處理。最后進(jìn)行圖片擴(kuò)展。實(shí)現(xiàn)盡可能多的圖片的互相關(guān)的注冊(cè)。具體步驟以樣題代碼為基礎(chǔ)，在matlab自帶guide平臺(tái)進(jìn)行圖形設(shè)計(jì)，建立xiangguan.flg文件并相應(yīng)產(chǎn)生xiangguan.m文件。在設(shè)計(jì)好大體的圖形界面后，對(duì)xiangguan.m文件進(jìn)行代碼編寫(xiě)。t最終guide平臺(tái)下的圖形界面如下：為優(yōu)化代碼使用全局變量，并在相應(yīng)的地方進(jìn)行聲明。t具體控件代碼如下：1)選圖按鈕Firstpicfunctionstart_Callback(hObject,eventdata,handles)%hObject%eventdataofMATLAB%handlesGUIDATA)subplot(3,3,2)[filename,pathname]=uigetfile({'*.jpg;*.tif;*.png;*.gif','AllImageFiles';...'*.*','AllFiles'},'open'); %打開(kāi)路徑下要檢索的圖像ifisequal([filename,pathname],[0,0])returnelse%讀取圖片pic=fullfile(pathname,filename);globalonion;onion=imread(pic);imshow(onion);end% Executesonbuttonpressinnext.functionnext_Callback(hObject,eventdata,handles)(seeGCBO)bedefinedinafutureversionofMATLABhandlesanduserdata(seeGUIDATA)[filename,pathname]=uigetfile({'*.jpg;*.tif;*.png;*.gif',,AllImageFiles''*.*','AllFiles'},'open'); %打開(kāi)路徑下要檢索的圖像ifisequal([filename,pathname],[0,0])returnelse%讀取圖片pict=fullfile(pathname,filename);globalpeppers;peppers=imread(pict);imshow(peppers);end% Executesonbuttonpressinrestart.functionrestart_Callback(hObject,eventdata,handles)%hObjecthandletorestart(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)相關(guān)比較按鈕Correlationsubplot(3,3,5)globalonion;globalpeppers;%non-interactivelyrect_onion=[111336558];rect_peppers=[16347143151];sub_onion=imcrop(onion,rect_onion);sub_peppers=imcrop(peppers,rect_peppers);c=normxcorr2(sub_onion(:,:,1),sub_peppers(:,:,1));surf(c),shadingflat%---Executesonbuttonpressinoverlay.functionoverlay_Callback(hObject,eventdata,handles)%hObjecthandletooverlay(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)按鈕Overlaysubplot(3,3,6)globalonion;globalpeppers;%non-interactivelyrect_onion=[111336558];rect_peppers=[16347143151];sub_onion=imcrop(onion,rect_onion);sub_peppers=imcrop(peppers,rect_peppers);c=normxcorr2(sub_onion(:,:,1),sub_peppers(:,:,1));%offsetfoundbycorrelation[max_c,imax]=max(abs(c(:)));[ypeak,xpeak]=ind2sub(size(c),imax(1));corr_offset=[(xpeak-size(sub_onion,2))(ypeak-size(sub_onion,1))];%relativeoffsetofpositionofsubimagesrect_offset=[(rect_peppers(1)-rect_onion(1))(rect_peppers(2)-rect_onion(2))];%totaloffsetoffset=corr_offset+rect_offset;xoffset=offset(1);yoffset=offset(2);xbegin=round(xoffset+1);xend=round(xoffset+size(onion,2));ybegin=round(yoffset+1);yend=round(yoffset+size(onion,1));%extractregionfrompeppersandcomparetoonionextracted_onion=peppers(ybegin:yend,xbegin:xend,:);ifisequal(onion,extracted_onion)disp('onion.pngwasextractedfrompeppers.png')endrecovered_onion=uint8(zeros(size(peppers)));recovered_onion(ybegin:yend,xbegin:xend,:)=onion;imshow(recovered_onion)% Executesonbuttonpressintransparent.functiontransparent_Callback(hObject,eventdata,handles)%hObjecthandletotransparent(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)5)透明處理按鈕Transparentsubplot(3,3,8)globalonion;globalpeppers;%non-interactivelyrect_onion=[111336558];rect_peppers=[16347143151];sub_onion=imcrop(onion,rect_onion);sub_peppers=imcrop(peppers,rect_peppers);c=normxcorr2(sub_onion(:,:,1),sub_peppers(:,:,1));%offsetfoundbycorrelation[max_c,imax]=max(abs(c(:)));[ypeak,xpeak]=ind2sub(size(c),imax(1));corr_offset=[(xpeak-size(sub_onion,2))(ypeak-size(sub_onion,1))];%relativeoffsetofpositionofsubimagesrect_offset=[(rect_peppers(1)-rect_onion(1))(rect_peppers(2)-rect_onion(2))];%totaloffsetoffset=corr_offset+rect_offset;xoffset=offset(1);yoffset=offset(2);xbegin=round(xoffset+1);xend=round(xoffset+size(onion,2));ybegin=round(yoffset+1);yend=round(yoffset+size(onion,1));%extractregionfrompeppersandcomparetoonionextracted_onion=peppers(ybegin:yend,xbegin:xend,:);ifisequal(onion,extracted_onion)disp('onion.pngwasextractedfrompeppers.png')endrecovered_onion=uint8(zeros(size(peppers)));recovered_onion(ybegin:yend,xbegin:xend,:)=onion;[m,n,p]=size(peppers);mask=ones(m,n);i=find(recovered_onion(:,:,1)==0);mask(i)=.2;%tryexperimentingwithdifferentlevelsof%transparency%overlayimageswithtransparencyimshow(peppers(:,:,1))%showonlyredplaneofpeppersholdonh=imshow(recovered_onion);%overlayrecovered_onionset(h,'AlphaData',mask)菜單quit% Executesonbuttonpressinquit.functionquit_Callback(hObject,eventdata,handles)%hObjecthandletoquit(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)button=questdlg('readytoquit?','exitdialog','保存退出','關(guān)閉運(yùn)行窗口','取消')switchbuttoncase'保存退出'

saveexitcase,取消'quitcase'關(guān)閉運(yùn)行窗口’closeend為了在同一