Physics in medicine and biology最新文献

Objective: T2-weighted 2D fast spin echo sequence serves as the standard sequence in clinical pelvic MR imaging protocols. However, motion artifacts and blurring caused by peristalsis present significant challenges. Patient preparation such as administering antiperistaltic agents is often required before examination to reduce artifacts, which discomfort the patients. This work introduce a novel dynamic approach for T2 weighted pelvic imaging to address peristalsis-induced motion issue without any patient preparation. Approach: A rapid dynamic data acquisition strategy with complementary sampling trajectory is designed to enable highly undersampled motion-resistant data sampling, and an unrolling method based on deep equilibrium model is leveraged to reconstruct images from the dynamic sampled k-space data. Moreover, the fix-point convergence of the equilibrium model ensures the stability of the reconstruction. The high acceleration factor in each temporal phase, which is much higher than that in traditional static imaging, has the potential to effectively freeze pelvic motion, thereby transforming the imaging problem from conventional motion prevention or removal to motion reconstruction. Main results: Experiments on both retrospective and prospective data have demonstrated the superior performance of the proposed dynamic approach in reducing motion artifacts and accurately depicting structural details compared to standard static imaging. Significance: The proposed dynamic approach effectively captures motion states through dynamic data acquisition and deep learning-based reconstruction, addressing motion-related challenges in pelvic imaging.

Objective: Blood flow sensitivity is a crucial metric for appraising the effectiveness of color Doppler flow imaging (CDFI). Color Doppler velocity maps based on classic autocorrelation techniques are widely used in clinical practice. However, these techniques often produce twinkling artifacts in noisy regions due to the inherent randomness of noise phases. To mitigate artifacts and improve image quality, Power Mask (PoM) technology becomes imperative. Nevertheless, PoM technology unintentionally filters out small flow signals that have similar power and frequency characteristics to noise signals, thereby reducing the imaging system's sensitivity to flow. Approach: To address this issue, a novel Flow Recycling Algorithm (FRA) based on phase anomaly is introduced in this study. This algorithm, excavating small flow signals from noise, aims to enhance the small flow signals with low-velocity by the phase characteristics of the color Doppler flow information. Main results: Experiments in multi-organ imaging have shown that the FRA-CDFI approach is more effective in suppressing twinkling artifacts in noisy regions, preserving intricate small flow signals, and markedly improving small blood flow sensitivity. This novel approach provides adequate technical support for clinical ultrasound imaging of organs with dense small blood vessels, such as the brain, kidneys, liver, and more. Significance: As a novel post-processing method, FRA-CDFI holds significant potential for future deployment in clinical high-frame-rate ultrasound imaging devices.

Objective.To develop a 2D MR acceleration method utilizing principal component analysis (PCA) in a hybrid fashion for rapid real-time applications.Approach.Retrospective testing was performed on 10 lung, 10 liver and 10 prostate 3T MRI data sets for image quality and target contourability. Sampling of k-space is performed by acquiring central (low-frequency) data in every frame while the high-frequency data is incoherently undersampled such that all of k-space is acquired in a pre-determined number of frames. Firstly, principal components (PCs) representative of intra-frame correlations between central and outer k-space data are used to estimate unsampled data in the frame of interest. Then to add further stability, PCs representative of time-domain fluctuations within a reconstruction window of the most recent frames are fit to outer k-space data (including above estimations) to obtain final estimates in the frame of interest. Accelerated reconstructions between 3x and 8x were tested for image quality and contourability along with the optimal number of PCs for fitting.Main results.It was found that at higher acceleration rates, image quality did not deteriorate significantly. Similarly, it was found that the images were of sufficient quality to contour a target using auto-contouring software at all tested acceleration rates and sites. SSIM values were found to be ⩾0.91 at all accelerations tested. Similarly dice coefficients at the different sites were found to be ⩾0.89 even at 8x accelerations which is on par with or better than intra-observer variation.Significance.This method appears to produce improved image quality and contourability compared to previous PCA methods while also allowing a greater number of PCs to be used in reconstruction. The method can be run using a simple single-channel coil and does not require significant computing power to meet real-time interventional standards (reconstruction times ∼60 ms/frame on Intel i5 CPU).

We read with great interest the paper by Limet al(2018Phys. Med. Biol.63035042) on bias reduction in Y-90 PET imaging. In particular, they proposed a new formulation of the tomographic reconstruction problem that enforces non-negativity in projection space as opposed to image space. We comment on the algorithm they derived from this formulation and bring some clarifications on the constraint that this algorithm respects.

Objective. Real-time adaptive particle therapy is being investigated as a means to maximize the treatment delivery accuracy. To react to dosimetric errors, a system for fast and reliable verification of the agreement between planned and delivered doses is essential. This study presents a clinically feasible, real-time 4D-dose reconstruction system, synchronized with the treatment delivery and motion of the patient, which can provide the necessary feedback on the quality of the delivery.Approach. A GPU-based analytical dose engine capable of millisecond dose calculation for carbon ion therapy has been developed and interfaced with the next generation of the dose delivery system (DDS) in use at Centro Nazionale di Adroterapia Oncologica (CNAO). The system receives the spot parameters and the motion information of the patient during the treatment and performs the reconstruction of the planned and delivered 4D-doses. After each iso-energy layer, the results are displayed on a graphical user interface by the end of the spill pause of the synchrotron, permitting verification against the reference dose. The framework has been verified experimentally at CNAO for a lung cancer case based on a virtual phantom 4DCT. The patient's motion was mimicked by a moving Ionization Chamber (IC) 2D-array.Mainresults. For the investigated static and 4D-optimized treatment delivery cases, real-time dose reconstruction was achieved with an average pencil beam dose calculation speed up to more than one order of magnitude smaller than the spot delivery. The reconstructed doses have been benchmarked against offline log-file based dose reconstruction with the TRiP98 treatment planning system, as well as QA measurements with the IC 2D-array, where an average gamma-index passing rate (3%/3 mm) of 99.8% and 98.3%, respectively, were achieved.Significance. This work provides the first real-time 4D-dose reconstruction engine for carbon ion therapy. The framework integration with the CNAO DDS paves the way for a swift transition to the clinics.

Objective.¹²N, having a half-life of 11 ms, is a highly effective positron emitter that can potentially provide near real-time feedback in proton therapy. There is currently no framework for comparing and validating positron emission imaging of¹²N. This work describes the development and validation of a Monte Carlo (MC) framework to calculate the images of¹²N, as well as long-lived isotopes, originating from activation by protons.Approach. The available dual-panel Biograph mCT PET scanner was modeled in GATE and validated by comparing the simulated sensitivity map with the measured one. The distributions of¹²N and long-lived isotopes were calculated by RayStation and used as the input of GATE simulations. The RayStation/GATE combination was verified using proton beam irradiations of homogeneous phantoms. A 120 MeV pulsed pencil beam with 10⁸protons per pulse was used. Two-dimensional images were created from the GATE output and compared with the images based on the measurements and the 1D longitudinal projection of the full 2D image was used to calculate the¹²N activity range.Main results. The simulated sensitivity in the center of the FoV (5.44%) agrees well with the measured one (5.41%). The simulated and measured 2D sensitivity maps agree in good detail. The relative difference between the measured and simulated positron activity range for both¹²N and long-lived isotopes is less than 1%. The broadening of the¹²N images relative to those of the longer-lived isotopes can be understood in terms of the large positron range of¹²N.Significance. We developed and validated a MC framework based on RayStation/GATE to support the in-beam PET method for quality assurance of proton therapy. The inclusion of the very short-lived isotope¹²N makes the framework useful for developing near real-time verification. This represents a significant step towards translating¹²N real-time in vivo verification to the clinic.

Objective.Photon-counting x-ray detectors (PCDs) can produce dual-energy (DE) x-ray images of lung cancer in a single x-ray exposure. It is important to understand the factors that affect contrast, noise and the contrast-to-noise ratio (CNR). This study quantifies the dependence of CNR on tube voltage, energy threshold and patient thickness in single exposure, DE, bone-suppressed thoracic imaging with PCDs, and elucidates how the fundamental processes inherent in x-ray detection by PCDs contribute to CNR degradation.Approach.We modeled the DE CNR for five theoretical PCDs, ranging from an ideal PCD that detects every primary photon in the correct energy bin while rejecting all scattered radiation to a non-ideal PCD that suffers from charge-sharing and electronic noise, and detects scatter. CNR was computed as a function of tube voltage and high energy threshold for average and larger-than-average patients. Model predictions were compared with experimental data extracted from images acquired using a cadmium telluride (CdTe) PCD with two energy bins and analog charge summing for charge-sharing suppression. The imaging phantom simulated attenuation, scatter and contrast in lung nodule imaging. We quantified CNR improvements achievable with anti-correlated noise reduction (ACNR) and measured the range of exposure rates over which pulse pile-up is negligible.Main Results.The realistic model predicted overall trends observed in the experimental data. CNR improvements with ACNR were approximately five-fold, and modeled CNR-enhancements were on average within 10% of experiment. CNR increased modestly (i.e.<20%) when increasing the tube voltage from 90 kV to 130 kV. Optimal energy thresholds ranged from 50 keV to 70 keV across all tube voltages and patient thicknesses with and without ACNR. Quantum efficiency, electronic noise, charge sharing and scatter degraded CNR by ~50%. Charge sharing and scatter had the largest effect on CNR, degrading it by ~30% and ~15% respectively. Dead-time losses were less than 5% for patient exposure rates within the range of clinical exposure rates.Significance.In this study, we (1) employed analytical and computational models to assess the impact of different factors on CNR in single-exposure DE imaging with PCDs, (2) evaluated the accuracy of these models in predicting experimental trends, (3) quantified improvements in CNR achievable through ACNR and (4) determined the range of patient exposure rates at which pulse pile-up can be considered negligible. To the best of our knowledge, this study represents the first systematic investigation of single-exposure DE imaging of lung nodules with PCDs.

Objective: Photon-counting detectors (PCDs) for CT imaging use energy thresholds to simultaneously acquire projections at multiple energies, making them suitable for spectral imaging and material decomposition. Unfortunately, setting multiple energy thresholds results in noisy analytical reconstructions due to low photon counts in high-energy bins. Iterative reconstruction provides high quality photon-counting CT (PCCT) images but requires enormous computation time for 5D (3D + energy + time) in vivo cardiac imaging. Approach. We recently introduced UnetU, a deep learning (DL) approach that accurately denoises axial slices from 4D (3D + energy) PCCT reconstructions at various acquisition settings. In this study, we explore UnetU configurations for 5D cardiac PCCT denoising, focusing on singular value decomposition (SVD) modifications along the energy and time dimensions and alternate network architectures such as 3D U-net, FastDVDNet, and Swin Transformer UNet. We compare our networks to multi-energy non-local means (ME NLM), an established PCCT denoising algorithm. Main results. Our evaluation, using real mouse data and the digital MOBY phantom, revealed that all DL methods were more than 16 times faster than iterative reconstruction. DL denoising with SVD along the energy dimension was most effective, consistently providing low root mean square error and spatio-temporal reduced reference entropic difference, alongside strong qualitative agreement with iterative reconstruction. This superiority was attributed to lower effective rank along the energy dimension than the time dimension in 5D cardiac PCCT reconstructions. ME NLM sometimes outperformed DL with time SVD or time and energy SVD, but lagged behind iterative reconstruction and DL with energy SVD. Among alternate DL architectures with energy SVD, none consistently outperformed UnetU Energy (2D). Significance. Our study establishes UnetU Energy as an accurate and efficient method for 5D cardiac PCCT denoising, offering a 32-fold speed increase from iterative reconstruction. This advancement sets a new benchmark for DL applications in cardiovascular imaging.

Objective:The formation of functional vasculature in solid tumours enables delivery of oxygen and nutrients, and is vital for effective treatment with chemotherapeutic agents. Longitudinal characterisation of vascular networks can be enabled using mesoscopic photoacoustic imaging, but requires accurate image co-registration to precisely assess local changes across disease development or in response to therapy. Co-registration in photoacoustic imaging is challenging due to the complex nature of the generated signal, including the sparsity of data, artefacts related to the illumination/detection geometry, scan-to-scan technical variability, and biological variability, such as transient changes in perfusion. To better inform the choice of co-registration algorithms, we compared five open-source methods, in physiological and pathological tissues, with the aim of aligning evolving vascular networks in tumours imaged over growth at different time-points.Approach:Co-registration techniques were applied to 3D vascular images acquired with photoacoustic mesoscopy from murine ears and breast cancer patient-derived xenografts, at a fixed time-point and longitudinally. Images were pre-processed and segmented using an unsupervised generative adversarial network. To compare co-registration quality in different settings, pairs of fixed and moving intensity images and/or segmentations were fed into five methods split into the following categories: affine intensity-based using 1)mutual information (MI) or 2)normalised cross-correlation (NCC) as optimisation metrics, affine shape-based using 3)NCC applied to distance-transformed segmentations or 4)iterative closest point algorithm, and deformable weakly supervised deep learning-based using 5)LocalNet co-registration. Percent-changes in Dice coefficients, surface distances, MI, structural similarity index measure and target registration errors were evaluated.Main results:Co-registration using MI or NCC provided similar alignment performance, better than shape-based methods. LocalNet provided accurate co-registration of substructures by optimising subfield deformation throughout the volumes, outperforming other methods, especially in the longitudinal breast cancer xenograft dataset by minimising target registration errors.Significance:We showed the feasibility of co-registering repeatedly or longitudinally imaged vascular networks in photoacoustic mesoscopy, taking a step towards longitudinal quantitative characterisation of these complex structures. These tools open new outlooks for monitoring tumour angiogenesis at the meso-scale and for quantifying treatment-induced co-localised alterations in the vasculature.

Objective.The energy deposition of photons and protons differs. It depends on the position in the proton Bragg peak (BP) and the linear energy transfer (LET) leading to a variable relative biological effectiveness (RBE). Here, we investigate LET dependent alterations on metabolic viability and proliferation of sarcoma and endothelium cell lines following proton irradiation in comparison to photon exposure.Approach.Using a multi-step range shifter, each column of a 96-well plate was positioned in a different depth along four BP curves with increasing intensities. The high-throughput experimental setup covers dose, LET, and RBE changes seen in a treatment field. Photon irradiation was performed to calculate the RBE along the BP curve. Two biological information out of one experiment were extracted allowing a correlation between metabolic viability and proliferation of the cells.Main results.The metabolic viability and cellular proliferation were column-wise altered showing a depth-dose profile. Endothelium cell viability recovers within 96 h post BP irradiation while sarcoma cell viability remains reduced. Highest RBE values were observed at the BP distal fall-off regarding proliferation of the sarcoma and endothelial cells.Significance.The high-throughput experimental setup introduced here (I) covers dose, LET, and RBE changes seen in a treatment field, (II) measures short-term effects within 48 h to 96 h post irradiation, and (III) can additionally be transferred to various cell types without time consuming experimental adaptations. Traditionally, RBE values are calculated from clonogenic cell survival. Measured RBE profiles strongly depend on physical characteristics such as dose and LET and biological characteristics for example cell type and time point. Metabolic viability and proliferation proofed to be in a similar effect range compared to clonogenic survival results. Based on limited data of combined irradiation with doxorubicin, future experiments will test combined treatment with systemic therapies applied in clinics e.g. cyclin-dependent inhibitors.