gaussfilter.m

this code is skin color detection algorithm

function y = GaussFilter(RGBimg)
%% HSV
% Mean = [0.08;0.56;0.75]; Variance = diag([0.03;0.1;0.1]);
% imgH=rgb2hsv(RGBimg/255);
% [N,M,t] = size(RGBimg);
% imgHShaped = reshape(imgH,N*M,[])';
% t1 = imgHShaped - repmat(Mean,1,N*M); t2 = inv(Variance)*t1;
% HSVimgFiltered = exp(-0.5*sum(t1.*t2,1))/((2*pi)^3/2 * det(Variance)^1/2);
% y = reshape(HSVimgFiltered,N,M);

%% RGB
load 'skinrgb.mat';
Mean        = mean (skinrgb);
Variance    = cov(skinrgb);
RGBimg = RGBimg/255;

[N,M,dummy] = size(RGBimg);
RGBimgShaped = reshape(RGBimg,N*M,[])';
t1 = RGBimgShaped - repmat(Mean',1,N*M); t2 = inv(Variance)*t1;
RGBimgFiltered = exp(-0.5*sum(t1.*t2,1))/((2*pi)^3/2 * det(Variance)^1/2);
y = reshape(RGBimgFiltered',N,M,[]);
end