divergencer.m

3D shape reconstruction matlab code. It used shape from defocus technique with diffusion. You can re

function [divF] = divergenceR(F,coeff)
% compute the divergence of the gradient of f
% multiplied by a diffusion coefficient
% backward difference scheme
%
% Copyright 2006 Paolo Favaro ([e-mail])
% 
% School of Engineering and Physical Sciences
% Heriot-Watt University, Edinburgh, UK
% 
% Last revision: August 2006
%
% This program can be used only for research purposes.
% This program is distributed WITHOUT ANY WARRANTY; 
% without even the implied warranty of MERCHANTABILITY 
% or FITNESS FOR A PARTICULAR PURPOSE.

[nri,nrj]=size(F);
Gx = zeros(nri,nrj);
Gy = zeros(nri,nrj);
Gx(:,1:nrj-1) = diff(F,1,2);
Gy(1:nri-1,:) = diff(F,1,1);
% diffusion coefficient
Gxcoeff = Gx.*coeff;
Gycoeff = Gy.*coeff;
Gxx(:,2:nrj) = diff(Gxcoeff,1,2);
Gyy(2:nri,:) = diff(Gycoeff,1,1);
divF = Gxx+Gyy;
return