--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/classes/@data3D/applymethod.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,89 @@
+% APPLYMETHOD applys the given method to the input 3D data.
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+% DESCRIPTION: APPLYMETHOD applys the given method to the input 3D data.
+%
+% CALL:        d = applymethod(d, pl)
+%
+% INPUTS:      d      - a 3D data object (xyzdata)
+%              pl     - a plist of configuration options
+%
+% PARAMETERS:
+%
+%       'method' - the method to apply to the data
+%       'axis'   - which axis vector to apply the method to. Possible values
+%                  are: 'X', 'Y', 'Z', 'XYZ' [default: 'Z']
+%       'option' - any additional option to pass to the method.
+%
+% VERSION:     $Id: applymethod.m,v 1.3 2011/04/17 09:19:17 hewitson Exp $
+%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+function pl = applymethod(ds, pl, method, getDefaultPlist, varargin)
+  
+  % Get function handles
+  dxFcn = {};
+  for jj = 1:numel(varargin)
+    if isa(varargin{jj}, 'function_handle')
+      dxFcn = varargin{jj};
+    end
+    if iscell(varargin{jj})
+      list = varargin{jj};
+      if ~isempty(list) && isa(list{1}, 'function_handle')
+        dxFcn = list{1};
+      end
+    end
+  end
+  
+  pl = applyDefaults(getDefaultPlist('3D'), pl);
+
+  % Get the axis we are dealing with
+  axis = find(pl, 'axis');
+  % Get the dimension to operate along
+  dim = find(pl, 'dim');
+  % Get any additional option
+  opt = find(pl, 'option');
+
+  % Loop over data objects
+  for jj=1:numel(ds)
+    switch lower(axis)
+      case 'x'
+        ds(jj).x  = apply(ds(jj).getX, method, dim, opt);
+      case 'y'
+        ds(jj).y  = apply(ds(jj).y, method, dim, opt);
+      case 'z'
+        ds(jj).z  = apply(ds(jj).z, method, dim, opt);
+      case 'xyz'
+        ds(jj).x  = apply(ds(jj).getX, method, dim, opt);
+        ds(jj).y  = apply(ds(jj).y, method, dim, opt);
+        ds(jj).z  = apply(ds(jj).z, method, dim, opt);
+      otherwise
+        error('### Unknown axis to operate on.');
+    end
+  end
+
+end
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%                               Local Functions                               %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+%-----------------------------------------------
+% Apply method to the vector v
+%-----------------------------------------------
+function v = apply(v, method, dim, opt)
+  if ~isempty(dim) && ~isempty(opt)
+    % User supplied a dimension and an option
+    v = feval(method, v, dim, opt);
+  elseif ~isempty(dim)
+    % User supplied only a dimension
+    v = feval(method, v, dim);
+  elseif ~isempty(opt)
+    % User supplied only an option
+    v = feval(method, v, opt);
+  else
+    % User supplied only a method
+    v = feval(method, v);
+  end
+end
+