bemd 二维经验模态分解

二维经验模态分解，分解图像等二维信号，可用于算法验证与学习 the process of sifting. It is dependent on the function SIFT.)

%%
function [ imf_matrix ] = bemd（ input_image ）

% BEMD This program calculates the Bidimensional EMD of a 2-d signal using
% the process of sifting. It is dependent on the function SIFT.

tic

% Make a ‘double‘ copy of the image

input_image_d = cast（input_image‘double‘）;

% Assigning the initial values of ‘h‘ （data for sifting） and ‘k‘ （imf index）

h_func = input_image_d;
k=1;

% The process of Sifting

while（k<4）
    [imf_temp residue_temp] = sift（h_func）; 
    imf_matrix（::k） = imf_temp; %#ok
    k = k+1;
    h_func = residue_temp;
end

% Assigning the final residue to the last IMF index

imf_matrix（::k） = residue_temp;

% End of BEMD Computation

toc

end

%%

function [ h_imf residue ] = sift（ input_image ）

% This function sifts for a single IMF of the given 2D signal i

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件       45810  2010-09-19 14:29  bemd.m
     文件        1364  2010-09-19 14:29  license.txt