Seed misplacement and stabilizing needles in transperineal permanent prostate implants

A review on tissue-needle interaction and path planning models for bevel tip type flexible needle minimal intervention

A flexible needle has emerged as a crucial clinical technique in contemporary medical practices, particularly for minimally invasive interventions. Its applicability spans diverse surgical domains such as brachytherapy, cardiovascular surgery, neurosurgery and others. Notably, flexible needles find utility in biopsies requiring deep skin penetration to access infected areas. Despite its minimally invasive advantages, the precise guidance of the needle to its intended target, while avoiding damage to bones, blood vessels, organs and tissues, remains a significant challenge for researchers. Consequently, extensive research has been dedicated to enhancing the steering and accuracy of flexible needles. Here, we aim to elucidate the recent advancements, trends and perspectives in flexible needle steering models and path planning over the last 15 years. The discussed models encompass various types, including symmetric-tip needles, curved-tip needles, tendon-actuated needles, programmable needles and the innovative fracture-directed waterjet needles. Moreover, the paper offers a comprehensive analysis, comparing the trajectories followed by these needle models to attain the desired target with minimal tissue damage. By delving into these aspects, the paper contributes to a deeper understanding of the current landscape of flexible needle technology and guides future research directions in this dynamic field.

Study of the surgical needle and biological soft tissue interaction phenomenon during insertion process for medical application: A Survey

The insertion of the surgical needle in soft tissue has involved significant interest in the current time because of its purpose in minimally invasive surgery (MIS) and percutaneous events like biopsies, PCNL, and brachytherapy. This study represents a review of the existing condition of investigation on insertion of a surgical needle in biological living soft tissue material. As observes the issue from numerous phases, like, analysis of the cutting forces modeling (insertion), tissue material deformation, analysis of the needle deflection for the period of the needle insertion, and the robot-controlled insertion procedures. All analysis confirms that the total needle insertion force is the total of dissimilar forces spread sideways the shaft of the insertion needle for example cutting force, stiffness force, and frictional force. Various investigations have analyzed all these kinds of forces during the needle insertion process. The force data in several measures are applied for recognizing the biological tissue materials as the needle is penetrated or for path planning. The deflection of the needle during insertion and tissue material deformation is the main trouble for defined needle placing and efforts have been prepared to model them. Applying existing models numerous insertion methods are established that are discussed in this review.

Feasibility of extracting tissue material properties via cohesive elements: a finite element approach to probe insertion procedures in non-invasive spine surgeries

Modeling the mechanical behavior of soft tissue probe insertion remains a challenging endeavor due to involved interdependent phenomena comprising tissue nonlinear deformation, contact between the probe and the tissue, crack propagation, and viscoelastic effects. To that matter, cohesive elements allow simulating crack formation and propagation, which provides a promising path to modeling the mechanical behavior of probe insertion in soft tissues. As such, the aim of the present study was to investigate the feasibility of devising and integrating an algorithm in a finite element (FE) case study in efforts of reverse engineering the material properties of non-homogeneous soft tissues. A layered nonlinear tissue model with a cohesive zone was created in the commercial software ABAQUS. Material properties were iteratively modified via a hybrid gradient descent optimization algorithm: minimizing the resultant error to first find optimum Ogden's hyperelastic parameters, followed by obtaining the damage parameters. Perceived material properties were then compared to those obtained via experimental human cadaver testing. Under the investigated four-layered muscle model, numerical results overlapped, to a great extent, with six different force-insertion experimental profiles with an average error of [Formula: see text] 15%. The best profile fit was realized when the highest sudden force drop was less than 60% of the peak force. Lastly, the FE analysis revealed an increase in stiffness as the probe advanced inside the tissue. The optimization algorithm demonstrated its capability to reverse engineer the material parameters required for the FE analysis of real, non-homogeneous, soft tissues. The significance of this procedure lies within its ability to extract tissue material parameters, in real time, with little to no intervention or invasive experimental tests. This could potentially further serve as a database for different muscle layers and force-insertion profiles, used for surgeon and physician clinical training purposes.

Intraoperative optimization of seed implantation plan in breast brachytherapy

PURPOSE

Low-dose-rate permanent-seed (LDR-PS) brachytherapy has shown a great potential for treating breast cancer. An implantation scheme indicating the template pose and needle trajectories is determined before the operation. However, when performing the pre-planned scheme intraoperatively, a change of the patient's posture will cause seed placements away from the desired locations. Hence, the implantation scheme should update based on the current patient's posture.

METHODS

A numerical method of optimizing the implantation scheme for the LDR-PS breast brachytherapy is presented here. The proposed algorithm determines the fewest needle trajectories and template poses for delivering the seeds to the intraoperative desired locations. The clinical demand, such as the minimum distance between the chest wall and the needle, is considered in the optimization process.

RESULTS

The method was simulated for a given LDR-PS brachytherapy procedure to evaluate the optimal scheme as the number of the template poses changing. The optimization parameters of the needles' number and the implantation errors are used to adjust the algorithm outcome. The results show that the implantation schemes obtained by our method have a satisfactory accuracy in the cases of 2 or 3 template poses. The computation time is about 76s to 150s according to the number of the template poses from 1 to 3.

CONCLUSION

The proposed method can find the optimal implantation scheme corresponding to the current desired seed locations immediately once there is a change of patient's posture. This work can be applied to the robot-assisted LDR-PS breast brachytherapy for improving the operation accuracy and efficiency.

Gels, jets, mosquitoes, and magnets: a review of implantation strategies for soft neural probes

Implantable neuroelectronic interfaces have enabled breakthrough advances in the clinical diagnosis and treatment of neurological disorders, as well as in fundamental studies of brain function, behavior, and disease. Intracranial electroencephalography (EEG) mapping with stereo-EEG (sEEG) depth electrodes is routinely adopted for precise epilepsy diagnostics and surgical treatment, while deep brain stimulation has become the standard of care for managing movement disorders. Intracortical microelectrode arrays for high-fidelity recordings of neural spiking activity have led to impressive demonstrations of the power of brain-machine interfaces for motor and sensory functional recovery. Yet, despite the rapid pace of technology development, the issue of establishing a safe, long-term, stable, and functional interface between neuroelectronic devices and the host brain tissue still remains largely unresolved. A body of work spanning at least the last 15 years suggests that safe, chronic integration between invasive electrodes and the brain requires a close match between the mechanical properties of man-made components and the neural tissue. In other words, the next generation of invasive electrodes should be soft and compliant, without sacrificing biological and chemical stability. Soft neuroelectronic interfaces, however, pose a new and significant surgical challenge: bending and buckling during implantation that can preclude accurate and safe device placement. In this topical review, we describe the next generation of soft electrodes and the surgical implantation methods for safe and precise insertion into brain structures. We provide an overview of the most recent innovations in the field of insertion strategies for flexible neural electrodes such as dissolvable or biodegradable carriers, microactuators, biologically-inspired support structures, and electromagnetic drives. In our analysis, we also highlight approaches developed in different fields, such as robotic surgery, which could be potentially adapted and translated to the insertion of flexible neural probes.

Characterization and registration of 3D ultrasound for use in permanent breast seed implant brachytherapy treatment planning

PURPOSE

Permanent breast seed implant (PBSI) brachytherapy is a novel technique for early-stage breast cancer. Computed tomography (CT) images are used for treatment planning and freehand 2D ultrasound for implant guidance. The multimodality imaging approach leads to discrepancies in target identification. To address this, a prototype 3D ultrasound (3DUS) system was recently developed for PBSI. In this study, we characterize the 3DUS system performance, establish QA baselines, and develop and test a method to register 3DUS images to CT images for PBSI planning.

METHODS AND MATERIALS

3DUS system performance was characterized by testing distance and volume measurement accuracy, and needle template alignment accuracy. 3DUS-CT registration was achieved through point-based registration using a 3D-printed model designed and constructed to provide visible landmarks on both images and tested on an in-house made gel breast phantom.

RESULTS

The 3DUS system mean distance measurement accuracy was within 1% in axial, lateral, and elevational directions. A volumetric error of 3% was observed. The mean needle template alignment error was 1.0° ± 0.3 ° and 1.3 ± 0.5 mm. The mean 3DUS-CT registration error was within 3 mm when imaging at the breast centre or across all breast quadrants.

CONCLUSIONS

This study provided baseline data to characterize the performance of a prototype 3DUS system for PBSI planning and developed and tested a method to obtain accurate 3DUS-CT image registration for PBSI planning. Future work will focus on system validation and characterization in a clinical context as well as the assessment of impact on treatment plans.

Establishing a simulation-based education program for radiation oncology learners in permanent seed implant brachytherapy: Building validation evidence

PURPOSE

The purpose of this study was to establish a simulation-based education program for radiation oncology learners in permanent seed implant brachytherapy. The first step in formalizing any education program is a validation process that builds evidence-based verification that the learning environment is appropriate.

METHODS AND MATERIALS

The primary education task allowed practitioners to use an anthropomorphic breast phantom to simulate a permanent seed implant brachytherapy delivery. Validation evidence is built by generating data to assess learner and expert cohorts according to their proficiency. Each practitioner's performance during the simulation was evaluated by seed placement accuracy, procedural time-to-complete, and two qualitative evaluation tools-a global rating scale and procedural checklist.

RESULTS

The average seed placement accuracy (±SD) was 8.1 ± 3.5 mm compared to 6.1 ± 2.6 mm for the learner and expert cohort, respectively. The median (range) procedural time-to-complete was 64 (60-77) minutes and 43 (41-50) minutes for the learner and expert cohort, respectively. Seed placement accuracy (student t-test, p < 0.05) and procedural time-to-complete (Mann-Whitney U-test, p < 0.05) were statistically different between the cohorts. In both the global rating scale and procedural checklist, the expert cohort demonstrated improved proficiency compared to the learner cohort.

CONCLUSIONS

This validation evidence supports the utilization of this simulation environment toward appropriately capturing the delivery experience of practitioners. The results demonstrate that, in all areas of evaluation, expert cohort proficiency was superior to learner cohort proficiency. This methodology will be used to establish a simulation-based education program for radiation oncology learners in permanent seed implant brachytherapy.

Challenges in MR-only seed localization for postimplant dosimetry in permanent prostate brachytherapy

PURPOSE

An MR-only postimplant dosimetry workflow for low dose rate (LDR) brachytherapy could reduce patient burden, improve accuracy, and improve cost efficiency. However, localization of brachytherapy seeds on MRI scans remains a major challenge for this type of workflow. In this study, we propose and validate an MR-only seed localization method and identify remaining challenges.

METHODS AND MATERIALS

The localization method was based on template matching of simulations of complex-valued imaging artifacts around metal brachytherapy seeds. The method was applied to MRI scans of 25 prostate cancer patients who underwent LDR brachytherapy and for whom postimplant dosimetry was performed after 4 weeks. The seed locations found with the MR-only method were validated against the seed locations found on CT. The circumstances in which detection errors were made were classified to gain an insight in the nature of the errors.

RESULTS

A total of 1490 of 1557 (96%) seeds were correctly detected, while 67 false-positive errors were made. The correctly detected seed locations had a high spatial accuracy with an average error of 0.8 mm compared with CT. A majority of the false positives occurred near other seeds. Most false negatives were found in either stranded configurations without spacers or near other seeds.

CONCLUSIONS

The low detection error rate and high localization accuracy obtained by the complex-valued template matching approach are promising for future clinical application of MR-only dosimetry. The most important remaining challenge is robustness with regard to configurations of multiple seeds in close vicinity, such as in strands of seeds without spacers. This issue could potentially be resolved by simulating specific configurations of multiple seeds or by constraining the treatment planning to avoid these configurations, which could make the proposed method competitive with CT-based seed localization.

Haptic Feedback in Needle Insertion Modeling and Simulation

Needle insertion is the most basic skill in medical care, and training has to be imparted not only for physicians but also for nurses and paramedics. In most needle insertion procedures, haptic feedback from the needle is the main stimulus in which novices need training. For better patient safety, the classical methods of training the haptic skills have to be replaced with simulators based on new robotic and graphics technologies. This paper reviews the current advances in needle insertion modeling, classified into three sections: needle insertion models, tissue deformation models, and needle-tissue interaction models. Although understated in the literature, the classical and dynamic friction models, which are critical for needle insertion modeling, are also discussed. The experimental setup or the needle simulators that have been developed to validate the models are described. The need of psychophysics for needle simulators and psychophysical parameter analysis of human perception in needle insertion are discussed, which are completely ignored in the literature.

Quantifying 125I placement accuracy in prostate brachytherapy using postimplant transrectal ultrasound images

PURPOSE

The quality of a prostate brachytherapy implant depends on the accurate placement of sources. This study quantifies the misplacement of ¹²⁵I sources from the intended location using intraoperative ultrasound images.

METHODS AND MATERIALS

¹²⁵I sources were manually identified in the postimplant ultrasound images and compared to the preoperative plan. Due to the subjective nature of the identifying sources, only sources identified with high confidence were included in the analysis. Misplacements from the original intended coordinate were measured along the X, Y, and Z axes and were stratified between overall misplacements and regions of the prostate gland.

RESULTS

A total of 1619 ¹²⁵I sources using 357 strands were implanted in 15 patients' prostate glands, with 1197 (74%) confidently identified for misplacement analysis. The overall mean displacement was 0.49 cm and in the X, Y, and Z direction was 0.13, 0.15, and 0.38 cm, respectively. Greater source misplacement occurred in the anterior part of the prostate gland than the posterior part of the prostate gland by a factor 1.33 (p < 0.0001). Comparing sources in the lateral vs. medial regions of the prostate, no statistically significant differences on source misplacement were observed. Comparing misplacement in the base vs. midgland vs. apex identified the greatest difference between the base and midgland by a factor of 1.29 (p < 0.0001).

CONCLUSIONS

This study has identified significant misplacement of ¹²⁵I sources from their intended locations with the greatest error misplacement occurring in the Z direction. Source misplacement tends to occur more commonly in the anterior gland and in the base of the prostate.

Accuracy evaluation of a 3D-printed individual template for needle guidance in head and neck brachytherapy

To transfer the preplan for the head and neck brachytherapy to the clinical implantation procedure, a preplan-based 3D-printed individual template for needle insertion guidance had previously been designed and used. The accuracy of needle insertion using this kind template was assessed in vivo In the study, 25 patients with head and neck tumors were implanted with ¹²⁵I radioactive seeds under the guidance of the 3D-printed individual template. Patients were divided into four groups based on the site of needle insertion: the parotid and masseter region group (nine patients); the maxillary and paranasal region group (eight patients); the submandibular and upper neck area group (five patients); and the retromandibular region group (six patients). The distance and angular deviations between the preplanned and placed needles were compared, and the complications and time required for needle insertion were assessed. The mean entrance point distance deviation for all 619 needles was 1.18 ± 0.81 mm, varying from 0.857 ± 0.545 to 1.930 ± 0.843 mm at different sites. The mean angular deviation was 2.08 ± 1.07 degrees, varying from 1.85 ± 0.93 to 2.73 ± 1.18 degrees at different sites. All needles were manually inserted to their preplanned positions in a single attempt, and the mean time to insert one needle was 7.5 s. No anatomical complications related to inaccurately placed implants were observed. Using the 3D-printed individual template for the implantation of ¹²⁵I radioactive seeds in the head and neck region can accurately transfer a CT-based preplan to the brachytherapy needle insertion procedure. Moreover, the addition of individual template guidance can reduce the time required for implantation and minimize the damage to normal tissues.

Needle Tracking and Deflection Prediction for Robot-Assisted Needle Insertion Using 2D Ultrasound Images

In many types of percutaneous needle insertion surgeries, tissue deformation and needle deflection can create significant difficulties for accurate needle placement. In this paper, we present a method for automatic needle tracking in 2D ultrasound (US) images, which is used in a needle–tissue interaction model to estimate current and future needle tip deflection. This is demonstrated using a semi-automatic needle steering system. The US probe can be controlled to follow the needle tip or it can be stopped at an appropriate position to avoid tissue deformation of the target area. US images are used to fully parameterize the needle-tissue model. Once the needle deflection reaches a pre-determined threshold, the robot rotates the needle to correct the tip’s trajectory. Experimental results show that the final needle tip deflection can be estimated with average accuracies between 0.7[Formula: see text]mm and 1.0[Formula: see text]mm for insertions with and without rotation. The proposed method provides surgeons with improved US feedback of the needle tip deflection and minimizes the motion of the US probe to reduce tissue deformation of the target area.

Ultrasound-Guided Model Predictive Control of Needle Steering in Biological Tissue

In needle-based medical procedures, beveled tip flexible needles are steered inside soft tissue to reach the desired target locations. In this paper, we have developed an autonomous image-guided needle steering system that enhances targeting accuracy in needle insertion while minimizing tissue trauma. The system has three main components. First is a novel mechanics-based needle steering model that predicts needle deflection and accepts needle tip rotation as an input for needle steering. The second is a needle tip tracking system that determines needle deflection from the ultrasound images. The needle steering model employs the estimated needle deflection at the present time to predict needle tip trajectory in the future steps. The third component is a nonlinear model predictive controller (NMPC) that steers the needle inside the tissue by rotating the needle beveled tip. The MPC controller calculates control decisions based on iterative optimization of the predictions of the needle steering model. To validate the proposed ultrasound-guided needle steering system, needle insertion experiments in biological tissue phantoms are performed in two cases–with and without obstacle. The results demonstrate that our needle steering strategy guides the needle to the desired targets with the maximum error of 2.85[Formula: see text]mm.

AAPM and GEC-ESTRO guidelines for image-guided robotic brachytherapy: report of Task Group 192

In the last decade, there have been significant developments into integration of robots and automation tools with brachytherapy delivery systems. These systems aim to improve the current paradigm by executing higher precision and accuracy in seed placement, improving calculation of optimal seed locations, minimizing surgical trauma, and reducing radiation exposure to medical staff. Most of the applications of this technology have been in the implantation of seeds in patients with early-stage prostate cancer. Nevertheless, the techniques apply to any clinical site where interstitial brachytherapy is appropriate. In consideration of the rapid developments in this area, the American Association of Physicists in Medicine (AAPM) commissioned Task Group 192 to review the state-of-the-art in the field of robotic interstitial brachytherapy. This is a joint Task Group with the Groupe Européen de Curiethérapie-European Society for Radiotherapy & Oncology (GEC-ESTRO). All developed and reported robotic brachytherapy systems were reviewed. Commissioning and quality assurance procedures for the safe and consistent use of these systems are also provided. Manual seed placement techniques with a rigid template have an estimated in vivo accuracy of 3-6 mm. In addition to the placement accuracy, factors such as tissue deformation, needle deviation, and edema may result in a delivered dose distribution that differs from the preimplant or intraoperative plan. However, real-time needle tracking and seed identification for dynamic updating of dosimetry may improve the quality of seed implantation. The AAPM and GEC-ESTRO recommend that robotic systems should demonstrate a spatial accuracy of seed placement ≤1.0 mm in a phantom. This recommendation is based on the current performance of existing robotic brachytherapy systems and propagation of uncertainties. During clinical commissioning, tests should be conducted to ensure that this level of accuracy is achieved. These tests should mimic the real operating procedure as closely as possible. Additional recommendations on robotic brachytherapy systems include display of the operational state; capability of manual override; documented policies for independent check and data verification; intuitive interface displaying the implantation plan and visualization of needle positions and seed locations relative to the target anatomy; needle insertion in a sequential order; robot-clinician and robot-patient interactions robustness, reliability, and safety while delivering the correct dose at the correct site for the correct patient; avoidance of excessive force on radioactive sources; delivery confirmation of the required number or position of seeds; incorporation of a collision avoidance system; system cleaning, decontamination, and sterilization procedures. These recommendations are applicable to end users and manufacturers of robotic brachytherapy systems.

Establishing implantation uncertainties for focal brachytherapy with I-125 seeds for the treatment of localized prostate cancer

BACKGROUND

The efficacy of focal continuous low dose-rate brachytherapy (CLDR-BT) for prostate cancer requires that appropriate margins are applied to ensure robust target coverage. In this study we propose a method to establish such margins by emulating a focal treatment in patients treated with CLDR-BT to the entire gland.

MATERIAL AND METHODS

In 15 patients with localized prostate cancer, prostate volumes and dominant intra-prostatic lesions were delineated on pre-treatment magnetic resonance imaging (MRI). Delineations and MRI were registered to trans-rectal ultrasound images in the operating theater. The patients received CLDR-BT treatment to the total prostate volume. The implantation consisted of two parts: an experimental focal plan covering the dominant intra-prostatic lesion (F-GTV), followed by a plan containing additional seeds to achieve entire prostate coverage. Isodose surfaces were reconstructed using follow-up computed tomography (CT). The focal dose was emulated by reconstructing seeds from the focal plan only. The distance to agreement between planned and delivered isodose surfaces and F-GTV coverage was determined to calculate the margin required for robust treatment.

RESULTS

If patients had been treated only focally, the target volume would have been reduced from an average of 40.9 cm3 for the entire prostate to 5.8 cm3 for the focal plan. The D90 for the F-GTV in the focal plan was 195±60 Gy, the V100 was 94% [range 71-100%]. The maximum distance (cd95) between the planned and delivered isodose contours was 0.48 cm.

CONCLUSIONS

This study provides an estimate of 0.5 cm for the margin required for robust coverage of a focal target volume prior to actually implementing a focal treatment protocol.

Autonomous real-time interventional scan plane control with a 3-D shape-sensing needle

This study demonstrates real-time scan plane control dependent on three-dimensional needle bending, as measured from magnetic resonance imaging (MRI)-compatible optical strain sensors. A biopsy needle with embedded fiber Bragg grating (FBG) sensors to measure surface strains is used to estimate its full 3-D shape and control the imaging plane of an MR scanner in real-time, based on the needle's estimated profile. The needle and scanner coordinate frames are registered to each other via miniature radio-frequency (RF) tracking coils, and the scan planes autonomously track the needle as it is deflected, keeping its tip in view. A 3-D needle annotation is superimposed over MR-images presented in a 3-D environment with the scanner's frame of reference. Scan planes calculated based on the FBG sensors successfully follow the tip of the needle. Experiments using the FBG sensors and RF coils to track the needle shape and location in real-time had an average root mean square error of 4.2 mm when comparing the estimated shape to the needle profile as seen in high resolution MR images. This positional variance is less than the image artifact caused by the needle in high resolution SPGR (spoiled gradient recalled) images. Optical fiber strain sensors can estimate a needle's profile in real-time and be used for MRI scan plane control to potentially enable faster and more accurate physician response.

Needle path planning and steering in a three-dimensional non-static environment using two-dimensional ultrasound images

Needle insertion is commonly performed in minimally invasive medical procedures such as biopsy and radiation cancer treatment. During such procedures, accurate needle tip placement is critical for correct diagnosis or successful treatment. Accurate placement of the needle tip inside tissue is challenging, especially when the target moves and anatomical obstacles must be avoided. We develop a needle steering system capable of autonomously and accurately guiding a steerable needle using two-dimensional (2D) ultrasound images. The needle is steered to a moving target while avoiding moving obstacles in a three-dimensional (3D) non-static environment. Using a 2D ultrasound imaging device, our system accurately tracks the needle tip motion in 3D space in order to estimate the tip pose. The needle tip pose is used by a rapidly exploring random tree-based motion planner to compute a feasible needle path to the target. The motion planner is sufficiently fast such that replanning can be performed repeatedly in a closed-loop manner. This enables the system to correct for perturbations in needle motion, and movement in obstacle and target locations. Our needle steering experiments in a soft-tissue phantom achieves maximum targeting errors of 0.86 ± 0.35 mm (without obstacles) and 2.16 ± 0.88 mm (with a moving obstacle).

Motion planning with sequential convex optimization and convex collision checking

We present a new optimization-based approach for robotic motion planning among obstacles. Like CHOMP (Covariant Hamiltonian Optimization for Motion Planning), our algorithm can be used to find collision-free trajectories from naïve, straight-line initializations that might be in collision. At the core of our approach are (a) a sequential convex optimization procedure, which penalizes collisions with a hinge loss and increases the penalty coefficients in an outer loop as necessary, and (b) an efficient formulation of the no-collisions constraint that directly considers continuous-time safety Our algorithm is implemented in a software package called TrajOpt.

We report results from a series of experiments comparing TrajOpt with CHOMP and randomized planners from OMPL, with regard to planning time and path quality. We consider motion planning for 7 DOF robot arms, 18 DOF full-body robots, statically stable walking motion for the 34 DOF Atlas humanoid robot, and physical experiments with the 18 DOF PR2. We also apply TrajOpt to plan curvature-constrained steerable needle trajectories in the SE(3) configuration space and multiple non-intersecting curved channels within 3D-printed implants for intracavitary brachytherapy. Details, videos, and source code are freely available at: http://rll.berkeley.edu/trajopt/ijrr .

THE POTENTIAL FIELD-BASED TRAJECTORY PLANNING FOR NEEDLE INSERTION IN A SOFT-TISSUE MODEL

This paper presents a 3D dynamic trajectory planning method for the insertion of a rigid needle into soft tissue. The optimal needle tip orientation was calculated in which applying an artificial potential field method to determine the 3D distribution of repulsive and attractive forces surrounding the target object and adjacent obstacles, e.g. bones, nerves, or arteries. Soft-tissue deformation occurs dynamically and continuously during the needle insertion. The trajectory planning was therefore temporally discretized, and the compartment searching method used in each time step. This trajectory planning method was then validated by a dynamic finite element method (FEM) simulation. The dynamic finite element model is built for the important displacement parameters of deformation node in needle insertion process. The Mooney–Rivlin material model combined with solid cubic element and an explicit center differencing scheme was used to compute the soft-tissue deformation at each time step and dynamically identify the target and obstacle positions. The proposed trajectory planning method can optimize the insertion path to achieve the target position while avoiding obstacles.

Development of a novel robot for transperineal needle based interventions: focal therapy, brachytherapy and prostate biopsies

PURPOSE

We report what is to our knowledge the initial experience with a new 3-dimensional ultrasound robotic system for prostate brachytherapy assistance, focal therapy and prostate biopsies. Its ability to track prostate motion intraoperatively allows it to manage motions and guide needles to predefined targets.

MATERIALS AND METHODS

A robotic system was created for transrectal ultrasound guided needle implantation combined with intraoperative prostate tracking. Experiments were done on 90 targets embedded in a total of 9 mobile, deformable, synthetic prostate phantoms. Experiments involved trying to insert glass beads as close as possible to targets in multimodal anthropomorphic imaging phantoms. Results were measured by segmenting the inserted beads in computerized tomography volumes of the phantoms.

RESULTS

The robot reached the chosen targets in phantoms with a median accuracy of 2.73 mm and a median prostate motion of 5.46 mm. Accuracy was better at the apex than at the base (2.28 vs 3.83 mm, p <0.001), and similar for horizontal and angled needle inclinations (2.7 vs 2.82 mm, p = 0.18).

CONCLUSIONS

To our knowledge this robot for prostate focal therapy, brachytherapy and targeted prostate biopsies is the first system to use intraoperative prostate motion tracking to guide needles into the prostate. Preliminary experiments show its ability to reach targets despite prostate motion.

A realistic deformable prostate phantom for multimodal imaging and needle-insertion procedures

PURPOSE

Phantoms are a vital step for the preliminary validation of new image-guided procedures. In this paper, the authors present a deformable prostate phantom for use with multimodal imaging (end-fire or side-fire ultrasound, CT and MRI) and more specifically for transperineal or transrectal needle-insertion procedures. It is made of soft polyvinyl chloride (PVC) plastic and includes a prostate, a perineum, a rectum, a soft periprostatic surrounding and embedded targets for image registration and needle-targeting. Its main particularity is its realistic deformability upon manipulation.

METHODS

After a detailed manufacturing description, the imaging and mechanical characteristics of the phantom are described and evaluated. First, the speed of sound and stress-strain relationship of the PVC material used in the phantom are described, followed by an analysis of its storage, imaging, needle-insertion force, and deformability characteristics.

RESULTS

The average speed of sound in the phantom was measured to be 1380 ± 20 m/s, while the stress-strain relationship was found to be viscoelastic and in the range of typical prostatic tissues. The mechanical and imaging characteristics of the phantom were found to remain stable at cooler storage temperatures. The phantom had clearly distinguishable morphology in all three imaging modalities, with embedded targets that could be precisely segmented, resulting in an average US-CT rigid registration error of 0.66 mm. The mobility of the phantom prostate upon needle insertion was between 2 and 4 mm, with rotations between 0° and 2°, about the US probe head.

CONCLUSION

The phantom's characteristics compare favorably with in vitro and in vivo measurements found in the literature. The authors believe that this realistic phantom could be of use to researchers studying new needle-based prostate diagnosis and therapy techniques.

Intra-operative 3D guidance and edema detection in prostate brachytherapy using a non-isocentric C-arm

PURPOSE

Brachytherapy (radioactive seed insertion) has emerged as one of the most effective treatment options for patients with prostate cancer, with the added benefit of a convenient outpatient procedure. The main limitation in contemporary brachytherapy is faulty seed placement, predominantly due to the presence of intra-operative edema (tissue expansion). Though currently not available, the capability to intra-operatively monitor the seed distribution, can make a significant improvement in cancer control. We present such a system here.

METHODS

Intra-operative measurement of edema in prostate brachytherapy requires localization of inserted radioactive seeds relative to the prostate. Seeds were reconstructed using a typical non-isocentric C-arm, and exported to a commercial brachytherapy treatment planning system. Technical obstacles for 3D reconstruction on a non-isocentric C-arm include pose-dependent C-arm calibration; distortion correction; pose estimation of C-arm images; seed reconstruction; and C-arm to TRUS registration.

RESULTS

In precision-machined hard phantoms with 40-100 seeds and soft tissue phantoms with 45-87 seeds, we correctly reconstructed the seed implant shape with an average 3D precision of 0.35 mm and 0.24 mm, respectively. In a DoD Phase-1 clinical trial on six patients with 48-82 planned seeds, we achieved intra-operative monitoring of seed distribution and dosimetry, correcting for dose inhomogeneities by inserting an average of over four additional seeds in the six enrolled patients (minimum 1; maximum 9). Additionally, in each patient, the system automatically detected intra-operative seed migration induced due to edema (mean 3.84 mm, STD 2.13 mm, Max 16.19 mm).

CONCLUSIONS

The proposed system is the first of a kind that makes intra-operative detection of edema (and subsequent re-optimization) possible on any typical non-isocentric C-arm, at negligible additional cost to the existing clinical installation. It achieves a significantly more homogeneous seed distribution, and has the potential to affect a paradigm shift in clinical practice. Large scale studies and commercialization are currently underway.

Kinematic analysis of a 5-DOF hybrid-driven MR compatible robot for minimally invasive prostatic interventions

This paper introduces the design and kinematic analysis of a 5-DOF (multiple degree of freedom) hybrid-driven MR (Magnetic Resonance) compatible robot for prostate brachytherapy. It can slip the leash of template and rely on the high precise of MR imaging. After a brief introduction on design requirements of MR compatible robot, a description of our robot structure, material selection, hybrid-driven, and control architecture are presented. Secondly, the forward kinematics equations are obtained according to the equivalent diagram of this robot, and the actual workspace can be outlined. This will help the designer to determine whether this robot can be operated in the MR core without intervention with patient. And then, the inverse kinematics equations combined with trajectory planning are used to calculate the actuators movement. This will help the control system to manipulate the robotic accurately. Finally, vision based experiments on phantoms are used to verify the mechanism precision. As the results shown, the needle tip precision of mechanism is 0.9 mm in the general lab environment.

Dose escalation to the dominant intraprostatic lesion defined by sextant biopsy in a permanent prostate I-125 implant: a prospective comparative toxicity analysis

PURPOSE

Using real-time intraoperative inverse-planned permanent seed prostate implant (RTIOP/PSI), multiple core biopsy maps, and three-dimensional ultrasound guidance, we planned a boost volume (BV) within the prostate to which hyperdosage was delivered selectively. The aim of this study was to investigate the potential negative effects of such a procedure.

METHODS AND MATERIALS

Patients treated with RTIOP/PSI for localized prostate cancer with topographic biopsy results received an intraprostatic boost (boost group [BG]). They were compared with patients treated with a standard plan (reference group [RG]). Plans were generated using a simulated annealing inverse planning algorithm. Prospectively recorded urinary, rectal, and sexual toxicities and dosimetric parameters were compared between groups.

RESULTS

The study included 120 patients treated with boost technique who were compared with 70 patients treated with a standard plan. Boost technique did not significantly change the number of seeds (55.1/RG vs. 53.6/BG). The intraoperative prostate V150 was slightly higher in BG (75.2/RG vs. 77.2/BG, p = 0.039). Urethra V100, urethra D90, and rectal D50 were significantly lower in the BG. No significant differences were seen in acute or late urinary, rectal, or sexual toxicities.

CONCLUSIONS

Because there were no differences between the groups in acute and late toxicities, we believe that BV can be planned and delivered to the dominant intraprostatic lesion without increasing toxicity. It is too soon to say whether a boost technique will ultimately increase local control.

Sensorless motion planning for medical needle insertion in deformable tissues

Minimally invasive medical procedures such as biopsies, anesthesia drug injections, and brachytherapy cancer treatments require inserting a needle to a specific target inside soft tissues. This is difficult because needle insertion displaces and deforms the surrounding soft tissues causing the target to move during the procedure. To facilitate physician training and preoperative planning for these procedures, we develop a needle insertion motion planning system based on an interactive simulation of needle insertion in deformable tissues and numerical optimization to reduce placement error. We describe a 2-D physically based, dynamic simulation of needle insertion that uses a finite-element model of deformable soft tissues and models needle cutting and frictional forces along the needle shaft. The simulation offers guarantees on simulation stability for mesh modifications and achieves interactive, real-time performance on a standard PC. Using texture mapping, the simulation provides visualization comparable to ultrasound images that the physician would see during the procedure. We use the simulation as a component of a sensorless planning algorithm that uses numerical optimization to compute needle insertion offsets that compensate for tissue deformations. We apply the method to radioactive seed implantation during permanent seed prostate brachytherapy to minimize seed placement error.

Robotic assistance for ultrasound-guided prostate brachytherapy

We present a robotically assisted prostate brachytherapy system and test results in training phantoms and Phase-I clinical trials. The system consists of a transrectal ultrasound (TRUS) and a spatially co-registered robot, fully integrated with an FDA-approved commercial treatment planning system. The salient feature of the system is a small parallel robot affixed to the mounting posts of the template. The robot replaces the template interchangeably, using the same coordinate system. Established clinical hardware, workflow and calibration remain intact. In all phantom experiments, we recorded the first insertion attempt without adjustment. All clinically relevant locations in the prostate were reached. Non-parallel needle trajectories were achieved. The pre-insertion transverse and rotational errors (measured with a Polaris optical tracker relative to the template's coordinate frame) were 0.25 mm (STD=0.17 mm) and 0.75 degrees (STD=0.37 degrees). In phantoms, needle tip placement errors measured in TRUS were 1.04 mm (STD=0.50mm). A Phase-I clinical feasibility and safety trial has been successfully completed with the system. We encountered needle tip positioning errors of a magnitude greater than 4mm in only 2 of 179 robotically guided needles, in contrast to manual template guidance where errors of this magnitude are much more common. Further clinical trials are necessary to determine whether the apparent benefits of the robotic assistant will lead to improvements in clinical efficacy and outcomes.

Parameter identification for a needle-tissue interaction model

In this work, needle-tissue interaction forces are modeled by a three parameter force distribution composed of two step functions with variable amplitudes and spacing. A finite element based simulation is used to adjust the parameters and fit the simulation results to the experimental data. In experiments, needle displacements and needle base forces were measured along with tissue displacements. A real-time version of the time-domain cross-correlation method was employed in this study to estimate the tissue displacements from ultrasound radio-frequency data, as done in elastography. In addition to the force model parameters, the elastic parameters of the tissue were adjusted to match the simulated and measured displacements.

3D simulation of needle-tissue interaction with application to prostate brachytherapy

This paper presents a needle-tissue interaction model that is a 3D extension of prior work based on needle and tissue models discretized using the Finite Element Method. The use of flexible needles necessitates remeshing the tissue during insertion, since simple mesh-node snapping to the tip can be detrimental to the simulation. In this paper, node repositioning and node addition are the two methods of mesh modification examined for coarse meshes. Our focus is on numerical approaches for fast implementation of these techniques. Although the two approaches compared, namely the Woodbury formula (matrix inversion lemma) and the boundary condition switches, have the same computational complexity, the Woodbury formula is shown to perform faster due to its cache-efficient order of operations. Furthermore, node addition is applied in constant time for both approaches, whereas node repositioning requires longer and variable computational times. A method for rendering the needle forces during simulated insertions into a 3D prostate model has been implemented. Combined with a detailed anatomical segmentation, this will be useful in teaching the practice of prostate brachytherapy. Issues related to discretization of such coupled (e.g., needle-tissue) models are also discussed.

Brachytherapy--an example of a urological minimally invasive robotic procedure

A review of the application of urological robots to the prostate is given, together with an examination of the conventional brachytherapy procedure for insertion of radioactive seeds to treat prostate cancer. The specification and design of a robotic system is provided, which can position a series of needles and radioactive pellets in accordance with a pre-operative plan. The needles can be withdrawn automatically, leaving the seeds in position.A separate motorised system is used to position a trans-rectal ultrasound (U/S) probe, which can be used to continually monitor the seed placement. The robot program can be updated intra-operatively if the U/S image shows this to be necessary. The demonstrator system has been demonstrated in-vitro using a variety of gel and animal tissue phantoms. The resulting robot performance shows this to be a viable approach.

Permanent prostate implant using high activity seeds and inverse planning with fast simulated annealing algorithm: A 12-year Canadian experience

PURPOSE

To report outcomes and toxicity of the first Canadian permanent prostate implant program.

METHODS AND MATERIALS

396 consecutive patients (Gleason < or =6, initial prostate specific antigen (PSA) < or =10 and stage T1-T2a disease) were implanted between June 1994 and December 2001. The median follow-up is of 60 months (maximum, 136 months). All patients were planned with fast-simulated annealing inverse planning algorithm with high activity seeds ([gt] 0.76 U). Acute and late toxicity is reported for the first 213 patients using a modified RTOG toxicity scale. The Kaplan-Meier biochemical failure-free survival (bFFS) is reported according to the ASTRO and Houston definitions.

RESULTS

The bFFS at 60 months was of 88.5% (90.5%) according to the ASTRO (Houston) definition and, of 91.4% (94.6%) in the low risk group (initial PSA < or =10 and Gleason < or =6 and Stage < or =T2a). Risk factors statistically associated with bFFS were: initial PSA >10, a Gleason score of 7-8, and stage T2b-T3. The mean D90 was of 151 +/- 36.1 Gy. The mean V100 was of 85.4 +/- 8.5% with a mean V150 of 60.1 +/- 12.3%. Overall, the implants were well tolerated. In the first 6 months, 31.5% of the patients were free of genitourinary symptoms (GUs), 12.7% had Grade 3 GUs; 91.6% were free of gastrointestinal symptoms (GIs). After 6 months, 54.0% were GUs free, 1.4% had Grade 3 GUs; 95.8% were GIs free.

CONCLUSION

The inverse planning with fast simulated annealing and high activity seeds gives a 5-year bFFS, which is comparable with the best published series with a low toxicity profile.

Robotic assistance for ultrasound guided prostate brachytherapy

We present a robotically assisted prostate brachytherapy system and test results in training phantoms. The system consists of a transrectal ultrasound (TRUS) and a spatially co-registered robot integrated with an FDA-approved commercial treatment planning system. The salient feature of the system is a small parallel robot affixed to the mounting posts of the template. The robot replaces the template interchangeably and uses the same coordinate system. Established clinical hardware, workflow and calibration are left intact. In these experiments, we recorded the first insertion attempt without adjustment. All clinically relevant locations were reached. Non-parallel needle trajectories were achieved. The pre-insertion transverse and rotational errors (measured with Polaris optical tracker relative to the template's coordinate frame) were 0.25mm (STD = 0.17mm) and 0.75 degrees (STD = 0.37 degrees). The needle tip placement errors measured in TRUS were 1.04mm (STD = 0.50mm). The system is in Phase-I clinical feasibility and safety trials, under Institutional Review Board approval.

Modeling of needle-tissue interaction using ultrasound-based motion estimation

A needle-tissue interaction model is an essential part of every needle insertion simulator. In this paper, a new experimental method for the modeling of needle-tissue interaction is presented. The method consists of measuring needle and tissue displacements with ultrasound, measuring needle base forces, and using a deformation simulation model to identify the parameters of a needle-tissue interaction model. The feasibility of this non-invasive approach was demonstrated in an experiment in which a brachytherapy needle was inserted into a prostate phantom. Ultrasound radio-frequency data and the time-domain cross-correlation method, often used in ultrasound elastography, were used to generate the tissue displacement field during needle insertion. A three-parameter force density model was assumed for the needle-tissue interaction. With the needle displacement, tissue displacement and needle base forces as input data, finite element simulations were carried out to adjust the model parameters to achieve a good fit between simulated and measured data.