--- 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) + + + + + +