Please send an email tp LINDAemail@example.com to register to the mailing list and receive notifications of new releases.
What is LINDA?
Is a neuroimaging toolkit for automatic segmentation of chronic stroke lesions. The method is described in Hum Brain Mapp. 2016 Apr;37(4):1405-21.
- Linux or Mac (since Oct 2016 ANTsR can run in Windows 10)
- R v3.0 or above
- The ANTsR package installed in R
- A T1-weighted MRI of a patient with (left hemispheric) stroke
Download the latest release (v0.2.6, 580Mb) and unzip in a local folder.
IMPORTANT: Only the above link contains the full release with trained models, templates, etc. Do not use the “Download” button on Github. It contains only smaller scripts because Github does not accept files above 100mb in main repositories.
Open R and source the file linda_predict.R
A file dialog will allow you to select the T1 nifti file of the patient with left hemispheric brain damage.
LINDA will run and you will see the timestamp of the various steps.
Results will be saved in a folder named “linda” in the same path where the T1 is located.
Note, LINDA will stop if you don’t have
ANTsR installed, and will automatically install the
Available prediction models:
Currently a model trained on 60 patients from Penn is used. The internal prediction engine works with 2mm voxel resolution. This does not mean your images need to be 2mm, any resolution should work.
14:47 Starting LINDA v0.2.6 …
14:47 Loading file: subjA_T1.nii
14:47 Creating folder: /mnt/c/Users/TESTsubject/linda
14:47 Loading template…
14:47 Skull stripping… (long process)
14:54 Saving skull stripped files
14:54 Loading LINDA model
14:54 Computing asymmetry mask…
14:55 Saving asymmetry mask…
14:55 Running 1st registration…
14:56 Feature calculation…
14:58 Running 1st prediction…
14:58 Saving 1st prediction…
14:58 Backprojecting 1st prediction…
14:58 Running 2nd registration…
15:00 Feature calculation…
15:01 Running 2nd prediction…
15:02 Saving 2nd prediction…
15:02 Backprojecting 2nd prediction…
15:02 Running 3rd registration… (long process)
15:35 Saving 3rd registration results…
15:36 Feature calculation…
15:37 Running 3rd prediction…
15:38 Saving 3rd final prediction…
15:38 Saving 3rd final prediction in template space…
15:38 Saving 3rd final prediction in native space…
15:38 Saving probabilistic prediction in template space…
15:38 Saving probabilistic prediction in native space…
15:38 Transferring data in MNI (ch2) space…
15:38 Saving subject in MNI (ch2) space…
15:38 Saving lesion in MNI (ch2) space…
BrainMask.nii.gz - Brain mask used to skull strip (native space)
N4corrected.nii.gz - Bias corrected, full image (native space)
N4corrected_Brain.nii.gz - Bias corrected, brain only (native space)
N4corrected_Brain_LRflipped.nii.gz - Flipped image used to compute asymmetry mask (native space)
Mask.lesion(1)(2)(3).nii.gz - masks used for registrations (native space)
Prediction(1)(2)(3).nii.gz - lesion predictions after each iteration (template space, but 2mm)
Prediction3template.nii.gz__ - final prediction (template space)
__Prediction3_native.nii.gz__ - final prediction (native space)
__Prediction3_probability_template.nii.gz__ - final graded probability (template space)
__Prediction3_probability_native.nii.gz__ - final graded probability (native space)
__Reg3_registered_to_template.nii.gz__ - Subject’s MRI, skull stripped, bias corrected, registered to template (template space)
__Reg3_sub_to_template(affine)(warp) - transformation matrices to register subject to template
Reg3template_to_sub(affine)(warp) - transformation matrices to backproject template to subject
Subject_in_MNI.nii.gz - Subject’s MRI, skull stripped, bias corrected, transformed in MNI space
Lesion_in_MNI.nii.gz - Lesion mask in MNI space
MNI is a space, not a template. There are many templates in MNI space, most of which do not have the same gyri or sulci at the same coordinates. LINDA uses the CH2 template which is included in LINDA. Please don’t use other MNI templates to overlay results. They may look good from a first look, but they will be wrong. Use the CH2 template that comes with LINDA. In alternative, you can register the CH2 template to another template (i.e., ICBM 2009a) and transform back and forth the results as necessary.
Frequently Asked Questions
- What file formats are accepted?
Nifti images (.nii and .nii.gz) are accepted. The platform can read many other formats, but we have limited the script to Nifti in order to avoid confusion with some formats (i.e., Analyze) that have unclear left-right orientation.
- Can I use it with right hemispheric lesions?
Yes, but you need to flip the L-R orientation before. After that, the prediction will work just as well. We plan to extend the script in the future for use in both left and right lesions.
- Can I use it with bilateral lesions?
It will likely be less accurate. One of the features we use is related to the left-right signal asymmetry
- Can I use it to segment acute and subacute stroke lesions?
No, the current model is trained only on chronic stroke patients. It might be possible to segment acute stroke with models trained on acute data (let us know if you want to contribute with those data).
- Can I use other images: FLAIR, T2, DWI?
The existing model uses only a T1, but we can adapt it with additional T2, FLAIR, DWI. We are open to collaborations with groups that have multimodal data and want to train LINDA with those. Having ~30 subjects is a good start.
- Can I train a model with my own data?
This is perfectly doable, but we haven’t made available the script online yet (needs some work to adapt it for widepsread use). If you are in a hurry, contact us and we will send the current version of the training script.
- Can I use LINDA to obtain registrations in MNI?
The transfer in MNI is obtained by concatenating transformations “Subject” -> “Penn Template” -> “ch2 MNI template”. Thus there are two sources of potential error. The second source of error is fixed for all subjects because our template has always the same registration on MNI. However, a fixed error is always an error. If you really want precise registration on MNI, we advise to register directly the subject to an MNI template (possibly using a group MNI template rather than the ch2).
The best way to keep track of bugs or failures is to open a New Issue on the Github system. You can also contact the author via email: dorian dot pustina at uphs dot upenn dot edu (translate from english). If the algorithm proceeds without errors but the automatic segmentation is erroneous, please send (i) your T1 image and (ii) the segmentation produced by LINDA in native space. Try also overlaying
Mask.lesion*.nii.gz files on the T1 to check whether the brain mask is wrong somewhere.