颜色模型的应用——利用HSI模型特性消除红眼

软件版本Matlab R2015a

在网上找了张汪的照片，jpg格式，下面直接上干活：

clear all;clc;                             %内存清理
[file, path] = uigetfile('*.jpg', 'Image');%获取图像名及存储位置
if isequal(file,0)                         %未选择图像，退出
    disp('取消选择');
else
    file=[path,file];                      %获取图像路径
    [filename,map]=imread(file);           %读取图像
    if  ndims(filename)~=3                 %是否为灰度图像,是退出并显示图像
        disp('非彩色图像');
        imshow(filename);
    else                                   %去红眼操作
       goagain=1;                          %操作控制变量，1进行操作
       change=2;                           %参数更改控制变量，1更改预设参数
       h1=-pi/4;                           %预设对rol区域内-pi/4<H<pi/4,S>0.3的像素点进行去红眼操作
       h2=pi/4;
       s=0.3;
       fprintf('预设对rol区域内-pi/4<H<pi/4,S>0.3的像素点进行去红眼操作，即令S=0\n')
       while goagain==1                           
           fprintf('是否修改H、S预设值：1、是；2、否\n')
           change=input('请选择：');       %是否更改预设参数
           while change~=1&change~=2         
               fprintf('错误命令！请重新输入~~\n');
               change=input('请选择：');
           end
           if change==1                    %预设参数更改
               h1=input('请输入H下限值：')/pi/2;
               while (h1>0)||(h1<-1/3)
                   fprintf('错误命令！请重新输入~~\n');
                   h1=input('请输入H下限值：')/pi/2;
               end;
               h2=input('请输入H上限值：')/pi/2;
               while (h2<0)||(h2>1/3)
                   fprintf('错误命令！请重新输入~~\n');
                   h2=input('请输入H上限值：')/pi/2
               end;
               s=input('请输入S值(0~1)：');
               while (s>1)||(s<0)
                   fprintf('错误命令！请重新输入~~\n');
                   s=input('请输入S值(0~1)：');
               end;
               change=2;
           end;
           B=roipoly(filename);            %绘制要去红眼的rol区域，返回与filename大小相同的逻辑矩阵B，rol内像素值为1，外为0
                                           %用户通过鼠标绘制任意多边形，点左键绘制顶点，点右键绘制最后的一
                                           %个顶点

           J=rgb2hsi(filename);            %rgb->hsi
           H=J(:,:,1);
           S=J(:,:,2);
           I=J(:,:,3);

           rol=find(B~=0);                 %获取rol内像素数；find 函数——找出矩阵X中的所有非零元素
           [L,M]=size(rol);                %size函数 当有两个输出参数时，函数将矩阵的行数返回到第一个输出变量L，将矩阵的列数返回到第二个输出变量M
           for i=1:L                       %遍历rol内像素，判定-pi/4<h<pi/4且s>0.3为红眼，并令s=0以去除红眼
               if (((H(rol(i))>0)&(H(rol(i))<h2))|((H(rol(i))>(1-abs(h1)))&(H(rol(i))<1)))
                   if S(rol(i))>s
                       S(rol(i))=0;
                   end;
               end;
           end;
           J=cat(3,H,S,I);                 %生成去红眼后的hsi图像
           rgb=hsi2rgb(J);                 %hsi->rgb
           imshow(rgb);                    %显示去红眼后的rgb图像
           fprintf('是否继续操作：1、是；2、否\n');              %提示是否继续操作
           goagain=input('请选择：');
           while goagain~=1&goagain~=2
              fprintf('错误命令！请重新输入~~\n');
              goagain=input('请选择：');
           end
           if goagain==1                   
               filename=rgb;
               clf;
           end;
       end;
    end
end;

程序是从网上搜的，关键词 [红眼消除-多媒体作业(附matlab程序)]，运行不出来，就一点一点改成功了。

下面是效果对比图：