归一化互相关

实训题目用归一化互相关注册一幅图像指导老师盛利元专业班级电信0901姓名耿宝印学号 1404090114handletostart(seeGCBO)reservedhandletostart(seeGCBO)reserved-tobedefinedinafutureversionstructurewithhandlesanduserdata(seeRegisteringanImageUsingNormalizedCross-Correlation这次实训，是在参考matlab样题后进行的。刚接触这个题目时，自己对什么是互相关也不知道，经过查资料，相关说简单点，就是相似性。构思；由一及多。选择两张特定的图片，大致选定区域，然后用normxcorr2判定相关性，找出补偿点，对图片进行相似处理。最后进行图片扩展。实现尽可能多的图片的互相关的注册。具体步骤以样题代码为基础，在matlab自带guide平台进行图形设计，建立xiangguan.flg文件并相应产生xiangguan.m文件。在设计好大体的图形界面后，对xiangguan.m文件进行代码编写。t最终guide平台下的图形界面如下：为优化代码使用全局变量，并在相应的地方进行声明。t具体控件代码如下：1)选图按钮Firstpicfunctionstart_Callback(hObject,eventdata,handles)%hObject%eventdataofMATLAB%handlesGUIDATA)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_Callback(hObject,eventdata,handles)(seeGCBO)bedefinedinafutureversionofMATLABhandlesanduserdata(seeGUIDATA)[filename,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% Executesonbuttonpressinrestart.functionrestart_Callback(hObject,eventdata,handles)%hObjecthandletorestart(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)相关比较按钮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','保存退出','关闭运行窗口','取消')switchbuttoncase'保存退出'

saveexitcase,取消'quitcase'关闭运行窗口’closeend为了在同一