%%
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