From 201711061219 at mail.bnu.edu.cn Mon Feb 8 08:10:14 2021 From: 201711061219 at mail.bnu.edu.cn (ZhengxinGong) Date: Mon, 8 Feb 2021 21:10:14 +0800 (GMT+08:00) Subject: [Neuroimaging] =?utf-8?q?=5Bfmriprep=5D=5Bciftify=5D_WHY_are_the?= =?utf-8?q?re_significant_differences_between_data_of_dtseries=2Enii_files?= =?utf-8?q?_derivated_from_two_piplines_respectively=3F?= Message-ID: Dear experts, We compared the data in the files ('~.dtseries.nii') proccessed respectively from the fmriprep pipeline and the ciftify pipeline, and found that the differences is intolerable. However, we had difficulties to further explore the causes. Following is the important details of our investigation. 1. We first preprocessed data by fmriprep20.2.1: [ codes in pdf ] Because we set on '--cifti-otput', fmriprep will produced .dtseries.nii file for each run, like '~space-fsLR_den-91k_bold.dtseries.nii'. 2. For some reasons in my project, we use ciftify to transfer outputs of fmriprep into HCP format: [ codes in pdf ] The second command takes input of the '~space-T1w_desc-prepro_bold.nii.gz' in fmriprep outputs, and produce the '~Atlas_s0.dtseries.nii'. 3. We then compare the time-series data in two .dtseries.nii. : 3.1 For each position(cortical vertex or subcortical voxel) in brain image, compute correlation of time series from two files. [ codes in pdf ] The mean of r is around 0.34 for a randomly selected scan, and we repeat the similar result on another scan. Combined with the histogram of r, indicating that there are even thousands of negative values, we were very astonished and curious about the cause of this probelm. 3.2 Visulization We tried to visulize two dtseries and the r-map, but find little cue on the cause. r map on a inflated template visualization at TR=0, two time-series data flash on the inflated template 4. Some additional analysis We thought that one plausible reason may lies on the the usage of '--use-aroma' in fmriprep. For this possibility, we check the log to determine how fmriprep generate .dtseries.nii, but the log information is unclear about it. However, we rerun fmriprep without aroma, and the comparison remain the same conclusion. We are planning to use traditional analysis like funtional localizer to continue the comparison as a validation of outputs from two pipelines. If anyone has already known some hint on this problem, please contact us! Thank you for reading. Best regards -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: dtseries_diff.gif Type: image/gif Size: 505485 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: r_map.png Type: image/png Size: 740236 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: hist.png Type: image/png Size: 15867 bytes Desc: not available URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: [fmriprep][ciftify] WHY are there significant differences between data.pdf Type: application/pdf Size: 1059217 bytes Desc: not available URL: