Graphics Processing Unit-Based Bioheat Simulation to Facilitate Rapid Decision Making Associated with Cryosurgery Training

Hybrid Simulation and Planning Platform for Cryosurgery with Microsoft HoloLens

Cryosurgery is a technique of growing popularity involving tissue ablation under controlled freezing. Technological advancement of devices along with surgical technique improvements have turned cryosurgery from an experimental to an established option for treating several diseases. However, cryosurgery is still limited by inaccurate planning based primarily on 2D visualization of the patient’s preoperative images. Several works have been aimed at modelling cryoablation through heat transfer simulations; however, most software applications do not meet some key requirements for clinical routine use, such as high computational speed and user-friendliness. This work aims to develop an intuitive platform for anatomical understanding and pre-operative planning by integrating the information content of radiological images and cryoprobe specifications either in a 3D virtual environment (desktop application) or in a hybrid simulator, which exploits the potential of the 3D printing and augmented reality functionalities of Microsoft HoloLens. The proposed platform was preliminarily validated for the retrospective planning/simulation of two surgical cases. Results suggest that the platform is easy and quick to learn and could be used in clinical practice to improve anatomical understanding, to make surgical planning easier than the traditional method, and to strengthen the memorization of surgical planning.

Evaluation of a numerical simulation for cryoablation - comparison with bench data, clinical kidney and lung cases

PURPOSE

The accuracy of a numerical simulation of cryoablation ice balls was evaluated in gel phantom data as well as clinical kidney and lung cases.

MATERIALS AND METHODS

To evaluate the accuracy, 64 experimental single-needle cryoablations and 12 multi-needle cryoablations in gel phantoms were re-simulated with the corresponding freeze-thaw-freeze cycles. The simulated temperatures were compared over time with the measurements of thermocouples. For single needles, temperature values were compared at each thermocouple location. For multiple needles, Euclidean distances between simulated and measured isotherms (10 °C, 0 °C, -20 °C, -40 °C) were computed. Furthermore, surface and volume of simulated 0 °C isotherms were compared to cryoablation-induced ice balls in 14 kidney and 13 lung patients. For this purpose, needle positions and relevant anatomical structures defining material parameters (kidney/lung, tumor) were reconstructed from pre-ablation CT images and fused with postablation CT images (from which ice balls were extracted by manual delineation).

RESULTS

The single-needle gel phantom cases showed less than 5 °C prediction error on average. Over all multiple needle experiments in gel, the mean and maximum isotherm distance were less than 2.3 mm and 4.1 mm, respectively. Average Dice coefficients of 0.82/0.63 (kidney/lung) and mean surface distances of 2.59/3.12 mm quantify the prediction performance of the numerical simulation. However, maximum surface distances of 10.57/10.8 mm indicate that locally larger errors have to be expected.

CONCLUSION

A very good agreement of the numerical simulations for gel experiments was measured and a satisfactory agreement of the numerical simulations with measured ice balls in patient data was shown.

A Novel, Efficient, Unit Circle-Based, Method for Positioning and Operating Cryo-Surgical Probes in Convex Target Areas

A novel method for positioning and operating needle-like cryo-surgical probes in 2D convex target areas is presented. The method is based on the recorded dynamic performance of a single probe, termed "unit circle," (UC) embedded in a semi-infinite, tissue-like medium. Up to 15 cryo-probes, inserted into the same depth, are operated uniformly for 2-5 min. A predetermined number of probes are rearranged inside the target area until a "tight configuration" is obtained. The probes are initially arranged inside the target area such that the "lethal temperature" circles produced by them are tangent to its contour and to both adjacent lethal temperature circles. Subsequently, all probes are repositioned inwardly, each at a specific distance that depends on the local radius of curvature of the target area. Resulting total "defect areas"-internal and external-for a number of demonstrated cases, amounted to between 2.5% and 7.6% of the target area. The lower values of the defect areas were obtained with increasing numbers of inserted probes coupled with shorter operating times. Possible freezing damages to regions beyond the target area were reduced by up to about 30% for these cases. Similar results were obtained for a case of combined convex-concave target area, treated with additional, externally inserted, heating probes.

GPU-based 3D iceball modeling for fast cryoablation simulation and planning

PURPOSE

The elimination of abdominal tumors by percutaneous cryoablation has been shown to be an effective and less invasive alternative to open surgery. Cryoablation destroys malignant cells by freezing them with one or more cryoprobes inserted into the tumor through the skin. Alternating cycles of freezing and thawing produce an enveloping iceball that causes the tumor necrosis. Planning such a procedure is difficult and time-consuming, as it is necessary to plan the number and cryoprobe locations and predict the iceball shape which is also influenced by the presence of heating sources, e.g., major blood vessels and warm saline solution, injected to protect surrounding structures from the cold.

METHODS

This paper describes a method for fast GPU-based iceball modeling based on the simulation of thermal propagation in the tissue. Our algorithm solves the heat equation within a cube around the cryoprobes tips and accounts for the presence of heating sources around the iceball.

RESULTS

Experimental results of two studies have been obtained: an ex vivo warm gel setup and simulation on five retrospective patient cases of kidney tumors cryoablation with various levels of complexity of the vascular structure and warm saline solution around the tumor tissue. The experiments have been conducted in various conditions of cube size and algorithm implementations. Results show that it is possible to obtain an accurate result within seconds.

CONCLUSION

The promising results indicate that our method yields accurate iceball shape predictions in a short time and is suitable for surgical planning.

Defeating Cancers' Adaptive Defensive Strategies Using Thermal Therapies: Examining Cancer's Therapeutic Resistance, Ablative, and Computational Modeling Strategies as a means for Improving Therapeutic Outcome

BACKGROUND

Diverse thermal ablative therapies are currently in use for the treatment of cancer. Commonly applied with the intent to cure, these ablative therapies are providing promising success rates similar to and often exceeding "gold standard" approaches. Cancer-curing prospects may be enhanced by deeper understanding of thermal effects on cancer cells and the hosting tissue, including the molecular mechanisms of cancer cell mutations, which enable resistance to therapy. Furthermore, thermal ablative therapies may benefit from recent developments in computer hardware and computation tools for planning, monitoring, visualization, and education.

METHODS

Recent discoveries in cancer cell resistance to destruction by apoptosis, autophagy, and necrosis are now providing an understanding of the strategies used by cancer cells to avoid destruction by immunologic surveillance. Further, these discoveries are now providing insight into the success of the diverse types of ablative therapies utilized in the clinical arena today and into how they directly and indirectly overcome many of the cancers' defensive strategies. Additionally, the manner in which minimally invasive thermal therapy is enabled by imaging, which facilitates anatomical features reconstruction, insertion guidance of thermal probes, and strategic placement of thermal sensors, plays a critical role in the delivery of effective ablative treatment.

RESULTS

The thermal techniques discussed include radiofrequency, microwave, high-intensity focused ultrasound, laser, and cryosurgery. Also discussed is the development of thermal adjunctive therapies-the combination of drug and thermal treatments-which provide new and more effective combinatorial physical and molecular-based approaches for treating various cancers. Finally, advanced computational and planning tools are also discussed.

CONCLUSION

This review lays out the various molecular adaptive mechanisms-the hallmarks of cancer-responsible for therapeutic resistance, on one hand, and how various ablative therapies, including both heating- and freezing-based strategies, overcome many of cancer's defenses, on the other hand, thereby enhancing the potential for curative approaches for various cancers.

The role of exposure time in computerized training of prostate cryosurgery: performance comparison of surgical residents with engineering students

PURPOSE

This study aims at the evaluation of a prototype of a computerized trainer for cryosurgery-the controlled destruction of cancer tumors by freezing. The hypothesis in this study is that computer-based cryosurgery training for an optimal cryoprobe layout is essentially a matter of exposure time, rather than trainee background or the specific computer-generated planning target. Key geometric features under considerations are associated with spatial limitations on cryoprobes placement and the match between the resulted thermal field and the unique anatomy of the prostate.

METHODS

All experiments in this study were performed on the cryosurgery trainer-a prototype platform for computerized cryosurgery training, which has been presented previously. Among its key features, the cryosurgery trainer displays the prostate shape and its contours and provides a distance measurement tool on demand, in order to address spatial constraints during ultrasound imaging guidance. Another unique feature of the cryosurgery trainer is an output movie, displaying the simulated thermal field at the end of the cryoprocedure.

RESULTS

The current study was performed on graduate engineering students having no formal background in medicine, and the results were benchmarked against data obtained on surgical residents having no experience with cryosurgery. Despite fundamental differences in background and experience, neither group displayed superior performance when it comes to cryoprobe layout planning.

CONCLUSIONS

This study demonstrates that computer-based training of an optimal cryoprobe layout is feasible. This study demonstrates that the training quality is essentially related to the training exposure time, rather than to a specific planning strategy from those investigated.

Computerized Planning of Prostate Cryosurgery and Shape Considerations

The current study aims to explore possible relationships between various prostate shapes and the difficulty in creating a computer-based plan for cryosurgery. This research effort is a part of an ongoing study to develop computational means in order to improve cryosurgery training and education. This study uses a computerized planner-a key building block of a recently developed prototype for cryosurgery training. The quality of planning is measured by the overall defect volume, a proprietary concept which refers to undercooled areas internal to the target region and overcooled areas external to it. Results of this study numerically confirm that the overall defect volume decreases with an increasing number of cryoprobes, regardless of the geometry of the prostate. However, the number of cryoprobes required to achieve the smallest possible defect may be unrealistically high (<30). Results of this study also demonstrate that the optimal cryoprobe layout is associated with a smaller defect for symmetric prostate geometries and, independently, for prostate models that better resemble a sphere. Furthermore, a smaller defect is typically achieved when the urethra passes through the center of the prostate model. This study proposes to create a cryoprobe convex hull for the purpose of initial planning, which is a subdomain similar in shape to the prostate but at a reduced size. Parametric studies indicate that a cryoprobe convex hull contracted by 7 to 9 mm in all directions from the prostate capsule serves as a quasi-optimal initial condition for planning, that is, a preselected number of cryoprobes placed in the cryoprobe convex hull yields favorable results for optimization. The cryoprobe convex hull could accelerate computer-based planning, while also being adopted as a concept for traditional cryosurgery training, when computerized means are absent.

A Computerized Tutor Prototype for Prostate Cryotherapy: Key Building Blocks and System Evaluation

This paper focuses on the evaluation of a prototype for a computer-based tutoring system for prostate cryosurgery, while reviewing its key building blocks and their benchmark performance. The tutoring system lists geometrical constraints of cryoprobe placement, displays a rendered shape of the prostate, simulates cryoprobe insertion, enables distance measurements, simulates the corresponding thermal history, and evaluates the mismatch between the target region shape and a pre-selected planning isotherm. The quality of trainee planning is measured in comparison with a computer-generated plan, created for each case study by a previously developed planning algorithm, known as bubble-packing. While the tutoring level in this study aims only at geometrical constraints on cryoprobe placement and the resulting thermal history, it creates a unique opportunity to gain insight into the process outside of the operation room. System validation of the tutor has been performed by collecting training data from surgical residents, having no prior experience or advanced knowledge of cryotherapy. Furthermore, the system has been evaluated by graduate engineering students having no formal education in medicine. In terms of match between a planning isotherm and the target region shape, results demonstrate medical residents' performance improved from 4.4% in a pretest to 37.8% in a posttest over a course of 50 minutes of training (within 10% margins from a computer-optimized plan). Comparing those results with the performance of engineering students indicates similar results, suggesting that planning of the cryoprobe layout essentially revolves around geometric considerations.

GPU-Based Simulation of Ultrasound Imaging Artifacts for Cryosurgery Training

This study presents an efficient computational technique for the simulation of ultrasound imaging artifacts associated with cryosurgery based on nonlinear ray tracing. This study is part of an ongoing effort to develop computerized training tools for cryosurgery, with prostate cryosurgery as a development model. The capability of performing virtual cryosurgical procedures on a variety of test cases is essential for effective surgical training. Simulated ultrasound imaging artifacts include reverberation and reflection of the cryoprobes in the unfrozen tissue, reflections caused by the freezing front, shadowing caused by the frozen region, and tissue property changes in repeated freeze-thaw cycles procedures. The simulated artifacts appear to preserve the key features observed in a clinical setting. This study displays an example of how training may benefit from toggling between the undisturbed ultrasound image, the simulated temperature field, the simulated imaging artifacts, and an augmented hybrid presentation of the temperature field superimposed on the ultrasound image. The proposed method is demonstrated on a graphic processing unit at 100 frames per second, on a mid-range personal workstation, at two orders of magnitude faster than a typical cryoprocedure. This performance is based on computation with C++ accelerated massive parallelism and its interoperability with the DirectX-rendering application programming interface.

A new method for temperature-field reconstruction during ultrasound-monitored cryosurgery using potential-field analogy

The current study aims at developing computational tools in order to gain information about the thermal history in areas invisible to ultrasound imaging during cryosurgery. This invisibility results from the high absorption rate of the ultrasound energy by the frozen region, which leads to an apparent opacity in the cryotreated area and a shadow behind it. A proof-of-concept for freezing-front estimation is demonstrated in the current study, using the new potential-field analogy method (PFAM). This method is further integrated with a recently developed temperature-field reconstruction method (TFRM) to estimate the temperature distribution within the frozen region. This study uses prostate cryosurgery as a developmental model and trans-rectal ultrasound imaging as a choice of practice. Results of this study indicate that the proposed PFAM is a viable and computationally inexpensive solution to estimate the extent of freezing in the acoustic shadow region. Comparison of PFAM estimations and experimental data shows an average mismatch of less than 2 mm in freezing-front location, which is comparable to the uncertainty in ultrasound imaging. Comparison of the integrated PFAM + TFRM scheme with a full-scale finite-elements analysis (FEA) indicates an average mismatch of 0.9 mm for the freezing front location and 0.1 mm for the lethal temperature isotherm of -45 °C. Comparison of the integrated PFAM + TFRM scheme with experimental temperature measurements show a difference in the range of 2 °C and 6 °C for selected points of measurement. Results of this study demonstrate the integrated PFAM + TFRM scheme as a viable and computationally inexpensive means to gain information about the thermal history in the frozen region during ultrasound-monitored cryosurgery.

Simulation-Based Cryosurgery Training: Variable Insertion Depth Planning in Prostate Cryosurgery

A proof-of-concept for an advanced-level computerized training tool for cryosurgery is demonstrated, based on three-dimensional cryosurgery simulations and a variable insertion depth strategy for cryoprobes. The objective for system development is two-fold: to identify a cryoprobe layout in order to best match a planning isotherm with the target region shape and to verify that cryoprobe placement does not violate accepted geometric constraints. System validation has been performed by collecting training data from 17 surgical residents having no prior experience or advanced knowledge of cryosurgery. This advanced-level study includes an improved training session design in order to enhance knowledge dissemination and elevate participant motivation to excel. In terms of match between a planning isotherm and the target region shape, results of this demonstrate trainee performance improvement from 4.4% in a pretest to 44.4% in a posttest over a course of 50 minutes of training. In terms of combined performance, including the above-mentioned geometrical match and constraints on cryoprobe placement, this study demonstrates trainee performance improvement from 2.2% in the pretest to 31.1% in the posttest. Given the relatively short training session and the lack of prior knowledge, these improvements are significant and encouraging. These results are of particular significance, as they have been obtained from a surgical resident population which are exposed to the typical stress and constraints in advanced surgical education.