Difference between revisions of "Meningioma growth modeling with TumorSim"

From NAMIC Wiki
Jump to: navigation, search
 
(12 intermediate revisions by the same user not shown)
Line 1: Line 1:
 
This is a follow-up project for [http://wiki.na-mic.org/Wiki/index.php/2009_Summer_Project_Week_Meningioma_growth_simulation Meningioma growth simulation project] within 2009 Boston programming week project.
 
This is a follow-up project for [http://wiki.na-mic.org/Wiki/index.php/2009_Summer_Project_Week_Meningioma_growth_simulation Meningioma growth simulation project] within 2009 Boston programming week project.
  
=Seed 1 (default)=
+
=tumorsim1.1rc2 experiments=
  
<flash>file=Seed_1_simulation.swf|width=1000|height=800|quality=best</flash>
+
'''NOTE: red outlines the boundary of the label corresponding to the initial tumor seed!'''
  
 +
{|
 +
|Seed 1 (default seed from Marcel)
 +
<flash>file=Seed_1_simulation.swf|width=500|height=400|quality=best</flash>
 +
|Seed 3 (ventricle)
 +
<flash>file=Seed_3_simulation.swf|width=500|height=400|quality=best</flash>
 +
|}
  
=Seed 3 (ventricle)=
+
=tumorsim1.1.rc4 experiments=
  
<flash>file=Seed_3_simulation.swf|width=1000|height=800|quality=best</flash>
+
{|
 +
|All the following experiments are based on the "thin dura" input labels, and the tumor seed image shown on the right. The tumor label is shown in green color. Note that the tumor is touching the gray matter label.
 +
 
 +
'''The seed image can be downloaded [[File:seed3.nrrd|here]]'''
 +
|<gallery perrow="3" heights="300px" widths="400px">
 +
Image:tumor_seed3.jpg
 +
Image:tumor_seed3_3dview.jpg
 +
</gallery>
 +
|}
 +
 
 +
Also, all of the experiments use the following configuration file, where '''POISSON''' value is substituted with 0.5, 0.4, and 0.2, with the simulation results shown below.
 +
 
 +
<pre>
 +
<?xml version="1.0"?>
 +
<tumor-simulation-parameters>
 +
 
 +
<dataset-name>SimMeningioma001</dataset-name>
 +
 
 +
<input-directory>../TumorSimInputWithThinDura</input-directory>
 +
<output-directory>./0.05</output-directory>
 +
 
 +
<brain-young-modulus>694</brain-young-modulus>
 +
<brain-poisson-ratio>POISSON</brain-poisson-ratio>
 +
 
 +
<deformation-seed>../SlicerScenes/NewSeed3-1.nrrd</deformation-seed>
 +
<deformation-iterations>10</deformation-iterations>
 +
<deformation-initial-pressure>3.0</deformation-initial-pressure>
 +
<deformation-kappa>1e+06</deformation-kappa>
 +
<deformation-damping>1</deformation-damping>
 +
 
 +
<deformation-solver-iterations>8</deformation-solver-iterations>
 +
<infiltration-solver-iterations>2</infiltration-solver-iterations>
 +
 
 +
<infiltration-iterations>2</infiltration-iterations>
 +
<infiltration-time-step>1e-6</infiltration-time-step>
 +
<infiltration-early-time>0.5e-6</infiltration-early-time>
 +
<infiltration-body-force-iterations>2</infiltration-body-force-iterations>
 +
<infiltration-body-force-coefficient>10.0</infiltration-body-force-coefficient>
 +
 
 +
<contrast-enhancement-type>uniform</contrast-enhancement-type>
 +
 
 +
<number-of-threads>16</number-of-threads>
 +
</pre>
 +
 
 +
The following screenshots capture the same axial slice (#119) of the following images in order: (1) image containing the tumor probabilities, and (2) the simulated Gad enhanced image, both overlayed with the warped tumor/brain tissues labels.
 +
 
 +
{|
 +
|'''Poisson ratio 0.5'''
 +
 
 +
(initial value)
 +
|<gallery perrow="3" heights="300px" widths="400px">
 +
Image:ts_rc4_p0.5_prob.jpg
 +
Image:ts_rc4_p0.5_gad.jpg
 +
</gallery>
 +
 
 +
Note very very minimal deformation of the surrounding brain tissue, although the tumor volume increased, for Poisson=0.5.
 +
|}
 +
 
 +
{|
 +
|'''Poisson ratio 0.4'''
 +
|<gallery perrow="3" heights="300px" widths="400px">
 +
Image:ts_rc4_p0.4_prob.jpg
 +
Image:ts_rc4_p0.4_gad.jpg
 +
</gallery>
 +
|}
 +
 
 +
{|
 +
|'''Poisson ratio 0.2'''
 +
|<gallery perrow="3" heights="300px" widths="400px">
 +
Image:ts_rc4_p0.2_prob.jpg
 +
Image:ts_rc4_p0.2_gad.jpg
 +
</gallery>
 +
|}
 +
 
 +
{|
 +
|<gallery perrow="3" heights="300px" widths="400px">
 +
Image:ts_rc4_p0.5.jpg|'''Poisson 0.5'''
 +
Image:ts_rc4_p0.49.jpg|'''Poisson 0.49'''
 +
Image:ts_rc4_p0.45.jpg|'''Poisson 0.45'''
 +
</gallery>
 +
|}
 +
 
 +
= Standing questions/items =
 +
 
 +
* limited tumor growth/surrounding brain deformation with 10 iterations -- try more iterations?
 +
* reduction of Poisson ratio makes results very unrealistic
 +
* rerun the previous version of tumorsim, which seemed to show nice more realistic deformation of the brain tissue (with Poisson=0.5)

Latest revision as of 13:42, 5 October 2009

Home < Meningioma growth modeling with TumorSim

This is a follow-up project for Meningioma growth simulation project within 2009 Boston programming week project.

tumorsim1.1rc2 experiments

NOTE: red outlines the boundary of the label corresponding to the initial tumor seed!

Seed 1 (default seed from Marcel)

<flash>file=Seed_1_simulation.swf|width=500|height=400|quality=best</flash>

Seed 3 (ventricle)

<flash>file=Seed_3_simulation.swf|width=500|height=400|quality=best</flash>

tumorsim1.1.rc4 experiments

All the following experiments are based on the "thin dura" input labels, and the tumor seed image shown on the right. The tumor label is shown in green color. Note that the tumor is touching the gray matter label.

The seed image can be downloaded File:Seed3.nrrd

Also, all of the experiments use the following configuration file, where POISSON value is substituted with 0.5, 0.4, and 0.2, with the simulation results shown below.

<?xml version="1.0"?>
<tumor-simulation-parameters>

<dataset-name>SimMeningioma001</dataset-name>

<input-directory>../TumorSimInputWithThinDura</input-directory>
<output-directory>./0.05</output-directory>

<brain-young-modulus>694</brain-young-modulus>
<brain-poisson-ratio>POISSON</brain-poisson-ratio>

<deformation-seed>../SlicerScenes/NewSeed3-1.nrrd</deformation-seed>
<deformation-iterations>10</deformation-iterations>
<deformation-initial-pressure>3.0</deformation-initial-pressure>
<deformation-kappa>1e+06</deformation-kappa>
<deformation-damping>1</deformation-damping>

<deformation-solver-iterations>8</deformation-solver-iterations>
<infiltration-solver-iterations>2</infiltration-solver-iterations>

<infiltration-iterations>2</infiltration-iterations>
<infiltration-time-step>1e-6</infiltration-time-step>
<infiltration-early-time>0.5e-6</infiltration-early-time>
<infiltration-body-force-iterations>2</infiltration-body-force-iterations>
<infiltration-body-force-coefficient>10.0</infiltration-body-force-coefficient>

<contrast-enhancement-type>uniform</contrast-enhancement-type>

<number-of-threads>16</number-of-threads>

The following screenshots capture the same axial slice (#119) of the following images in order: (1) image containing the tumor probabilities, and (2) the simulated Gad enhanced image, both overlayed with the warped tumor/brain tissues labels.

Poisson ratio 0.5

(initial value)

Note very very minimal deformation of the surrounding brain tissue, although the tumor volume increased, for Poisson=0.5.

Poisson ratio 0.4
Poisson ratio 0.2

Standing questions/items

  • limited tumor growth/surrounding brain deformation with 10 iterations -- try more iterations?
  • reduction of Poisson ratio makes results very unrealistic
  • rerun the previous version of tumorsim, which seemed to show nice more realistic deformation of the brain tissue (with Poisson=0.5)