biac:analysis:resting_pipeline
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionNext revisionBoth sides next revision | ||
biac:analysis:resting_pipeline [2012/07/09 19:02] – [Step 4] admin | biac:analysis:resting_pipeline [2013/10/02 15:29] – cmp12 | ||
---|---|---|---|
Line 17: | Line 17: | ||
4 - normalize data | 4 - normalize data | ||
5 - regress out WM/CSF | 5 - regress out WM/CSF | ||
- | 6 - lowpass | + | 6 - bandpass |
7 - do parcellation and produce correlation matrix from label file | 7 - do parcellation and produce correlation matrix from label file | ||
* or split it up: | * or split it up: | ||
Line 25: | Line 25: | ||
--corrts overrides default location for input parcellation | --corrts overrides default location for input parcellation | ||
results (outputpath/ | results (outputpath/ | ||
+ | 8 - functional connectivity density mapping | ||
+ | |||
Line 33: | Line 35: | ||
--throwaway=4 | --throwaway=4 | ||
run | run | ||
+ | --t1=/ | ||
-p func, --prefix=func | -p func, --prefix=func | ||
- | prefix for all resulting images, defaults to name of | + | prefix for all resulting images, defaults to name of |
input | input | ||
- | --t1=/ | ||
-s 0,1,2,3, --steps=0, | -s 0,1,2,3, --steps=0, | ||
comma seperated string of steps. ' | comma seperated string of steps. ' | ||
Line 54: | Line 56: | ||
standard brain | standard brain | ||
--refcsf=FILE | --refcsf=FILE | ||
+ | standard brain | ||
+ | --refgm=FILE | ||
standard brain | standard brain | ||
--refbrainmask=FILE | --refbrainmask=FILE | ||
Line 61: | Line 65: | ||
brain | brain | ||
--betfval=0.4 | --betfval=0.4 | ||
- | --anatbetfval=0.5 | + | --anatbetfval=0.5 |
+ | | ||
--lpfreq=0.08 | --lpfreq=0.08 | ||
- | is .08hz | + | is .08hz. highpass is fixed at .001hz. |
--corrlabel=FILE | --corrlabel=FILE | ||
correlation search. default is the 116 region AAL | correlation search. default is the 116 region AAL | ||
Line 92: | Line 97: | ||
If --dvarsthreshold is specified, then | If --dvarsthreshold is specified, then | ||
--dvarsnumnumneighbors specifies how many neighboring | --dvarsnumnumneighbors specifies how many neighboring | ||
+ | volumes, before and after the initially excluded | ||
+ | volumes, should also be excluded. | ||
+ | --fdthreshold=THRESH | ||
+ | Any volume contributing to a FD value greater than | ||
+ | this threshold will be excluded (" | ||
+ | (final) correlation step. FD calculation is performed | ||
+ | on the results of the last pre-processing step, and is | ||
+ | calculated as described by Power, J.D., et al., | ||
+ | " | ||
+ | connectivity MRI networks arise from subject motion", | ||
+ | NeuroImage(2011). | ||
+ | the final correlation, | ||
+ | operations that require the full signal, like | ||
+ | regression, etc. | ||
+ | --fdnumneighbors=NUMNEIGHBORS | ||
+ | If --fdthreshold is specified, then | ||
+ | --fdnumnumneighbors specifies how many neighboring | ||
volumes, before and after the initially excluded | volumes, before and after the initially excluded | ||
volumes, should also be excluded. | volumes, should also be excluded. | ||
Line 119: | Line 141: | ||
perform motion correction, then this option is | perform motion correction, then this option is | ||
ignored. | ignored. | ||
+ | --scrubop=SCRUBOP | ||
+ | --fdthreshold are specified, then --scrubop specifies | ||
+ | the aggregation operator used to determine the final | ||
+ | list of excluded volumes. | ||
+ | means a volume will be excluded if *any* of its | ||
+ | thresholds are exceeded, whereas ' | ||
+ | thresholds must be exceeded to be excluded. | ||
+ | --powerscrub | ||
+ | --fdnumneighbors=0 --dvarsthreshold=0.5% | ||
+ | --dvarsnumneigbhors=0 --scrubop=' | ||
+ | method used in the Power et al. article. | ||
+ | conflicting options specified before or after this | ||
+ | will override these. | ||
--scrubkeepminvols=NUMVOLS | --scrubkeepminvols=NUMVOLS | ||
- | If --motionthreshold | + | If --motionthreshold, --dvarsthreshold, or |
- | specified, then --scrubminvols specifies the minimum | + | |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
+ | --fcdmthresh=THRESH | ||
+ | connectivity density mapping ( step8 ). Default is set | ||
+ | to 0.6. Algorithm from Tomasi et al, PNAS(2010), vol. | ||
+ | 107, no. 21. Calculates the fcdm of functional data | ||
+ | from last completed step, inside a dilated gray matter | ||
+ | mask | ||
--cleanup | --cleanup | ||
</ | </ | ||
Line 175: | Line 216: | ||
==== Step 6 ==== | ==== Step 6 ==== | ||
* this step will band-pass filter data to remove high-frequency noise using custom python code | * this step will band-pass filter data to remove high-frequency noise using custom python code | ||
- | * the default is 0.08 HZ | + | * the default |
+ | * highpass is fixed at .001 HZ | ||
* if you'd like to chose a different frequency, please use ** --lpfreq ** | * if you'd like to chose a different frequency, please use ** --lpfreq ** | ||
==== Step 7 ==== | ==== Step 7 ==== | ||
Line 202: | Line 244: | ||
* that is the zr_values from the subject.graphml displayed in matplotlib. | * that is the zr_values from the subject.graphml displayed in matplotlib. | ||
- | For a beta version of the viewer: | + | For the rs-pipeline |
< | < | ||
- | / | + | / |
Usage: | Usage: | ||
- | beta_viewer.py --graphml / | + | rspipe_viewer.py --graphml / |
Program to display graphml output from resting_pipeline | Program to display graphml output from resting_pipeline | ||
Line 223: | Line 265: | ||
</ | </ | ||
+ | |||
+ | ==== Step 8 ==== | ||
+ | * Functional connectivity density mapping | ||
+ | * Takes functional data from last step and calculates how connected they are to the voxels around them | ||
+ | * uses ( --fcdmthresh and --refgm ) as the pearson r-value and gray matter mask | ||
+ | * if defaults are used, then a dilated gray matter mask is used from FAST segmentation of MNI brain and a pearson r value of 0.6 | ||
+ | * Iteratively goes to all neighboring voxels and counts the number that have correlated signal until they are under the r threshold | ||
+ | * adapted from Dardo Tomasi, PNAS(2010), vol. 107, no. 21. 9885–9890 | ||
+ | * resulting file with be " | ||
+ | |||
+ | {{: | ||
===== Things to consider ===== | ===== Things to consider ===== | ||
* this was designed to be modular, so that you only need to run the steps you need | * this was designed to be modular, so that you only need to run the steps you need | ||
Line 362: | Line 415: | ||
- | Still under-development | ||
3D VTK: | 3D VTK: | ||
biac/analysis/resting_pipeline.txt · Last modified: 2024/06/21 15:44 by 127.0.0.1