數(shù)字圖像處理代碼Ch5《彩色圖像處理》

./例5.1函數(shù)說明。clc,clear,closeall;%函數(shù)demof=imread<'Fig0604<a><iris>.tif'>;subplot131;imshow<f>;title<'<a>RGB圖像'>;[X1,map1]=rgb2ind<f,8,'nodither'>;subplot132;imshow<X1,map1>;title<'<b>未經(jīng)抖動(dòng)處理的顏色數(shù)減少到8的圖像'>;[X2,map2]=rgb2ind<f,8,'dither'>;subplot133;imshow<X2,map2>;title<'<c>經(jīng)抖動(dòng)處理的顏色數(shù)減少到8的圖像'>;g=rgb2gray<f>;g1=dither<g>;figure;subplot121;imshow<g>;title<'<d>使用函數(shù)rgb2gray得到的圖<a>的灰度圖像'>;subplot122;imshow<g1>;title<'<e>抖動(dòng)處理后的灰度圖像〔二值圖像'>;運(yùn)行結(jié)果：例5.2RGB轉(zhuǎn)化為HSI。clc,clear,closeall;%從RGB轉(zhuǎn)化到HSIf=imread<'Fig0602<b><RGB_color_cube>.tif'>;subplot221;imshow<f>;title<'<a>RGB圖像'>;f1=rgb2hsi<f>;H=f1<:,:,1>;S=f1<:,:,2>;I=f1<:,:,3>;subplot222;imshow<H>;title<'<b>色調(diào)圖像'>;subplot223;imshow<S>;title<'<c>飽和度圖像'>;subplot224;imshow<I>;title<'<d>亮度圖像'>;運(yùn)行結(jié)果：例5.3基于L*a*b彩色空間創(chuàng)建一個(gè)感覺上一致的彩色空間。clc,clear,closeall;L=linspace<40,80,1024>;%創(chuàng)建40到80之間等分1024個(gè)值的斜坡radius=70;%明顯不同的相鄰顏色,通過改變彩色的色調(diào)滿足,theta=linspace<0,pi,1024>;%色調(diào)對(duì)應(yīng)于a*b*平面中彩色坐標(biāo)的極角a=radius*cos<theta>;b=radius*sin<theta>;L=repmat<L,100,1>;a=repmat<a,100,1>;b=repmat<b,100,1>;%創(chuàng)建該L*a*b*彩色標(biāo)尺的一幅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彩色空間的感覺上一致的標(biāo)尺'>;運(yùn)行結(jié)果：例5.4ICC彩色剖面的軟件證明。clc,clear,closeall;f=imread<'Fig0604<a><iris>.tif'>;fp=padarray<f,[4040],255,'both'>;%圍繞圖像加一較粗的白色邊框和一個(gè)較細(xì)灰色邊框fp=padarray<fp,[44],230,'both'>;%便于查看模擬的新聞報(bào)紙的白度subplot121;imshow<fp>;title<'<a>帶有白色邊框的原始圖像'>;p_srgb=iccread<'sRGB.icm'>;%讀入兩個(gè)剖面,使它們把虹膜圖像從sRGB轉(zhuǎn)為新聞紙顏色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>;subplot122;imshow<fp_proof>;title<'<b>圖像打印到新聞紙上后的模擬外觀'>;運(yùn)行結(jié)果：例5.5單色負(fù)片和彩色分量的反映射。clc,clear,closeall;f=imread<'Fig0304<a><breast_digital_Xray>.tif'>;g=ice<'image',f,'space','rgb'>;f1=imread<'Fig0614<a><ChalkOriginal>.tif.tif'>;g1=<'image',f,'space','rgb'>;%得到圖像的負(fù)片運(yùn)行結(jié)果：例5.6單色和彩色對(duì)比度增強(qiáng)。f2=imread<'Fig0615<a><AerialOriginal>.tif.tif'>;g2=ice<'image',f2,'space','rgb'>;f3=imread<'Fig0615<d><IrisOriginal>.tif.tif'>;g3=ice<'image',f3,'space','rgb'>;運(yùn)行結(jié)果：例5.7偽彩色映射。f4=imread<'Fig0616<a><WeldOriginal>.tif.tif'>;g4=ice<'image',f4>;運(yùn)行結(jié)果：例5.8彩色平衡。f5=imread<'Fig0617<a><JLKMagenta>.tif.tif'>;g5=ice<'image',f5,'space','CMY'>;運(yùn)行結(jié)果：例5.9基于直方圖的映射。f5=imread<'Fig0618<a><CasterOriginal>.tif.tif'>;imshow<f5>;g5=ice<'image',f5,'space','hsi'>;運(yùn)行結(jié)果：例5.10彩色圖像平滑。clc,clear,closeall;f=imread<'Fig0619<a><RGB_iris>.tif'>;R=f<:,:,1>;G=f<:,:,2>;B=f<:,:,3>;subplot141;imshow<f>;title<'<a>原圖'>;subplot142;imshow<R>;title<'<a>R'>;subplot143;imshow<G>;title<'<b>G'>;subplot144;imshow<B>;title<'<c>B'>;h=rgb2hsi<f>;H=h<:,:,1>;S=h<:,:,2>;I=h<:,:,3>;figure;subplot131;imshow<H>;title<'<a>H'>;subplot132;imshow<S>;title<'<b>S'>;subplot133;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>;subplot131;imshow<f1>;title<'<a>分別平滑RGB分量得到平滑后RGB圖像'>;subplot132;imshow<F>;title<'<b>HSI中亮度分量I濾波后'>;subplot133;imshow<f2>;title<'<c>平滑所有三個(gè)HSI分量的結(jié)果'>;運(yùn)行結(jié)果：例5.11彩色圖像銳化。clc,clear,closeall;f=imread<'Fig0619<a><RGB_iris>.tif'>;w=fspecial<'average',5>;fb=imfilter<f,w,'replicate'>;subplot121;imshow<fb>;title<'<a>模糊圖像'>;lapmask=[111;1-81;111];fb=tofloat<fb>;fen=fb-imfilter<fb,lapmask,'replicate'>;subplot122;imshow<fen>;title<'<b>使用拉普拉斯算子增強(qiáng)后的圖像'>;運(yùn)行結(jié)果：例5.11使用函數(shù)colorgrad檢測RGB圖像的邊緣。clc,clear,closeall;R=imread<'Fig0624<a>RGB2-red.tif'>;G=imread<'Fig0624<b>RGB2-green.tif'>;B=imread<'Fig0624<c><RGB2-blue>.tif'>;f=imread<'Fig0624<d><RGB2-fullcolor>.tif'>;subplot221;imshow<R>;title<'<a>RGB圖像R分量'>;subplot222;imshow<G>;title<'<b>RGB圖像G分量'>;subplot223;imshow<B>;title<'<c>RGB圖像B分量'>;subplot224;imshow<f>;title<'<d>RGB圖像'>;figure;[VG,A,PPG]=colorgrad<f>;subplot121;imshow<VG>;title<'<e>在RGB向量空間中直接計(jì)算的梯度'>;subplot122;imshow<PPG>;title<'<f>分別計(jì)算每幅RGB分量圖像的二維梯度并將結(jié)果相加得到的合成梯度'>;figure;f1=imread<'Fig0604<a><iris>.tif'>;[VG1,A1,PPG1]=colorgrad<f1>;sub=abs<VG1-PPG1>;m1=max<sub>;%絕對(duì)差映射到[0,1]區(qū)間m2=max<m1>;sub=sub*1/m2;subplot141;imshow<f1>;title<'<a>RGB圖像'>;subplot142;imshow<VG1>;title<'<b>在RGB向量空間中直接計(jì)算的梯度'>;subplot143;imshow<PPG1>;title<'<c>通過分別計(jì)算每幅RGB分量圖像的二維梯度并將結(jié)果相加得到的合成梯度'>;subplot144;imshow<sub>;title<'<d>圖<b><c>的絕對(duì)差,標(biāo)定區(qū)間[0,1]'>;運(yùn)行結(jié)果：例5.11RGB彩色圖像分割。clc,clear,closeall;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;subplot121;imshow<f>;title<'<a>木衛(wèi)Io表面區(qū)域的偽彩色圖像'>;subplot122;imshow<g>;title<'<b>使用函數(shù)roipoly交互提取的感興趣區(qū)域'>;[M,N,K]=size<g>;%提取ROI區(qū)域中的點(diǎn)的坐標(biāo)I=reshape<g,M*N,3>;%g中彩色像素重新排列為I的行idx=find<mask>;%找出非黑彩色分量的標(biāo)準(zhǔn)差的倍數(shù)I=double<I<idx,1:3>>;[C,m]=covmatrix<I>;%C的主對(duì)角線上包含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;subplot241;imshow<E25>;title<'<a>T為25時(shí)‘