Hi,
I have been given a set of DICOM data and I am using heudiconv
to convert them to NIfTI and BIDS. Unfortunately, I have very limited experience converting DICOM data.
I have a custon heuristic file that I have developped based on the reproin
example available on the heudiconv
repository. The data I have contains T1w data, dMRI data (AP and PA) and SE field maps.
I am using a singularity image built on Oct 29 from the latest version. My call to heudiconv
is simple:
singularity exec \
--bind ${data_dirname} \
${heudiconv_singularity_fname} \
heudiconv \
-f ${config_fname} \
--bids \
-o ${out_dirname} \
--files ${in_dirname}
where the variables have been appropriately set to the folders of interest.
The DICOM files are located in three separate folders, T1
, DWI_AP
, DWI_PA
, Spin_Echo_Maps
, each one containing one folder containing all the DICOM files per participant.
When running the conversion, the task raises an assertion error because when checking that the study session that is being parsed is not already in the list,
it finds that the session is already in the list of study sessions.
The data belongs to the same study (at least for my purposes), but the StudySessionInfo
instances show that they have no session
information, and the locator
(study) is different, e.g.
StudySessionInfo(locator='Investigators/MyStudySZPain10820', session=None, subject='092743')
(...)
StudySessionInfo(locator='Investigators/MyStudySZPain8232021', session=None, subject='197668')
(...)
There seems to be 4 different locators (MyStudy
is a substitute for the real name):
Investigators/MyStudy
Investigators/MyStudyFDNeuro08242021
Investigators/MyStudyFDPain08242021
Investigators/MyStudySZPain
Investigators/MyStudySZPain10820
Investigators/MyStudySZPain8232021
If I keep only the DWI folders in my in_dirname
, heudiconv
is able to complete the conversion, and organizes the data into an Investigators
main folder with the rest of MyStudy*
subfolders containing a subset of the 68 participants contained in the dataset.
If I keep the T1
and DWI_AP
folder in the in_dirname
, I get the error mentioned above.
I found this PR
ENH: grouping by mgxd · Pull Request #359 · nipy/heudiconv · GitHub
and thought that maybe the --grouping all
flag could be helpful to solve this. When using the flag I get an assertion error from my config file, because the seq infos are not unique (as said, inherited from the reproin
example):
I do not know how to go about this or how to debug this.
Any help is highly appreciated.
Thanks.