DPABISurf: A Surface-Based Resting-State fMRI Data Analysis Toolbox

Submitted by YAN Chao-Gan on

The toolbox for Data Processing & Analysis for Brain Imaging on Surface (DPABISurf) (Yan et al., 2021) is based on fMRIPrep (Esteban et al., 2019), FreeSurfer (Dale et al., 1999), ANTs (Avants et al., 2008), FSL (Jenkinson et al., 2002), AFNI (Cox, 1996), SPM (Ashburner, 2012), dcm2niix (Li et al., 2016), PALM (Winkler et al., 2016), GNU Parallel (Tange, 2011), MATLAB (The MathWorks Inc., Natick, MA, US), Docker (https://docker.com) and DPABI (Yan et al., 2016). DPABISurf provides user-friendly graphical user interface (GUI) for pipeline surface-based preprocessing, statistical analyses and results viewing, while requires no programming/scripting skills from the users.

The DPABISurf pipeline first converts the user specified data into BIDS format (Gorgolewski et al., 2016), and then calls fMRIPprep docker to preprocess the structural and functional MRI data, which integrates FreeSurfer, ANTs, FSL and AFNI. With fMRIPprep, the data is processed into FreeSurfer fsaverage5 surface space and MNI volume space. DPABISurf further performs nuisance covariates regression (including ICA-AROMA) on the surface-based data (volume-based data is processed as well), and then calculate the commonly used R-fMRI metrics: amplitude of low frequency fluctuation (ALFF) (Zang et al., 2007), fractional ALFF (Zou et al., 2008), regional homogeneity (Zang et al., 2004), degree centrality (Zuo and Xing, 2014), and seed-based functional connectivity. DPABISurf also performs surface-based smoothing by calling FreeSurfer’s mri_surf2surf command. These processed metrics then enters surfaced-based statistical analyses within DPABISurf, which could perform surfaced-based permutation test with TFCE by integrating PALM. Finally, the corrected results could be viewed by the convenient surface viewer DPABISurf_VIEW, which is derived from spm_mesh_render.m.

DPABISurf is designed to make surface-based data analysis require minimum manual operations and almost no programming/scripting experience. We anticipate this open-source toolbox will assist novices and expert users alike and continue to support advancing R-fMRI methodology and its application to clinical translational studies.

DPABISurf is open-source and distributed under GNU/GPL, available with DPABI at http://www.rfmri.org/dpabi. It supports Windows 10 Pro, MacOS and Linux operating systems. You can run it with or without MATLAB.

1. With MATLAB.   
1.1. Please go to http://www.rfmri.org/dpabi to download DPABI.    
1.2. Add with subfolders for DPABI in MATLAB's path setting.   
1.3. Input 'dpabi' and then follow the instructions of the "Install" Button on DPABISurf.   
2. Without MATLAB.   
2.1. Install Docker.   
2.2. Terminal: docker pull cgyan/dpabi   
2.3. Terminal: docker run -d --rm -v /My/FreeSurferLicense/Path/license.txt:/opt/freesurfer/license.txt -v /My/Data/Path:/data -p 5925:5925 cgyan/dpabi x11vnc -forever -shared -usepw -create -rfbport 5925   
/My/FreeSurferLicense/Path/license.txt: Where you stored the FreeSurferLicense got from https://surfer.nmr.mgh.harvard.edu/registration.html.   
/My/Data/Path: This is where you stored your data. In Docker, the path is /data.   
2.4. Open VNC Viewer, connect to localhost:5925, the password is 'dpabi'.   
2.5. In the terminal within the VNC Viewer, input "bash", and then input:   
/opt/DPABI/DPABI_StandAlone/run_DPABI_StandAlone.sh ${MCRPath}


Now please enjoy the StandAlone version of DPABISurf with GUI!


If you don't want to run with GUI, you can also call the compiled version of DPABISurf_run. E.g.,   
docker run -it --rm -v /My/FreeSurferLicense/Path/license.txt:/opt/freesurfer/license.txt -v /My/Data/Path:/data cgyan/dpabi /bin/bash   
/opt/DPABI/DPABI_StandAlone/run_DPABISurf_run_StandAlone.sh ${MCRPath} /data/DPABISurf_Cfg.mat

 
New features of DPABISurf_V3.0_231111 within DPABI_V8.0_231111 (download at http://rfmri.org/dpabi).

1. DPABISurfSlurm: High Performance Computing Version of DPABISurf was released. Please read here for how to use it: http://rfmri.org/DPABISurfSlurm  
2. Back engine updated to fMRIPrep V23.1.4, thus using FreeSurfer 7.3.2.  
3. Segment subregions (Thalamus, Hippocampus, Amygdala and Brainstem) with FreeSurfer 7.3.2 was set to the default option now.  
4. Added freesurfer atlases: {DPABISurf}/SurfTemplates/fsaverage5_lh_aparc_a2009s_annot.label.gii and fsaverage5_lh_aparc_annot.label.gii.  
5. Could make use of DPABI Harmonization module.

 
Tips: 
1) For Linux or Mac OS, please start matlab from terminal in order to reach docker in DPABI (e.g., Linux: matlab; Mac: open /Applications/MATLAB_R2018a.app/).
2) Before running DPABISurf_Pipeline, you can test the docker environment by running DPABI->DPABISurf->Utilities->Volume-Surface Projector. If the file can be successfully projected to surface, then the software is all set. 
 
New features of DPABISurf_V1.8_220915 within DPABI_V6.2_220915 (download at http://rfmri.org/dpabi).
1. Fixed a bug in displaying negative values in DPABISurf_VIEW.
2. Added area and curv in DPABISurf templates.
 
 
New features of DPABISurf_V1.7_220101 within DPABI_V6.1_220101 (download at http://rfmri.org/dpabi).
1. fMRIPrep backend updated to version 20.2.5 (Long-Term Support).
2. If the TR info in NIfTI header is different from what specified by the user in DPABISurf_Pipeline GUI, then rewrite as user specified.
 
 
New features of DPABISurf_V1.6_210501 within DPABI_V6.0_210501 (download at http://rfmri.org/dpabi).
1. Can use individual ROI definition files for each subject. The first line of the ROI definition file should be: 'Seed_ROI_List:'. An example SeedROI.txt:
Seed_ROI_List:
/Processing/fmriprep/sub-Sub001/func/sub-Sub001_task-rest_space-MNI152NLin2009cAsym_desc-aseg_dseg.nii.gz
/Processing/fmriprep/sub-Sub002/func/sub-Sub002_task-rest_space-MNI152NLin2009cAsym_desc-aseg_dseg.nii.gz
/Processing/fmriprep/sub-Sub003/func/sub-Sub003_task-rest_space-MNI152NLin2009cAsym_desc-aseg_dseg.nii.gz
...
2. Fixed a bug in dealing with multiple sessions of FieldMap.
3. Allow skipping subjects in TRInfo.tsv.
4. The surface for ReHo calculation were changed from pial to white.
 
 
New features of DPABISurf_V1.5_201201 within DPABI_V5.1_201201 (download at http://rfmri.org/dpabi, please also update the docker file by: docker pull cgyan/dpabi):
1. fMRIPrep backend updated to version 20.2.1 (Long-Term Support).
2. Fixed a bug during processing data with multiple sessions.
3. DPABISurf_VIEW. Updated Cluster Report function. Please see here for usage.
4. DPABISurf_VIEW. Fixed a bug when setting p = 0.001, but still display 1. Fixed a bug to input index. Added save colorbar.
 
New features of DPABISurf_V1.3_200401 within DPABI_V4.3_200401 (download at http://rfmri.org/dpabi, please also update the docker file by: docker pull cgyan/dpabi):
Fixed the errors of "--ignore-aroma-denoising-errors" and "--template-resampling-grid". These were caused by the input format changing in fmriprep V20.
 
 
New features of DPABISurf_V1.4_201001 within DPABI_V5.0_201001 (download at http://rfmri.org/dpabi, please also update the docker file by: docker pull cgyan/dpabi):
1. fMRIPrep backend updated to version 20.2.0 (Long-Term Support).
2. For surface based multiple comparison correction, the method of FDR correction was added to DPABISurf. See DPABISurf_VIEW->Cluster...->Apply FDR Correction.
3. For surface based multiple comparison correction, the method of Monte Carlo simulation was added to DPABISurf. See DPABISurf_VIEW->Cluster...->Apply FWE (Monte Carlo Simulation) Correction.
4. Pre-calculated Monte Carlo Simulation tables were added for most used situations. Please see {DPABI}/StatisticalAnalysis/DPABISurf_MonteCarlo/MonteCarloTable/. Users can index the minimum cluster area for an estimated smoothness and set at DPABISurf_VIEW->Cluster...->Set Cluster Size.
5. Docker File: MCR files extracted for easy use in Singularity.
 
New features of DPABISurf_V1.3_200301 within DPABI_V4.3_200301 (download at http://rfmri.org/dpabi, please also update the docker file by: docker pull cgyan/dpabi):
1. Stability Analysis module was added. You can calculate volume-based and surface-based stability from DPABI->Dynamic & Stability Analyses. You can read our recent work for the details of stability measure: Li, L., Lu, B., Yan, C.G. (2019). Stability of dynamic functional architecture differs between brain networks and states. Neuroimage, 116230, doi:10.1016/j.neuroimage.2019.116230.
2. Field map correction was added, both for DPARSF and DPABISurf. If you want to perform field map Correction, you need to arrange each subject's field map DICOM files in one directory, and then put them in "FieldMap" directory under the working directory. i.e.:';...
    '{Working Directory}\FieldMap\PhaseDiffRaw\Subject001\xxxxx001.dcm';...
    '{Working Directory}\FieldMap\PhaseDiffRaw\Subject001\xxxxx002.dcm';...
    '...';...
    '{Working Directory}\FieldMap\PhaseDiffRaw\Subject002\xxxxx001.dcm';...
    '{Working Directory}\FieldMap\PhaseDiffRaw\Subject002\xxxxx002.dcm';...
    '...';...
    '{Working Directory}\FieldMap\Magnitude1Raw\Subject001\xxxxx001.dcm';...
    '{Working Directory}\FieldMap\Magnitude1Raw\Subject001\xxxxx002.dcm';...
    '...';...
    '{Working Directory}\FieldMap\Magnitude1Raw\Subject002\xxxxx001.dcm';...
    '{Working Directory}\FieldMap\Magnitude1Raw\Subject002\xxxxx002.dcm';...
    '...';...
    '...';...
Then you can click the button of “FieldMap” button to set field map correction parameters. In most cases, you can use the default  “0” value to let the program read the parameters (e.g., echo times) from the DICOM files.
3. Check data organization function added. For the new users of DPARSF and DPABISurf, most of the errors were caused by data organization! Please use DPABI->Utilities->Check Data Organization to check your data organization before running DPARSF or DPABISurf. This program will lead you to organize your data correctly with prompting messages!
4. Slice timing information read from DICOM files. If you are starting with DICOM files, you no longer need to set the slice timing correction parameters. Just leave it as default (slice number: 0), then DPARSF or DPABISurf will read the parameters from DICOM files. This new feature thanks to Dr. Chris Rorden's new dcm2niiX program (version v1.0.20190902).
5. DPARSF V5.0 now is compatible with BIDS format. You can start with BIDS format data by checking checkbox “BIDS to DPARSF” and setting  “Starting Directory Name” to “BIDS”.
6. DPABISurf V1.3. Check and re-run fmriprep failed subjects. If for any reason, the program failed fmriprep running in DPABISurf, you just need to re-run starting with the step “Preprocessing with fmriprep” and set the “Starting Directory Name” to “BIDS” in DPABISurf_Pipeline. Alternatively, you can run a single checking step from DPABI->DPABISurf->Utilities->Re-Run fmriprep Failed Subjects.
7. DPABI->Utilities->DICOM Sorter. Now DICOM Sorter will remove illegal characters for file names in DICOM sorter. E.g., a subject id of “Wang’#$#’s” will no longer cause a problem in sorting DICOM files.
8. y_Call_DPABISurf_VIEW_FromVolume.m and y_Call_DPABISurf_VIEW.m were added. Now you can use script to call DPABISurf_VIEW to generate surface maps (with higher quality now) in batch. Remember to add “close all” in the for loop to prevent too many windows.
 
New features of DPABISurf_V1.2_190919 within DPABI_V4.2_190919 (download at http://rfmri.org/dpabi, please also update the docker file by: docker pull cgyan/dpabi):
1. DPABISurf_V1.2_190919 updated.
1.1. A quality control module was added to DPABISurf. Now users can quality control surface reconstruction, EPI to T1 registration and T1 to MNI registration for all the subjects in one HTML file, respectively (based on fmriprep 1.5.0). For volume-based analysis, users can also generate group mask for DPABISurf, and exclude subjects by thresholding coverage and head motion.
1.2. DPABISurf now also output sulcus depth and volume in fsaverage and fsaverage5 spaces for statistical analysis.
1.3. In results organizer of DPABISurf, the redundant files would not be organized now. In addition, the fmriprep and freesurfer files were backed up, while excluding T1 image that may have privacy information such as face.
2. DPABI_VIEW has a new function "Surface View with DPABISurf_VIEW" now. This function will convert the files to fsaverage surface using freesurfer's mri_vol2surf command. Then the results were displayed by calling DPABISurf_VIEW to generate surface-based picture. 
 
New features of DPABISurf_V1.1_190725  (download at http://rfmri.org/dpabi):
1. New module for Surface-Based Temporal Dynamic Analysis (DPABI_TDA_Surf) was added. Dynamic regional indices (ALFF, fALFF, ReHo, Degree Centrality and Global Signal Correlation) and dynamic functional connectivity could be automatically calculated by one click through DPABI_TDA_Surf (with DPABISurf preprocessed data). The statistics maps (CV, Mean and SD) of the dynamic regional indices would also be generated by DPABI_TDA_Surf. A neuroimaging index which measures the concordance of the dynamic regional indices is incorporated into DPABI_TDA_Surf. Please see more details at: Yan, C.-G., Yang, Z., Colcombe, S.J., Zuo, X.-N., Milham, M.P., 2017. Concordance among indices of intrinsic brain function: insights from inter-individual variation and temporal dynamics. Sci Bull 62, 1572-1584.
2. The calculation of degree centrality now considers a vertex's correlation to both left and right hemispheres.
3. Standardization considers bilateral hemispheres. 
4. Smooth function after Standardization was added.
5. If ICA-AROMA was chosen, no head motion realign parameters would be regressed out.
6. Fixed compatibility issues with old matlab versions.
7. The default surface-based smoothing kernel changed to 6mm instead of 10mm.
8. The DPABISurf results organizing function was added to the R-fMRI Maps Project.
9. Output an excel table for the volume of subcortical structures (calculated by freesurfer): {WorkingDir}/Results/AnatVolu/Anat_Segment_Volume.tsv.
10. DPABISurf_VIEW, the surface-based viewer now has a function to yoke between different viewers.
11. Docker updated basing on fMRIPrep 1.4.1. Besides pull from docker hub, the docker file can be also downloaded form baidu (extract code: enmn).
12. Besides the stand alone version of DPABI (with GUI), the compiled version of DPABISurf_run was also added to docker. Users can run DPABISurf_run with scripting. E.g.,
docker run -it --rm -v /My/FreeSurferLicense/Path/license.txt:/opt/freesurfer/license.txt -v /My/Data/Path:/data cgyan/dpabi /bin/bash
/opt/DPABI/DPABI_StandAlone/run_DPABISurf_run_StandAlone.sh ${MCRPath} /data/DPABISurf_Cfg.mat
 
Tips for Linux or Mac OS: please start matlab from terminal in order to reach docker in DPABI (e.g., Linux: matlab; Mac: open /Applications/MATLAB_R2018a.app/).

 

References:

·                Ashburner, J. (2012). SPM: a history. Neuroimage, 62(2), 791-800, doi:10.1016/j.neuroimage.2011.10.025.

·                Avants, B.B., Epstein, C.L., Grossman, M., Gee, J.C. (2008). Symmetric diffeomorphic image registration with cross-correlation: evaluating automated labeling of elderly and neurodegenerative brain. Med Image Anal, 12(1), 26-41, doi:10.1016/j.media.2007.06.004.

·                Cox, R.W. (1996). AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. Comput Biomed Res, 29(3), 162-173.

·                Dale, A.M., Fischl, B., Sereno, M.I. (1999). Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage, 9(2), 179-194, doi:10.1006/nimg.1998.0395.

·                Esteban, O., Markiewicz, C.J., Blair, R.W., Moodie, C.A., Isik, A.I., Erramuzpe, A., Kent, J.D., Goncalves, M., DuPre, E., Snyder, M., Oya, H., Ghosh, S.S., Wright, J., Durnez, J., Poldrack, R.A., Gorgolewski, K.J. (2019). fMRIPrep: a robust preprocessing pipeline for functional MRI. Nat Methods, 16, 111-116, doi:10.1038/s41592-018-0235-4.

·                Gorgolewski, K.J., Auer, T., Calhoun, V.D., Craddock, R.C., Das, S., Duff, E.P., Flandin, G., Ghosh, S.S., Glatard, T., Halchenko, Y.O., Handwerker, D.A., Hanke, M., Keator, D., Li, X., Michael, Z., Maumet, C., Nichols, B.N., Nichols, T.E., Pellman, J., Poline, J.B., Rokem, A., Schaefer, G., Sochat, V., Triplett, W., Turner, J.A., Varoquaux, G., Poldrack, R.A. (2016). The brain imaging data structure, a format for organizing and describing outputs of neuroimaging experiments. Sci Data, 3, 160044, doi:10.1038/sdata.2016.44.

·                Jenkinson, M., Bannister, P., Brady, M., Smith, S. (2002). Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage, 17(2), 825-841.

·                Li, X., Morgan, P.S., Ashburner, J., Smith, J., Rorden, C. (2016). The first step for neuroimaging data analysis: DICOM to NIfTI conversion. J Neurosci Methods, 264, 47-56, doi:10.1016/j.jneumeth.2016.03.001.

·                Tange, O. (2011). Gnu parallel-the command-line power tool. The USENIX Magazine, 36(1), 42-47.

·                Winkler, A.M., Ridgway, G.R., Douaud, G., Nichols, T.E., Smith, S.M. (2016). Faster permutation inference in brain imaging. Neuroimage, 141, 502-516, doi:10.1016/j.neuroimage.2016.05.068.

·                Yan, C.-G., Wang, X.-D., Lu, B. (2021). DPABISurf: data processing & analysis for brain imaging on surface. Sci Bull, 66(24), 2453-2455, doi: https://doi.org/10.1016/j.scib.2021.09.016.

·                Yan, C.G., Wang, X.D., Zuo, X.N., Zang, Y.F. (2016). DPABI: Data Processing & Analysis for (Resting-State) Brain Imaging. Neuroinformatics, 14(3), 339-351, doi:10.1007/s12021-016-9299-4.

·                Zang, Y., Jiang, T., Lu, Y., He, Y., Tian, L. (2004). Regional homogeneity approach to fMRI data analysis. Neuroimage, 22(1), 394-400, doi:http://dx.doi.org/10.1016/j.neuroimage.2003.12.030.

·                Zang, Y.F., He, Y., Zhu, C.Z., Cao, Q.J., Sui, M.Q., Liang, M., Tian, L.X., Jiang, T.Z., Wang, Y.F. (2007). Altered baseline brain activity in children with ADHD revealed by resting-state functional MRI. Brain Dev, 29(2), 83-91, doi:10.1016/j.braindev.2006.07.002.

·                Zou, Q.-H., Zhu, C.-Z., Yang, Y., Zuo, X.-N., Long, X.-Y., Cao, Q.-J., Wang, Y.-F., Zang, Y.-F. (2008). An improved approach to detection of amplitude of low-frequency fluctuation (ALFF) for resting-state fMRI: Fractional ALFF. Journal of Neuroscience Methods, 172(1), 137-141, doi:http://dx.doi.org/10.1016/j.jneumeth.2008.04.012.

·                Zuo, X.-N., Xing, X.-X. (2014). Test-retest reliabilities of resting-state FMRI measurements in human brain functional connectomics: A systems neuroscience perspective. Neuroscience & Biobehavioral Reviews, 45, 100-118, doi:http://dx.doi.org/10.1016/j.neubiorev.2014.05.009.

严老师您好,我想研究一下海马尾部和全脑的功能连接,但是不知道如何设置海马尾部的ROI,也找不到相关的mask,希望得到老师的帮助!

Hello,
I am using DPABISurf  4.3 200301 and the Demo data with 3 subjects.
 
I got this error below-
fmriprep: error: unrecognized arguments: --ignore-aroma-denoising-errors --template-resampling-grid 2
After re-run fmriprep, these subjects have not run fmriprep yet:
    'sub-Sub001'
    'sub-Sub002'
    'sub-Sub003'
 
I tried re-running, starting with preprocessing with fmriprep and giving BIDS as my starting directory. It still fails witht he error/s below.
Appreciate any suggestions. Thank you.
 
fmriprep: error: unrecognized arguments: --ignore-aroma-denoising-errors --template-resampling-grid 2
 
You are using fMRIPrep-20.0.1, and a newer version of fMRIPrep is available: 20.0.5.
Please check out our documentation about how and when to upgrade:
https://fmriprep.readthedocs.io/en/latest/faq.html#upgrading
WARNING: Version 20.0.1 of fMRIPrep (current) has been FLAGGED
(reason: fsLR resampling error (problematic only when using the --cifti-output flag)).
That means some severe flaw was found in it and we strongly discourage its usage.
 
Error using DPABISurf_run (line 518)
Error detected during running fmriprep, please check!
 
Error in DPABISurf_Pipeline>pushbuttonDPABISurfRun_Callback (line 1693)
[Error, Cfg]=DPABISurf_run(handles.Cfg);
 
Error in gui_mainfcn (line 95)
        feval(varargin{:});
 
Error in DPABISurf_Pipeline (line 42)
    gui_mainfcn(gui_State, varargin{:});
 
Error in
matlab.graphics.internal.figfile.FigFile/read>@(hObject,eventdata)DPABISurf_Pipeline('pushbuttonDPABISurfRun_Callback',hObject,eventdata,guidata(hObject)) 
Error while evaluating UIControl Callback.
 

 

There is a minor bug when using ICA-AROMA since fMRIPrep-20.0.1 changed the input format.

For "Version 20.0.1 of fMRIPrep (current) has been FLAGGED", you don't need to worry because DPABISurf doesn't use cifti-output.

As for now, please download our development version from https://github.com/Chaogan-Yan/DPABI

Please let me know if it works.

Thank you, with the GitHub version fMRIprep started running. 

After few hours, the pipeline has stopped with subject 002. Though I don’t think it finished 001 successfully either.

Matlab window has following error. Please let me know how to fix this. Thank you.

-------------------------------------------------------------

[Node] Running "lta2itk_fwd" ("niworkflows.interfaces.freesurfer.PatchedLTAConvert"), a CommandLine Interface with command:

lta_convert --inlta /data/fmriprepwork/sub-Sub002/fmriprep_wf/single_subject_Sub002_wf/anat_preproc_wf/surface_recon_wf/t1w2fsnative_xfm/out.lta --outitk /data/fmriprepwork/sub-Sub002/fmriprep_wf/single_subject_Sub002_wf/anat_preproc_wf/anat_derivatives_wf/lta2itk_fwd/out.txt

200327-06:36:15,538 nipype.workflow INFO:

                 [Node] Finished "fmriprep_wf.single_subject_Sub002_wf.anat_preproc_wf.anat_derivatives_wf.lta2itk_fwd".

200327-06:36:19,0 nipype.workflow ERROR:

                 could not run node: fmriprep_wf.single_subject_Sub002_wf.anat_preproc_wf.brain_extraction_wf.norm

You are using fMRIPrep-20.0.1, and a newer version of fMRIPrep is available: 20.0.5. Please check out our documentation about how and when to upgrade:

https://fmriprep.readthedocs.io/en/latest/faq.html#upgrading

WARNING: Version 20.0.1 of fMRIPrep (current) has been FLAGGED

(reason: fsLR resampling error (problematic only when using the --cifti-output flag)).

That means some severe flaw was found in it and we strongly

discourage its usage.

fMRIPrep failed: Workflow did not execute cleanly. Check log for details

Preprocessing did not finish successfully. Errors occurred while processing data from participants: Sub002 (1). Check the HTML reports for details.

 

/usr/local/miniconda/lib/python3.7/site-packages/scipy/fftpack/basic.py:160: FutureWarning: Using a non-tuple sequence for multidimensional indexing is deprecated; use `arr[tuple(seq)]` instead of `arr[seq]`. In the future this will be interpreted as an array index, `arr[np.array(seq)]`, which will result either in an error or a different result.

 

  z[index] = x

 

老师好按照您说的,在github上下载了最新的工具包。仍然无法跑通这一步。用的是windows系统。

请问老师应该怎么办

 

错误使用 DPABISurf_run (line 591)
运行 fmriprep 时检测到错误,请检查!
 
出错DPABISurf_Pipeline>pushbuttonDPABISurfRun_Callback (line 1712)
[错误,Cfg]=DPABISurf_run(handles.Cfg);
 
出错gui_mainfcn(第95行)
        feval(varargin{:});
 
出错DPABISurf_Pipeline(第42行)
    gui_mainfcn(gui_State, varargin{:});
 
出错
matlab.graphics.internal.figfile.FigFile/read>@(hObject,eventdata)DPABISurf_Pipeline('pushbuttonDPABISurfRun_Callback',hObject,eventdata,guidata(hObject)) 
UIControl Callback 时计算出错。

Hetingting

Wed, 06/10/2020 - 08:10

严老师,

您好!我在用dpabi做预处理的时候总是出现以下报错,想请问您是什么原因?应该怎么解决?

Running DARTEL: Normalize to MNI space for VBM. Modulated version With smooth kernel [8 8 8].
 
 
------------------------------------------------------------------------
Running job #1
------------------------------------------------------------------------
Running 'Normalise to MNI Space'
Failed  'Normalise to MNI Space'
Error using read_hdr (line 32)
Error reading header file "D:\MRIAnalysissoftware\spm12\tpm\TPM.nii".
In file "D:\MRIAnalysissoftware\spm12\@nifti\private\read_hdr.m" (v4967), function "read_hdr" at line 32.
In file "D:\MRIAnalysissoftware\spm12\@nifti\nifti.m" (v4986), function "nifti" at line 26.
In file "D:\MRIAnalysissoftware\spm12\spm_get_space.m" (v6379), function "spm_get_space" at line 18.
In file "D:\MRIAnalysissoftware\spm12\toolbox\DARTEL\spm_dartel_norm_fun.m" (v5668), function "spm_dartel_norm_fun" at line 43.
 
The following modules did not run:
Failed: Normalise to MNI Space
 
Error using MATLABbatch system
Job execution failed. The full log of this run can be found in MATLAB command
window, starting with the lines (look for the line showing the exact #job as
displayed in this error message)
------------------
Running job #1
------------------
 
Error while evaluating UIControl Callback
 
>> 
期待您的回复,在此感谢

严老师您好,

        我在用DPABISurf进行到preprocessing with fmriprep过程中matlab突然被killed

-VM:~/sharefolder/BOLD/T2DM$ matlab
MATLAB is selecting SOFTWARE OPENGL rendering.
Killed
         然后程序就直接闪退了,一直是这样,您有什么解决办法吗?
 

yamazaki020602

Fri, 06/12/2020 - 06:38

        严老师好!我用台式机windows10-pro系统(Intel i3-4160U 3.6GHz  双核4线程,16G内存,硬盘2T),Matlab2013b  运行DPABI_V4.3_200401中DPABISurf 1.3 【最新版docker2.3.0.3(45519),docker分了2个核,10G memory,load from local file <dpabi43_200401docker.tar.gz,licence.txt均顺利安装】,从DPABISurf_DemoData已做好fmriprep的步骤,照着视频设置往下跑(parallel:1#):

Regressing Out Covariates...
.................................................................................................
Regressing Out Covariates finished.
Warning: Several fields match intent type. Using first.
> In @gifti\private\isintent at 56
  In gifti.subsref at 19
  In y_RegressOutCovariates_OnSurf at 31
  In DPABISurf_run>(parfor body) at 920
  In parallel_function>make_general_channel/channel_general at 922
  In remoteParallelFunction at 28 
 
Load mask "".
Error using y_RegressOutCovariates_OnSurf (line 63)
Dimensions of matrices being concatenated are not consistent.
 
Error in DPABISurf_run>(parfor body) (line 920)
                y_RegressOutCovariates_OnSurf(InFile,CovariablesDef,OutFile,'');
 
Error in DPABISurf_run (line 780)
        parfor i=1:Cfg.SubjectNum
 
Error in DPABISurf_Pipeline>pushbuttonDPABISurfRun_Callback (line 1693)
[Error, Cfg]=DPABISurf_run(handles.Cfg);
 
Error in gui_mainfcn (line 96)
        feval(varargin{:});
 
Error in DPABISurf_Pipeline (line 42)
    gui_mainfcn(gui_State, varargin{:});
 
Error in @(hObject,eventdata)DPABISurf_Pipeline('pushbuttonDPABISurfRun_Callback',hObject,eventdata,guidata(hObject))
 
 
Caused by:
    Error using horzcat
    Dimensions of matrices being concatenated are not consistent.
 
Error while evaluating uicontrol Callback
失败〒_〒 
 
没有气馁 
尝试用VNC viewer 运行DPABISurf_StandAlone GUI version,parallel 1#刚开始run时,滚动一串error,速度比matlab慢但还是能跑下去,到了smoothing the results...                                   ______________________________
parallel:warning:input is read from the terminal...             \              DPABISurf   Pipeline                        /
parallel:warning: are donging(in which case: YOU ARE... \    working DIR  /data/...                               /
parallel:Warning: consider going through the tutorial...    \         sub_Sub_001                                         /
parallel:Warning:Press CTRL-D to exit.                                 \                                                                           /
 
在这停住了,也没有报错提示 ,CTRL-D也无反应,只能docker stop o(╯□╰)o
 
是不是并行设置问题?盼严老师指教,谢谢!
                      
 

Huiling Li

Wed, 06/24/2020 - 08:43

老师,您好,在处理静息态数据的时候,我想用自己的mask,但是我的mask是91x109x91,然后就没办法开展了,所以我想问一下如何将维度为91x109x91 的脑图转成61x73x61?

Huiling Li

Tue, 06/30/2020 - 08:03

老师,您好,我从HCP数据库下载了他们预处理后的静息态数据,现在想用dpabi分析一下功能连接。问题是每个被试的静息态数据包括两个扫描方向的结果一个是从左向右扫描的静息态数据,如“rfMRI_REST1_LR”,一个是从右向左扫描的静息态数据“rfMRI_REST1_RL”。那如果使用dpabi分析的话,是先将这两种数据合并成一个数据吗?还是说先分开分析,最后将脑结果平均?麻烦了,期待回复!

黄浦江

Mon, 10/05/2020 - 14:01

严老师您好,我叫黄浦江

最近在使用dpabi 预处理T1的数据 Reorient这个过程完成之后会报错,想请问下您该如何解决,问题如下:

Reorienting T1 Image Interactively for sub60: OK
 
Reference to non-existent field 'TRSet'.
 
Error in DPARSFA_run>(parfor body) (line 3918)
                TR = AutoDataProcessParameter.TRSet(i,iFunSession);
 
Error in parallel_function (line 466)
            F(base, limit);
 
Error in DPARSFA_run (line 3915)
        parfor i=1:AutoDataProcessParameter.SubjectNum
 
Error in DPARSFA>pushbuttonRun_Callback (line 1845)
    [Error]=DPARSFA_run(handles.Cfg);
 
Error in gui_mainfcn (line 96)
        feval(varargin{:});
 
Error in DPARSFA (line 30)
    gui_mainfcn(gui_State, varargin{:});
 
Error while evaluating uicontrol Callback
 
期待您的回复

严老师您好:

我在使用如下的afni指令进行图像的3D重构时,总是会遇到问题。

脚本指令:

uniq_images IMG*.dcm > uniq_image_list.txt

Dimon -infile_list uniq_image_list.txt \
      -gert_create_dataset             \
      -gert_write_as_nifti             \
      -gert_to3d_prefix T1_01.3D          \
      -gert_outdir ..                  \
      -dicom_org                       \
      -use_last_elem                   \
      -save_details Dimon.details      \
      -gert_quit_on_err
\rm uniq_image_list.txt
 
 
遇到的问题:
 
** ERROR: DICOM file IMG-0015-00179.dcm: not MONOCHROME!
** failed to read file 'IMG-0015-00179.dcm' as dicom
** ERROR: DICOM file IMG-0015-00180.dcm: not MONOCHROME!
** failed to read file 'IMG-0015-00180.dcm' as dicom
 
** no_wait: no volume found in 180 files
-- writing file list to Dimon.details...
no_vol : 
 
请问该如何解决呢?

 

老师您好:

我使用DPARSF进行任务态数据的预处理,发现一开始在FunRaw同样数量,同样大小的dcm文件,跑完预处理之后文件的大小在不同的run里不一样,大约差2mb,想请问一下发生这种情况的原因以及这样的数据会对后期的数据分析有影响吗?谢谢老师!

Dear, Prof. Yan

After re-run thefmriprep failed subjects, how should I choice the options?

'''

ERROR: could not determine type of /data/freesurfer/sub_1001Post/surf/lh.volume
ERROR: could not determine type of /data/freesurfer/sub_1002Post/surf/lh.volume
ERROR: could not determine type of /data/freesurfer/sub_1003Post/surf/lh.volume
ERROR: could not determine type of /data/freesurfer/sub_1004Post/surf/lh.volume
ERROR: could not determine type of /data/freesurfer/sub_1005Post/surf/lh.volume
ERROR: could not determine type of /data/freesurfer/sub_1008Post/surf/lh.volume
ERROR: could not determine type of /data/freesurfer/sub_1009Post/surf/lh.volume
ERROR: could not determine type of /data/freesurfer/sub_1013Post/surf/lh.volume
ERROR: could not determine type of /data/freesurfer/sub_1015Post/surf/lh.volume
ERROR: could not determine type of /data/freesurfer/sub_1001Post/surf/rh.volume
ERROR: could not determine type of /data/freesurfer/sub_1002Post/surf/rh.volume
ERROR: could not determine type of /data/freesurfer/sub_1003Post/surf/rh.volume
...
 
Error in DPABISurf_run (line 536)
    y_Organize_fmriprep(Cfg);
 
Error in DPABISurf_Pipeline>pushbuttonDPABISurfRun_Callback (line 1693)
[Error, Cfg]=DPABISurf_run(handles.Cfg);
 
Error in gui_mainfcn (line 95)
        feval(varargin{:});
 
Error in DPABISurf_Pipeline (line 42)
    gui_mainfcn(gui_State, varargin{:});
 
Error in
matlab.graphics.internal.figfile.FigFile/read>@(hObject,eventdata)DPABISurf_Pipeline('pushbuttonDPABISurfRun_Callback',hObject,eventdata,guidata(hObject)) 
Error while evaluating UIControl Callback.

'''

 

best,

Hui Zheng

马丽华

Fri, 04/16/2021 - 09:24

[environment]
cpu_count = 4
exec_env = "docker"
free_mem = 2.7
overcommit_policy = "heuristic"
overcommit_limit = "50%"
nipype_version = "1.5.1"
templateflow_version = "0.6.3"
version = "20.2.1"

[execution]
bids_dir = "/data/BIDS"
bids_database_dir = "/data/fmriprepwork/sub-sub0983/20210414-115444_7e921fc7-a77e-450c-8825-95139b8682d7/bids_db"
bids_description_hash = "ceb55f9ef8d5658e7418760e3be9884a4199d08a3fe8d1097cc0ad9faf5ae0a2"
boilerplate_only = false
sloppy = false
debug = []
fmriprep_dir = "/data/fmriprep"
fs_license_file = "/opt/freesurfer/license.txt"
fs_subjects_dir = "/data/freesurfer"
layout = "BIDS Layout: .../data/BIDS | Subjects: 1 | Sessions: 0 | Runs: 0"
log_dir = "/data/fmriprep/logs"
log_level = 25
low_mem = true
md_only_boilerplate = false
notrack = false
output_dir = "/data"
output_layout = "legacy"
output_spaces = "fsaverage:den-10k MNI152NLin2009cAsym:res-native"
reports_only = false
run_uuid = "20210414-115444_7e921fc7-a77e-450c-8825-95139b8682d7"
participant_label = [ "sub0983",]
templateflow_home = "/home/fmriprep/.cache/templateflow"
work_dir = "/data/fmriprepwork/sub-sub0983"
write_graph = false

[workflow]
anat_only = false
aroma_err_on_warn = false
aroma_melodic_dim = -200
bold2t1w_dof = 6
bold2t1w_init = "register"
cifti_output = false
fmap_bspline = false
force_syn = false
hires = true
ignore = []
longitudinal = false
medial_surface_nan = false
regressors_all_comps = false
regressors_dvars_th = 1.5
regressors_fd_th = 0.5
run_reconall = true
skull_strip_fixed_seed = false
skull_strip_template = "OASIS30ANTs"
skull_strip_t1w = "force"
spaces = "fsaverage:den-10k MNI152NLin2009cAsym:res-native"
use_aroma = false
use_syn_sdc = false

[nipype]
crashfile_format = "txt"
get_linked_libs = false
nprocs = 4
omp_nthreads = 3
plugin = "MultiProc"
resource_monitor = true
stop_on_first_crash = false

[seeds]
master = 536
ants = 25065

[nipype.plugin_args]
maxtasksperchild = 1
raise_insufficient = false

mazrooyisebd

Sat, 06/05/2021 - 22:03

Hello,

I am new to DPABI Surf and wanted to use its ICA-AROMA noise remove method for my own analysis. I have installed the Docker as it described in the website and got the varification for it as well. However I still get error which I hope someone here can help me fixing it. I got the following lines of errors:

ERROR: could not determine type of /data/freesurfer/sub-NV163V1/surf/rh.area

ERROR: could not determine type of /data/freesurfer/sub-NV163V1/surf/ lh.curv

ERROR: could not determine type of /data/freesurfer/sub-NV163V1/surf/ rh.curv

ERROR: could not determine type of /data/freesurfer/sub-NV163V1/surf/ lh.sulc

ERROR: could not determine type of /data/freesurfer/sub-NV163V1/surf/ rh.sulc

ERROR: could not determine type of /data/freesurfer/sub-NV163V1/surf/ lh.volume

ERROR: could not determine type of /data/freesurfer/sub-NV163V1/surf/ rh.volume

Error using y_Organize_fmriprep>(parfor body) (line 118)

cp: cannot stat

'/mnt/MXMVPGroupShare/JXH108/NV_Sanders/Original_Data/NoDelyn/NoDelyn_12_16_2020/fmriprep/sub-NV004V1/anat/*_space-MNI152NLin2009cAsym_*':

No such file or directory

 

Error in y_Organize_fmriprep (line 116)

parfor i=1:Cfg.SubjectNum

 

Error in DPABISurf_run (line 479)

    y_Organize_fmriprep(Cfg);

 

Error in DPABISurf_Pipeline>pushbuttonDPABISurfRun_Callback (line 1693)

[Error, Cfg]=DPABISurf_run(handles.Cfg);

 

Error in gui_mainfcn (line 95)

        feval(varargin{:});

 

Error in DPABISurf_Pipeline (line 42)

    gui_mainfcn(gui_State, varargin{:});

 

Error in

@(hObject,eventdata)DPABISurf_Pipeline('pushbuttonDPABISurfRun_Callback',hObject,eventdata,guidata(hObject))

 

Error while evaluating UIControl Callback

Thanks in advance for your helps!

three_brains

Sun, 06/20/2021 - 07:15

老师好!
我在预处理的时候点完reorient的图像后,便会出现以下报错,请问该怎么解决呢?谢谢老师!

Bet finished
Copying T1 image Files:Sub_001 OK
Copying T1 image Files:Sub_002 OK
Copying T1 image Files:Sub_003 OK
 
Error using DPARSFA_run (line 1620)
Index exceeds matrix dimensions.
 
Error in DPARSFA>pushbuttonRun_Callback (line 1795)
    [Error]=DPARSFA_run(handles.Cfg);
 
Error in gui_mainfcn (line 95)
        feval(varargin{:});
 
Error in DPARSFA (line 30)
    gui_mainfcn(gui_State, varargin{:});
 
Error while evaluating UIControl Callback.

asus_learner

Wed, 06/23/2021 - 09:38

严老师您好,我在网上下载了其他的数据想进行数据处理,但是跑了很多次都在同一个地方报错,我自己不太可能解决这个问题,希望得到老师的帮助!以下是报错的内容

错误使用 gunzip>gunzipwrite (line 234)
尝试对文件“sub-25629_ses-2_T1w.nii”执行GUNZIP时突然终止输入流。
 
出错 gunzip>gunzipEntries (line 154)
    名称{k} = gunzipwrite(entries(k).file, outputDir, baseName, streamCopier);
 
出错 gunzip(第 90 行)
   名称 = gunzipEntries(entries, outputDir);
 
出错 DPARSFA_run (line 1584)
                            gunzip(DirImg(j).name);
 
出错 DPARSFA>pushbuttonRun_Callback (line 1847)
    [错误,Cfg]=DPARSFA_run(handles.Cfg);
 
出错gui_mainfcn(第95行)
        feval(varargin{:});
 
出错 DPARSFA(第 30 行)
    gui_mainfcn(gui_State, varargin{:});
 
UIControl Callback 时计算出错。

严老师您好,我最近跑DPABISurf会出现报错,有几个地方我觉得很疑惑

(1)fmriprep这一步被报告没跑成功的被试在其对应的网页上却显示:

“Errors
 
No errors to report!”
 
(2)被试文件夹里的log文件会显示:
[environment]
cpu_count = 48
exec_env = "docker"
free_mem = 18.7
overcommit_policy = "heuristic"
overcommit_limit = "50%"
nipype_version = "1.5.1"
templateflow_version = "0.6.3"
version = "20.2.0"
 
[execution]
bids_dir = "/data/BIDS"
bids_database_dir = "/data/fmriprepwork/sub-hyr/20210908-073726_1b31ac36-a899-4bb1-94d9-f8fce44b11cf/bids_db"
bids_description_hash = "ceb55f9ef8d5658e7418760e3be9884a4199d08a3fe8d1097cc0ad9faf5ae0a2"
boilerplate_only = false
sloppy = false
debug = []
fs_license_file = "/opt/freesurfer/license.txt"
fs_subjects_dir = "/data/freesurfer"
layout = "BIDS Layout: .../data/BIDS | Subjects: 1 | Sessions: 0 | Runs: 0"
log_dir = "/data/fmriprep/logs"
log_level = 25
low_mem = true
md_only_boilerplate = false
notrack = false
output_dir = "/data"
output_spaces = "fsaverage:den-10k MNI152NLin2009cAsym:res-2"
reports_only = false
run_uuid = "20210908-073726_1b31ac36-a899-4bb1-94d9-f8fce44b11cf"
participant_label = [ "hyr",]
templateflow_home = "/home/fmriprep/.cache/templateflow"
work_dir = "/data/fmriprepwork/sub-hyr"
write_graph = false
 
[workflow]
anat_only = false
aroma_err_on_warn = false
aroma_melodic_dim = -200
bold2t1w_dof = 6
bold2t1w_init = "register"
cifti_output = false
fmap_bspline = false
force_syn = false
hires = true
ignore = [ "slicetiming",]
longitudinal = false
medial_surface_nan = false
regressors_all_comps = false
regressors_dvars_th = 1.5
regressors_fd_th = 0.5
run_reconall = true
skull_strip_fixed_seed = false
skull_strip_template = "OASIS30ANTs"
skull_strip_t1w = "force"
spaces = "fsaverage:den-10k MNI152NLin2009cAsym:res-2 MNI152NLin6Asym:res-2"
use_aroma = true
use_syn_sdc = false
 
[nipype]
crashfile_format = "txt"
get_linked_libs = false
nprocs = 48
omp_nthreads = 8
plugin = "MultiProc"
resource_monitor = true
stop_on_first_crash = false
 
[seeds]
master = 30561
ants = 13279
 
[nipype.plugin_args]
maxtasksperchild = 1
raise_insufficient = false
我看有类似提问,您回答说应该是没有错,可以继续往下跑,但是我这边跑到nipype这一步,matlab就会报告
nipype.workflow WARNING:
 
Could not retrieve profiling information for node "fmriprep_wf.single_subject_hyr_wf.func_preproc_task_rest_wf.final_boldref_wf.enhance_and_skullstrip_bold_wf.n4_correct" (mapflow 1/1).
 
(3)matlab最后会报告:fMRIPrep finished successfully!;同时报告These subjects have failed run fmriprep(文件夹里每个被试都会被报告跑失败了);最后会报错,如下:
Array indices must be positive integers or logical values.
 
Error in y_ReRunfmriprepFailedSubjects (line 164)
        if strcmpi(Cfg.Normalize.VoxelSize(end-1:end),'mm')
 
Error in DPABISurf_run (line 509)
    [FailedID WaitingID SuccessID]=y_ReRunfmriprepFailedSubjects(Cfg,Cfg.WorkingDir,Cfg.SubjectID);
 
Error in DPABISurf_Pipeline>pushbuttonDPABISurfRun_Callback (line 1693)
[Error, Cfg]=DPABISurf_run(handles.Cfg);
 
Error in gui_mainfcn (line 95)
        feval(varargin{:});
 
Error in DPABISurf_Pipeline (line 42)
    gui_mainfcn(gui_State, varargin{:});
 
Error in
matlab.graphics.internal.figfile.FigFile/read>@(hObject,eventdata)DPABISurf_Pipeline('pushbuttonDPABISurfRun_Callback',hObject,eventdata,guidata(hObject)) 
Error while evaluating UIControl Callback.
 
希望可以得到老师的帮助,谢谢!

 

Yapeng Qi

Thu, 09/30/2021 - 13:49

老师您好!

 
我在用dpabi做数据预处理时候遇到这样的问题。我的数据是并采加速采集的,共有62层,TR为2s,全脑分为前后两部分、各31层同时采样,采集顺序如下:
   前半部分:[2:2:30 1:2:31]
   后半部分:[33:2:61 32:2:62]
   每部分的采集时刻为:0, 62.5, 127.5, 190, 255, 317.5, 382.5, 445, 510, 575, 637.5, 702.5, 765, 830, 892.5, 957.5, 1020, 1085, 1150, 1212.5, 1277.5, 1340, 1405, 1467.5, 1532.5, 1595, 1660, 1722.5, 1787.5, 1852.5, 1915。
 
 
请问这种情况可以使用dpabi进行时间层校正吗?如果可以的话slice order应该怎样填写呢?是否需要修改脚本?
 
期待您的回复,谢谢老师!

xieht0123

Wed, 10/13/2021 - 06:49

严老师您好:

用西门子3.0T核磁采集的fmri序列,同时采集了bold_filed_maping序列,如果想做Field map correction,原始数据该怎么设置?

您的示范数据—FieldMap中放了3个文件夹,谢谢!

请参阅:

https://github.com/Puyunfashi/Suzhou_rumination/blob/master/sort_data/Organize_FieldMap.m

https://github.com/Puyunfashi/Suzhou_rumination/blob/master/sort_data/c_DivideFieldMapMagnitudeFiles.m

严老师:

您好!

我目前在做一组儿童的静息态fMRI预处理,因此在https://www.nitrc.org/projects/chn-pd这里下载了贺永老师团队创建的一个中国儿童标准脑模板,利用其TPM先验模板来做segment和normalize,但得到的灰质、白质分割结果是一个这样的图,因此我想请问老师为什么会出现这样的结果,有什么方法可以解决吗?如果我做更换模板这样的操作,需要注意哪些问题呢?

严老师,您好!

以上这个结果是我做了替换后得到的分割结果,而且同时也用spm尝试了一下,发现也是这样的,目前没有找到问题出在哪,麻烦老师帮忙解答一下,谢谢!

enjoylife

Fri, 11/19/2021 - 10:16

 

严老师您好,我有个问题想请教您,您在DPABISurf输出文件夹与计算结果概述 的课程里提到-----存在Data>fmriprep>sub-Sub001>log文件夹的话,一般是跑fmriprep可能出问题了。老师我想问您,我还要继续运行DPABIsurf吗?

Parallel pool using the 'local' profile is shutting down.
Starting parallel pool (parpool) using the 'local' profile ...
Connected to the parallel pool (number of workers: 12).
Preprocessing with fmriprep, this process is very time consuming, please be patient..