Merge pull request #17 from pnlbwh/individual-fa

Individual fa

Author: tashrifbillah

TUTORIAL.md

@@ -265,6 +265,10 @@ as follows:
     --labelMap atlas.nii.gz --space MNI.nii.gz
 
 
+Unlike [original TBSS](https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/TBSS/UserGuide) approach, we use the [ENIGMA](http://enigma.ini.usc.edu/wp-content/uploads/DTI_Protocols/ENIGMA_ROI_protocol_USC.pdf) approach 
+that identifies the direction of projection onto the skeleton based on the individual FA maps rather than on the mean FA map.
+
+
 ## 3. --studyTemplate
 
 With this branch, a study-specific template is created using `antsMultivariateTemplateConstruction2.sh`. 
@@ -502,7 +506,33 @@ corresponding input images are FA.
 
 However, if you wish to run FA first and then other modalities in future, use option (iv) from above. There, you should 
 provide the directory containing warp/affine obtained during registration of subject FA to the template/standard space. 
-This way, we bypass doing the same non-linear registration once again.
+This way, we bypass doing the same non-linear registration once again. In addition, provide any templates created 
+during FA TBSS.Here are a few sample commands for running separate nonFA TBSS:
+
+    
+**partial** `--engima` TBSS
+    
+    $libDir/tbss_all -i MD/origdata,RD/origdata \
+    -c $CASELIST \                                      # same caselist that was used for FA TBSS
+    --xfrmDir $testDir/enigmaTemplateOutput/transform \ # transform files are obtained from here
+    --modality MD,RD --enigma \                         # --enigma tells to use enigma templates
+    --avg                                               # -o is not required
+
+
+
+**partial** `--studyTemplate` TBSS
+
+    $libDir/tbss_all -i AD/origdata \
+    -c $CASELIST \
+    --xfrmDir $testDir/studyTemplateOutput/template \
+    --modality AD \
+    --template $testDir/studyTemplateOutput/template/template0.nii.gz \     # provide created templates 
+    --templateMask $testDir/studyTemplateOutput/stats/mean_FA_mask.nii.gz \
+    --skeleton $testDir/studyTemplateOutput/stats/mean_FA_skeleton.nii.gz \
+    --skeletonMask $testDir/studyTemplateOutput/stats/mean_FA_skeleton_mask.nii.gz \
+    --skeletonMaskDst $testDir/studyTemplateOutput/stats/mean_FA_skeleton_mask_dst.nii.gz \
+    -s $FSLDIR/data/standard/FMRIB58_FA_1mm.nii.gz \                        # provide space defining image if wanted
+    -l $FSLDIR/data/atlases/JHU/JHU-ICBM-labels-1mm.nii.gz \                # provide atlas if wanted
 
 
 # List creation
@@ -616,7 +646,7 @@ Finally, if needed, you can copy the transform files in the [transform](#ii-tran
 Processing can be multi-threaded over the cases. Besides, `antsMultivariateTemplateConstruction2.sh` utilizes 
 multiple threads to speed-up template construction. 
 
-    --nproc 8 # default is 4, use -1 for all available
+    --ncpu 8 # default is 4, use -1 for all available
 
 However, multi-threading comes with a price of slowing down other processes that may be running in your system. So, it 
 is advisable to leave out at least two cores for other processes to run smoothly.

_version.py

@@ -1 +1 @@
-__version__ = '0.1.1'
+__version__ = '0.1.2'

lib/antsTemplate.py

@@ -35,6 +35,11 @@ def antsMult(caselist, outPrefix, logDir, N_proc, verbose):
                            caselist]), shell= True, stdout= f, stderr= sys.stdout)
 
 
+    if f.name!='<sys.stdout>':
+        f.close()
+
+
+
 def antsReg(fixedImg, movingImg, outPrefix, logDir, verbose):
 
     if verbose:
@@ -54,3 +59,6 @@ def antsReg(fixedImg, movingImg, outPrefix, logDir, verbose):
                            '-o', outPrefix,
                            '-e', '123456']), shell=True, stdout= f, stderr= sys.stdout)
 
+
+    if f.name!='<sys.stdout>':
+        f.close()

lib/project_skeleton.py

@@ -14,7 +14,7 @@
 from tbssUtil import pjoin, load, check_call
 
 
-def project_skeleton(c, imgPath, args, meanFA, skelDir):
+def project_skeleton(c, imgPath, args, skelDir):
 
     '''
     Part of FSL (ID: 5.0.11)
@@ -46,7 +46,7 @@ def project_skeleton(c, imgPath, args, meanFA, skelDir):
     if args.modality == 'FA':
 
         check_call((' ').join(['tbss_skeleton',
-                              '-i', meanFA,
+                              '-i', imgPath,
                               '-p', args.SKEL_THRESH, args.skeletonMaskDst, args.SEARCH_RULE_MASK,
                               imgPath, modImgSkel,
                               '-s', args.skeletonMask]),
@@ -55,7 +55,7 @@ def project_skeleton(c, imgPath, args, meanFA, skelDir):
     else:
 
         check_call((' ').join(['tbss_skeleton',
-                              '-i', meanFA,
+                              '-i', imgPath,
                               '-p', args.SKEL_THRESH, args.skeletonMaskDst, args.SEARCH_RULE_MASK,
                               pjoin(args.outDir, 'FA', 'warped', f'{c}_FA_to_target.nii.gz'),
                               modImgSkel, '-a', imgPath,

lib/skeletonize.py

@@ -11,7 +11,7 @@
 # View LICENSE at https://github.com/pnlbwh/tbss/blob/master/LICENSE
 # ===============================================================================
 
-from tbssUtil import load, save_nifti, pjoin, check_call, environ, basename, Pool
+from tbssUtil import load, save_nifti, pjoin, check_call, environ, basename, Pool, psplit
 import numpy as np
 from project_skeleton import project_skeleton
 
@@ -212,8 +212,7 @@ def skeletonize(imgs, cases, args, statsDir, skelDir, miFile):
     allFAskeletonizedData= np.zeros((len(imgs), X, Y, Z), dtype= 'float32')
     pool= Pool(args.ncpu)
     for c, imgPath in zip(cases, imgs):
-
-        res.append(pool.apply_async(project_skeleton, (c, imgPath, args, meanFA, skelDir)))
+            res.append(pool.apply_async(project_skeleton, (c, imgPath, args, skelDir)))
 
 
     for i,r in enumerate(res):

lib/tbss_single.py

@@ -19,7 +19,6 @@
 from glob import glob
 from plumbum.cmd import antsApplyTransforms, fslmaths
 from measureSimilarity import measureSimilarity
-from plumbum import FG
 from skeletonize import skeletonize
 from roi_analysis import roi_analysis
 from antsTemplate import antsMult

lib/tests/data/enigma_res/AD_combined_roi_avg.csv

@@ -1,5 +1,5 @@
 Cases,ACR,ACR-L,ACR-R,ALIC,ALIC-L,ALIC-R,AverageAD,BCC,CGC,CGC-L,CGC-R,CGH,CGH-L,CGH-R,CP,CP-L,CP-R,CST,CST-L,CST-R,EC,EC-L,EC-R,FX,FX_ST,FX_ST-L,FX_ST-R,GCC,ICP,ICP-L,ICP-R,IFO,IFO-L,IFO-R,ML,ML-L,ML-R,PCR,PCR-L,PCR-R,PLIC,PLIC-L,PLIC-R,PTR,PTR-L,PTR-R,RLIC,RLIC-L,RLIC-R,SCC,SCP,SCP-L,SCP-R,SCR,SCR-L,SCR-R,SFO,SFO-L,SFO-R,SLF,SLF-L,SLF-R,SS,SS-L,SS-R,UNC,UNC-L,UNC-R
-dti_tra_dir16_AP_9_134431,0.001176,0.001170,0.001182,0.001250,0.001257,0.001243,0.001184,0.001658,0.001286,0.001327,0.001242,0.000655,0.000728,0.000523,0.001127,0.001211,0.001030,0,0,0,0.001060,0.001128,0.000984,0.002487,0.001167,0.001214,0.001099,0.001572,0,0,0,0.000600,0.000779,0.000128,0,0,0,0.001254,0.001249,0.001258,0.001380,0.001394,0.001367,0.001372,0.001376,0.001368,0.001335,0.001381,0.001291,0.001602,0,0,0,0.001186,0.001215,0.001158,0.001147,0.001163,0.001132,0.001190,0.001210,0.001171,0.001115,0.001276,0.000877,0.001577,0.001535,0.001606
-dti_tra_dir16_PA_6_134431,0.001196,0.001196,0.001196,0.001295,0.001275,0.001314,0.001201,0.001679,0.001321,0.001343,0.001298,0.000615,0.000675,0.000509,0.001176,0.001191,0.001159,0,0,0,0.001051,0.001110,0.000986,0.002797,0.001202,0.001216,0.001183,0.001625,0,0,0,0.000608,0.000791,0.000126,0,0,0,0.001255,0.001242,0.001267,0.001428,0.001390,0.001466,0.001381,0.001377,0.001384,0.001373,0.001361,0.001385,0.001624,0,0,0,0.001212,0.001223,0.001202,0.001210,0.001175,0.001244,0.001216,0.001195,0.001236,0.001116,0.001269,0.000868,0.001566,0.001535,0.001588
-dti_tra_dir16_PA_rot_12_132222,0.001207,0.001193,0.001222,0.001301,0.001265,0.001333,0.001245,0.001681,0.001419,0.001502,0.001330,0.000986,0.000904,0.001078,0.001373,0.001347,0.001397,0,0,0,0.001181,0.001150,0.001212,0.002723,0.001281,0.001248,0.001324,0.001647,0,0,0,0.001373,0.001361,0.001383,0,0,0,0.001259,0.001234,0.001281,0.001473,0.001471,0.001474,0.001428,0.001437,0.001419,0.001367,0.001372,0.001363,0.001664,0,0,0,0.001219,0.001221,0.001218,0.001156,0.001142,0.001169,0.001216,0.001190,0.001242,0.001416,0.001384,0.001449,0.001603,0.001558,0.001634
-dti_tra_dir16_PA_rot_12_134431,0.001192,0.001190,0.001194,0.001288,0.001269,0.001305,0.001194,0.001672,0.001308,0.001342,0.001271,0.000556,0.000612,0.000374,0.001108,0.001211,0.000995,0,0,0,0.001057,0.001096,0.001013,0.002826,0.001154,0.001202,0.001081,0.001605,0,0,0,0.000596,0.000613,0.000024,0,0,0,0.001258,0.001240,0.001274,0.001421,0.001398,0.001443,0.001360,0.001376,0.001345,0.001358,0.001347,0.001368,0.001623,0,0,0,0.001207,0.001216,0.001198,0.001197,0.001155,0.001238,0.001210,0.001187,0.001232,0.001067,0.001170,0.000862,0.001553,0.001539,0.001563
+dti_tra_dir16_AP_9_134431,0.001162,0.001160,0.001165,0.001242,0.001251,0.001234,0.001171,0.001661,0.001164,0.001172,0.001156,0.000659,0.000734,0.000523,0.001208,0.001277,0.001108,0,0,0,0.001063,0.001096,0.001023,0.002445,0.001173,0.001224,0.001093,0.001553,0,0,0,0.000973,0.000973,0,0,0,0,0.001240,0.001233,0.001247,0.001374,0.001392,0.001356,0.001356,0.001360,0.001353,0.001313,0.001357,0.001271,0.001583,0,0,0,0.001176,0.001199,0.001154,0.001076,0.001064,0.001087,0.001158,0.001179,0.001138,0.001241,0.001248,0.001218,0.001460,0.001398,0.001502
+dti_tra_dir16_PA_6_134431,0.001184,0.001185,0.001182,0.001288,0.001271,0.001304,0.001189,0.001678,0.001201,0.001189,0.001214,0.000621,0.000683,0.000511,0.001243,0.001255,0.001227,0,0,0,0.001054,0.001076,0.001026,0.002701,0.001213,0.001201,0.001231,0.001622,0,0,0,0.000993,0.000993,0,0,0,0,0.001244,0.001230,0.001256,0.001417,0.001386,0.001448,0.001387,0.001360,0.001413,0.001352,0.001337,0.001366,0.001595,0,0,0,0.001203,0.001207,0.001199,0.001139,0.001094,0.001182,0.001178,0.001155,0.001201,0.001216,0.001240,0.001074,0.001493,0.001398,0.001559
+dti_tra_dir16_PA_rot_12_132222,0.001188,0.001178,0.001199,0.001291,0.001261,0.001319,0.001210,0.001689,0.001295,0.001349,0.001237,0.000982,0.000898,0.001076,0.001348,0.001338,0.001357,0,0,0,0.001145,0.001121,0.001168,0.002636,0.001255,0.001255,0.001254,0.001619,0,0,0,0.001351,0.001354,0.001348,0,0,0,0.001248,0.001227,0.001267,0.001467,0.001466,0.001468,0.001405,0.001397,0.001413,0.001348,0.001359,0.001338,0.001626,0,0,0,0.001205,0.001203,0.001207,0.001149,0.001132,0.001165,0.001166,0.001146,0.001186,0.001402,0.001372,0.001432,0.001531,0.001510,0.001545
+dti_tra_dir16_PA_rot_12_134431,0.001189,0.001184,0.001194,0.001275,0.001259,0.001290,0.001183,0.001670,0.001190,0.001203,0.001177,0.000560,0.000619,0.000369,0.001233,0.001278,0.001172,0,0,0,0.001054,0.001067,0.001038,0.002722,0.001188,0.001221,0.001122,0.001606,0,0,0,0.000930,0.000930,0,0,0,0,0.001244,0.001230,0.001256,0.001413,0.001393,0.001433,0.001378,0.001360,0.001395,0.001337,0.001317,0.001355,0.001592,0,0,0,0.001196,0.001208,0.001185,0.001161,0.001120,0.001201,0.001170,0.001144,0.001195,0.001182,0.001200,0.000712,0.001525,0.001478,0.001557

lib/tests/data/enigma_res/FA_combined_roi_avg.csv

@@ -1,5 +1,5 @@
 Cases,ACR,ACR-L,ACR-R,ALIC,ALIC-L,ALIC-R,AverageFA,BCC,CGC,CGC-L,CGC-R,CGH,CGH-L,CGH-R,CP,CP-L,CP-R,CST,CST-L,CST-R,EC,EC-L,EC-R,FX,FX_ST,FX_ST-L,FX_ST-R,GCC,ICP,ICP-L,ICP-R,IFO,IFO-L,IFO-R,ML,ML-L,ML-R,PCR,PCR-L,PCR-R,PLIC,PLIC-L,PLIC-R,PTR,PTR-L,PTR-R,RLIC,RLIC-L,RLIC-R,SCC,SCP,SCP-L,SCP-R,SCR,SCR-L,SCR-R,SFO,SFO-L,SFO-R,SLF,SLF-L,SLF-R,SS,SS-L,SS-R,UNC,UNC-L,UNC-R
-dti_tra_dir16_AP_9_134431,0.419876,0.425731,0.413988,0.556162,0.562564,0.550342,0.419650,0.616682,0.512402,0.544969,0.477565,0.260001,0.291410,0.203539,0.504261,0.552324,0.448920,0,0,0,0.389659,0.421413,0.354173,0.419737,0.479476,0.455101,0.515359,0.689340,0,0,0,0.143168,0.183610,0.036286,0,0,0,0.478090,0.469669,0.485592,0.660847,0.653474,0.668146,0.519267,0.523969,0.514710,0.609203,0.615799,0.602983,0.714921,0,0,0,0.489787,0.496939,0.482907,0.479324,0.470568,0.487811,0.473810,0.473450,0.474162,0.388581,0.444175,0.306143,0.570173,0.533605,0.595375
-dti_tra_dir16_PA_6_134431,0.422007,0.425961,0.418031,0.574196,0.576491,0.572110,0.417941,0.618775,0.519960,0.550291,0.487516,0.238570,0.270217,0.182609,0.505573,0.522886,0.485738,0,0,0,0.374179,0.398238,0.347251,0.421627,0.467459,0.468745,0.465620,0.657946,0,0,0,0.139294,0.179535,0.032944,0,0,0,0.471736,0.469447,0.473775,0.662038,0.662409,0.661671,0.518990,0.531467,0.506896,0.617020,0.617719,0.616361,0.727670,0,0,0,0.486069,0.502513,0.470252,0.478499,0.466305,0.490320,0.480248,0.481307,0.479212,0.378176,0.441090,0.275999,0.554327,0.538903,0.564957
-dti_tra_dir16_PA_rot_12_132222,0.414567,0.414075,0.415061,0.556832,0.552736,0.560555,0.438677,0.627485,0.561131,0.614705,0.503823,0.371355,0.354715,0.390156,0.605077,0.619010,0.592075,0,0,0,0.409992,0.408171,0.411806,0.385335,0.505064,0.497541,0.514674,0.696609,0,0,0,0.467126,0.489761,0.448264,0,0,0,0.477984,0.484206,0.472441,0.654561,0.667458,0.641795,0.530505,0.541331,0.520011,0.588465,0.599607,0.577960,0.738051,0,0,0,0.472109,0.481474,0.463101,0.435697,0.427541,0.443604,0.478935,0.468329,0.489320,0.531611,0.522392,0.540859,0.551498,0.568008,0.540120
-dti_tra_dir16_PA_rot_12_134431,0.421583,0.423827,0.419327,0.566532,0.574298,0.559473,0.417577,0.614816,0.517959,0.556796,0.476416,0.211582,0.239383,0.121713,0.457369,0.527764,0.379650,0,0,0,0.366979,0.384583,0.346822,0.412827,0.462908,0.473975,0.446095,0.653552,0,0,0,0.111028,0.114201,0.003142,0,0,0,0.480895,0.480755,0.481020,0.657673,0.664191,0.651220,0.509762,0.526007,0.494016,0.607179,0.614069,0.600682,0.725085,0,0,0,0.480810,0.496148,0.466056,0.473443,0.454241,0.492057,0.478223,0.481382,0.475131,0.353617,0.399767,0.262304,0.562773,0.539770,0.578626
+dti_tra_dir16_AP_9_134431,0.422245,0.426844,0.417613,0.557781,0.562671,0.553337,0.437611,0.619984,0.517415,0.547055,0.485709,0.262096,0.294670,0.203539,0.562475,0.597456,0.512301,0,0,0,0.400146,0.422143,0.373681,0.433788,0.489969,0.464499,0.529791,0.691827,0,0,0,0.256093,0.256093,0,0,0,0,0.483600,0.473751,0.492373,0.661646,0.653863,0.669351,0.530328,0.533059,0.527681,0.608848,0.614758,0.603276,0.716856,0,0,0,0.491489,0.498141,0.485090,0.483796,0.491110,0.476705,0.478385,0.474977,0.481721,0.457305,0.458056,0.454900,0.572257,0.539978,0.594504
+dti_tra_dir16_PA_6_134431,0.423806,0.425870,0.421712,0.576617,0.577226,0.576064,0.436355,0.622903,0.527081,0.549085,0.503544,0.239735,0.272041,0.182609,0.565398,0.587727,0.533848,0,0,0,0.382721,0.395559,0.366933,0.442159,0.481613,0.472675,0.495587,0.668124,0,0,0,0.249813,0.249813,0,0,0,0,0.473331,0.470660,0.475710,0.663306,0.663703,0.662914,0.528718,0.545530,0.512423,0.616793,0.617271,0.616342,0.732847,0,0,0,0.486520,0.501285,0.472318,0.482910,0.450523,0.514305,0.482480,0.482910,0.482059,0.441440,0.451299,0.384415,0.558384,0.548723,0.565043
+dti_tra_dir16_PA_rot_12_132222,0.413580,0.415421,0.411728,0.563548,0.555872,0.570526,0.448449,0.629895,0.574804,0.634598,0.510844,0.369402,0.351851,0.389232,0.594242,0.616545,0.573430,0,0,0,0.410802,0.413133,0.408481,0.391234,0.495661,0.491090,0.501500,0.698609,0,0,0,0.467704,0.491033,0.448264,0,0,0,0.476752,0.484017,0.470281,0.654710,0.668341,0.641216,0.526425,0.532505,0.520531,0.587848,0.599558,0.576808,0.739195,0,0,0,0.464604,0.474205,0.455370,0.450097,0.444213,0.455800,0.476587,0.466908,0.486063,0.532118,0.525666,0.538590,0.549543,0.562751,0.540440
+dti_tra_dir16_PA_rot_12_134431,0.425340,0.425052,0.425633,0.569071,0.577920,0.561027,0.436782,0.618339,0.515120,0.553615,0.473942,0.211370,0.240172,0.118264,0.538410,0.579724,0.481540,0,0,0,0.375081,0.382567,0.365659,0.418468,0.477262,0.482824,0.466337,0.660991,0,0,0,0.231260,0.231260,0,0,0,0,0.486896,0.485122,0.488476,0.658720,0.665745,0.651766,0.520291,0.535497,0.505146,0.606193,0.614361,0.598472,0.728759,0,0,0,0.482329,0.497878,0.467372,0.478855,0.469337,0.488081,0.480986,0.481527,0.480457,0.416426,0.426970,0.145948,0.564925,0.544112,0.579269

lib/tests/data/enigma_res/MD_combined_roi_avg.csv

@@ -1,5 +1,5 @@
 Cases,ACR,ACR-L,ACR-R,ALIC,ALIC-L,ALIC-R,AverageMD,BCC,CGC,CGC-L,CGC-R,CGH,CGH-L,CGH-R,CP,CP-L,CP-R,CST,CST-L,CST-R,EC,EC-L,EC-R,FX,FX_ST,FX_ST-L,FX_ST-R,GCC,ICP,ICP-L,ICP-R,IFO,IFO-L,IFO-R,ML,ML-L,ML-R,PCR,PCR-L,PCR-R,PLIC,PLIC-L,PLIC-R,PTR,PTR-L,PTR-R,RLIC,RLIC-L,RLIC-R,SCC,SCP,SCP-L,SCP-R,SCR,SCR-L,SCR-R,SFO,SFO-L,SFO-R,SLF,SLF-L,SLF-R,SS,SS-L,SS-R,UNC,UNC-L,UNC-R
-dti_tra_dir16_AP_9_134431,0.000796,0.000789,0.000803,0.000724,0.000724,0.000724,0.000789,0.000917,0.000785,0.000786,0.000784,0.000405,0.000443,0.000338,0.000666,0.000692,0.000637,0,0,0,0.000723,0.000757,0.000684,0.001711,0.000741,0.000788,0.000673,0.000811,0,0,0,0.000488,0.000634,0.000102,0,0,0,0.000801,0.000805,0.000798,0.000726,0.000739,0.000714,0.000842,0.000840,0.000844,0.000753,0.000769,0.000737,0.000799,0,0,0,0.000748,0.000760,0.000737,0.000720,0.000732,0.000708,0.000769,0.000781,0.000757,0.000737,0.000830,0.000599,0.000923,0.000927,0.000920
-dti_tra_dir16_PA_6_134431,0.000804,0.000803,0.000806,0.000740,0.000729,0.000750,0.000800,0.000928,0.000799,0.000788,0.000810,0.000388,0.000416,0.000338,0.000696,0.000699,0.000693,0,0,0,0.000729,0.000760,0.000694,0.001906,0.000758,0.000780,0.000726,0.000860,0,0,0,0.000504,0.000655,0.000105,0,0,0,0.000807,0.000800,0.000813,0.000749,0.000730,0.000768,0.000845,0.000834,0.000856,0.000763,0.000754,0.000771,0.000801,0,0,0,0.000766,0.000761,0.000771,0.000757,0.000746,0.000768,0.000779,0.000766,0.000792,0.000745,0.000827,0.000611,0.000929,0.000926,0.000931
-dti_tra_dir16_PA_rot_12_132222,0.000814,0.000807,0.000822,0.000756,0.000740,0.000770,0.000823,0.000916,0.000819,0.000821,0.000817,0.000624,0.000565,0.000690,0.000751,0.000733,0.000768,0,0,0,0.000796,0.000776,0.000815,0.001938,0.000801,0.000787,0.000818,0.000830,0,0,0,0.000875,0.000848,0.000898,0,0,0,0.000805,0.000785,0.000823,0.000781,0.000770,0.000792,0.000859,0.000857,0.000861,0.000787,0.000776,0.000797,0.000804,0,0,0,0.000781,0.000774,0.000787,0.000757,0.000753,0.000761,0.000786,0.000776,0.000796,0.000852,0.000839,0.000865,0.000953,0.000907,0.000985
-dti_tra_dir16_PA_rot_12_134431,0.000800,0.000798,0.000802,0.000742,0.000725,0.000757,0.000795,0.000927,0.000796,0.000784,0.000809,0.000350,0.000380,0.000254,0.000672,0.000692,0.000649,0,0,0,0.000737,0.000756,0.000715,0.001934,0.000727,0.000764,0.000670,0.000850,0,0,0,0.000514,0.000528,0.000021,0,0,0,0.000802,0.000791,0.000812,0.000749,0.000733,0.000764,0.000836,0.000838,0.000834,0.000763,0.000754,0.000772,0.000802,0,0,0,0.000765,0.000760,0.000770,0.000753,0.000739,0.000767,0.000778,0.000763,0.000793,0.000722,0.000775,0.000617,0.000921,0.000940,0.000908
+dti_tra_dir16_AP_9_134431,0.000778,0.000776,0.000781,0.000711,0.000714,0.000709,0.000779,0.000905,0.000769,0.000768,0.000770,0.000405,0.000443,0.000338,0.000673,0.000693,0.000645,0,0,0,0.000732,0.000748,0.000713,0.001651,0.000725,0.000766,0.000660,0.000784,0,0,0,0.000772,0.000772,0,0,0,0,0.000791,0.000792,0.000790,0.000718,0.000732,0.000704,0.000826,0.000829,0.000823,0.000737,0.000752,0.000723,0.000784,0,0,0,0.000741,0.000750,0.000733,0.000699,0.000698,0.000700,0.000756,0.000768,0.000744,0.000825,0.000832,0.000801,0.000896,0.000911,0.000886
+dti_tra_dir16_PA_6_134431,0.000794,0.000792,0.000796,0.000725,0.000717,0.000733,0.000792,0.000913,0.000784,0.000772,0.000797,0.000390,0.000419,0.000338,0.000698,0.000698,0.000697,0,0,0,0.000743,0.000754,0.000730,0.001767,0.000752,0.000757,0.000744,0.000841,0,0,0,0.000805,0.000805,0,0,0,0,0.000793,0.000788,0.000797,0.000741,0.000724,0.000757,0.000854,0.000822,0.000886,0.000751,0.000746,0.000755,0.000779,0,0,0,0.000759,0.000750,0.000767,0.000755,0.000739,0.000771,0.000768,0.000757,0.000779,0.000818,0.000831,0.000746,0.000908,0.000906,0.000909
+dti_tra_dir16_PA_rot_12_132222,0.000800,0.000793,0.000807,0.000734,0.000726,0.000741,0.000798,0.000902,0.000799,0.000805,0.000792,0.000622,0.000562,0.000689,0.000743,0.000723,0.000762,0,0,0,0.000782,0.000764,0.000799,0.001802,0.000789,0.000792,0.000786,0.000799,0,0,0,0.000864,0.000843,0.000882,0,0,0,0.000796,0.000776,0.000814,0.000774,0.000765,0.000783,0.000850,0.000840,0.000860,0.000770,0.000769,0.000771,0.000783,0,0,0,0.000773,0.000766,0.000779,0.000757,0.000746,0.000768,0.000770,0.000762,0.000778,0.000837,0.000822,0.000852,0.000924,0.000896,0.000943
+dti_tra_dir16_PA_rot_12_134431,0.000795,0.000793,0.000798,0.000723,0.000707,0.000737,0.000786,0.000909,0.000782,0.000778,0.000787,0.000352,0.000382,0.000254,0.000711,0.000706,0.000719,0,0,0,0.000746,0.000754,0.000737,0.001787,0.000746,0.000763,0.000714,0.000834,0,0,0,0.000763,0.000763,0,0,0,0,0.000789,0.000779,0.000798,0.000740,0.000722,0.000758,0.000843,0.000820,0.000866,0.000748,0.000738,0.000758,0.000779,0,0,0,0.000756,0.000753,0.000758,0.000749,0.000729,0.000768,0.000767,0.000752,0.000782,0.000805,0.000817,0.000495,0.000903,0.000900,0.000905