spinalcordtoolbox / PAM50

https://github.com/neuropoly/spinalcordtoolbox
2 stars 1 forks source link

Linear scaling : potentially inaccurate positioning of individual spinal levels #10

Open sandrinebedard opened 1 year ago

sandrinebedard commented 1 year ago

Description

We are currently working to update the spinal levels files in the PAM50 (see #3 ).

Philips Lab created updated spinal levels files : Philips Lab spinal level files from the Mendez (2021).

The measures used from Mendez: Segment length at dorsal column entry

To bring the spinal level positioning in the PAM50 space, linear scaling between length of Mendez and the length of pixels between C2 and L5 in the PAM50 was done.

However, the spinal segments length doesn't vary linearly along the spinal cord. Using the length between C2-L5 won't necessarly be accurate for all levels to get the position in the PAM50 space.

New spinal segment positions vs before on PAM50

Here are the different positions of the center of each segment with the updated spinal levels vs before:

<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:x="urn:schemas-microsoft-com:office:excel" xmlns="http://www.w3.org/TR/REC-html40">

Level | Mendez (Slices in I-->S) | Old (Slices in I-->S) | Difference -- | -- | -- | -- C2 | 954 | 945 | 9 C3 | 927 | 927 | 0 C4 | 899 | 891 | 8 C5 | 868 | 857 | 11 C6 | 839 | 826 | 13 C7 | 812 | 794 | 18 C8 | 787 | 757 | 30 T1 | 760 | 713 | 47 T2 | 729 | 668 | 61 T3 | 695 | 622 | 73 T4 | 655 | 572 | 83 T5 | 608 | 522 | 86 T6 | 558 | 471 | 87 T7 | 507 | 417 | 90 T8 | 454 | 365 | 89 T9 | 403 | 310 | 93 T10 | 351 | 251 | 100 T11 | 301 | 194 | 107 T12 | 259 | 169 | 90 L1 | 221 |   |   L2 | 188 |   |   L3 | 162 |   |   L4 | 137 |   |   L5 | 114 |   |   ### Examples Here are some examples: **Mendez --> yellow; Old --> blue** #### Segment C3 ![c3](https://user-images.githubusercontent.com/71230552/229629141-db448adf-c96a-4240-9c78-79c21018d73c.gif) ![image](https://user-images.githubusercontent.com/71230552/229597609-2356161c-5853-4188-a48e-59eec66d61f6.png) The center position is the same. We see the rootlet in the axial view. #### Segment C4 ![c4](https://user-images.githubusercontent.com/71230552/229629251-244784b6-f7b3-4cea-8e9e-92fcc07eee38.gif) **Axial view at the center of the segment for the new and old spinal segment** ![image](https://user-images.githubusercontent.com/71230552/229626454-e2e2ea2a-d22d-4b1e-aa87-79cf0616e5e6.png) ![image](https://user-images.githubusercontent.com/71230552/229626563-56adbb9f-ee22-4565-aaaf-6de301fd801c.png) --> Not much difference of the appearance of the spinal rootlet #### Segment C5 ![c5](https://user-images.githubusercontent.com/71230552/229629511-d021bf18-de1f-4b77-8805-35bf57f19d85.gif) **Axial view at the center of the segment for the new and old spinal segment** ![image](https://user-images.githubusercontent.com/71230552/229626662-7bafbef0-d4bc-4f9a-8134-a10354f37cde.png) ![image](https://user-images.githubusercontent.com/71230552/229626750-457a5f49-a48f-47fa-aecd-433e8c277d26.png) --> We can better see the spinal rootlet in the old spinal segment file #### Segment T6 ![T6](https://user-images.githubusercontent.com/71230552/229629594-b7ad15ac-bc39-44ea-8bdf-c702aa64c977.gif) #### Segment T10 ![T10](https://user-images.githubusercontent.com/71230552/229629617-01e346f4-c7b8-454f-bf9b-6a5d20f3b740.gif) ### Suggestion We need a surrogate measure to get the position of each spinal segment that we can relate to the PAM50. Measures that relate the vertebrae to the spinal rootlets could be relevant here since we can position it for each level in the PAM50. - Intervertebral foramen to rostral rootlet ![image](https://user-images.githubusercontent.com/71230552/229630859-816e8366-608a-4225-b893-e36d04e49793.png) - Intervertebral foramen to caudal rootlet ![image](https://user-images.githubusercontent.com/71230552/229631409-83b61cb3-0324-4f98-bbb0-f8e3f949e201.png) This would require identifying the intervertebral foramen on the PAM50 template