www.pudn.com > colorseg.zip > convProd.m, change:2003-03-07,size:284b

function Y = convProd(X,H)
% CONVPROD Product, not sum, convolution.  Note that H must be 3x3.
 
  
[xrows,xcols] = size(X);

Y = zeros(xrows,xcols);

for i=2:(xrows-1)
  for j=2:(xcols-1)
    [junk1,junk2,v] = find(H.*X((i-1):(i+1),(j-1):(j+1)));
    Y(i,j) = prod(v);
  end
end