Matlab圖像顏色空間轉(zhuǎn)換

1、精選文檔Matlab圖像顏色空間轉(zhuǎn)換試驗內(nèi)容用matlab軟件編程實現(xiàn)下述任務(wù)：讀入彩色圖像，提取其中的R、G、B顏色重量，并呈現(xiàn)出來。我們學習了多種表示圖像的顏色空間，請編寫程序?qū)D像轉(zhuǎn)換到Y(jié)UV、YIQ、YCrCb、HIS、CMY等顏色空間，并呈現(xiàn)出來。顏色空間的轉(zhuǎn)化關(guān)系參考以下公式：原始圖片三個色調(diào)重量YUV與RGB之間的轉(zhuǎn)換Y0.229R0.587G0.114BU0.147R0.289G+0.436BV=0.615R0.515G0.100BYIQ與RGB之間的轉(zhuǎn)換Y0.299R0.587G0.114BI0.596R0.275G0.321BQ=0.212R0.523G0.311BYCrC

2、b與RGB之間的轉(zhuǎn)換Y  = 0.2990R + 0.5870G + 0.1140BCr = 0.5000R - 0.4187G - 0.0813B + 128Cb = -0.1687R - 0.3313G + 0.5000B + 128HSI與RGB之間的轉(zhuǎn)換I（RGB）/3Harccos 0.5*(R-G)+(R-B) / (R-G)2 + (R-B)(G-B)0.5 S=1-min（R，G，B）/ I CMY與RGB之間的轉(zhuǎn)換心得體會查閱了很多資料，并且學習了關(guān)于matlab實現(xiàn)圖像顏色空間轉(zhuǎn)換的過程。不同的顏色空間在描述圖像的顏色時側(cè)重點不同。如RGB（紅、綠、藍三原色）顏色

3、空間適用于彩色監(jiān)視器和彩色攝像機，HSI（色調(diào)、飽和度、亮度）更符合人描述和解釋顏色的方式（或稱為HSV，色調(diào)、飽和度、亮度），CMY（青、深紅、黃）、CMYK（青、深紅、黃、黑。）主要針對彩色打印機、復(fù)印機等，YIQ（亮度、色差、色差）是用于NTSC規(guī)定的電視系統(tǒng)格式，YUV（亮度、色差、色差）是用于PAL規(guī)定的電視系統(tǒng)格式，YCbCr（亮度單一要素、藍色與參考值的差值、紅色與參考值的差值）在數(shù)字影像中廣泛應(yīng)用。近年來消滅了另一種顏色空間l，由于其把亮度和顏色信息最大限度的分別，在該顏色空間可以分別處理亮度或顏色而不相互影響。通過這次試驗，實現(xiàn)了五種顏色空間的轉(zhuǎn)換，看到了不同的絢麗結(jié)果,把握

4、了一些基本的學問。程序clearrgb=imread('G:LearningMultiMedia666.jpg');rgb2hsi(rgb);rgb_r=rgb(:,:,1);rgb_g=rgb(:,:,2);rgb_b=rgb(:,:,3);n, m = size(rgb);zero=zeros(n,m/3);Y = 0.229 * rgb_r + 0.587 * rgb_g + 0.114 * rgb_b;U = -0.147 * rgb_r - 0.289 * rgb_g + 0.436 * rgb_b;V = 0.615 * rgb_r - 0.515 * rgb_g

5、- 0.100 * rgb_b;I = 0.596 * rgb_r - 0.275 * rgb_g - 0.321 * rgb_b;Q = 0.212 * rgb_r - 0.523 * rgb_g + 0.311 * rgb_b;Cr = 0.5 * rgb_r - 0.4187 * rgb_g - 0.0813 * rgb_b + 128;Cb = -0.1687 * rgb_r - 0.3313 * rgb_g + 0.5 * rgb_b + 128;I = (rgb_r + rgb_g + rgb_b) / 3;R=cat(3,Y,zero,zero);G=cat(3,zero,U,z

6、ero);B=cat(3,zero,zero,V);RGB=cat(3, Y, Cr, Cb);%通過修改參數(shù)的值可以顯示各種顏色空間的效果imshow(RGB);subplot(2,2,1),imshow(R),title('紅色重量');subplot(2,2,2),imshow(G),title('綠色重量');subplot(2,2,3),imshow(B),title('藍色重量');subplot(2,2,4),imshow(RGB);HIS:function hsi=rgb2hsi(rgb)%提取單通道重量rgb=im2double

7、(rgb);r=rgb(:,:,1);g=rgb(:,:,2);b=rgb(:,:,3);%實現(xiàn)轉(zhuǎn)換num=0.5*(r-g)+(r-b);den=sqrt(r-g).2+(r-b).*(g-b);theta=acos(num./(den+eps);H=theta;H(b>g)=2*pi-H(b>g);H=H/(2*pi);num=min(min(r,g),b);den=r+g+b;den(den=0)=eps;S=1-3.*num./den;H(S=0)=0;I=(r+g+b)/3;hsi=cat(3,H,S,I);imshow(hsi);CMY:function hsi=rgb2CMY(rgb)rgb=im2