1
|
Liu H, McKenzie E, Xu D, Xu Q, Chin RK, Ruan D, Sheng K. MUsculo-Skeleton-Aware (MUSA) deep learning for anatomically guided head-and-neck CT deformable registration. Med Image Anal 2024; 99:103351. [PMID: 39388843 DOI: 10.1016/j.media.2024.103351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Revised: 06/05/2024] [Accepted: 09/16/2024] [Indexed: 10/12/2024]
Abstract
Deep-learning-based deformable image registration (DL-DIR) has demonstrated improved accuracy compared to time-consuming non-DL methods across various anatomical sites. However, DL-DIR is still challenging in heterogeneous tissue regions with large deformation. In fact, several state-of-the-art DL-DIR methods fail to capture the large, anatomically plausible deformation when tested on head-and-neck computed tomography (CT) images. These results allude to the possibility that such complex head-and-neck deformation may be beyond the capacity of a single network structure or a homogeneous smoothness regularization. To address the challenge of combined multi-scale musculoskeletal motion and soft tissue deformation in the head-and-neck region, we propose a MUsculo-Skeleton-Aware (MUSA) framework to anatomically guide DL-DIR by leveraging the explicit multiresolution strategy and the inhomogeneous deformation constraints between the bony structures and soft tissue. The proposed method decomposes the complex deformation into a bulk posture change and residual fine deformation. It can accommodate both inter- and intra- subject registration. Our results show that the MUSA framework can consistently improve registration accuracy and, more importantly, the plausibility of deformation for various network architectures. The code will be publicly available at https://github.com/HengjieLiu/DIR-MUSA.
Collapse
Affiliation(s)
- Hengjie Liu
- Physics and Biology in Medicine Graduate Program, University of California Los Angeles, Los Angeles, CA, USA; Department of Radiation Oncology, University of California Los Angeles, Los Angeles, CA, USA
| | - Elizabeth McKenzie
- Department of Radiation Oncology, Cedars-Sinai Medical Center, Los Angeles, CA, USA
| | - Di Xu
- UCSF/UC Berkeley Graduate Program in Bioengineering, University of California San Francisco, San Francisco, CA, USA; Department of Radiation Oncology, University of California San Francisco, San Francisco, CA, USA
| | - Qifan Xu
- UCSF/UC Berkeley Graduate Program in Bioengineering, University of California San Francisco, San Francisco, CA, USA; Department of Radiation Oncology, University of California San Francisco, San Francisco, CA, USA
| | - Robert K Chin
- Department of Radiation Oncology, University of California Los Angeles, Los Angeles, CA, USA
| | - Dan Ruan
- Physics and Biology in Medicine Graduate Program, University of California Los Angeles, Los Angeles, CA, USA; Department of Radiation Oncology, University of California Los Angeles, Los Angeles, CA, USA
| | - Ke Sheng
- UCSF/UC Berkeley Graduate Program in Bioengineering, University of California San Francisco, San Francisco, CA, USA; Department of Radiation Oncology, University of California San Francisco, San Francisco, CA, USA.
| |
Collapse
|
2
|
Lobov GI, Kosareva ME. Comparative Characterization of Capsule Mechanical Properties in Mesenteric Lymph Nodes of Young and Aging Bulls. J EVOL BIOCHEM PHYS+ 2022. [DOI: 10.1134/s0022093022050076] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
3
|
Qin A, Chen S, Liang J, Snyder M, Yan D. Evaluation of DIR schemes on tumor/organ with progressive shrinkage: accuracy of tumor/organ internal tissue tracking during the radiation treatment. Radiother Oncol 2022; 173:170-178. [PMID: 35667570 DOI: 10.1016/j.radonc.2022.05.039] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Revised: 05/31/2022] [Accepted: 05/31/2022] [Indexed: 11/19/2022]
Abstract
PURPOSE Accuracy of intratumoral treatment dose accumulation and response assessment highly depends on the accuracy of a DIR method. However, achievable accuracy of the existing DIR methods for tumor/organ with large and progressive shrinkage during the radiotherapy course have not been explored. This study aimed to use a bio-tissue phantom to quantify the achievable accuracy of different DIR schemes. MATERIALS /METHODS A fresh porcine liver was used for phantom material. Sixty gold markers were implanted on the surface and inside of the liver. To simulate the progressive radiation-induced tumor/organ shrinkage, the phantom was heated using a microwave oven incrementally from 30s to 200s in 8 phases. For each phase, the phantom was scanned by CT. Two extra image sets were generated from the original images: 1) the image set with overriding the high-density gold markers (feature image); 2) the image set with overriding the entire phantom to the mean soft tissue intensity (featureless image). Ten DIR schemes were evaluated to mimic clinical treatment situations of tumor/critical organ with respect to their surface and internal condition of image features, availability of intermediate feedback images and DIR methods. The internal marker's positions were utilized to evaluate DIR accuracy quantified by target registration error (TRE). RESULTS Volume reduction was about 20% ∼ 40% of the initial volume after 90s ∼ 200s of the heating. Without image features on the surface and inside of the phantom, the hybrid-DIR (image-based DIR followed by biomechanical model-based refinement) with the surface constraint achieved the registration TRE from 2.6 ± 1.2mm to 5.3 ± 2.6mm proportional to the %volume shrinkage. Meanwhile, the hybrid-DIR with the surface-marker constraint achieved the TRE from 2.4 ± 1.2mm to 2.6 ± 1.0mm. If both the surface and internal image features would be viable on the feedback images, the achievable accuracy could be minimal with the TRE from 1.6±0.9mm to 1.9 ± 1.2mm. CONCLUSIONS Standard DIR methods cannot guarantee intratumoral tissue registration accuracy for tumor/organ with large progressive shrinkage. Achievable accuracy with using the hybrid DIR method is highly dependent on the surface registration accuracy. If the surface registration mean TRE can be controlled within 2mm, the mean TRE of internal tissue can be controlled within 3mm.
Collapse
Affiliation(s)
- An Qin
- Dept. of Radiation Oncology, Beaumont Health System, Royal Oak, United States
| | - Shupeng Chen
- Dept. of Radiation Oncology, Beaumont Health System, Royal Oak, United States
| | - Jian Liang
- Dept. of Radiation Oncology, Beaumont Health System, Royal Oak, United States
| | - Michael Snyder
- Dept. of Radiation Oncology, Beaumont Health System, Royal Oak, United States
| | - Di Yan
- Dept. of Radiation Oncology, Beaumont Health System, Royal Oak, United States; Radiation Oncology, Huaxi Hospitals & Medical School, Chengdu, China.
| |
Collapse
|
4
|
Qin A, Ionascu D, Liang J, Han X, O’Connell N, Yan D. The evaluation of a hybrid biomechanical deformable registration method on a multistage physical phantom with reproducible deformation. Radiat Oncol 2018; 13:240. [PMID: 30514348 PMCID: PMC6280462 DOI: 10.1186/s13014-018-1192-x] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2018] [Accepted: 11/23/2018] [Indexed: 12/22/2022] Open
Abstract
BACKGROUND Advanced clinical applications, such as dose accumulation and adaptive radiation therapy, require deformable image registration (DIR) algorithms capable of voxel-wise accurate mapping of treatment dose or functional imaging. By utilizing a multistage deformable phantom, the authors investigated scenarios where biomechanical refinement method (BM-DIR) may be better than the pure image intensity based deformable registration (IM-DIR). METHODS The authors developed a biomechanical-model based DIR refinement method (BM-DIR) to refine the deformable vector field (DVF) from any initial intensity-based DIR (IM-DIR). The BM-DIR method was quantitatively evaluated on a novel phantom capable of ten reproducible gradually-increasing deformation stages using the urethra tube as a surrogate. The internal DIR accuracy was inspected in term of the Dice similarity coefficient (DSC), Hausdorff and mean surface distance as defined in of the urethra structure inside the phantom and compared with that of the initial IM-DIR under various stages of deformation. Voxel-wise deformation vector discrepancy and Jacobian regularity were also inspected to evaluate the output DVFs. In addition to phantom, two pairs of Head&Neck patient MR images with expert-defined landmarks inside parotids were utilized to evaluate the BM-DIR accuracy with target registration error (TRE). RESULTS The DSC and surface distance measures of the inner urethra tube indicated the BM-DIR method can improve the internal DVF accuracy on masked MR images for the phases of a large degree of deformation. The smoother Jacobian distribution from the BM-DIR suggests more physically-plausible internal deformation. For H&N cancer patients, the BM-DIR improved the TRE from 0.339 cm to 0.210 cm for the landmarks inside parotid on the masked MR images. CONCLUSIONS We have quantitatively demonstrated on a multi-stage physical phantom and limited patient data that the proposed BM-DIR can improve the accuracy inside solid organs with large deformation where distinctive image features are absent.
Collapse
Affiliation(s)
- An Qin
- Department of Radiation Oncology, Beaumont Health System, Royal Oak, MI USA
| | - Dan Ionascu
- Department of Radiation Oncology, College of Medicine, University of Cincinnati, Cincinnati, OH USA
| | - Jian Liang
- Department of Radiation Oncology, Beaumont Health System, Royal Oak, MI USA
| | - Xiao Han
- Elekta Inc., Maryland Heights, MO USA
| | | | - Di Yan
- Department of Radiation Oncology, Beaumont Health System, Royal Oak, MI USA
| |
Collapse
|
5
|
Qin A, Liang J, Han X, O'Connell N, Yan D. Technical Note: The impact of deformable image registration methods on dose warping. Med Phys 2018; 45:1287-1294. [PMID: 29297939 DOI: 10.1002/mp.12741] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2017] [Revised: 11/15/2017] [Accepted: 12/12/2017] [Indexed: 11/08/2022] Open
Abstract
PURPOSE The purpose of this study was to investigate the clinical-relevant discrepancy between doses warped by pure image based deformable image registration (IM-DIR) and by biomechanical model based DIR (BM-DIR) on intensity-homogeneous organs. METHODS AND MATERIALS Ten patients (5Head&Neck, 5Prostate) were included. A research DIR tool (ADMRIE_v1.12) was utilized for IM-DIR. After IM-DIR, BM-DIR was carried out for organs (parotids, bladder, and rectum) which often encompass sharp dose gradient. Briefly, high-quality tetrahedron meshes were generated and deformable vector fields (DVF) from IM-DIR were interpolated to the surface nodes of the volume meshes as boundary condition. Then, a FEM solver (ABAQUS_v6.14) was used to simulate the displacement of internal nodes, which were then interpolated to image-voxel grids to get the more physically plausible DVF. Both geometrical and subsequent dose warping discrepancies were quantified between the two DIR methods. Target registration discrepancy(TRD) was evaluated to show the geometry difference. The re-calculated doses on second CT were warped to the pre-treatment CT via two DIR. Clinical-relevant dose parameters and γ passing rate were compared between two types of warped dose. The correlation was evaluated between parotid shrinkage and TRD/dose discrepancy. RESULT The parotid shrunk to 75.7% ± 9% of its pre-treatment volume and the percentage of volume with TRD>1.5 mm) was 6.5% ± 4.7%. The normalized mean-dose difference (NMDD) of IM-DIR and BM-DIR was -0.8% ± 1.5%, with range (-4.7% to 1.5%). 2 mm/2% passing rate was 99.0% ± 1.4%. A moderate correlation was found between parotid shrinkage and TRD and NMDD. The bladder had a NMDD of -9.9% ± 9.7%, with BM-DIR warped dose systematically higher. Only minor deviation was observed for rectum NMDD (0.5% ± 1.1%). CONCLUSION Impact of DIR method on treatment dose warping is patient and organ-specific. Generally, intensity-homogeneous organs, which undergo larger deformation/shrinkage during treatment and encompass sharp dose gradient, will have greater dose warping uncertainty. For these organs, BM-DIR could be beneficial to the evaluation of DIR/dose-warping uncertainty.
Collapse
Affiliation(s)
- An Qin
- Dept. of Radiation Oncology, Beaumont Health System, Royal Oak, MI, USA
| | - Jian Liang
- Dept. of Radiation Oncology, Beaumont Health System, Royal Oak, MI, USA
| | - Xiao Han
- Elekta Inc., Maryland Heights, MO, 63043, USA
| | | | - Di Yan
- Dept. of Radiation Oncology, Beaumont Health System, Royal Oak, MI, USA
| |
Collapse
|
6
|
McGrath DM, Lee J, Foltz WD, Samavati N, van der Kwast T, Jewett MAS, Chung P, Ménard C, Brock KK. MR elastography to measure the effects of cancer and pathology fixation on prostate biomechanics, and comparison with T 1, T 2 and ADC. Phys Med Biol 2017; 62:1126-1148. [PMID: 28092638 DOI: 10.1088/1361-6560/aa52f4] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
MRI is under evaluation for image-guided intervention for prostate cancer. The sensitivity and specificity of MRI parameters is determined via correlation with the gold-standard of histopathology. Whole-mount histopathology of prostatectomy specimens can be digitally registered to in vivo imaging for correlation. When biomechanical-based deformable registration is employed to account for deformation during histopathology processing, the ex vivo biomechanical properties are required. However, these properties are altered by pathology fixation, and vary with disease. Hence, this study employs magnetic resonance elastography (MRE) to measure ex vivo prostate biomechanical properties before and after fixation. A quasi-static MRE method was employed to measure high resolution maps of Young's modulus (E) before and after fixation of canine prostate and prostatectomy specimens (n = 4) from prostate cancer patients who had previously received radiation therapy. For comparison, T 1, T 2 and apparent diffusion coefficient (ADC) were measured in parallel. E (kPa) varied across clinical anatomy and for histopathology-identified tumor: peripheral zone: 99(±22), central gland: 48(±37), tumor: 85(±53), and increased consistently with fixation (factor of 11 ± 5; p < 0.02). T 2 decreased consistently with fixation, while changes in T 1 and ADC were more complex and inconsistent. The biomechanics of the clinical prostate specimens varied greatly with fixation, and to a lesser extent with disease and anatomy. The data obtained will improve the precision of prostate pathology correlation, leading to more accurate disease detection and targeting.
Collapse
Affiliation(s)
- Deirdre M McGrath
- Radiation Medicine Program, Princess Margaret Hospital, University Health Network, Toronto, Ontario M5G 2M9, Canada
| | | | | | | | | | | | | | | | | |
Collapse
|
7
|
Velec M, Moseley JL, Svensson S, Hårdemark B, Jaffray DA, Brock KK. Validation of biomechanical deformable image registration in the abdomen, thorax, and pelvis in a commercial radiotherapy treatment planning system. Med Phys 2017; 44:3407-3417. [PMID: 28453911 DOI: 10.1002/mp.12307] [Citation(s) in RCA: 49] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2016] [Revised: 03/20/2017] [Accepted: 04/20/2017] [Indexed: 01/13/2023] Open
Abstract
PURPOSE The accuracy of deformable image registration tools can vary widely between imaging modalities and specific implementations of the same algorithms. A biomechanical model-based algorithm initially developed in-house at an academic institution was translated into a commercial radiotherapy treatment planning system and validated for multiple imaging modalities and anatomic sites. METHODS Biomechanical deformable registration (Morfeus) is a geometry-driven algorithm based on the finite element method. Boundary conditions are derived from the model-based segmentation of controlling structures in each image which establishes a point-to-point surface correspondence. For each controlling structure, material properties and fixed or sliding interfaces are assigned. The displacements of internal volumes for controlling structures and other structures implicitly deformed are solved with finite element analysis. Registration was performed for 74 patients with images (mean vector resolution) of thoracic and abdominal 4DCT (2.8 mm) and MR (5.3 mm), liver CT-MR (4.5 mm), and prostate MR (2.6 mm). Accuracy was quantified between deformed and actual target images using distance-to-agreement (DTA) for structure surfaces and the target registration error (TRE) for internal point landmarks. RESULTS The results of the commercial implementation were as follows. The mean DTA was ≤ 1.0 mm for controlling structures and 1.0-3.5 mm for implicitly deformed structures on average. TRE ranged from 2.0 mm on prostate MR to 5.1 mm on lung MR on average, within 0.1 mm or lower than the image voxel sizes. Accuracy was not overly sensitive to changes in the material properties or variability in structure segmentations, as changing these inputs affected DTA and TRE by ≤ 0.8 mm. Maximum DTA > 5 mm occurred for 88% of the structures evaluated although these were within the inherent segmentation uncertainty for 82% of structures. Differences in accuracy between the commercial and in-house research implementations were ≤ 0.5 mm for mean DTA and ≤ 0.7 mm for mean TRE. CONCLUSIONS Accuracy of biomechanical deformable registration evaluated on a large cohort of images in the thorax, abdomen and prostate was similar to the image voxel resolution on average across multiple modalities. Validation of this treatment planning system implementation supports biomechanical deformable registration as a versatile clinical tool to enable accurate target delineation at planning and treatment adaptation.
Collapse
Affiliation(s)
- Michael Velec
- Techna Institute and Princess Margaret Cancer Centre, University Health Network, Toronto, M5G 2M9, Canada
| | - Joanne L Moseley
- Techna Institute and Princess Margaret Cancer Centre, University Health Network, Toronto, M5G 2M9, Canada
| | - Stina Svensson
- RaySearch Laboratories AB, Sveavägen 44, SE-103 65, Stockholm, Sweden
| | - Björn Hårdemark
- RaySearch Laboratories AB, Sveavägen 44, SE-103 65, Stockholm, Sweden
| | - David A Jaffray
- Techna Institute and Princess Margaret Cancer Centre, University Health Network, Toronto, M5G 2M9, Canada.,Department of Radiation Oncology, Medical Biophysics, and Institute of Biomaterials and Biomedical Engineering, University of Toronto, Toronto, M5S 3E2, Canada
| | - Kristy K Brock
- Department of Radiation Oncology, University of Michigan, Ann Arbor, MI, 48109, USA
| |
Collapse
|
8
|
Brock KK, Mutic S, McNutt TR, Li H, Kessler ML. Use of image registration and fusion algorithms and techniques in radiotherapy: Report of the AAPM Radiation Therapy Committee Task Group No. 132. Med Phys 2017; 44:e43-e76. [PMID: 28376237 DOI: 10.1002/mp.12256] [Citation(s) in RCA: 518] [Impact Index Per Article: 74.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2016] [Revised: 02/13/2017] [Accepted: 02/19/2017] [Indexed: 11/07/2022] Open
Abstract
Image registration and fusion algorithms exist in almost every software system that creates or uses images in radiotherapy. Most treatment planning systems support some form of image registration and fusion to allow the use of multimodality and time-series image data and even anatomical atlases to assist in target volume and normal tissue delineation. Treatment delivery systems perform registration and fusion between the planning images and the in-room images acquired during the treatment to assist patient positioning. Advanced applications are beginning to support daily dose assessment and enable adaptive radiotherapy using image registration and fusion to propagate contours and accumulate dose between image data taken over the course of therapy to provide up-to-date estimates of anatomical changes and delivered dose. This information aids in the detection of anatomical and functional changes that might elicit changes in the treatment plan or prescription. As the output of the image registration process is always used as the input of another process for planning or delivery, it is important to understand and communicate the uncertainty associated with the software in general and the result of a specific registration. Unfortunately, there is no standard mathematical formalism to perform this for real-world situations where noise, distortion, and complex anatomical variations can occur. Validation of the software systems performance is also complicated by the lack of documentation available from commercial systems leading to use of these systems in undesirable 'black-box' fashion. In view of this situation and the central role that image registration and fusion play in treatment planning and delivery, the Therapy Physics Committee of the American Association of Physicists in Medicine commissioned Task Group 132 to review current approaches and solutions for image registration (both rigid and deformable) in radiotherapy and to provide recommendations for quality assurance and quality control of these clinical processes.
Collapse
Affiliation(s)
- Kristy K Brock
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, 1400 Pressler St, FCT 14.6048, Houston, TX, 77030, USA
| | - Sasa Mutic
- Department of Radiation Oncology, Washington University School of Medicine, St. Louis, MO, USA
| | - Todd R McNutt
- Department of Radiation Oncology, Johns Hopkins Medical Institute, Baltimore, MD, USA
| | - Hua Li
- Department of Radiation Oncology, Washington University School of Medicine, St. Louis, MO, USA
| | - Marc L Kessler
- Department of Radiation Oncology, University of Michigan, Ann Arbor, MI, USA
| |
Collapse
|
9
|
Patient-specific Deformation Modelling via Elastography: Application to Image-guided Prostate Interventions. Sci Rep 2016; 6:27386. [PMID: 27272239 PMCID: PMC4895338 DOI: 10.1038/srep27386] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2016] [Accepted: 05/16/2016] [Indexed: 01/22/2023] Open
Abstract
Image-guided prostate interventions often require the registration of preoperative magnetic resonance (MR) images to real-time transrectal ultrasound (TRUS) images to provide high-quality guidance. One of the main challenges for registering MR images to TRUS images is how to estimate the TRUS-probe-induced prostate deformation that occurs during TRUS imaging. The combined statistical and biomechanical modeling approach shows promise for the adequate estimation of prostate deformation. However, the right setting of the biomechanical parameters is very crucial for realistic deformation modeling. We propose a patient-specific deformation model equipped with personalized biomechanical parameters obtained from shear wave elastography to reliably predict the prostate deformation during image-guided interventions. Using data acquired from a prostate phantom and twelve patients with suspected prostate cancer, we compared the prostate deformation model with and without patient-specific biomechanical parameters in terms of deformation estimation accuracy. The results show that the patient-specific deformation model possesses favorable model ability, and outperforms the model without patient-specific biomechanical parameters. The employment of the patient-specific biomechanical parameters obtained from elastography for deformation modeling shows promise for providing more precise deformation estimation in applications that use computer-assisted image-guided intervention systems.
Collapse
|
10
|
Toth R, Sperling D, Madabhushi A. Quantifying Post- Laser Ablation Prostate Therapy Changes on MRI via a Domain-Specific Biomechanical Model: Preliminary Findings. PLoS One 2016; 11:e0150016. [PMID: 27088600 PMCID: PMC4835053 DOI: 10.1371/journal.pone.0150016] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2015] [Accepted: 02/08/2016] [Indexed: 11/18/2022] Open
Abstract
Focal laser ablation destroys cancerous cells via thermal destruction of tissue by a laser. Heat is absorbed, causing thermal necrosis of the target region. It combines the aggressive benefits of radiation treatment (destroying cancer cells) without the harmful side effects (due to its precise localization). MRI is typically used pre-treatment to determine the targeted area, and post-treatment to determine efficacy by detecting necrotic tissue, or tumor recurrence. However, no system exists to quantitatively evaluate the post-treatment effects on the morphology and structure via MRI. To quantify these changes, the pre- and post-treatment MR images must first be spatially aligned. The goal is to quantify (a) laser-induced shape-based changes, and (b) changes in MRI parameters post-treatment. The shape-based changes may be correlated with treatment efficacy, and the quantitative effects of laser treatment over time is currently poorly understood. This work attempts to model changes in gland morphology following laser treatment due to (1) patient alignment, (2) changes due to surrounding organs such as the bladder and rectum, and (3) changes due to the treatment itself. To isolate the treatment-induced shape-based changes, the changes from (1) and (2) are first modeled and removed using a finite element model (FEM). A FEM models the physical properties of tissue. The use of a physical biomechanical model is important since a stated goal of this work is to determine the physical shape-based changes to the prostate from the treatment, and therefore only physical real deformations are to be allowed. A second FEM is then used to isolate the physical, shape-based, treatment-induced changes. We applied and evaluated our model in capturing the laser induced changes to the prostate morphology on eight patients with 3.0 Tesla, T2-weighted MRI, acquired approximately six months following treatment. Our results suggest the laser treatment causes a decrease in prostate volume, which appears to manifest predominantly at the site of ablation. After spatially aligning the images, changes to MRI intensity values are clearly visible at the site of ablation. Our results suggest that our new methodology is able to capture and quantify the degree of laser-induced changes to the prostate. The quantitative measurements reflecting of the deformation changes can be used to track treatment response over time.
Collapse
Affiliation(s)
- Robert Toth
- Department of Biomedical Engineering, Case Western Reserve University, Cleveland, OH, United States of America
- Toth Technology LLC, Long Valley, NJ, United States of America
- * E-mail:
| | - Dan Sperling
- Sperling Prostate Center, Manhattan, NY, United States of America
| | - Anant Madabhushi
- Department of Biomedical Engineering, Case Western Reserve University, Cleveland, OH, United States of America
| |
Collapse
|
11
|
Rosewall T, Wheat J, Currie G, Kong V, Bayley AJ, Moseley J, Chung P, Catton C, Craig T, Milosevic M. Planned versus 'delivered' bladder dose reconstructed using solid and hollow organ models during prostate cancer IMRT. Radiother Oncol 2016; 119:417-22. [PMID: 27072936 DOI: 10.1016/j.radonc.2016.03.022] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2016] [Revised: 03/11/2016] [Accepted: 03/13/2016] [Indexed: 10/22/2022]
Abstract
BACKGROUND AND PURPOSE All studies to date have evaluated the dosimetric effect of bladder deformation using an organ model that includes the dose to the urine. This research reconstructed bladder dose using both hollow and solid organ models, to determine if dose/volume differences exist. MATERIALS AND METHODS 35 prostate IMRT patients were selected, who had received 78Gy in 39 fractions and full bladder instructions. Biomechanical modelling and finite element analysis were used to reconstruct bladder dose (solid and hollow organ model) using every third CBCT throughout the treatment course. RESULTS Reconstructed dose (ReconDose) was 11.3Gy greater than planned dose (planDose) with a hollow bladder model (p<0.001) and 12.3Gy greater with a solid bladder model (p<0.0001). Median reconstructed volumes within the 30Gy, 65Gy and 78Gy isodoses were 3-4 times larger with the solid organ model (p<0.0001). The difference between planning bladder volume and median treatment volume was associated with the difference between the planDose and reconDose below 78Gy (R(2)>0.61). CONCLUSIONS Substantial differences exist between planned and reconstructed bladder dose, associated with the differences in bladder filling between planning and treatment. Dose reconstructed using a solid bladder model over-reports the volume of bladder within key isodose levels and overestimates the differences between planned and reconstructed dose. Dose reconstruction with a hollow organ model is recommended if the goal is to associate that dose with toxicity.
Collapse
Affiliation(s)
- Tara Rosewall
- Radiation Medicine Program, Princess Margaret Cancer Centre, Toronto, Canada; Department of Radiation Oncology, University of Toronto, Canada; Faculty of Science, Charles Sturt University, Australia.
| | - Janelle Wheat
- Faculty of Science, Charles Sturt University, Australia; Faculty of Medicine and Health Sciences, Macquarie University, Australia
| | - Geoffrey Currie
- Faculty of Science, Charles Sturt University, Australia; Faculty of Medicine and Health Sciences, Macquarie University, Australia
| | - Vickie Kong
- Radiation Medicine Program, Princess Margaret Cancer Centre, Toronto, Canada
| | - Andrew J Bayley
- Radiation Medicine Program, Princess Margaret Cancer Centre, Toronto, Canada; Department of Radiation Oncology, University of Toronto, Canada
| | - Joanne Moseley
- Radiation Medicine Program, Princess Margaret Cancer Centre, Toronto, Canada
| | - Peter Chung
- Radiation Medicine Program, Princess Margaret Cancer Centre, Toronto, Canada; Department of Radiation Oncology, University of Toronto, Canada
| | - Charles Catton
- Radiation Medicine Program, Princess Margaret Cancer Centre, Toronto, Canada; Department of Radiation Oncology, University of Toronto, Canada
| | - Tim Craig
- Radiation Medicine Program, Princess Margaret Cancer Centre, Toronto, Canada; Department of Radiation Oncology, University of Toronto, Canada
| | - Michael Milosevic
- Radiation Medicine Program, Princess Margaret Cancer Centre, Toronto, Canada; Department of Radiation Oncology, University of Toronto, Canada
| |
Collapse
|
12
|
van de Ven WJM, Hu Y, Barentsz JO, Karssemeijer N, Barratt D, Huisman HJ. Biomechanical modeling constrained surface-based image registration for prostate MR guided TRUS biopsy. Med Phys 2016; 42:2470-81. [PMID: 25979040 DOI: 10.1118/1.4917481] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
PURPOSE Adding magnetic resonance (MR)-derived information to standard transrectal ultrasound (TRUS) images for guiding prostate biopsy is of substantial clinical interest. A tumor visible on MR images can be projected on ultrasound (US) by using MR-US registration. A common approach is to use surface-based registration. The authors hypothesize that biomechanical modeling will better control deformation inside the prostate than a regular nonrigid surface-based registration method. The authors developed a novel method by extending a nonrigid surface-based registration algorithm with biomechanical finite element (FE) modeling to better predict internal deformations of the prostate. METHODS Data were collected from ten patients and the MR and TRUS images were rigidly registered to anatomically align prostate orientations. The prostate was manually segmented in both images and corresponding surface meshes were generated. Next, a tetrahedral volume mesh was generated from the MR image. Prostate deformations due to the TRUS probe were simulated using the surface displacements as the boundary condition. A three-dimensional thin-plate spline deformation field was calculated by registering the mesh vertices. The target registration errors (TREs) of 35 reference landmarks determined by surface and volume mesh registrations were compared. RESULTS The median TRE of a surface-based registration with biomechanical regularization was 2.76 (0.81-7.96) mm. This was significantly different than the median TRE of 3.47 (1.05-7.80) mm for regular surface-based registration without biomechanical regularization. CONCLUSIONS Biomechanical FE modeling has the potential to improve the accuracy of multimodal prostate registration when comparing it to a regular nonrigid surface-based registration algorithm and can help to improve the effectiveness of MR guided TRUS biopsy procedures.
Collapse
Affiliation(s)
- Wendy J M van de Ven
- Department of Radiology and Nuclear Medicine, Radboud University Medical Center, Nijmegen 6525 GA, The Netherlands
| | - Yipeng Hu
- Centre for Medical Image Computing, University College London, London WC1E 6BT, United Kingdom
| | - Jelle O Barentsz
- Department of Radiology and Nuclear Medicine, Radboud University Medical Center, Nijmegen 6525 GA, The Netherlands
| | - Nico Karssemeijer
- Department of Radiology and Nuclear Medicine, Radboud University Medical Center, Nijmegen 6525 GA, The Netherlands
| | - Dean Barratt
- Centre for Medical Image Computing, University College London, London WC1E 6BT, United Kingdom
| | - Henkjan J Huisman
- Department of Radiology and Nuclear Medicine, Radboud University Medical Center, Nijmegen 6525 GA, The Netherlands
| |
Collapse
|
13
|
Wang Y, Cheng JZ, Ni D, Lin M, Qin J, Luo X, Xu M, Xie X, Heng PA. Towards Personalized Statistical Deformable Model and Hybrid Point Matching for Robust MR-TRUS Registration. IEEE TRANSACTIONS ON MEDICAL IMAGING 2016; 35:589-604. [PMID: 26441446 DOI: 10.1109/tmi.2015.2485299] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Registration and fusion of magnetic resonance (MR) and 3D transrectal ultrasound (TRUS) images of the prostate gland can provide high-quality guidance for prostate interventions. However, accurate MR-TRUS registration remains a challenging task, due to the great intensity variation between two modalities, the lack of intrinsic fiducials within the prostate, the large gland deformation caused by the TRUS probe insertion, and distinctive biomechanical properties in patients and prostate zones. To address these challenges, a personalized model-to-surface registration approach is proposed in this study. The main contributions of this paper can be threefold. First, a new personalized statistical deformable model (PSDM) is proposed with the finite element analysis and the patient-specific tissue parameters measured from the ultrasound elastography. Second, a hybrid point matching method is developed by introducing the modality independent neighborhood descriptor (MIND) to weight the Euclidean distance between points to establish reliable surface point correspondence. Third, the hybrid point matching is further guided by the PSDM for more physically plausible deformation estimation. Eighteen sets of patient data are included to test the efficacy of the proposed method. The experimental results demonstrate that our approach provides more accurate and robust MR-TRUS registration than state-of-the-art methods do. The averaged target registration error is 1.44 mm, which meets the clinical requirement of 1.9 mm for the accurate tumor volume detection. It can be concluded that the presented method can effectively fuse the heterogeneous image information in the elastography, MR, and TRUS to attain satisfactory image alignment performance.
Collapse
|
14
|
Multi-axis dose accumulation of noninvasive image-guided breast brachytherapy through biomechanical modeling of tissue deformation using the finite element method. J Contemp Brachytherapy 2015; 7:55-71. [PMID: 25829938 PMCID: PMC4371066 DOI: 10.5114/jcb.2015.49355] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2015] [Revised: 01/29/2015] [Accepted: 02/01/2015] [Indexed: 01/08/2023] Open
Abstract
PURPOSE Noninvasive image-guided breast brachytherapy delivers conformal HDR (192)Ir brachytherapy treatments with the breast compressed, and treated in the cranial-caudal and medial-lateral directions. This technique subjects breast tissue to extreme deformations not observed for other disease sites. Given that, commercially-available software for deformable image registration cannot accurately co-register image sets obtained in these two states, a finite element analysis based on a biomechanical model was developed to deform dose distributions for each compression circumstance for dose summation. MATERIAL AND METHODS The model assumed the breast was under planar stress with values of 30 kPa for Young's modulus and 0.3 for Poisson's ratio. Dose distributions from round and skin-dose optimized applicators in cranial-caudal and medial-lateral compressions were deformed using 0.1 cm planar resolution. Dose distributions, skin doses, and dose-volume histograms were generated. Results were examined as a function of breast thickness, applicator size, target size, and offset distance from the center. RESULTS Over the range of examined thicknesses, target size increased several millimeters as compression thickness decreased. This trend increased with increasing offset distances. Applicator size minimally affected target coverage, until applicator size was less than the compressed target size. In all cases, with an applicator larger or equal to the compressed target size, > 90% of the target covered by > 90% of the prescription dose. In all cases, dose coverage became less uniform as offset distance increased and average dose increased. This effect was more pronounced for smaller target-applicator combinations. CONCLUSIONS The model exhibited skin dose trends that matched MC-generated benchmarking results within 2% and clinical observations over a similar range of breast thicknesses and target sizes. The model provided quantitative insight on dosimetric treatment variables over a range of clinical circumstances. These findings highlight the need for careful target localization and accurate identification of compression thickness and target offset.
Collapse
|
15
|
Samavati N, McGrath DM, Jewett MA, van der Kwast T, Ménard C, Brock KK. Effect of material property heterogeneity on biomechanical modeling of prostate under deformation. Phys Med Biol 2015; 60:195-209. [PMID: 25489840 PMCID: PMC4443715 DOI: 10.1088/0031-9155/60/1/195] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Biomechanical model based deformable image registration has been widely used to account for prostate deformation in various medical imaging procedures. Biomechanical material properties are important components of a biomechanical model. In this study, the effect of incorporating tumor-specific material properties in the prostate biomechanical model was investigated to provide insight into the potential impact of material heterogeneity on the prostate deformation calculations. First, a simple spherical prostate and tumor model was used to analytically describe the deformations and demonstrate the fundamental effect of changes in the tumor volume and stiffness in the modeled deformation. Next, using a clinical prostate model, a parametric approach was used to describe the variations in the heterogeneous prostate model by changing tumor volume, stiffness, and location, to show the differences in the modeled deformation between heterogeneous and homogeneous prostate models. Finally, five clinical prostatectomy examples were used in separately performed homogeneous and heterogeneous biomechanical model based registrations to describe the deformations between 3D reconstructed histopathology images and ex vivo magnetic resonance imaging, and examine the potential clinical impact of modeling biomechanical heterogeneity of the prostate. The analytical formulation showed that increasing the tumor volume and stiffness could significantly increase the impact of the heterogeneous prostate model in the calculated displacement differences compared to the homogeneous model. The parametric approach using a single prostate model indicated up to 4.8 mm of displacement difference at the tumor boundary compared to a homogeneous model. Such differences in the deformation of the prostate could be potentially clinically significant given the voxel size of the ex vivo MR images (0.3 × 0.3 × 0.3 mm). However, no significant changes in the registration accuracy were observed using heterogeneous models for the limited number of clinical prostatectomy patients modeled and evaluated in this study.
Collapse
Affiliation(s)
- Navid Samavati
- Institute of Biomaterials and Biomedical Engineering, University of Toronto, Toronto, Ontario M5S 3G9
| | - Deirdre M. McGrath
- Radiation Medicine Program, Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario M5G 2M9, Canada
| | - Michael A.S. Jewett
- Division of Urology, Department of Surgery and Surgical Oncology, Princess Margaret Cancer Centre, University of Toronto, Toronto, Ontario M5G 2C4, Canada
| | - Theo van der Kwast
- Department of Pathology, University Health Network, Toronto, Ontario M5G 2C4, Canada
| | - Cynthia Ménard
- Radiation Medicine Program, Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario M5G 2M9, Canada
- Department of Radiation Oncology, University of Toronto
| | - Kristy K. Brock
- Department of Radiation Oncology, University of Michigan, Ann Arbor, MI, USA, 48109
| |
Collapse
|
16
|
Toth R, Traughber B, Ellis R, Kurhanewicz J, Madabhushi A. A Domain Constrained Deformable (DoCD) Model for Co-registration of Pre- and Post-Radiated Prostate MRI. Neurocomputing 2014; 114:3-12. [PMID: 25267873 DOI: 10.1016/j.neucom.2014.01.058] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
External beam radiation treatment (EBRT) is a popular method for treating prostate cancer (CaP) involving destroying tumor cells with ionizing radiation. Following EBRT, biochemical failure has been linked with disease recurrence. However, there is a need for methods for evaluating early treatment related changes to allow for an early intervention in case of incomplete disease response. One method for looking at treatment evaluation is to detect changes in MRI markers on a voxel-by-voxel basis following treatment. Changes in MRI markers may be correlated with disease recurrence and complete or partial response. In order to facilitate voxel-by-voxel imaging related treatment changes, and also to evaluate morphologic changes in the gland post treatment, the pre- and post-radiated MRI must first be brought into spatial alignment via image registration. However, EBRT induces changes in the prostate volume and distortion to the internal anatomy of the prostate following radiation treatment. The internal substructures of the prostate, the central gland (CG) and peripheral zone (PZ), may respond to radiation differently, and their resulting shapes may change drastically. Biomechanical models of the prostate that have been previously proposed tend to focus on how external forces affect the surface of the prostate (not the internals), and assume that the prostate is a volume-preserving entity. In this work we present DoCD, a biomechanical model for automatically registering pre-, post-EBRT MRI with the aim of expressly modeling the (1) changes in volume, and (2) changes to the CG and PZ. DoCD was applied to a cohort of 30 patients and achieved a root mean square error of 2.994 mm, which was statistically significantly better a traditional biomechanical model which did not consider changes to the internal anatomy of the prostate (mean of 5.071 mm).
Collapse
Affiliation(s)
- Robert Toth
- Rutgers, The State University of New Jersey. New Brunswick, NJ ; Case Western Reserve University, Cleveland, OH
| | | | | | - John Kurhanewicz
- Department of Radiology, University of California, San Francisco, CA
| | | |
Collapse
|
17
|
Reaungamornrat S, Liu WP, Wang AS, Otake Y, Nithiananthan S, Uneri A, Schafer S, Tryggestad E, Richmon J, Sorger JM, Siewerdsen JH, Taylor RH. Deformable image registration for cone-beam CT guided transoral robotic base-of-tongue surgery. Phys Med Biol 2013; 58:4951-79. [PMID: 23807549 PMCID: PMC3990286 DOI: 10.1088/0031-9155/58/14/4951] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Transoral robotic surgery (TORS) offers a minimally invasive approach to resection of base-of-tongue tumors. However, precise localization of the surgical target and adjacent critical structures can be challenged by the highly deformed intraoperative setup. We propose a deformable registration method using intraoperative cone-beam computed tomography (CBCT) to accurately align preoperative CT or MR images with the intraoperative scene. The registration method combines a Gaussian mixture (GM) model followed by a variation of the Demons algorithm. First, following segmentation of the volume of interest (i.e. volume of the tongue extending to the hyoid), a GM model is applied to surface point clouds for rigid initialization (GM rigid) followed by nonrigid deformation (GM nonrigid). Second, the registration is refined using the Demons algorithm applied to distance map transforms of the (GM-registered) preoperative image and intraoperative CBCT. Performance was evaluated in repeat cadaver studies (25 image pairs) in terms of target registration error (TRE), entropy correlation coefficient (ECC) and normalized pointwise mutual information (NPMI). Retraction of the tongue in the TORS operative setup induced gross deformation >30 mm. The mean TRE following the GM rigid, GM nonrigid and Demons steps was 4.6, 2.1 and 1.7 mm, respectively. The respective ECC was 0.57, 0.70 and 0.73, and NPMI was 0.46, 0.57 and 0.60. Registration accuracy was best across the superior aspect of the tongue and in proximity to the hyoid (by virtue of GM registration of surface points on these structures). The Demons step refined registration primarily in deeper portions of the tongue further from the surface and hyoid bone. Since the method does not use image intensities directly, it is suitable to multi-modality registration of preoperative CT or MR with intraoperative CBCT. Extending the 3D image registration to the fusion of image and planning data in stereo-endoscopic video is anticipated to support safer, high-precision base-of-tongue robotic surgery.
Collapse
Affiliation(s)
- S Reaungamornrat
- Department of Computer Science, Johns Hopkins University, Baltimore, MD 21218, USA
| | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
18
|
Mori S, Inaniwa T, Furukawa T, Zenklusen S, Shirai T, Noda K. Effects of a difference in respiratory cycle between treatment planning and irradiation for phase-controlled rescanning and carbon pencil beam scanning. Br J Radiol 2013; 86:20130163. [PMID: 23833034 DOI: 10.1259/bjr.20130163] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
OBJECTIVE To evaluate the impact of variation in respiratory cycle between treatment planning and irradiation for pencil beam scanning and phase-controlled rescanning (PCR) on the resulting dose distribution, we conducted a simulation study based on four-dimensional CT (4DCT) data for lung cancer patients. METHODS 4DCT data were acquired for seven patients with lung tumours. Treatment planning was designed to ensure the delivery of 95% of the prescribed dose to the clinical target volume in respective phases of the 4DCT by taking account of intrafractional beam range variations. Carbon ion pencil beam scanning dose distributions were calculated for various respiratory cycles that differed from the reference respiration (=4.4 s) but which stayed regular during irradiation. The number of rescannings was changed to 1, 4 or 8 times. PCR was correlated with the gating window in treatment planning to calculate the beam weighting map. RESULTS 8×PCR improved dose conformation to the target for all irradiation respiratory cycles. Minimum dose (Dmin) and lowest dose encompassing 95% of the target (D95) values with 4×PCR were decreased from 94.1% and 98.1% to 88.4% and 93.5% with an altered irradiation respiratory cycle of 2.4 s. However, these values were improved with 8×PCR to over 94.9% for Dmin and 98.6% for D95 for respective irradiation respiratory cycles. CONCLUSION Pencil beam scanning treatment with eight or more PCRs consistently improved dose conformation for moving lung targets even when different respiratory cycles were used for treatment planning and irradiation. ADVANCES IN KNOWLEDGE Scanning treatment with eight or more rescannings consistently improved dose homogeneity to a moving target even though respiratory cycles varied during treatment.
Collapse
Affiliation(s)
- S Mori
- Medical Physics Research Group, Research Center for Charged Particle Therapy, National Institute of Radiological Sciences, Chiba, Japan.
| | | | | | | | | | | |
Collapse
|
19
|
Spirka T, Kenton K, Brubaker L, Damaser MS. Effect of material properties on predicted vesical pressure during a cough in a simplified computational model of the bladder and urethra. Ann Biomed Eng 2013; 41:185-94. [PMID: 22907256 PMCID: PMC3677772 DOI: 10.1007/s10439-012-0637-x] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2012] [Accepted: 08/03/2012] [Indexed: 10/28/2022]
Abstract
Stress urinary incontinence is a condition that affects mainly women and is characterized by the involuntary loss of urine in conjunction with an increase in abdominal pressure but in the absence of a bladder contraction. In spite of the large number of women affected by this condition, little is known regarding the mechanics associated with the maintenance of continence in women. Urodynamic measurements of the pressure acting on the bladder and the pressures developed within the bladder and the urethra offer a potential starting point for constructing computational models of the bladder and urethra during stress events. The measured pressures can be utilized in these models to provide information to specify loads and validate the models. The main goals of this study were to investigate the feasibility of incorporating human urodynamic pressure data into a computational model of the bladder and the urethra during a cough and determine if the resulting model could be validated through comparison of predicted and measured vesical pressure. The results of this study indicated that simplified models can predict vesical pressures that differ by less than 5 cmH(2)O (<10%) compared to urodynamic pressure measurements. In addition, varying material properties had a minimal impact on the vesical pressure and displacements predicted by the model. The latter finding limits the use of vesical pressure as a validation criterion since different parameters can yield similar results in the same model. However, the insensitivity of vesical pressure predictions to material properties ensures that the outcome of our models is not highly sensitive to tissue material properties, which are not well characterized.
Collapse
Affiliation(s)
- Thomas Spirka
- Department of Biomedical Engineering, Cleveland Clinic Lerner Research Institute, Cleveland Clinic Main Campus, Mail Code ND20, 9500 Euclid Avenue, Cleveland, OH 44195, USA
| | | | | | | |
Collapse
|
20
|
Fedorov A, Tuncali K, Fennessy FM, Tokuda J, Hata N, Wells WM, Kikinis R, Tempany CM. Image registration for targeted MRI-guided transperineal prostate biopsy. J Magn Reson Imaging 2012; 36:987-92. [PMID: 22645031 PMCID: PMC3434292 DOI: 10.1002/jmri.23688] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2011] [Accepted: 03/28/2012] [Indexed: 11/09/2022] Open
Abstract
PURPOSE To develop and evaluate image registration methodology for automated re-identification of tumor-suspicious foci from preprocedural MR exams during MR-guided transperineal prostate core biopsy. MATERIALS AND METHODS A hierarchical approach for automated registration between planning and intra-procedural T2-weighted prostate MRI was developed and evaluated on the images acquired during 10 consecutive MR-guided biopsies. Registration accuracy was quantified at image-based landmarks and by evaluating spatial overlap for the manually segmented prostate and sub-structures. Registration reliability was evaluated by simulating initial mis-registration and analyzing the convergence behavior. Registration precision was characterized at the planned biopsy targets. RESULTS The total computation time was compatible with a clinical setting, being at most 2 min. Deformable registration led to a significant improvement in spatial overlap of the prostate and peripheral zone contours compared with both rigid and affine registration. Average in-slice landmark registration error was 1.3 ± 0.5 mm. Experiments simulating initial mis-registration resulted in an estimated average capture range of 6 mm and an average in-slice registration precision of ±0.3 mm. CONCLUSION Our registration approach requires minimum user interaction and is compatible with the time constraints of our interventional clinical workflow. The initial evaluation shows acceptable accuracy, reliability and consistency of the method.
Collapse
Affiliation(s)
- Andriy Fedorov
- Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, USA.
| | | | | | | | | | | | | | | |
Collapse
|
21
|
Zhong H, Kim J, Li H, Nurushev T, Movsas B, Chetty IJ. A finite element method to correct deformable image registration errors in low-contrast regions. Phys Med Biol 2012; 57:3499-515. [PMID: 22581269 DOI: 10.1088/0031-9155/57/11/3499] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Image-guided adaptive radiotherapy requires deformable image registration to map radiation dose back and forth between images. The purpose of this study is to develop a novel method to improve the accuracy of an intensity-based image registration algorithm in low-contrast regions. A computational framework has been developed in this study to improve the quality of the 'demons' registration. For each voxel in the registration's target image, the standard deviation of image intensity in a neighborhood of this voxel was calculated. A mask for high-contrast regions was generated based on their standard deviations. In the masked regions, a tetrahedral mesh was refined recursively so that a sufficient number of tetrahedral nodes in these regions can be selected as driving nodes. An elastic system driven by the displacements of the selected nodes was formulated using a finite element method (FEM) and implemented on the refined mesh. The displacements of these driving nodes were generated with the 'demons' algorithm. The solution of the system was derived using a conjugated gradient method, and interpolated to generate a displacement vector field for the registered images. The FEM correction method was compared with the 'demons' algorithm on the computed tomography (CT) images of lung and prostate patients. The performance of the FEM correction relating to the 'demons' registration was analyzed based on the physical property of their deformation maps, and quantitatively evaluated through a benchmark model developed specifically for this study. Compared to the benchmark model, the 'demons' registration has the maximum error of 1.2 cm, which can be corrected by the FEM to 0.4 cm, and the average error of the 'demons' registration is reduced from 0.17 to 0.11 cm. For the CT images of lung and prostate patients, the deformation maps generated by the 'demons' algorithm were found unrealistic at several places. In these places, the displacement differences between the 'demons' registrations and their FEM corrections were found in the range of 0.4 and 1.1 cm. The mesh refinement and FEM simulation were implemented in a single thread application which requires about 45 min of computation time on a 2.6 GHz computer. This study has demonstrated that the FEM can be integrated with intensity-based image registration algorithms to improve their registration accuracy, especially in low-contrast regions.
Collapse
Affiliation(s)
- Hualiang Zhong
- Department of Radiation Oncology, Henry Ford Health System, Detroit, MI, USA.
| | | | | | | | | | | |
Collapse
|
22
|
Liu H, Wu Q. Evaluations of an adaptive planning technique incorporating dose feedback in image-guided radiotherapy of prostate cancer. Med Phys 2012; 38:6362-70. [PMID: 22149819 DOI: 10.1118/1.3658567] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE Online image guidance (IG) has been used to effectively correct the setup error and inter-fraction rigid organ motion for prostate cancer. However, planning margins are still necessary to account for uncertainties such as deformation and intra-fraction motion. The purpose of this study is to investigate the effectiveness of an adaptive planning technique incorporating offline dose feedback to manage inter-fraction motion and residuals from online correction. METHODS Repeated helical CT scans from 28 patients were included in the study. The contours of prostate and organs-at-risk (OARs) were delineated on each CT, and online IG was simulated by matching center-of-mass of prostate between treatment CTs and planning CT. A seven beam intensity modulated radiation therapy (IMRT) plan was designed for each patient on planning CT for a total of 15 fractions. Dose distribution at each fraction was evaluated based on actual contours of the target and OARs from that fraction. Cumulative dose up to each fraction was calculated by tracking each voxel based on a deformable registration algorithm. The cumulative dose was compared with the dose from initial plan. If the deviation exceeded the pre-defined threshold, such as 2% of the D₉₉ to the prostate, an adaptive planning technique called dose compensation was invoked, in which the cumulative dose distribution was fed back to the treatment planning system and the dose deficit was made up through boost radiation in future treatment fractions. The dose compensation was achieved by IMRT inverse planning. Two weekly compensation delivery strategies were simulated: one intended to deliver the boost dose in all future fractions (schedule A) and the other in the following week only (schedule B). The D₉₉ to prostate and generalized equivalent uniform dose (gEUD) to rectal wall and bladder were computed and compared with those without the dose compensation. RESULTS If only 2% underdose is allowed at the end of the treatment course, then 11 patients fail. If the same criteria is assessed at the end of each week (every five fractions), then 14 patients fail, with three patients failing the 1st or 2nd week but passing at the end. The average dose deficit from these 14 patients was 4.4%. They improved to 2% after the weekly compensation. Out of these 14 patients who needed dose compensation, ten passed the dose criterion after weekly dose compensation, three patients failed marginally, and one patient still failed the criterion significantly (10% deficit), representing 3.6% of the patient population. A more aggressive compensation frequency (every three fractions) could successfully reduce the dose deficit to the acceptable level for this patient. The average number of required dose compensation re-planning per patient was 0.82 (0.79) per patient for schedule A (B) delivery strategy. The doses to OARs were not significantly different from the online IG only plans without dose compensation. CONCLUSIONS We have demonstrated the effectiveness of offline dose compensation technique in image-guided radiotherapy for prostate cancer. It can effectively account for residual uncertainties which cannot be corrected through online IG. Dose compensation allows further margin reduction and critical organs sparing.
Collapse
Affiliation(s)
- Han Liu
- Department of Radiation Oncology, Duke University Medical Center, Durham, North Carolina 27710, USA
| | | |
Collapse
|
23
|
McGrath DM, Foltz WD, Al-Mayah A, Niu CJ, Brock KK. Quasi-static magnetic resonance elastography at 7 T to measure the effect of pathology before and after fixation on tissue biomechanical properties. Magn Reson Med 2011; 68:152-65. [DOI: 10.1002/mrm.23223] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2011] [Revised: 08/15/2011] [Accepted: 08/29/2011] [Indexed: 01/22/2023]
|
24
|
Hu Y, Carter TJ, Ahmed HU, Emberton M, Allen C, Hawkes DJ, Barratt DC. Modelling prostate motion for data fusion during image-guided interventions. IEEE TRANSACTIONS ON MEDICAL IMAGING 2011; 30:1887-1900. [PMID: 21632296 DOI: 10.1109/tmi.2011.2158235] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
There is growing clinical demand for image registration techniques that allow multimodal data fusion for accurate targeting of needle biopsy and ablative prostate cancer treatments. However, during procedures where transrectal ultrasound (TRUS) guidance is used, substantial gland deformation can occur due to TRUS probe pressure. In this paper, the ability of a statistical shape/motion model, trained using finite element simulations, to predict and compensate for this source of motion is investigated. Three-dimensional ultrasound images acquired on five patient prostates, before and after TRUS-probe-induced deformation, were registered using a nonrigid, surface-based method, and the accuracy of different deformation models compared. Registration using a statistical motion model was found to outperform alternative elastic deformation methods in terms of accuracy and robustness, and required substantially fewer target surface points to achieve a successful registration. The mean final target registration error (based on anatomical landmarks) using this method was 1.8 mm. We conclude that a statistical model of prostate deformation provides an accurate, rapid and robust means of predicting prostate deformation from sparse surface data, and is therefore well-suited to a number of interventional applications where there is a need for deformation compensation.
Collapse
Affiliation(s)
- Yipeng Hu
- UCL Centre for Medical Image Computing, the Departmentof Medical Physics and Bioengineering, and the Department of ComputerScience, University College London, UK.
| | | | | | | | | | | | | |
Collapse
|
25
|
Chai X, van Herk M, van de Kamer JB, Hulshof MCCM, Remeijer P, Lotz HT, Bel A. Finite element based bladder modeling for image-guided radiotherapy of bladder cancer. Med Phys 2010; 38:142-50. [DOI: 10.1118/1.3523624] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
26
|
Vásquez Osorio EM, Hoogeman MS, Teguh DN, Al-Mamgani A, Kolkman-Deurloo IKK, Bondar L, Levendag PC, Heijmen BJM. Three-dimensional dose addition of external beam radiotherapy and brachytherapy for oropharyngeal patients using nonrigid registration. Int J Radiat Oncol Biol Phys 2010; 80:1268-77. [PMID: 21129854 DOI: 10.1016/j.ijrobp.2010.10.006] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2010] [Revised: 09/30/2010] [Accepted: 10/02/2010] [Indexed: 10/18/2022]
Abstract
PURPOSE To develop and evaluate a method for adding dose distributions of combined external beam radiotherapy (EBRT) and brachytherapy (BT) for oropharyngeal patients. METHODS AND MATERIALS Two computed tomography (CT) scans were used for 5 patients: the EBRT CT, used for EBRT planning, and the BT CT, acquired after catheter implantation. For each scan, the salivary glands and the chewing and swallowing muscles were contoured, and a dose distribution was calculated. A nonrigid transformation was obtained by registering the organs' surfaces. Then the BT dose distribution was mapped onto the EBRT dose distribution by applying the transformation obtained. To account for differences in fractionation, the physical doses were converted to equivalent dose in 2 Gy (EQD(2)), and the total dose was found by adding dose voxel by voxel. The robustness of the dose addition was investigated by varying delineations and input parameters of the registration method and by varying the α/β parameter for EQD(2). The effect of the perturbations was quantified using dose-volume histograms (DVH) and gamma analyses (distance-to-agreement/dose-difference = 1 mm/1 Gy). RESULTS The variations in input parameters and delineations caused only small perturbations in the DVH of the added dose distributions. For most organs the gamma index was low, and it was moderately elevated for organs lying in areas with a steep gradient (median gamma index ≤ 2.3 for constrictor muscles, ≤ 0.7 for all other organs). CONCLUSIONS The presented method allows adding dose distributions of combined EBRT and BT for oropharyngeal patients. In general, the method is reliable and robust with respect to uncertainties in organ delineation, perturbations in input parameters of the method, and α/β values.
Collapse
Affiliation(s)
- Eliana M Vásquez Osorio
- Department of Radiation Oncology, Daniel den Hoed Cancer Center, Erasmus Medical Center, Rotterdam, The Netherlands.
| | | | | | | | | | | | | | | |
Collapse
|
27
|
Al-Mayah A, Moseley J, Velec M, Hunter S, Brock K. Deformable image registration of heterogeneous human lung incorporating the bronchial tree. Med Phys 2010; 37:4560-71. [PMID: 20964173 DOI: 10.1118/1.3471020] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE To investigate the effect of the bronchial tree on the accuracy of biomechanical-based deformable image registration of human lungs. METHODS Three dimensional finite element models have been developed using four dimensional computed tomography image data of ten lung cancer patients. Each model is built of a body, left and right lungs, tumor, and bronchial trees. Triangular shell elements are used for the bronchial trees while tetrahedral elements are used for other components. Hyperelastic material properties based on experimental investigation on human lungs are used for the lung parenchyma. Different material properties are assigned for the bronchial tree using five values for the modulus of elasticity of 0.01, 0.12, 0.5, 10, and 18 MPa. Lungs are modeled to slide inside chest cavities by applying frictionless contact surfaces between each lung and corresponding chest cavity. The accuracy of the models is examined using an average of 40 bronchial bifurcation points identified on inhale and exhale images. Relative accuracy is evaluated by comparing the displacement of all nodes within the lungs as well as the dosimetric difference at the exhale position predicted by the model. RESULTS There is no significant effect of bronchial tree on the model accuracy based on the bifurcation points analysis. However, on the local level, using an average of 38 000 nodes, there is a maximum difference of 8.5 mm in the deformation of the bronchial trees, as the modulus of elasticity of the bronchial trees increases from 0.01 to 18 MPa; however, more than 96% of nodes are within a 2.5 mm difference in each direction. The average dose difference at the predicted exhale position is less than 35 cGy between the models. CONCLUSIONS The bronchial tree has little effect on the global deformation and the accuracy of deformable image registration of lungs. Hence, the homogenous model is a reasonable assumption. Since there are some local deformation differences between nodes as the material properties of the bronchial tree change that may affect the accuracy of dosimetric results, heterogeneity may be required for a smaller scale modeling of lungs.
Collapse
Affiliation(s)
- Adil Al-Mayah
- Radiation Medicine Program, Princess Margaret Hospital, 610 University Avenue, Toronto, Ontario M5G 2M9, Canada.
| | | | | | | | | |
Collapse
|
28
|
Hu Y, van den Boom R, Carter T, Taylor Z, Hawkes D, Ahmed HU, Emberton M, Allen C, Barratt D. A comparison of the accuracy of statistical models of prostate motion trained using data from biomechanical simulations. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2010; 103:262-72. [PMID: 20869389 DOI: 10.1016/j.pbiomolbio.2010.09.009] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/12/2010] [Revised: 08/31/2010] [Accepted: 09/15/2010] [Indexed: 11/18/2022]
Abstract
Statistical shape models (SSM) are widely used in medical image analysis to represent variability in organ shape. However, representing subject-specific soft-tissue motion using this technique is problematic for applications where imaging organ changes in an individual is not possible or impractical. One solution is to synthesise training data by using biomechanical modelling. However, for many clinical applications, generating a biomechanical model of the organ(s) of interest is a non-trivial task that requires a significant amount of user-interaction to segment an image and create a finite element mesh. In this study, we investigate the impact of reducing the effort required to generate SSMs and the accuracy with which such models can predict tissue displacements within the prostate gland due to transrectal ultrasound probe pressure. In this approach, the finite element mesh is based on a simplified geometric representation of the organs. For example, the pelvic bone is represented by planar surfaces, or the number of distinct tissue compartments is reduced. Such representations are much easier to generate from images than a geometrically accurate mesh. The difference in the median root-mean-square displacement error between different SSMs of prostate was <0.2 mm. We conclude that reducing the geometric complexity of the training model in this way made little difference to the absolute accuracy of SSMs to recover tissue displacements. The implication is that SSMs of organ motion based on simulated training data may be generated using simplified geometric representations, which are much more compatible with the time constraints of clinical workflows.
Collapse
Affiliation(s)
- Yipeng Hu
- Centre for Medical Image Computing, Department of Medical Physics and Bioengineering, University College London, London, UK.
| | | | | | | | | | | | | | | | | |
Collapse
|
29
|
Kuo HC, Chuang KS, Mah D, Wu A, Hong L, Yaparpalvi R, Kalnicki S. Multi-scale regularization approaches of non-parametric deformable registrations. J Digit Imaging 2010; 24:586-97. [PMID: 20574767 DOI: 10.1007/s10278-010-9313-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
Abstract
Most deformation algorithms use a single-value smoother during optimization. We investigate multi-scale regularizations (smoothers) during the multi-resolution iteration of two non-parametric deformable registrations (demons and diffeomorphic algorithms) and compare them to a conventional single-value smoother. Our results show that as smoothers increase, their convergence rate decreases; however, smaller smoothers also have a large negative value of the Jacobian determinant suggesting that the one-to-one mapping has been lost; i.e., image morphology is not preserved. A better one-to-one mapping of the multi-scale scheme has also been established by the residual vector field measures. In the demons method, the multi-scale smoother calculates faster than the large single-value smoother (Gaussian kernel width larger than 0.5) and is equivalent to the smallest single-value smoother (Gaussian kernel width equals to 0.5 in this study). For the diffeomorphic algorithm, since our multi-scale smoothers were implemented at the deformation field and the update field, calculation times are longer. For the deformed images in this study, the similarity measured by mean square error, normal correlation, and visual comparisons show that the multi-scale implementation has better results than large single-value smoothers, and better or equivalent for smallest single-value smoother. Between the two deformable registrations, diffeormophic method constructs better coherence space of the deformation field while the deformation is large between images.
Collapse
Affiliation(s)
- Hsiang-Chi Kuo
- Department of Radiation Oncology, Montefiore Medical Center, Bronx, NY 10461, USA.
| | | | | | | | | | | | | |
Collapse
|
30
|
Zhou J, Kim S, Jabbour S, Goyal S, Haffty B, Chen T, Levinson L, Metaxas D, Yue NJ. A 3D global-to-local deformable mesh model based registration and anatomy-constrained segmentation method for image guided prostate radiotherapy. Med Phys 2010; 37:1298-308. [DOI: 10.1118/1.3298374] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
31
|
Misra S, Macura KJ, Ramesh KT, Okamura AM. The importance of organ geometry and boundary constraints for planning of medical interventions. Med Eng Phys 2008; 31:195-206. [PMID: 18815068 DOI: 10.1016/j.medengphy.2008.08.002] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2008] [Revised: 07/01/2008] [Accepted: 08/04/2008] [Indexed: 10/21/2022]
Abstract
Realistic modeling of medical interventions involving tool-tissue interactions has been considered to be a key requirement in the development of high-fidelity simulators and planners. Organ geometry, soft-tissue constitutive laws, and boundary conditions imposed by the connective tissues surrounding the organ are some of the factors that govern the accuracy of medical intervention planning. In this study it is demonstrated that, for needle path planning, the organ geometry and boundary constraints surrounding the organ are the most important factors influencing the deformation. As an example, the procedure of needle insertion into the prostate (e.g. for biopsy or brachytherapy) is considered. Image segmentation is used to extract the anatomical details from magnetic resonance images, while object-oriented finite element analysis (OOF) software is used to generate finite element (FE) meshes from the segmented images. Two-dimensional FE simulations that account for complex anatomical details along with relative motion between the prostate and its surrounding structure using cohesive zone models are compared with traditional simulation models having simple organ geometry and boundary constraints. Nodal displacements for these simpler models were observed to be up to 14 times larger than those obtained from the anatomically accurate models.
Collapse
Affiliation(s)
- S Misra
- Department of Mechanical Engineering, The Johns Hopkins University, Baltimore, MD 21218, USA.
| | | | | | | |
Collapse
|
32
|
Lin L, Shi C, Liu Y, Swanson G, Papanikolaou N. Development of a novel post-processing treatment planning platform for 4D radiotherapy. Technol Cancer Res Treat 2008; 7:125-32. [PMID: 18345701 DOI: 10.1177/153303460800700205] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
The aim of this study is to develop an Automatic Post-processing Tool for four-dimensional (4D) treatment planning (APT4D) that enables the user to perform some necessary procedures related to 4D treatment planning, such as automated image registration, automatic propagation of regions of interest, and dose distribution transformation. Demons-based deformable registrations were performed to map the moving phase images (such as the end-inhalation phase or 0% phase) to the reference phase (typically the end-exhalation fixed phase or 50% phase). Contours were automatically propagated into the moving phase using the image registration results. The dose distribution of each moving phase was transformed to the fixed phase and subsequently was summed as an average with equal weighting factor. To validate the application of APT4D utility, the 4D computed tomography (CT) images of a lung cancer patient and an abdominal cancer patient were acquired and resorted into ten respiratory phases. 4D plans based on the 4D CT images were developed. The correlation coefficient ranged from 0.992 to 0.999 for the re-sampled deformed moving phase image against the fixed phase image for the lung patient plan and from 0.977 to 0.999 for the abdominal patient plan. For all the organs, the match indices between the manual contours and automatic contour propagation results were around 0.92 to 0.95. The 4D composite dose-volume histogram showed dosimetric reductions for liver and kidneys in the high dose region. The APT4D adds automation, efficiency, and functionality, while integrating the whole process of 4D treatment planning.
Collapse
Affiliation(s)
- Lan Lin
- Department of Medical Physics, Cancer Therapy and Research Center, San Antonio, TX 78229, USA.
| | | | | | | | | |
Collapse
|
33
|
Söhn M, Birkner M, Chi Y, Wang J, Di Y, Berger B, Alber M. Model-independent, multimodality deformable image registration by local matching of anatomical features and minimization of elastic energy. Med Phys 2008; 35:866-78. [PMID: 18404923 DOI: 10.1118/1.2836951] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
With respect to the demands of adaptive and 4D-radiotherapy applications, an algorithm is proposed for a fully automatic, multimodality deformable registration that follows the concept of translational relocation of regularly distributed image subvolumes governed by local anatomical features. Thereby, the problem of global deformable registration is broken down to multiple independent local registration steps which allows for straightforward parallelization of the algorithm. In a subsequent step, possible local misregistrations are corrected for by minimization of the elastic energy of the displacement field under consideration of image information. The final displacement field results from interpolation of the subvolume shift vectors. The algorithm can employ as a similarity measure both the correlation coefficient and mutual information. The latter allows the application to intermodality deformable registration problems. The typical calculation time on a modern multiprocessor PC is well below 1 min, which facilitates almost-interactive, "online" usage. CT-to-MRI and CT-to-cone-beam-CT registrations of head-and-neck data sets are presented, as well as inhale-to-exhale registrations of lung CT data sets. For quantitative evaluation of registration accuracy, a virtual thorax phantom was developed; additionally, a landmark-based evaluation on four lung respiratory-correlated CT data sets was performed. This consistently resulted in average registration residuals on the order of the voxel size or less (3D-residuals approximately 1-2 mm). Summarizing, the presented algorithm allows an accurate multimodality deformable registration with calculation times well below 1 min, and thus bears promise as a versatile basic tool in adaptive and 4D-radiotherapy applications.
Collapse
Affiliation(s)
- Matthias Söhn
- Section for Biomedical Physics, University Hospital for Radiation Oncology, Hoppe-Seyler-Strasse 3, 72076 Tiibingen, Germany.
| | | | | | | | | | | | | |
Collapse
|
34
|
Sharpe M, Brock KK. Quality Assurance of Serial 3D Image Registration, Fusion, and Segmentation. Int J Radiat Oncol Biol Phys 2008; 71:S33-7. [DOI: 10.1016/j.ijrobp.2007.06.087] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2007] [Revised: 06/19/2007] [Accepted: 06/20/2007] [Indexed: 11/28/2022]
|
35
|
Zhong H, Weiss E, Siebers JV. Assessment of dose reconstruction errors in image-guided radiation therapy. Phys Med Biol 2008; 53:719-36. [PMID: 18199911 PMCID: PMC2819061 DOI: 10.1088/0031-9155/53/3/013] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Dose reconstruction can be used to improve the accuracy of dose evaluation throughout a treatment course. Its working mechanism is based on deformable image registration (DIR). The purpose of this paper is to develop a method to estimate the dose reconstruction error associated with the inaccuracy of DIR algorithms. To reach this goal, we quantified dominant errors in DIR in terms of unbalanced energy (UE), which were compared with the standard displacement error (SDE). Their high similarity, characterized by Pearson correlation coefficient, was verified through nine 'demons' registration instances performed within simulated reference frames. Based on the similarity, the dose-warping discrepancy at each voxel was defined as a line integral of the dose gradient within the voxel's neighborhood whose boundary was determined by the voxel's UE value. From this definition, the dose reconstruction error was then calculated at each voxel on nine prostate computed tomography images, obtained from a patient treatment course. The average of the Pearson correlation coefficients between UE and SDE over the simulated registration instances was above 70%. The mean value of the dose reconstruction errors in a target volume was calculated for each of nine treatment fractions. The averaged percentage of these mean values with respect to the prescribed dose on the target volume was 1.68%. These results are consistent with contour-based mean dose error evaluations. This paper has established a relation between a registration error and its induced dose reconstruction discrepancy. It allows an automatic validation method to be developed to estimate the dose accumulation error at each voxel in clinical settings.
Collapse
Affiliation(s)
- Hualiang Zhong
- Virginia Commonwealth University, Richmond, VA 23298, USA.
| | | | | |
Collapse
|
36
|
Hu Y, Morgan D, Ahmed HU, Pendsé D, Sahu M, Allen C, Emberton M, Hawkes D, Barratt D. A Statistical Motion Model Based on Biomechanical Simulations for Data Fusion during Image-Guided Prostate Interventions. ACTA ACUST UNITED AC 2008; 11:737-44. [DOI: 10.1007/978-3-540-85988-8_88] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/05/2023]
|
37
|
|
38
|
Hensel JM, Ménard C, Chung PWM, Milosevic MF, Kirilova A, Moseley JL, Haider MA, Brock KK. Development of Multiorgan Finite Element-Based Prostate Deformation Model Enabling Registration of Endorectal Coil Magnetic Resonance Imaging for Radiotherapy Planning. Int J Radiat Oncol Biol Phys 2007; 68:1522-8. [PMID: 17674983 DOI: 10.1016/j.ijrobp.2007.04.004] [Citation(s) in RCA: 66] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2006] [Revised: 04/03/2007] [Accepted: 04/03/2007] [Indexed: 11/25/2022]
Abstract
PURPOSE Endorectal coil (ERC) magnetic resonance imaging (MRI) provides superior visualization of the prostate compared with computed tomography at the expense of deformation. This study aimed to develop a multiorgan finite element deformable method, Morfeus, to accurately co-register these images for radiotherapy planning. METHODS Patients with prostate cancer underwent fiducial marker implantation and computed tomography simulation for radiotherapy planning. A series of axial MRI scans were acquired with and without an ERC. The prostate, bladder, rectum, and pubic bones were manually segmented and assigned linear elastic material properties. Morfeus mapped the surface of the bladder and rectum between two imaged states, calculating the deformation of the prostate through biomechanical properties. The accuracy of deformation was measured as fiducial marker error and residual surface deformation between the inferred and actual prostate. The deformation map was inverted to deform from 100 cm(3) to no coil. RESULTS The data from 19 patients were analyzed. Significant prostate deformation occurred with the ERC (mean intrapatient range, 0.88 +/- 0.25 cm). The mean vector error in fiducial marker position (n = 57) was 0.22 +/- 0.09 cm, and the mean vector residual surface deformation (n = 19) was 0.15 +/- 0.06 cm for deformation from no coil to 100-cm(3) ERC, with an image vector resolution of 0.22 cm. Accurately deformed MRI scans improved soft-tissue resolution of the anatomy for radiotherapy planning. CONCLUSIONS This method of multiorgan deformable registration enabled accurate co-registration of ERC-MRI scans with computed tomography treatment planning images. Superior structural detail was visible on ERC-MRI, which has potential for improving target delineation.
Collapse
Affiliation(s)
- Jennifer M Hensel
- Radiation Medicine Program, Princess Margaret Hospital, Toronto, ON, Canada
| | | | | | | | | | | | | | | |
Collapse
|
39
|
Zhong H, Peters T, Siebers JV. FEM-based evaluation of deformable image registration for radiation therapy. Phys Med Biol 2007; 52:4721-38. [PMID: 17671331 DOI: 10.1088/0031-9155/52/16/001] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
This paper presents a new concept to automatically detect the neighborhood in an image where deformable registration is mis-performing. Specifically, the displacement vector field (DVF) from a deformable image registration is substituted into a finite-element-based elastic framework to calculate unbalanced energy in each element. The value of the derived energy indicates the quality of the DVF in its neighborhood. The new voxel-based evaluation approach is compared with three other validation criteria: landmark measurement, a finite element approach and visual comparison, for deformable registrations performed with the optical-flow-based 'demons' algorithm as well as thin-plate spline interpolation. This analysis was performed on three pairs of prostate CT images. The results of the analysis show that the four criteria give mutually comparable quantitative assessments on the six registration instances. As an objective concept, the unbalanced energy presents no requirement on boundary constraints in its calculation, different from traditional mechanical modeling. This method is automatic, and at voxel level suitable to evaluate deformable registration in a clinical setting.
Collapse
Affiliation(s)
- Hualiang Zhong
- Department of Radiation Oncology, Virginia Commonwealth University, Richmond, VA 23298, USA.
| | | | | |
Collapse
|
40
|
Sodee DB, Sodee AE, Bakale G. Synergistic Value of Single-Photon Emission Computed Tomography/Computed Tomography Fusion to Radioimmunoscintigraphic Imaging of Prostate Cancer. Semin Nucl Med 2007; 37:17-28. [PMID: 17161036 DOI: 10.1053/j.semnuclmed.2006.07.002] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
The rationale on which positron emission tomography/computed tomography (PET/CT) imaging is based, combining the functional features of PET with the anatomic detail of CT, provides many advantages that are easily transferable to single-photon emission computed tomography (SPECT)/CT imaging. Our efforts have focused on applying fused SPECT/CT imaging to identify prostate cancer and its metastasis and recurrence through radioimmunoscintigraphy (RIS). This application of RIS to imaging prostate cancer requires 2 key components: (1) a well-defined target associated with the cancer and (2) a "magic bullet" to seek that target. A well-characterized RIS target for prostate cancer is prostate-specific membrane antigen, or PSMA, and finding the bullet to seek this target with high sensitivity and specificity has been the focus of intensive study for nearly two decades. One of the candidate bullets developed is capromab pendetide, which is a monoclonal antibody that seeks PSMA. This antibody is commercially available as ProstaScint, which can be labeled with indium-111 to localize prostate cancer via SPECT imaging. In the course of applying fused SPECT/CT ProstaScint imaging to more than 800 prostate cancer cases, numerous refinements to our protocol have evolved that are aimed at staging the cancer with utmost accuracy. In addition to optimizing the localization of prostate cancer and its metastasis, these refinements also have been extended toward guiding both the implantation of radioactive seeds in brachytherapy and in other types of radiation therapy which is illustrated through 5 case reports. Progress in the therapeutic targeting of PSMA is also being actively explored, which has more universal ramifications because PSMA is found in the neovasculature of other types of cancers.
Collapse
Affiliation(s)
- D Bruce Sodee
- Department of Radiology, Division of Nuclear Medicine, University Hospitals of Cleveland and Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA.
| | | | | |
Collapse
|
41
|
Abstract
Ths paper examines several applications of deformable registration algorithms in the field of image-guided radiotherapy. The first part focuses on the description of input and output of deformable registration algorithms, with a brief review of conventional and most current methods. The typical applications of deformable registration are then reviewed on the basis of four practical examples. The first two sets of examples deal with the fusion of images obtained from the same patient (inter-fraction registration), with time intervals of several days between each image. The other two examples deal with the fusion of images obtained in immediate sequence (intra-fraction registration); in this case, the focus is the displacement during image acquisition or patient treatment (mainly due to respiratory movement), with time intervals in the order of magnitude of tenths of seconds. Finally, the registration of images of different patients (inter-patient registration) is also discussed. In conclusion, deformable registration has become a fundamental tool for image analysis in radiotherapy. Although extensive validation of the numerous existing methods is required before extending its clinical use, deformable registration is expected to become a standard methodology in the treatment planning systems in the near future.
Collapse
Affiliation(s)
- David Sarrut
- Radiotherapy Department, Centre Léon Bérard, 28 rue Laënnec, 69373 Lyon, France.
| |
Collapse
|