New functions: ribbon.m spinmap.m

[Kai Habel]

Hello all,

I have two new graphic functions implemented: ribbon and spinmap.


Examples
* ribbon:

[x,y,Z] = sombrero;
[X,Y] = meshgrid(x,y);
ribbon(Y,Z);


* spinmap:

image;spinmap

For more complex graphics like surf(peaks) spinmap is not very fast, but
might be useful nevertheless.

Kai

ChangeLog:

2007-11-12  Kai Habel  <address@hidden>

  * plot/spinmap.m, plot/ribbon.m: New functions.
  * plot/Makefile.in (SOURCES): Add them to the list.

--- Makefile.in 2007-11-12 20:56:20.000000000 +0100
+++ Makefile.in.new     2007-11-12 20:55:43.000000000 +0100
@@ -108,6 +108,7 @@
   print.m \
   quiver.m \
   replot.m \
+  ribbon.m \
   semilogx.m \
   semilogxerr.m \
   semilogy.m \
@@ -115,6 +116,7 @@
   shading.m \
   shg.m \
   sombrero.m \
+  spinmap.m \
   stairs.m \
   stem.m \
   subplot.m \

## Copyright (C) 2007 Kai Habel
##
## This program is free software; you can redistribute it and/or modify it
## under the terms of the GNU General Public License as published by
## the Free Software Foundation; either version 2, or (at your option)
## any later version.
##
## OctPlot is distributed in the hope that it will be useful, but
## WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
## General Public License for more details.
##
## You should have received a copy of the GNU General Public License
## along with OctPlot; see the file COPYING.  If not, write to the Free
## Software Foundation, 59 Temple Place - Suite 330, Boston, MA
## 02111-1307, USA.

## -*- texinfo -*-
## @deftypefn {Function File}  ribbon (@var{X}, @var{Y}, @var{WIDTH})
## @deftypefnx {Function File}  ribbon (@var{X}, @var{Y})
## @deftypefnx {Function File}  ribbon (@var{Y})
## @deftypefnx {Function File}  @var{h} = ribbon (...)
## Plots ribbon plot for the columns of @var{Y} vs. @var{X}. The optional 
## parameter @var{WIDTH} specifies the width of a single ribbon (default is 
0.75).
## If @var{X} is omitted, a vector containing the row numbers is assumed 
(1:rows(Y)).
## If requested a vector @var{h} of the handles to the surface objects is 
returned.
## @end deftypefn
## @seealso{gca, colorbar}

## Author: Kai Habel <kai.habel at gmx.de>

function h = ribbon(X, Y, W)

  newplot ();

  if (nargin == 1)
    Y = X;
    if (isvector(Y))
      Y = Y(:);
    endif
    [nr, nc] = size(Y);
    X = repmat((1 : nr)', 1, nc);
    W = 0.75;
  elseif (nargin == 2)
    W = 0.75;
  elseif (nargin == 3)
  else
    print_usage();
  end

  if (isvector(X) && isvector(Y))
    if (length(X) != length(Y))
      error("In case of vectors, X and Y must have same length")
    else
      [X, Y] = meshgrid(X, Y);
    endif
  else
    if (!all(size(X) == size(Y)))
      error("In case of matrices, X and Y must have same size")
    endif
  endif

  [nr,nc] = size(Y);
  tmp = zeros(1,nc);

  for c = nc:-1:1
    ZZ = [Y(:,c) Y(:,c)];
    y = X(:,c);
    x = [c - W / 2, c + W / 2];
    [XX,YY] = meshgrid(x,y);
    CC = ones(size(ZZ))*c;
    tmp(c) = surface(XX,YY,ZZ,CC);
  endfor

  ax = get (tmp(c), "parent");

  if (!ishold ())
    set (ax, "view", [-37.5, 30], 
"box","off","xgrid","on","ygrid","on","zgrid","on");
  endif

  if (nargout > 0)
    h = tmp;
  endif
end

## Copyright (C) 2007 Kai Habel
##
## This program is free software; you can redistribute it and/or modify it
## under the terms of the GNU General Public License as published by
## the Free Software Foundation; either version 2, or (at your option)
## any later version.
##
## OctPlot is distributed in the hope that it will be useful, but
## WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
## General Public License for more details.
##
## You should have received a copy of the GNU General Public License
## along with OctPlot; see the file COPYING.  If not, write to the Free
## Software Foundation, 59 Temple Place - Suite 330, Boston, MA
## 02111-1307, USA.

## -*- texinfo -*-
## @deftypefn {Function File}  spinmap (@var{T}, @var{inc})
## @deftypefnx {Function File}  spinmap (@var{T})
## @deftypefnx {Function File}  spinmap ()
## Cycles the colormap for @var{T} seconds with an increment
## of @var{inc}. Both parameter are optional. In that case 5s and an increment 
## of 2 is taken.
## A higher @var{inc} causes a faster cycle through the colormap.
## @end deftypefn
## @seealso{gca, colorbar}

## Author: Kai Habel <kai.habel at gmx.de>

function spinmap(T, inc)

  if (nargin == 0)
    inc = 2;
    T = 5;
  elseif (nargin == 1)
    inc = 2;
  endif

  cmap = get(gcf,"Colormap");
  clen = rows(cmap);

  t0 = clock;

  while (etime(clock, t0) < T)
    for n = 1 : inc : clen
      newmap = shift(cmap, n, 1);
      set(gcf,"Colormap",newmap)
      drawnow
    endfor
  end
  set(gcf,"Colormap",cmap)
end