歸一化互相關(guān)

1、實(shí)訓(xùn)題目用歸一化互相關(guān)注冊一幅圖像指導(dǎo)老師盛利元專業(yè)班級電信0901姓名耿寶印學(xué)號1404090114RegisteringanImageUsingNormalizedCross-Correlation這次實(shí)訓(xùn)，是在參考matlab樣題后進(jìn)行的。剛接觸這個題目時，自己對什么是互相關(guān)也不知道，經(jīng)過查資料，相關(guān)說簡單點(diǎn)，就是相似性。構(gòu)思；由一及多。選擇兩張?zhí)囟ǖ膱D片，大致選定區(qū)域，然后用normxcorr2判定相關(guān)性，找出補(bǔ)償點(diǎn)，對圖片進(jìn)行相似處理。最后進(jìn)行圖片擴(kuò)展。實(shí)現(xiàn)盡可能多的圖片的互相關(guān)的注冊。具體步驟以樣題代碼為基礎(chǔ)，在matlab自帶guide平臺進(jìn)行圖形設(shè)計，建立xiangguan.f

2、lg文件并相應(yīng)產(chǎn)生xiangguan.m文件。在設(shè)計好大體的圖形界面后，對xiangguan.m文件進(jìn)行代碼編寫。最終guide平臺下的圖形界面如下：和Hiangguan口皿為優(yōu)化代碼使用全局變量，并在相應(yīng)的地方進(jìn)行聲明。具體控件代碼如下：1)選圖按鈕Firstpicfunctionstart_Callback(hObject,eventdata,handles)%hObjecthandletostart(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesandus

3、erdata(seeGUIDATA)subplot(3,3,2)filename,pathname=uigetfile(*.jpg;*.tif;*.png;*.gif,AllImageFiles;.*.*,AllFiles,open);%打開路徑下要檢索的圖像ifisequal(filename,pathname,0,0)returnelse%讀取圖片pic=fullfile(pathname,filename);globalonion;onion=imread(pic);imshow(onion);end%-Executesonbuttonpressinnext.functionnext_C

4、allback(hObject,eventdata,handles)t(seeGCBO)bedefinedinafutureversionofMATLABhhandlesanduserdata(seeGUIDATA)操作方法廿別點(diǎn)擊firstpic按鈕和secondpic按紀(jì)選擇圖片。點(diǎn)擊correlation按鈕進(jìn)行相關(guān)性比較。依次點(diǎn)擊overlay按鈕與transparent按鈕,講行重臺分析。L操作箱firstpicsecondpiccorrelationbjectoverhandletonex%hO%eventdatareserved-to%handle二2)選圖按鈕Sesubplot

5、(3,3,3);structureransparenisecondpicwitfilename,pathname=uigetfile(*.jpg;*.tif;*.png;*.gif,AllImageFiles*.*,AllFiles,open);%打開路徑下要檢索的圖像ifisequal(filename,pathname,0,0)returnelse%讀取圖片pict=fullfile(pathname,filename);globalpeppers;peppers=imread(pict);imshow(peppers);end%-Executesonbuttonpressinrestar

6、t.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=1

7、6347143151;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-tobedefine

8、dinafutureversionofMATLAB%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_p

9、eppers(:,:,1);%offsetfoundbycorrelationmax_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+

10、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.pn

11、gwasextractedfrompeppers.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-tobedef

12、inedinafutureversionofMATLAB%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

13、(:,:,1),sub_peppers(:,:,1);%offsetfoundbycorrelationmax_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

14、=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)

15、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%overlayimageswithtransparencyimsh

16、ow(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為了在