数字图像处理代码Ch5《彩色图像处理》

例5.1 函数说明。

clc,clear,close all;

%函数demo

f=imread('Fig0604(a)(iris).tif');

subplot 131;imshow(f);title('(a)RGB图像');

[X1,map1]=rgb2ind(f,8,'nodither');

subplot 132;imshow(X1,map1);title('(b)未经抖动处理的颜色数减少到8的图像');

[X2,map2]=rgb2ind(f,8,'dither');

subplot 133;imshow(X2,map2);title('(c)经抖动处理的颜色数减少到8的图像');

g=rgb2gray(f); g1=dither(g);

figure;subplot 121;imshow(g);

title('(d)使用函数rgb2gray得到的图(a)的灰度图像');

subplot 122;imshow(g1);

title('(e)抖动处理后的灰度图像（二值图像）');

运行结果：

例5.2 RGB转化为HSI。

clc,clear,close all;

%从RGB转化到HSI

f=imread('Fig0602(b)(RGB_color_cube).tif'); subplot 221;imshow(f);

title('(a)RGB图像');

f1=rgb2hsi(f);

H=f1(:,:,1);

S=f1(:,:,2);

I=f1(:,:,3);

subplot 222;imshow(H);title('(b)色调图像'); subplot 223;imshow(S);title('(c)饱和度图像'); subplot 224;imshow(I);title('(d)亮度图像');

运行结果：

例5.3 基于L*a*b彩色空间创建一个感觉上一致的彩色空间。

clc,clear,close all;

L=linspace(40,80,1024); %创建40到80之间等分1024个值的斜坡

radius=70; %明显不同的相邻颜色，通过改变彩色的色调满足，

theta=linspace(0,pi,1024); %色调对应于a*b*平面中彩色坐标的极角

a=radius*cos(theta); b=radius*sin(theta);

L=repmat(L,100,1);a=repmat(a,100,1);b=repmat(b,100,1);

%创建该L*a*b*彩色标尺的一幅100×1024×3 的图像

lab_scale=cat(3,L,a,b);

cform=makecform('lab2srgb');

rgb_scale=applycform(lab_scale,cform);

imshow(rgb_scale);title('基于L*a*b彩色空间的感觉上一致的标尺');

运行结果：

例5.4 ICC彩色剖面的软件证明。

clc,clear,close all;

f=imread('Fig0604(a)(iris).tif');

fp=padarray(f,[40 40],255,'both');%围绕图像加一较粗的白色边框和一个较细灰色边框fp=padarray(fp,[4 4],230,'both'); %便于查看模拟的新闻报纸的白度

subplot 121;imshow(fp);title('(a)带有白色边框的原始图像');

p_srgb=iccread('sRGB.icm'); %读入两个剖面，使它们把虹膜图像从sRGB转为新闻纸颜色p_snap=iccread('SNAP2007.icc');

cform1=makecform('icc',p_srgb,p_snap); fp_newsprint=applycform(fp,cform1); cform2=makecform('icc',p_snap,p_srgb,...

'SourceRenderingIntent','AbsoluteColorimetric',..

'DestRenderingIntent','AbsoluteColorimetric'); fp_proof=applycform(fp_newsprint,cform2);

subplot 122; imshow(fp_proof);title('(b)图像打印到新闻纸上后的模拟外观');

运行结果：

例5.5 单色负片和彩色分量的反映射。

clc,clear,close all;

f=imread('Fig0304(a)(breast_digital_Xray).tif');

g=ice('image',f,'space','rgb');

f1=imread('Fig0614(a)(Chalk Original).tif.tif');

g1=('image',f,'space','rgb'); %得到图像的负片

运行结果：

例5.6 单色和彩色对比度增强。

f2=imread('Fig0615(a)(Aerial Original).tif.tif');

g2=ice('image',f2,'space','rgb');

f3=imread('Fig0615(d)(Iris Original).tif.tif');

g3=ice('image',f3,'space','rgb');

运行结果：