emd（2） 在MATLAB环境下实现二维emd分解的实现

在MATLAB环境下实现二维emd分解的实现,效果还不错

function imf=bemd（im） 
if nargin~=1  
 [filename pathname] = uigetfile（ {‘*.bmp‘ ‘Bitmap Files （*.bmp）‘;‘*.jpg‘ ‘Jpeg Files （*.jpg）‘;‘*.jpeg‘ ‘Jpeg Files （*.jpeg）‘} ‘Open a image file‘）; 
 fn=strcat（pathnamefilename）; 
 im=imread（fn）; 
end 
x=1:120; 
y=1:110; 
[xy]=meshgrid（yx）; 
im=40*sin（2*pi/40*x）+40*sin（2*pi/60*y）+30*sin（2*pi/14*x）+50*sin（2*pi/17*y）+190*sin（2*pi/0.2*x）+190*sin（2*pi/0.3*y）; 
%im=5*sin（2*pi/40*x+2*pi/60*y）+6*sin（2*pi/14*x+2*pi/17*y）+50*sin（2*pi/0.2*x+2*pi/0.3*y）; 
dim=size（im）; 
if size（dim2）>2 
    im=rgb2gray（im）; 
end 
figure;imshow（fn）; 
im=imresize（im[128108]‘bicubic‘）; 
i=0; 
im=double（im）;  
figure;surf（im）; 
title（‘original mesh‘）; 
ch=1.4; 
cw=1.4; 
rim=im; 
while 1
    if max（max（im））<7
        i=i+1; 
        imf（i::）=i

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1384  2013-01-04 20:50  emd（2）\bemd.m

     文件      66614  2013-01-04 18:05  emd（2）\lena.BMP

     文件       2469  2012-12-25 10:50  emd（2）\maxsurf.m

     文件       2493  2012-12-25 10:51  emd（2）\minsurf.m

     文件       1341  2012-12-25 10:51  emd（2）\wkl.m

     目录          0  2013-01-05 11:21  emd（2）

----------- ---------  ---------- -----  ----

                74301                    6