aded support for custom IRF in FSL models

Author: TheChymera

nipype/interfaces/fsl/model.py

@@ -42,6 +42,7 @@ class Level1DesignInputSpec(BaseInterfaceInputSpec):
             'dgamma'), traits.Dict(traits.Enum('derivs'), traits.Bool)),
         traits.Dict(traits.Enum('gamma'), traits.Dict(
                     traits.Enum('derivs', 'gammasigma', 'gammadelay'))),
+        traits.Dict(traits.Enum('custom'), traits.Dict(traits.Dict(traits.Enum('bfcustompath'), traits.Str))),
         traits.Dict(traits.Enum('none'), traits.Dict()),
         traits.Dict(traits.Enum('none'), traits.Enum(None)),
         mandatory=True,
@@ -193,10 +194,11 @@ def _create_ev_files(
                             evinfo.insert(j, [onset, cond['duration'][j], amp])
                         else:
                             evinfo.insert(j, [onset, cond['duration'][0], amp])
+                    ev_parameters['cond_file'] = evfname
                     ev_parameters['ev_num'] = num_evs[0]
                     ev_parameters['ev_name'] = name
+                    ev_parameters['fsldir'] = os.environ['FSLDIR']
                     ev_parameters['tempfilt_yn'] = do_tempfilter
-                    ev_parameters['cond_file'] = evfname
                     try:
                         ev_parameters['temporalderiv'] = int(bool(ev_parameters.pop('derivs')))
                     except KeyError:

nipype/interfaces/script_templates/feat_ev_custom.tcl

@@ -0,0 +1,38 @@
+# EV title
+set fmri(evtitle$ev_num) "$ev_name"
+
+# Basic waveform shape (EV $ev_num)
+# 0 : Square
+# 1 : Sinusoid
+# 2 : Custom (1 entry per volume)
+# 3 : Custom (3 column format)
+# 4 : Interaction
+# 10 : Empty (all zeros)
+set fmri(shape$ev_num) 3
+
+# Convolution (EV $ev_num)
+# 0 : None
+# 1 : Gaussian
+# 2 : Gamma
+# 3 : Double-Gamma HRF
+# 4 : Gamma basis functions
+# 5 : Sine basis functions
+# 6 : FIR basis functions
+# 7 : Optimal/custom basis functions
+set fmri(convolve$ev_num) 7
+
+# Convolve phase (EV $ev_num)
+set fmri(convolve_phase$ev_num) 0
+
+# Apply temporal filtering (EV $ev_num)
+set fmri(tempfilt_yn$ev_num) $tempfilt_yn
+
+# Add temporal derivative (EV $ev_num)
+set fmri(deriv_yn$ev_num) $temporalderiv
+
+# Custom EV file (EV $ev_num)
+set fmri(custom$ev_num) "$cond_file"
+
+# Optimal/custom HRF convolution file
+set fmri(default_bfcustom) "${fsldir}etc/default_flobs.flobs/hrfbasisfns.txt"
+set fmri(bfcustom) "$bfcustompath"

nipype/interfaces/script_templates/feat_ev_gamma.tcl

@@ -18,6 +18,7 @@ set fmri(shape$ev_num) 3
 # 4 : Gamma basis functions
 # 5 : Sine basis functions
 # 6 : FIR basis functions
+# 7 : Optimal/custom basis functions
 set fmri(convolve$ev_num) 2
 
 # Convolve phase

nipype/interfaces/script_templates/feat_ev_hrf.tcl

@@ -18,6 +18,7 @@ set fmri(shape$ev_num) 3
 # 4 : Gamma basis functions
 # 5 : Sine basis functions
 # 6 : FIR basis functions
+# 7 : Optimal/custom basis functions
 set fmri(convolve$ev_num) 3
 
 # Convolve phase (EV $ev_num)

nipype/interfaces/script_templates/feat_ev_none.tcl

@@ -18,6 +18,7 @@ set fmri(shape$ev_num) 2
 # 4 : Gamma basis functions
 # 5 : Sine basis functions
 # 6 : FIR basis functions
+# 7 : Optimal/custom basis functions
 set fmri(convolve$ev_num) 0
 
 # Apply temporal filtering