Substitute CT generation from a single ultra short time echo MRI sequence: Preliminary study

Soumya Ghose, Jason A. Dowling, Robba Rai, Gary P. Liney

Research output: Contribution to journalOther Journal Articlepeer-review

11 Citations (Scopus)


In MR guided radiation therapy planning both MR and CT images for a patient are acquired and co-registered to obtain a tissue specific HU map. Generation of the HU map directly from the MRI would eliminate the CT acquisition and may improve radiation therapy planning. In this preliminary study of substitute CT (sCT) generation, two porcine leg phantoms were scanned using a 3D ultrashort echo time (PETRA) sequence and co-registered to corresponding CT images to build tissue specific regression models. The model was created from one co-registered CT-PETRA pair to generate the sCT for the other PETRA image. An expectation maximization based clustering was performed on the co-registered PETRA image to identify the soft tissues, dense bone and air class membership probabilities. A tissue specific non linear regression model was built from one registered CT-PETRA pair dataset to predict the sCT of the second PETRA image in a two-fold cross validation schema. A complete substitute CT is generated in 3 min. The mean absolute HU error for air was 0.3 HU, bone was 95 HU, fat was 30 HU and for muscle it was 10 HU. The mean surface reconstruction error for the bone was 1.3 mm. The PETRA sequence enabled a low mean absolute surface distance for the bone and a low HU error for other classes. The sCT generated from a single PETRA sequence shows promise for the generation of fast sCT for MRI based radiation therapy planning.

Original languageEnglish
Pages (from-to)2950-2960
Number of pages11
JournalPhysics in Medicine and Biology
Issue number8
Publication statusPublished - 17 Mar 2017
Externally publishedYes


Dive into the research topics of 'Substitute CT generation from a single ultra short time echo MRI sequence: Preliminary study'. Together they form a unique fingerprint.

Cite this