biac:analysis:physiological
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:physiological [2011/12/08 20:46] – [AFNI details] petty | biac:analysis:physiological [2022/01/20 14:46] – cmp12 | ||
---|---|---|---|
Line 1: | Line 1: | ||
====== Physiological Correction Methods ====== | ====== Physiological Correction Methods ====== | ||
- | Currently this assumes that you have recorded your physiological responses through Cigal. | + | ===== Generate resampled TSV from raw GE physio files automatically copied from the scanners |
- | ===== physio_run.py ===== | + | Stock GE physio |
- | There is currently a tool for running AFNI-based correction of functional image data using physiological data. | + | |
- | < | ||
- | physio_run.py --help | ||
- | </ | ||
< | < | ||
- | Usage: | + | [cmp12@blade17 20211218_03399]$ |
- | physio_run.py | + | usage: convert_GEphysio.py |
+ | | ||
+ | | ||
- | Program to produce formatted physiological | + | create 100Hz resampled TSV from raw GE physio |
- | Data will be corrected | + | |
- | BXH and PDIGM required | + | optional arguments: |
- | | + | -h, --help |
- | | + | |
+ | | ||
+ | --physdir PHYSDIR | ||
+ | | ||
+ | | ||
+ | --delay DELAYTIME | ||
+ | | ||
- | Options: | ||
- | -h, --help | ||
- | -f string, --format=string | ||
- | output format ( afni ) | ||
</ | </ | ||
- | '' | + | Provide a valid BXH ( and JSON if possible ) and a path to your subject physio folder to convert_GEphysio.py function. |
- | The outputs | + | This function |
- | ===== physio_create.py ===== | + | If matching respiration and pulse files are found theoretical scan triggers on each TR will be created and the respiration and pulse data will be resampled |
- | + | ||
- | The underlying tool used by '' | + | |
- | + | ||
- | This tool uses various fields from the bxh and pdigm file to output " | + | |
- | + | ||
- | The python tool will call cigal functions directly to do the corrections, | + | |
< | < | ||
- | physio_create.py --help | + | / |
</ | </ | ||
- | < | + | By default you will get a 100Hz resampled TSV with headers and timing. |
- | physio_create.py --bxh / | + | |
- | Program to produce formatted physiological data: | ||
- | Data will be corrected and recreated from the pdigm file based on information from run data | ||
- | BXH and PDIGM required | ||
- | OUTPATH defaults to PWD | ||
- | HZ sampling rate | ||
- | FORMAT: defaults to fsl | ||
- | - fsl = cardio, | ||
- | - npm = time, | ||
- | - afni = cardiac.txt, | ||
+ | The **--bidsout** flag will output a TSV without headers and a BIDS sidecar with appropriate column labels and sampling frequency. | ||
- | Options: | ||
- | -h, --help | ||
- | -b FILE, --bxh=FILE | ||
- | -p FILE, --pdigm=FILE | ||
- | pdigm file for run | ||
- | -f string, --format=string | ||
- | output format ( fsl, | ||
- | --hz sampling rate is hz ( 100 ) | ||
- | -o PATH, --outpath=PATH | ||
- | location to store output files | ||
- | </ | + | The **--afnislicetime** flag will output two 1D slicetiming files to be used with RetroTS.py and afni_proc.py directly. |
- | ===== AFNI details ===== | + | === Using afni === |
- | This section details the steps that the tool '' | + | You can create regressors for afni preprocessing from output files above using RetroTS.py |
- | Create the AFNI formated physiological text data. | + | < |
- | < | + | python3 |
- | This would create my " | + | -n #slices |
+ | -v TR in seconds | ||
- | <code matlab> | ||
- | addpath / | ||
- | Opt.Respfile = '/ | ||
- | Opt.Cardfile = '/ | ||
- | Opt.VolTR = 1.5 | ||
- | Opt.Nslices = 34 | ||
- | Opt.PhysFS = 100 | ||
- | Opt.SliceOrder = ' | ||
- | %please choose the correct order if not running interleaved: | ||
- | |||
- | RetroTS(Opt) | ||
</ | </ | ||
- | The output of runnning the MATLAB script | + | RetroTS |
- | + | ||
- | Convert your data to 4D nifti if you haven' | + | |
- | <code bash> | + | |
- | + | ||
- | Create the afni script which you will use to actually run the 3dretroicor functionality: | + | |
- | <code bash> | + | |
- | ss_method ' | + | |
- | + | ||
- | The result of the above command would be a tcsh script, which would run the despiking and 3dretroicor correction only. Run it: | + | |
- | <code bash> | + | |
- | + | ||
- | There will be a resulting folder | + | |
- | <code bash> | + | |
- | # | + | |
- | bxh2analyze --niigz -s 12345.results/ | + | |
- | </ | + | |
- | + | ||
- | Run the resting state pipeline: | + | |
- | <code bash> | + | |
- | ===== FSL Feat ===== | + | |
- | We've gotten the FSL beta functions for physiological correction, which we'll test. | ||
- | They are based on the same algorithms used in afni with retroicor. | ||
- | More information can be seen here: | ||
- | [[http:// | ||
- | The beta options will allow you to plug in the resulting regressors directly into Feat. After pre-processing you can use the **res4d.nii.gz** files for further analysis ( filtered_func_data would not have the physio regressor removed ). | ||
+ | ====== CIGAL pdigm specific ====== | ||
+ | [[biac: | ||
biac/analysis/physiological.txt · Last modified: 2024/06/21 15:44 by 127.0.0.1