We propose a new single image super resolution (SR) algorithm via Bayesian modeling with a natural image prior modeled by a high-order Markov random field (MRF). SR is one of the long-standing and active topics in image processing community. It is of great use in many practical applications, such as astronomical observation, medical imaging, and the adaptation of low-resolution contents onto high-resolution displays. One category of the conventional approaches for image SR is formulating the problem with Bayesian modeling techniques and then obtaining its maximum-a-posteriori solution, which actually boils down to a regularized regression task. Although straightforward, this approach cannot exploit the full potential offered by the probabilistic modeling, as only the posterior mode is sought. On the other hand, current Bayesian SR approaches using the posterior mean estimation typically use very simple prior models for natural images to ensure the computational tractability. In this paper, we present a Bayesian image SR approach with a flexible high-order MRF model as the prior for natural images. The minimum mean square error (MMSE) criteria are used for estimating the HR image. A Markov chain Monte Carlo-based sampling algorithm is presented for obtaining the MMSE solution. The proposed method cannot only enjoy the benefits offered by the flexible prior, but also has the advantage of making use of the probabilistic modeling to perform a posterior mean estimation, thus is less sensitive to the local minima problem as the MAP solution. Experimental results indicate that the proposed method can generate competitive or better results than state-of-the-art SR algorithms.

No abstract

Image-guided adaptive radiotherapy requires deformable image registration to map radiation dose back and forth between images. The purpose of this study is to develop a novel method to improve the accuracy of an intensity-based image registration algorithm in low-contrast regions. A computational framework has been developed in this study to improve the quality of the “demons” registration. For each voxel in the registration’s target image, the standard deviation of image intensity in a neighborhood of this voxel was calculated. A mask for high-contrast regions was generated based on their standard deviations. In the masked regions, a tetrahedral mesh was refined recursively so that a sufficient number of tetrahedral nodes in these regions can be selected as driving nodes. An elastic system driven by the displacements of the selected nodes was formulated using a finite element method (FEM) and implemented on the refined mesh. The displacements of these driving nodes were generated with the “demons” algorithm. The solution of the system was derived using a conjugated gradient method, and interpolated to generate a displacement vector field for the registered images. The FEM correction method was compared with the “demons” algorithm on the CT images of lung and prostate patients. The performance of the FEM correction relating to the “demons” registration was analyzed based on the physical property of their deformation maps, and quantitatively evaluated through a benchmark model developed specifically for this study. Compared to the benchmark model, the “demons” registration has the maximum error of 1.2 cm, which can be corrected by the FEM method to 0.4 cm, and the average error of the “demons” registration is reduced from 0.17 cm to 0.11 cm. For the CT images of lung and prostate patients, the deformation maps generated by the “demons” algorithm were found unrealistic at several places. In these places, the displacement differences between the “demons” registrations and their FEM corrections were found in the range of 0.4 cm and 1.1cm. The mesh refinement and FEM simulation were implemented in a single thread application which requires about 45 minutes of computation time on a 2.6 GH computer. This study has demonstrated that the finite element method can be integrated with intensity-based image registration algorithms to improve their registration accuracy, especially in low-contrast regions.

The purpose of this study is to characterize the dosimetric properties and accuracy of a novel treatment platform (Edge radiosurgery system) for localizing and treating patients with frameless, image‐guided stereotactic radiosurgery (SRS) and stereotactic body radiotherapy (SBRT). Initial measurements of various components of the system, such as a comprehensive assessment of the dosimetric properties of the flattening filter‐free (FFF) beams for both high definition (HD120) MLC and conical cone‐based treatment, positioning accuracy and beam attenuation of a six degree of freedom (6DoF) couch, treatment head leakage test, and integrated end‐to‐end accuracy tests, have been performed. The end‐to‐end test of the system was performed by CT imaging a phantom and registering hidden targets on the treatment couch to determine the localization accuracy of the optical surface monitoring system (OSMS), cone‐beam CT (CBCT), and MV imaging systems, as well as the radiation isocenter targeting accuracy. The deviations between the percent depth‐dose curves acquired on the new linac‐based system (Edge), and the previously published machine with FFF beams (TrueBeam) beyond Dmax were within 1.0% for both energies. The maximum deviation of output factors between the Edge and TrueBeam was 1.6%. The optimized dosimetric leaf gap values, which were fitted using Eclipse dose calculations and measurements based on representative spine radiosurgery plans, were 0.700 mm and 1.000 mm, respectively. For the conical cones, 6X FFF has sharper penumbra ranging from 1.2−1.8 mm (80%‐20%) and 1.9−3.8 mm (90%‐10%) relative to 10X FFF, which has 1.2−2.2 mm and 2.3−5.1 mm, respectively. The relative attenuation measurements of the couch for PA, PA (rails‐in), oblique, oblique (rails‐out), oblique (rails‐in) were: −2.0%, −2.5%, −15.6%, −2.5%, −5.0% for 6X FFF and −1.4%, −1.5%, −12.2%, −2.5%, −5.0% for 10X FFF, respectively, with a slight decrease in attenuation versus field size. The systematic deviation between the OSMS and CBCT was −0.4±0.2 mm, 0.1±0.3 mm, and 0.0±0.1 mm in the vertical, longitudinal, and lateral directions. The mean values and standard deviations of the average deviation and maximum deviation of the daily Winston‐Lutz tests over three months are 0.20±0.03 mm and 0.66±0.18 mm, respectively. Initial testing of this novel system demonstrates the technology to be highly accurate and suitable for frameless, linac‐based SRS and SBRT treatment.PACS number: 87.56.J‐

The purpose of this study was to perform comprehensive measurements and testing of a Novalis Tx linear accelerator, and to develop technical guidelines for commissioning from the time of acceptance testing to the first clinical treatment. The Novalis Tx (NTX) linear accelerator is equipped with, among other features, a high‐definition MLC (HD120 MLC) with 2.5 mm central leaves, a 6D robotic couch, an optical guidance positioning system, as well as X‐ray‐based image guidance tools to provide high accuracy radiation delivery for stereotactic radiosurgery and stereotactic body radiation therapy procedures. We have performed extensive tests for each of the components, and analyzed the clinical data collected in our clinic. We present technical guidelines in this report focusing on methods for: (1) efficient and accurate beam data collection for commissioning treatment planning systems, including small field output measurements conducted using a wide range of detectors; (2) commissioning tests for the HD120 MLC; (3) data collection for the baseline characteristics of the on‐board imager (OBI) and ExacTrac X‐ray (ETX) image guidance systems in conjunction with the 6D robotic couch; and (4) end‐to‐end testing of the entire clinical process. Established from our clinical experience thus far, recommendations are provided for accurate and efficient use of the OBI and ETX localization systems for intra‐ and extracranial treatment sites. Four results are presented. (1) Basic beam data measurements: Our measurements confirmed the necessity of using small detectors for small fields. Total scatter factors varied significantly (30% to approximately 62%) for small field measurements among detectors. Unshielded stereotactic field diode (SFD) overestimated dose by ~ 2% for large field sizes. Ion chambers with active diameters of 6 mm suffered from significant volume averaging. The sharpest profile penumbra was observed for the SFD because of its small active diameter (0.6 mm). (2) MLC commissioning: Winston Lutz test, light/radiation field congruence, and Picket Fence tests were performed and were within criteria established by the relevant task group reports. The measured mean MLC transmission and dynamic leaf gap of 6 MV SRS beam were 1.17% and 0.36 mm, respectively. (3) Baseline characteristics of OBI and ETX: The isocenter localization errors in the left/right, posterior/anterior, and superior/inferior directions were, respectively, −0.2±0.2 mm, −0.8±0.2 mm, and −0.8±0.4 mm for ETX, and 0.5±0.7 mm, 0.6±0.5 mm, and 0.0±0.5 mm for OBI cone‐beam computed tomography. The registration angular discrepancy was 0.1±0.2°, and the maximum robotic couch error was 0.2°. (4) End‐to‐end tests: The measured isocenter dose differences from the planned values were 0.8% and 0.4%, measured respectively by an ion chamber and film. The gamma pass rate, measured by EBT2 film, was 95% (3% DD and 1 mm DTA). Through a systematic series of quantitative commissioning experiments and end‐to‐end tests and our initial clinical experience, described in this ...