用MATLAB怎么把YCbCr颜色空间的亮度分离啊？

/**

*颜色选择器颜色转换算法类

*<p>Copyright: langzisoft 13/10/2004</p>

* 山东大学计算机系

* @author 邵建

* @version 1.0

**/

import java.lang.Math

import java.awt.*

public class ColorConverter{

public ColorConverter(){

}

//RGB转换CMYk

public int[] rgb2cmyk(int R,int G,int B){

int cmyk []= new int[4]

cmyk[3]=(int)(Math.min(Math.min(255-R,255-G),255-B)/2.55)//cmykK

int MyR = (int)(R/2.55)

int Div = 100-cmyk[3]

if (Div == 0)Div = 1

cmyk[0] = ((100-MyR-cmyk[3])/Div)*100//cmykC

int MyG = (int)(G/2.55)

cmyk[1] = ((100-MyG-cmyk[3])/Div)*100

int MyB = (int)(B/2.55)

cmyk[2] = ((100-MyB-cmyk[3])/Div)*100

return cmyk

}

public Color cmyk2rgb(int C,int M,int Y,int K){

float MyC = C/100

float MyM = M/100

float MyY = Y/100

float MyK = K/100

int R = (int)((1-(MyC*(1-MyK)+MyK))*255)

int G = (int)((1-(MyM*(1-MyK)+MyK))*255)

int B = (int)((1-(MyY*(1-MyK)+MyK))*255)

if (R<0) R=0

if (G<0) G=0

if (B<0) B=0

if (R>255) R=255

if (G>255) G=255

if (B>255) B=255

Color rgb = new Color(R,G,B)

return rgb

}

public int[] rgb2yuv(int R,int G,int B){

int yuv[] = new int[3]

yuv[0] =(int)(0.299*R+0.587*G+0.114*B)

yuv[1] = (int)(-0.147*R-0.289*G+0.437*B)

yuv[2] = (int)(0.615*R-0.515*G-0.1*B)

return yuv

}

public Color yuv2rgb(int y,int u,int v){

int R = (int)(y+1.14*v)

int G = (int)(y-0.394*u-0.581*v)

int B = (int)(y+2.028*u)

if (R<0) R=0

if (G<0) G=0

if (B<0) B=0

if (R>255) R=255

if (G>255) G=255

if (B>255) B=255

Color rgb=new Color(R,G,B)

return rgb

}

public int[] rgb2yiq(int R,int G,int B){

int yiq[] = new int[3]

yiq[0] = (int)(0.299*R+0.587*G+0.114*B)

yiq[1] = (int)(0.596*R-0.274*G-0.322*B)

yiq[2] = (int)(0.212*R-0.523*G+0.311*B)

return yiq

}

public Color yiq2rgb(float y,float i,float q){

int R = (int)(y+0.956*i+0.621*q)

int G = (int)(y-0.272*i-0.647*q)

int B = (int)(y-1.105*i+1.702*q)

if (R<0) R=0

if (G<0) G=0

if (B<0) B=0

if (R>255) R=255

if (G>255) G=255

if (B>255) B=255

Color rgb = new Color(R,G,B)

return rgb

}

//RGB转换YCbCr

public int[] rgb2ycbcr(int R, int G,int B){

int ycbcr[] = new int[3]

ycbcr[0] = (int)(0.299*R+0.587*G+0.114*B)

ycbcr[1] = (int)(-0.1687*R-0.3313*G+0.5*B+128)

if (ycbcr[1] >255) ycbcr[1] = 255

ycbcr[2] = (int)(0.5*R-0.4187*G-0.0813*B+128)

if (ycbcr[2] >255) ycbcr[2] = 255

return ycbcr

}

public Color ycbcr2rgb(int Y,int Cb,int Cr){

int MyR = (int)(Y + 1.402 * (Cr -128))

int MyG = (int)(Y - 0.34414 * (Cb - 128) - 0.71414 * (Cr - 128))

int MyB = (int)(Y + 1.772 * (Cb -128))

if (MyR >255) MyR = 255

if (MyG >255) MyG = 255

if (MyB >255) MyB = 255

if (MyR <0) MyR =0

if (MyG <0) MyG =0

if (MyB <0) MyB =0

Color rgb = new Color(MyR,MyG,MyB)

return rgb

}

}

close allclearclc

%读取并显示原图

image_rgb=imread('woman1.jpg')

image_ycbcr=rgb2ycbcr(image_rgb)%将RGB色彩值变换为YcbCr色彩空间（将RGB真彩色图像转化为YcbCr色彩空间中相等的图像）

y=image_ycbcr(:,:,1) %这三行分别是Ycbcr空间的y,cb,cr对应矩阵

cr=image_ycbcr(:,:,2)

cb=image_ycbcr(:,:,3)

[W,H]=size(y)% y对应矩阵大小

image_face=zeros(W,H)%创建和原图等大的矩阵用来存放接下来二值化人脸数据

%以下循环是根据设定的阈值把人脸转换成二值图像

for x1=1:W %%%%%%%%(cb=( 100,140),cr=(70,160)这是阈值，这个阈值应该是人脸肤色的范围

for y1=1:H

if (cb(x1,y1)>=70)&&(cb(x1,y1)<=140)&&(cr(x1,y1)>=80)&&(cr(x1,y1)<=170)

image_face(x1,y1)=0

else

image_face(x1,y1)=255

end

end

end

figure(1),imshow(image_face) %显示初步转换成二值图像的人脸

%%%%%%%%% 找到人脸区域并设为白色%%%% (r=35)

se1=strel('square',35) %这个函数可以创建边长35的方形元素

fo=imclose(image_face,se1) %利用上面创建的方形元素，弥补figure（1）中人脸中以及其他部分残留的小块黑色

figure(2),imshow(fo) %显示修补漏洞后的人脸

log=logical(fo) %生成一个逻辑索引矩阵，为非零元素时返回1，否则返回0。用于把figure（2）中黑色部分记录为0，这样对非人脸部分做修改也不会影响人脸

r=image_rgb(:,:,1) %这个是rgb空间三维的数据

g=image_rgb(:,:,2)

b=image_rgb(:,:,3)

r(log==0)=255 %根据前面的逻辑数据把之前的黑色部分对应位置设为255，以后显示RGB图时，非人脸部分就是白色啦

g(log==0)=255

b(log==0)=255

face_rgb=cat(3,r,g,b)%把修改后的RGB矩阵重组

figure(3),imshow(face_rgb),title('face_rgb')% 显示非人脸部分已经转白的RGB图

afm = rgb2gray(face_rgb) % rgb图转灰度图

figure(4), imshow(afm), title('surface image') % 显示人脸灰度图

这个是我个人加的注释 不知道你能不能理解 不能理解可以追问 我们一起学习下 哈哈 我也是新手

---