X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?p=CreaPhase.git;a=blobdiff_plain;f=octave_packages%2Fgeometry-1.5.0%2Fgeom2d%2FdrawEllipseArc.m;fp=octave_packages%2Fgeometry-1.5.0%2Fgeom2d%2FdrawEllipseArc.m;h=bd111e84c1f10f67033b67f46371cc447363e7ca;hp=0000000000000000000000000000000000000000;hb=f5f7a74bd8a4900f0b797da6783be80e11a68d86;hpb=1705066eceaaea976f010f669ce8e972f3734b05

diff --git a/octave_packages/geometry-1.5.0/geom2d/drawEllipseArc.m b/octave_packages/geometry-1.5.0/geom2d/drawEllipseArc.m
new file mode 100644
index 0000000..bd111e8
--- /dev/null
+++ b/octave_packages/geometry-1.5.0/geom2d/drawEllipseArc.m
@@ -0,0 +1,161 @@
+%% Copyright (c) 2011, INRA
+%% 2003-2011, David Legland <david.legland@grignon.inra.fr>
+%% 2011 Adapted to Octave by Juan Pablo Carbajal <carbajal@ifi.uzh.ch>
+%%
+%% All rights reserved.
+%% (simplified BSD License)
+%%
+%% Redistribution and use in source and binary forms, with or without
+%% modification, are permitted provided that the following conditions are met:
+%%
+%% 1. Redistributions of source code must retain the above copyright notice, this
+%%    list of conditions and the following disclaimer.
+%%     
+%% 2. Redistributions in binary form must reproduce the above copyright notice, 
+%%    this list of conditions and the following disclaimer in the documentation
+%%    and/or other materials provided with the distribution.
+%%
+%% THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+%% AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+%% IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+%% ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
+%% LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR 
+%% CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+%% SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS 
+%% INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+%% CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+%% ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+%% POSSIBILITY OF SUCH DAMAGE.
+%%
+%% The views and conclusions contained in the software and documentation are
+%% those of the authors and should not be interpreted as representing official
+%% policies, either expressed or implied, of copyright holder.
+
+%% -*- texinfo -*-
+%% @deftypefn {Function File} {@var{h} = } drawEllipseArc (@var{arc})
+%% Draw an ellipse arc on the current axis.
+%%
+%%   drawEllipseArc(ARC) 
+%%   draw ellipse arc specified by ARC. ARC has the format:
+%%     ARC = [XC YC A B THETA T1 T2]
+%%   or:
+%%     ARC = [XC YC A B T1 T2] (isothetic ellipse)
+%%   with center (XC, YC), main axis of half-length A, second axis of
+%%   half-length B, and ellipse arc running from t1 to t2 (both in degrees,
+%%   in Counter-Clockwise orientation).
+%%
+%%   Parameters can also be arrays. In this case, all arrays are suposed to
+%%   have the same size...
+%%
+%%   @example
+%%     % draw an ellipse arc: center = [10 20], radii = 50 and 30, theta = 45
+%%     arc = [10 20 50 30 45 -90 270];
+%%     figure;
+%%     axis([-50 100 -50 100]); axis equal;
+%%     hold on
+%%     drawEllipseArc(arc, 'color', 'r')
+%%
+%%     % draw another ellipse arc, between angles -60 and 70
+%%     arc = [10 20 50 30 45 -60 (60+70)];
+%%     figure;
+%%     axis([-50 100 -50 100]); axis equal;
+%%     hold on
+%%     drawEllipseArc(arc, 'LineWidth', 2);
+%%     ray1 = createRay([10 20], deg2rad(-60+45));
+%%     drawRay(ray1)
+%%     ray2 = createRay([10 20], deg2rad(70+45));
+%%     drawRay(ray2)
+%% @end example
+%%
+%%   @seealso{ellipses2d, drawEllipse, drawCircleArc}
+%% @end deftypefn
+
+function varargout = drawEllipseArc(varargin)
+
+  %% Extract input arguments
+
+  % extract dawing style strings
+  styles = {};
+  for i = 1:length(varargin)
+      if ischar(varargin{i})
+          styles = varargin(i:end);
+          varargin(i:end) = [];
+          break;
+      end
+  end
+
+  if length(varargin)==1
+      ellipse = varargin{1};
+      x0 = ellipse(1);
+      y0 = ellipse(2);
+      a  = ellipse(3);
+      b  = ellipse(4);
+      if size(ellipse, 2)>6
+          theta   = ellipse(5);
+          start   = ellipse(6);
+          extent  = ellipse(7);
+      else
+          theta   = zeros(size(x0));
+          start   = ellipse(5);
+          extent  = ellipse(6);
+      end
+      
+  elseif length(varargin)>=6
+      x0 = varargin{1};
+      y0 = varargin{2};
+      a  = varargin{3};
+      b  = varargin{4};
+      if length(varargin)>6
+          theta   = varargin{5};
+          start   = varargin{6};
+          extent  = varargin{7};
+      else
+          theta   = zeros(size(x0));
+          start   = varargin{5};
+          extent  = varargin{6};
+      end
+      
+  else
+      error('drawellipse: please specify center x, center y and radii a and b');
+  end
+
+
+  %% Drawing
+
+  % allocate memory for handles
+  h = zeros(size(x0));
+
+  for i = 1:length(x0)
+      % start and end angles
+      t1 = deg2rad(start);
+      t2 = t1 + deg2rad(extent);
+      
+      % vertices of ellipse
+      t = linspace(t1, t2, 60);
+      
+      % convert angles to ellipse parametrisation
+      sup = cos(t) > 0;
+      t(sup)  = atan(a(i) / b(i) * tan(t(sup)));
+      t(~sup) = atan2(a(i) / b(i) * tan(2*pi - t(~sup)), -1);
+      t = mod(t, 2*pi);
+      
+      % precompute cos and sin of theta (given in degrees)
+      cot = cosd(theta(i));
+      sit = sind(theta(i));
+
+      % compute position of points
+      xt = x0(i) + a(i)*cos(t)*cot - b(i)*sin(t)*sit;
+      yt = y0(i) + a(i)*cos(t)*sit + b(i)*sin(t)*cot;
+      
+      h(i) = plot(xt, yt, styles{:});
+  end
+
+
+  %% Process output arguments
+
+  if nargout > 0
+      varargout = {h};
+  end
+
+endfunction
+