lesion_sci_t2

Intramedullary SCI lesion and cord segmentation in T2w MRI

This segmentation model for spinal cord injury segmentation uses a 3D U-Net architecture, and was trained with the nnUNetV2 framework. It is a multiclass model, outputting segmentations for both the hyperintense SCI lesions and spinal cord. Training data consisted of T2w images from 7 sites with traumatic (acute pre-operative, intermediate, chronic), non-traumatic (DCM) and ischemic SCI lesions spanning numerous resolutions, orientations, as well as multiple scanner manufacturers and field strengths.

Reference

@InProceedings{10.1007/978-3-031-82007-6_19,
               author="Karthik, Enamundram Naga and Valo{{s}}ek, Jan and Farner, Lynn and Pfyffer, Dario and Schading-Sassenhausen, Simon and Lebret, Anna and David, Gergely and Smith, Andrew C. and Weber II, Kenneth A. and Seif, Maryam and Freund, Patrick and Cohen-Adad, Julien",
               editor="Wu, Shandong and Shabestari, Behrouz and Xing, Lei",
               title="SCIsegV2: A Universal Tool for Segmentation of Intramedullary Lesions in Spinal Cord Injury",
               booktitle="Applications of Medical Artificial Intelligence",
               year="2025",
               publisher="Springer Nature Switzerland",
               address="Cham",
               pages="198--209",
               abstract="Spinal cord injury (SCI) is a devastating incidence leading to permanent paralysis and loss of sensory-motor functions potentially resulting in the formation of lesions within the spinal cord. Imaging biomarkers obtained from magnetic resonance imaging (MRI) scans can predict the functional recovery of individuals with SCI and help choose the optimal treatment strategy. Currently, most studies employ manual quantification of these MRI-derived biomarkers, which is a subjective and tedious task. In this work, we propose (i) a universal tool for the automatic segmentation of intramedullary SCI lesions, dubbed SCIsegV2, and (ii) a method to automatically compute the width of the tissue bridges from the segmented lesion. Tissue bridges represent the spared spinal tissue adjacent to the lesion, which is associated with functional recovery in SCI patients. The tool was trained and validated on a heterogeneous dataset from 7 sites comprising patients from different SCI phases (acute, sub-acute, and chronic) and etiologies (traumatic SCI, ischemic SCI, and degenerative cervical myelopathy). Tissue bridges quantified automatically did not significantly differ from those computed manually, suggesting that the proposed automatic tool can be used to derive relevant MRI biomarkers. SCIsegV2 and the automatic tissue bridges computation are open-source and available in Spinal Cord Toolbox (v6.4 and above) via the sct{\_}deepseg -task seg{\_}sc{\_}lesion{\_}t2w{\_}sci and sct{\_}analyze{\_}lesion functions, respectively.",
               isbn="978-3-031-82007-6"
}

Project URL: https://github.com/ivadomed/model_seg_sci

usage: sct_deepseg lesion_sci_t2 [-i <file> [<file> ...]] [-o <str>] [-install]
                                 [-custom-url CUSTOM_URL [CUSTOM_URL ...]]
                                 [-largest {0,1}] [-fill-holes {0,1}]
                                 [-remove-small REMOVE_SMALL [REMOVE_SMALL ...]]
                                 [-qc <folder>] [-qc-dataset <str>]
                                 [-qc-subject <str>] [-qc-plane <str>]
                                 [-qc-seg <file>] [-h] [-v <int>] [-r {0,1}]
                                 [-test-time-aug]

INPUT/OUTPUT

-i

Image filename(s) to segment. If segmenting multiple files, separate filenames with a space.

-o

Output file name. The chosen filename will be used as a base name, and model-specific suffixes will be added to the end depending on the type of output (e.g. ‘_cord.nii.gz’, ‘_gm.nii.gz’, etc.).

TASKS

-install

Install models that are required for specified task.

Default: False

-custom-url

URL(s) pointing to the .zip asset for a model release. This option can be used with -install to install a specific version of a model. To use this option, navigate to the ‘Releases’ page of the model, find release you wish to install, and right-click + copy the URL of the .zip listed under ‘Assets’. Example: sct_deepseg lesion_sci_t2 -install -custom-url CUSTOM_URL sct_deepseg lesion_sci_t2 -i t2.nii.gz

PARAMETERS

-largest

Possible choices: 0, 1

Keep the largest connected object from each output segmentation; if not set, all objects are kept.

Default: 0

-fill-holes

Possible choices: 0, 1

If set, small holes in the segmentation will be filled in automatically.

Default: 0

-remove-small

Minimal object size to keep with unit (mm3 or vox). A single value can be provided or one value per prediction class. Single value example: 1mm3, 5vox. Multiple values example: 10 20 10vox (remove objects smaller than 10 voxels for class 1 and 3, and smaller than 20 voxels for class 2).

-test-time-aug

Perform test-time augmentation (TTA) by flipping the input image along all axes and averaging the resulting predictions. Note: The time it takes to run the model will increase due to the additional predictions.

Default: False

MISC ARGUMENTS

-qc

The path where the quality control generated content will be saved.

-qc-dataset

If provided, this string will be mentioned in the QC report as the dataset the process was run on.

-qc-subject

If provided, this string will be mentioned in the QC report as the subject the process was run on.

-qc-plane

Possible choices: Axial, Sagittal

Plane of the output QC. If Sagittal, it is highly recommended to provide the -qc-seg option, as it will ensure the output QC is cropped to a reasonable field of view.

Default: 'Axial'

-qc-seg

Segmentation file to use for cropping the QC. This option is useful when you want to QC a region that is different from the output segmentation. For example, it might be useful to provide a dilated cord segmentation to expand the QC field of view.

If -qc-seg is not provided, the default behavior will depend on the value of -qc-plane:

• ‘Axial’: Without ‘-qc-seg’, a sensible crop radius between 15-40 vox will be automatically used, depending on the resolution and segmentation type.

• ‘Sagittal’: Without ‘-qc-seg’, the full image will be displayed by default. (For very large images, this may cause a crash, so using -qc-seg is highly recommended.)

-v

Possible choices: 0, 1, 2

Verbosity. 0: Display only errors/warnings, 1: Errors/warnings + info messages, 2: Debug mode.

Default: 1

-r

Possible choices: 0, 1

Remove temporary files.

Default: 1