enh: add SHORE model

Resolves: #38.

setup.cfg

@@ -25,6 +25,7 @@ setup_requires =
     setuptools_scm >= 6.2
     wheel
 install_requires =
+    cvxpy
     dipy>=1.3.0
     joblib
     nipype>= 1.5.1, < 2.0

src/eddymotion/model.py

@@ -44,7 +44,7 @@ def init(gtab, model="DTI", **kwargs):
             return AverageDWModel(gtab=gtab, **kwargs)
 
         # Generate a GradientTable object for DIPY
-        if model.lower() in ("dti", "dki"):
+        if model.lower() in ("dti", "dki", "shore"):
             Model = globals()[f"{model.upper()}Model"]
             return Model(gtab, **kwargs)
 
@@ -285,13 +285,30 @@ class DTIModel(BaseModel):
     _model_class = "dipy.reconst.dti.TensorModel"
 
 
-class DKIModel(BaseModel):
+class DKIModel(DTIModel):
     """A wrapper of :obj:`dipy.reconst.dki.DiffusionKurtosisModel`."""
 
-    _modelargs = DTIModel._modelargs
     _model_class = "dipy.reconst.dki.DiffusionKurtosisModel"
 
 
+class SHOREModel(BaseModel):
+    """A wrapper of :obj:`dipy.reconst.shore.ShoreModel`."""
+
+    _modelargs = (
+        "radial_order",
+        "zeta",
+        "lambdaN",
+        "lambdaL",
+        "tau",
+        "constrain_e0",
+        "positive_constraint",
+        "pos_grid",
+        "pos_radius",
+        "cvxpy_solver",
+    )
+    _model_class = "dipy.reconst.shore.ShoreModel"
+
+
 def _rasb2dipy(gradient):
     gradient = np.asanyarray(gradient)
     if gradient.ndim == 1: