z2.asv

这是一般都会有的作业答案 请认真阅读您的文件包然后写出其具体功能(至少要20个字)。尽量不要让站长把时间都花费在为您修正说明上。压缩包解压时不能有密码。系统会自动删除debug和release目录

lear;
A=imread('Water lilies','JPEG');
subplot(1,2,1);
subimage(A);
title('变化前的图');
A=im2double(A);
A=A/256;
R=A(:,:,1);
G=A(:,:,2);
B=A(:,:,3);
I=(R+G+B)/3;
S=1-min(R,min(G,B))./I;
H=zeros(size(S));%为H分配空间
t=find(S==0.0);%当S为时，H=0.0
H(t)=0.0;
t=find(S~=0.0);
H(t)=((R(t)-G(t))+(R(t)-B(t)))/2.0;
H(t)=H(t)./sqrt((R(t)-G(t)).^2+(R(t)-B(t)).*(G(t)-B(t)));
H(t)=acos(H(t));
t=find(B>G);
H(t)=2*pi-H(t);%转换成HSI 
H=H+40/180*pi;%色调加40度
t=find(H>2*pi);
H(t)=H(t)-2*pi;
find(0<=H&H<=2*pi/3);
B(t)=I(t).*(1-S(t));
R(t)=I(t).*(1+S(t).*cos(H(t))./cos(pi/3-H(t)));
G(t)=3*I(t)-(B(t)+R(t));
t=find(2*pi/3<H&H<=4*pi/3);
R(t)=I(t).*(1-S(t));
G(t)=I(t).*(1+S(t).*cos(H(t)-2*pi/3)./cos(pi-H(t)));
B(t)=3*I(t)-(R(t)+G(t));
t=find(4*pi/3<H&H<=2*pi);
G(t)=I(t).*(1-S(t));
B(t)=I(t).*(1+S(t).*cos(H(t)-4*pi/3)./cos(5*pi/3-H(t)));
R(t)=3*I(t)-(G(t)+B(t));%转换成RGB
A(:,:,1)=R;
A(:,:,2)=G;
A(:,:,3)=B;
A=A*256;
A=im2uint8(A);
subplot(1,2,2);
subimage(A);
title('转换后的图');