elec_fit_ellipse_optim - Optimization for elec_fit_ellipse.m Called from elec_fit_ellipse.m
0001 function [f] = elec_fit_ellipse_optim(r, X, Y, Z, xo, yo, zo) 0002 0003 % elec_fit_ellipse_optim - Optimization for elec_fit_ellipse.m 0004 % 0005 % Called from elec_fit_ellipse.m 0006 % 0007 0008 % $Revision: 1.1 $ $Date: 2004/11/12 01:32:33 $ 0009 0010 % Licence: GNU GPL, no implied or express warranties 0011 % History: 02/2002, Darren.Weber_at_radiology.ucsf.edu 0012 % 0013 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0014 0015 % r is a 3x1 vector of radius values for each x,y,z axis component of ellipse 0016 % 0017 % equation of ellipsoid with center (xo,yo,zo) and radius for each axis (x,y,z) = (a,b,c): 0018 % (( x - xo )^2 / a^2) + (( y - yo )^2 / b^2) + (( z - zo )^2 / c^2) = 1 0019 % 0020 % This function below creates a scalar value to 0021 % return to the fminsearch function in elec_fit_sphere. 0022 0023 E = ( (X-xo).^2 )/r(1).^2 + ((Y-yo).^2)/r(2).^2 + ((Z-zo).^2)/r(3).^2 - 1; 0024 0025 f = sum( E .* E ); % sum of squares returned 0026 0027 return