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

例5.1 函数说明。

clc,clear,close all;%函数demof=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转化到HSIf=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');运行结果：例5.7 伪彩色映射。

f4=imread('Fig0616(a)(Weld Original).tif.tif');g4=ice('image',f4);运行结果：例5.8 彩色平衡。

f5=imread('Fig0617(a)(JLK Magenta).tif.tif');g5=ice('image',f5,'space','CMY');运行结果：例5.9 基于直方图的映射。

f5=imread('Fig0618(a)(Caster Original).tif.tif');imshow(f5);g5=ice('image',f5,'space','hsi');运行结果：例5.10 彩色图像平滑。

clc,clear,close all;f=imread('Fig0619(a)(RGB_iris).tif');R=f(:,:,1); G=f(:,:,2); B=f(:,:,3);subplot 141;imshow(f);title('(a)原图');subplot 142;imshow(R);title('(a)R');subplot 143;imshow(G);title('(b)G');subplot 144;imshow(B);title('(c)B');h=rgb2hsi(f);H=h(:,:,1);S=h(:,:,2);I=h(:,:,3);figure;subplot 131;imshow(H);title('(a)H');subplot 132;imshow(S);title('(b)S');subplot 133;imshow(I);title('(c)I');w=fspecial('average',25); I_filtered=imfilter(I,w,'replicate');H_filtered=imfilter(H,w,'replicate'); S_filtered=imfilter(S,w,'replicate');h=cat(3,H,S,I_filtered);F=hsi2rgb(h);h1=cat(3,H_filtered,S_filtered,I_filtered);f2=hsi2rgb(h1);figure;r_filtered=imfilter(R,w,'replicate'); g_filtered=imfilter(G,w,'replicate');b_filtered=imfilter(B,w,'replicate');f1=cat(3,r_filtered,g_filtered,b_filtered);subplot 131;imshow(f1);title('(a)分别平滑RGB分量得到平滑后RGB图像');subplot 132;imshow(F);title('(b)HSI中亮度分量I滤波后');subplot 133;imshow(f2);title('(c)平滑所有三个HSI分量的结果');运行结果：例5.11 彩色图像锐化。

clc,clear,close all;f=imread('Fig0619(a)(RGB_iris).tif');w=fspecial('average',5);fb=imfilter(f,w,'replicate');subplot 121;imshow(fb);title('(a)模糊图像');lapmask=[1 1 1;1 -8 1;1 1 1];fb=tofloat(fb); fen=fb-imfilter(fb,lapmask,'replicate');subplot 122;imshow(fen);title('(b)使用拉普拉斯算子增强后的图像');运行结果：例5.11 使用函数colorgrad检测RGB图像的边缘。