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

1、例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);ti

2、tle('(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');

3、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=li

4、nspace(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(

5、9;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'); %便

6、于查看模拟的新闻报纸的白度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,

7、. '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('Fig03

8、04(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(

9、9;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)(J

10、LK 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).

11、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);ti

12、tle('(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_filtere

13、d);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)分

14、别平滑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 12

15、1;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图像的边缘。clc,clear,close all;R=imread('Fig0624(a)RGB2-red.tif');G=imread('

16、;Fig0624(b)RGB2-green.tif');B=imread('Fig0624(c)(RGB2-blue).tif');f=imread('Fig0624(d)(RGB2-fullcolor).tif');subplot 221; imshow(R); title('(a)RGB图像R分量');subplot 222; imshow(G); title('(b)RGB图像G分量');subplot 223; imshow(B); title('(c)RGB图像B分量');subplot 224;

17、 imshow(f); title('(d)RGB图像'); figure; VG,A,PPG=colorgrad(f);subplot 121;imshow(VG);title('(e)在RGB向量空间中直接计算的梯度');subplot 122;imshow(PPG);title('(f)分别计算每幅RGB分量图像的二维梯度并将结果相加得到的合成梯度');figure;f1=imread('Fig0604(a)(iris).tif');VG1,A1,PPG1=colorgrad(f1);sub=abs(VG1-PPG1);m1

18、=max(sub); %绝对差映射到0,1区间m2=max(m1);sub=sub*1/m2;subplot 141; imshow(f1); title('(a)RGB图像');subplot 142;imshow(VG1); title('(b)在RGB向量空间中直接计算的梯度');subplot 143;imshow(PPG1);title('(c)通过分别计算每幅RGB分量图像的二维梯度并将结果相加得到的合成梯度');subplot 144;imshow(sub); title('(d)图(b)(c)的绝对差，标定区间0,1

19、9;);运行结果：例5.11 RGB彩色图像分割。clc,clear,close all;f=imread('Fig0627(a)(jupitermoon_original).tif');mask=roipoly(f);red=immultiply(mask, f(:,:,1);green=immultiply(mask, f(:,:,2);blue=immultiply(mask, f(:,:,3);g=cat(3,red,green,blue);figure;subplot 121;imshow(f);title('(a)木卫Io表面区域的伪彩色图像');s

20、ubplot 122;imshow(g);title('(b)使用函数roipoly交互提取的感兴趣区域'); M,N,K=size(g); %提取ROI区域中的点的坐标I=reshape(g,M*N,3); %g中彩色像素重新排列为I的行idx=find(mask); %找出非黑彩色分量的标准差的倍数I=double(I(idx,1:3);C,m=covmatrix(I); %C的主对角线上包含RGB分量的方差d=diag(C);sd=sqrt(d)' E25=colorseg('euclidean',f,25,m);E50=colorseg('euclidean',f,50,m);E75=colorseg('euclidean',f,75,m);E100=colorseg('euclidean',f,100,m);figure;subplot 241;imshow(E25);title('(a)T为25时euclidean选项分割的结果