软件版本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程序)],运行不出来,就一点一点改成功了。
下面是效果对比图: