[RESOLVED] Hcp_fmri_volume : Data requirements

Hi Jure,

I’m having problems with the hcp_fmri_volume step when I’m testing the effect of parameter tuning with parallel times using my own data.

I’ll start by providing you with mapping.txt and batch.txt used in parallel testing.

mapping.txt (136 Bytes)

batch.txt (8.0 KB)

In the mapping.txt file I only used T1w,T2w,Bold_rest from all the sequences, perhaps because I only provided the bold image and not the corresponding spin-echo field maps, which made it problematic to do the hcp_fmri_volume processing. The terminal output is as follows:

(base) jiayf@EyeOpener:~/qunex$ qunex_container hcp_fmri_volume    --sessionsfolder="${STUDY_FOLDER}/sessions"     --batchfile="${STUDY_FOLDER}/processing/batch.txt"     --parelements=1  --container="${QUNEX_CONTAINER}"

---> QuNex will run the command over 4 sessions. It will utilize:

    Maximum sessions run in parallel for a job: 1.
    Maximum elements run in parallel for a session: 1.
    Up to 1 processes will be utilized for a job.

    Job #1 will run sessions: 31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004
(base) jiayf@EyeOpener:~/qunex$ ---> unsetting the following environment variables: PATH MATLABPATH PYTHONPATH QUNEXVer TOOLS QUNEXREPO QUNEXPATH QUNEXEXTENSIONS QUNEXLIBRARY QUNEXLIBRARYETC TemplateFolder FSL_FIXDIR FREESURFERDIR FREESURFER_HOME FREESURFER_SCHEDULER FreeSurferSchedulerDIR WORKBENCHDIR DCMNIIDIR DICMNIIDIR MATLABDIR MATLABBINDIR OCTAVEDIR OCTAVEPKGDIR OCTAVEBINDIR RDIR HCPWBDIR AFNIDIR PYLIBDIR FSLDIR FSLBINDIR PALMDIR QUNEXMCOMMAND HCPPIPEDIR CARET7DIR GRADUNWARPDIR HCPPIPEDIR_Templates HCPPIPEDIR_Bin HCPPIPEDIR_Config HCPPIPEDIR_PreFS HCPPIPEDIR_FS HCPPIPEDIR_FS_CUSTOM HCPPIPEDIR_PostFS HCPPIPEDIR_fMRISurf HCPPIPEDIR_fMRIVol HCPPIPEDIR_tfMRI HCPPIPEDIR_dMRI HCPPIPEDIR_dMRITract HCPPIPEDIR_Global HCPPIPEDIR_tfMRIAnalysis HCPCIFTIRWDIR MSMBin HCPPIPEDIR_dMRITractFull HCPPIPEDIR_dMRILegacy AutoPtxFolder EDDYCUDA USEOCTAVE QUNEXENV CONDADIR MSMBINDIR MSMCONFIGDIR R_LIBS FSL_FIX_CIFTIRW FSFAST_HOME SUBJECTS_DIR MINC_BIN_DIR MNI_DIR MINC_LIB_DIR MNI_DATAPATH FSF_OUTPUT_FORMAT ANTSDIR CUDIMOT

========================================================================
Generated by QuNex
------------------------------------------------------------------------
Version: 1.1.0 [QIO]
User: root
System: a3ea1a2f307e
OS: Debian Linux #1 SMP Tue Nov 5 00:21:55 UTC 2024
------------------------------------------------------------------------

        ██████\                  ║      ██\   ██\
       ██  __██\                 ║      ███\  ██ |
       ██ /  ██ |██\   ██\       ║      ████\ ██ | ██████\ ██\   ██\
       ██ |  ██ |██ |  ██ |      ║      ██ ██\██ |██  __██\\██\ ██  |
       ██ |  ██ |██ |  ██ |      ║      ██ \████ |████████ |\████  /
       ██ ██\██ |██ |  ██ |      ║      ██ |\███ |██   ____|██  ██\
       \██████ / \██████  |      ║      ██ | \██ |\███████\██  /\██\
        \___███\  \______/       ║      \__|  \__| \_______\__/  \__|
            \___|                ║


                       DEVELOPED & MAINTAINED BY:

               Mind & Brain Lab, University of Ljubljana
                       Cho Lab, Yale University

                      COPYRIGHT & LICENSE NOTICE:

Use of this software is subject to the terms and conditions defined in
QuNex LICENSES which can be found in the LICENSES folder of the QuNex
repository or at https://qunex.yale.edu/qunex-registration
========================================================================

---> Setting up Octave


.......................... Running QuNex v1.1.0 [QIO] ..........................


--- Full QuNex call for command: hcp_fmri_volume

qunex hcp_fmri_volume --sessionsfolder="/home/jiayf/qunex/parallel/sessions" --parelements="1" --batchfile="/home/jiayf/qunex/parallel/processing/batch.txt" --sessions="31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004"

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


# Generated by QuNex 1.1.0 [QIO] on 2025-03-17_01.30.02.136063#
=================================================================
qunex hcp_fmri_volume \
  --sessionsfolder="/home/jiayf/qunex/parallel/sessions" \
  --parelements="1" \
  --sessions="/home/jiayf/qunex/parallel/processing/batch.txt" \
  --sessionids="31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004" \
=================================================================

Starting multiprocessing sessions in /home/jiayf/qunex/parallel/processing/batch.txt with a pool of 4 concurrent processes


Adding processing of session 31-JA-JNLA-0001 to the pool at Monday, 17. March 2025 01:30:02

Adding processing of session 31-JA-JNLA-0002 to the pool at Monday, 17. March 2025 01:30:02

Adding processing of session 31-JA-JNLA-0003 to the pool at Monday, 17. March 2025 01:30:02

Adding processing of session 31-JA-JNLA-0004 to the pool at Monday, 17. March 2025 01:30:02
=======================================================================================================================================
=======================================================================================================================================
=======================================================================================================================================
=======================================================================================================================================

------------------------------------------------------------
Session id: 31-JA-JNLA-0001
[started on Monday, 17. March 2025 01:30:02]
Running HCP fMRI Volume pipeline [LegacyStyleData] ...
---> PreFS results present.
---> FS results present.
---> PostFS results present.
---> Looking for spin echo fieldmap set images [NONE/NONE].
     ... TOPUP configuration file present

---> Preprocessing settings (unwarpdir, refimage, moveref, seimage) for BOLD 1
     ... phase encoding direction: PA
     ... unwarp direction: j
     ... using image specific EchoSpacing: 0.000530002 s
     ... ERROR: No spin echo fieldmap set images present!
     ... bold image present
     ... reference image not used
     ... using the HCPpipelines default BOLD mask

Running 1 BOLD images in parallel
---> ERROR: something missing, skipping this BOLD!

HCP fMRIVolume completed on Monday, 17. March 2025 01:30:02
------------------------------------------------------------

------------------------------------------------------------
Session id: 31-JA-JNLA-0002
[started on Monday, 17. March 2025 01:30:02]
Running HCP fMRI Volume pipeline [LegacyStyleData] ...
---> PreFS results present.
---> FS results present.
---> PostFS results present.
---> Looking for spin echo fieldmap set images [NONE/NONE].
     ... TOPUP configuration file present

---> Preprocessing settings (unwarpdir, refimage, moveref, seimage) for BOLD 1
     ... phase encoding direction: PA
     ... unwarp direction: j
     ... using image specific EchoSpacing: 0.000530002 s
     ... ERROR: No spin echo fieldmap set images present!
     ... bold image present
     ... reference image not used
     ... using the HCPpipelines default BOLD mask

Running 1 BOLD images in parallel
---> ERROR: something missing, skipping this BOLD!

HCP fMRIVolume completed on Monday, 17. March 2025 01:30:02
------------------------------------------------------------

------------------------------------------------------------
Session id: 31-JA-JNLA-0003
[started on Monday, 17. March 2025 01:30:02]
Running HCP fMRI Volume pipeline [LegacyStyleData] ...
---> PreFS results present.
---> FS results present.
---> PostFS results present.
---> Looking for spin echo fieldmap set images [NONE/NONE].
     ... TOPUP configuration file present

---> Preprocessing settings (unwarpdir, refimage, moveref, seimage) for BOLD 1
     ... phase encoding direction: PA
     ... unwarp direction: j
     ... using image specific EchoSpacing: 0.000530002 s
     ... ERROR: No spin echo fieldmap set images present!
     ... bold image present
     ... reference image not used
     ... using the HCPpipelines default BOLD mask

Running 1 BOLD images in parallel
---> ERROR: something missing, skipping this BOLD!

HCP fMRIVolume completed on Monday, 17. March 2025 01:30:02
------------------------------------------------------------

------------------------------------------------------------
Session id: 31-JA-JNLA-0004
[started on Monday, 17. March 2025 01:30:02]
Running HCP fMRI Volume pipeline [LegacyStyleData] ...
---> PreFS results present.
---> FS results present.
---> PostFS results present.
---> Looking for spin echo fieldmap set images [NONE/NONE].
     ... TOPUP configuration file present

---> Preprocessing settings (unwarpdir, refimage, moveref, seimage) for BOLD 1
     ... phase encoding direction: PA
     ... unwarp direction: j
     ... using image specific EchoSpacing: 0.000530002 s
     ... ERROR: No spin echo fieldmap set images present!
     ... bold image present
     ... reference image not used
     ... using the HCPpipelines default BOLD mask

Running 1 BOLD images in parallel
---> ERROR: something missing, skipping this BOLD!

HCP fMRIVolume completed on Monday, 17. March 2025 01:30:02
------------------------------------------------------------


---> Final report for command hcp_fmri_volume
... 31-JA-JNLA-0001 ---> HCP fMRI Volume: bolds 1 not ready
... 31-JA-JNLA-0002 ---> HCP fMRI Volume: bolds 1 not ready
... 31-JA-JNLA-0003 ---> HCP fMRI Volume: bolds 1 not ready
... 31-JA-JNLA-0004 ---> HCP fMRI Volume: bolds 1 not ready
---> Not all tasks completed fully!

The information I’ve reviewed indicates that HCPpipeline requires paired spin echo fieldmaps (AP+PA) to generate the aberration field (topup) for processing, and that a lack of fieldmaps will result in an inability to compute the aberration field, with the error No spin echo fieldmap set images present. Some alternatives state that a GRE fieldmap (gradient echo fieldmap) can be used instead, I see that I have MAG_IMAGES_0022 and PHA_IMAGES_0023 in my sequences folder, does HCPpipeline support this substitution?

I can provide you with the session.txt of the input data, which contains sequence information that may help you in your judgment.
session.txt (3.2 KB)

Addendum: I realized that my dmri data may not meet the HCPpipeline data requirements either, as one batch of data is 64 directions, while the other batch only even has 48 directions!!! I hope my fmri data processing is still saved!:sob:

Best,
YJia

Hi YJia,

Do not worry, with QuNex you can process pretty much any kind of data :slight_smile: , if your acquisition protocol deviates from HCP’s then preprocessed data quality will be most likely a bit lower. Most data I worked with and I would also say most data from our users does not exactly follow HCP acquisition, meaning that we often work with data similar as yours. Based on the attached files, your acquisition protocol looks OK and I think there is nothing to worry about.

Are images 221 : Mag_Images and 231 : Pha_Images in your session.txt fieldmap magnitude and phase? If that is the case, then you need to add the following to the mapping file:

Mag_Images => FM-Magnitude
Pha_Images => FM-Phase

See HCP Pipeline Mapping in QuNex — QuNex documentation for additional details. After doing this, you need to rerun relevant QuNex commands (create_session_info, setup_hcp, create_batch) with --overwrite=yes. Note that your structural processing results (HCP FreeSurfer chain) will be preserved so nothing to worry here. Note though that the quality of structural processing can be also improved with fieldmaps through the --hcp_avgrdcmethod parameter, so you might want to reprocess that as well.

After that you need to set the relevant hcp_fmri_volume parameters:

qunex_container hcp_fmri_volume \
  ...
  --hcp_bold_dcmethod="SiemensFieldMap" \
  ...

I wrote this from the top of my head. I currently do not have access to our main processing system where I could double check this, so there is a chance that this might need some tweaks, but the general logic should be sound.

Best, Jure

Hi Jure.

I’m going to try this right away and look forward to bringing you good news!

Best,
YJia

Hi, Jure

I have reworked the data in the .zip and the mapping.txt as you suggested, and the revised mapping.txt is as follows:

#  parallel mapping file
#  ----------------

41                  => T1w
51                  => T2w
121                => bold:rest
161                => DWI:dir98_AP
221                => FM-Magnitude
231                => FM-Phase

After that I executed create_session_info, setup_hcp, create_batch) with --overwrite=yes (the setup_hcp directive doesn’t seem to have an overwrite parameter) as you requested. Because I encountered the following error while executing the setup_hcp` command:

qunex_container setup_hcp     --sessionsfolder="${STUDY_FOLDER}/sessions"
--batchfile="${STUDY_FOLDER}/processing/batch.txt"   --overwrite=yes  --container="${QUNEX_CONTAINER}"

---> QuNex will run the command over 4 sessions. It will utilize:

    Maximum sessions run in parallel for a job: 1.
    Maximum elements run in parallel for a session: 1.
    Up to 1 processes will be utilized for a job.

    Job #1 will run sessions: 31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004
(base) jiayf@EyeOpener:~/qunex$ ---> unsetting the following environment variables: PATH MATLABPATH PYTHONPATH QUNEXVer TOOLS QUNEXREPO QUNEXPATH QUNEXEXTENSIONS QUNEXLIBRARY QUNEXLIBRARYETC TemplateFolder FSL_FIXDIR FREESURFERDIR FREESURFER_HOME FREESURFER_SCHEDULER FreeSurferSchedulerDIR WORKBENCHDIR DCMNIIDIR DICMNIIDIR MATLABDIR MATLABBINDIR OCTAVEDIR OCTAVEPKGDIR OCTAVEBINDIR RDIR HCPWBDIR AFNIDIR PYLIBDIR FSLDIR FSLBINDIR PALMDIR QUNEXMCOMMAND HCPPIPEDIR CARET7DIR GRADUNWARPDIR HCPPIPEDIR_Templates HCPPIPEDIR_Bin HCPPIPEDIR_Config HCPPIPEDIR_PreFS HCPPIPEDIR_FS HCPPIPEDIR_FS_CUSTOM HCPPIPEDIR_PostFS HCPPIPEDIR_fMRISurf HCPPIPEDIR_fMRIVol HCPPIPEDIR_tfMRI HCPPIPEDIR_dMRI HCPPIPEDIR_dMRITract HCPPIPEDIR_Global HCPPIPEDIR_tfMRIAnalysis HCPCIFTIRWDIR MSMBin HCPPIPEDIR_dMRITractFull HCPPIPEDIR_dMRILegacy AutoPtxFolder EDDYCUDA USEOCTAVE QUNEXENV CONDADIR MSMBINDIR MSMCONFIGDIR R_LIBS FSL_FIX_CIFTIRW FSFAST_HOME SUBJECTS_DIR MINC_BIN_DIR MNI_DIR MINC_LIB_DIR MNI_DATAPATH FSF_OUTPUT_FORMAT ANTSDIR CUDIMOT

========================================================================
Generated by QuNex
------------------------------------------------------------------------
Version: 1.1.0 [QIO]
User: root
System: da69f526ff52
OS: Debian Linux #1 SMP Tue Nov 5 00:21:55 UTC 2024
------------------------------------------------------------------------

        ██████\                  ║      ██\   ██\
       ██  __██\                 ║      ███\  ██ |
       ██ /  ██ |██\   ██\       ║      ████\ ██ | ██████\ ██\   ██\
       ██ |  ██ |██ |  ██ |      ║      ██ ██\██ |██  __██\\██\ ██  |
       ██ |  ██ |██ |  ██ |      ║      ██ \████ |████████ |\████  /
       ██ ██\██ |██ |  ██ |      ║      ██ |\███ |██   ____|██  ██\
       \██████ / \██████  |      ║      ██ | \██ |\███████\██  /\██\
        \___███\  \______/       ║      \__|  \__| \_______\__/  \__|
            \___|                ║


                       DEVELOPED & MAINTAINED BY:

               Mind & Brain Lab, University of Ljubljana
                       Cho Lab, Yale University

                      COPYRIGHT & LICENSE NOTICE:

Use of this software is subject to the terms and conditions defined in
QuNex LICENSES which can be found in the LICENSES folder of the QuNex
repository or at https://qunex.yale.edu/qunex-registration
========================================================================

---> Setting up Octave


.......................... Running QuNex v1.1.0 [QIO] ..........................


--- Full QuNex call for command: setup_hcp

qunex setup_hcp --sessionsfolder="/home/jiayf/qunex/parallel/sessions" --overwrite="yes" --batchfile="/home/jiayf/qunex/parallel/processing/batch.txt" --sessions="31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004"

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


ERROR: Extra argument overwrite is not valid! Please check your command!

ERROR in running setup_hcp at gmri:
Invalid arguments
Provided unknown arguments: [overwrite]!

The call received was:
(please note that when run through scheduler, all possible parameters,
even non relevant ones are passed)

qunex setup_hcp \
    --sessionsfolder=/home/jiayf/qunex/parallel/sessions \
    --overwrite=yes \
    --batchfile=/home/jiayf/qunex/parallel/processing/batch.txt \
    --sessions=31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004

Please run `qunex setup_hcp --help` to get help for the failed command.

After that I encountered the following when executing the hcp_fmri_volume command:

(base) jiayf@EyeOpener:~/qunex$ qunex_container hcp_fmri_volume    --sessionsfolder="${STUDY_FOLDER}/sessions"     --batchfile="${STUDY_FOLDER}/processing/batch.txt"  --hcp_bold_dcmethod="SiemensFieldMap"   --parelements
=1  --container="${QUNEX_CONTAINER}"

---> QuNex will run the command over 4 sessions. It will utilize:

    Maximum sessions run in parallel for a job: 1.
    Maximum elements run in parallel for a session: 1.
    Up to 1 processes will be utilized for a job.

    Job #1 will run sessions: 31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004
(base) jiayf@EyeOpener:~/qunex$ ---> unsetting the following environment variables: PATH MATLABPATH PYTHONPATH QUNEXVer TOOLS QUNEXREPO QUNEXPATH QUNEXEXTENSIONS QUNEXLIBRARY QUNEXLIBRARYETC TemplateFolder FSL_FIXDIR FREESURFERDIR FREESURFER_HOME FREESURFER_SCHEDULER FreeSurferSchedulerDIR WORKBENCHDIR DCMNIIDIR DICMNIIDIR MATLABDIR MATLABBINDIR OCTAVEDIR OCTAVEPKGDIR OCTAVEBINDIR RDIR HCPWBDIR AFNIDIR PYLIBDIR FSLDIR FSLBINDIR PALMDIR QUNEXMCOMMAND HCPPIPEDIR CARET7DIR GRADUNWARPDIR HCPPIPEDIR_Templates HCPPIPEDIR_Bin HCPPIPEDIR_Config HCPPIPEDIR_PreFS HCPPIPEDIR_FS HCPPIPEDIR_FS_CUSTOM HCPPIPEDIR_PostFS HCPPIPEDIR_fMRISurf HCPPIPEDIR_fMRIVol HCPPIPEDIR_tfMRI HCPPIPEDIR_dMRI HCPPIPEDIR_dMRITract HCPPIPEDIR_Global HCPPIPEDIR_tfMRIAnalysis HCPCIFTIRWDIR MSMBin HCPPIPEDIR_dMRITractFull HCPPIPEDIR_dMRILegacy AutoPtxFolder EDDYCUDA USEOCTAVE QUNEXENV CONDADIR MSMBINDIR MSMCONFIGDIR R_LIBS FSL_FIX_CIFTIRW FSFAST_HOME SUBJECTS_DIR MINC_BIN_DIR MNI_DIR MINC_LIB_DIR MNI_DATAPATH FSF_OUTPUT_FORMAT ANTSDIR CUDIMOT

========================================================================
Generated by QuNex
------------------------------------------------------------------------
Version: 1.1.0 [QIO]
User: root
System: df90ddb8d837
OS: Debian Linux #1 SMP Tue Nov 5 00:21:55 UTC 2024
------------------------------------------------------------------------

        ██████\                  ║      ██\   ██\
       ██  __██\                 ║      ███\  ██ |
       ██ /  ██ |██\   ██\       ║      ████\ ██ | ██████\ ██\   ██\
       ██ |  ██ |██ |  ██ |      ║      ██ ██\██ |██  __██\\██\ ██  |
       ██ |  ██ |██ |  ██ |      ║      ██ \████ |████████ |\████  /
       ██ ██\██ |██ |  ██ |      ║      ██ |\███ |██   ____|██  ██\
       \██████ / \██████  |      ║      ██ | \██ |\███████\██  /\██\
        \___███\  \______/       ║      \__|  \__| \_______\__/  \__|
            \___|                ║


                       DEVELOPED & MAINTAINED BY:

               Mind & Brain Lab, University of Ljubljana
                       Cho Lab, Yale University

                      COPYRIGHT & LICENSE NOTICE:

Use of this software is subject to the terms and conditions defined in
QuNex LICENSES which can be found in the LICENSES folder of the QuNex
repository or at https://qunex.yale.edu/qunex-registration
========================================================================

---> Setting up Octave


.......................... Running QuNex v1.1.0 [QIO] ..........................


--- Full QuNex call for command: hcp_fmri_volume

qunex hcp_fmri_volume --sessionsfolder="/home/jiayf/qunex/parallel/sessions" --hcp_bold_dcmethod="SiemensFieldMap" --parelements="1" --batchfile="/home/jiayf/qunex/parallel/processing/batch.txt" --sessions="31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004"

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


# Generated by QuNex 1.1.0 [QIO] on 2025-03-17_07.09.56.514061#
=================================================================
qunex hcp_fmri_volume \
  --sessionsfolder="/home/jiayf/qunex/parallel/sessions" \
  --hcp_bold_dcmethod="SiemensFieldMap" \
  --parelements="1" \
  --sessions="/home/jiayf/qunex/parallel/processing/batch.txt" \
  --sessionids="31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004" \
=================================================================

Starting multiprocessing sessions in /home/jiayf/qunex/parallel/processing/batch.txt with a pool of 4 concurrent processes


Adding processing of session 31-JA-JNLA-0001 to the pool at Monday, 17. March 2025 07:09:56

Adding processing of session 31-JA-JNLA-0002 to the pool at Monday, 17. March 2025 07:09:56

Adding processing of session 31-JA-JNLA-0003 to the pool at Monday, 17. March 2025 07:09:56

Adding processing of session 31-JA-JNLA-0004 to the pool at Monday, 17. March 2025 07:09:56
=======================================================================================================================================

------------------------------------------------------------
Session id: 31-JA-JNLA-0001
[started on Monday, 17. March 2025 07:09:56]
Running HCP fMRI Volume pipeline [HCPStyleData] ...
---> PreFS results present.
---> FS results present.
---> PostFS results present.

---> Preprocessing settings (unwarpdir, refimage, moveref, seimage) for BOLD 1
     ... phase encoding direction: PA
     ... unwarp direction: j
     ... using image specific EchoSpacing: 0.000530002 s
---> ERROR: No fieldmap number specified for the BOLD image!
     ... bold image present
     ... reference image not used
     ... using the HCPpipelines default BOLD mask

Running 1 BOLD images in parallel

------------------------------------------------------------
Running HCP Pipelines command via QuNex:

/opt/HCP/HCPpipelines/fMRIVolume/GenericfMRIVolumeProcessingPipeline.sh
    --path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/hcp"
    --session="31-JA-JNLA-0001"
    --fmriname="1"
    --fmritcs="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/hcp/31-JA-JNLA-0001/unprocessed/BOLD_1_PA/31-JA-JNLA-0001_BOLD_1_PA.nii.gz"
    --fmriscout="NONE"
    --fmapmag="NONE"
    --fmapphase="NONE"
    --fmapcombined="NONE"
    --echospacing="0.000530002"
    --echodiff="NONE"
    --unwarpdir="j"
    --fmrires="2"
    --dcmethod="SiemensFieldMap"
    --biascorrection="NONE"
    --gdcoeffs="NONE"
    --mctype="MCFLIRT"
    --processing-mode="HCPStyleData"
    --doslicetime="FALSE"
------------------------------------------------------------

---> ERROR: something missing, skipping this BOLD!

HCP fMRIVolume completed on Monday, 17. March 2025 07:09:56
------------------------------------------------------------
=======================================================================================================================================
=======================================================================================================================================
=======================================================================================================================================

------------------------------------------------------------
Session id: 31-JA-JNLA-0003
[started on Monday, 17. March 2025 07:09:56]
Running HCP fMRI Volume pipeline [HCPStyleData] ...
---> PreFS results present.
---> FS results present.
---> PostFS results present.

---> Preprocessing settings (unwarpdir, refimage, moveref, seimage) for BOLD 1
     ... phase encoding direction: PA
     ... unwarp direction: j
     ... using image specific EchoSpacing: 0.000530002 s
---> ERROR: No fieldmap number specified for the BOLD image!
     ... bold image present
     ... reference image not used
     ... using the HCPpipelines default BOLD mask

Running 1 BOLD images in parallel

------------------------------------------------------------
Running HCP Pipelines command via QuNex:

/opt/HCP/HCPpipelines/fMRIVolume/GenericfMRIVolumeProcessingPipeline.sh
    --path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0003/hcp"
    --session="31-JA-JNLA-0003"
    --fmriname="1"
    --fmritcs="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0003/hcp/31-JA-JNLA-0003/unprocessed/BOLD_1_PA/31-JA-JNLA-0003_BOLD_1_PA.nii.gz"
    --fmriscout="NONE"
    --fmapmag="NONE"
    --fmapphase="NONE"
    --fmapcombined="NONE"
    --echospacing="0.000530002"
    --echodiff="NONE"
    --unwarpdir="j"
    --fmrires="2"
    --dcmethod="SiemensFieldMap"
    --biascorrection="NONE"
    --gdcoeffs="NONE"
    --mctype="MCFLIRT"
    --processing-mode="HCPStyleData"
    --doslicetime="FALSE"
------------------------------------------------------------

---> ERROR: something missing, skipping this BOLD!

HCP fMRIVolume completed on Monday, 17. March 2025 07:09:56
------------------------------------------------------------

------------------------------------------------------------
Session id: 31-JA-JNLA-0002
[started on Monday, 17. March 2025 07:09:56]
Running HCP fMRI Volume pipeline [HCPStyleData] ...
---> PreFS results present.
---> FS results present.
---> PostFS results present.

---> Preprocessing settings (unwarpdir, refimage, moveref, seimage) for BOLD 1
     ... phase encoding direction: PA
     ... unwarp direction: j
     ... using image specific EchoSpacing: 0.000530002 s
---> ERROR: No fieldmap number specified for the BOLD image!
     ... bold image present
     ... reference image not used
     ... using the HCPpipelines default BOLD mask

Running 1 BOLD images in parallel

------------------------------------------------------------
Running HCP Pipelines command via QuNex:

/opt/HCP/HCPpipelines/fMRIVolume/GenericfMRIVolumeProcessingPipeline.sh
    --path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0002/hcp"
    --session="31-JA-JNLA-0002"
    --fmriname="1"
    --fmritcs="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0002/hcp/31-JA-JNLA-0002/unprocessed/BOLD_1_PA/31-JA-JNLA-0002_BOLD_1_PA.nii.gz"
    --fmriscout="NONE"
    --fmapmag="NONE"
    --fmapphase="NONE"
    --fmapcombined="NONE"
    --echospacing="0.000530002"
    --echodiff="NONE"
    --unwarpdir="j"
    --fmrires="2"
    --dcmethod="SiemensFieldMap"
    --biascorrection="NONE"
    --gdcoeffs="NONE"
    --mctype="MCFLIRT"
    --processing-mode="HCPStyleData"
    --doslicetime="FALSE"
------------------------------------------------------------

---> ERROR: something missing, skipping this BOLD!

HCP fMRIVolume completed on Monday, 17. March 2025 07:09:56
------------------------------------------------------------

------------------------------------------------------------
Session id: 31-JA-JNLA-0004
[started on Monday, 17. March 2025 07:09:56]
Running HCP fMRI Volume pipeline [HCPStyleData] ...
---> PreFS results present.
---> FS results present.
---> PostFS results present.

---> Preprocessing settings (unwarpdir, refimage, moveref, seimage) for BOLD 1
     ... phase encoding direction: PA
     ... unwarp direction: j
     ... using image specific EchoSpacing: 0.000530002 s
---> ERROR: No fieldmap number specified for the BOLD image!
     ... bold image present
     ... reference image not used
     ... using the HCPpipelines default BOLD mask

Running 1 BOLD images in parallel

------------------------------------------------------------
Running HCP Pipelines command via QuNex:

/opt/HCP/HCPpipelines/fMRIVolume/GenericfMRIVolumeProcessingPipeline.sh
    --path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0004/hcp"
    --session="31-JA-JNLA-0004"
    --fmriname="1"
    --fmritcs="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0004/hcp/31-JA-JNLA-0004/unprocessed/BOLD_1_PA/31-JA-JNLA-0004_BOLD_1_PA.nii.gz"
    --fmriscout="NONE"
    --fmapmag="NONE"
    --fmapphase="NONE"
    --fmapcombined="NONE"
    --echospacing="0.000530002"
    --echodiff="NONE"
    --unwarpdir="j"
    --fmrires="2"
    --dcmethod="SiemensFieldMap"
    --biascorrection="NONE"
    --gdcoeffs="NONE"
    --mctype="MCFLIRT"
    --processing-mode="HCPStyleData"
    --doslicetime="FALSE"
------------------------------------------------------------

---> ERROR: something missing, skipping this BOLD!

HCP fMRIVolume completed on Monday, 17. March 2025 07:09:56
------------------------------------------------------------


---> Final report for command hcp_fmri_volume
... 31-JA-JNLA-0001 ---> HCP fMRI Volume: bolds 1 not ready
... 31-JA-JNLA-0003 ---> HCP fMRI Volume: bolds 1 not ready
... 31-JA-JNLA-0002 ---> HCP fMRI Volume: bolds 1 not ready
... 31-JA-JNLA-0004 ---> HCP fMRI Volume: bolds 1 not ready
---> Not all tasks completed fully!

So I changed your suggested Mag_Images => FM-Magnitude to 221 => FM-Magnitude and after this change ERROR: No fieldmap number specified for the BOLD image! disappeared but other errors appeared!

(base) jiayf@EyeOpener:~/qunex$ qunex_container hcp_fmri_volume    --sessionsfolder="${STUDY_FOLDER}/sessions"     --batchfile="${STUDY_FOLDER}/processing/batch.txt"  --hcp_bold_dcmethod="SiemensFieldMap"   --parelements=1  --container="${QUNEX_CONTAINER}"

---> QuNex will run the command over 4 sessions. It will utilize:

    Maximum sessions run in parallel for a job: 1.
    Maximum elements run in parallel for a session: 1.
    Up to 1 processes will be utilized for a job.

    Job #1 will run sessions: 31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004
(base) jiayf@EyeOpener:~/qunex$ ---> unsetting the following environment variables: PATH MATLABPATH PYTHONPATH QUNEXVer TOOLS QUNEXREPO QUNEXPATH QUNEXEXTENSIONS QUNEXLIBRARY QUNEXLIBRARYETC TemplateFolder FSL_FIXDIR FREESURFERDIR FREESURFER_HOME FREESURFER_SCHEDULER FreeSurferSchedulerDIR WORKBENCHDIR DCMNIIDIR DICMNIIDIR MATLABDIR MATLABBINDIR OCTAVEDIR OCTAVEPKGDIR OCTAVEBINDIR RDIR HCPWBDIR AFNIDIR PYLIBDIR FSLDIR FSLBINDIR PALMDIR QUNEXMCOMMAND HCPPIPEDIR CARET7DIR GRADUNWARPDIR HCPPIPEDIR_Templates HCPPIPEDIR_Bin HCPPIPEDIR_Config HCPPIPEDIR_PreFS HCPPIPEDIR_FS HCPPIPEDIR_FS_CUSTOM HCPPIPEDIR_PostFS HCPPIPEDIR_fMRISurf HCPPIPEDIR_fMRIVol HCPPIPEDIR_tfMRI HCPPIPEDIR_dMRI HCPPIPEDIR_dMRITract HCPPIPEDIR_Global HCPPIPEDIR_tfMRIAnalysis HCPCIFTIRWDIR MSMBin HCPPIPEDIR_dMRITractFull HCPPIPEDIR_dMRILegacy AutoPtxFolder EDDYCUDA USEOCTAVE QUNEXENV CONDADIR MSMBINDIR MSMCONFIGDIR R_LIBS FSL_FIX_CIFTIRW FSFAST_HOME SUBJECTS_DIR MINC_BIN_DIR MNI_DIR MINC_LIB_DIR MNI_DATAPATH FSF_OUTPUT_FORMAT ANTSDIR CUDIMOT

========================================================================
Generated by QuNex
------------------------------------------------------------------------
Version: 1.1.0 [QIO]
User: root
System: cc3e1ddbd0da
OS: Debian Linux #1 SMP Tue Nov 5 00:21:55 UTC 2024
------------------------------------------------------------------------

        ██████\                  ║      ██\   ██\
       ██  __██\                 ║      ███\  ██ |
       ██ /  ██ |██\   ██\       ║      ████\ ██ | ██████\ ██\   ██\
       ██ |  ██ |██ |  ██ |      ║      ██ ██\██ |██  __██\\██\ ██  |
       ██ |  ██ |██ |  ██ |      ║      ██ \████ |████████ |\████  /
       ██ ██\██ |██ |  ██ |      ║      ██ |\███ |██   ____|██  ██\
       \██████ / \██████  |      ║      ██ | \██ |\███████\██  /\██\
        \___███\  \______/       ║      \__|  \__| \_______\__/  \__|
            \___|                ║


                       DEVELOPED & MAINTAINED BY:

               Mind & Brain Lab, University of Ljubljana
                       Cho Lab, Yale University

                      COPYRIGHT & LICENSE NOTICE:

Use of this software is subject to the terms and conditions defined in
QuNex LICENSES which can be found in the LICENSES folder of the QuNex
repository or at https://qunex.yale.edu/qunex-registration
========================================================================

---> Setting up Octave


.......................... Running QuNex v1.1.0 [QIO] ..........................


--- Full QuNex call for command: hcp_fmri_volume

qunex hcp_fmri_volume --sessionsfolder="/home/jiayf/qunex/parallel/sessions" --hcp_bold_dcmethod="SiemensFieldMap" --parelements="1" --batchfile="/home/jiayf/qunex/parallel/processing/batch.txt" --sessions="31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004"

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


# Generated by QuNex 1.1.0 [QIO] on 2025-03-17_07.11.59.372562#
=================================================================
qunex hcp_fmri_volume \
  --sessionsfolder="/home/jiayf/qunex/parallel/sessions" \
  --hcp_bold_dcmethod="SiemensFieldMap" \
  --parelements="1" \
  --sessions="/home/jiayf/qunex/parallel/processing/batch.txt" \
  --sessionids="31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004" \
=================================================================

Starting multiprocessing sessions in /home/jiayf/qunex/parallel/processing/batch.txt with a pool of 4 concurrent processes


Adding processing of session 31-JA-JNLA-0001 to the pool at Monday, 17. March 2025 07:11:59

Adding processing of session 31-JA-JNLA-0002 to the pool at Monday, 17. March 2025 07:11:59

Adding processing of session 31-JA-JNLA-0003 to the pool at Monday, 17. March 2025 07:11:59

Adding processing of session 31-JA-JNLA-0004 to the pool at Monday, 17. March 2025 07:11:59
=======================================================================================================================================
=======================================================================================================================================

------------------------------------------------------------
Session id: 31-JA-JNLA-0001
[started on Monday, 17. March 2025 07:11:59]
Running HCP fMRI Volume pipeline [HCPStyleData] ...
---> PreFS results present.
---> FS results present.
---> PostFS results present.

---> Preprocessing settings (unwarpdir, refimage, moveref, seimage) for BOLD 1
     ... phase encoding direction: PA
     ... unwarp direction: j
     ... using image specific EchoSpacing: 0.000530002 s
     ... Siemens fieldmap magnitude image 1 present
     ... Siemens fieldmap phase image 1 present
     ... ERROR: hcp_bold_echodiff not defined correctly: "NONE"!
     ... bold image present
     ... reference image not used
     ... using the HCPpipelines default BOLD mask

Running 1 BOLD images in parallel
---> ERROR: something missing, skipping this BOLD!

HCP fMRIVolume completed on Monday, 17. March 2025 07:11:59
------------------------------------------------------------
=======================================================================================================================================

------------------------------------------------------------
Session id: 31-JA-JNLA-0002
[started on Monday, 17. March 2025 07:11:59]
Running HCP fMRI Volume pipeline [HCPStyleData] ...
---> PreFS results present.
---> FS results present.
---> PostFS results present.

---> Preprocessing settings (unwarpdir, refimage, moveref, seimage) for BOLD 1
     ... phase encoding direction: PA
     ... unwarp direction: j
     ... using image specific EchoSpacing: 0.000530002 s
     ... Siemens fieldmap magnitude image 1 present
     ... Siemens fieldmap phase image 1 present
     ... ERROR: hcp_bold_echodiff not defined correctly: "NONE"!
     ... bold image present
     ... reference image not used
     ... using the HCPpipelines default BOLD mask

Running 1 BOLD images in parallel
---> ERROR: something missing, skipping this BOLD!

HCP fMRIVolume completed on Monday, 17. March 2025 07:11:59
------------------------------------------------------------
=======================================================================================================================================

------------------------------------------------------------
Session id: 31-JA-JNLA-0003
[started on Monday, 17. March 2025 07:11:59]
Running HCP fMRI Volume pipeline [HCPStyleData] ...
---> PreFS results present.
---> FS results present.
---> PostFS results present.

---> Preprocessing settings (unwarpdir, refimage, moveref, seimage) for BOLD 1
     ... phase encoding direction: PA
     ... unwarp direction: j
     ... using image specific EchoSpacing: 0.000530002 s
     ... Siemens fieldmap magnitude image 1 present
     ... Siemens fieldmap phase image 1 present
     ... ERROR: hcp_bold_echodiff not defined correctly: "NONE"!
     ... bold image present
     ... reference image not used
     ... using the HCPpipelines default BOLD mask

Running 1 BOLD images in parallel
---> ERROR: something missing, skipping this BOLD!

HCP fMRIVolume completed on Monday, 17. March 2025 07:11:59
------------------------------------------------------------

------------------------------------------------------------
Session id: 31-JA-JNLA-0004
[started on Monday, 17. March 2025 07:11:59]
Running HCP fMRI Volume pipeline [HCPStyleData] ...
---> PreFS results present.
---> FS results present.
---> PostFS results present.

---> Preprocessing settings (unwarpdir, refimage, moveref, seimage) for BOLD 1
     ... phase encoding direction: PA
     ... unwarp direction: j
     ... using image specific EchoSpacing: 0.000530002 s
     ... Siemens fieldmap magnitude image 1 present
     ... Siemens fieldmap phase image 1 present
     ... ERROR: hcp_bold_echodiff not defined correctly: "NONE"!
     ... bold image present
     ... reference image not used
     ... using the HCPpipelines default BOLD mask

Running 1 BOLD images in parallel
---> ERROR: something missing, skipping this BOLD!

HCP fMRIVolume completed on Monday, 17. March 2025 07:11:59
------------------------------------------------------------


---> Final report for command hcp_fmri_volume
... 31-JA-JNLA-0001 ---> HCP fMRI Volume: bolds 1 not ready
... 31-JA-JNLA-0002 ---> HCP fMRI Volume: bolds 1 not ready
... 31-JA-JNLA-0003 ---> HCP fMRI Volume: bolds 1 not ready
... 31-JA-JNLA-0004 ---> HCP fMRI Volume: bolds 1 not ready
---> Not all tasks completed fully!

This time it was reporting the error ERROR: hcp_bold_echodiff not defined correctly: “NONE”!, so I consulted MAG_IMAGES_0022.json and PHA_IMAGES_0023_ph.json. There are two parameters in PHA_IMAGES_0023_ph.json: “EchoTime1”: 0.02, “EchoTime2”: 0.015, and there is one parameter in MAG_IMAGES_0022.json: “EchoTime”: 0.02, so I calculated and made hcp_ bold_echodiff to be 20-15=5 (ms). Now hcp_fmri_volume seems to work, expect it to turn out correctly.

These are the problems I encountered during the tuning process and the corresponding solutions, so I hope I didn’t make any mistakes in my parameter modifications. In the meantime, I have followed your advice and re-included the dMRI data, even though it only has 64 directions (which does not meet the direction requirement of the HCPpipeline), and I have also named it DWI:dir98_AP, I don’t know if the result of this process is feasible.

Mr. Jure, please give me strength and give me some good advice! :heartbeat:

Best,
YJia

Hi,

OK, there is a lot to unpack here. :slight_smile:

  1. You are correct, setup_hcp does not have overwrite as the situation there is a bit more complex (setup_hcp — QuNex documentation). Instead of --overwrite=yes, you should set --existing=backup. This is the safest option as it creates a backup of the current hcp folder before creating anything, so even if something goes wrong, you have the old state saved. Without completed setup_hcp, the files are perhaps not setup correctly, and your processing might still not be optimal. So best to have every command complete successfully before moving on.

  2. Again speaking a bit from the top of my head. The logic for automatically assigning fieldmaps to BOLDs is that the last acquired fieldmap before BOLD is assigned to it. As your fieldmap was acquired after BOLD, I think you might have to do this manually. The webpage I linked above shows you how do to this. I think in your case this is the case:

#  parallel mapping file
#  ----------------

41  => T1w
51  => T2w
121 => bold:rest : fm(1)
161 => DWI:dir98_AP
221 => FM-Magnitude : fm(1)
231 => FM-Phase : fm(1)

Note that mapping from names is usually more robust. For example, 41 might not always be your T1w, but ABI1_t1iso_mprage will be. Not sure what the issue before is, I see there is a weird danlging ph in there, try:

Mag_Images    => FM-Magnitude
Pha_Images ph => FM-Phase
  1. ... ERROR: hcp_bold_echodiff not defined correctly: "NONE". QuNex tries to automatically calculate this based on provided JSON sidecars. Because setup_hcp failed, it might be that it cannot find them. I also see that you have _hcp_bold_echodiff : NONE in your batch file, remove this from it and from the parameters template that you used. There is a slight bug for this parameter and “NONE” right now, it will be fixed in the next version.

  2. Sice you do not have an AP/PA diffusion pair, you will not be able to use the hcp_diffusion pipelines. For such cases we have our own diffusion preprocessing pipeline, see dwi_legacy_gpu — QuNex documentation for details.

Best, Jure

Hi, :blush:

Thanks for the explanation, I think I understand the use of the setup_hcp command better.

Regarding the second point you mentioned - As your fieldmap was acquired after BOLD, I think you might have to do this manually, what does do this manually mean here? Does it mean that I need to fill in the mapping of fieldmap (FM-Magnitude and FM-Phase) in mapping.txt before each processing? Also the : fm(1) in 121 => bold:rest : fm(1) is just to mark the relationship between BOLD and fieldmap, right? There should be no need to include the : fm(1) tag in the official mapping.txt, right?

As for your point about using names as mapping relationships being more robust, I’ll study your suggestion! You also mentioned a weird danlging ph in mapping.txt, I checked session.txt and the raw data naming. The original DICOM folder for the sequence is named PHA_IMAGES_0023 and after conversion to NIFTI format the name is PHA_IMAGES_0023_ph.nii. The session.txt is as follows:

41  : ABI1_t1iso_mprage         : TR(2.53): DwellTime(8.2e-06)
51  : t2_tse_tra_320_p2         : TR(4.21): PEDirection(i): DwellTime(7.1e-06)
121 : ABI1_bold_rest            : TR(2): PEDirection(j): EchoSpacing(0.000530002): DwellTime(3.3e-06)
161 : ABI1_diff_hardi           : TR(13.7): PEDirection(j): EchoSpacing(0.000344999): DwellTime(2.6e-06)
221 : Mag_Images                : TR(0.028): DwellTime(1.3e-05)
231 : Pha_Images ph             : TR(0.028): DwellTime(1.3e-05)

You also mentioned the _hcp_bold_echodiff : NONE setting. In fact, there were other errors reported before that one, and it was only after resolving the others and re-completing the create_session_info, setup_hcp, create_batch command that I encountered the ERROR: hcp_bold_echodiff not defined correctly: ” NONE", for which I manually checked the batch.txt and manually calculated and modified the echodiff based on the information I consulted (this process was mentioned at the end of the previous reply, I’m not sure if the values are correct). After this modification, hcp_fmri_volume works fine and it took 50min to complete this step. But you mentioned remove this from it and from the parameters template that you used, do you mean I don’t need to use this parameter and I should remove it?

Lastly, about dMRI processing. I currently have two kinds of dMRI data in my hands, one has only one DICOM folder named DTI_b0_b1k_iso2mm_SaveBySlc_601, and the other one has more data related to the dMRI, and its folder naming screenshot is as follows:

The relevant information for these sequences is as follows, you can see from this that there is a lot of dMRI related data, I’m not sure if I have the AP/PA diffusion pair that you mentioned, if not I’ll take your advice and try to use the dwi_legacy_gpu command.

11  :                :ABI1_localizer [1/3] i00001: TR(0.0086): PEDirection(i): DwellTime(6.1e-06)
12  :                :ABI1_localizer [2/3] i00002: TR(0.0086): PEDirection(i): DwellTime(6.1e-06)
13  :                :ABI1_localizer [3/3] i00003: TR(0.0086): PEDirection(j-): DwellTime(6.1e-06)
21  :                :ABI1_localizer2 [1/3] i00001: TR(0.0086): PEDirection(i): DwellTime(6.1e-06)
22  :                :ABI1_localizer2 [2/3] i00002: TR(0.0086): PEDirection(i): DwellTime(6.1e-06)
23  :                :ABI1_localizer2 [3/3] i00003: TR(0.0086): PEDirection(j-): DwellTime(6.1e-06)
31  :                :t1_tirm_sag: TR(1.23): PEDirection(i): DwellTime(6e-06)
41  :T1w             :ABI1_t1iso_mprage: TR(2.53): DwellTime(8.2e-06)
51  :T2w             :t2_tse_tra_320_p2: TR(4.21): PEDirection(i): DwellTime(7.1e-06)
61  :                :t1_tirm_tra: TR(1.53): PEDirection(i): DwellTime(6e-06)
71  :                :ep2d_diff_3scan_trace_p2: TR(16.265): PEDirection(j-): EchoSpacing(0.000545015): DwellTime(3.1e-06)
81  :                :ep2d_diff_3scan_trace_p2_ADC: TR(5.4): PEDirection(j-): EchoSpacing(0.000545015): DwellTime(3.1e-06)
91  :                :t2_tirm_tra_dark-fluid: TR(5): PEDirection(i): DwellTime(6.9e-06)
101 :                :ABI1_fieldmap_rest [1/2] e1: TR(0.4): PEDirection(j): DwellTime(2.69e-05)
102 :                :ABI1_fieldmap_rest [2/2] e2: TR(0.4): PEDirection(j): DwellTime(2.69e-05)
111 :                :ABI1_fieldmap_rest e2_ph: TR(0.4): PEDirection(j): DwellTime(2.69e-05)
121 :bold1:rest      :ABI1_bold_rest: TR(2): PEDirection(j): EchoSpacing(0.000530002): DwellTime(3.3e-06)
131 :                :ABI1_fieldmap_hardi [1/2] e1: TR(0.727): PEDirection(j): DwellTime(7.9e-06)
132 :                :ABI1_fieldmap_hardi [2/2] e2: TR(0.727): PEDirection(j): DwellTime(7.9e-06)
141 :                :ABI1_fieldmap_hardi e2_ph: TR(0.727): PEDirection(j): DwellTime(7.9e-06)
151 :                :ABI1_diff_hardi_b0: TR(13.755): PEDirection(j): EchoSpacing(0.000475): DwellTime(2.6e-06)
161 :                :ABI1_diff_hardi: TR(13.7): PEDirection(j): EchoSpacing(0.000344999): DwellTime(2.6e-06)
171 :                :ABI1_diff_hardi_ADC: TR(13.7): PEDirection(j): EchoSpacing(0.000344999): DwellTime(2.6e-06)
181 :                :ABI1_diff_hardi_TRACEW: TR(13.7): PEDirection(j): EchoSpacing(0.000344999): DwellTime(2.6e-06)
191 :                :ABI1_diff_hardi_FA: TR(13.7): PEDirection(j): EchoSpacing(0.000344999): DwellTime(2.6e-06)
201 :                :ABI1_diff_hardi_ColFA: EchoSpacing(0.000344999)
221 :                :Mag_Images: TR(0.028): DwellTime(1.3e-05)
231 :                :Pha_Images ph: TR(0.028): DwellTime(1.3e-05)
241 :                :mIP_Images(SW): TR(0.028): DwellTime(1.3e-05)
251 :                :SWI_Images: TR(0.028): DwellTime(1.3e-05)
261 :                :ep2d_tra_pasl_p2: TR(3.5): PEDirection(j): EchoSpacing(0.000265001): DwellTime(3.4e-06)
271 :                :Perfusion_Weighted: TR(3.5): PEDirection(j): EchoSpacing(0.000265001): DwellTime(3.4e-06)
281 :                :relCBF: TR(3.5): PEDirection(j): EchoSpacing(0.000265001): DwellTime(3.4e-06)

Thanks again for your patient response!
Best, YJia

Hi, Jure

I was using dwi_legacy_gpu to process the test data at hand and then tried to process it using hcp_diffusion and it ran into reporting errors. From the error report, it was observed that its problem occurs here.

—> The following pos direction files were found: EMPTY
----> The following neg direction files were found:
31-JA-JNLA-0001_DWI_dir98_AP.nii.gz
----> Using image specific EchoSpacing: 0.000344999 s

By looking at the terminal output of other posters, I have a rough idea of what you mean by AP/PA diffusion pair.

—> The following pos direction files were found:
HCPA001_DWI_dir98_PA.nii.gz
HCPA001_DWI_dir99_PA.nii.gz
—> The following neg direction files were found:
HCPA001_DWI_dir98_AP.nii.gz
HCPA001_DWI_dir99_AP.nii.gz
—> Using image specific EchoSpacing: 0.000689998 s

If that’s the case, then both of my dMRI data may not qualify for hcp_diffusion. Maybe I do have to use dwi_legacy_gpu only for the data I have in hand.

(base) jiayf@EyeOpener:~/qunex$ qunex_container hcp_diffusion \
>  --sessionsfolder="${STUDY_FOLDER}/sessions" \
>  --batchfile="${STUDY_FOLDER}/processing/batch.txt" \
> --overwrite="yes" \
> --container="${QUNEX_CONTAINER}" \
> --hcp_nogpu

---> QuNex will run the command over 4 sessions. It will utilize:

    Maximum sessions run in parallel for a job: 1.
    Maximum elements run in parallel for a session: 1.
    Up to 1 processes will be utilized for a job.

    Job #1 will run sessions: 31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004
(base) jiayf@EyeOpener:~/qunex$ ---> unsetting the following environment variables: PATH MATLABPATH PYTHONPATH QUNEXVer TOOLS QUNEXREPO QUNEXPATH QUNEXEXTENSIONS QUNEXLIBRARY QUNEXLIBRARYETC TemplateFolder FSL_FIXDIR FREESURFERDIR FREESURFER_HOME FREESURFER_SCHEDULER FreeSurferSchedulerDIR WORKBENCHDIR DCMNIIDIR DICMNIIDIR MATLABDIR MATLABBINDIR OCTAVEDIR OCTAVEPKGDIR OCTAVEBINDIR RDIR HCPWBDIR AFNIDIR PYLIBDIR FSLDIR FSLBINDIR PALMDIR QUNEXMCOMMAND HCPPIPEDIR CARET7DIR GRADUNWARPDIR HCPPIPEDIR_Templates HCPPIPEDIR_Bin HCPPIPEDIR_Config HCPPIPEDIR_PreFS HCPPIPEDIR_FS HCPPIPEDIR_FS_CUSTOM HCPPIPEDIR_PostFS HCPPIPEDIR_fMRISurf HCPPIPEDIR_fMRIVol HCPPIPEDIR_tfMRI HCPPIPEDIR_dMRI HCPPIPEDIR_dMRITract HCPPIPEDIR_Global HCPPIPEDIR_tfMRIAnalysis HCPCIFTIRWDIR MSMBin HCPPIPEDIR_dMRITractFull HCPPIPEDIR_dMRILegacy AutoPtxFolder EDDYCUDA USEOCTAVE QUNEXENV CONDADIR MSMBINDIR MSMCONFIGDIR R_LIBS FSL_FIX_CIFTIRW FSFAST_HOME SUBJECTS_DIR MINC_BIN_DIR MNI_DIR MINC_LIB_DIR MNI_DATAPATH FSF_OUTPUT_FORMAT ANTSDIR CUDIMOT

========================================================================
Generated by QuNex
------------------------------------------------------------------------
Version: 1.1.0 [QIO]
User: root
System: 8aa6afb6439e
OS: Debian Linux #1 SMP Tue Nov 5 00:21:55 UTC 2024
------------------------------------------------------------------------

        ██████\                  ║      ██\   ██\
       ██  __██\                 ║      ███\  ██ |
       ██ /  ██ |██\   ██\       ║      ████\ ██ | ██████\ ██\   ██\
       ██ |  ██ |██ |  ██ |      ║      ██ ██\██ |██  __██\\██\ ██  |
       ██ |  ██ |██ |  ██ |      ║      ██ \████ |████████ |\████  /
       ██ ██\██ |██ |  ██ |      ║      ██ |\███ |██   ____|██  ██\
       \██████ / \██████  |      ║      ██ | \██ |\███████\██  /\██\
        \___███\  \______/       ║      \__|  \__| \_______\__/  \__|
            \___|                ║


                       DEVELOPED & MAINTAINED BY:

               Mind & Brain Lab, University of Ljubljana
                       Cho Lab, Yale University

                      COPYRIGHT & LICENSE NOTICE:

Use of this software is subject to the terms and conditions defined in
QuNex LICENSES which can be found in the LICENSES folder of the QuNex
repository or at https://qunex.yale.edu/qunex-registration
========================================================================

---> Setting up Octave


.......................... Running QuNex v1.1.0 [QIO] ..........................


--- Full QuNex call for command: hcp_diffusion

qunex hcp_diffusion --sessionsfolder="/home/jiayf/qunex/parallel/sessions" --overwrite="yes" --hcp_nogpu --batchfile="/home/jiayf/qunex/parallel/processing/batch.txt" --sessions="31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004"

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


# Generated by QuNex 1.1.0 [QIO] on 2025-03-18_01.05.36.758297#
=================================================================
qunex hcp_diffusion \
  --sessionsfolder="/home/jiayf/qunex/parallel/sessions" \
  --overwrite="yes" \
  --hcp_nogpu="True" \
  --sessions="/home/jiayf/qunex/parallel/processing/batch.txt" \
  --sessionids="31-JA-JNLA-0001,31-JA-JNLA-0002,31-JA-JNLA-0003,31-JA-JNLA-0004" \
=================================================================

Starting multiprocessing sessions in /home/jiayf/qunex/parallel/processing/batch.txt with a pool of 4 concurrent processes


Adding processing of session 31-JA-JNLA-0001 to the pool at Tuesday, 18. March 2025 01:05:36

Adding processing of session 31-JA-JNLA-0002 to the pool at Tuesday, 18. March 2025 01:05:36

Adding processing of session 31-JA-JNLA-0003 to the pool at Tuesday, 18. March 2025 01:05:36

Adding processing of session 31-JA-JNLA-0004 to the pool at Tuesday, 18. March 2025 01:05:36
Running external command: /opt/HCP/HCPpipelines/DiffusionPreprocessing/DiffPreprocPipeline.sh
--path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0004/hcp"                 --session="31-JA-JNLA-0004"                 --PEdir=2                 --posData="EMPTY"                 --negData="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0004/hcp/31-JA-JNLA-0004/unprocessed/Diffusion/31-JA-JNLA-0004_DWI_dir98_AP.nii.gz"                 --echospacing-seconds="0.000344999"                 --gdcoeffs="NONE"                 --combine-data-flag="1"                 --printcom=""                --gpu=False

You can follow command's progress in:
/home/jiayf/qunex/parallel/processing/logs/comlogs/tmp_hcp_diffusion_31-JA-JNLA-0004_2025-03-18_01.05.36.774928.log
------------------------------------------------------------
Running external command: /opt/HCP/HCPpipelines/DiffusionPreprocessing/DiffPreprocPipeline.sh
--path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0003/hcp"                 --session="31-JA-JNLA-0003"                 --PEdir=2                 --posData="EMPTY"                 --negData="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0003/hcp/31-JA-JNLA-0003/unprocessed/Diffusion/31-JA-JNLA-0003_DWI_dir98_AP.nii.gz"                 --echospacing-seconds="0.000344999"                 --gdcoeffs="NONE"                 --combine-data-flag="1"                 --printcom=""                --gpu=False

You can follow command's progress in:
/home/jiayf/qunex/parallel/processing/logs/comlogs/tmp_hcp_diffusion_31-JA-JNLA-0003_2025-03-18_01.05.36.775014.log
------------------------------------------------------------
Running external command: /opt/HCP/HCPpipelines/DiffusionPreprocessing/DiffPreprocPipeline.sh
--path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/hcp"                 --session="31-JA-JNLA-0001"                 --PEdir=2                 --posData="EMPTY"                 --negData="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/hcp/31-JA-JNLA-0001/unprocessed/Diffusion/31-JA-JNLA-0001_DWI_dir98_AP.nii.gz"                 --echospacing-seconds="0.000344999"                 --gdcoeffs="NONE"                 --combine-data-flag="1"                 --printcom=""                --gpu=False

You can follow command's progress in:
/home/jiayf/qunex/parallel/processing/logs/comlogs/tmp_hcp_diffusion_31-JA-JNLA-0001_2025-03-18_01.05.36.775107.log
------------------------------------------------------------
Running external command: /opt/HCP/HCPpipelines/DiffusionPreprocessing/DiffPreprocPipeline.sh
--path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0002/hcp"                 --session="31-JA-JNLA-0002"                 --PEdir=2                 --posData="EMPTY"                 --negData="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0002/hcp/31-JA-JNLA-0002/unprocessed/Diffusion/31-JA-JNLA-0002_DWI_dir98_AP.nii.gz"                 --echospacing-seconds="0.000344999"                 --gdcoeffs="NONE"                 --combine-data-flag="1"                 --printcom=""                --gpu=False

You can follow command's progress in:
/home/jiayf/qunex/parallel/processing/logs/comlogs/tmp_hcp_diffusion_31-JA-JNLA-0002_2025-03-18_01.05.36.777324.log
------------------------------------------------------------

------------------------------------------------------------
Session id: 31-JA-JNLA-0004
[started on Tuesday, 18. March 2025 01:05:36]
Running HCP DiffusionPreprocessing Pipeline [HCPStyleData] ...
---> The following pos direction files were found:
     EMPTY
---> The following neg direction files were found:
     31-JA-JNLA-0004_DWI_dir98_AP.nii.gz
---> Using image specific EchoSpacing: 0.000344999 s

------------------------------------------------------------
Running HCP Pipelines command via QuNex:

/opt/HCP/HCPpipelines/DiffusionPreprocessing/DiffPreprocPipeline.sh
    --path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0004/hcp"
    --session="31-JA-JNLA-0004"
    --PEdir=2
    --posData="EMPTY"
    --negData="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0004/hcp/31-JA-JNLA-0004/unprocessed/Diffusion/31-JA-JNLA-0004_DWI_dir98_AP.nii.gz"
    --echospacing-seconds="0.000344999"
    --gdcoeffs="NONE"
    --combine-data-flag="1"
    --printcom=""
    --gpu=False
------------------------------------------------------------


Running HCP Diffusion Preprocessing

ERROR: Running HCP Diffusion Preprocessing failed with error 1
...
command executed:
/opt/HCP/HCPpipelines/DiffusionPreprocessing/DiffPreprocPipeline.sh                 --path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0004/hcp"                 --session="31-JA-JNLA-0004"                 --PEdir=2                 --posData="EMPTY"                 --negData="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0004/hcp/31-JA-JNLA-0004/unprocessed/Diffusion/31-JA-JNLA-0004_DWI_dir98_AP.nii.gz"                 --echospacing-seconds="0.000344999"                 --gdcoeffs="NONE"                 --combine-data-flag="1"
   --printcom=""                --gpu=False

---> logfile: /home/jiayf/qunex/parallel/processing/logs/comlogs/error_hcp_diffusion_31-JA-JNLA-0004_2025-03-18_01.05.36.774928.log

HCP Diffusion Preprocessing completed on Tuesday, 18. March 2025 01:05:39
------------------------------------------------------------

------------------------------------------------------------
Session id: 31-JA-JNLA-0002
[started on Tuesday, 18. March 2025 01:05:36]
Running HCP DiffusionPreprocessing Pipeline [HCPStyleData] ...
---> The following pos direction files were found:
     EMPTY
---> The following neg direction files were found:
     31-JA-JNLA-0002_DWI_dir98_AP.nii.gz
---> Using image specific EchoSpacing: 0.000344999 s

------------------------------------------------------------
Running HCP Pipelines command via QuNex:

/opt/HCP/HCPpipelines/DiffusionPreprocessing/DiffPreprocPipeline.sh
    --path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0002/hcp"
    --session="31-JA-JNLA-0002"
    --PEdir=2
    --posData="EMPTY"
    --negData="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0002/hcp/31-JA-JNLA-0002/unprocessed/Diffusion/31-JA-JNLA-0002_DWI_dir98_AP.nii.gz"
    --echospacing-seconds="0.000344999"
    --gdcoeffs="NONE"
    --combine-data-flag="1"
    --printcom=""
    --gpu=False
------------------------------------------------------------


Running HCP Diffusion Preprocessing

ERROR: Running HCP Diffusion Preprocessing failed with error 1
...
command executed:
/opt/HCP/HCPpipelines/DiffusionPreprocessing/DiffPreprocPipeline.sh                 --path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0002/hcp"                 --session="31-JA-JNLA-0002"                 --PEdir=2                 --posData="EMPTY"                 --negData="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0002/hcp/31-JA-JNLA-0002/unprocessed/Diffusion/31-JA-JNLA-0002_DWI_dir98_AP.nii.gz"                 --echospacing-seconds="0.000344999"                 --gdcoeffs="NONE"                 --combine-data-flag="1"
   --printcom=""                --gpu=False

---> logfile: /home/jiayf/qunex/parallel/processing/logs/comlogs/error_hcp_diffusion_31-JA-JNLA-0002_2025-03-18_01.05.36.777324.log

HCP Diffusion Preprocessing completed on Tuesday, 18. March 2025 01:05:39
------------------------------------------------------------

------------------------------------------------------------
Session id: 31-JA-JNLA-0001
[started on Tuesday, 18. March 2025 01:05:36]
Running HCP DiffusionPreprocessing Pipeline [HCPStyleData] ...
---> The following pos direction files were found:
     EMPTY
---> The following neg direction files were found:
     31-JA-JNLA-0001_DWI_dir98_AP.nii.gz
---> Using image specific EchoSpacing: 0.000344999 s

------------------------------------------------------------
Running HCP Pipelines command via QuNex:

/opt/HCP/HCPpipelines/DiffusionPreprocessing/DiffPreprocPipeline.sh
    --path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/hcp"
    --session="31-JA-JNLA-0001"
    --PEdir=2
    --posData="EMPTY"
    --negData="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/hcp/31-JA-JNLA-0001/unprocessed/Diffusion/31-JA-JNLA-0001_DWI_dir98_AP.nii.gz"
    --echospacing-seconds="0.000344999"
    --gdcoeffs="NONE"
    --combine-data-flag="1"
    --printcom=""
    --gpu=False
------------------------------------------------------------


Running HCP Diffusion Preprocessing

ERROR: Running HCP Diffusion Preprocessing failed with error 1
...
command executed:
/opt/HCP/HCPpipelines/DiffusionPreprocessing/DiffPreprocPipeline.sh                 --path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/hcp"                 --session="31-JA-JNLA-0001"                 --PEdir=2                 --posData="EMPTY"                 --negData="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/hcp/31-JA-JNLA-0001/unprocessed/Diffusion/31-JA-JNLA-0001_DWI_dir98_AP.nii.gz"                 --echospacing-seconds="0.000344999"                 --gdcoeffs="NONE"                 --combine-data-flag="1"
   --printcom=""                --gpu=False

---> logfile: /home/jiayf/qunex/parallel/processing/logs/comlogs/error_hcp_diffusion_31-JA-JNLA-0001_2025-03-18_01.05.36.775107.log

HCP Diffusion Preprocessing completed on Tuesday, 18. March 2025 01:05:39
------------------------------------------------------------

------------------------------------------------------------
Session id: 31-JA-JNLA-0003
[started on Tuesday, 18. March 2025 01:05:36]
Running HCP DiffusionPreprocessing Pipeline [HCPStyleData] ...
---> The following pos direction files were found:
     EMPTY
---> The following neg direction files were found:
     31-JA-JNLA-0003_DWI_dir98_AP.nii.gz
---> Using image specific EchoSpacing: 0.000344999 s

------------------------------------------------------------
Running HCP Pipelines command via QuNex:

/opt/HCP/HCPpipelines/DiffusionPreprocessing/DiffPreprocPipeline.sh
    --path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0003/hcp"
    --session="31-JA-JNLA-0003"
    --PEdir=2
    --posData="EMPTY"
    --negData="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0003/hcp/31-JA-JNLA-0003/unprocessed/Diffusion/31-JA-JNLA-0003_DWI_dir98_AP.nii.gz"
    --echospacing-seconds="0.000344999"
    --gdcoeffs="NONE"
    --combine-data-flag="1"
    --printcom=""
    --gpu=False
------------------------------------------------------------


Running HCP Diffusion Preprocessing

ERROR: Running HCP Diffusion Preprocessing failed with error 1
...
command executed:
/opt/HCP/HCPpipelines/DiffusionPreprocessing/DiffPreprocPipeline.sh                 --path="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0003/hcp"                 --session="31-JA-JNLA-0003"                 --PEdir=2                 --posData="EMPTY"                 --negData="/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0003/hcp/31-JA-JNLA-0003/unprocessed/Diffusion/31-JA-JNLA-0003_DWI_dir98_AP.nii.gz"                 --echospacing-seconds="0.000344999"                 --gdcoeffs="NONE"                 --combine-data-flag="1"
   --printcom=""                --gpu=False

---> logfile: /home/jiayf/qunex/parallel/processing/logs/comlogs/error_hcp_diffusion_31-JA-JNLA-0003_2025-03-18_01.05.36.775014.log

HCP Diffusion Preprocessing completed on Tuesday, 18. March 2025 01:05:39
------------------------------------------------------------


---> Final report for command hcp_diffusion
... 31-JA-JNLA-0004 ---> Error
... 31-JA-JNLA-0002 ---> Error
... 31-JA-JNLA-0001 ---> Error
... 31-JA-JNLA-0003 ---> Error
---> Not all tasks completed fully!

Best, YJia

Hi, Jure

After initially completing the three modal preprocessing of my test data, I need its QC result, where should I get it?

Should I use run_qc to get the relevant QC pictures?

Or should I get it from sessions\31-JA-JNLA-0001\hcp\31-JA-JNLA-0001\MNINonLinear\StructuralQC\snapshots and sessions\31-JA-JNLA-0001\hcp\31-JA-JNLA- 0001\MNINonLinear\Results\1\fMRIQC\snapshots in obtaining relevant quality control images?

Or should I get the relevant QC pictures from sessions\QC? (This folder is empty)

Or should I not check the effects from these QC results plots and should I check the processed .nii data directly?

Best, YJia

Hi,

1. Volume processing

Based on the information I have, I think you are not using fieldmaps for hcp_fmri_volume which will lead to suboptimal results. To use fieldmaps you need to execute the whole create_session_info => create_batch => setup_hcp chain correctly. Fieldmaps need bo assigned to other modalities (BOLD, struct …). In your case they are not. The traditional acquisition is to acquire fieldmap which is then immediately followed by say a BOLD scan. Here, QuNex will automatically link a field map preceding the BOLD scan. In your case the fieldmap was acquired after the BOLD, which is not the most common way, as a result there are no fieldmaps preceding the BOLD and QuNex will not assign the fieldmap to your BOLD. Because of this you need to do this assignment manually in the mapping file. Try:

ABI1_t1iso_mprage  => T1w
t2_tse_tra_320_p2  => T2w
ABI1_bold_rest     => bold:rest : fm(1)
ABI1_diff_hardi    => DWI:dir98_AP
Mag_Images         => FM-Magnitude
Pha_Images ph      => FM-Phase

Whether fieldmaps were assigned to bolds can be seen from either a file called session_hcp.txt or the batch.txt file. Since you have a single BOLD and a single fieldmap the BOLD should be linked with fm(1), the BOLD will have this tag in its line.

2. QC

HCP prepares some of its own QC outputs (StructuralQC, fMRIQC). QuNex has a specialized run_qc command that works across multiple modalities and generates two types of outputs. PNG images for a quick review, and Connectome Workbench scenes for a deeper inspection when the PNGs are not enough. This is what we use on our end. The folder sessions\QC is empty because you did not execute run_qc, this is the default location where the results will be stored.

3. Diffusion

Unfortunately, I am not a sequence expert, so I cannot really tell you what all the DIFF folders are. If I had to take I guess, I would say that ABI1_diff_hardi is the “real” diffusion data. So, you do not have an AP/PA pair, meaning that you cannot use hcp_diffusion and any other diffusion things you might have (b0). If this is indeed the case, you need to use dwi_legacy_gpu.

Best, Jure

Thanks, Jure

Actually I executed the whole create_session_info => create_batch => setup_hcp process, maybe the versions of mapping.txt and batch.txt I provided you before were not correct (I saved several versions of the process results folder). The current session_hcp.txt I’m using is as follows, and I notice that these mappings are followed by fm(1) links, does this mean that fieldmaps are being used properly?

hcpready: true

41  :T1w             :ABI1_t1iso_mprage: TR(2.53): DwellTime(8.2e-06): fm(1)
51  :T2w             :t2_tse_tra_320_p2: TR(4.21): PEDirection(i): DwellTime(7.1e-06): fm(1)
121 :bold1:rest      :ABI1_bold_rest: TR(2): PEDirection(j): EchoSpacing(0.000530002): DwellTime(3.3e-06): fm(1)
161 :DWI:dir98_AP    :ABI1_diff_hardi: TR(13.7): PEDirection(j): EchoSpacing(0.000344999): DwellTime(2.6e-06): fm(1)
221 :FM-Magnitude    :Mag_Images: TR(0.028): DwellTime(1.3e-05): fm(1)
231 :FM-Phase        :Pha_Images ph: TR(0.028): DwellTime(1.3e-05): fm(1)

Regarding QC, I first looked at HCPpipeline’s own output of QC result images, and then later by checking the forums using the run_qc command to get QC result images for the modal and raw data used. If I subsequently need to use Connectome Workbench to view the detail information of the .scene file, do I need to download and use Connectome Workbench separately to view it? (QuNex doesn’t seem to have a command to support viewing).

Regarding the Diffusion data, I did preprocess it using ABI1_diff_hardi, and I took your earlier suggestion of using the dwi_legacy_gpu command instead of hcp_diffusion (although I named ABI1_diff_hardi as DWI:dir98_AP)

Here is a graph of my BOLD, T1w and DWI partial QC results after executing the run_qc command.



Right now I’m referencing https://forum.qunex.yale.edu/t/resolved-setting-up-connectivity-analysis/1233/18 and within the forums https://qunex.readthedocs.io/en/ latest/wiki/UsageDocs/BOLDPreprocessing.html for other processing.

Best, YJia

Hi,

In terms of processing everything looks good (session_hcp.txt) looks correct, and the fieldmaps are properly linked.

QuNex is traditionally used on compute systems or clusters and does not have a graphical interface. Even if it had it, you usually cannot use a GUI on such systems. How we do QC is that we first review the run_qc PNGs. Based on this review, the sessions that look good pass the QC step. For sessions that might have some issues, we download the zipped workbench scenes that QuNex prepared and review them in detail on our personal computers with Connectome Workbench (Connectome - Using Connectome Workbench). We decide whether a session is OK or not based on this review. Furthermore, this is usually done by 2 persons in parallel and then they discuss sessions where their QC conclusions differ.

Best, Jure

Hi,

I completed the BOLD process by referring to this post, which is below:
BOLD task-evoked and resting-state functional connectivity preprocessing — QuNex documentation

map_hcp_data
run_QC
create_bold_brain_masks
compute_bold_stats
create_stats_report
extract_nuisance_signal
preprocess_bold

After completing the above steps, I would like to use the fc_compute_roifc command for connectivity analysis. I referenced this post during that time:
[已解决]设置连通性分析 - 使用问题 - QuNex Discourse — [RESOLVED] Setting up connectivity analysis - Use Questions - QuNex Discourse

I followed the post and used the map2hcp command to create an AAL mapping in CIFTI format. But I still get the following error when running fc_compute_roifc:

(base) jiayf@EyeOpener:~/qunex$ qunex_container fc_compute_roifc --sessionsfolder="${STUDY_FOLDER}/sessions" --flist="listname:31-JA-JNLA-0001|session id:31-JA-JNLA-0001|file:/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/images/functional/bold1_Atlas_s_hpss_res-mVWMWB1d.dtseries.nii" --roiinfo="/home/jiayf/qunex/AAL3v1.names" --targetf="/home/jiayf/qunex" --options="sessions=all|roimethod=mean|eventdata=all|ignore=use,fidl|badevents=use|fcmeasure=r|debug=false|saveind=mat,long,wide_single|itargetf:gfolder" --frames_info="5" --container="${QUNEX_CONTAINER}"
(base) jiayf@EyeOpener:~/qunex$ ---> unsetting the following environment variables: PATH MATLABPATH PYTHONPATH QUNEXVer TOOLS QUNEXREPO QUNEXPATH QUNEXEXTENSIONS QUNEXLIBRARY QUNEXLIBRARYETC TemplateFolder FSL_FIXDIR FREESURFERDIR FREESURFER_HOME FREESURFER_SCHEDULER FreeSurferSchedulerDIR WORKBENCHDIR DCMNIIDIR DICMNIIDIR MATLABDIR MATLABBINDIR OCTAVEDIR OCTAVEPKGDIR OCTAVEBINDIR RDIR HCPWBDIR AFNIDIR PYLIBDIR FSLDIR FSLBINDIR PALMDIR QUNEXMCOMMAND HCPPIPEDIR CARET7DIR GRADUNWARPDIR HCPPIPEDIR_Templates HCPPIPEDIR_Bin HCPPIPEDIR_Config HCPPIPEDIR_PreFS HCPPIPEDIR_FS HCPPIPEDIR_FS_CUSTOM HCPPIPEDIR_PostFS HCPPIPEDIR_fMRISurf HCPPIPEDIR_fMRIVol HCPPIPEDIR_tfMRI HCPPIPEDIR_dMRI HCPPIPEDIR_dMRITract HCPPIPEDIR_Global HCPPIPEDIR_tfMRIAnalysis HCPCIFTIRWDIR MSMBin HCPPIPEDIR_dMRITractFull HCPPIPEDIR_dMRILegacy AutoPtxFolder EDDYCUDA USEOCTAVE QUNEXENV CONDADIR MSMBINDIR MSMCONFIGDIR R_LIBS FSL_FIX_CIFTIRW FSFAST_HOME SUBJECTS_DIR MINC_BIN_DIR MNI_DIR MINC_LIB_DIR MNI_DATAPATH FSF_OUTPUT_FORMAT ANTSDIR CUDIMOT

========================================================================
Generated by QuNex
------------------------------------------------------------------------
Version: 1.1.0 [QIO]
User: root
System: e3bc6bc0ec08
OS: Debian Linux #1 SMP Tue Nov 5 00:21:55 UTC 2024
------------------------------------------------------------------------

        ██████\                  ║      ██\   ██\
       ██  __██\                 ║      ███\  ██ |
       ██ /  ██ |██\   ██\       ║      ████\ ██ | ██████\ ██\   ██\
       ██ |  ██ |██ |  ██ |      ║      ██ ██\██ |██  __██\\██\ ██  |
       ██ |  ██ |██ |  ██ |      ║      ██ \████ |████████ |\████  /
       ██ ██\██ |██ |  ██ |      ║      ██ |\███ |██   ____|██  ██\
       \██████ / \██████  |      ║      ██ | \██ |\███████\██  /\██\
        \___███\  \______/       ║      \__|  \__| \_______\__/  \__|
            \___|                ║


                       DEVELOPED & MAINTAINED BY:

               Mind & Brain Lab, University of Ljubljana
                       Cho Lab, Yale University

                      COPYRIGHT & LICENSE NOTICE:

Use of this software is subject to the terms and conditions defined in
QuNex LICENSES which can be found in the LICENSES folder of the QuNex
repository or at https://qunex.yale.edu/qunex-registration
========================================================================

---> Setting up Octave


.......................... Running QuNex v1.1.0 [QIO] ..........................


--- Full QuNex call for command: fc_compute_roifc

qunex fc_compute_roifc --sessionsfolder="/home/jiayf/qunex/parallel/sessions" --flist="listname:31-JA-JNLA-0001|session id:31-JA-JNLA-0001|file:/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/images/functional/bold1_Atlas_s_hpss_res-mVWMWB1d.dtseries.nii" --roiinfo="/home/jiayf/qunex/AAL3v1.names" --targetf="/home/jiayf/qunex" --options="sessions=all|roimethod=mean|eventdata=all|ignore=use,fidl|badevents=use|fcmeasure=r|debug=false|saveind=mat,long,wide_single|itargetf:gfolder" --frames_info="5"

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



Running:
>>> fc_compute_roifc('listname:31-JA-JNLA-0001|session id:31-JA-JNLA-0001|file:/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/images/functional/bold1_Atlas_s_hpss_res-mVWMWB1d.dtseries.nii', '/home/jiayf/qunex/AAL3v1.names', '', '/home/jiayf/qunex', 'sessions=all|roimethod=mean|eventdata=all|ignore=use,fidl|badevents=use|fcmeasure=r|debug=false|saveind=mat,long,wide_single|itargetf:gfolder')


Checking ...
 ... listing files to process
 ... reading file list:
     - session id: 31-JA-JNLA-0001
       ... found image file (/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/images/functional/bold1_Atlas_s_hpss_res-mVWMWB1d.dtseries.nii)

 ... done ... done.

---------------------------------
Processing session 31-JA-JNLA-0001     ... reading image file(s) ... 240 frames read, done.
     ... creating ROI mask

Matlab Error! Processing Failed!

ERROR: The roi file [/home/jiayf/qunex/Precentral_L|1|] specified in the .names file [/home/jiayf/qunex/AAL3v1.names] is not a valid image file!


ERROR: fc_compute_roifc failed! Please check output / log!

It looks like my roiinfo file is incorrect, but I don’t think it’s any different than the aal.names in the post.

Precentral_L|1|
Precentral_R|2|
Frontal_Sup_2_L|3|
Frontal_Sup_2_R|4|
Frontal_Mid_2_L|5|
Frontal_Mid_2_R|6|
Frontal_Inf_Oper_L|7|
Frontal_Inf_Oper_R|8|
Frontal_Inf_Tri_L|9|
Frontal_Inf_Tri_R|10|
Frontal_Inf_Orb_2_L|11|
Frontal_Inf_Orb_2_R|12|
Rolandic_Oper_L|13|
Rolandic_Oper_R|14|
Supp_Motor_Area_L|15|
Supp_Motor_Area_R|16|
Olfactory_L|17|
Olfactory_R|18|
Frontal_Sup_Medial_L|19|
Frontal_Sup_Medial_R|20|
Frontal_Med_Orb_L|21|
Frontal_Med_Orb_R|22|
Rectus_L|23|
Rectus_R|24|
OFCmed_L|25|
OFCmed_R|26|
OFCant_L|27|
OFCant_R|28|
OFCpost_L|29|
OFCpost_R|30|
OFClat_L|31|
OFClat_R|32|
Insula_L|33|
Insula_R|34|
Cingulate_Ant_L|35|
Cingulate_Ant_R|36|
Cingulate_Mid_L|37|
Cingulate_Mid_R|38|
Cingulate_Post_L|39|
Cingulate_Post_R|40|
Hippocampus_L|41|
Hippocampus_R|42|
ParaHippocampal_L|43|
ParaHippocampal_R|44|
Amygdala_L|45|
Amygdala_R|46|
Calcarine_L|47|
Calcarine_R|48|
Cuneus_L|49|
Cuneus_R|50|
Lingual_L|51|
Lingual_R|52|
Occipital_Sup_L|53|
Occipital_Sup_R|54|
Occipital_Mid_L|55|
Occipital_Mid_R|56|
Occipital_Inf_L|57|
Occipital_Inf_R|58|
Fusiform_L|59|
Fusiform_R|60|
Postcentral_L|61|
Postcentral_R|62|
Parietal_Sup_L|63|
Parietal_Sup_R|64|
Parietal_Inf_L|65|
Parietal_Inf_R|66|
SupraMarginal_L|67|
SupraMarginal_R|68|
Angular_L|69|
Angular_R|70|
Precuneus_L|71|
Precuneus_R|72|
Paracentral_Lobule_L|73|
Paracentral_Lobule_R|74|
Caudate_L|75|
Caudate_R|76|
Putamen_L|77|
Putamen_R|78|
Pallidum_L|79|
Pallidum_R|80|
Thalamus_L|81|
Thalamus_R|82|
Heschl_L|83|
Heschl_R|84|
Temporal_Sup_L|85|
Temporal_Sup_R|86|
Temporal_Pole_Sup_L|87|
Temporal_Pole_Sup_R|88|
Temporal_Mid_L|89|
Temporal_Mid_R|90|
Temporal_Pole_Mid_L|91|
Temporal_Pole_Mid_R|92|
Temporal_Inf_L|93|
Temporal_Inf_R|94|
Cerebellum_Crus1_L|95|
Cerebellum_Crus1_R|96|
Cerebellum_Crus2_L|97|
Cerebellum_Crus2_R|98|
Cerebellum_3_L|99|
Cerebellum_3_R|100|
Cerebellum_4_5_L|101|
Cerebellum_4_5_R|102|
Cerebellum_6_L|103|
Cerebellum_6_R|104|
Cerebellum_7b_L|105|
Cerebellum_7b_R|106|
Cerebellum_8_L|107|
Cerebellum_8_R|108|
Cerebellum_9_L|109|
Cerebellum_9_R|110|
Cerebellum_10_L|111|
Cerebellum_10_R|112|
Vermis_1_2|113|
Vermis_3|114|
Vermis_4_5|115|
Vermis_6|116|
Vermis_7|117|
Vermis_8|118|
Vermis_9|119|
Vermis_10|120|
Thal_AV_L|121|
Thal_AV_R|122|
Thal_LP_L|123|
Thal_LP_R|124|
Thal_VA_L|125|
Thal_VA_R|126|
Thal_VL_L|127|
Thal_VL_R|128|
Thal_VPL_L|129|
Thal_VPL_R|130|
Thal_IL_L|131|
Thal_IL_R|132|
Thal_Re_L|133|
Thal_Re_R|134|
Thal_MDm_L|135|
Thal_MDm_R|136|
Thal_MDl_L|137|
Thal_MDl_R|138|
Thal_LGN_L|139|
Thal_LGN_R|140|
Thal_MGN_L|141|
Thal_MGN_R|142|
Thal_PuI_L|143|
Thal_PuI_R|144|
Thal_PuM_L|145|
Thal_PuM_R|146|
Thal_PuA_L|147|
Thal_PuA_R|148|
Thal_PuL_L|149|
Thal_PuL_R|150|
ACC_sub_L|151|
ACC_sub_R|152|
ACC_pre_L|153|
ACC_pre_R|154|
ACC_sup_L|155|
ACC_sup_R|156|
N_Acc_L|157|
N_Acc_R|158|
VTA_L|159|
VTA_R|160|
SN_pc_L|161|
SN_pc_R|162|
SN_pr_L|163|
SN_pr_R|164|
Red_N_L|165|
Red_N_R|166|
LC_L|167|
LC_R|168|
Raphe_D|169|
Raphe_M|170|

Finally, I am a little confused about the BOLD task-evoked and resting-state functional connectivity preprocessing. From my cursory understanding, it seems that the fmrivolume step in the HCPpipeline has already done some denoising, so why is there a need for additional preprocessing of the BOLD data through multiple steps before the resting-state functional connectivity analysis? (e.g. extract_nuisance_signal etc.)

Best, YJia

Your roiinfo file is indeed incorrect, if you take a look at the forum example you pasted, you can see that the first row needs to be the path to the .dscalar.nii file.

Best, Jure

Hi, Jure

You’re right! I tried again, it works now! Thanks for your generous help!

Here are the contents of my AAL3v1.names and the output of the fc_compute_roifc command, hopefully the stupid mistakes I made will help others!

/home/jiayf/qunex/AAL3v1.dscalar.nii
Precentral_L|1|
Precentral_R|2|
Frontal_Sup_2_L|3|
Frontal_Sup_2_R|4|
Frontal_Mid_2_L|5|
Frontal_Mid_2_R|6|
Frontal_Inf_Oper_L|7|
Frontal_Inf_Oper_R|8|
Frontal_Inf_Tri_L|9|
Frontal_Inf_Tri_R|10|
Frontal_Inf_Orb_2_L|11|
Frontal_Inf_Orb_2_R|12|
Rolandic_Oper_L|13|
Rolandic_Oper_R|14|
Supp_Motor_Area_L|15|
Supp_Motor_Area_R|16|
Olfactory_L|17|
Olfactory_R|18|
Frontal_Sup_Medial_L|19|
Frontal_Sup_Medial_R|20|
Frontal_Med_Orb_L|21|
Frontal_Med_Orb_R|22|
Rectus_L|23|
Rectus_R|24|
OFCmed_L|25|
OFCmed_R|26|
OFCant_L|27|
OFCant_R|28|
OFCpost_L|29|
OFCpost_R|30|
OFClat_L|31|
OFClat_R|32|
Insula_L|33|
Insula_R|34|
Cingulate_Ant_L|35|
Cingulate_Ant_R|36|
Cingulate_Mid_L|37|
Cingulate_Mid_R|38|
Cingulate_Post_L|39|
Cingulate_Post_R|40|
Hippocampus_L|41|
Hippocampus_R|42|
ParaHippocampal_L|43|
ParaHippocampal_R|44|
Amygdala_L|45|
Amygdala_R|46|
Calcarine_L|47|
Calcarine_R|48|
Cuneus_L|49|
Cuneus_R|50|
Lingual_L|51|
Lingual_R|52|
Occipital_Sup_L|53|
Occipital_Sup_R|54|
Occipital_Mid_L|55|
Occipital_Mid_R|56|
Occipital_Inf_L|57|
Occipital_Inf_R|58|
Fusiform_L|59|
Fusiform_R|60|
Postcentral_L|61|
Postcentral_R|62|
Parietal_Sup_L|63|
Parietal_Sup_R|64|
Parietal_Inf_L|65|
Parietal_Inf_R|66|
SupraMarginal_L|67|
SupraMarginal_R|68|
Angular_L|69|
Angular_R|70|
Precuneus_L|71|
Precuneus_R|72|
Paracentral_Lobule_L|73|
Paracentral_Lobule_R|74|
Caudate_L|75|
Caudate_R|76|
Putamen_L|77|
Putamen_R|78|
Pallidum_L|79|
Pallidum_R|80|
Thalamus_L|81|
Thalamus_R|82|
Heschl_L|83|
Heschl_R|84|
Temporal_Sup_L|85|
Temporal_Sup_R|86|
Temporal_Pole_Sup_L|87|
Temporal_Pole_Sup_R|88|
Temporal_Mid_L|89|
Temporal_Mid_R|90|
Temporal_Pole_Mid_L|91|
Temporal_Pole_Mid_R|92|
Temporal_Inf_L|93|
Temporal_Inf_R|94|
Cerebellum_Crus1_L|95|
Cerebellum_Crus1_R|96|
Cerebellum_Crus2_L|97|
Cerebellum_Crus2_R|98|
Cerebellum_3_L|99|
Cerebellum_3_R|100|
Cerebellum_4_5_L|101|
Cerebellum_4_5_R|102|
Cerebellum_6_L|103|
Cerebellum_6_R|104|
Cerebellum_7b_L|105|
Cerebellum_7b_R|106|
Cerebellum_8_L|107|
Cerebellum_8_R|108|
Cerebellum_9_L|109|
Cerebellum_9_R|110|
Cerebellum_10_L|111|
Cerebellum_10_R|112|
Vermis_1_2|113|
Vermis_3|114|
Vermis_4_5|115|
Vermis_6|116|
Vermis_7|117|
Vermis_8|118|
Vermis_9|119|
Vermis_10|120|
Thal_AV_L|121|
Thal_AV_R|122|
Thal_LP_L|123|
Thal_LP_R|124|
Thal_VA_L|125|
Thal_VA_R|126|
Thal_VL_L|127|
Thal_VL_R|128|
Thal_VPL_L|129|
Thal_VPL_R|130|
Thal_IL_L|131|
Thal_IL_R|132|
Thal_Re_L|133|
Thal_Re_R|134|
Thal_MDm_L|135|
Thal_MDm_R|136|
Thal_MDl_L|137|
Thal_MDl_R|138|
Thal_LGN_L|139|
Thal_LGN_R|140|
Thal_MGN_L|141|
Thal_MGN_R|142|
Thal_PuI_L|143|
Thal_PuI_R|144|
Thal_PuM_L|145|
Thal_PuM_R|146|
Thal_PuA_L|147|
Thal_PuA_R|148|
Thal_PuL_L|149|
Thal_PuL_R|150|
ACC_sub_L|151|
ACC_sub_R|152|
ACC_pre_L|153|
ACC_pre_R|154|
ACC_sup_L|155|
ACC_sup_R|156|
N_Acc_L|157|
N_Acc_R|158|
VTA_L|159|
VTA_R|160|
SN_pc_L|161|
SN_pc_R|162|
SN_pr_L|163|
SN_pr_R|164|
Red_N_L|165|
Red_N_R|166|
LC_L|167|
LC_R|168|
Raphe_D|169|
Raphe_M|170|
(base) jiayf@EyeOpener:~/qunex$ qunex_container fc_compute_roifc --sessionsfolder="${STUDY_FOLDER}/sessions" --flist="listname:31-JA-JNLA-0001|session id:31-JA-JNLA-0001|file:/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/images/functional/bold1_Atlas_s_hpss_res-mVWMWB1d.dtseries.nii" --roiinfo="/home/jiayf/qunex/AAL3v1.names" --targetf="/home/jiayf/qunex" --options="sessions=all|roimethod=mean|eventdata=all|ignore=use,fidl|badevents=use|fcmeasure=r|debug=false|saveind=mat,long,wide_single|itargetf:gfolder" --frames_info="5" --container="${QUNEX_CONTAINER}"
(base) jiayf@EyeOpener:~/qunex$ ---> unsetting the following environment variables: PATH MATLABPATH PYTHONPATH QUNEXVer TOOLS QUNEXREPO QUNEXPATH QUNEXEXTENSIONS QUNEXLIBRARY QUNEXLIBRARYETC TemplateFolder FSL_FIXDIR FREESURFERDIR FREESURFER_HOME FREESURFER_SCHEDULER FreeSurferSchedulerDIR WORKBENCHDIR DCMNIIDIR DICMNIIDIR MATLABDIR MATLABBINDIR OCTAVEDIR OCTAVEPKGDIR OCTAVEBINDIR RDIR HCPWBDIR AFNIDIR PYLIBDIR FSLDIR FSLBINDIR PALMDIR QUNEXMCOMMAND HCPPIPEDIR CARET7DIR GRADUNWARPDIR HCPPIPEDIR_Templates HCPPIPEDIR_Bin HCPPIPEDIR_Config HCPPIPEDIR_PreFS HCPPIPEDIR_FS HCPPIPEDIR_FS_CUSTOM HCPPIPEDIR_PostFS HCPPIPEDIR_fMRISurf HCPPIPEDIR_fMRIVol HCPPIPEDIR_tfMRI HCPPIPEDIR_dMRI HCPPIPEDIR_dMRITract HCPPIPEDIR_Global HCPPIPEDIR_tfMRIAnalysis HCPCIFTIRWDIR MSMBin HCPPIPEDIR_dMRITractFull HCPPIPEDIR_dMRILegacy AutoPtxFolder EDDYCUDA USEOCTAVE QUNEXENV CONDADIR MSMBINDIR MSMCONFIGDIR R_LIBS FSL_FIX_CIFTIRW FSFAST_HOME SUBJECTS_DIR MINC_BIN_DIR MNI_DIR MINC_LIB_DIR MNI_DATAPATH FSF_OUTPUT_FORMAT ANTSDIR CUDIMOT

========================================================================
Generated by QuNex
------------------------------------------------------------------------
Version: 1.1.0 [QIO]
User: root
System: 1c9a614548d6
OS: Debian Linux #1 SMP Tue Nov 5 00:21:55 UTC 2024
------------------------------------------------------------------------

        ██████\                  ║      ██\   ██\
       ██  __██\                 ║      ███\  ██ |
       ██ /  ██ |██\   ██\       ║      ████\ ██ | ██████\ ██\   ██\
       ██ |  ██ |██ |  ██ |      ║      ██ ██\██ |██  __██\\██\ ██  |
       ██ |  ██ |██ |  ██ |      ║      ██ \████ |████████ |\████  /
       ██ ██\██ |██ |  ██ |      ║      ██ |\███ |██   ____|██  ██\
       \██████ / \██████  |      ║      ██ | \██ |\███████\██  /\██\
        \___███\  \______/       ║      \__|  \__| \_______\__/  \__|
            \___|                ║


                       DEVELOPED & MAINTAINED BY:

               Mind & Brain Lab, University of Ljubljana
                       Cho Lab, Yale University

                      COPYRIGHT & LICENSE NOTICE:

Use of this software is subject to the terms and conditions defined in
QuNex LICENSES which can be found in the LICENSES folder of the QuNex
repository or at https://qunex.yale.edu/qunex-registration
========================================================================

---> Setting up Octave


.......................... Running QuNex v1.1.0 [QIO] ..........................


--- Full QuNex call for command: fc_compute_roifc

qunex fc_compute_roifc --sessionsfolder="/home/jiayf/qunex/parallel/sessions" --flist="listname:31-JA-JNLA-0001|session id:31-JA-JNLA-0001|file:/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/images/functional/bold1_Atlas_s_hpss_res-mVWMWB1d.dtseries.nii" --roiinfo="/home/jiayf/qunex/AAL3v1.names" --targetf="/home/jiayf/qunex" --options="sessions=all|roimethod=mean|eventdata=all|ignore=use,fidl|badevents=use|fcmeasure=r|debug=false|saveind=mat,long,wide_single|itargetf:gfolder" --frames_info="5"

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



Running:
>>> fc_compute_roifc('listname:31-JA-JNLA-0001|session id:31-JA-JNLA-0001|file:/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/images/functional/bold1_Atlas_s_hpss_res-mVWMWB1d.dtseries.nii', '/home/jiayf/qunex/AAL3v1.names', '', '/home/jiayf/qunex', 'sessions=all|roimethod=mean|eventdata=all|ignore=use,fidl|badevents=use|fcmeasure=r|debug=false|saveind=mat,long,wide_single|itargetf:gfolder')


Checking ...
 ... listing files to process
 ... reading file list:
     - session id: 31-JA-JNLA-0001
       ... found image file (/home/jiayf/qunex/parallel/sessions/31-JA-JNLA-0001/images/functional/bold1_Atlas_s_hpss_res-mVWMWB1d.dtseries.nii)

 ... done ... done.

---------------------------------
Processing session 31-JA-JNLA-0001     ... reading image file(s) ... 240 frames read, done.
     ... creating ROI mask
     ... generating extraction sets
         -> timeseries: 1 good events, 231 good frames
     ... computing fc matrices
         ... set timeseries
         ... extracted ts
         ... computed fc matrix
     ... saving results
         ... saving mat file ... done
         ... saving long tsv file ... done
         ... saving wide tsv file ... done


---> Successful completion of task

Best, YJia

Great! Good luck with your connectome analyses. :slight_smile:

Best, Jure