Front-tracking image reconstruction algorithm for EIT-monitored cryosurgery using the boundary element method

Experimental observation of cavity-free ice-free isochoric vitrification via combined pressure measurements and photon counting x-ray computed tomography

Isochoric (constant-volume or volumetrically confined) vitrification has shown potential as an alternative cryopreservation-by-vitrification technique, but the complex processes at play within the chamber are yet poorly characterized, and recent investigations have prompted significant debate around whether a truly isochoric vitrification process (in which the liquid remains completely confined by solid boundaries) is indeed feasible. Based on a recent thermomechanical simulation of a high-concentration Me2SO solution, Solanki and Rabin (Cryobiology, 2023, 111, 9-15.) argue that isochoric vitrification is not feasible, because differential thermal contraction of the solution and container will necessarily drive generation of a cavity, corrupting the rigid confinement of the liquid. Here, we provide direct experimental evidence to the contrary, demonstrating cavity-free isochoric vitrification of a ∼3.5 M vitrification solution by combined isochoric pressure measurement (IPM) and photon-counting x-ray computed tomography (PC-CT). We hypothesize that the absence of a cavity is due to the minimal thermal contraction of the solution, which we support with additional volumetric analysis of the PC-CT reconstructions. In total, this study provides experimental evidence both demonstrating the feasibility of isochoric vitrification and highlighting the potential of designing vitrification solutions that exhibit minimal thermal contraction.

Imaging cryosurgery with EIT: tracking the ice front and post-thaw tissue viability

Cryosurgery employs freezing for targeted destruction of undesirable tissues such as cancer. Ice front imaging has made controlled treatment of deep body tumors possible. One promising method, recently explored for this task, is EIT, which recovers images of electrical impedance from measurements made at boundary electrodes. However, since frozen tissue near the ice front survives, ice front imaging is insufficient. Monitoring treatment effect would enable iterative cryosurgery, where extents of ablation and need for further treatment are assessed upon thawing. Since lipid bilayers are strong barriers to low frequency electrical current and cell destruction implies impaired membranes, EIT should be able to detect the desired effect of cryosurgery: cell death. Previous work has tested EIT for ice front imaging with tank studies while others have simulated EIT in detecting cryoablation, but in vivo tests have not been reported in either case. To address this, we report 3D images of differential conductivity throughout the freeze-thaw cycle in a rat liver model in vivo with histological validation, first testing our system for ice front imaging in a gel and for viability imaging post-thaw in a raw potato slice.

A boundary element approach for image-guided near-infrared absorption and scatter estimation

Multimodality NIR spectroscopy systems offer the possibility of region-based vascular and molecular characterization of tissue in vivo. However, computationally efficient 3D image reconstruction algorithms specific to these image-guided systems currently do not exist. Image reconstruction is often based on finite-element methods (FEMs), which require volume discretization. Here, a boundary element method (BEM) is presented using only surface discretization to recover the optical properties in an image-guided setting. The reconstruction of optical properties using BEM was evaluated in a domain containing a 30 mm inclusion embedded in two layer media with different noise levels and initial estimates. For 5% noise in measurements, and background starting values for reconstruction, the optical properties were recovered to within a mean error of 6.8%. When compared with FEM for this case, BEM showed a 28% improvement in computational time. BEM was also applied to experimental data collected from a gelatin phantom with a 25 mm inclusion and could recover the true absorption to within 6% of expected values using less time for computation compared with FEM. When applied to a patient-specific breast mesh generated using MRI, with a 2 cm ductal carcinoma, BEM showed successful recovery of optical properties with less than 5% error in absorption and 1% error in scattering, using measurements with 1% noise. With simpler and faster meshing schemes required for surface grids as compared with volume grids, BEM offers a powerful and potentially more feasible alternative for high-resolution 3D image-guided NIR spectroscopy.

Finite-element analysis of ex vivo and in vivo hepatic cryoablation

Cryoablation is a widely used method for the treatment of nonresectable primary and metastatic liver tumors. A model that can accurately predict the size of a cryolesion may allow more effective treatment of tumor, while sparing normal liver tissue. We generated a computer model of tissue cryoablation using the finite-element method (FEM). In our model, we considered the heat transfer mechanism inside the cryoprobe and also cryoprobe surfaces so our model could incorporate the effect of heat transfer along the cryoprobe from the environment at room temperature. The modeling of the phase shift from liquid to solid was a key factor in the accurate development of this model. The model was verified initially in an ex vivo liver model. Temperature history at three locations around one cryoprobe and between two cryoprobes was measured. The comparison between the ex vivo result and the FEM modeling result at each location showed a good match, where the maximum difference was within the error range acquired in the experiment (< 5 degrees C). The FEM model prediction of the lesion size was within 0.7 mm of experimental results. We then validated our FEM in an in vivo experimental porcine model. We considered blood perfusion in conjunction with blood viscosity depending on temperature. The in vivo iceball size was smaller than the ex vivo iceball size due to blood perfusion as predicted in our model. The FEM results predicted this size within 0.1-mm error. The FEM model we report can accurately predict the extent of cryoablation in the liver.

A narrow-band level set method applied to EIT in brain for cryosurgery monitoring

In this paper, we investigate the feasibility of applying a novel level set reconstruction technique to electrical imaging of the human brain. We focus particularly on the potential application of electrical impedance tomography (EIT) to cryosurgery monitoring. In this application, cancerous tissue is treated by a local freezing technique using a small needle-like cryosurgery probe. The interface between frozen and nonfrozen tissue can be expected to have a relatively high contrast in conductivity and we treat the inverse problem of locating and monitoring this interface during the treatment. A level set method is used as a powerful and flexible tool for tracking the propagating interfaces during the monitoring process. For calculating sensitivities and the Jacobian when deforming the interfaces we employ an adjoint formula rather than a direct differentiation technique. In particular, we are using a narrow-band technique for this procedure. This combination of an adjoint technique and a narrow-band technique for calculating Jacobians results in a computationally efficient and extremely fast method for solving the inverse problem. Moreover, due to the reduced number of unknowns in each step of the narrow-band approach compared to a pixel- or voxel-based technique, our reconstruction scheme tends to be much more stable. We demonstrate that our new method also outperforms its pixel-/voxel-based counterparts in terms of image quality in this application.

Electrical impedance tomography imaging using a priori ultrasound data

Background

Different imaging systems (e.g. electrical, magnetic, and ultrasound) rely on a wide variety of physical properties, and the datasets obtained from such systems provide only partial information about the unknown true state. One approach is to choose complementary imaging systems, and to combine the information to achieve a better representation.

Methods

This paper discusses the combination of ultrasound and electrical impedance tomography (EIT) information. Ultrasound reflection signals are good at locating sharp acoustic density changes associated with the boundaries of objects. Some boundaries, however, may be indeterminable due to masking from intermediate boundaries or because they are outside the ultrasound beam. Conversely, the EIT data contains relatively low-quality information, but it includes the whole region enclosed by the electrodes.

Results

Results are shown from a narrowband level-set method applied to 2D and 3D EIT incorporating limited angle ultrasound time of flight data.

Conclusion

The EIT reconstruction is shown to be faster and more accurate using the additional edge information from both one and four transducer ultrasound systems.