4D SPECT/CT: A Hybrid Approach to Primary Hyperparathyroidism.

Our rationale was to review the imaging options for patients with primary hyperparathyroidism and to advocate for judicious use of 4-dimensional (4D) SPECT/CT to visualize diseased parathyroid glands in patients with complex medical profiles or in whom other imaging modalities fail. We review the advantages and disadvantages of traditional imaging modalities used in preoperative assessment of patients with primary hyperparathyroidism: ultrasound, SPECT, and 4D CT. We describe a scheme for optimizing and individualizing preoperative imaging of patients with hyperfunctioning parathyroid glands using traditional modalities in tandem with 4D SPECT/CT. Using the input from radiologists, endocrinologists, and surgeons, we apply patient criteria such as large body habitus, concomitant multiglandular disease, multinodular thyroid disease, confusing previous imaging, and unsuccessful previous surgery to create an imaging paradigm that uses 4D SPECT/CT yet is cost-effective, accurate, and limits extraneous radiation exposure. 4D SPECT/CT capitalizes on the strengths of SPECT and 4D CT and addresses limitations that exist when these modalities are used in isolation. In select patients with complicated clinical parameters, preoperative imaging with 4D SPECT/CT can improve accuracy yet remain cost-effective.

Establishing a 4D-CT lung function related volumetric dose model to reduce radiation pneumonia.

In order to study how to use pulmonary functional imaging obtained through 4D-CT fusion for radiotherapy planning, and transform traditional dose volume parameters into functional dose volume parameters, a functional dose volume parameter model that may reduce level 2 and above radiation pneumonia was obtained. 41 pulmonary tumor patients who underwent 4D-CT in our department from 2020 to 2023 were included. MIM Software (MIM 7.0.7; MIM Software Inc., Cleveland, OH, USA) was used to register adjacent phase CT images in the 4D-CT series. The three-dimensional displacement vector of CT pixels was obtained when changing from one respiratory state to another respiratory state, and this three-dimensional vector was quantitatively analyzed. Thus, a color schematic diagram reflecting the degree of changes in lung CT pixels during the breathing process, namely the distribution of ventilation function strength, is obtained. Finally, this diagram is fused with the localization CT image. Select areas with Jacobi > 1.2 as high lung function areas and outline them as fLung. Import the patient's DVH image again, fuse the lung ventilation image with the localization CT image, and obtain the volume of fLung different doses (V60, V55, V50, V45, V40, V35, V30, V25, V20, V15, V10, V5). Analyze the functional dose volume parameters related to the risk of level 2 and above radiation pneumonia using R language and create a predictive model. By using stepwise regression and optimal subset method to screen for independent variables V35, V30, V25, V20, V15, and V10, the prediction formula was obtained as follows: Risk = 0.23656-0.13784 * V35 + 0.37445 * V30-0.38317 * V25 + 0.21341 * V20-0.10209 * V15 + 0.03815 * V10. These six independent variables were analyzed using a column chart, and a calibration curve was drawn using the calibrate function. It was found that the Bias corrected line and the Apparent line were very close to the Ideal line, The consistency between the predicted value and the actual value is very good. By using the ROC function to plot the ROC curve and calculating the area under the curve: 0.8475, 95% CI 0.7237-0.9713, it can also be determined that the accuracy of the model is very high. In addition, we also used Lasso method and random forest method to filter out independent variables with different results, but the calibration curve drawn by the calibration function confirmed poor prediction performance. The function dose volume parameters V35, V30, V25, V20, V15, and V10 obtained through 4D-CT are key factors affecting radiation pneumonia. Establishing a predictive model can provide more accurate lung restriction basis for clinical radiotherapy planning.

Detection of scapholunate interosseous ligament injury using dynamic computed tomography-derived arthrokinematics: A prospective clinical trial.

Scapholunate interosseous ligament injuries are a major cause of wrist instability and can be difficult to diagnose radiographically. To improve early diagnosis of scapholunate ligament injuries, we compared injury detection between bilateral routine clinical radiographs, static CT, and dynamic four-dimensional CT (4DCT) during wrist flexion-extension and radioulnar deviation. Participants with unilateral scapholunate ligament injuries were recruited to a prospective clinical trial investigating the diagnostic utility of 4DCT imaging for ligamentous wrist injury. Twenty-one participants underwent arthroscopic surgery to confirm scapholunate ligament injury. Arthrokinematics, defined as distributions of interosseous proximities across radioscaphoid and scapholunate articular surfaces at different positions within the motion cycle, were used as CT-derived biomarkers. Preoperative radiographs, static CT, and extrema of 4DCT were compared between uninjured and injured wrists using Wilcoxon signed rank or Kolmogorov-Smirnov tests. Median interosseous proximities at the scapholunate interval were significantly greater in the injured versus the uninjured wrists at static-neutral and maximum flexion, extension, radial deviation, and ulnar deviation. Mean cumulative distribution functions at the radioscaphoid joint were not significantly different between wrists but were significantly shifted at the scapholunate interval towards increased interosseous proximities in injured versus uninjured wrists in all positions. Median and cumulative distribution scapholunate proximities from static-neutral and 4DCT-derived extrema reflect injury status.

A back propagation neural network based respiratory motion modelling method.

BACKGROUND: This study presents the development of a backpropagation neural network-based respiratory motion modelling method (BP-RMM) for precisely tracking arbitrary points within lung tissue throughout free respiration, encompassing deep inspiration and expiration phases. METHODS: Internal and external respiratory data from four-dimensional computed tomography (4DCT) are processed using various artificial intelligence algorithms. Data augmentation through polynomial interpolation is employed to enhance dataset robustness. A BP neural network is then constructed to comprehensively track lung tissue movement. RESULTS: The BP-RMM demonstrates promising accuracy. In cases from the public 4DCT dataset, the average target registration error (TRE) between authentic deep respiration phases and those forecasted by BP-RMM for 75 marked points is 1.819 mm. Notably, TRE for normal respiration phases is significantly lower, with a minimum error of 0.511 mm. CONCLUSIONS: The proposed method is validated for its high accuracy and robustness, establishing it as a promising tool for surgical navigation within the lung.

Dosimetric Effect of Target Position Accuracy on Single-Isocenter Multiple Liver Metastases SBRT.

Purpose: To evaluate the dosimetric effects of intrafraction baseline shifts combined with rotational errors on Four-dimensional computed tomography-guided stereotactic body radiotherapy for multiple liver metastases (MLMs). Methods: A total of 10 patients with MLM (2 or 3 lesions) were selected for this retrospective study. Baseline shift errors of 0.5, 1.0, and 2.0 mm; and rotational errors of 0.5°, 1°, and 1.5°, were simulated about all axes. All of the baseline shifts and rotation errors were simulated around the planned isocenter using a matrix transformation of 6° of freedom. The coverage degradation of baseline shifts and rotational errors were analyzed according to the dose to 95% of the planning target volume (D95) and the volume covered by 95% of the prescribed dose (V95), and related changes in gross tumor volume were also analyzed. Results: At the rotation error of 0.5° and the baseline offset of less than 0.5 mm, the D95 and V95 values of all targets were >95%. For rotational errors of 1.0° (combined with all baseline shift errors), 36.3% of targets had D95 and V95 values of <95%. Coverage worsened substantially when the baseline shift errors were increased to 1.0 mm. D95 and V95 values were >95% for about 77.3% of the targets. Only 11.4% of the D95 and V95 values were >95% when the baseline shift errors were increased to 2.0 mm. When the rotational error was increased to 1.5° and baseline shift errors increased to 1.0 mm, the D95 and V95 values were >95% in only 3 cases. Conclusions: The multivariate regression model analysis in this study showed that the coverage of the target decreased further with reduced target volume, increasing the baseline drift, the rotation error, and the distance to the target.

Optimized scoring of end-to-end dosimetry audits for passive motion management - A simulation study using the IROC thorax phantom.

Dosimetry audits for passive motion management require dynamically-acquired measurements in a moving phantom to be compared to statically calculated planned doses. This study aimed to characterise the relationship between planning and delivery errors, and the measured dose in the Imaging and Radiation Oncology Core (IROC) thorax phantom, to assess different audit scoring approaches. Treatment plans were created using a 4DCT scan of the IROC phantom, equipped with film and thermoluminescent dosimeters (TLDs). Plans were created on the average intensity projection from all bins. Three levels of aperture complexity were explored: dynamic conformal arcs (DCAT), low-, and high-complexity volumetric modulated arcs (VMATLo, VMATHi). Simulated-measured doses were generated by modelling motion using isocenter shifts. Various errors were introduced including incorrect setup position and target delineation. Simulated-measured film doses were scored using gamma analysis and compared within specific regions of interest (ROIs) as well as the entire film plane. Positional offsets were estimated based on isodoses on the film planes, and point doses within TLD contours were compared. Motion-induced differences between planned and simulated-measured doses were evident even without introduced errors Gamma passing rates within target-centred ROIs correlated well with error-induced dose differences, while whole film passing rates did not. Isodose-based setup position measurements demonstrated high sensitivity to errors. Simulated point doses at TLD locations yielded erratic responses to introduced errors. ROI gamma analysis demonstrated enhanced sensitivity to simulated errors compared to whole film analysis. Gamma results may be further contextualized by other metrics such as setup position or maximum gamma.

4D-Precise: Learning-based 3D motion estimation and high temporal resolution 4DCT reconstruction from treatment 2D+t X-ray projections.

BACKGROUND AND OBJECTIVE: In radiotherapy treatment planning, respiration-induced motion introduces uncertainty that, if not appropriately considered, could result in dose delivery problems. 4D cone-beam computed tomography (4D-CBCT) has been developed to provide imaging guidance by reconstructing a pseudo-motion sequence of CBCT volumes through binning projection data into breathing phases. However, it suffers from artefacts and erroneously characterizes the averaged breathing motion. Furthermore, conventional 4D-CBCT can only be generated post-hoc using the full sequence of kV projections after the treatment is complete, limiting its utility. Hence, our purpose is to develop a deep-learning motion model for estimating 3D+t CT images from treatment kV projection series. METHODS: We propose an end-to-end learning-based 3D motion modelling and 4DCT reconstruction model named 4D-Precise, abbreviated from Probabilistic reconstruction of image sequences from CBCT kV projections. The model estimates voxel-wise motion fields and simultaneously reconstructs a 3DCT volume at any arbitrary time point of the input projections by transforming a reference CT volume. Developing a Torch-DRR module, it enables end-to-end training by computing Digitally Reconstructed Radiographs (DRRs) in PyTorch. During training, DRRs with matching projection angles to the input kVs are automatically extracted from reconstructed volumes and their structural dissimilarity to inputs is penalised. We introduced a novel loss function to regulate spatio-temporal motion field variations across the CT scan, leveraging planning 4DCT for prior motion distribution estimation. RESULTS: The model is trained patient-specifically using three kV scan series, each including over 1200 angular/temporal projections, and tested on three other scan series. Imaging data from five patients are analysed here. Also, the model is validated on a simulated paired 4DCT-DRR dataset created using the Surrogate Parametrised Respiratory Motion Modelling (SuPReMo). The results demonstrate that the reconstructed volumes by 4D-Precise closely resemble the ground-truth volumes in terms of Dice, volume similarity, mean contour distance, and Hausdorff distance, whereas 4D-Precise achieves smoother deformations and fewer negative Jacobian determinants compared to SuPReMo. CONCLUSIONS: Unlike conventional 4DCT reconstruction techniques that ignore breath inter-cycle motion variations, the proposed model computes both intra-cycle and inter-cycle motions. It represents motion over an extended timeframe, covering several minutes of kV scan series.

Dosimetric Impact of Delineation and Motion Uncertainties on the Heart and Substructures in Lung Cancer Radiotherapy.

AIMS: Delineation variations and organ motion produce difficult-to-quantify uncertainties in planned radiation doses to targets and organs at risk. Similar to manual contouring, most automatic segmentation tools generate single delineations per structure; however, this does not indicate the range of clinically acceptable delineations. This study develops a method to generate a range of automatic cardiac structure segmentations, incorporating motion and delineation uncertainty, and evaluates the dosimetric impact in lung cancer. MATERIALS AND METHODS: Eighteen cardiac structures were delineated using a locally developed auto-segmentation tool. It was applied to lung cancer planning CTs for 27 curative (planned dose ≥50 Gy) cases, and delineation variations were estimated by using ten mapping-atlases to provide separate substructure segmentations. Motion-related cardiac segmentation variations were estimated by auto-contouring structures on ten respiratory phases for 9/27 cases that had 4D-planning CTs. Dose volume histograms (DVHs) incorporating these variations were generated for comparison. RESULTS: Variations in mean doses (Dmean), defined as the range in values across ten feasible auto-segmentations, were calculated for each cardiac substructure. Over the study cohort the median variations for delineation uncertainty and motion were 2.20-11.09 Gy and 0.72-4.06 Gy, respectively. As relative values, variations in Dmean were between 18.7%-65.3% and 7.8%-32.5% for delineation uncertainty and motion, respectively. Doses vary depending on the individual planned dose distribution, not simply on segmentation differences, with larger dose variations to cardiac structures lying within areas of steep dose gradient. CONCLUSION: Radiotherapy dose uncertainties from delineation variations and respiratory-related heart motion were quantified using a cardiac substructure automatic segmentation tool. This predicts the 'dose range' where doses to structures are most likely to fall, rather than single DVH curves. This enables consideration of these uncertainties in cardiotoxicity research and for future plan optimisation. The tool was designed for cardiac structures, but similar methods are potentially applicable to other OARs.

A novel approach for estimating lung tumor motion based on dynamic features in 4D-CT.

Due to the high expenses involved, 4D-CT data for certain patients may only include five respiratory phases (0%, 20%, 40%, 60%, and 80%). This limitation can affect the subsequent planning of radiotherapy due to the absence of lung tumor information for the remaining five respiratory phases (10%, 30%, 50%, 70%, and 90%). This study aims to develop an interpolation method that can automatically derive tumor boundary contours for the five omitted phases using the available 5-phase 4D-CT data. The dynamic mode decomposition (DMD) method is a data-driven and model-free technique that can extract dynamic information from high-dimensional data. It enables the reconstruction of long-term dynamic patterns using only a limited number of time snapshots. The quasi-periodic motion of a deformable lung tumor caused by respiratory motion makes it suitable for treatment using DMD. The direct application of the DMD method to analyze the respiratory motion of the tumor is impractical because the tumor is three-dimensional and spans multiple CT slices. To predict the respiratory movement of lung tumors, a method called uniform angular interval (UAI) sampling was developed to generate snapshot vectors of equal length, which are suitable for DMD analysis. The effectiveness of this approach was confirmed by applying the UAI-DMD method to the 4D-CT data of ten patients with lung cancer. The results indicate that the UAI-DMD method effectively approximates the lung tumor's deformable boundary surface and nonlinear motion trajectories. The estimated tumor centroid is within 2 mm of the manually delineated centroid, a smaller margin of error compared to the traditional BSpline interpolation method, which has a margin of 3 mm. This methodology has the potential to be extended to reconstruct the 20-phase respiratory movement of a lung tumor based on dynamic features from 10-phase 4D-CT data, thereby enabling more accurate estimation of the planned target volume (PTV).

Comparing interplay effects in scanned proton therapy of lung cancer: Free breathing with various layer and volume rescanning versus respiratory gating with different gate widths.

PURPOSE: We investigated interplay effects and treatment time (TT) in scanned proton therapy for lung cancer patients. We compared free-breathing (FB) approaches with multiple rescanning strategies and respiratory-gating (RG) methods with various gating widths to identify the superior irradiation technique. METHODS: Plans were created with 4/1, 2/2, and 1/4 layered/volume rescans of FB (L4V1, L2V2, and L1V4), and 50%, 30%, and 10% gating widths of the total respiratory curves (G50, G30, and G10) of the RG plans with L4V1. We calculated 4-dimensional dynamic doses assuming a constant sinusoidal curve for six irradiation methods. The reconstructed doses per fraction were compared with planned doses in terms of dose differences in 99% clinical-target-volume (CTV) (ΔD99%), near-maximum dose differences (ΔD2%) at organs-at-risk (OARs), and TT. RESULTS: The mean/minimum CTV ΔD99% values for FB were -1.0%/-4.9%, -0.8%/-4.3%, and -0.1%/-1.0% for L4V1, L2V2, and L1V4, respectively. Those for RG were -0.3%/-1.7%, -0.1%/-1.0%, and 0.0%/-0.5% for G50, G30, and G10, respectively. The CTV ΔD99% of the RGs with less than 50% gate width and the FBs of L1V4 were within the desired tolerance (±3.0%), and the OARs ΔD2% for RG were lower than those for FB. The mean TTs were 90, 326, 824, 158, 203, and 422 s for L4V1, L2V2, L1V4, G50, G30, and G10, respectively. CONCLUSIONS: FB (L4V1) is the most efficient treatment, but not necessarily the optimal choice due to interplay effects. To satisfy both TT extensions and interplay, RG with a gate width as large as possible within safety limits is desirable.

CT respiratory motion synthesis using joint supervised and adversarial learning.

Objective.Four-dimensional computed tomography (4DCT) imaging consists in reconstructing a CT acquisition into multiple phases to track internal organ and tumor motion. It is commonly used in radiotherapy treatment planning to establish planning target volumes. However, 4DCT increases protocol complexity, may not align with patient breathing during treatment, and lead to higher radiation delivery.Approach.In this study, we propose a deep synthesis method to generate pseudo respiratory CT phases from static images for motion-aware treatment planning. The model produces patient-specific deformation vector fields (DVFs) by conditioning synthesis on external patient surface-based estimation, mimicking respiratory monitoring devices. A key methodological contribution is to encourage DVF realism through supervised DVF training while using an adversarial term jointly not only on the warped image but also on the magnitude of the DVF itself. This way, we avoid excessive smoothness typically obtained through deep unsupervised learning, and encourage correlations with the respiratory amplitude.Main results.Performance is evaluated using real 4DCT acquisitions with smaller tumor volumes than previously reported. Results demonstrate for the first time that the generated pseudo-respiratory CT phases can capture organ and tumor motion with similar accuracy to repeated 4DCT scans of the same patient. Mean inter-scans tumor center-of-mass distances and Dice similarity coefficients were 1.97 mm and 0.63, respectively, for real 4DCT phases and 2.35 mm and 0.71 for synthetic phases, and compares favorably to a state-of-the-art technique (RMSim).Significance.This study presents a deep image synthesis method that addresses the limitations of conventional 4DCT by generating pseudo-respiratory CT phases from static images. Although further studies are needed to assess the dosimetric impact of the proposed method, this approach has the potential to reduce radiation exposure in radiotherapy treatment planning while maintaining accurate motion representation. Our training and testing code can be found athttps://github.com/cyiheng/Dynagan.

The impact of motion on onboard MRI-guided pencil beam scanned proton therapy treatments.

Objective.Online magnetic resonance imaging (MRI) guidance could be especially beneficial for pencil beam scanned (PBS) proton therapy of tumours affected by respiratory motion. For the first time to our knowledge, we investigate the dosimetric impact of respiratory motion on MRI-guided proton therapy compared to the scenario without magnetic field.Approach.A previously developed analytical proton dose calculation algorithm accounting for perpendicular magnetic fields was extended to enable 4D dose calculations. For two geometrical phantoms and three liver and two lung patient cases, static treatment plans were optimised with and without magnetic field (0, 0.5 and 1.5 T). Furthermore, plans were optimised using gantry angle corrections (0.5 T +5° and 1.5 T +15°) to reproduce similar beam trajectories compared to the 0 T reference plans. The effect of motion was then considered using 4D dose calculations without any motion mitigation and simulating 8-times volumetric rescanning, with motion for the patient cases provided by 4DCT(MRI) data sets. Each 4D dose calculation was performed for different starting phases and the CTV dose coverageV95%and homogeneityD5%-D95%were analysed.Main results.For the geometrical phantoms with rigid motion perpendicular to the beam and parallel to the magnetic field, a comparable dosimetric effect was observed independent of the magnetic field. Also for the five 4DCT(MRI) cases, the influence of motion was comparable for all magnetic field strengths with and without gantry angle correction. On average, the motion-induced decrease in CTVV95%from the static plan was 17.0% and 18.9% for 1.5 T and 0.5 T, respectively, and 19.9% without magnetic field.Significance.For the first time, this study investigates the combined impact of magnetic fields and respiratory motion on MR-guided proton therapy. The comparable dosimetric effects irrespective of magnetic field strength indicate that the effects of motion for future MR-guided proton therapy may not be worse than for conventional PBS proton therapy.

Morphological and dynamic analysis of the normal aortic valve with 4D computed tomography.

OBJECTIVES: To evaluate the precise dimensions of the normal aortic root, especially the true aortic annulus, during the cardiac cycle using an innovative reconstruction method based on multiphase cardiac computed tomography and to assess the feasibility and the reproducibility of this method for aortic root analysis. METHODS: Between January 2019 and June 2021, 30 optimal consecutive ECG-gated multiphase cardiac computed tomography of patients with normal tricuspid aortic valve were analysed using an in-house software. Aortic annulus border was pinpointed on 9 reconstructed planes and the 3D coordinates of the 18 consecutive points were interpolated into a 3D curve using a cubic spline. Three additional planes were generated at the level of the left ventricular outflow tract, the level of the Valsalva sinus and the level of the sinotubular junction. This procedure was repeated for all the 10 temporal phases of the RR interval. RESULTS: The aortic annulus mean 3D and 2D areas were 7.67 ± 1.51 and 5.16 ± 1.40 cm2, respectively. The mean 2D diameter was 2.51 ± 0.23 cm. The mean global area expansion was 11.8 ± 3.5% and the mean perimeter expansion of 7.1 ± 2.6%. During the cardiac cycle, the left ventricle outflow tract expands, reaching its maximum surface at the end of diastole, followed by the aortic annulus, the Valsalva sinuses and the sinotubular junction. The aorta changes from a clover-shaped cone during diastole to more cylindrical shape during systole. Compared to the 3D measurements, the analysis of the virtual basal ring significantly underestimates the annulus area, perimeter, and mean diameter. CONCLUSIONS: 4D morphometric analysis enables to have a precise and reproducible evaluation of the aortic annulus. The aortic annulus and root are deformable structures that undergo a unique expansion sequence during the cardiac cycle which should be considered for procedural planning.

Motion compensated cone-beam CT reconstruction using ana priorimotion model from CT simulation: a pilot study.

Objective. To combat the motion artifacts present in traditional 4D-CBCT reconstruction, an iterative technique known as the motion-compensated simultaneous algebraic reconstruction technique (MC-SART) was previously developed. MC-SART employs a 4D-CBCT reconstruction to obtain an initial model, which suffers from a lack of sufficient projections in each bin. The purpose of this study is to demonstrate the feasibility of introducing a motion model acquired during CT simulation to MC-SART, coined model-based CBCT (MB-CBCT).Approach. For each of 5 patients, we acquired 5DCTs during simulation and pre-treatment CBCTs with a simultaneous breathing surrogate. We cross-calibrated the 5DCT and CBCT breathing waveforms by matching the diaphragms and employed the 5DCT motion model parameters for MC-SART. We introduced the Amplitude Reassignment Motion Modeling technique, which measures the ability of the model to control diaphragm sharpness by reassigning projection amplitudes with varying resolution. We evaluated the sharpness of tumors and compared them between MB-CBCT and 4D-CBCT. We quantified sharpness by fitting an error function across anatomical boundaries. Furthermore, we compared our MB-CBCT approach to the traditional MC-SART approach. We evaluated MB-CBCT's robustness over time by reconstructing multiple fractions for each patient and measuring consistency in tumor centroid locations between 4D-CBCT and MB-CBCT.Main results. We found that the diaphragm sharpness rose consistently with increasing amplitude resolution for 4/5 patients. We observed consistently high image quality across multiple fractions, and observed stable tumor centroids with an average 0.74 ± 0.31 mm difference between the 4D-CBCT and MB-CBCT. Overall, vast improvements over 3D-CBCT and 4D-CBCT were demonstrated by our MB-CBCT technique in terms of both diaphragm sharpness and overall image quality.Significance. This work is an important extension of the MC-SART technique. We demonstrated the ability ofa priori5DCT models to provide motion compensation for CBCT reconstruction. We showed improvements in image quality over both 4D-CBCT and the traditional MC-SART approach.

PNMC: Four-dimensional conebeam CT reconstruction combining prior network and motion compensation.

Four-dimensional conebeam computed tomography (4D CBCT) is an efficient technique to overcome motion artifacts caused by organ motion during breathing. 4D CBCT reconstruction in a single scan usually divides projections into different groups of sparsely sampled data based on the respiratory phases. The reconstructed images within each group present poor image quality due to the limited number of projections. To improve the image quality of 4D CBCT in a single scan, we propose a novel reconstruction scheme that combines prior knowledge with motion compensation. We apply the reconstructed images of the full projections within a single routine as prior knowledge, providing structural information for the network to enhance the restoration structure. The prior network (PN-Net) is proposed to extract features of prior knowledge and fuse them with the sparsely sampled data using an attention mechanism. The prior knowledge guides the reconstruction process to restore the approximate organ structure and alleviates severe streaking artifacts. The deformation vector field (DVF) extracted using deformable image registration among different phases is then applied in the motion-compensated ordered-subset simultaneous algebraic reconstruction algorithm to generate 4D CBCT images. Proposed method has been evaluated using simulated and clinical datasets and has shown promising results by comparative experiment. Compared with previous methods, our approach exhibits significant improvements across various evaluation metrics.

Imaging Characteristics and Interpretation Strategy for Renal Hyperparathyroidism: A Retrospective 4-Dimensional Computed Tomography Study.

OBJECTIVE: Parathyroidectomy treats uncontrolled renal hyperparathyroidism (RHPT), requiring identification of all glands. Three types of enhancement are proposed. Type A lesions have higher arterial phase attenuation than the thyroid, type B lesions lack higher arterial phase attenuation but have lower venous phase attenuation, and type C lesions have neither higher arterial phase attenuation nor lower venous phase attenuation than the thyroid. We aimed to outline the image features of problematic parathyroid glands in RHPT and propose a 4-dimensional computed tomography (4DCT) interpretation algorithm. METHODS: This retrospective study involved data collection from patients with RHPT who underwent preoperative 4DCT for parathyroidectomy between January and November 2022. Pathologically confirmed parathyroid lesions were retrospectively identified on 4DCT according to the location and size described in the surgical notes. The attenuation of parathyroid lesions and the thyroid glands was assessed in 3 phases, and demographic data of the patients were collected. RESULTS: Ninety-seven pathology-proven parathyroid glands from 27 patients were obtained, with 86 retrospectively detected on 4DCT. In the arterial phase, the attenuation of parathyroid lesions in RHPT did not exceed that of the thyroid gland (P < .001). In the venous phase, parathyroid lesions demonstrated lower attenuation than the thyroid gland (P < .001). A total of 81 parathyroid lesions (94.2%) exhibited type B patterns. CONCLUSION: Unlike primary hyperparathyroidism, lesions in RHPT exhibited more type B enhancement, making them less readily identifiable in the arterial phase. Therefore, we propose a distinct imaging interpretation strategy to locate these problematic glands more efficiently.

A model that predicts a real-time tumour surface using intra-treatment skin surface and end-of-expiration and end-of-inhalation planning CT images.

OBJECTIVES: To develop a mapping model between skin surface motion and internal tumour motion and deformation using end-of-exhalation (EOE) and end-of-inhalation (EOI) 3D CT images for tracking lung tumours during respiration. METHODS: Before treatment, skin and tumour surfaces were segmented and reconstructed from the EOE and the EOI 3D CT images. A non-rigid registration algorithm was used to register the EOE skin and tumour surfaces to the EOI, resulting in a displacement vector field that was then used to construct a mapping model. During treatment, the EOE skin surface was registered to the real-time, yielding a real-time skin surface displacement vector field. Using the mapping model generated, the input of a real-time skin surface can be used to calculate the real-time tumour surface. The proposed method was validated with and without simulated noise on 4D CT images from 15 patients at Léon Bérard Cancer Center and the 4D-lung dataset. RESULTS: The average centre position error, dice similarity coefficient (DSC), 95%-Hausdorff distance and mean distance to agreement of the tumour surfaces were 1.29 mm, 0.924, 2.76 mm, and 1.13 mm without simulated noise, respectively. With simulated noise, these values were 1.33 mm, 0.920, 2.79 mm, and 1.15 mm, respectively. CONCLUSIONS: A patient-specific model was proposed and validated that was constructed using only EOE and EOI 3D CT images and real-time skin surface images to predict internal tumour motion and deformation during respiratory motion. ADVANCES IN KNOWLEDGE: The proposed method achieves comparable accuracy to state-of-the-art methods with fewer pre-treatment planning CT images, which holds potential for application in precise image-guided radiation therapy.

Benchmarking machine learning-based real-time respiratory signal predictors in 4D SBRT.

BACKGROUND: Stereotactic body radiotherapy of thoracic and abdominal tumors has to account for respiratory intrafractional tumor motion. Commonly, an external breathing signal is continuously acquired that serves as a surrogate of the tumor motion and forms the basis of strategies like breathing-guided imaging and gated dose delivery. However, due to inherent system latencies, there exists a temporal lag between the acquired respiratory signal and the system response. Respiratory signal prediction models aim to compensate for the time delays and to improve imaging and dose delivery. PURPOSE: The present study explores and compares six state-of-the-art machine and deep learning-based prediction models, focusing on real-time and real-world applicability. All models and data are provided as open source and data to ensure reproducibility of the results and foster reuse. METHODS: The study was based on 2502 breathing signals ( t t o t a l ≈ 90 $t_{total} \approx 90$  h) acquired during clinical routine, split into independent training (50%), validation (20%), and test sets (30%). Input signal values were sampled from noisy signals, and the target signal values were selected from corresponding denoised signals. A standard linear prediction model (Linear), two state-of-the-art models in general univariate signal prediction (Dlinear, Xgboost), and three deep learning models (Lstm, Trans-Enc, Trans-TSF) were chosen. The prediction performance was evaluated for three different prediction horizons (480, 680, and 920 ms). Moreover, the robustness of the different models when applied to atypical, that is, out-of-distribution (OOD) signals, was analyzed. RESULTS: The Lstm model achieved the lowest normalized root mean square error for all prediction horizons. The prediction errors only slightly increased for longer horizons. However, a substantial spread of the error values across the test signals was observed. Compared to typical, that is, in-distribution test signals, the prediction accuracy of all models decreased when applied to OOD signals. The more complex deep learning models Lstm and Trans-Enc showed the least performance loss, while the performance of simpler models like Linear dropped the most. Except for Trans-Enc, inference times for the different models allowed for real-time application. CONCLUSION: The application of the Lstm model achieved the lowest prediction errors. Simpler prediction filters suffer from limited signal history access, resulting in a drop in performance for OOD signals.

Annotation-free prediction of treatment-specific tissue outcome from 4D CT perfusion imaging in acute ischemic stroke.

Acute ischemic stroke is a critical health condition that requires timely intervention. Following admission, clinicians typically use perfusion imaging to facilitate treatment decision-making. While deep learning models leveraging perfusion data have demonstrated the ability to predict post-treatment tissue infarction for individual patients, predictions are often represented as binary or probabilistic masks that are not straightforward to interpret or easy to obtain. Moreover, these models typically rely on large amounts of subjectively segmented data and non-standard perfusion analysis techniques. To address these challenges, we propose a novel deep learning approach that directly predicts follow-up computed tomography images from full spatio-temporal 4D perfusion scans through a temporal compression. The results show that this method leads to realistic follow-up image predictions containing the infarcted tissue outcomes. The proposed compression method achieves comparable prediction results to using perfusion maps as inputs but without the need for perfusion analysis or arterial input function selection. Additionally, separate models trained on 45 patients treated with thrombolysis and 102 treated with thrombectomy showed that each model correctly captured the different patient-specific treatment effects as shown by image difference maps. The findings of this work clearly highlight the potential of our method to provide interpretable stroke treatment decision support without requiring manual annotations.

Myocardial analysis from routine 4D cardiac-CT to predict reverse remodeling and clinical outcomes after transcatheter aortic valve implantation.

PURPOSE: Our study aimed to determine whether 4D cardiac computed tomography (4DCCT) based quantitative myocardial analysis may improve risk stratification and can predict reverse remodeling (RRM) and mortality after transcatheter aortic valve implantation (TAVI). METHODS: Consecutive patients undergoing clinically indicated 4DCCT prior to TAVI were prospectively enrolled. 4DCCT-derived left- (LV) and right ventricular (RV), and left atrial (LA) dimensions, mass, ejection fraction (EF) and myocardial strain were evaluated to predict RRM and survival. RRM was defined by either relative increase in LVEF by 5% or relative decline in LV end diastolic diameter (LVEDD) by 5% assessed by transthoracic echocardiography prior TAVI, at discharge, and at 12-month follow-up compared to baseline prior to TAVI. RESULTS: Among 608 patients included in this study (55 % males, age 81 ± 6.6 years), RRM was observed in 279 (54 %) of 519 patients at discharge and in 218 (48 %) of 453 patients at 12-month echocardiography. While no CCT based measurements predicted RRM at discharge, CCT based LV mass index and LVEF independently predicted RRM at 12-month (ORadj = 1.012; 95 %CI:1.001-1.024; p = 0.046 and ORadj = 0.969; 95 %CI:0.943-0.996; p = 0.024, respectively). The most pronounced changes in LVEF and LVEDD were observed in patients with impaired LV function at baseline. In multivariable analysis age (HRadj = 1.037; 95 %CI:1.005-1.070; p = 0.022) and CCT-based LVEF (HRadj = 0.972; 95 %CI:0.945-0.999; p = 0.048) and LAEF (HRadj = 0.982; 95 %CI:0.968-0.996; p = 0.011) independently predicted survival. CONCLUSION: Comprehensive myocardial functional information derived from routine 4DCCT in patients with severe aortic stenosis undergoing TAVI could predict reverse remodeling and clinical outcomes at 12-month following TAVI.