Weird outputs from the CAT segmentation process

Hi @Francois & all,
I have noticed some weird output from the CAT segmentation in some of my data,
here are some views of the outputs:


image

I got this in two dataset, the MRI looks fine but the final output is not.

I'm trying to find the issue, so any help is appreciated.
I'm sharing also one MRI where this issue is observed,

Thanks for you help :slight_smile:

Takfarinas

I could reproduce this behavior by running CAT12 from Brainstorm.

The execution report shows that there is something off with the definition of the volume coordinates (the correspondence between voxels and scanner coordinates).

The .nii file has both SFORM and QFORM matrices defined, and one is probably unusual or wrong.
image

@tmedani
Where is the volume coming from?
How was it pre-processed?
Can you try processing the volume with CAT12 without this extra pre-processing?

@CGaser
Do you have any idea of the cause of the problem?
How can we ensure a .nii file like this gets processed correctly by CAT?

Here is the full CAT output folder:
https://www.dropbox.com/s/s1e2t8pjnsmscux/cat12_badmri.zip?dl=0

1 Like

from one of our collaborator that works in hospital.

As far as I know they just anonymize the data and convert from dicom to nifti.
but Im going to check.

Yes, I will ask for the "original" data and I will let you know.

Thank you

@tmedani
Francois is absolutely right. There is something odd with your SFORM and QFORM matrices. If you visualize the image in SPM12 you see that there is some shearing applied and the voxel sizes are not correct. I have tried to reorient, but there the shearing couldn't be repaired. Please ask for the dicom images of that subject or whether the nifti files are converted using a different tool. This seems to be not nifti-conform even if it looks fine with other tools.

1 Like

@tmedani We don't see the shearing in Brainstorm, BrainSuite or MRIcron because all these tools display slices of voxels from the cube of values saved in the .nii file. The SFORM/QFORM transformation do not alter the visualization, they only change the computation of the coordinates.
At the contrary, SPM/CAT compute slices through the volume along arbitrary axes (not necessarily aligned with the voxel axes), and including the SFORM or QFORM transformation. Errors in these matrices become visible.

Please ask for the dicom images of that subject or whether the nifti files are converted using a different tool.

I hope the answer is not going to be that the conversion was done with Brainstorm... :slight_smile:

1 Like

Hello @Francois & @CGaser

Thanks again for your help.

This is an update regarding this issue.
I was able to get he original MRIs from the hospital without any processing and the CAT/SPM/SIMNIBS pipelines worked fine.

I want also to highlight that when I used Brainsuite to segment these "BadMRIs", I was able to get good outputs (correct shape and alignement). This mean that Brainsuite get/recompute the missing value from somewhere else at the NiftiFile.

The remaining issue with Brainsuite, in some cases, the WM and GM surface are intersecting, and thus its impossible to generate correct tetramesh with the current tools.

no, we are safe :slight_smile:

Thanks :slight_smile:

Or that BrainSuite ignores the SFORM/QFORM information for the segmentation.