R(ind)=I(ind).*(1.0+S(ind).*cos(H(ind))./cos(pi/3.0-H(ind))); G(ind)=1.0-(R(ind)+B(ind)); ind=find((H>2*pi/3)&(H<4*pi/3)); H(ind)=H(ind)-pi*2/3; R(ind)=I(ind).*(1.0-S(ind));
G(ind)=I(ind).*(1.0+S(ind).*cos(H(ind))./cos(pi/3.0-H(ind))); B(ind)=1.0-(R(ind)+G(ind)); ind=find((H>=4*pi/3)&(H<2*pi)); H(ind)=H(ind)-pi*4/3; G(ind)=I(ind).*(1.0-S(ind));
B(ind)=I(ind).*(1.0+S(ind).*cos(H(ind))./cos(pi/3.0-H(ind))); R(ind)=1.0-(G(ind)+B(ind)); RGB=cat(3,R,G,B); figure();
imshow(rgb1);%RGB原图像 figure();
imshow(RGB);%对I分量均衡化后的彩色图像转换成拥有显示的RGB图像
实验结果：
H分量 I分量
I分量直方图均衡化后结果 I分量直方图
S分量 I分量直方图均衡化后直方图
RGB原图 对I分量均衡化后的色彩图像转化成拥有显示的RGB图像
实验结果分析：
其强度分量I的范围值归一化后为[0,1]。正如在处理前强度分量I直方图中看到的，强度分布还是比较宽的。只对强度I均衡化处理而不改变图像的色调H和饱和度值S，将其结果转换到RGB空间，从显示的图像上可以看出，它的确影响啦整体图像的彩色感官。 6.RGB彩色空间分割：
rgb=imread('006.jpg'); rgb1=im2double(rgb); r=rgb1(:,:,1); g=rgb1(:,:,2); b=rgb1(:,:,3); figure(); imshow(rgb1); figure(); imshow(r); figure(); imshow(g); figure(); imshow(b);
r1=r(129:256,86:170); r1_u=mean(mean(r1(:))); [m,n]=size(r1); sd1=0.0; for i=1:m for j=1:n
sd1=sd1+(r1(i,j)-r1_u)*(r1(i,j)-r1_u); end end
r1_d=sqrt(sd1/(m*n));
r2=zeros(size(rgb1,1),size(rgb1,2));
ind=find((r>r1_u-1.25*r1_d)&(r<r1_u)+1.25*r1_d); r2(ind)=1; figure(); imshow(r2);
实验分析：
参考文献：
(1)数字图像处理及MATLAB实现，电子工业出版社，杨杰。
(2)数字图像处理及MATLAB实现—学习与实验指导，电子工业出版
社，杨杰。
(3)刘维一，于德月，王肇圻等.用迭代法消除数字图像放大后 的模糊[J].光电子激光，2002，13（4）：398~400. (4)陶洪.数字共焦显微技术及其图像复原算法研究四川大学 硕士学位论文[D].CNKI::CDMD:10610.2.2003.6632. (5)罗军辉，冯平等.MATLAB7.0在图像处理中的应用[M].北 京：机械工业出版社，2006：202~204，257~278.