消除红眼matlab程序.zip_begantgn_matlab_去红眼_去红眼 matlab_图像去红颜

1 %本程序的功能是去除用户指定区域的红眼 clear all;clc;%内存清理 [file, path] = uigetfile('*.bmp', '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=1/6; %预设对rol区域内-pi/3<H<pi/3,S>0.2的像素点进行去红眼操作 h2=1/6; s=0.2; fprintf('预设对rol区域内-pi/3<H<pi/3,S>0.2的像素点进行去红眼操作，即令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下限值（-2*pi/3~0）：')/pi/2; while (h1>0)||(h1<-1/3) fprintf('错误命令！请重新输入~~\n'); h1=input('请输入H下限值（-2*pi/3~0）：')/pi/2; end; h2=input('请输入H上限值（0~2*pi/3）：')/pi/2; while (h2<0)||(h2>1/3) fprintf('错误命令！请重新输入~~\n'); h2=input('请输入H上限值（0~2*pi/3）：')/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 %用户通过鼠标绘制任意多边形，点左键绘制顶点，点右键绘制最后的一 %个顶点 [hsi,H,S,I]=rgb2hsi(filename); %rgb->hsi rol=find(B~=0); %获取rol内像素数 [L,M]=size(rol); for i=1:L %遍历rol内像素，判定-1/6<h<1/6且s>0.2为红眼，并令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; I(rol(i))=0; end; end; end; hsi=cat(3,H,S,I); %生成去红眼后的hsi图像 rgb=hsi2rgb(hsi); %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;