Realistic simulation of cardiac magnetic resonance studies modeling anatomical variability, trabeculae, and papillary muscles.

Simulated magnetic resonance imaging brain studies have been generated for over a decade. Despite their useful potential, simulated cardiac studies are only emerging. This article focuses on the realistic simulation of cardiac magnetic resonance imaging datasets. The methodology is based on the XCAT phantom, which is modified to increase realism of the simulated images. Modifications include the modeling of trabeculae and papillary muscles based on clinical measurements and published data. To develop and evaluate our approach, the clinical database included 40 patients for anatomical measurements, 10 patients for papillary muscle modeling, and 10 patients for local gray value statistics. The virtual database consisted of 40 digital voxel phantoms. Histograms from different tissues were obtained from the real datasets and compared with histograms of the simulated datasets with the Chi-square dissimilarity metric (χ(2)) and Kullback-Leibler divergence. For the original phantom, χ(2) values averaged 0.65 ± 0.06 and Kullboek-Leibler values averaged 0.69 ± 0.38. For the modified phantom, χ(2) values averaged 0.34 ± 0.12 and Kullboek-Leibler values averaged 0.32 ± 0.15. The proposed approach demonstrated a noticeable improvement of the local appearance of the simulated images with respect to the ones obtained originally.

Simulation of late gadolinium enhancement cardiac magnetic resonance studies.

In this study we propose a pipeline for simulation of late gadolinium enhancement images. We used a modified version of the XCAT phantom to improve simulation realism. Modifications included the modeling of trabeculae and papillary muscles, and the increase of sublabels to resemble tissue intensity variability. Magnetic properties for each body tissue were sampled in three settings: from Gaussian distributions, combining Rayleigh-Gaussian distributions, and from Rayleigh distributions. Thirty-two simulated datasets were compared with 32 clinical datasets from infarcted patients. Histograms were obtained for five tissues: lung, pericardium, myocardium, blood and hyper-enhanced area. Real and simulated histograms were compared with the Chi-square dissimilarity metric (χ(2)) and Kullback-Leibler divergence (KL). The generated simulated images look similar to real images according to both metrics. Rayleigh and the Rayleigh-Gaussian models obtained comparable average results (respectively: χ(2)= 0.16 ± 0.12 and 0.18 ± 0.11; KL=0.15 ± 0.17 and 0.16 ± 0.18).