/Modules/Python/StochasticTractographyServerMP/autoTracto2regions.py
Python | 76 lines | 48 code | 22 blank | 6 comment | 12 complexity | 5ca82468da547d280635400565afa876 MD5 | raw file
- import os, sys, glob, time, shutil
- if __name__ == '__main__':
- rootA = os.path.abspath(os.path.dirname(sys.argv[0]))
- print 'Input : ', rootA
- cont = os.walk('.')
- d = cont.next()
- dirs = d[1]
- print 'Folders to visit : ', dirs
- for i in range(len(dirs)):
- print 'Go into folder : ', dirs[i]
- os.chdir(dirs[i])
-
- contRM = os.walk('.')
- dRM = contRM.next()
- dirsRM = dRM[1]
- if len(dirsRM)>0:
- for s in range(len(dirsRM)):
- os.system('rm -rf ' + dirsRM[s])
- os.chdir('..')
- print 'Cleaning done'
- for i in range(len(dirs)):
- print 'Go into folder : ', dirs[i]
- os.chdir(dirs[i])
-
- fullDir = os.getcwd()
-
- pars = glob.glob('*.in')
- dwis = glob.glob('*.dwi')
- rois = glob.glob('*.roi')
- wmas = glob.glob('*.wm')
- #try:
- # ret = os.system('python ' + rootA + '/TractoCell0/ModifyHeader.py ' + pars[0] + ' stepSize 0.9' )
- # ret = os.system('python ' + rootA + '/TractoCell0/ModifyHeader.py ' + pars[0] + ' totalTracts 40' )
- # ret = os.system('python ' + rootA + '/TractoCell0/ModifyHeader.py ' + pars[0] + ' maxLength 200')
- #except:
- # print "Exception: continue anyway!"
- if len(pars)==1 and len(dwis)==1:
-
- print 'ROIS : ', rois
- if len(rois)!=2:
- print 'No consistent number of ROIs'
- else:
- try:
- for s in range(len(rois)/2):
- tmpD = rois[s].split('.')[0] + '_' + rois[s+1].split('.')[0]
- isDir = os.access(tmpD, os.F_OK)
- if not isDir:
- os.mkdir(tmpD)
- tmpF = './' + tmpD + '/'
- os.chdir(tmpF)
- print sys.executable + ' ' + rootA + '/' + 'StochasticTractoGraphyClusterFiles.py ' + fullDir + '/' + pars[0] + ' ' + fullDir + '/' + dwis[0] + ' ' + fullDir + '/' + rois[s] + ' ' + fullDir + '/' + rois[s+1]
- ret = os.system(sys.executable + ' ' + rootA + '/' + 'StochasticTractoGraphyClusterFiles.py ' + fullDir + '/' + pars[0] + ' ' + fullDir + '/' + dwis[0] + ' ' + fullDir + '/' + rois[s] + ' ' + fullDir + '/' + rois[s+1])
- os.chdir('..')
- except:
- print "Exception: continue anyway!"
- os.chdir('..')
- time.sleep(4) # breathe
-
- print 'Completed'