Towards guidelines to harmonize textural features in PET: Haralick textural features vary with image noise, but exposure-invariant domains enable comparable PET radiomics.

Prenosil, George Amadeus; Weitzel, Thilo; Fürstner, Markus; Hentschel, Michael; Krause, Thomas; Cumming, Paul; Rominger, Axel; Klaeser, Bernd (2020). Towards guidelines to harmonize textural features in PET: Haralick textural features vary with image noise, but exposure-invariant domains enable comparable PET radiomics. PLoS ONE, 15(3), e0229560. Public Library of Science 10.1371/journal.pone.0229560

[img]
Preview
Text
pone.0229560.pdf - Published Version
Available under License Creative Commons: Attribution (CC-BY).

Download (3MB) | Preview

PURPOSE

Image texture is increasingly used to discriminate tissues and lesions in PET/CT. For quantification or in computer-aided diagnosis, textural feature analysis must produce robust and comparable values. Because statistical feature values depend on image count statistics, we investigated in depth the stability of Haralick features values as functions of acquisition duration, and for common image resolutions and reconstructions.

METHODS

A homogeneous cylindrical phantom containing 9.6 kBq/ml Ge-68 was repeatedly imaged on a Siemens Biograph mCT, with acquisition durations ranging from three seconds to three hours. Images with 1.5, 2, and 4 mm isometrically spaced voxels were reconstructed with filtered back-projection (FBP), ordered subset expectation maximization (OSEM), and the Siemens TrueX algorithm. We analysed Haralick features derived from differently quantized (3 to 8-bit) grey level co-occurrence matrices (GLCMs) as functions of exposure E, which we defined as the product of activity concentration in a volume of interest (VOI) and acquisition duration. The VOI was a 50 mm wide cube at the centre of the phantom. Feature stability was defined for df/dE → 0.

RESULTS

The most stable feature values occurred in low resolution FBPs, whereas some feature values from 1.5 mm TrueX reconstructions ranged over two orders of magnitude. Within the same reconstructions, most feature value-exposure curves reached stable plateaus at similar exposures, regardless of GLCM quantization. With 8-bit GLCM, median time to stability was 16 s and 22 s for FBPs, 18 s and 125 s for OSEM, and 23 s, 45 s, and 76 s for PSF reconstructions, with longer durations for higher resolutions. Stable exposures coincided in OSEM and TrueX reconstructions with image noise distributions converging to a Gaussian. In FBP, the occurrence of stable values coincided the disappearance of negatives image values in the VOI.

CONCLUSIONS

Haralick feature values depend strongly on exposure, but invariance exists within defined domains of exposure. Here, we present an easily replicable procedure to identify said stable exposure domains, where image noise does not substantially add to textural feature values. Only by imaging at predetermined feature-invariant exposure levels and by adjusting exposure to expected activity concentrations, can textural features have a quantitative use in PET/CT. The necessary exposure levels are attainable by modern PET/CT systems in clinical routine.

Item Type:

Journal Article (Review Article)

Division/Institute:

04 Faculty of Medicine > Department of Radiology, Neuroradiology and Nuclear Medicine (DRNN) > Clinic of Nuclear Medicine

UniBE Contributor:

Prenosil, George; Weitzel, Thilo; Fürstner, Markus; Hentschel, Michael; Krause, Thomas Michael; Cumming, Paul; Rominger, Axel Oliver and Klaeser, Bernd

Subjects:

600 Technology > 610 Medicine & health

ISSN:

1932-6203

Publisher:

Public Library of Science

Language:

English

Submitter:

Sabine Lanz

Date Deposited:

29 Dec 2020 15:10

Last Modified:

03 Jan 2021 02:53

Publisher DOI:

10.1371/journal.pone.0229560

PubMed ID:

32176698

BORIS DOI:

10.48350/149138

URI:

https://boris.unibe.ch/id/eprint/149138

Actions (login required)

Edit item Edit item
Provide Feedback