imextrem.m

这个是EMD分解算法,相对比较简单,十分适合初学者,我开始也是从个程序入手的,觉得十分不错,与大家共享.

function [Imax,Imin] = Imextrem(Image,conn)
% Image 是要检测极大值的图像
%Imextrem 函数用来计算灰度图像的区域极值.
%其中Imax、Imin分别代表图像Image的区域极大值图像和区域极小值图像；Image为输入的原图像；conn为区域邻域大小，通常取conn=4或conn=8
%I=double(Image);
%---计算Image的区域极大值---%
X=Image;
if conn==4
    Ire=imregionalmax(X,4); % 在检测到极大值点的位置处用1来表示，其余点用0来表示
    %将检测到的极大值位置映射回原图像，在极大值点处的像素值用原图像像素值表示，其余点用0表示
    Imax=Ire.*X;
end
if conn==8
    Ire=imregionalmax(X,8); % 在检测到极大值点的位置处用1来表示，其余点用0来表示
    %将检测到的极大值位置映射回原图像，在极大值点处的像素值用原图像像素值表示，其余点用0表示
    Imax=Ire.*X;
end

%---计算Image的区域极小值---%
XX=Image;
if conn==4
    Iree=imregionalmin(XX,4); % 在检测到极小值点的位置处用1来表示，其余点用0来表示
    %将检测到的极大值位置映射回原图像，在极小值点处的像素值用原图像像素值表示，其余点用0表示
    Imin=Iree.*XX;
end
if conn==8
    Iree=imregionalmin(XX,8); % 在检测到极小值点的位置处用1来表示，其余点用0来表示
    %将检测到的极大值位置映射回原图像，在极小值点处的像素值用原图像像素值表示，其余点用0表示
    Imin=Iree.*XX;
end

%figure;imshow(Image,[])
%figure;imshow(Imax,[])
%figure;imshow(Imin,[])