biac:analysis:physiological
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionLast revisionBoth sides next revision | ||
biac:analysis:physiological [2011/10/25 19:37] – [Physiological Correction Methods] petty | biac:analysis:physiological [2022/01/20 17:17] – 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 |
- | There' | + | Stock GE physio is hardcoded by GE to be 25Hz for respiration and 100Hz for pulse. |
- | 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 | + | [cmp12@blade17 20211218_03399]$ / |
+ | usage: convert_GEphysio.py [-h] --xmlh XMLFNAME [--json JSONFNAME] --physdir | ||
+ | | ||
+ | | ||
+ | |||
+ | create 100Hz resampled TSV from raw GE physio data | ||
+ | |||
+ | optional arguments: | ||
+ | -h, --help | ||
+ | --xmlh XMLFNAME | ||
+ | --json JSONFNAME | ||
+ | --physdir PHYSDIR | ||
+ | --bidsout | ||
+ | --afnislicetime | ||
+ | --delay DELAYTIME | ||
+ | | ||
</ | </ | ||
+ | |||
+ | Provide a valid BXH ( and JSON if possible ) and a path to your subject physio folder to convert_GEphysio.py function. | ||
+ | |||
+ | This function will search the physio directory for physio files within a specified time range based on the functional scantime + timepoints + a default delay in writing. | ||
+ | |||
+ | 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 to 100Hz and locked to the onset of the scan. If disdaqs are detected in the BXH the corresponding physio time will be discarded. | ||
< | < | ||
- | physio_create.py --bxh /path/to/run4.bxh --pdigm /path/to/ | + | / |
+ | </code> | ||
- | Program to produce formatted physiological data: | + | By default you will get a 100Hz resampled TSV with headers |
- | Data will be corrected | + | |
- | BXH and PDIGM required | + | |
- | OUTPATH defaults to PWD | + | |
- | FORMAT: defaults to fsl | + | |
- | - fsl = cardio, | + | |
- | - npm = time, | + | |
- | - afni = cardiac.txt, | + | |
- | Options: | + | The **--bidsout** flag will output a TSV without headers |
- | | + | |
- | -b FILE, --bxh=FILE | + | |
- | -p FILE, --pdigm=FILE | + | |
- | pdigm file for run | + | |
- | -f string, --format=string | + | |
- | output format ( fsl, | + | |
- | -o PATH, --outpath=PATH | + | |
- | location to store output files | + | |
- | </ | ||
- | ===== AFNI ===== | + | The **--afnislicetime** flag will output two 1D slicetiming files to be used with RetroTS.py and afni_proc.py directly. |
- | Create the AFNI formated physiological text data. | + | === Using afni === |
- | <code bash> | + | |
- | This would create | + | You can create |
- | Now set-up build the regressors that will be used in afni with their matlab tool. | + | < |
+ | python3 / | ||
- | <code matlab> | + | -n #slices |
- | addpath / | + | -v TR in seconds |
- | Opt.Respfile = '/ | + | |
- | Opt.Cardfile = '/ | + | |
- | Opt.VolTR = 1.5 | + | |
- | Opt.Nslices = 34 | + | |
- | Opt.PhysFS = 100 | + | |
- | RetroTS(Opt) | ||
</ | </ | ||
- | The output | + | 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: | + | === setup afni despike, retroicor and slicetime correction === |
- | <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 with all the afni data inside. | + | < |
- | < | + | python3 / |
- | # | + | |
- | 3dAFNItoNIFTI 12345.results/pb02.12345.r01.ricor+orig.BRIK | + | |
- | fslwrapbxh pb02.12345.r01.ricor.BRIK | + | -dsets ../../../ |
+ | -tcat_remove_first_trs 2 \ | ||
+ | | ||
+ | | ||
+ | | ||
+ | -verb 1 | ||
</ | </ | ||
- | Run the resting state pipeline: | + | The above will produce an afni script to throw out 2 timepoints, perform a despike, retroicor physio regression and slice time correction based on the RetroTS output and the physio TSV/json above. |
- | <code bash> | + | |
+ | Convert the afni corrected output back to nii.gz to be used for other pre-processing if not using AFNI for all steps | ||
- | ===== FSL Feat ===== | + | < |
+ | bxh2analyze --niigz -s pb03.sub-00104_run-4_ricor_phys.r01.tshift+orig.HEAD sub-00104_run-4_physiocorrected | ||
+ | </ | ||
+ | |||
- | 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 processing ====== | ||
+ | [[biac: | ||
biac/analysis/physiological.txt · Last modified: 2024/06/21 15:44 by 127.0.0.1