Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 14, Article number: 5599 (2024 ) Cite this article Portable Radiation
![X-ray dark-field computed tomography for monitoring of tissue freezing | Scientific Reports](http://cdn.globalso.com/ergodi-radiation/15a6ba395.png)
Accurately monitoring the extent of freezing in biological tissue is an important requirement for cryoablation, a minimally invasive cancer treatment that induces cell death by freezing tissue with a cryoprobe. During the procedure, monitoring is required to avoid unnecessary harm to the surrounding healthy tissue and to ensure the tumor is properly encapsulated. One commonly used monitoring method is attenuation-based computed tomography (CT), which visualizes the ice ball by utilizing its hypoattenuating properties compared to unfrozen tissue. However, the contrast between frozen and unfrozen tissue remains low. In a proof-of-principle experiment, we show that the contrast between frozen and unfrozen parts of a porcine phantom mimicking breast tissue can be greatly enhanced by acquiring X-ray dark-field images that capture the increasing small-angle scattering caused by the ice crystals formed during the procedure. Our results show that, compared to X-ray attenuation, the frozen region is detected significantly better in dark-field radiographs and CT scans of the phantom. These findings demonstrate that X-ray dark-field imaging could be a potential candidate for improved monitoring of cryoablation procedures.
Cryoablation, a local ablative technique used in interventional oncology mainly for malignancies of the liver, kidney, prostate, lung and breast, has recently received increased interest1. In this procedure, cell death is induced by freezing tissue using a cryoprobe inserted percutaneously near the tumor. Compared to conventional surgical resection, which can lead to chronic pain, infection, hematoma and seroma formation, cryoablation has several advantages. Due to the natural analgesic effect of cold, it is usually painless and can be performed with only local anesthesia. Scarring is reduced due to the small needle hole and the patient can be discharged the same day2. However, the extent of the treated regions must be monitored during the procedure to minimize damage to healthy tissue while freezing the tumor region with a safety margin of 8–10 mm3,4. Unnecessarily large frozen regions have been shown to lead to fat necrosis, infection, and skin injury5,6.
For the treatment of breast cancer, the procedure has been repeatedly shown to be an effective option for tumors smaller than 15 mm in size7,8,9. The most popular methods for monitoring breast cryoablation procedures are ultrasound (US) and attenuation-based computed tomography (CT)2,10, although magnetic resonance imaging (MRI)-guided procedures have also been successfully performed11. While US-guided monitoring allows visualization of the anterior surface of the frozen region, the posterior margin cannot be resolved due to acoustic shadowing8,12. For this reason, CT or MRI guidance is recommended in addition to US to provide a safe distance from critical structures in the thorax11. Real-time US imaging has been reported as the most operator-dependent technique, as it requires considerable practice and knowledge of the artifacts that occur10.
Because soft tissue contrast is inherently low for X-ray attenuation imaging, imaging techniques that use the phase shift and/or scattering of the X-ray wavefront passing through the tissue as an additional contrast mechanism are the subject of ongoing research. Such approaches include propagation-based imaging13,14,15, which has been successfully applied for breast radiographies and CTs at synchrotron sources, as well as grating-based imaging (GBI)16,17,18,19 and coded apertures20,21. The latter two approaches have been demonstrated to be compatible with conventional X-ray sources and detectors. GBI has the additional advantage of simultaneously providing a dark-field, phase-contrast, and attenuation signal in a single scan by modulating the X-ray beam using up to three gratings with absorbing or phase-shifting properties22. In this context, dark-field is defined as a signal associated with the local visibility reduction of the interferometer due to local small-angle scattering caused by spatially random unresolved microstructure.
In recent years, it has been shown that the additional information provided by phase-contrast and dark-field modalities is beneficial for the detection of lesions and microcalcifications in the breast23,24,25,26,27. However, the dark-field signal may not only be useful for diagnosis, but also for the monitoring of freezing processes: We propose that frozen tissue, causing strong X-ray scattering due to the presence of ice crystals, should be well distinguishable from unfrozen tissue, which causes comparatively little scattering. This seems plausible as the effect has been previously used for distinguishing between frozen and unfrozen fruit28. In the currently used attenuation-based CT-guidance, only the slight change in the density of the frozen tissue, which leads to a hypoattenuating region, is used as a contrast mechanism29.
The potential of X-ray dark-field monitoring for the freezing of tissue models has so far not been investigated in either radiography or computed tomography. In this work, we show that the extent of the frozen region in a porcine phantom mimicking breast tissue can be monitored more accurately with a grating interferometer at a laboratory X-ray source than with attenuation-based imaging alone. For this purpose, the phantom is actively frozen using a pipe connected to a liquid nitrogen reservoir. The imaging results are first demonstrated in radiographs and then transferred to a CT setup.
First, a piece of pork neck is frozen using a copper rod cooled by liquid nitrogen, the setup being illustrated in Fig. 1. The process is monitored using grating-based radiographs acquired using 11 grating-phase steps, each with an exposure time of 6 s, resulting in a total acquisition duration of 66 s per radiograph. Three of the resulting pairs of attenuation and dark-field images at different times are shown along with their accompanying line plots in Fig. 2. It may be observed that the attenuation images show no change over time while the dark-field images clearly show the expansion of the frozen region around the needle. This is most likely caused by the formation of ice crystals in the tissue, leading to increased small-angle scattering and therefore increased dark-field signal. The signal shows a decreasing intensity from the rod outwards, corresponding to the density of ice crystals. The line plot of the attenuation shows no significant change. A small lateral shift of the intensity values corresponding to the copper rod is most likely caused by a slight contraction of the sample during the freezing. In the dark-field line plot, the expansion of the frozen region is easy to detect as well. The average dose per radiograph is 4.4 mGy—each radiograph simultaneously providing information for both a dark-field and attenuation image.
Sketch of the experimental setup used for the measurements. The X-ray source is followed by the source grating G\(_0\) and the phase-shifting grating G\(_1\) . Next, the beam passes through the sample tube, which is either cooled in the center by liquid nitrogen flowing from a basin through a thin pipe (CT experiments) or using a copper rod connected to a liquid nitrogen reservoir on one end (radiography experiments). The X-rays then reach the analyzer grating G\(_2\) and the beam is absorbed by a detector. The distances l and d mark the distances G\(_0\) –G\(_1\) and G\(_1\) –G\(_2\) respectively.
Series of radiographs (left) and line plots (right) of a pork neck slab frozen using a copper rod connected to a liquid nitrogen reservoir on one end. While the extent of the frozen region is clearly visible in the dark-field images, it is not visually detectable in the attenuation images. The line plots show the signal change over time in the dotted white rectangular region, averaged in the vertical direction.
Next, the process is transferred from radiography to three-dimensional imaging. For this purpose, ten consecutive grating-based CT scans are performed while a piece of pork belly is being frozen. Compared to the radiography, a significantly lower exposure time of 0.12 s per phase step is chosen to enable sufficient time resolution. For each CT scan, the sample completes a 360\(^\circ \) rotation within 2.4 min; the mean glandular dose is estimated at 22 mGy. Freezing is induced by the flow of liquid nitrogen through a thin tube pierced through the center of the sample. Figure 3 shows an axial slice of the sample over time in the attenuation and dark-field modalities together with plots of the intensity values along the dotted white line. In the dark-field, a radially symmetric increase in signal intensity over time is observed near the wall of the central cooling tube. A change in signal intensity is also observed in the attenuation images, although it is more difficult to detect. Here, a hypoattenuating region extends radially outward from the cooling tube. As indicated by the red arrows, the boundary of this region is less sharply defined than the boundary seen in the corresponding dark-field images. Histograms of the image slices in the first and third column are provided in Fig. 4 to show that the dark-field signal experiences a larger change in contrast compared to the attenuation signal, regardless of the specific windowing chosen in Fig. 3.
Series of axial slices (left) and line plots (right) near the top of the sample in dark-field and attenuation images at different time points. A radially symmetric increase in signal strength is observed in the dark-field while the same region appears slightly darker in the attenuation images. The red arrows mark the boundary of the frozen region, which is easier to see in dark-field. The acquisition time for each multimodal CT scan is 2.4 min. The line plots show the intensity along the dotted white line, averaged over 10 axial slices extending beyond the depicted image plane.
Overlayed histograms for the attenuation image slice (left) and extinction coefficient image slice (right) from Fig. 3. The image slices correspond to the leftmost (orange) and rightmost (blue) column of Fig. 3 respectively. While the attenuation signal shows a slight shift of values towards the left, the extinction values are shifted by a larger magnitude towards the right (red arrow).
The dark circle in the center of all images corresponds to the inside of the polytetrafluoroethylene (PTFE) cooling tube, which is mostly filled with evaporating liquid nitrogen. Its walls produce a visible signal in both the attenuation and dark-field images. In terms of artifacts, the dark-field images show streaks caused by the strong scattering signal present in the frozen region. They are approximately radially symmetric and brighten the background outside the frozen region over time. The attenuation images display streaking artifacts in the form of two faint streaks running horizontally from the longer sides of the elliptically compressed rod.
To examine how easily freezing is detected in the two imaging modalities, a region of interest as shown in Fig. 5a is selected and analyzed. The maximum radius corresponds to the largest distance at which signs of freezing were observed during the time window of the experiment. The region excludes the elliptical region bounded by the edge of the cooling tube, which is not relevant for the freezing analysis. To create a region with mostly homogeneous attenuation coefficients, the fat region (dark) is also excluded. This allows tracking changes in one specific tissue without creating an overlap with results from other tissues. Figure 5b shows a scatter plot of the attenuation and dark-field signal of each voxel in the region of interest at time \(t = {0}\,{\hbox {min}}\) and at \(t = {7.2}\, {\hbox {min}}\) . The mean extinction coefficient value at the beginning of the experiment is \(-\,0.002 \pm 0.002\) a.u. During the freezing process, higher extinction coefficient values are reached and the distribution widens in the vertical direction. Overall, the center of mass of the distribution is significantly (\(> 2 \sigma )\) shifted upward, where \(\sigma \) corresponds to the standard deviation. A shift of the center of mass towards lower attenuation coefficients is also observed. However, this shift is within \(2\sigma \) of the original distribution. These differences in the magnitude of the shift underline numerically that the freezing process is easier to detect in the dark-field.
Analysis of the change in attenuation and extinction coefficients over time in muscle tissue. (a) Region of interest (red) used for the analysis depicted in (b). The cooling tube and a neighboring fat region (dark) remain unselected to create a homogeneous region. It extends over 10 axial slices beyond the depicted image plane. (b) Plot of the attenuation and extinction coefficient values of all voxels in the region of interest defined in (a) at time \(t = {0}\,{\hbox {min}}\) (orange) and after 7.2 min (blue). The dashed gray lines indicate the positions of the mean signal in each channel \(\pm 2\sigma \) at \(t = {0}\,{\hbox {min}}\) , \(\sigma \) being the standard deviation. During the freezing process, the original center of mass (white cross) moves to a new position (red cross) by a vertical distance greater than \(2\sigma \) . This means that the average extinction coefficient experiences a significant increase. In comparison, the average attenuation coefficient decreases insignificantly and remains below \(2\sigma \) .
This result is then used as a method to classify a given voxel as frozen depending on its dark-field value. For this purpose, a threshold is defined by adding \(2\sigma \) to the mean of the extinction coefficients of all voxels at time \(t = {0}\,{\hbox {min}}\) , excluding the sample tube edge and the cooling pipe. In Fig. 6, a binary map (blue) of all dark-field voxels that satisfy this condition at \(t = {7.2}\,{\hbox {min}}\) is superimposed on the corresponding attenuation image. The result is a single image that combines anatomical (attenuation) and functional (dark-field) information.
Fusion image of dark-field and attenuation for an example image slice. The binary map (blue) marks voxels as frozen if they experienced an increase in the extinction coefficient of more than \(2\sigma \) compared to the mean of the extinction coefficients of all voxels at the beginning (excluding the tube border and the cooling pipe). A morphological closing with a (3,3) kernel was performed to reduce the impact of noise on the visual appearance. Dark-field information from inside the cooling pipe and the edge of the sample tube was discarded using circular masks, as it is not caused by freezing, but by the scattering properties of the materials.
In this proof-of-concept work, we have examined the suitability of X-ray attenuation and dark-field signals for tracking an active freezing process in a porcine phantom mimicking breast tissue. We have come to the conclusion that dark-field imaging allows for detecting the frozen region with higher significance than attenuation alone. Compared to attenuation, which shows little change during the freezing process, the dark-field signal has a contrast mechanism better suited to the monitoring task: At room temperature, muscle and fat tissue do not cause significant small-angle scattering and therefore provide little dark-field signal, which changes when ice crystals form during freezing. The method was exemplarily tested on tissue mimicking the composition of the human breast, but may offer similar advantages for the monitoring of freezing in other soft tissue areas (e.g. liver, kidney, prostate). It is likely less suitable for tissue types that already create strong small-angle scattering in the unfrozen state, as is the case for lung tissue30. A limitation of our study is the studied sample thickness, which required at most 3 cm depth of X-ray penetration. For a successful translation to the monitoring of real cryoablation procedures, larger sample sizes will need to be studied for their impact on image quality and measurement times.
Moreover, we showed that the dark-field information can be superimposed on the simultaneously acquired attenuation image to provide an image that combines anatomical and functional information, similar to the color overlays used in Doppler ultrasonography31. In this experiment, a threshold of \(2\sigma \) was used to classify a voxel as frozen. A possible improvement to the visualization would be to display images at multiple thresholds depending on the application: A lower threshold image may show more frozen voxels than are actually present, to avoid accidentally freezing adjacent healthy tissue. A higher threshold image provides even more confidence that the voxels containing the tumor are sufficiently frozen.
In the radiography experiment, the freezing occurred more rapidly on the left side of the needle than on the right (see Fig. 2). This effect is most likely explained by the higher proportion of fat tissue on the right side (visible as a darker region in the attenuation channel), which has a thermal conductivity half that of muscle tissue32. This is consistent with the results shown for the CT freezing experiment in Fig. 6, where the frozen region extends further out in muscle tissue compared to fat tissue, as would be expected for the lower thermal conductivity of fat.
It is to be noted that the dark-field CT measurements yield slightly negative extinction coefficients for some pixels. This can be attributed in part to the high amount of noise in the measured projections: If the visibility of a certain pixel in the reference image is decreased due to noise, the visibility may increase by chance in the corresponding sample pixel. This then leads to a positive \(V_\text {sample}/V_\text {ref}\) quotient, interpreted as a negative extinction coefficient. Another contributing factor is the fact that extinction coefficients obtained via least-squares fitting of the stepping curve are biased: They are underestimated in measurements with low photon counts33, as is the case in our experiment.
A disadvantage of dark-field CT-guided imaging and CT-guided imaging in general compared to ultrasound and MRI is the ionizing radiation to which patients are exposed. In our experiment, each CT scan yielded a dose of 22 mGy, but multiple scans are needed to monitor the procedure over time. As our experiment only constitutes a proof of principle, a number of steps still need to be taken to decrease the dose. Firstly, the exposure time could be reduced by increasing the effective pixel size via additional binning. Since the medically relevant distances for safety margins to organs are in the range of 8–10 mm3,4, a doubling of the effective pixel size from 290 to 580 \(\upmu {\hbox {m}}\) may still be acceptable. Additionally, the total dose can be lowered by increasing the time in between CT scans, sacrificing temporal resolution, or by decreasing the exposure time further while implementing self-supervised deep-learning denoising to differentiate the frozen border even in very high-noise settings34. A limitation of the dose estimate used for this experiment is that the estimate is calculated using a dosimetry method developed for the human breast, although the sample is only similar in composition.
Recently, the feasibility of a human-scale dark-field CT scanner has been demonstrated by adding phase-shifting gratings to the source and absorption gratings near the detector of an existing commercial CT35. Measurement times of few seconds have been achieved in this approach by using the intrinsic vibrations of a commercial CT scanner to sample the phase-stepping curve. This demonstrates that the time resolution of 2.4 min of our CT experiment, which was limited due to the applied measurement protocol, can in principle be drastically improved. Factors helping with the adoption of the proposed cryoablation monitoring method would be the equipping of existing, conventional CT scanners with gratings or the development of a dark-field CT arm. Attenuation-based breast CTs, which have been shown to be feasible for breast cancer screening36, could also be equipped with gratings to allow intraoperative use for monitoring.
The experimental setup is shown in Fig. 1 and uses a rotating anode (MicroMax HF007, Rigaku Corporation, Tokyo, Japan) with a molybdenum target as an X-ray source. The variables G\(_i\) refer to the i-th grating in the beam direction. The gratings for both experiments are detailed in Table 1. The geometries of the setups used for radiograph and CT measurements differ slightly due to an upgrade performed to the setup to enable fast CT imaging. The distances and resulting geometric magnifications are provided in Table 2.
For the radiography experiment, the source is operated with a voltage of 33 kVp and a current of 24 mA. A scintillation-based energy integrating detector (Dexela 1512, PerkinElmer, Inc., Waltham, United States) with a pixel pitch of 74.8 \(\upmu {\hbox {m}}\) and 1944 \(\times \) 1536 pixels is used. The effective pixel size of the setup in this configuration is 53.2 \(\upmu {\hbox {m}}\) .
For the CT experiment, the source is set to a voltage of 50 kVp and a current of 24 mA. The projections are acquired using a prototype photon-counting detector (Santis 0808 GaAs HR, DECTRIS AG, Baden, Switzerland) with a module thickness of 500 \(\upmu {\hbox {m}}\) and 1030 \(\times \) 1024 pixels with a pitch of 75 \(\upmu {\hbox {m}}\) . Due to the geometric magnification of the setup, the effective pixel size is 58 \(\upmu {\hbox {m}}\) . The grating is moved using a P-620.1CD piezoelectric actuator controlled by an E-709.1C1L digital controller (both Physik Instrumente (PI) GmbH & Co. KG, Karlsruhe, Germany).
Instead of the sample tube shown in Fig. 1, a piece of pork neck is held in place by two polycarbonate plates. To begin the experiment, a copper rod is inserted into the sample and cooled using a polystyrene bath filled with liquid nitrogen. Grating-based radiographs are then taken at regular intervals: Each radiograph consists of 11 phase steps of 6 s exposure time each. The exposure time is chosen to result in a mean glandular dose estimate on the order of the mean glandular dose in two-view digital mammography37. The frame time \(T_\text {F}\) is approximately equal to the exposure time \(T_\text {E} = {6} \,{\hbox {s}}\) , because the idle time \(T_\text {I}\) of the detector is chosen to be negligibly small. This differs from the CT measurement protocol described in detail in the next subsection. A total of 20 radiographs are taken over a period of 22 min. The time marks in the Fig. 2 refer to the starting time of the acquisition of a new image.
In order to capture the changes in the sample over time in CT scans, a measurement protocol similar to Zanette et al.38 is used. The piezoelectric actuator that controls the position of the \(\text {G}_1\) grating is moved continuously in a ramp-like pattern over time while synchronized to the detector. Per 360\(^\circ \) rotation of the sample, the detector records a total of 650 frames with a frame time of \(T_\text {F}\) = 0.22 s each, resulting in a full sample rotation duration of 2.4 min. Four additional frames from the next round are used in the processing step to obtain the last projection. Each frame consists of \(T_\text {E}\) = 0.12 s exposure time and \(T_\text {I}\) = 0.10 s idle time. The idle time, during which no photons are recorded, is present to allow the actuator to return to its original position after the fourth exposure without imaging the resulting motion artifacts. Each frame contains the same idle time to ensure equal spacing of the recorded images. The value of the idle time is chosen to be as small as possible considering the maximum motion speed of the actuator. The exposure time is chosen as a trade-off between the available image statistics and the need for fast rotations to resolve tissue changes over time. Two actuator ramps and the corresponding detector actions are shown in Fig. 7. Due to the periodicity of the phase grating, stepper position #5 produces the same pattern as position #0.
Illustration of the measurement protocol used for the CT scans. The graph shows the movement of the piezoelectric actuator over time while the sample is being rotated. The orange dot indicates the start of a new detector frame of duration \(T_\text {F}\) , which consists of an exposure time \(T_\text {E}\) and an idle time \(T_\text {I}\) . Photons are only counted during the exposure time. Per 360\(^\circ \) rotation of the sample, 130 ramps are produced by the actuator. The stepping curve provided by 5 consecutive frames is later used in the processing step to calculate one projection. By shifting the starting point by one frame and using the consecutive frames again, the next projection is calculated in a sliding-window fashion.
Before beginning the recording, liquid nitrogen is filled into the polystyrene basin at the top of the FalconTM tube. It then flows through a hollow polytetrafluoroethylene (PTFE) tube with an outer diameter of 3 mm and an inner diameter of 2 mm, which is placed inside the sample and sealed at the bottom with silicone to prevent direct contact of the liquid nitrogen with the tissue. This tube material was chosen to avoid streaking artifacts that would occur with a metal rod. Once the reservoir is filled, the measurement is initiated by recording 10 full ramp movements of the actuator without the sample in the beam to be used for the flat field correction later on. The sample is then moved into the beam path, the rotation is started and the measurement is performed as described earlier. In our experiment, the sample is rotated continuously for 10 rounds of 360\(^\circ \) .
After the flat-field correction39, a weighted least-squares fit of sine curves to the measurements in each pixel is performed to extract the attenuation, dark-field, and phase-contrast signals. The dark-field signal D is calculated as
where \(V_\text {sample}\) is the visibility in a given pixel with the sample in the beam and \(V_\text {ref}\) is the visibility without the sample. The attenuation signal is obtained analogously by substituting the visibilities for the measured photon counts. For the CT measurements, the phase retrieval is performed in a sliding-window fashion: The recorded detector frames numbered #0 to #4 are used to obtain projection #0, then the frames #1 to #5 yield the next processed projection #1, as shown in Fig. 7. The phase retrieval is additionally performed in a patch-wise manner, with the patches consisting of the center pixel and its 4 direct neighbors, to mitigate the effects of noise on the phase retrieval40. After this step, the data is binned by a factor of 5, resulting in a new effective pixel size of 290 \(\upmu {\hbox {m}}\) . A ring artifact correction is applied to all signals in the projection space. The dark-field signal is corrected as described in a subsequent subsection to mitigate the effects of beam hardening. A bilateral filter as described by Allner et al. is applied to both signals to reduce the amount of noise present due to the low dose while preserving most of the edges41. For the attenuation image, the phase-contrast data is provided to the bilateral filter algorithm as additional edge information.
Because the sensor of the energy-integrating detector used for the radiographs is affected by electronic noise, a dark-current correction is applied in addition to the flat-field correction. The radiographs are neither binned nor filtered, i.e., the resolution corresponds to the effective pixel size.
For the radiographs, an unpreserved piece of pork neck approximately \(8 \times 4 \times 3\) cm in dimensions is used. It is held in place using a holder made of two polycarbonate plates with an inner distance of 3 cm. The inserted copper needle has a diameter of 3 mm. The sample is shown in measurement position with the inserted needle in Fig. 8a.
Photographs of the samples used in the experiments. (a) Piece of pork neck in measurement position, as used for the radiography experiment. The inserted copper needle is visible in the middle of the sample; it is cooled using the polystyrene basin on top of the polycarbonate sample holder. (b) Piece of pork belly used for the CT experiment after removal of the cooling tube. The orange arrow indicates the top-facing part of the sample.
A similar sample is chosen for the CT scans: A piece of unpreserved pork belly is cut to fit into a FalconTM tube (Fisher Scientific GmbH, Schwerte, Germany) 3 cm in diameter. It comprises comparably large sections of adipose tissue and muscle and additionally includes a layer of skin. The sample is shown in Fig. 8b.
In order to remove the spectrally induced dark-field signal caused by beam-hardening from the measurements, a calibration as described by Pelzer et al. is performed for the CT measurements42. For this purpose, different Polyoxymethylene (POM) blocks ranging in thickness from 7 to 75 mm are placed in the setup, replacing the sample tube in Fig. 1. POM is chosen to mimic the composition of the sample and is comparable to a 1:1-mixture of muscle and fat43. Radiographs with an exposure time of 4 s are acquired for each thickness, yielding a transmission and dark-field value. The measured dark-field values D and transmission values T in each the image are meaned and fitted using the empirical relation
to obtain a look-up table. a and b are the fitting coefficients, for which the values 1.1 and 0.33 are obtained respectively when averaging over the used detector module for more reliable fitting. These values are used to correct the dark-field data in the CT scan using43
where \(D^\text {(SAS)}\) is the component of the dark-field signal representing small-angle scattering and D is the measured dark-field signal, \(D^\text {(sp)}\) is the dark-field due to spectral effects, as determined in the calibration experiment. Due to the homogeneity of POM at micrometer length scales, it does not cause strong small-angle scattering—for this reason, it is assumed that all recorded signal stems from beam-hardening effects43. The dependence on the phase of the signal \(\varphi \) is neglected due to the low visibility of \(< 25\%\) of the interferometer43.
Because the porcine sample is similar in composition to a human breast, the dose in this work is given as the mean glandular dose (MGD). The MGD is calculated using the monoenergetic normalized glandular dose coefficients DgN(E) (with E the X-ray energy) as tabulated by Boone et al. for different compressed breast thicknesses44. Taking into account the spectrum of the rotating anode and summing over all energy bins E, this results in the following formula for the MGD45:
Here, K(E) is the air kerma, the total sum of kinetic energy of all charged particles per unit mass in the first step of interaction. The air kerma is measured separately for the radiograph and CT setup by placing the clinical dosimeter (Diados T60005 MAM, PTW Freiburg GmbH, Freiburg, Germany) at the sample position and leaving all gratings in the beam path. Since the holder used in the radiography experiment contains a polycarbonate plate of 5 mm in thickness in front of the sample, such a plate is also placed in front of the dosimeter. The comparatively thin walls of the FalconTM tube are deemed negligible for the dose calculation; as such, no object is placed in front of the dosimeter for the air kerma measurement in the CT setup. The resulting air kerma values are \({0.166}\,\upmu \hbox {Gys}^{-1}\) and \({0.566}\,\upmu \hbox {Gys}^{-1}\) for the radiograph and CT setup respectively. A breast glandularity of 50% is assumed for all dose calculations since approximately half of each sample consists of fat.
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
Mahnken, AH, König, AM & Figiel, JH Current technique and application of percutaneous cryotherapy.In RöFo advances in the field of X-rays and imaging methods Vol. 190 (eds Mahnken, AH et al.) 836–846 (Georg Thieme Verlag KG, 2018).
Pusceddu, C., Paliogiannis, P., Nigri, G. & Fancellu, A. Cryoablation in the management of breast cancer: Evidence to date. Breast Cancer Targets Ther. 1, 283–292 (2019).
Rewcastle, J. C., Sandison, G. A., Muldrew, K., Saliken, J. C. & Donnelly, B. J. A model for the time dependent three-dimensional thermal distribution within iceballs surrounding multiple cryoprobes. Med. Phys. 28, 1125–1137 (2001).
Article CAS PubMed Google Scholar
Niu, L.-Z., Li, J.-L. & Xu, K.-C. Percutaneous cryoablation for liver cancer. J. Clin. Transl. Hepatol. 2, 182 (2014).
PubMed PubMed Central Google Scholar
Poplack, S. P. et al. A pilot study of ultrasound-guided cryoablation of invasive ductal carcinomas up to 15 mm with mri follow-up and subsequent surgical resection. Am. J. Roentgenol. 204, 1100 (2015).
Bland, K. L., Gass, J. & Klimberg, V. S. Radiofrequency, cryoablation, and other modalities for breast cancer ablation. Surg. Clin. N. Am. 87, 539–550 (2007).
Sabel, M. S. et al. Cryoablation of early-stage breast cancer: Work-in-progress report of a multi-institutional trial. Ann. Surg. Oncol. 11, 542–549 (2004).
Pfleiderer, S. O., Marx, C., Camara, O., Gajda, M. & Kaiser, W. A. Ultrasound-guided, percutaneous cryotherapy of small (\(\le \) 15 mm) breast cancers. Investig. Radiol. 40, 472–477 (2005).
Sabel, M. S. Cryoablation as a replacement for surgical resection in early stage breast cancer. Curr. Breast Cancer Rep. 3, 109–116 (2011).
Fornage, B. D. & Hwang, R. F. Current status of imaging-guided percutaneous ablation of breast cancer. Am. J. Roentgenol. 203, 442–448 (2014).
Lanza, E. et al. Percutaneous image-guided cryoablation of breast cancer: A systematic review. J. Vasc. Interv. Radiol. 26, 1652–1657 (2015).
Pfleiderer, S. O. et al. Cryotherapy of breast cancer under ultrasound guidance: Initial results and limitations. Eur. Radiol. 12, 3009–3014 (2002).
Article MathSciNet PubMed Google Scholar
Arfelli, F. et al. Mammography with synchrotron radiation: Phase-detection techniques. Radiology 215, 286–293 (2000).
Article CAS PubMed Google Scholar
Dreossi, D. et al. The mammography project at the syrmep beamline. Eur. J. Radiol. 68, S58–S62 (2008).
Article CAS PubMed Google Scholar
Longo, R. et al. Clinical study in phase-contrast mammography: Image-quality analysis. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 372, 20130025 (2014).
David, C., Nöhammer, B., Solak, H. & Ziegler, E. Differential X-ray phase contrast imaging using a shearing interferometer. Appl. Phys. Lett. 81, 3287–3289 (2002).
Article ADS CAS Google Scholar
Momose, A. et al. Demonstration of X-ray talbot interferometry. Jpn. J. Appl. Phys. 42, L866 (2003).
Article ADS CAS Google Scholar
Momose, A. Recent advances in X-ray phase imaging. Jpn. J. Appl. Phys. 44, 6355 (2005).
Article ADS CAS Google Scholar
Pfeiffer, F., Weitkamp, T., Bunk, O. & David, C. Phase retrieval and differential phase-contrast imaging with low-brilliance X-ray sources. Nat. Phys. 2, 258–261 (2006).
Modregger, P. et al. Small angle X-ray scattering with edge-illumination. Sci. Rep. 6, 30940 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Endrizzi, M., Astolfo, A., Vittoria, F. A., Millard, T. P. & Olivo, A. Asymmetric masks for laboratory-based X-ray phase-contrast imaging with edge illumination. Sci. Rep. 6, 25466 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Pfeiffer, F. et al. Hard-X-ray dark-field imaging using a grating interferometer. Nat. Mater. 7, 134–137 (2008).
Article ADS CAS PubMed Google Scholar
Stampanoni, M. et al. The first analysis and clinical evaluation of native breast tissue using differential phase-contrast mammography. Investig. Radiol. 46, 801–806 (2011).
Grandl, S. et al. Improved visualization of breast cancer features in multifocal carcinoma using phase-contrast and dark-field mammography: An ex vivo study. Eur. Radiol. 25, 3659–3668 (2015).
Article PubMed PubMed Central Google Scholar
Scherer, K. et al. Improved diagnostics by assessing the micromorphology of breast calcifications via X-ray dark-field radiography. Sci. Rep. 6, 36991 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Arboleda, C. et al. Towards clinical grating-interferometry mammography. Eur. Radiol. 30, 1419–1425 (2020).
Heck, L. & Herzen, J. Recent advances in X-ray imaging of breast tissue: From two- to three-dimensional imaging. Phys. Med. 79, 69–79 (2020).
Article CAS PubMed Google Scholar
Nielsen, M. S., Christensen, L. B. & Feidenhans, R. Frozen and defrosted fruit revealed with X-ray dark-field radiography. Food Control 39, 222–226 (2014).
Cazzato, R. L. et al. Percutaneous image-guided cryoablation: Current applications and results in the oncologic field. Med. Oncol. 33, 1–16 (2016).
Willer, K. et al. X-ray dark-field imaging of the human lung—A feasibility study on a deceased body. PLoS ONE 13, e0204565 (2018).
Article PubMed PubMed Central Google Scholar
Frates, M. C., Benson, C. B., Doubilet, P. M., Cibas, E. S. & Marqusee, E. Can color doppler sonography aid in the prediction of malignancy of thyroid nodules? J. Ultrasound Med. 22, 127–131 (2003).
El-Brawany, M. et al. Measurement of thermal and ultrasonic properties of some biological tissues. J. Med. Eng. Technol. 33, 249–256 (2009).
Article CAS PubMed Google Scholar
Ji, X., Ge, Y., Zhang, R., Li, K. & Chen, G.-H. Studies of signal estimation bias in grating-based X-ray multicontrast imaging. Med. Phys. 44, 2453–2465 (2017).
Article CAS PubMed PubMed Central Google Scholar
Hendriksen, A. A., Pelt, D. M. & Batenburg, K. J. Noise2inverse: Self-supervised deep convolutional denoising for tomography. IEEE Trans. Comput. Imaging 6, 1320–1335 (2020).
Viermetz, M. et al. Dark-field computed tomography reaches the human scale. Proc. Natl. Acad. Sci. 119, e2118799119 (2022).
Article CAS PubMed PubMed Central Google Scholar
Kalender, W. A. et al. Technical feasibility proof for high-resolution low-dose photon-counting ct of the breast. Eur. Radiol. 27, 1081–1086 (2017).
Chevalier, M. et al. Patient dose in digital mammography. Med. Phys. 31, 2471–2479 (2004).
Zanette, I., Bech, M., Pfeiffer, F. & Weitkamp, T. Interlaced phase stepping in phase-contrast X-ray tomography. Appl. Phys. Lett. 98, 094101 (2011).
Seibert, J. A., Boone, J. M. & Lindfors, K. K. Flat-field correction technique for digital detectors. In Medical Imaging 1998: Physics of Medical Imaging Vol. 3336 (eds Seibert, J. A. et al.) 348–354 (SPIE, 1998).
Haeusele, J. et al. Advanced phase-retrieval for stepping-free X-ray dark-field computed tomography. IEEE Trans. Med. Imaging 42, 2876–2885. https://doi.org/10.1109/TMI.2023.3271413 (2023).
Allner, S. et al. Bilateral filtering using the full noise covariance matrix applied to X-ray phase-contrast computed tomography. Phys. Med. Biol. 61, 3867 (2016).
Article CAS PubMed Google Scholar
Pelzer, G. et al. A beam hardening and dispersion correction for X-ray dark-field radiography. Med. Phys. 43, 2774–2779 (2016).
De Marco, F. Image Reconstruction, Pre-clinical Studies, and Signal Formation Investigations at a Dark-Field Chest Radiography Setup. Ph.D. thesis, Technical University Munich (2021).
Boone, J. Normalized glandular dose (dgn) coefficients for arbitrary X-ray spectra in mammography: Computer-fit values of Monte Carlo derived data. Med. Phys. 29, 869–875 (2002).
Eggl, E. et al. Dose-compatible grating-based phase-contrast mammography on mastectomy specimens using a compact synchrotron source. Sci. Rep. 8, 15700 (2018).
Article ADS PubMed PubMed Central Google Scholar
The authors thank Andreas H. Mahnken for helpful input on the state of the art of cryoablation procedures. They would also like to thank Nico Huber for discussions on the development of the freezing device and Lara Szeimies for her help with the fabrication of the device.
Open Access funding enabled and organized by Projekt DEAL.
Research Group Biomedical Imaging Physics, Department of Physics, TUM School of Natural Sciences, Technical University of Munich, 85748, Garching, Germany
Dominik John, Wolfgang Gottwald, Daniel Berthe, Sami Wirtensohn, Julia Hickler, Lisa Heck & Julia Herzen
Chair of Biomedical Physics, Department of Physics, TUM School of Natural Sciences, Technical University of Munich, 85748, Garching, Germany
Dominik John, Wolfgang Gottwald, Daniel Berthe, Sami Wirtensohn, Julia Hickler, Lisa Heck & Julia Herzen
Munich Institute of Biomedical Engineering, Technical University of Munich, 85748, Garching, Germany
Dominik John, Wolfgang Gottwald, Daniel Berthe, Sami Wirtensohn, Julia Hickler, Lisa Heck & Julia Herzen
Institute of Materials Physics, Helmholtz-Zentrum hereon, 21502, Geesthacht, Germany
Dominik John & Sami Wirtensohn
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
W.G., D.J. and J. He. devised the experiments. D.J. conducted the CT measurements and processed the data. D.J., D.B. and J. Hi. developed the CT measurement acquisition protocol. W.G. and D.B. performed and processed the radiograph measurements. S.W. developed the scatter plot visualization. D.J., J. He. and L.H. wrote the manuscript with input from all authors. All authors reviewed the manuscript.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
John, D., Gottwald, W., Berthe, D. et al. X-ray dark-field computed tomography for monitoring of tissue freezing. Sci Rep 14, 5599 (2024). https://doi.org/10.1038/s41598-024-56201-3
DOI: https://doi.org/10.1038/s41598-024-56201-3
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
![X-ray dark-field computed tomography for monitoring of tissue freezing | Scientific Reports](http://cdn.globalso.com/ergodi-radiation/14f207c94.png)
Gm Survey Meter Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.