Source Code for Module mvpa.misc.fsl.melodic

 1  # emacs: -*- mode: python; py-indent-offset: 4; indent-tabs-mode: nil -*- 
 2  # vi: set ft=python sts=4 ts=4 sw=4 et: 
 3  ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ## 
 4  # 
 5  #   See COPYING file distributed along with the PyMVPA package for the 
 6  #   copyright and license terms. 
 7  # 
 8  ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ## 
 9  """Wrapper around the output of MELODIC (part of FSL)""" 
10   
11  __docformat__ = 'restructuredtext' 
12   
13  import os 
14  import numpy as N 
15   
16  from mvpa.base import externals 
17  if externals.exists('nifti', raiseException=True): 
18      import nifti 
19   
20   
21 -class MelodicResults( object ): 
22      """Easy access to MELODIC output. 
23   
24      Only important information is available (important as judged by the 
25      author). 
26      """ 
27 -    def __init__( self, path ): 
28          """Reads all information from the given MELODIC output path. 
29          """ 
30          self.__outputpath = path 
31          self.__icapath = os.path.join( path, 'filtered_func_data.ica' ) 
32          self.__ic = \ 
33              nifti.NiftiImage( os.path.join( self.__icapath, 
34                                              'melodic_IC' ) ) 
35          self.__funcdata = \ 
36              nifti.NiftiImage( os.path.join( self.__outputpath, 
37                                              'filtered_func_data' ) ) 
38          self.__tmodes = N.fromfile( os.path.join( self.__icapath, 
39                                                    'melodic_Tmodes' ), 
40                                      sep = ' ' ).reshape( self.tr, self.nic ) 
41          self.__smodes = N.fromfile( os.path.join( self.__icapath, 
42                                                    'melodic_Smodes' ), 
43                                      sep = ' ' ) 
44          self.__icstats = N.fromfile( os.path.join( self.__icapath, 
45                                                     'melodic_ICstats' ), 
46                                       sep = ' ' ).reshape( self.nic, 4 ) 
47   
48      # properties 
49      path     = property( fget=lambda self: self.__respath ) 
50      ic       = property( fget=lambda self: self.__ic ) 
51      nic      = property( fget=lambda self: self.ic.extent[3] ) 
52      funcdata = property( fget=lambda self: self.__funcdata ) 
53      tr       = property( fget=lambda self: self.funcdata.extent[3] ) 
54      tmodes   = property( fget=lambda self: self.__tmodes ) 
55      smodes   = property( fget=lambda self: self.__smodes ) 
56      icastats = property( fget=lambda self: self.__icstats ) 
57      relvar_per_ic  = property( fget=lambda self: self.__icstats[:, 0] ) 
58      truevar_per_ic = property( fget=lambda self: self.__icstats[:, 1] ) 
59