前言

其实数字图像处理并不是我的研究方向，仅仅是因为在开题中简单提到了基于图像的设备安全监控，就这样被导师安排去学习图像处理的相关知识。

Orz……好吧，假期回家还不到两周就这样匆匆忙忙的提前赶回了学校，和实验室三维重建方向的的筒子们一起去上美国罗文大学唐瑛（Ying (Gina) Tang）教授开设的数字图像处理入门课程。这周从周一到周三一共上了三节大课，算是简单了解了图像处理的相关内容吧，课程结束后估计也就是浅尝辄止。

在唐教授的安排下阅读了Oliva, Torrallaba, 和Schyns所发表的论文《Hybrid images》，今天主要写写如何用matlab将论文所提到的图像混合进行代码的实现（估计这是图像处理小白级的内容了吧，大神请无视哈~）

实现的整体思想比较简单，就是对两幅图像分别进行高通和低通滤波，从而将两幅图像融合到一起，同时通过修改截止频率还可以改善融合效果，这样就可以实现在近处看是其中一副图像的效果，而在远处看又是另一幅图像的效果，如下所示：

1. 待混合图像1
   
2. 待混合图像2
   
3. 混合完成图像
   

matlab实现hybrid images

hybrid_images函数的编写

对两幅图像分别进行高通和低通滤波，其中涉及到傅里叶变换和反变换，图像混合函数所在文件hybrid_images.m如下：

function re=hybrid_images(ImageFile1,ImageFile2,radius_input)
    
    I1=ImageFile1;
    I2=ImageFile2;
    ImgFileOut='C:\Users\dell\Desktop\hybird_images/out.jpg';
    
    radius = radius_input;      
    I1_ = fftshift(fft2(double(I1)));  %进行傅里叶变换，并将零点移动到中心
    I2_ = fftshift(fft2(double(I2)));
    [m, n, z]=size(I1);
    h = fspecial('gaussian', [m,n], radius);
    h = h./max(max(h));
    for colorI = 1:3
       J_(:,:,colorI) = I1_(:,:,colorI).*(1-h) + I2_(:,:,colorI).*h;
    end
    J = uint8(real(ifft2(ifftshift(J_))));          %反变换
    imwrite(J,ImgFileOut);    %图片保存到本地
    re=J;
end

GUI设计

GUI关键代码

% --- Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
global im1;
[filename,pathname]=uigetfile({'*.*';'*.bmp';'*.jpg';'*.tif';'*.jpg'},'选择图像');
if isequal(filename,0)||isequal(pathname,0)
  errordlg('您还没有选取图片！！','温馨提示');%如果没有输入，则创建错误对话框 
  return;
else
    image=[pathname,filename];%合成路径+文件名
    im1=imread(image);%读取图像
    set(handles.axes1,'HandleVisibility','ON');%打开坐标，方便操作
    axes(handles.axes1);%%使用图像，操作在坐标1
    imshow(im1);%在坐标axes1显示原图像 
end


% --- Executes on button press in pushbutton2.
function pushbutton2_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton2 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
global im2;
[filename,pathname]=uigetfile({'*.*';'*.bmp';'*.jpg';'*.tif';'*.jpg'},'选择图像');
if isequal(filename,0)||isequal(pathname,0)
  errordlg('您还没有选取图片！！','温馨提示');%如果没有输入，则创建错误对话框 
  return;
else
    image=[pathname,filename];%合成路径+文件名
    im2=imread(image);%读取图像
    set(handles.axes2,'HandleVisibility','ON');%打开坐标，方便操作
    axes(handles.axes2);%%使用图像，操作在坐标1
    imshow(im2);%在坐标axes1显示原图像 
end


% --- Executes on slider movement.
function slider1_Callback(hObject, eventdata, handles)
% hObject    handle to slider1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'Value') returns position of slider
%        get(hObject,'Min') and get(hObject,'Max') to determine range of slider
slide_value=get(hObject,'Value');%获取滑块当前值
global im1;
global im2;
radius=slide_value;
set(handles.edit1,'string',num2str(slide_value));
x=hybrid_images(im1,im2,radius);
set(handles.axes5,'HandleVisibility','ON');%打开坐标，方便操作
axes(handles.axes5)
imshow(x);

% --- Executes during object creation, after setting all properties.
function slider1_CreateFcn(hObject, eventdata, handles)
% hObject    handle to slider1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: slider controls usually have a light gray background.
if isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor',[.9 .9 .9]);
end

运行程序

paper链接

http://cvcl.mit.edu/hybrid/OlivaTorralb_Hybrid_Siggraph06.pdf