彩色图像处理部分程序.pdfVIP

例 9.1 clear rgb_R=zeros(128,128); rgb_R(1:64,1:64)=1; rgb_G=zeros(128,128); rgb_G(1:64,65:128)=1; rgb_B=zeros(128,128); rgb_B(65:128,1:64)=1; rgb=cat(3,rgb_R,rgb_G,rgb_B); figure,imshow(rgb),title 例 9.2 rgb=imread(i_flower673.jpg); rgb1=im2double(rgb); r=rgb1(:, :,1); g=rgb1(:, :,2); b=rgb1(:,:,3); I=(r+g+b)/3; tmp1=min(min(r,g),b); tmp2=r+g+b; tmp2(tmp2==0)=eps; S=1-3.*tmp1./tmp2; tmp1=0.5*((r-g)+(r-b)); tmp2=sqrt((r-g).^2+(r-b).*(g-b)); theta=acos(tmp1./(tmp2+eps)); H=theta; H(bg)=2*pi-H(bg); H=H/(2*pi); H(S==0)=0; 例 9.3 I=imread(i_peppers_gray.bmp); GS8=graysl