Difference between revisions of "Documentation/4.8/Modules/GeodesicSlicer"

From Slicer Wiki
Jump to: navigation, search
 
(37 intermediate revisions by the same user not shown)
Line 3: Line 3:
 
{{documentation/{{documentation/version}}/module-header}}
 
{{documentation/{{documentation/version}}/module-header}}
 
<!-- ---------------------------- -->
 
<!-- ---------------------------- -->
[[File:GeodesicSlicer logo.png|128x128px|thumb|center|GeodesicSlicer logo]]
+
[[File:GeodesicSlicer logo.png|128x128px|thumb|left|GeodesicSlicer logo]]
 +
[[File:Screen-shot of the GeodesicSlicer program. .png|thumb|right|The users enters 1) the T1-weighted whole-brain anatomical image 2) Place four points: the nasion, the inion, the left tragus and the right tragus. The program make a 3D mesh morphed to the structural MRI data of a participant and calculates the 10-20 system EEG with T3P3, and outputs the distance between the anatomical target and the T3 electrode.]]
 
<!-- ---------------------------- -->
 
<!-- ---------------------------- -->
 
{{documentation/{{documentation/version}}/module-section|Introduction and Acknowledgements}}
 
{{documentation/{{documentation/version}}/module-section|Introduction and Acknowledgements}}
 
{{documentation/{{documentation/version}}/module-introduction-start|{{documentation/modulename}}}}
 
{{documentation/{{documentation/version}}/module-introduction-start|{{documentation/modulename}}}}
 
{{documentation/{{documentation/version}}/module-introduction-row}}
 
{{documentation/{{documentation/version}}/module-introduction-row}}
: '''Author(s)/Contributor(s):''' Frederic Briend (ISTS, UNICAEN), Antoine Nourry (UMS 3408)<br>
+
: '''Author(s)/Contributor(s):''' Frederic Briend (ISTS EA 7466, UNICAEN), Antoine Nourry (UMS 3408)<br>
: '''Acknowledgements:''' This work was supported by CHU Caen, Region Normandie and UNICAEN.<br>
+
: '''Acknowledgements:''' This work was supported by a Perceneige-Fondamental prize, CHU Caen, Region Normandie and UNICAEN.<br>
 
: '''Contact:''' Frederic Briend, <email>briend@cyceron.fr</email><br>
 
: '''Contact:''' Frederic Briend, <email>briend@cyceron.fr</email><br>
[[File:Logo ISTS Elise final trans petit.jpg|146x122px|thumb|center|ISTS]]
 
 
{{documentation/{{documentation/version}}/module-introduction-row}}
 
{{documentation/{{documentation/version}}/module-introduction-row}}
 
{{documentation/{{documentation/version}}/module-introduction-end}}
 
{{documentation/{{documentation/version}}/module-introduction-end}}
  
The module has been developed based on ideas and feedback from the community. We would like to especially thank:
+
 
 +
The module has been developed based on ideas and feedbacks from the community. We would like to especially thank:
 
* Dr. Olivier Etard, M.D., Ph.D., CHU de Caen.
 
* Dr. Olivier Etard, M.D., Ph.D., CHU de Caen.
 +
* Dr. Clément Nathou, M.D., Ph.D., CHU de Caen.
 +
* Dr. Nicolas Delcroix, Ph.D., UMS 3408.
 
* Dr. Sonia Dollfus, M.D., Ph.D., CHU de Caen, header of [http://www.ists.cyceron.fr/spip.php?rubrique17 ISTS].
 
* Dr. Sonia Dollfus, M.D., Ph.D., CHU de Caen, header of [http://www.ists.cyceron.fr/spip.php?rubrique17 ISTS].
 
* Dr. Csaba Pinter, MSc, Queen's University.
 
* Dr. Csaba Pinter, MSc, Queen's University.
 
* Dr. Andras Lasso, Ph.D., Queen's University.
 
* Dr. Andras Lasso, Ph.D., Queen's University.
 +
 +
''If you use this module, please cite the following article: <ref name="Briend 2018">Briend F. et al., Repetitive transcranial magnetic stimulation (rTMS) treatment for auditory hallucinations: personalized or standardized targets? Brain Stimulation, submitted</ref>.''
  
 
__TOC__
 
__TOC__
Line 25: Line 30:
 
<!-- ---------------------------- -->
 
<!-- ---------------------------- -->
 
{{documentation/{{documentation/version}}/module-section|Module Description}}
 
{{documentation/{{documentation/version}}/module-section|Module Description}}
This module calculates geodesic path in 3D structure. Thanks to this geodesic path, this module could draw an EEG 10-20 system.
+
This module calculates geodesic path in 3D structure. Thanks to this geodesic path, this module can draw an EEG 10-20 system, determine the projected scalp stimulation site (MRI guided brain stimulation without the use of a neuronavigation System) and correct the rTMS resting motor threshold by correction factor.
  
 
'''Terminology'''
 
'''Terminology'''
 
*'''''Mesh''''' A mesh or polygon mesh is a collection of vertices, edges and faces that defines the shape of a polyhedral object in 3D computer graphics and solid modeling.
 
*'''''Mesh''''' A mesh or polygon mesh is a collection of vertices, edges and faces that defines the shape of a polyhedral object in 3D computer graphics and solid modeling.
 
*'''''Shortest path''''' In graph theory, the shortest path problem is the problem of finding a path between two vertices (or nodes) in a graph such that the sum of the weights of its constituent edges is minimized.
 
*'''''Shortest path''''' In graph theory, the shortest path problem is the problem of finding a path between two vertices (or nodes) in a graph such that the sum of the weights of its constituent edges is minimized.
*'''''10-20 EEG system''''' The International 10-20 system is commonly used for EEG electrode placement and for correlating external skull locations to underlying cortical areas.<ref name="Jasper 1958">Jasper, H.H. The ten-twenty electrode system of the International Federation. Electroencephalogr. Clin. Neurophysiol., 1958, 10: 370-375</ref>  
+
*'''''10-20 EEG system''''' The International 10-20 system is commonly used for electroencephalogram (EEG) electrode placement and to correlate external skull locations with underlying cortical areas.<ref name="Jasper 1958">Jasper, H. (1958). The ten twenty electrode system of the international federation. Electroencephalography and Clinical Neurophysiology, 10, 371‑375.</ref>  
 +
 
 +
<!-- ---------------------------- -->
 +
{{documentation/{{documentation/version}}/module-section|Installation (in progress)}}
 +
# First, open 3D Slicer
 +
# Open the Slicer Extensions from the icon on the menu bar
 +
# Choose "Geodesic Slicer" module from the list of extensions and click "INSTALL" button.
 +
# Once you restart 3D Slicer, the Geodesic Slicer module should show up on the Modules menu (under Informatics->Geodesic Slicer)
  
 
<!-- ---------------------------- -->
 
<!-- ---------------------------- -->
 
{{documentation/{{documentation/version}}/module-section|Use Cases}}
 
{{documentation/{{documentation/version}}/module-section|Use Cases}}
 
The overall goal is to allow users to find the shortest paths between nodes in a graph and via the Dijkstra's algorithm to make  10-20 system.  This module can be used for:
 
The overall goal is to allow users to find the shortest paths between nodes in a graph and via the Dijkstra's algorithm to make  10-20 system.  This module can be used for:
*Stimulation in psychiatry
+
*Stimulation in psychiatry: MRI guided brain stimulation without the use of a neuronavigation system.
*Surgery measurement
+
*Surgery measurement.
*3D printing
+
*3D printing.
  
 
<!-- ---------------------------- -->
 
<!-- ---------------------------- -->
 
{{documentation/{{documentation/version}}/module-section|Panels and their use}}
 
{{documentation/{{documentation/version}}/module-section|Panels and their use}}
 
[[File:Wiki main picture.png|thumb|center|Geodesic Slicer. Panel set]]
 
  
 
==== Create a mesh ====
 
==== Create a mesh ====
 
+
[[File:Create mesh.png|thumb|right]]
Image à insérer: "wiki_main_picture2.png"
 
  
 
A typical straightforward Geodesic Slicer workflow for consists of the following steps:
 
A typical straightforward Geodesic Slicer workflow for consists of the following steps:
  
# Load a volume.nii.
+
# Load a volume.nii (by Drag & Drop or the Add Data dialogue).
# Enter the Geodesic Slicer module using either the toolbar or the Modules menu button.
+
# Enter in the Geodesic Slicer module using either the toolbar or the Modules menu button.
# Create a mesh.
+
# Press the button "Create a quick mesh" or "Create a mesh" (with filling holes smoothing, better for the next part but longer).
#*Wait a moment
+
#*Wait a moment.
#Use this mesh.
 
#*If your background was named 't1' the label will be called 't1-label'.
 
 
#Go to '''Parameters to find the shortest path''' or '''Make 10-20 EEG system electrode''' section.
 
#Go to '''Parameters to find the shortest path''' or '''Make 10-20 EEG system electrode''' section.
  
 
==== Parameters to find the shortest path ====
 
==== Parameters to find the shortest path ====
# Source points.
+
[[File:Shortest past.png|thumb|right]]
# Input STL model.
+
 
#* Calculate in centimeter the shortest path via the Dijkstra's algorithm.
+
# Source points: The list of fiducial points on the curve, since the "Create-and-place Fiducial" button (in green in the figure above).  
#* Draw the Dijkstra's algorithm shortest path.
+
# Input STL model: The model you use (after "use this mesh", the T1.stl created).
#** Length (cm)
+
#* Find the shortest path: Calculate in centimeter the geodesic (shortest) path via the Dijkstra's algorithm.
 +
#* Draw the shortest path: Draw the Dijkstra's algorithm shortest path.
 +
#** Length (cm): The length of the current curve is shown in centimeter.
 +
 
 +
====10-20 system electrode====
 +
*Run the Dijkstra's algorithm to '''make the 10-20 system electrode'''.
 +
 
 +
[[File:4 landmarks.png|thumb|right|Four anatomical landmarks are used for the essential positioning of the electrodes: the nasion, the inion, the pre auricular to the left ear and the pre auricular to the right ear. ]]
 +
# 4 anatomical landmarks: (Sources Points) The list of fiducial points on the curve, since the "Create-and-place Fiducial" button (in green in the figure above). Four anatomical landmarks are used for the essential positioning of the electrodes (in this order!):
 +
#* 1/ The nasion
 +
#* 2/ The inion
 +
#* 3/ The pre auricular to the left ear
 +
#* 4/ The pre auricular to the right ear
 +
# Input STL model: The model you use (after "use this mesh", the T1.stl created).
 +
# Press the button "Make 10-20 EEG system electrode" to draw the 10-20 EEG system via the Dijkstra's algorithm.
 +
#* The traditional T3P3 site according to the International 10–20 system of electroencephalogram was identified.
 +
 
 +
*'''Project the stimulation site''' on the 10-20 system electrode distances and characterize it.
 +
# Stimulation Site placed: Place on the T1-weighted anatomical image the stimulation point that you want since the "Create-and-place Fiducial" button. Once this point given, click on 'Yes'.
 +
# Press the button "Project the stimulation site" to project the stimulation point on the scalp and find the 3 nearest electrodes around it.
 +
#* Nearest electrode 1: The distance in centimeter between the first nearest electrode and the projected stimulation site.
 +
#* Nearest electrode 2: The distance in centimeter between the second nearest electrode and the projected stimulation site.
 +
#* Nearest electrode 3: The distance in centimeter between the third nearest electrode and the projected stimulation site.
  
==10-20 system electrode==
 
Run the Dijkstra's algorithm to make the 10-20 system electrode distances
 
  
===4 anatomical landmarks===
+
====rTMS resting motor threshold- Correction factor====
Four anatomical landmarks are used for the essential positioning of the electrodes:
+
Calculate correction factors to adjust the rTMS dose for the treatment (according to the depth of the stimulation site).
* 1/the nasion
 
* 2/the inion
 
* 3/the pre auricular to the left ear
 
* 4/the pre auricular to the right ear
 
  
===Input STL model===
+
[[File:M1.png|thumb|Localization of the motor hand area via a knob on the precentral gyrus]]
This model required to contain a dense set of points.
 
  
===Make 10-20 EEG system electrode===
+
# M1 Point Placed: Place on the T1-weighted anatomical image a point targeting the human [https://en.wikipedia.org/wiki/Motor_cortex|motor motor cortex] since the "Create-and-place Fiducial" button. Once this point given, click on 'Yes'. Help via the [https://pdfs.semanticscholar.org/ba38/045e9f01ec4d128c5fbe5a46dc209fccaac4.pdf Yousry's method].
Make 10-20 EEG system electrode via the Dijkstra's algorithm.
+
# Set the stimulation intensity of the  resting motor threshold.
* Distance T3-Stim in the T3P3 axis are give in centimeters.
+
# Press the button "Correct the motor threshold" to correct the unadjusted motor threshold (rMT) in % stimulator output.
 +
#* Two adjusted motor threshold (AdjMT%) in % stimulator output are given where SCDx is the scalp-to-cortex distance between the scalp and and the Stimulation Site, SCDm is the scalp-to-cortex distance between the scalp and M1.
 +
#* 1/ The first according to Stokes et al. Clin  Neurophysiol 2007 <ref name="Stokes 2007">Stokes, M. G., Chambers, C. D., Gould, I. C., English, T., McNaught, E., McDonald, O., & Mattingley, J. B. (2007). Distance-adjusted motor threshold for transcranial magnetic stimulation. Clinical Neurophysiology, 118(7), 1617‑1625.</ref> , where [AdjMT% = 2,7*(SCDx - SCDm) + rMT]
 +
#* 2/ The second according to Hoffman et al. Biol Psychiatry 2013 <ref name="Hoffman 2013">Hoffman, R. E., Wu, K., Pittman, B., Cahill, J. D., Hawkins, K. A., Fernandez, T., & Hannestad, J. (2013). Transcranial magnetic stimulation of Wernicke’s and Right homologous sites to curtail « voices »: a randomized trial. Biological Psychiatry, 73(10), 1008‑1014. </ref> , where [AdjMT% = 0.90*rMT*e0.036*(SCDx-SCDm)]
  
 
<!-- ---------------------------- -->
 
<!-- ---------------------------- -->

Latest revision as of 18:12, 7 January 2019

Home < Documentation < 4.8 < Modules < GeodesicSlicer


For the latest Slicer documentation, visit the read-the-docs.


GeodesicSlicer logo
The users enters 1) the T1-weighted whole-brain anatomical image 2) Place four points: the nasion, the inion, the left tragus and the right tragus. The program make a 3D mesh morphed to the structural MRI data of a participant and calculates the 10-20 system EEG with T3P3, and outputs the distance between the anatomical target and the T3 electrode.

Introduction and Acknowledgements

Author(s)/Contributor(s): Frederic Briend (ISTS EA 7466, UNICAEN), Antoine Nourry (UMS 3408)
Acknowledgements: This work was supported by a Perceneige-Fondamental prize, CHU Caen, Region Normandie and UNICAEN.
Contact: Frederic Briend, <email>briend@cyceron.fr</email>


The module has been developed based on ideas and feedbacks from the community. We would like to especially thank:

  • Dr. Olivier Etard, M.D., Ph.D., CHU de Caen.
  • Dr. Clément Nathou, M.D., Ph.D., CHU de Caen.
  • Dr. Nicolas Delcroix, Ph.D., UMS 3408.
  • Dr. Sonia Dollfus, M.D., Ph.D., CHU de Caen, header of ISTS.
  • Dr. Csaba Pinter, MSc, Queen's University.
  • Dr. Andras Lasso, Ph.D., Queen's University.

If you use this module, please cite the following article: [1].

Module Description

This module calculates geodesic path in 3D structure. Thanks to this geodesic path, this module can draw an EEG 10-20 system, determine the projected scalp stimulation site (MRI guided brain stimulation without the use of a neuronavigation System) and correct the rTMS resting motor threshold by correction factor.

Terminology

  • Mesh A mesh or polygon mesh is a collection of vertices, edges and faces that defines the shape of a polyhedral object in 3D computer graphics and solid modeling.
  • Shortest path In graph theory, the shortest path problem is the problem of finding a path between two vertices (or nodes) in a graph such that the sum of the weights of its constituent edges is minimized.
  • 10-20 EEG system The International 10-20 system is commonly used for electroencephalogram (EEG) electrode placement and to correlate external skull locations with underlying cortical areas.[2]

Installation (in progress)

  1. First, open 3D Slicer
  2. Open the Slicer Extensions from the icon on the menu bar
  3. Choose "Geodesic Slicer" module from the list of extensions and click "INSTALL" button.
  4. Once you restart 3D Slicer, the Geodesic Slicer module should show up on the Modules menu (under Informatics->Geodesic Slicer)

Use Cases

The overall goal is to allow users to find the shortest paths between nodes in a graph and via the Dijkstra's algorithm to make 10-20 system. This module can be used for:

  • Stimulation in psychiatry: MRI guided brain stimulation without the use of a neuronavigation system.
  • Surgery measurement.
  • 3D printing.

Panels and their use

Create a mesh

Create mesh.png

A typical straightforward Geodesic Slicer workflow for consists of the following steps:

  1. Load a volume.nii (by Drag & Drop or the Add Data dialogue).
  2. Enter in the Geodesic Slicer module using either the toolbar or the Modules menu button.
  3. Press the button "Create a quick mesh" or "Create a mesh" (with filling holes smoothing, better for the next part but longer).
    • Wait a moment.
  4. Go to Parameters to find the shortest path or Make 10-20 EEG system electrode section.

Parameters to find the shortest path

Shortest past.png
  1. Source points: The list of fiducial points on the curve, since the "Create-and-place Fiducial" button (in green in the figure above).
  2. Input STL model: The model you use (after "use this mesh", the T1.stl created).
    • Find the shortest path: Calculate in centimeter the geodesic (shortest) path via the Dijkstra's algorithm.
    • Draw the shortest path: Draw the Dijkstra's algorithm shortest path.
      • Length (cm): The length of the current curve is shown in centimeter.

10-20 system electrode

  • Run the Dijkstra's algorithm to make the 10-20 system electrode.
Four anatomical landmarks are used for the essential positioning of the electrodes: the nasion, the inion, the pre auricular to the left ear and the pre auricular to the right ear.
  1. 4 anatomical landmarks: (Sources Points) The list of fiducial points on the curve, since the "Create-and-place Fiducial" button (in green in the figure above). Four anatomical landmarks are used for the essential positioning of the electrodes (in this order!):
    • 1/ The nasion
    • 2/ The inion
    • 3/ The pre auricular to the left ear
    • 4/ The pre auricular to the right ear
  2. Input STL model: The model you use (after "use this mesh", the T1.stl created).
  3. Press the button "Make 10-20 EEG system electrode" to draw the 10-20 EEG system via the Dijkstra's algorithm.
    • The traditional T3P3 site according to the International 10–20 system of electroencephalogram was identified.
  • Project the stimulation site on the 10-20 system electrode distances and characterize it.
  1. Stimulation Site placed: Place on the T1-weighted anatomical image the stimulation point that you want since the "Create-and-place Fiducial" button. Once this point given, click on 'Yes'.
  2. Press the button "Project the stimulation site" to project the stimulation point on the scalp and find the 3 nearest electrodes around it.
    • Nearest electrode 1: The distance in centimeter between the first nearest electrode and the projected stimulation site.
    • Nearest electrode 2: The distance in centimeter between the second nearest electrode and the projected stimulation site.
    • Nearest electrode 3: The distance in centimeter between the third nearest electrode and the projected stimulation site.


rTMS resting motor threshold- Correction factor

Calculate correction factors to adjust the rTMS dose for the treatment (according to the depth of the stimulation site).

Localization of the motor hand area via a knob on the precentral gyrus
  1. M1 Point Placed: Place on the T1-weighted anatomical image a point targeting the human motor cortex since the "Create-and-place Fiducial" button. Once this point given, click on 'Yes'. Help via the Yousry's method.
  2. Set the stimulation intensity of the resting motor threshold.
  3. Press the button "Correct the motor threshold" to correct the unadjusted motor threshold (rMT) in % stimulator output.
    • Two adjusted motor threshold (AdjMT%) in % stimulator output are given where SCDx is the scalp-to-cortex distance between the scalp and and the Stimulation Site, SCDm is the scalp-to-cortex distance between the scalp and M1.
    • 1/ The first according to Stokes et al. Clin Neurophysiol 2007 [3] , where [AdjMT% = 2,7*(SCDx - SCDm) + rMT]
    • 2/ The second according to Hoffman et al. Biol Psychiatry 2013 [4] , where [AdjMT% = 0.90*rMT*e0.036*(SCDx-SCDm)]

Information for Developers

The code is available at Github.

References

  1. Briend F. et al., Repetitive transcranial magnetic stimulation (rTMS) treatment for auditory hallucinations: personalized or standardized targets? Brain Stimulation, submitted
  2. Jasper, H. (1958). The ten twenty electrode system of the international federation. Electroencephalography and Clinical Neurophysiology, 10, 371‑375.
  3. Stokes, M. G., Chambers, C. D., Gould, I. C., English, T., McNaught, E., McDonald, O., & Mattingley, J. B. (2007). Distance-adjusted motor threshold for transcranial magnetic stimulation. Clinical Neurophysiology, 118(7), 1617‑1625.
  4. Hoffman, R. E., Wu, K., Pittman, B., Cahill, J. D., Hawkins, K. A., Fernandez, T., & Hannestad, J. (2013). Transcranial magnetic stimulation of Wernicke’s and Right homologous sites to curtail « voices »: a randomized trial. Biological Psychiatry, 73(10), 1008‑1014.