No Arabic abstract
Statistical iterative reconstruction is expected to improve the image quality of megavoltage computed tomography (MVCT). However, one of the challenges of iterative reconstruction is its large computational cost. The purpose of this work is to develop a fast iterative reconstruction algorithm by combining several iterative techniques and by optimizing reconstruction parameters. Megavolt projection data was acquired from a TomoTherapy system and reconstructed using our statistical iterative reconstruction. Total variation was used as the regularization term and the weight of the regularization term was determined by evaluating signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), and visual assessment of spatial resolution using Gammex and Cheese phantoms. Gradient decent with an adaptive convergence parameter, ordered subset expectation maximization (OSEM), and CPU/GPU parallelization were applied in order to accelerate the present reconstruction algorithm. The SNR and CNR of the iterative reconstruction were several times better than that of filtered back projection (FBP). The GPU parallelization code combined with the OSEM algorithm reconstructed an image several hundred times faster than a CPU calculation. With 500 iterations, which provided good convergence, our method produced a 512$times$512 pixel image within a few seconds. The image quality of the present algorithm was much better than that of FBP for patient data. An image from the iterative reconstruction in TomoTherapy can be obtained within few seconds by fine-tuning the parameters. The iterative reconstruction with GPU was fast enough for clinical use, and largely improve the MVCT images.
In radial fast spin-echo MRI, a set of overlapping spokes with an inconsistent T2 weighting is acquired, which results in an averaged image contrast when employing conventional image reconstruction techniques. This work demonstrates that the problem may be overcome with the use of a dedicated reconstruction method that further allows for T2 quantification by extracting the embedded relaxation information. Thus, the proposed reconstruction method directly yields a spin-density and relaxivity map from only a single radial data set. The method is based on an inverse formulation of the problem and involves a modeling of the received MRI signal. Because the solution is found by numerical optimization, the approach exploits all data acquired. Further, it handles multi-coil data and optionally allows for the incorporation of additional prior knowledge. Simulations and experimental results for a phantom and human brain in vivo demonstrate that the method yields spin-density and relaxivity maps that are neither affected by the typical artifacts from TE mixing, nor by streaking artifacts from the incomplete k-space coverage at individual echo times.
Statistical event reconstruction techniques can give better results for gamma cameras than the traditional centroid method. However, implementation of such techniques requires detailed knowledge of the PMT light response functions. Here we describe an iterative technique which allows to obtain the response functions from flood irradiation data without imposing strict requirements on the spatial uniformity of the event distribution. A successful application of the technique for medical gamma cameras is demonstrated using both simulated and experimental data. We show that this technique can be used for monitoring of the photomultiplier gain variations. An implementation of the iterative reconstruction technique capable of operating in real-time is also presented.
Introduction: Treating pregnant women in the radiotherapy clinic is a rare occurrence. When it does occur, it is vital that the dose received by the developing embryo or foetus is understood as fully as possible. This study presents the first investigation of foetal doses delivered during helical tomotherapy treatments. Materials & Methods: Six treatment plans were delivered to an anthropomorphic phantom using a tomotherapy machine. These included treatments of the brain, unilateral and bilateral head-and-neck, chest wall, and upper lung. Measurements of foetal dose were made with an ionisation chamber positioned at various locations longitudinally within the phantom to simulate a variety of patient anatomies. Results: All measurements were below the established limit of 100 mGy for a high risk of damage during the first trimester. The largest dose encountered was 75 mGy (0.125% of prescription dose). The majority of treatments with measurement positions less than 30 cm fell into the range of uncertain risk (50 - 100 mGy). All treatments with measurement positions beyond 30 cm fell into the low risk category (< 50 mGy). Conclusions: For the cases in this study, tomotherapy resulted in foetal doses that are at least on par with, if not significantly lower than, similar 3D conformal or intensity-modulated treatments delivered with other devices. Recommendations were also provided for estimating foetal doses from tomotherapy plans.
Simulation-based image quality metrics are adapted and investigated for characterizing the parameter dependences of linear iterative image reconstruction for DBT. Three metrics based on 2D DBT simulation are investigated: (1) a root-mean-square-error (RMSE) between the test phantom and reconstructed image, (2) a gradient RMSE where the comparison is made after taking a spatial gradient of both image and phantom, and (3) a region-of-interest (ROI) Hotelling observer (HO) for signal-known-exactly/background-known-exactly (SKE/BKE) and signal-known-exactly/background-known-statistically (SKE/BKS) detection tasks. Two simulation studies are performed using the aforementioned metrics, varying voxel aspect ratio and regularization strength for two types of Tikhonov regularized least-squares optimization. The RMSE metrics are applied to a 2D test phantom and the ROI-HO metric is applied to two tasks relevant to DBT: large, low contrast lesion detection and small, high contrast microcalcification detection. The RMSE metric trends are compared with visual assessment of the reconstructed test phantom. The ROI-HO metric trends are compared with 3D reconstructed images from ACR phantom data acquired with a Hologic Selenia Dimensions DBT system. Sensitivity of image RMSE to mean pixel value is found to limit its applicability to the assessment of DBT image reconstruction. Image gradient RMSE is insensitive to mean pixel value and appears to track better with subjective visualization of the reconstructed bar-pattern phantom. The ROI-HO metric shows an increasing trend with regularization strength for both forms of Tikhonov-regularized least-squares; however, this metric saturates at intermediate regularization strength indicating a point of diminishing returns for signal detection. Visualization with reconstructed ACR phantom images appears to show a similar dependence with regularization strength.
A number of image-processing problems can be formulated as optimization problems. The objective function typically contains several terms specifically designed for different purposes. Parameters in front of these terms are used to control the relative weights among them. It is of critical importance to tune these parameters, as quality of the solution depends on their values. Tuning parameter is a relatively straightforward task for a human, as one can intelligently determine the direction of parameter adjustment based on the solution quality. Yet manual parameter tuning is not only tedious in many cases, but becomes impractical when a number of parameters exist in a problem. Aiming at solving this problem, this paper proposes an approach that employs deep reinforcement learning to train a system that can automatically adjust parameters in a human-like manner. We demonstrate our idea in an example problem of optimization-based iterative CT reconstruction with a pixel-wise total-variation regularization term. We set up a parameter tuning policy network (PTPN), which maps an CT image patch to an output that specifies the direction and amplitude by which the parameter at the patch center is adjusted. We train the PTPN via an end-to-end reinforcement learning procedure. We demonstrate that under the guidance of the trained PTPN for parameter tuning at each pixel, reconstructed CT images attain quality similar or better than in those reconstructed with manually tuned parameters.