Structured errors in reconstruction methods for Non-Cartesian MR data

OPERA: a novel method to reduce ghost and aliasing artifacts

OBJECTIVE

A method for Orthogonal Phase Encoding Reduction of Artifact (OPERA) was developed and tested.

MATERIALS AND METHODS

Because the position of ghosts and aliasing artifacts is predictable along columns or rows, OPERA combines the intensity values of two images acquired using the same parameters, but with swapped phase-encoding directions, to correct the artifacts. Simulations and phantom experiments were conducted to define the efficacy, robustness, and reproducibility. Clinical validation was performed on a total of 1003 images by comparing the OPERA-corrected images and the corresponding image standard in terms of Signal-to-Noise Ratio (SNR) and Contrast-to-Noise Ratio (CNR). The method efficacy was also rated using a Likert-type scale response by two experienced independent radiologists using a single-blinded procedure.

RESULTS

Simulations and phantom experiments demonstrated the robustness and effectiveness of OPERA in reducing artifacts strength. OPERA application did not significantly change the SNR [+ 4.16%; inter-quartile range (IQR): 2.72-5.01%] and CNR (+ 4.30%; IQR: 2.86-6.04%) values. The two radiologists observed a total of 893 original images with artifacts (89.03% of the total images), a reduction in the perceived artifacts of 82.0% and 83.9% (p < 0.0001), and an improvement in the perceived SNR (82.8% and 88.5%; K = 0.714) and perceived CNR (86.9-88.9%; K = 0.722).

DISCUSSION

The study demonstrated that OPERA reduces MR artifacts and improves the perceived image quality.

A Conway-Maxwell-Poisson (CMP) model to address data dispersion on positron emission tomography

Positron emission tomography (PET) in medicine exploits the properties of positron-emitting unstable nuclei. The pairs of γ- rays emitted after annihilation are revealed by coincidence detectors and stored as projections in a sinogram. It is well known that radioactive decay follows a Poisson distribution; however, deviation from Poisson statistics occurs on PET projection data prior to reconstruction due to physical effects, measurement errors, correction of deadtime, scatter, and random coincidences. A model that describes the statistical behavior of measured and corrected PET data can aid in understanding the statistical nature of the data: it is a prerequisite to develop efficient reconstruction and processing methods and to reduce noise. The deviation from Poisson statistics in PET data could be described by the Conway-Maxwell-Poisson (CMP) distribution model, which is characterized by the centring parameter λ and the dispersion parameter ν, the latter quantifying the deviation from a Poisson distribution model. In particular, the parameter ν allows quantifying over-dispersion (ν<1) or under-dispersion (ν>1) of data. A simple and efficient method for λ and ν parameters estimation is introduced and assessed using Monte Carlo simulation for a wide range of activity values. The application of the method to simulated and experimental PET phantom data demonstrated that the CMP distribution parameters could detect deviation from the Poisson distribution both in raw and corrected PET data. It may be usefully implemented in image reconstruction algorithms and quantitative PET data analysis, especially in low counting emission data, as in dynamic PET data, where the method demonstrated the best accuracy.