CSA Measurements and Lesion Measurements

Hello,

Three questions regarding the measurements by the SCT:

  1. We are generating csa.csv files using “sct_process_segmentation”. Can you please advice which of the columns are showing the spinal cord volume? Or what can be the equation to measure this?

  2. We are generating lesionseg_analyzis Excel files from Axial and Sagital images with an equal number of lesions, but the measurements are very different in volume - can you share, based on your experience the differences between the two?

  3. Is there a way to combine Sagital and Axial MRI images to improve the results? I tried converting Dicom to NII but the result was two files instead of a combined one (using dcim2niix).

Thank you!

hi,

  1. We are generating csa.csv files using “sct_process_segmentation”. Can you please advice which of the columns are showing the spinal cord volume? Or what can be the equation to measure this?

If you output the slicewise csv (i.e., one CSA per slice), you can compute the cord volume by summing the column “MEAN(area)” across the desired number of slices, and then multiply the result by the slice thickness.

  1. We are generating lesionseg_analyzis Excel files from Axial and Sagital images with an equal number of lesions, but the measurements are very different in volume - can you share, based on your experience the differences between the two?

The difference is related to the precision of your measure. With thick sagittal and axial slices, the delineation of the lesion is hampered by partial volume effect, so you will inevitably experience this loss of precision. One way to handle this is to compute lesion volume in both images and then average your results. For more details, see this article, which points to open-source code that deals with this situation.

  1. Is there a way to combine Sagital and Axial MRI images to improve the results? I tried converting Dicom to NII but the result was two files instead of a combined one (using dcim2niix).

See my answer above

best
Julien

Dear Julien,

Thank you so much for your answers, they were really helpful!

We are working now to produce segmentation of the white & gray matter with SCT 4.0.0.
We were able to segment the Gray Matter using: sct_deepseg_gm. Please see the following screenshot:
gm_1
Is it possible to measure (to Excel) the Gray Matter volume/area? We couldn’t use the tools mentioned in 2016 article since they have been changed.

In addition, we were unable to segment the White Matter - we couldn’t find an appropriate tool for that. Can you share the tools used, if any, to segment the White Matter as well?

Thank you very much,
Achiron

Hi,
To compute GM and WM area/volume, see the batch_processing.sh example file, as well as the SCT course (pages 80-84).
Here is an example syntax:

# Subtract GM segmentation from cord segmentation to obtain WM segmentation
sct_maths -i t2s_seg.nii.gz -sub t2s_gmseg.nii.gz -o t2s_wmseg.nii.gz
# Compute cross-sectional area of the gray and white matter between C2 and C5
sct_process_segmentation -i t2s_wmseg.nii.gz -vert 2:5 -perlevel 1 -o csa_wm.csv
sct_process_segmentation -i t2s_gmseg.nii.gz -vert 2:5 -perlevel 1 -o csa_gm.csv

Dear Julien
I hope you are well in these catastrophic times…
We have been working on our data of spinal cord MRI with the program. We found that if we repeat the analysis again we get different numbers.
Please advice

Anat Achiron

Hi Anat
I need a bit more specific information. Namely: what syntax and what SCT version produces variable results?

Dear Julien
It seems that the reason for the different results were related the implementation of white and grey matter measurements which we added.
When we re-tested the data results were the same!
Thanks for that.
The SCT version we use is 4.1.1

Anat

1 Like