Diff of /ants/utils/mni2tal.py [000000] .. [5d12a0]

Switch to side-by-side view

--- a
+++ b/ants/utils/mni2tal.py
@@ -0,0 +1,59 @@
+
+__all__ = ['mni2tal']
+
+def mni2tal(xin):
+    """
+    mni2tal for converting from ch2/mni space to tal - very approximate.
+
+    This is a standard approach but it's not very accurate.
+
+    ANTsR function: `mni2tal`
+
+    Arguments
+    ---------
+    xin : tuple
+        point in mni152 space.
+
+    Returns
+    -------
+    tuple
+
+    Example
+    -------
+    >>> import ants
+    >>> ants.mni2tal( (10,12,14) )
+
+    References
+    ----------
+    http://bioimagesuite.yale.edu/mni2tal/501_95733_More\\%20Accurate\\%20Talairach\\%20Coordinates\\%20SLIDES.pdf
+    http://imaging.mrc-cbu.cam.ac.uk/imaging/MniTalairach
+    """
+    if (not isinstance(xin, (tuple,list))) or (len(xin) != 3):
+        raise ValueError('xin must be tuple/list with 3 coordinates')
+
+    x = list(xin)
+    # The input image is in RAS coordinates but we use ITK which returns LPS
+    # coordinates.  So we need to flip the coordinates such that L => R and P => A to
+    # get RAS (MNI) coordinates
+    x[0] = x[0] * (-1)  # flip X
+    x[1] = x[1] * (-1)  # flip Y
+    
+    xout = x
+    
+    if (x[2] >= 0):
+        xout[0] = x[0] * 0.99
+        xout[1] = x[1] * 0.9688 + 0.046 * x[2]
+        xout[2] = x[1] * (-0.0485) + 0.9189 * x[2]
+    
+    if (x[2] < 0):
+        xout[0] = x[0] * 0.99
+        xout[1] = x[1] * 0.9688 + 0.042 * x[2]
+        xout[2] = x[1] * (-0.0485) + 0.839 * x[2]
+    
+    return(xout)
+
+
+
+
+
+