Segmentation, separation and pose estimation of prostate brachytherapy seeds in CT images.

International audienceIn this paper, we address the development of an automatic approach for the computation of pose information (position + orientation) of prostate brachytherapy loose seeds from 3D CT images. From an initial detection of a set of seed candidates in CT images using a threshold and connected component method, the orientation of each individual seed is estimated by using the principal components analysis (PCA) method. The main originality of this approach is the ability to classify the detected objects based on a priori intensity and volume information and to separate groups of closely spaced seeds using three competing clustering methods: the standard and a modified k-means method and a Gaussian mixture model with an Expectation-Maximization algorithm. Experiments were carried out on a series of CT images of two phantoms and patients. The fourteen patients correspond to a total of 1063 implanted seeds. Detections are compared to manual segmentation and to related work in terms of detection performance and calculation time. This automatic method has proved to be accurate and fast including the ability to separate groups of seeds in a reliable way and to determine the orientation of each seed. Such a method is mandatory to be able to compute precisely the real dose delivered to the patient post-operatively instead of assuming the alignment of seeds along the theoretical insertion direction of the brachytherapy needles

A 3D US Guidance System for Permanent Breast Seed Implantation: Development and Validation

Permanent breast seed implantation (PBSI) is a promising breast radiotherapy technique that suffers from operator dependence. We propose and have developed an intraoperative 3D ultrasound (US) guidance system for PBSI. A tracking arm mounted to a 3D US scanner registers a needle template to the image. Images were validated for linear and volumetric accuracy, and image quality in a volunteer. The tracking arm was calibrated, and the 3D image registered to the scanner. Tracked and imaged needle positions were compared to assess accuracy and a patient-specific phantom procedure guided with the system. Median/mean linear and volumetric error was ±1.1% and ±4.1%, respectively, with clinically suitable volunteer scans. Mean tracking arm error was 0.43mm and 3D US target registration error ≤0.87mm. Mean needle tip/trajectory error was 2.46mm/1.55°. Modelled mean phantom procedure seed displacement was 2.50mm. To our knowledge, this is the first reported PBSI phantom procedure with intraoperative 3D image guidance

Brachytherapy Seed and Applicator Localization via Iterative Forward Projection Matching Algorithm using Digital X-ray Projections

Interstitial and intracavitary brachytherapy plays an essential role in management of several malignancies. However, the achievable accuracy of brachytherapy treatment for prostate and cervical cancer is limited due to the lack of intraoperative planning and adaptive replanning. A major problem in implementing TRUS-based intraoperative planning is an inability of TRUS to accurately localize individual seed poses (positions and orientations) relative to the prostate volume during or after the implantation. For the locally advanced cervical cancer patient, manual drawing of the source positions on orthogonal films can not localize the full 3D intracavitary brachytherapy (ICB) applicator geometry. A new iterative forward projection matching (IFPM) algorithm can explicitly localize each individual seed/applicator by iteratively matching computed projections of the post-implant patient with the measured projections. This thesis describes adaptation and implementation of a novel IFPM algorithm that addresses hitherto unsolved problems in localization of brachytherapy seeds and applicators. The prototype implementation of 3-parameter point-seed IFPM algorithm was experimentally validated using a set of a few cone-beam CT (CBCT) projections of both the phantom and post-implant patient’s datasets. Geometric uncertainty due to gantry angle inaccuracy was incorporated. After this, IFPM algorithm was extended to 5-parameter elongated line-seed model which automatically reconstructs individual seed orientation as well as position. The accuracy of this algorithm was tested using both the synthetic-measured projections of clinically-realistic Model-6711 125I seed arrangements and measured projections of an in-house precision-machined prostate implant phantom that allows the orientations and locations of up to 100 seeds to be set to known values. The seed reconstruction error for simulation was less than 0.6 mm/3o. For the physical phantom experiments, IFPM absolute accuracy for position, polar angle, and azimuthal angel were (0.78 ± 0.57) mm, (5.8 ± 4.8)o, and (6.8 ± 4.0)o, respectively. It avoids the need to match corresponding seeds in each projection and accommodates incomplete data, overlapping seed clusters, and highly-migrated seeds. IFPM was further generalized from 5-parameter to 6-parameter model which was needed to reconstruct 3D pose of arbitrary-shape applicators. The voxelized 3D model of the applicator was obtained from external complex combinatorial geometric modeling. It is then integrated into the forward projection matching method for computing the 2D projections of the 3D ICB applicators, iteratively. The applicator reconstruction error for simulation was about 0.5 mm/2o. The residual 2D registration error (positional difference) between computed and actual measured applicator images was less than 1 mm for the intrauterine tandem and about 1.5 mm for the bilateral colpostats in each detector plane. By localizing the applicator’s internal structure and the sources, the effect of intra and inter-applicator attenuation can be included in the resultant dose distribution and CBCT metal streaking artifact mitigation. The localization accuracy of better than 1 mm and 6o has the potential to support more accurate Monte Carlo-based or 2D TG-43 dose calculations in clinical practice. It is hoped the clinical implementation of IFPM approach to localize elongated line-seed/applicator for intraoperative brachytherapy planning may have a positive impact on the treatment of prostate and cervical cancers

Deformable registration of X-ray and MRI for post-implant dosimetry in low-dose-rate prostate brachytherapy

Purpose Dosimetric assessment following permanent prostate brachytherapy (PPB) commonly involves seed localization using CT and prostate delineation using coregistered MRI. However, pelvic CT leads to additional imaging dose and requires significant resources to acquire and process both CT and MRI. In this study, we propose an automatic postimplant dosimetry approach that retains MRI for soft‐tissue contouring, but eliminates the need for CT and reduces imaging dose while overcoming the inconsistent appearance of seeds on MRI with three projection x rays acquired using a mobile C‐arm. Methods Implanted seeds are reconstructed using x rays by solving a combinatorial optimization problem and deformably registered to MRI. Candidate seeds are located in MR images using local hypointensity identification. X ray‐based seeds are registered to these candidate seeds in three steps: (a) rigid registration using a stochastic evolutionary optimizer, (b) affine registration using an iterative closest point optimizer, and (c) deformable registration using a local feature point search and nonrigid coherent point drift. The algorithm was evaluated using 20 PPB patients with x rays acquired immediately postimplant and T2‐weighted MR images acquired the next day at 1.5 T with mean 0.8 × 0.8 × 3.0 mmurn:x-wiley:00942405:media:mp13667:mp13667-math-0001 voxel dimensions. Target registration error (TRE) was computed based on the distance from algorithm results to manually identified seed locations using coregistered CT acquired the same day as the MRI. Dosimetric accuracy was determined by comparing prostate D90 determined using the algorithm and the ground truth CT‐based seed locations. Results The mean ± standard deviation TREs across 20 patients including 1774 seeds were 2.23 ± 0.52 mm (rigid), 1.99 ± 0.49 mm (rigid + affine), and 1.76 ± 0.43 mm (rigid + affine + deformable). The corresponding mean ± standard deviation D90 errors were 5.8 ± 4.8%, 3.4 ± 3.4%, and 2.3 ± 1.9%, respectively. The mean computation time of the registration algorithm was 6.1 s. Conclusion The registration algorithm accuracy and computation time are sufficient for clinical PPB postimplant dosimetry

Encapsulated Contrast Agent Markers for MRI-based Post-implant Dosimetry

Low-dose-rate prostate brachytherapy involves the implantation of tiny radioactive seeds into the prostate to treat prostate cancer. The current standard post-implant imaging modality is computed tomography (CT). On CT images, the radioactive seeds can be distinctively localized but delineation of the prostate and surrounding soft tissue is poor. Magnetic resonance imaging (MRI) provides better prostate and soft tissue delineation, but seed localization is difficult. To aid with seed localization, MRI markers with encapsulated contrast agent that provide positive-contrast on MRI images (Sirius MRI markers; C4 Imaging, Houston, TX) have been proposed to be placed adjacent to the negative-contrast seeds. This dissertation describes the development of the Sirius MRI markers for prostate post-implant dosimetry. First, I compared the dose-volume histogram and other dosimetry parameters generated by MIM Symphony (a brachytherapy treatment planning system that allow the use of MRI images for treatment planning; MIM Software Inc., Cleveland, OH) and VariSeed (a widely used brachytherapy treatment planning system; Varian Medical Systems, Inc., Palo Alto, CA), and found the dosimetry between both brachytherapy treatment planning systems to be comparable. To gain more insight into the MRI contrast characteristics of the Sirius MRI markers, I measured the Sirius MRI marker contrast agent\u27s spin-lattice and spin-spin relaxivities, and studied the relaxation characteristics\u27 dependence on MRI field strength, temperature, and orientation. From the Sirius MRI marker\u27s contrast agent relaxation characteristics, I systematically studied the effect of varying MRI scan parameters such as flip angle, number of excitations, bandwidth, field of view, slice thickness, and encoding steps, on the Sirius MRI markers\u27 signal and contrast, as well as image noise, artifact and scan time. On patients implanted with Sirius MRI markers, I evaluated the visibility of the Sirius MRI markers and image artifacts. Lastly, I semi-automated the localization of markers and seeds to more enable the efficient incorporation of Sirius MRI markers as part of the clinical post-implant workflow. Ultimately, the Sirius MRI markers may change the paradigm from CT-based to MRI-based post-implant dosimetry, for a more accurate understanding of dose-response relationships in patients undergoing low dose rate prostate brachytherapy

Automatic C-arm pose estimation via 2D/3D hybrid registration of a radiographic fiducial

ABSTRACT Motivation: In prostate brachytherapy, real-time dosimetry would be ideal to allow for rapid evaluation of the implant quality intra-operatively. However, such a mechanism requires an imaging system that is both real-time and which provides, via multiple C-arm fluoroscopy images, clear information describing the three-dimensional position of the seeds deposited within the prostate. Thus, accurate tracking of the C-arm poses proves to be of critical importance to the process. Methodology: We compute the pose of the C-arm relative to a stationary radiographic fiducial of known geometry by employing a hybrid registration framework. Firstly, by means of an ellipse segmentation algorithm and a 2D/3D feature based registration, we exploit known FTRAC geometry to recover an initial estimate of the C-arm pose. Using this estimate, we then initialize the intensity-based registration which serves to recover a refined and accurate estimation of the C-arm pose. Results: Ground-truth pose was established for each C-arm image through a published and clinically tested segmentation-based method. Using 169 clinical C-arm images and a ±10° and ±10 mm random perturbation of the ground-truth pose, the average rotation and translation errors were 0.68° (std = 0.06°) and 0.64 mm (std = 0.24 mm). Conclusion: Fully automated C-arm pose estimation using a 2D/3D hybrid registration scheme was found to be clinically robust based on human patient data

SEED LOCALIZATION IN IMAGE-GUIDED PROSTATE BRACHYTHERAPY INTRAOPERATIVE DOSIMETRY SYSTEMS

Prostate cancer is the most common cancer among men in the United States. Many treatments are available, but prostate brachytherapy is acknowledged as a standard treatment for patients with localized cancer. Prostate brachytherapy is a minimally invasive surgery involving the permanent implantation of approximately 100 grain-sized radioactive seeds into the prostate. While effective, contemporary practice of brachytherapy is suboptimal because it spreads the stages of planning, implant, and dosimetry over several weeks. Although brachytherapy is now moving towards intraoperative treatment planning (ITP) which integrates all three stages into a single day in the operating room,the American Brachytherapy Society states, “the major current limitation of ITP is the inability to localize the seeds in relation to the prostate.” While the procedure is traditionally guided by transrectal ultrasound (TRUS), poor image quality prevents TRUS from accurately localizing seeds to compute dosimetry intraoperatively. Alternative methods exist, but are generally impractical to implement in clinics worldwide. The subject of this dissertation is the development of two intraoperative dosimetry systems to practically solve the problem of seed localization in ITP. The first system fuses TRUS with X-ray fluoroscopy using the ubiquitous non-isocentric mobile C-arm.The primary contributions of this dissertation include an automatic fiducial and seed segmentation algorithm for fluoroscopic images, as well as a next generation intraoperative dosimetry system based on a fiducial with seed-like markers. Results from over 30 patients prove that both contributions are significant for localizing seeds with high accuracy and demonstrate the capability of detecting cold spots. The second intraoperative dosimetry system is based on photoacoustic imaging, and uses the already implemented TRUS probe to detect ultrasonic waves converted from electromagnetic waves generated by a laser. The primary contributions of this dissertation therefore also include a prototype benchtop photoacoustic system and an improved clinical version usable in the operating room. Results from gelatin phantoms, an ex vivo dog prostate, and an in vivo dog study reveal that multiple seeds are clearly visible with high contrast using photoacoustic imaging at clinically safe laser energies.Together, both systems significantly progress the latest technologies to provide optimal care to patients through ITP

Prostate volume delineation and seed localization using a 3-T magnetic resonance spectrometer

With approximately one in six men affected by prostate cancer at some point in their lives, effective treatment of the disease remains a focus of oncology research. Effective treatment using radiation requires the delivery of a significant dose to the prostate volume while sparing surrounding sensitive structures. Treatment success can then be determined by localization of the seeds following implantation and the calculation of a dose distribution across the target volume. Magnetic Resonance Imaging (MRI) yields images with soft tissue contrast that is superior to CT or ultrasound, but has been under-appreciated as a dosimetric tool due to the difficulty in localizing the implanted seeds; To optimize scan parameters for seed localization, a phantom was constructed of tissue-equivalent gelatin. Seeds were implanted during construction so various scan protocols could be tested for seed visualization and volume calculation prior to patient studies. Five healthy volunteers and five patients with permanently implanted seeds were then imaged to validate the phantom studies. Images were evaluated based on anatomical clarity and seed visualization rates; Optimization of the scan protocols for use with this equipment yields images with clearly defined anatomical boundaries as well as clearly defined seeds. Phantom volume measurements deviated from known values by less than 2.5% T2-weighted images provide superior anatomical delineation, but suffer from broad susceptibility artifacts that make determination of seed locations difficult. Proton density-weighted images clearly show seed locations and tissue margins. The selection of a 1 mm slice thickness and a 4 mm interstice gap allowed maximum seed visualization rates of 93.3%; Keywords: prostate, brachytherapy, dosimetry, magnetic resonance imaging

Prostate deformation from inflatable rectal probe cover and dosimetric effects in prostate seed implant brachytherapy: Deformation on prostate permanent brachytherapy dosimetry

Prostate brachytherapy is an important treatment technique for patients with localized prostate cancer. An inflatable rectal ultrasound probe cover is frequently utilized during the procedure to adjust for unfavorable prostate position relative to the implant grid. However, the inflated cover causes prostate deformation, which is not accounted for during dosimetric planning. Most of the therapeutic dose is delivered after the procedure when the prostate and surrounding organs-at-risk are less deformed. The aim of this study is to quantify the potential dosimetry changes between the initial plan (prostate deformed) and the more realistic dosimetry when the prostate is less deformed without the cover