图像处理作业第五次图像变暗处理

1、数字图像处理第五次作业汇报人：杨鹏斌1 1）设计内容： 将图像caomei.bmp分解成HIS模型，显示H、S、I三变量的灰度图像，分别改变三分量，显示改变后的图像。 2）设计要求：用matlab或VC+等语言编程，编程时，设计中除了读写、显示、数据类型转换等函数外以外，其它主要功能函数需自己编写。 ）制作PPT进行讲解，限时8分钟。一2二变换原理HSI 色彩空间是从人的视觉系统出发，用色调（Hue）、色饱 和度（Saturation或Chroma）和亮度 （Intensity或Brightness）来描述色彩。HSI 色彩空间可以用一个圆锥空间模型来描述。用这种 描述HIS 色彩空间的圆锥模

2、型相当复杂，但确能把色调、亮度和色饱和度的变化情形表现得很清楚。 通常把色调和饱和度通称为色度，用来表示 颜色的类别与深浅程度。由于人的视觉对 亮度的敏感 程度远强于对 颜色浓淡的敏感程度，为了便于色彩处理和识别，人的视觉系统经常采用HSI 色彩空间， 它比RGB色彩空间更符合人的视觉特性。31.读入图片并进行显示；2.进行R,G,B三通道显示；3.运用变换公式对图像进行颜色模式的转换；4.显示H,I,S三变量的灰度图像；5.分别对相应变量进行滤波处理；6.显示改变变量后的图像。三设计流程4原始图像四5改变RGB三通道效果图五6六H、S、I三变量的灰度图像7对比图：8最终对比图：9matlab

3、程序代码%数字图像处理:RGB与HSI模型转换Matlab程序% im1=imread(caomei.bmp); im3=im1; im1=im2double(im1); m,n,q=size(im1); % 获取图像的 RGB 3个通道 R=im1(:,:,1); G=im1(:,:,2); B=im1(:,:,3); % 改变通道来观察图像 figure,subplot(1,3,1),imshow(cat(3,R*0.5,G,B),title(R通道减半); set(gcf,outerposition,get(0,screensize); set(gcf,NumberTitle,off,N

4、ame,改变通道来观察图像); subplot(1,3,2),imshow(cat(3,R,G*2,B),title(G通道翻倍); subplot(1,3,3),imshow(cat(3,R,G,zeros(m,n),title(B通道为0); % 观察HSI通道图像 % 因系统没有rgb2hsi函数,只有rgb2hsv函数,所以自己编写程序实现 H=zeros(m,n); S=H; for i1=1:m for i2=1:n numerator=0.5*(R(i1,i2)-G(i1,i2)+R(i1,i2)-B(i1,i2); denominator=sqrt(R(i1,i2)-G(i1,

5、i2)2+(R(i1,i2)-B(i1,i2)*(G(i1,i2)-B(i1,i2); theta=acos(numerator/denominator)*180/pi; if(B(i1,i2)=G(i1,i2) H(i1,i2)=theta; else H(i1,i2)=360-theta; end min1=min(R(i1,i2),G(i1,i2); min1=min(B(i1,i2),min1); S(i1,i2)=1-3/(R(i1,i2)+G(i1,i2)+B(i1,i2)*min1; endendI=(R+G+B)/3; figure,subplot(1,3,1),imshow(

6、H,),title(色度图H); set(gcf,outerposition,get(0,screensize); set(gcf,NumberTitle,off,Name,HSI通道图像);subplot(1,3,2),imshow(S),title(饱和度图S); subplot(1,3,3),imshow(I),title(强度图I); 10% 从HSI转换回RGB模型并观察图像,系统也不自带hsi2rgb,仅自带hsv2rgb im2=zeros(m,n,3); for i1=1:m for i2=1:n value1=I(i1,i2)*(1+S(i1,i2)*cos(H(i1,i2)

7、/cos(pi/3-H(i1,i2)/180*pi); value2=I(i1,i2)*(1-S(i1,i2); switch (floor(H(i1,i2)/120) case 0 im2(i1,i2,1)=value1; im2(i1,i2,2)=3*I(i1,i2)-R(i1,i2)-B(i1,i2); im2(i1,i2,3)=value2; case 1 H(i1,i2)=H(i1,i2)-120; im2(i1,i2,1)=value2; im2(i1,i2,2)=value1; im2(i1,i2,3)=3*I(i1,i2)-R(i1,i2)-G(i1,i2); case 2 H(i1,i2)=H(i1,i2)-240; im2(i1,i2,1)=3*I(i1,i2)-G(i1,i2)-B(i1,i2); im2(i1,i2,2)=value2; im2(i1,i2,3)=value1; end endend% 模型转换完毕,显示图像figure,subplot(1,2,1),imshow(im1),title(原图); set(gcf,outerposition,g