shapesphere.m

level set matlab code

function data = shapeSphere(grid, center, radius)
% shapeSphere: implicit surface function for a sphere.
%
%   data = shapeSphere(grid, center, radius)
%
% Creates an implicit surface function (actually signed distance) 
%   for a sphere.
%
% Can be used to create circles in 2D or intervals in 1D.
%
% parameters:
%   grid	Grid structure (see processGrid.m for details).
%   center      Vector (length grid.dim) specifying point at the center 
%                 of the sphere (defaults to the origin).
%   radius	Scalar specifying the radius of the sphere
%                 (defaults to 1).
%
%   data	Output data array (of size grid.size) containing the
%                 implicit surface function.

% Copyright 2004 Ian M. Mitchell ([e-mail]).
% This software is used, copied and distributed under the licensing 
%   agreement contained in the file LICENSE in the top directory of 
%   the distribution.
%
% Ian Mitchell, 6/23/04

%---------------------------------------------------------------------------
% Default parameter values.
if(nargin < 2)
  center = zeros(grid.dim, 1);
end
if(nargin < 3)
  radius = 1;
end

%---------------------------------------------------------------------------
% Signed distance function calculation.
data = (grid.xs{1} - center(1)).^2;
for i = 2 : grid.dim
  data = data + (grid.xs{i} - center(i)).^2;
end
data = sqrt(data) - radius;

%---------------------------------------------------------------------------
% Warn the user if there is no sign change on the grid
%  (ie there will be no implicit surface to visualize).
if(all(data(:) < 0) || (all(data(:) > 0)))
  warning([ 'Implicit surface not visible because function has ' ...
            'single sign on grid' ]);
end