归一化互相关

《归一化互相关》由会员分享，可在线阅读，更多相关《归一化互相关（8页珍藏版）》请在装配图网上搜索。

1、实训题目用归一化互相关注册一幅图像指导老师盛利元 专业班级电信 0901 姓名耿宝印学号 1404090114Registering an Image Using Normalized Cross-Correlation这次实训，是在参考 matlab 样题后进行的。刚接触这个题目时，自己对什么是互相关 也不知道，经过查资料，相关说简单点，就是相似性。构思；由一及多。选择两张特定的图片，大致选定区域，然后用 normxcorr2 判定相关 性，找出补偿点，对图片进行相似处理。最后进行图片扩展。实现尽可能多的图片的互相关 的注册。具体步骤以样题代码为基础，在 matlab 自带 guide 平台

2、进行图形设计，建立 xiangguan.flg 文件并 相应产生 xiangguan.m 文件。在设计好大体的图形界面后，对 xiangguan.m 文件进行代码编 写。最终 guide 平台下的图形界面如下：和 Hiangguan口皿为优化代码使用全局变量，并在相应的地方进行声明。具体控件代码如下：1)选图按钮First pic functionstart_Callback(hObject, eventdata, handles)% hObject handle to start (see GCBO)% eventdata reserved - to be defined in a futu

3、re version of MATLAB% handles structure with handles and user data (see GUIDATA)subplot(3,3,2)filename,pathname=uigetfile( *.jpg;*.tif;*.png;*.gif ,All Image Files ;.*.* ,All Files ,open);%打开路径下要检索的图像ifisequal(filename,pathname,0,0)returnelse%读取图片pic = fullfile(pathname,filename);global onion;onion

4、= imread(pic);imshow(onion);end% - Executes on button press in next.function next_Callback(hObject, eventdata, handles)t (see GCBO)be defined in a future version of MATLABh handles and user data (see GUIDATA)操作方法廿别点击first pic 按钮和second pic 按紀选择图片。点击 correlation 按 钮进行相关性比较。 依次点击overlay按 钮与 transparen

5、t 按 钮,讲行重台分析。L操作箱first picsecond piccorrelatio nbjec t ove rhandle tonex% hO% eventdata reserved - to% handle二2)选图按钮Sesubplot(3,3,3);structureransparenisecon d picwitfilename,pathname=uigetfile( *.jpg;*.tif;*.png;*.gif ,All Image Files *.* ,All Files ,open);%打开路径下要检索的图像ifisequal(filename,pathname,0,

6、0)returnelse%读取图片pict = fullfile(pathname,filename);global peppers;peppers = imread(pict);imshow(peppers);end% - Executes on button press in restart. function restart_Callback(hObject, eventdata, handles)% hObject handle to restart (see GCBO)% eventdata reserved - to be defined in a future version o

7、f MATLAB% handles structure with handles and user data (see GUIDATA)3) 相关比较按钮Correlationsubplot(3,3,5)global onion; global peppers;% non-interactively rect_onion = 111 33 65 58; rect_peppers = 163 47 143 151; sub_onion = imcrop(onion,rect_onion); sub_peppers = imcrop(peppers,rect_peppers);c = normxc

8、orr2(sub_onion(:,:,1),sub_peppers(:,:,1); surf(c), shading flat% - Executes on button press in overlay.function overlay_Callback(hObject, eventdata, handles)% hObject handle to overlay (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and us

9、er data (see GUIDATA)4) 按钮 Overlay subplot(3,3,6) global onion; global peppers;% non-interactively rect_onion = 111 33 65 58; rect_peppers = 163 47 143 151; sub_onion = imcrop(onion,rect_onion);sub_peppers = imcrop(peppers,rect_peppers);c = normxcorr2(sub_onion(:,:,1),sub_peppers(:,:,1);% offset fou

10、nd by correlationmax_c, imax = max(abs(c(:);ypeak, xpeak = ind2sub(size(c),imax(1); corr_offset = (xpeak-size(sub_onion,2)(ypeak-size(sub_onion,1);% relative offset of position of subimages rect_offset = (rect_peppers(1)-rect_onion(1)(rect_peppers(2)-rect_onion(2);% total offsetoffset = corr_offset

11、+ 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);% extract region from peppers and compare to onion extracted_onion = peppers(ybegin:yend,xbegin:xend,:); ifisequal(onion,ex

12、tracted_onion)disp( onion.png was extracted from peppers.png ) endrecovered_onion = uint8(zeros(size(peppers); recovered_onion(ybegin:yend,xbegin:xend,:) = onion; imshow(recovered_onion)% - Executes on button press in transparent. function transparent_Callback(hObject, eventdata, handles)% hObject h

13、andle to transparent (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)5)透明处理按钮Transparen t subplot(3,3,8)global onion; global peppers;% non-interactively rect_onion = 111 33 65 58;rect_peppers = 163 47 143 151; s

14、ub_onion = imcrop(onion,rect_onion);sub_peppers = imcrop(peppers,rect_peppers);c = normxcorr2(sub_onion(:,:,1),sub_peppers(:,:,1);% offset found by correlationmax_c, imax = max(abs(c(:);ypeak, xpeak = ind2sub(size(c),imax(1); corr_offset = (xpeak-size(sub_onion,2)(ypeak-size(sub_onion,1);% relative

15、offset of position of subimages rect_offset = (rect_peppers(1)-rect_onion(1)(rect_peppers(2)-rect_onion(2);% total offsetoffset = 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(yof

16、fset+size(onion,1);% extract region from peppers and compare to onion extracted_onion = peppers(ybegin:yend,xbegin:xend,:); ifisequal(onion,extracted_onion)disp( onion.png was extracted from peppers.png ) endrecovered_onion = uint8(zeros(size(peppers); recovered_onion(ybegin:yend,xbegin:xend,:) = on

17、ion;m,n,p = size(peppers);mask = ones(m,n);i = find(recovered_onion(:,:,1)=0);mask(i) = .2; % try experimenting with different levels of% transparency% overlay images with transparency imshow(peppers(:,:,1) % show only red plane of peppers holdonh = imshow(recovered_onion); % overlay recovered_onion

18、set(h, AlphaData ,mask)菜单 quit% - Executes on button press in quit.function quit_Callback(hObject, eventdata, handles)% hObject handle to quit (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)butt on=ques tdlg(rea

19、dy to qui t? , exi t dialog,保存退出,关闭运行窗口 ,取消) switch button case 保存退出save exitcase取消quitcase 关闭运行窗口closeend为了在同一个图形界面显示出完整的图像，用subpl。t ()函数对图形界面进行分区。用 uigetfile ()函数指定图像的路径，以使程序得到更多的应用。程序操作：点选firs t pic和second pic按钮会出现如下对话框一次点击选项按钮后一个样例如下图所示1200200-1J0400点击透明处理按钮后则会得到最终影像在 guide 平台上验证无误后利用matlab自带的编译器进行最终的程序输出 代码如下mex-setup用 mcc 函数编译mcc -m xiangguan.m 编译成功后得到运行 xiangguan.exe 文件结果如下0 0