Elastic Registration of Prostate MR Images Based on State Estimation of Dynamical Systems
|
|
- Oswin Reed
- 5 years ago
- Views:
Transcription
1 Elastic Registration of Prostate MR Images Based on State Estimation of Dynamical Systems Bahram Marami a,b,c, Suha Ghoul a, Shahin Sirouspour c, David W. Capson d, Sean R. H. Davidson e, John Trachtenberg f and Aaron Fenster a,b,g a Imaging Research Laboratories, Robarts Research Institute, London, Ontario, Canada b Biomedical Engineering Graduate Program, The University of Western Ontario, Canada c Department of Electrical and Computer Engineering, McMaster University, Canada d Department of Electrical and Computer Engineering, University of Victoria, Canada e Ontario Cancer Institute, University Health Network, Toronto, Ontario, Canada f Department of Surgical Oncology, University Health Network, Toronto, Ontario, Canada g Department of Medical Biophysics, University of Western Ontario, London, Ontario, Canada ABSTRACT Magnetic resonance imaging (MRI) is being increasingly used for image-guided biopsy and focal therapy of prostate cancer. A combined rigid and deformable registration technique is proposed to register pre-treatment diagnostic 3 T magnetic resonance (MR) images, with the identified target tumor(s), to the intra-treatment 1.5 T MR images. The pre-treatment 3 T images are acquired with patients in strictly supine position using an endorectal coil, while 1.5 T images are obtained intra-operatively just before insertion of the ablation needle with patients in the lithotomy position. An intensity-based registration routine rigidly aligns two images in which the transformation parameters is initialized using three pairs of manually selected approximate corresponding points. The rigid registration is followed by a deformable registration algorithm employing a generic dynamic linear elastic deformation model discretized by the finite element method (FEM). The model is used in a classical state estimation framework to estimate the deformation of the prostate based on a similarity metric between pre- and intra-treatment images. Registration results using 1 sets of prostate MR images showed that the proposed method can significantly improve registration accuracy in terms of target registration error (TRE) for all prostate substructures. The root mean square (RMS) TRE of 46 manually identified fiducial points was found to be 2.4±1.2 mm, 2.51±1.2 mm, and 2.28±1.22 mm for the whole gland (WG), central gland (CG), and peripheral zone (PZ), respectively after deformable registration. These values are improved from 3.15±1.6 mm, 3.9±1.5 mm, and 3.2±1.73 mm in the WG, CG and PZ, respectively resulted from rigid registration. Registration results are also evaluated based on the Dice similarity coefficient (DSC), mean absolute surface distances (MAD) and maximum absolute surface distances (MAXD) of the WG and CG in the prostate images. Keywords: Prostate MR image registration, elastic model-based registration, motion compensation, focal ablation therapy, finite element method 1. INTRODUCTION Prostate cancer is the most frequently diagnosed non-skin cancer and the third leading cause of cancer death among males in the developed countries. 1 Early detection and treatment can significantly reduce morality rate from prostate cancer. Accurate delivery of treatment is extremely important in localized prostate cancer therapy. Due to the superior performance in differentiating soft tissue structures, MRI systems are employed in imageguided interventions such as prostate biopsy, and radiation and focal therapy of prostate cancer. 2, 3 High-field MRI systems with an endorectal coil can provide detailed images with ability to visualize substructures of the prostate, i.e., central gland (CG) and peripheral zone (PZ), and their margins. These images also allow accurate localizing of the tumor for use in guiding the biopsy needles 2 and focal ablation therapy techniques. 3, 4 Medical Imaging 214: Image Processing, edited by Sebastien Ourselin, Martin A. Styner, Proc. of SPIE Vol. 934, 934F 214 SPIE CCC code: /14/$18 doi: / Proc. of SPIE Vol F-1
2 START ˆx Time Update/Prediction J k no J k 1 ˆx k 1 J Measurement Update k ˆx k yes FINISH ˆx k Preoperative Images Deformation, Transformation and Interpolation T [u ] k _ Intra-operative Images + R Similarity Measure dx c Figure 1: Block diagram of the iterative registration algorithm based on the concept of dynamic state estimation. Focal ablation procedures performed within the bore of the MR scanner are time sensitive, requiring acquisition of the images in a short time. We acquired intra-treatment images using MRI systems with a field-strength of 1.5 T for ablation procedures. 3 Thus, the intra-treatment images acquired by our interventional MRI system had a lower signal-to-noise ratio (SNR), spatial resolution and contrast than those of diagnostic images acquired at 3 T. Accurate registration of the diagnostic pre-treatment images with delineated tumor(s) to intra-treatment images is instrumental in guiding the focal ablation applicators to the identified target and ensuring that the ablation zone covers the target tumor. Prostate tissue usually undergoes non-rigid deformation between imaging sessions due to the presence and absence of the endorectal coil, different inflation volumes of the endorectal coil balloon, patient s position and different bladder volume. Therefore, deformable registration methods are required to align pre-treatment and intra-treatment images. A 3D-3D deformable registration method based on the concept of state estimation for dynamical systems is proposed in this paper. The registration method employs a general linear elastic finite element (FE)-based deformation model. The registration is achieved through a Kalman-like filtering process, which incorporates information from the deformation model and an observation error computed from an intensity-based similarity metric, namely, modality independent neighborhood descriptor (MIND). 5 No specific geometry of the prostate and its substructures is used for constructing the deformation model, rather a cubic volume of tetrahedral finite elements is employed. In the state estimation framework, the estimation process takes into account modeling and observation uncertainties in the form of an unknown process and measurement disturbances. Furthermore, the registration method requires minimal manual intervention of the user and no pre-processing (feature extraction or segmentation) of the images is involved. 2. METHODS The proposed registration method in this paper involves both rigid and deformable registration. In the first step, the pre- and intra-treatment images are rigidly co-registered by minimizing the sum-of-absolute-differences (SAD) of the MINDs between two images. The MIND algorithm extracts distinctive structures in a local neighborhood based on the similarity of small image patches within the image. The extraction of structures using this method is independent from image modalities, contrast, noise and intensity levels of the images. 5 MIND can automatically deal with intensity inhomogeneities around the endorectal coil and no preprocessing of the images is required. In the optimization, the initial values for the rigid transformation parameters are computed based on aligning three manually identified approximate corresponding point pairs in two images. The registration algorithm is robust with respect to the selection of these approximate corresponding points and they are not used for the validation of the registrations. In the second step, the deformation of the prostate is approximated by a dynamic linear elastic model based on a matching criterion between pre- and intra-treatment images using a filtering approach. The dynamics for the tissue deformation are modeled by following discrete-time linear state-space equations x k = Ax k 1 + w k 1 (1a) z k = Hx k + v k (1b) Proc. of SPIE Vol F-2
3 In these equations, x k is the vector of deformation states at time step k, which is comprised of displacements and velocities of the nodal points of the FE mesh, projected to a lower dimensionality. The measurement vector z k is the vector of displacements at control points from their initial positions. Control points, i.e., x c, are evenly distributed inside the FE mesh and their spacing is larger than intra-treatment image voxels. A and H are the state transition and output matrices, respectively, which can be computed based on the dynamic linear elastic FE model 6 and the position of the control points. Moreover, w k and v k are process and measurement noise, and represent uncertainties/errors (model mismatches, unknown external forces, etc) in the deformation dynamics and observation model, respectively. Given any set of pre- and intra-treatment images, the deformable registration is achieved through an iterative routine. The flowchart of the proposed state estimation-based registration algorithm is given in Figure 1 and it involves following steps. 1. Time Update/Prediction: Starting from zero initial states, the vector of state estimates from the previous step ˆx k 1 is updated to ˆx k using ˆx k ˆx k = Aˆx k 1 (2) is an a priori estimate of the states,7 based on which the predicted displacements for the nodal points of the FE mesh, i.e., u k, are computed. 2. Deformation, Transformation and Interpolation: The deformation for the grid of control points x c and intratreatment image R are computed based on the deformation of the FE mesh using elements shape functions. 6 Deformed grids are transformed to the coordinates of the pre-treatment image T using the rigid transformation. The predicted deformed image T [u k ] can then be interpolated from the pre-treatment image volume. 3. Observation Prediction Error: The observation prediction error dx c is computed at control points based on the gradient of the similarity measure as dx c = z k Hˆx k = γ J(R, T [u k ], x c) (3) where J(R, T [u k ], x c) is the SSD between MINDs in two images at x c Measurement Update: The state estimates are updated based on the observation prediction error using ˆx k = ˆx k + Γdx c (4) where Γ is the steady-state Kalman gain, which is computed offline given process and measurement noise covariance matrices. 7 Γ is computed once and can be used in the registration of different image sets. The algorithm terminates when the relative change in the similarity metric between R and T [u] (i.e., J(R, T [u])) is less than a given small number ε or the number of iterations reaches a maximum value. It should be noted that both the deformation and observation models are independent from the image sets. Therefore, γ is the only parameter, which is needed to be tuned for the registration of different image sets. This parameter weights imaged-derived forces to provide a balance between internal elastic energy of the deformation model and external energy defined based on the image similarity measure. A hierarchical multi-level approach (three image resolution levels) is employed in the estimation of non-rigid deformation of the tissue. 3. EXPERIMENTS AND RESULTS The proposed registration algorithm was evaluated in the registration of pre-treatment 3 T T2-weighted MR images of the prostate to their intra-treatment 1.5 T T2-weighted MR images. Ten MR-guided prostate focal laser ablation therapy patients were identified and all patients underwent pre-treatment 3 T MR imaging using a GE Discovery MR75 scanner. Intra-treatment images were obtained using a GE Signa HDxt 1.5 T scanner just before inserting the ablation needle. An endorectal coil was used for acquiring both pre- and intra-treatment images. The pixel sizes in pre- and intra-treatment images were.4.4 mm and mm, respectively and the slice thickness was 3 mm in both image sets. Proc. of SPIE Vol F-3
4 z x 4 2 y 4 (a) (b) Figure 2: (a) Undeformed FE mesh employed in the registration (b) deformed FE mesh after registration of an example image pair. z x y A cubic mesh with 1,476 tetrahedral elements and 2,93 nodal points was employed to construct the deformation model, which is shown in Figure 2a. Using this mesh, an isotropic linear elastic model of deformation with the Young s elasticity modulus E = 15 kpa, Poisson s ratio ν =.49 and a mass density of ρ = 1 g/cm 3 was created. The constructed model is used to estimate the deformation of the prostate in all image pairs. For any given image set, the intra-treatment image volume is scaled and transformed to be placed inside the FE mesh. A regular grid of was also used as control points inside the cube of finite elements. For an example image pair, the deformed 3D FE mesh after registration are given in Figure 2b. Axial views of the registered prostate images with tumor outlines in two different example cases are given in Figure 3. In each row of this figure, pre-treatment images are rigidly (column (a)) and non-rigidly (column (c)) registered to the intra-treatment images in column (b). Tumor outlines were identified by a radiologist in the pre-treatment images. These outlines were transformed and overlaid on the intra-treatment (column (b)) and registered pre-treatment images (column (c)) based on the obtained rigid and non-rigid displacement field. As can be observed from these figures, most of the deformation occurs in the PZ caused by the endorectal coil. For most image cases the tumor cannot be differentiated from the surrounding healthy tissue in the intra-treatment images. These cases demonstrate the need for an algorithm that can accurately register pre-treatment target volumes to intra-treatment images in MRI-guided targeted biopsy or focal ablation therapy. Both pre-treatment and intra-treatment images were manually segmented to determine the volume of the WG and CG in the prostate images and used to evaluate the registration results. The volume of the PZ was computed by subtracting the volume of CG from the WG volume in each image set. Table 1 summarizes the results of both rigid and deformable registrations for the different prostate substructures in terms of average and standard deviations of the TRE, DSC (%), MAD and MAXD. Furthermore, the frequency distribution of TREs for 46 fiducial points (23 points in each region) in the initial rigid alignment and the deformable registration are shown in Figure 4. It is evident from this figure that the distribution of TREs in deformable registration has a higher peak at small values. In contrast, the TRE from the rigid registration is distributed with a wider range. The registration results given in Table 1 and Figure 4 show that deformable registration improves the quality of matching in terms of DSC and TRE especially in the PZ. The average TRE resulted for deformable registration in the PZ is decreased more than that in the CG. This is partly due to the fact that most of the deformation occur in the PZ of the prostate. Also, DSC improvement in the PZ and WG is greater than that in the CG. The average MAD and MAXD results given in Table 1 also show that deformable registration is more effective Proc. of SPIE Vol F-4
5 (a) (b) (c) Figure 3: Axial views of the prostate images with tumor outlines, in each row: (a) rigidly aligned pre-treatment image, (b) intra-treatment image and (c) non-rigidly registered pre-treatment image. Table 1: Summary of registration results (average ± std) based on TRE, DSC, MAD and MAXD TRE (mm) DSC (%) MAD (mm) MAXD (mm) PZ CG WG PZ CG WG CG WG CG WG Rigid 3.2± ± ± ± ± ± ± ± ± ±1.6 Deformable 2.28± ± ± ± ± ± ± ± ± ±.86 in the WG rather than the CG. 4. COMPUTATIONS Computational tasks involved in the proposed registration method can be categorized in two groups: common tasks before starting the registration, which are performed once for all image sets, and tasks for registering each set of images. We implemented all steps of the registration algorithm using MATLAB on a 3.5 GHz Intel(R) Core(TM) i7-397x processor with a 64. GB RAM. Common tasks required 34 min. For the registration of each image set, it required about 2 min to identify three approximate corresponding points in the prostate image pairs, determine the approximate prostate region in the intra-treatment images and rigidly register two images. The time for iterative elastic registration algorithm was about 2 min. After the iterations end, it required 26 sec to compute the deformed pre-treatment image (high-resolution) based on the estimated states of the deformation. Furthermore, the time for finding the deformation of pre-treatment tumor outline and overlaying on the intratreatment images was less than 5 sec. Thus, the total time required is about 5 min to perform both rigid and deformable registration algorithms and compute the deformed registered image. 5. BREAKTHROUG WORK TO BE PRESENTED A novel linear elastic registration method is presented for deformable co-registration of 3D images. An earlier version of the proposed method was evaluated in the registration of breast MR images. 8, 9 In this paper, the method is further developed and generalized. Also, a different similarity measure is employed in the registration. Moreover, the method is evaluated in the registration of prostate MR images. Proc. of SPIE Vol F-5
6 frequency PZ CG frequency PZ CG TRE (mm) TRE (mm) (a) (b) Figure 4: Frequency distribution of TREs in CG and PZ after (a) rigid alignment and (b) deformable registration. 6. CONCLUSIONS A combined rigid and deformable registration approach was proposed to match pre- and intra-treatment MR images in the focal ablation therapy of prostate cancer. The proposed elastic registration method estimates the tissue deformation using a dynamic linear elastic FE-based model through a filtering process. The registration approach involves minimal manual intervention of the user and does not require segmentation of the prostate boundaries. Experimental results showed that the deformable registration improves the quality of matching especially in the PZ where the endorectal coil deforms the prostate during different imaging sessions. ACKNOWLEDGMENTS The authors gratefully acknowledge funding from the Ontario Institute for Cancer Research (OICR) and Canadian Institutes of Health Research (CIHR). A. Fenster holds a Canada Research Chair in Biomedical Engineering, and acknowledges the support of the Canada Research Chair Program. The authors would like to thank Dr. Jeremy Cepek for collecting and analyzing prostate MR image data. REFERENCES [1] Jemal, A., Bray, F., Center, M. M., Ferlay, J., Ward, E., and Forman, D., Global cancer statistics, CA: a cancer journal for clinicians 61(2), 69 9 (211). [2] Hata, N., Jinzaki, M., Kacher, D., Cormak, R., Gering, D., Nabavi, A., Silverman, S. G., DAmico, A. V., Kikinis, R., Jolesz, F. A., et al., MR imaging-guided prostate biopsy with surgical navigation software: Device validation and feasibility, Radiology 22(1), (21). [3] Raz, O., Haider, M. A., Davidson, S. R., Lindner, U., Hlasny, E., Weersink, R., Gertner, M. R., Kucharcyzk, W., McCluskey, S. A., and Trachtenberg, J., Real-time magnetic resonance imaging guided focal laser therapy in patients with low-risk prostate cancer, European urology 58(1), (21). [4] Lindner, U., Lawrentschuk, N., Weersink, R. A., Davidson, S. R., Raz, O., Hlasny, E., Langer, D. L., Gertner, M. R., Van der Kwast, T., Haider, M. A., and Trachtenberg, J., Focal laser ablation for prostate cancer followed by radical prostatectomy: validation of focal therapy and imaging accuracy, European urology 57(6), (21). [5] Heinrich, M. P., Jenkinson, M., Matin, T., Brady, S. M., Schnabel, J. A., and et al., MIND: Modality independent neighbourhood descriptor for multi-modal deformable registration, Med. Im. Anal. (212). [6] Bathe, K., [Finite Element Procedures], Prentice Hall, Englewood Cliffs, NJ (1996). [7] Welch, G. and Bishop, G., An introduction to the kalman filter, (1995). [8] Marami, B., Sirouspour, S., and Capson, D. W., Model-based deformable registration of preoperative 3D to intraoperative low-resolution 3D and 2D sequences of MR images, in [Medical Image Computing and Computer-Assisted Intervention MICCAI 211], , Springer (211). [9] Marami, B., Elastic Registration of Medical Images Using Generic Dynamic Deformation Models, PhD Thesis, McMaster University Library (213). Proc. of SPIE Vol F-6
Non-rigid MRI-TRUS Registration in Targeted Prostate Biopsy
Non-rigid MRI-TRUS Registration in Targeted Prostate Biopsy Bahram Marami a,b, Shahin Sirouspour c, Suha Ghoul a, Shadi Emami Abarghouei c, Yue Sun a,b and Aaron Fenster a,b,d a Imaging Research Laboratories,
More informationIntraoperative Prostate Tracking with Slice-to-Volume Registration in MR
Intraoperative Prostate Tracking with Slice-to-Volume Registration in MR Sean Gill a, Purang Abolmaesumi a,b, Siddharth Vikal a, Parvin Mousavi a and Gabor Fichtinger a,b,* (a) School of Computing, Queen
More informationA fast breast nonlinear elastography reconstruction technique using the Veronda-Westman model
A fast breast nonlinear elastography reconstruction technique using the Veronda-Westman model Mohammadhosein Amooshahi a and Abbas Samani abc a Department of Electrical & Computer Engineering, University
More informationAtlas Based Segmentation of the prostate in MR images
Atlas Based Segmentation of the prostate in MR images Albert Gubern-Merida and Robert Marti Universitat de Girona, Computer Vision and Robotics Group, Girona, Spain {agubern,marly}@eia.udg.edu Abstract.
More informationMethodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion
Methodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion Mattias P. Heinrich Julia A. Schnabel, Mark Jenkinson, Sir Michael Brady 2 Clinical
More informationAn Automated Image-based Method for Multi-Leaf Collimator Positioning Verification in Intensity Modulated Radiation Therapy
An Automated Image-based Method for Multi-Leaf Collimator Positioning Verification in Intensity Modulated Radiation Therapy Chenyang Xu 1, Siemens Corporate Research, Inc., Princeton, NJ, USA Xiaolei Huang,
More informationA Non-Linear Image Registration Scheme for Real-Time Liver Ultrasound Tracking using Normalized Gradient Fields
A Non-Linear Image Registration Scheme for Real-Time Liver Ultrasound Tracking using Normalized Gradient Fields Lars König, Till Kipshagen and Jan Rühaak Fraunhofer MEVIS Project Group Image Registration,
More informationImproved Non-rigid Registration of Prostate MRI
Improved Non-rigid Registration of Prostate MRI Aloys du Bois d Aische 1,2,3, Mathieu De Craene 1,2,3, Steven Haker 2, Neil Weisenfeld 2,3, Clare Tempany 2, Benoit Macq 1, and Simon K. Warfield 2,3 1 Communications
More informationREAL-TIME ADAPTIVITY IN HEAD-AND-NECK AND LUNG CANCER RADIOTHERAPY IN A GPU ENVIRONMENT
REAL-TIME ADAPTIVITY IN HEAD-AND-NECK AND LUNG CANCER RADIOTHERAPY IN A GPU ENVIRONMENT Anand P Santhanam Assistant Professor, Department of Radiation Oncology OUTLINE Adaptive radiotherapy for head and
More informationOverview of Proposed TG-132 Recommendations
Overview of Proposed TG-132 Recommendations Kristy K Brock, Ph.D., DABR Associate Professor Department of Radiation Oncology, University of Michigan Chair, AAPM TG 132: Image Registration and Fusion Conflict
More informationA Combined Statistical and Biomechanical Model for Estimation of Intra-operative Prostate Deformation
A Combined Statistical and Biomechanical Model for Estimation of Intra-operative Prostate Deformation Ashraf Mohamed 1,2, Christos Davatzikos 1,2, and Russell Taylor 1 1 CISST NSF Engineering Research
More informationDEFORMATION MODELING OF BREAST BASED ON MRI AND FEM.
Sandomierz, 5 th 9 th September 2016 DEFORMATION MODELING OF BREAST BASED ON MRI AND FEM. Marta Danch-Wierzchowska 1,, Damian Borys 1, Andrzej Swierniak 1 1 Institute of Automatic Control, Silesian University
More informationdoi: /
Yiting Xie ; Anthony P. Reeves; Single 3D cell segmentation from optical CT microscope images. Proc. SPIE 934, Medical Imaging 214: Image Processing, 9343B (March 21, 214); doi:1.1117/12.243852. (214)
More informationMedicale Image Analysis
Medicale Image Analysis Registration Validation Prof. Dr. Philippe Cattin MIAC, University of Basel Prof. Dr. Philippe Cattin: Registration Validation Contents 1 Validation 1.1 Validation of Registration
More informationNon-rigid Registration of Preprocedural MRI and Intra-procedural CT for CT-guided Cryoablation Therapy of Liver Cancer
NA-MIC Non-rigid Registration of Preprocedural MRI and Intra-procedural CT for CT-guided Cryoablation Therapy of Liver Cancer Atsushi Yamada, Dominik S. Meier and Nobuhiko Hata Brigham and Women s Hospital
More informationBuilding tools for image-guided adaptive radiotherapy of bladder cancer Chai, X.
UvA-DARE (Digital Academic Repository) Building tools for image-guided adaptive radiotherapy of bladder cancer Chai, X. Link to publication Citation for published version (APA): Chai, X. (2012). Building
More information3D imaging with SPACE vs 2D TSE in MR guided prostate biopsy
3D imaging with SPACE vs 2D TSE in MR guided prostate biopsy Poster No.: C-2004 Congress: ECR 2011 Type: Authors: Keywords: DOI: Scientific Paper M. Garmer, S. Mateiescu, M. Busch, D. Groenemeyer; Bochum/
More informationNonrigid Liver Registration for Image-Guided Surgery using Partial Surface Data: A Novel Iterative Approach
Nonrigid Liver Registration for Image-Guided Surgery using Partial Surface Data: A Novel Iterative Approach D. Caleb Rucker a, Yifei Wu a, Janet E. Ondrake a, Thomas S. Pheiffer a, Amber L. Simpson a,
More informationMaterial Properties Estimation of Layered Soft Tissue Based on MR Observation and Iterative FE Simulation
Material Properties Estimation of Layered Soft Tissue Based on MR Observation and Iterative FE Simulation Mitsunori Tada 1,2, Noritaka Nagai 3, and Takashi Maeno 3 1 Digital Human Research Center, National
More informationMR to Ultrasound Registration for Image-Guided Prostate Biopsy
Western University Scholarship@Western Electronic Thesis and Dissertation Repository May 2014 MR to Ultrasound Registration for Image-Guided Prostate Biopsy Yue Sun The University of Western Ontario Supervisor
More informationAssociation between pathology and texture features of multi parametric MRI of the prostate
Association between pathology and texture features of multi parametric MRI of the prostate 1,2 Peter Kuess, 3 D. Nilsson, 1,2 P. Andrzejewski, 2,4 P. Georg, 1 J. Knoth, 5 M. Susani, 3 J. Trygg, 2,6 T.
More informationMulti-slice-to-volume registration for MRI-guided transperineal prostate biopsy
DOI 10.1007/s11548-014-1108-7 ORIGINAL ARTICLE Multi-slice-to-volume registration for MRI-guided transperineal prostate biopsy Helen Xu Andras Lasso Andriy Fedorov Kemal Tuncali Clare Tempany Gabor Fichtinger
More informationTG 132: Use of Image Registration and Fusion in RT
TG 132: Use of Image Registration and Fusion in RT Kristy K Brock, PhD, DABR, FAAPM Associate Professor Department of Radiation Oncology, University of Michigan Chair, AAPM TG 132: Image Registration and
More informationPredicting Tumour Location by Simulating Large Deformations of the Breast using a 3D Finite Element Model and Nonlinear Elasticity
Predicting Tumour Location by Simulating Large Deformations of the Breast using a 3D Finite Element Model and Nonlinear Elasticity P. Pathmanathan 1 D. Gavaghan 1 J. Whiteley 1 M. Brady 2 M. Nash 3 P.
More informationFINITE element analysis is a powerful computational tool
1 Automated Finite Element Analysis for Deformable Registration of Prostate Images Jessica R. Crouch, Member, IEEE, Stephen M. Pizer, Senior Member, IEEE, Edward L. Chaney, Yu-Chi Hu, Gig S. Mageras, Marco
More informationTarget Motion Tracking in MRI-guided Transrectal Robotic Prostate Biopsy
Target Motion Tracking in MRI-guided Transrectal Robotic Prostate Biopsy Hadi Tadayyon, Member, IEEE, Andras Lasso, Member, IEEE, Aradhana Kaushal, Peter Guion, Member, IEEE, and Gabor Fichtinger, Member,
More informationSoftware Strategy for Robotic Transperineal Prostate Therapy in Closed-Bore MRI
Software Strategy for Robotic Transperineal Prostate Therapy in Closed-Bore MRI Junichi Tokuda 1,GregoryS.Fischer 2,CsabaCsoma 2, Simon P. DiMaio 3, David G. Gobbi 4, Gabor Fichtinger 4,ClareM.Tempany
More informationBasic principles of MR image analysis. Basic principles of MR image analysis. Basic principles of MR image analysis
Basic principles of MR image analysis Basic principles of MR image analysis Julien Milles Leiden University Medical Center Terminology of fmri Brain extraction Registration Linear registration Non-linear
More informationRADIOMICS: potential role in the clinics and challenges
27 giugno 2018 Dipartimento di Fisica Università degli Studi di Milano RADIOMICS: potential role in the clinics and challenges Dr. Francesca Botta Medical Physicist Istituto Europeo di Oncologia (Milano)
More informationLandmark-Guided Surface Matching and Volumetric Warping for Improved Prostate Biopsy Targeting and Guidance
Landmark-Guided Surface Matching and Volumetric Warping for Improved Prostate Biopsy Targeting and Guidance Steven Haker, Simon K. Warfield, and Clare M.C. Tempany Surgical Planning Laboratory Harvard
More information2D Rigid Registration of MR Scans using the 1d Binary Projections
2D Rigid Registration of MR Scans using the 1d Binary Projections Panos D. Kotsas Abstract This paper presents the application of a signal intensity independent registration criterion for 2D rigid body
More informationDeformable Registration Using Scale Space Keypoints
Deformable Registration Using Scale Space Keypoints Mehdi Moradi a, Purang Abolmaesoumi a,b and Parvin Mousavi a a School of Computing, Queen s University, Kingston, Ontario, Canada K7L 3N6; b Department
More informationAbbie M. Diak, PhD Loyola University Medical Center Dept. of Radiation Oncology
Abbie M. Diak, PhD Loyola University Medical Center Dept. of Radiation Oncology Outline High Spectral and Spatial Resolution MR Imaging (HiSS) What it is How to do it Ways to use it HiSS for Radiation
More informationImplementation of Advanced Image Guided Radiation Therapy
Image Acquisition Course Outline Principles, characteristics& applications of the available modalities Image Processing in the T x room Image guided treatment delivery What can / can t we do in the room
More information3D Volume Mesh Generation of Human Organs Using Surface Geometries Created from the Visible Human Data Set
3D Volume Mesh Generation of Human Organs Using Surface Geometries Created from the Visible Human Data Set John M. Sullivan, Jr., Ziji Wu, and Anand Kulkarni Worcester Polytechnic Institute Worcester,
More informationFunctional MRI in Clinical Research and Practice Preprocessing
Functional MRI in Clinical Research and Practice Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization
More informationFinite Element Simulation of Moving Targets in Radio Therapy
Finite Element Simulation of Moving Targets in Radio Therapy Pan Li, Gregor Remmert, Jürgen Biederer, Rolf Bendl Medical Physics, German Cancer Research Center, 69120 Heidelberg Email: pan.li@dkfz.de Abstract.
More informationEPI Data Are Acquired Serially. EPI Data Are Acquired Serially 10/23/2011. Functional Connectivity Preprocessing. fmri Preprocessing
Functional Connectivity Preprocessing Geometric distortion Head motion Geometric distortion Head motion EPI Data Are Acquired Serially EPI Data Are Acquired Serially descending 1 EPI Data Are Acquired
More informationBasic fmri Design and Analysis. Preprocessing
Basic fmri Design and Analysis Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial filtering
More informationMedially Based Meshing with Finite Element Analysis of Prostate Deformation
Medially Based Meshing with Finite Element Analysis of Prostate Deformation Jessica R. Crouch 1, Stephen M. Pizer 1, Edward L. Chaney 1, and Marco Zaider 2 1 Medical Image Display & Analysis Group, University
More informationInteractive Differential Segmentation of the Prostate using Graph-Cuts with a Feature Detector-based Boundary Term
MOSCHIDIS, GRAHAM: GRAPH-CUTS WITH FEATURE DETECTORS 1 Interactive Differential Segmentation of the Prostate using Graph-Cuts with a Feature Detector-based Boundary Term Emmanouil Moschidis emmanouil.moschidis@postgrad.manchester.ac.uk
More informationAutomatic 3D Segmentation of Ultrasound Images Using Atlas Registration and Statistical Texture Prior
Cum Laude Poster Award Automatic 3D Segmentation of Ultrasound Images Using Atlas Registration and Statistical Texture Prior Xiaofeng Yang 1, David Schuster 1, Viraj Master 2, Peter Nieh 2, Aaron Fenster
More informationA high accuracy multi-image registration method for tracking MRIguided
A high accuracy multi-image registration method for tracking MRIguided robots Weijian Shang* a, Gregory S. Fischer* a a Worcester Polytechnic Institute, 100 Institute Road, Worcester, MA, USA 01609; ABSTRACT
More informationMRI-Guided Prostate Motion Tracking using Multislice-to-Volume Registration
MRI-Guided Prostate Motion Tracking using Multislice-to-Volume Registration by Hadi Tadayyon A thesis submitted to the Department of Electrical and Computer Engineering in conformity with the requirements
More informationSPM8 for Basic and Clinical Investigators. Preprocessing. fmri Preprocessing
SPM8 for Basic and Clinical Investigators Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial
More informationA Radiometry Tolerant Method for Direct 3D/2D Registration of Computed Tomography Data to X-ray Images
A Radiometry Tolerant Method for Direct 3D/2D Registration of Computed Tomography Data to X-ray Images Transfer Function Independent Registration Boris Peter Selby 1, Georgios Sakas 2, Stefan Walter 1,
More informationUsing Probability Maps for Multi organ Automatic Segmentation
Using Probability Maps for Multi organ Automatic Segmentation Ranveer Joyseeree 1,2, Óscar Jiménez del Toro1, and Henning Müller 1,3 1 University of Applied Sciences Western Switzerland (HES SO), Sierre,
More informationRIGID IMAGE REGISTRATION
RIGID IMAGE REGISTRATION Duygu Tosun-Turgut, Ph.D. Center for Imaging of Neurodegenerative Diseases Department of Radiology and Biomedical Imaging duygu.tosun@ucsf.edu What is registration? Image registration
More informationSimultaneous Model-based Segmentation of Multiple Objects
Simultaneous Model-based Segmentation of Multiple Objects Astrid Franz 1, Robin Wolz 1, Tobias Klinder 1,2, Cristian Lorenz 1, Hans Barschdorf 1, Thomas Blaffert 1, Sebastian P. M. Dries 1, Steffen Renisch
More informationWhole Body MRI Intensity Standardization
Whole Body MRI Intensity Standardization Florian Jäger 1, László Nyúl 1, Bernd Frericks 2, Frank Wacker 2 and Joachim Hornegger 1 1 Institute of Pattern Recognition, University of Erlangen, {jaeger,nyul,hornegger}@informatik.uni-erlangen.de
More informationMR-Guided Mixed Reality for Breast Conserving Surgical Planning
MR-Guided Mixed Reality for Breast Conserving Surgical Planning Suba Srinivasan (subashini7@gmail.com) March 30 th 2017 Mentors: Prof. Brian A. Hargreaves, Prof. Bruce L. Daniel MEDICINE MRI Guided Mixed
More informationHybrid Spline-based Multimodal Registration using a Local Measure for Mutual Information
Hybrid Spline-based Multimodal Registration using a Local Measure for Mutual Information Andreas Biesdorf 1, Stefan Wörz 1, Hans-Jürgen Kaiser 2, Karl Rohr 1 1 University of Heidelberg, BIOQUANT, IPMB,
More informationSPM8 for Basic and Clinical Investigators. Preprocessing
SPM8 for Basic and Clinical Investigators Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial
More informationFast CT-CT Fluoroscopy Registration with Respiratory Motion Compensation for Image-Guided Lung Intervention
Fast CT-CT Fluoroscopy Registration with Respiratory Motion Compensation for Image-Guided Lung Intervention Po Su a,b, Zhong Xue b*, Kongkuo Lu c, Jianhua Yang a, Stephen T. Wong b a School of Automation,
More informationMEDICAL IMAGE ANALYSIS
SECOND EDITION MEDICAL IMAGE ANALYSIS ATAM P. DHAWAN g, A B IEEE Engineering in Medicine and Biology Society, Sponsor IEEE Press Series in Biomedical Engineering Metin Akay, Series Editor +IEEE IEEE PRESS
More informationABSTRACT 1. INTRODUCTION
Image to physical space registration of supine breast MRI for image guided breast surgery Rebekah H. Conley a, Ingrid M. Meszoely b,thomas S. Pheiffer a, Jared A. Weis c, Thomas E. Yankeelov a,c,d,e, Michael
More informationFINITE element analysis is a powerful computational tool
1 Automated Finite Element Analysis for Deformable Registration of Prostate Images Jessica R. Crouch, Member, IEEE, Stephen M. Pizer, Senior Member, IEEE, Edward L. Chaney, Yu-Chi Hu, Gig S. Mageras, Marco
More informationNIH Public Access Author Manuscript Med Image Comput Comput Assist Interv. Author manuscript; available in PMC 2013 July 16.
NIH Public Access Author Manuscript Published in final edited form as: Med Image Comput Comput Assist Interv. 2012 ; 15(0 3): 107 114. Selection of Optimal Hyper-Parameters for Estimation of Uncertainty
More informationNIH Public Access Author Manuscript Proc Soc Photo Opt Instrum Eng. Author manuscript; available in PMC 2014 October 07.
NIH Public Access Author Manuscript Published in final edited form as: Proc Soc Photo Opt Instrum Eng. 2014 March 21; 9034: 903442. doi:10.1117/12.2042915. MRI Brain Tumor Segmentation and Necrosis Detection
More informationRigid and Deformable Vasculature-to-Image Registration : a Hierarchical Approach
Rigid and Deformable Vasculature-to-Image Registration : a Hierarchical Approach Julien Jomier and Stephen R. Aylward Computer-Aided Diagnosis and Display Lab The University of North Carolina at Chapel
More informationClosed-Loop Control in Fused MR-TRUS Image-Guided Prostate Biopsy
Closed-Loop Control in Fused MR-TRUS Image-Guided Prostate Biopsy Sheng Xu 1, Jochen Kruecker 1, Peter Guion 2, Neil Glossop 3, Ziv Neeman 2, Peter Choyke 2, Anurag K. Singh 2, and Bradford J. Wood 2 1
More informationDigital Volume Correlation for Materials Characterization
19 th World Conference on Non-Destructive Testing 2016 Digital Volume Correlation for Materials Characterization Enrico QUINTANA, Phillip REU, Edward JIMENEZ, Kyle THOMPSON, Sharlotte KRAMER Sandia National
More informationSegmentation of Bony Structures with Ligament Attachment Sites
Segmentation of Bony Structures with Ligament Attachment Sites Heiko Seim 1, Hans Lamecker 1, Markus Heller 2, Stefan Zachow 1 1 Visualisierung und Datenanalyse, Zuse-Institut Berlin (ZIB), 14195 Berlin
More informationPredicting Tumour Location by Modelling the Deformation of the Breast using Nonlinear Elasticity
Predicting Tumour Location by Modelling the Deformation of the Breast using Nonlinear Elasticity November 8th, 2006 Outline Motivation Motivation Motivation for Modelling Breast Deformation Mesh Generation
More informationSuper-resolution Reconstruction of Fetal Brain MRI
Super-resolution Reconstruction of Fetal Brain MRI Ali Gholipour and Simon K. Warfield Computational Radiology Laboratory Children s Hospital Boston, Harvard Medical School Worshop on Image Analysis for
More informationA multi-atlas approach for prostate segmentation in MR images
A multi-atlas approach for prostate segmentation in MR images Geert Litjens, Nico Karssemeijer, and Henkjan Huisman Diagnostic Image Analysis Group, Radboud University Nijmegen Medical Centre, Nijmegen,
More informationNavigation System for ACL Reconstruction Using Registration between Multi-Viewpoint X-ray Images and CT Images
Navigation System for ACL Reconstruction Using Registration between Multi-Viewpoint X-ray Images and CT Images Mamoru Kuga a*, Kazunori Yasuda b, Nobuhiko Hata a, Takeyoshi Dohi a a Graduate School of
More informationLocal Image Registration: An Adaptive Filtering Framework
Local Image Registration: An Adaptive Filtering Framework Gulcin Caner a,a.murattekalp a,b, Gaurav Sharma a and Wendi Heinzelman a a Electrical and Computer Engineering Dept.,University of Rochester, Rochester,
More informationLandmark-based 3D Elastic Registration of Pre- and Postoperative Liver CT Data
Landmark-based 3D Elastic Registration of Pre- and Postoperative Liver CT Data An Experimental Comparison Thomas Lange 1, Stefan Wörz 2, Karl Rohr 2, Peter M. Schlag 3 1 Experimental and Clinical Research
More informationFOREWORD TO THE SPECIAL ISSUE ON MOTION DETECTION AND COMPENSATION
Philips J. Res. 51 (1998) 197-201 FOREWORD TO THE SPECIAL ISSUE ON MOTION DETECTION AND COMPENSATION This special issue of Philips Journalof Research includes a number of papers presented at a Philips
More informationHybrid Formulation of the Model-Based Non-rigid Registration Problem to Improve Accuracy and Robustness
Hybrid Formulation of the Model-Based Non-rigid Registration Problem to Improve Accuracy and Robustness Olivier Clatz 1,2,Hervé Delingette 1, Ion-Florin Talos 2, Alexandra J. Golby 2,RonKikinis 2, Ferenc
More informationEdge-Preserving MRI Super Resolution Using a High Frequency Regularization Technique
Edge-Preserving MRI Super Resolution Using a High Frequency Regularization Technique Kaveh Ahmadi Department of EECS University of Toledo, Toledo, Ohio, USA 43606 Email: Kaveh.ahmadi@utoledo.edu Ezzatollah
More informationEstimating 3D Respiratory Motion from Orbiting Views
Estimating 3D Respiratory Motion from Orbiting Views Rongping Zeng, Jeffrey A. Fessler, James M. Balter The University of Michigan Oct. 2005 Funding provided by NIH Grant P01 CA59827 Motivation Free-breathing
More informationVALIDATION OF DIR. Raj Varadhan, PhD, DABMP Minneapolis Radiation Oncology
VALIDATION OF DIR Raj Varadhan, PhD, DABMP Minneapolis Radiation Oncology Overview Basics: Registration Framework, Theory Discuss Validation techniques Using Synthetic CT data & Phantoms What metrics to
More informationNon-Rigid Multimodal Medical Image Registration using Optical Flow and Gradient Orientation
M. HEINRICH et al.: MULTIMODAL REGISTRATION USING GRADIENT ORIENTATION 1 Non-Rigid Multimodal Medical Image Registration using Optical Flow and Gradient Orientation Mattias P. Heinrich 1 mattias.heinrich@eng.ox.ac.uk
More informationUse of MRI in Radiotherapy: Technical Consideration
Use of MRI in Radiotherapy: Technical Consideration Yanle Hu, PhD Department of Radiation Oncology, Mayo Clinic Arizona 04/07/2018 2015 MFMER slide-1 Conflict of Interest: None 2015 MFMER slide-2 Objectives
More informationBiomedical Image Processing
Biomedical Image Processing Jason Thong Gabriel Grant 1 2 Motivation from the Medical Perspective MRI, CT and other biomedical imaging devices were designed to assist doctors in their diagnosis and treatment
More informationIntroduction to Digitization Techniques for Surgical Guidance
Introduction to Digitization Techniques for Surgical Guidance Rebekah H. Conley, MS and Logan W. Clements, PhD Department of Biomedical Engineering Biomedical Modeling Laboratory Outline Overview of Tracking
More informationTowards Meshless Methods for Surgical Simulation
Ashley Horton, Adam Wittek, Karol Miller Intelligent Systems for Medicine Laboratory School of Mechanical Engineering, The University of Western Australia 35 Stirling Highway, Crawley, Perth WA 6009 Australia
More informationTHE EFFECTS OF THE PLANFORM SHAPE ON DRAG POLAR CURVES OF WINGS: FLUID-STRUCTURE INTERACTION ANALYSES RESULTS
March 18-20, 2013 THE EFFECTS OF THE PLANFORM SHAPE ON DRAG POLAR CURVES OF WINGS: FLUID-STRUCTURE INTERACTION ANALYSES RESULTS Authors: M.R. Chiarelli, M. Ciabattari, M. Cagnoni, G. Lombardi Speaker:
More informationLOCAL-GLOBAL OPTICAL FLOW FOR IMAGE REGISTRATION
LOCAL-GLOBAL OPTICAL FLOW FOR IMAGE REGISTRATION Ammar Zayouna Richard Comley Daming Shi Middlesex University School of Engineering and Information Sciences Middlesex University, London NW4 4BT, UK A.Zayouna@mdx.ac.uk
More informationSpatio-Temporal Registration of Biomedical Images by Computational Methods
Spatio-Temporal Registration of Biomedical Images by Computational Methods Francisco P. M. Oliveira, João Manuel R. S. Tavares tavares@fe.up.pt, www.fe.up.pt/~tavares Outline 1. Introduction 2. Spatial
More informationCoustyx Tutorial Indirect Model
Coustyx Tutorial Indirect Model 1 Introduction This tutorial is created to outline the steps required to compute radiated noise from a gearbox housing using Coustyx software. Detailed steps are given on
More informationPROSTATE CANCER DETECTION USING LABEL IMAGE CONSTRAINED MULTIATLAS SELECTION
PROSTATE CANCER DETECTION USING LABEL IMAGE CONSTRAINED MULTIATLAS SELECTION Ms. Vaibhavi Nandkumar Jagtap 1, Mr. Santosh D. Kale 2 1 PG Scholar, 2 Assistant Professor, Department of Electronics and Telecommunication,
More informationAlternate Biplanar MR Navigation for Microwave Ablation of Liver Tumors
Magnetic Resonance in Medical Sciences, Vol. 4, No. 2, p. 89 94, 2005 TECHNICAL NOTE Alternate Biplanar MR Navigation for Microwave Ablation of Liver Tumors Koichiro SATO 1 *, Shigehiro MORIKAWA 2,ToshiroINUBUSHI
More informationClassification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging
1 CS 9 Final Project Classification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging Feiyu Chen Department of Electrical Engineering ABSTRACT Subject motion is a significant
More informationLucy Phantom MR Grid Evaluation
Lucy Phantom MR Grid Evaluation Anil Sethi, PhD Loyola University Medical Center, Maywood, IL 60153 November 2015 I. Introduction: The MR distortion grid, used as an insert with Lucy 3D QA phantom, is
More information2 Michael E. Leventon and Sarah F. F. Gibson a b c d Fig. 1. (a, b) Two MR scans of a person's knee. Both images have high resolution in-plane, but ha
Model Generation from Multiple Volumes using Constrained Elastic SurfaceNets Michael E. Leventon and Sarah F. F. Gibson 1 MIT Artificial Intelligence Laboratory, Cambridge, MA 02139, USA leventon@ai.mit.edu
More informationABSTRACT 1. INTRODUCTION 2. METHODS
Finding Seeds for Segmentation Using Statistical Fusion Fangxu Xing *a, Andrew J. Asman b, Jerry L. Prince a,c, Bennett A. Landman b,c,d a Department of Electrical and Computer Engineering, Johns Hopkins
More informationMachine Learning for Medical Image Analysis. A. Criminisi
Machine Learning for Medical Image Analysis A. Criminisi Overview Introduction to machine learning Decision forests Applications in medical image analysis Anatomy localization in CT Scans Spine Detection
More informationA Transmission Line Matrix Model for Shielding Effects in Stents
A Transmission Line Matrix Model for hielding Effects in tents Razvan Ciocan (1), Nathan Ida (2) (1) Clemson University Physics and Astronomy Department Clemson University, C 29634-0978 ciocan@clemon.edu
More informationA MORPHOLOGY-BASED FILTER STRUCTURE FOR EDGE-ENHANCING SMOOTHING
Proceedings of the 1994 IEEE International Conference on Image Processing (ICIP-94), pp. 530-534. (Austin, Texas, 13-16 November 1994.) A MORPHOLOGY-BASED FILTER STRUCTURE FOR EDGE-ENHANCING SMOOTHING
More information4D Shape Registration for Dynamic Electrophysiological Cardiac Mapping
4D Shape Registration for Dynamic Electrophysiological Cardiac Mapping Kevin Wilson 1,5, Gerard Guiraudon 4,5, Doug Jones 3,4, and Terry M. Peters 1,2,5 1 Biomedical Engineering Program 2 Department of
More informationValidation of Image Segmentation and Expert Quality with an Expectation-Maximization Algorithm
Validation of Image Segmentation and Expert Quality with an Expectation-Maximization Algorithm Simon K. Warfield, Kelly H. Zou, and William M. Wells Computational Radiology Laboratory and Surgical Planning
More informationIs deformable image registration a solved problem?
Is deformable image registration a solved problem? Marcel van Herk On behalf of the imaging group of the RT department of NKI/AVL Amsterdam, the Netherlands DIR 1 Image registration Find translation.deformation
More informationThe Insight Toolkit. Image Registration Algorithms & Frameworks
The Insight Toolkit Image Registration Algorithms & Frameworks Registration in ITK Image Registration Framework Multi Resolution Registration Framework Components PDE Based Registration FEM Based Registration
More informationDavid Wagner, Kaan Divringi, Can Ozcan Ozen Engineering
Internal Forces of the Femur: An Automated Procedure for Applying Boundary Conditions Obtained From Inverse Dynamic Analysis to Finite Element Simulations David Wagner, Kaan Divringi, Can Ozcan Ozen Engineering
More informationOutline. follows the structure of the report
Outline follows the structure of the report Introduction Mesh-based Modeling of Cuts Finite Element Simulation for Virtual Cutting Numerical Solvers Meshfree Methods Summary & Application Study Discussion
More informationMagnetic Resonance Elastography (MRE) of Liver Disease
Magnetic Resonance Elastography (MRE) of Liver Disease Authored by: Jennifer Dolan Fox, PhD VirtualScopics Inc. jennifer_fox@virtualscopics.com 1-585-249-6231 1. Overview of MRE Imaging MRE is a magnetic
More informationAdditional file 1: Online Supplementary Material 1
Additional file 1: Online Supplementary Material 1 Calyn R Moulton and Michael J House School of Physics, University of Western Australia, Crawley, Western Australia. Victoria Lye, Colin I Tang, Michele
More information