数字图像处理结课论文(9)

figure(); imshow(rgb);

lapMatrix=[1 1 1;1 -8 1;1 1 1];

fR_filtered=imfilter(fR,lapMatrix,'replicate'); fG_filtered=imfilter(fG,lapMatrix,'replicate'); fB_filtered=imfilter(fB,lapMatrix,'replicate'); figure();

imshow(fR_filtered); figure();

imshow(fG_filtered); figure();

imshow(fB_filtered);

rgb_tmp=cat(3,fR_filtered,fG_filtered,fB_filtered); rgb_filtered=imsubtract(rgb,rgb_tmp); figure();

imshow(rgb_filtered); I1=(fR+fG+fB)/3;

tmp1=min(min(fR,fG),fB); tmp2=fR+fG+fB; tmp2(tmp2==0)=eps; S=1-3.*tmp1./tmp2;

tmp1=0.5*((fB-fG)+(fR-fB));

tmp2=sqrt((fR-fG).^2+(fR-fB).*(fG-fB)); theta=acos(tmp1./(tmp2+eps)); H1=theta;

H1(fB>fG)=2*pi-H1(fB>fG); H1=H1/(2*pi); H1(S==0)=0; figure(); imshow(H1); figure(); imshow(S); figure(); imshow(I1);

lapMatrix=[1 1 1;1 -8 1;1 1 1]; I=imfilter(I1,lapMatrix,'replicate'); figure();

imshow(I);%I分量锐化

实验结果：