stencilxgridfrvc.m

来自「用matlab实现图像融合」· M 代码 · 共 26 行

function F = stencilxgridfRVC(G, Stencil, center)
%------------------------------------------------------------------------------
% Multiplies gridfunction G with Stencil.
% Excess area is filled by Vertex-Centered (VC) Reflection across boundaries.
%
% Design and implementation by:
% Dr. Paul M. de Zeeuw <[e-mail]>  http://homepages.cwi.nl/~pauldz/
% Last Revision: June 26, 2000.
% Copyright 1999-2002 Stichting CWI, Amsterdam.
%------------------------------------------------------------------------------
o=[0 0];
if ~all(size(o) == size(center))
  error(' stencilxgridfRVC - unexpected dimensions of center ')
end
F=zeros(size(G));
[n, m] = size(Stencil);
for i=1:n
   for j=1:m
      sij = Stencil(i,j);
      if sij ~= 0   
        F = F + rayxgridfRVC(G, ([i j] - center), sij);
      end
   end
end   
%------------------------------------------------------------------------------