Difference between revisions of "Image Registration for MRg Prostate Biopsy"

From NAMIC Wiki
Jump to: navigation, search
m (moved Image registration with BRAINSFitIGT to Image Registration for MRg Prostate Biopsy: use a more descriptive title that correspond to the purpose of this page)
Line 1: Line 1:
== Intro and Disclaimer ==
+
=Summary=
  
BRAINSFitIGT code is a 3D Slicer module based on the [http://wiki.slicer.org/slicerWiki/index.php/Modules:BRAINSFit BRAINSFit] code originally developed by Hans Johnson and his team at U.Iowa. BRAINSFitIGT has been modified and customized to support registration of prostate preoperative MRI to intraoperative MRI during image-guided biopsy procedures. It has been used to produce reasonable registration results for the last 7 or so biopsy cases. Quantitative evaluation of the accuracy is in progress.
+
BRAINSFitIGT code is a 3D Slicer module based on the [http://wiki.slicer.org/slicerWiki/index.php/Modules:BRAINSFit BRAINSFit] code originally developed by Hans Johnson and his team at U.Iowa. BRAINSFitIGT has been modified and customized to support registration of prostate preoperative MRI to intraoperative MRI during image-guided biopsy procedures. This module is NOT included in the binary distribution of 3D Slicer: you need to compile Slicer version 3 from source code to access it. This page explains how to compile and use this module.
  
While using this code, please remember that it can change at any time, as I (Andrey Fedorov) am working to evaluate and improve it. It has not been evaluated for organs other than prostate. I will be interested to learn about your experience using it for new applications.
+
BRAINSFitIGT module has been developed specifically to register pre-procedural T2w MRI that are typically using endorectal coil to intra-procedural MRI that do not use endorectal coil.  
  
'''Please treat all the training data provided to you for testing as non-public. You MUST check with me or with Dr. Noby Hata before you copy the data outside the SPL network or use it in any publication!'''
+
=Build instructions=
 +
{|style="text-align:left;" border="0"
 +
|
 +
You must compile 3D Slicer version 3 from source. General build instructions and prerequisites can be found here: http://www.slicer.org/slicerWiki/index.php/Slicer3:Build_Instructions. You should also enable BRAINSFit module in the CMake options, as shown in the screnshot.
  
== General instructions ==
+
|[[Image:BRAINSFitIGT_build_option.png]]
 +
|-
 +
|}
  
* The software components used in this pipeline are not available in the standard Slicer distribution. Instead, the following binary '''available on opel.bwh.harvard.edu''' should be used:
+
=Usage instructions=
<pre>
 
/workspace/fedorov/Slicer/Stable-Release/Slicer3-superbuild/Slicer3-build/Slicer3
 
</pre>
 
* log on to opel with X forwarding enabled:
 
<pre>
 
ssh -XY opel
 
</pre>
 
* the code is also available on the cart workstation at SPL. Log in using "develop" account (email me fedorov AT bwh) if you need the password. Use the Slicer icon on the desktop to launch Slicer with BRAINSFitIGT.
 
* You are supposed to be familiar with 3D Slicer. If this is not the case, familiarize yourself at least with the tutorials 1, 2, 4 and 5 from [http://wiki.slicer.org/slicerWiki/index.php/Slicer_3.6:Training Slicer 3.6 training page]
 
  
== WARNINGS ==
+
You are supposed to be familiar with 3D Slicer. If this is not the case, familiarize yourself at least with the tutorials 1, 2, 4 and 5 from [http://wiki.slicer.org/slicerWiki/index.php/Slicer_3.6:Training Slicer 3.6 training page].
  
Do not's ...
+
Before registering pre-procedural prostate T2w MRI, the following workflow is advised:
 +
# prepare approximate contouring of the prostate gland in the pre-procedural T2w image
 +
# apply N4 inhomogeneity correction (using [http://www.slicer.org/slicerWiki/index.php/Modules:N4ITKBiasFieldCorrection-Documentation-3.6 this module of Slicer]) to the pre-procedural image, in case it was acquired using endorectal coil. While performing inhomogeneity correction, use the prostate mask prepared in the previous step. The recommended number of iterations is "500,400,300".
  
* do not make slice visible in 3D viewer: this will crash Slicer (this is a known problem specific to SPL fat nodes, there is no solution to fix it right now)
+
Next, proceed with the registration steps as described below.
* do not use your home directory for Slicer temporary storage -- it will fill up quickly, and Slicer will not work correctly
 
 
 
== Processing Steps ==
 
  
 
{|style="text-align:left;" border="0"
 
{|style="text-align:left;" border="0"
Line 37: Line 32:
 
## create a new label image
 
## create a new label image
 
## trace the approximate boundary of the prostate capsule in each image using "Draw" tool (let's say, they are called "preop-label" "intraop-label")
 
## trace the approximate boundary of the prostate capsule in each image using "Draw" tool (let's say, they are called "preop-label" "intraop-label")
# Perform bias correction on the intraop image
 
## go to "N4ITK MRI Bias Correction" module
 
## ''Parameters: Input image = "intraop" or "preop" image, Mask image = "intraop-label" or "preop-label", Output volume = create new volume and name it "preop-N4"/"intraop-N4"''
 
  
 
'''2. Registration step'''
 
'''2. Registration step'''
Line 49: Line 41:
 
|-
 
|-
 
|}
 
|}
 
== Testing Data ==
 
Example datasets of prostate MRI acquired before and during Mr-guided biopsy procedures are available on the cart workstation at SPL:
 
* load a training dataset by clicking on File --> Load Scene... --> Prostate_registration_training
 
* there should be three training datasets: Case003, Case 004. Case005
 
* choose e.g. Case004/RegistrationScene.mrml
 
* nomenclature:
 
** t2ax = preoperative prostate image data
 
** t2ax-intraop = intraoperative prostate image data
 
The same training data is also available on SPL network in /projects/igtcases/ProstateBx_cases/Training/
 

Revision as of 15:59, 28 March 2012

Home < Image Registration for MRg Prostate Biopsy

Summary

BRAINSFitIGT code is a 3D Slicer module based on the BRAINSFit code originally developed by Hans Johnson and his team at U.Iowa. BRAINSFitIGT has been modified and customized to support registration of prostate preoperative MRI to intraoperative MRI during image-guided biopsy procedures. This module is NOT included in the binary distribution of 3D Slicer: you need to compile Slicer version 3 from source code to access it. This page explains how to compile and use this module.

BRAINSFitIGT module has been developed specifically to register pre-procedural T2w MRI that are typically using endorectal coil to intra-procedural MRI that do not use endorectal coil.

Build instructions

You must compile 3D Slicer version 3 from source. General build instructions and prerequisites can be found here: http://www.slicer.org/slicerWiki/index.php/Slicer3:Build_Instructions. You should also enable BRAINSFit module in the CMake options, as shown in the screnshot.

BRAINSFitIGT build option.png

Usage instructions

You are supposed to be familiar with 3D Slicer. If this is not the case, familiarize yourself at least with the tutorials 1, 2, 4 and 5 from Slicer 3.6 training page.

Before registering pre-procedural prostate T2w MRI, the following workflow is advised:

  1. prepare approximate contouring of the prostate gland in the pre-procedural T2w image
  2. apply N4 inhomogeneity correction (using this module of Slicer) to the pre-procedural image, in case it was acquired using endorectal coil. While performing inhomogeneity correction, use the prostate mask prepared in the previous step. The recommended number of iterations is "500,400,300".

Next, proceed with the registration steps as described below.

1. Preparation step: start Slicer, load the images (let's say, they are called "preop.nrrd" and "intraop.nrrd"). Do the following:

  1. Contour the object you want to register in both images
    1. go to Slicer Editor module
    2. create a new label image
    3. trace the approximate boundary of the prostate capsule in each image using "Draw" tool (let's say, they are called "preop-label" "intraop-label")

2. Registration step

  1. Choose Registration -> BRAINSFitIGT module
    1. Parameters: Fixed image volume = "intraop-N4", Moving image volume = "preop-N4", Input fixed mask = "intraop-label", Input moving mask = "preop-label", Slicer BSpline transform = create a new transform and name it "bspline-REG-tfm", Output image volume = create a new volume and name it "preop-REG"
    2. Upon completion of registration, inspect the result and check if it makes sense
Fields that need to be selected for the 'Registration step of the workflow are highlighted