小波变换分割彩色图像

close all
clc
clear
img=imread（‘F:\图像\沉船滤波图像.jpg‘）;
img=rgb2gray（img）;
img=im2double（img）;
LUV=zeros（size（img））;
h = ones（55）/25;
MELH=zeros（size（img1）size（img2））;
MEHL=zeros（size（img1）size（img2））;
LUV = RGB2ULV（img）;
img_gray=rgb2gray（img）;
MELH_HL=wavlet_decomposition（img_gray） ;
Enl=0.00001;
err=0.0001;
Entropy_ns=0;
alpha=0.5;
X1=[];
%%%%%%%%%TL%%%%%%%%%%
while 1
[Entropy_nsalpha1 TLIL] = ns（LUV（::1）Enl）;
fprintf（‘Enl= %f  Entropy_ns= %f  error= %f\n‘ EnlEntropy_ns（Entropy_ns-Enl）/Enl）;
if （（Entropy_ns-Enl）/Enl）< err
    break
else
    Enl=Entropy_ns ;
end
end
disp（‘---------------------‘）;
X=zeros（size（img1）size（img2））;
TL_mean = imfilter（TLh‘same‘）;
X（ILX（IL>=alpha1）=TL_mean（IL>=alpha1）; %required
imshow（X）
X1（::1）=X;
%%%%%%%%%TU%%%%%%%%%%
Enl=0.00001;
while 1
[Entropy_nsalpha2 TUIU] = ns（LUV（::2）Enl）;
fprintf（‘Enl= %f  Entropy_ns= %f  error= %f\n‘ EnlEntropy_ns（Entropy_ns-Enl）/Enl）;
if （（Entropy_ns-Enl）/Enl）< err
    break
else
    Enl=Entropy_ns ;
end
end
disp（‘---------------------‘）;
X=zeros（size（img1）size（img2））;
TU_mean = imfilter（TUh‘same‘）;
X（IUX（IU>=alpha2）=TU_mean（IU>=alpha2）; %required
imshow（X）
X1（::2）=X;
%%%%%%%%%%TV%%%%%%%%%
Enl=0.00001;
while 1
[Entropy_nsalpha3 TVIV] = ns（LUV（::3）Enl）;
fprintf（‘Enl= %f  Entropy_ns= %f  error= %f\n‘ EnlEntropy_ns（Entropy_ns-Enl）/Enl）;
if （（Entropy_ns-Enl）/Enl）< err
    break
else
    Enl=Entropy_ns ;
end
end
disp（‘----------------