Question

求这段MATLAB代码注释！！！

clear;
A=imread('Water lilies','JPEG');
subplot(1,2,1);
subimage(A);
title('before');
A=im2double(A);%converts the intensity image A to double
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);%将A转换成uint8型
subplot(1,2,2);
subimage(A);
title('after');

Answer 1

这段程序用于增加图片的色调，要结合HSI空间的公式来看
clear;
A=imread('1','JPEG');%%读取名为Water lilies的图片
subplot(1,2,1);%%显示1行2列的图片组的第一幅
subimage(A);%%第一幅图片A
title('before');%%变换前的图片标题“title”
A=im2double(A);%转换数据格式，unchar到double
A=A/256;%%转换为0~1之间的小数
R=A(:,:,1);%%原始输入图片为RGB空间的，共3维分别表示R,G,B，这里就是取第一维R分量放入R矩阵中
G=A(:,:,2);
B=A(:,:,3);
%%%以下开始从RGB空间转化到HSI空间，具体的有相关公式参考
I=(R+G+B)/3;%%公式计算亮度矩阵I
S=1-min(R,min(G,B))./I;%%公式计算饱和度S
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));%%反余弦函数，H值为0~pi之间
t=find(B>G);
H(t)=2*pi-H(t);%转换成HSI
%%%%转换HSI结束
H=H+40/180*pi;%色调加40度，目的就在于这里,增加色调
%%%%%从HSI空间转化到RGB空间的反变换
t=find(H>2*pi);%%从这里开始就相当于上面操作的逆变换,只不过H已经增加了40
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;%%RGB三维赋值给图片矩阵A
A(:,:,2)=G;
A(:,:,3)=B;
A=A*256;
A=im2uint8(A);%将A转换成char型
subplot(1,2,2);
subimage(A);
title('after');

Answer 2

电脑:‘开始 所有程序 附件 命令提示符’
你操作一下