mexGetSPFBasisMatrix
====================

.. code-block:: none

  % mexGetSPFBasisMatrix: get SPF basis matrix with given spherical and radial orders, orientationMatrix, bVector
  %
  % Usage:
  %   spfMatrix = mexGetSPFBasisMatrix(shOrder, radialOrder, orientationMatrix, bVector)
  %   spfMatrix = mexGetSPFBasisMatrix(shOrder, radialOrder, orientationMatrix, bVector, params)
  % 
  % Inputs:
  %   shOrder:               SH order
  %   radialOrder:           order in radial part
  %   orientationMatrix :    orientation matrix where each row is an orientation (x,y,z), caresian format
  %   bVector:               b values in a vector, each row is a b value
  %   params.scale :         (optional) scale value for the SPF basis in radial part
  %   params.tau:            (optional) tau value to calculate qVector from bVector, b=4\pi^2 \tau q^2
  %   params.original:       (optional) true: (default) original SPF basis matrix 
  %                                     false: independent SPF basis matrix, considering E(0)=1
  %
  % Outputs:
  %   spfMatrix:              each column is a SPF basis
  %
  %
  % Reference: 
  %  1. "Model-Free, Regularized, Fast, and Robust Analytical Orientation Distribution Function Estimation",
  %     Jian Cheng, Aurobrata Ghosh, Rachid Deriche, Tianzi Jiang, Medical Image Computing and Computer-Assisted Intervention (MICCAI'10), vol. 6361, pp. 648–656, sep, 2010.
  %  2. "Model-free and Analytical EAP Reconstruction via Spherical Polar Fourier Diffusion MRI", 
  %     Jian Cheng, Aurobrata Ghosh, Tianzi Jiang, Rachid Deriche, Medical Image Computing and Computer-Assisted Intervention (MICCAI'10), vol. 6361, pp. 590–597, sep, 2010. 
  %
  % Copyright (c) 2013, Jian Cheng <jian.cheng.1983@gmail.com>
  %