clear;close?all;
rgb_image=imread（‘E:\9.jpg‘）;
z=~（~rgb_image）;
ycbcr_image=double（rgb2ycbcr（rgb_image））;
ycbcr_image=ycbcr_image.*z;
[rcz]=size（ycbcr_image）;
hsv_image=rgb2hsv（rgb_image）;
a=25.39;b=14.03;ecx=1.60;ecy=2.41;cx=109.38;cy=152.02;x=0;y=0;
sinr=sin（2.53）;cosr=cos（2.53）;
for?i=1:1:r
????for?j=1:1:c???????
????????imgy=ycbcr_image（ij1）;
????????cb=ycbcr_image（ij2）;
????????cr=ycbcr_image（ij3）;
????????x=cosr*（cb-cx）+sinr*（cr-cy）;
????????y=-sinr*（cb-cx）+cosr*（cr-cy）;
????????if?imgy<60?ycbcr_image（ij:）=0;
????????elseif?imgy>=60&&imgy<=230
????????????if?（（x-ecx）^2）/（a^2）+（（y-ecy）^2）/（b^2）>1
????????????????ycbcr_image（ij:）=0;
????????????end
????????else?
????????????if?（（x-ecx）^2）/（（1.2*a）^2）+（（y-ecy）^2）/（（1.2*b）^2）>1
????????????????ycbcr_image（ij:）=0;
????????????end
????????end??????????
????end
end
subplot（121）;imshow（rgb_image）;
subplot（122）;imshow（ycbcr_image）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

????.CA....???????937??2010-04-15?20:40??CbCrEllipse\CbCrEllipse.asv

????.CA....???????938??2010-04-15?20:59??CbCrEllipse\CbCrEllipse.m

????.C.D...?????????0??2010-04-21?10:49??CbCrEllipse

-----------?---------??----------?-----??----

?????????????????1875????????????????????3