Estimation of Hounsfield unit conversion parameters for pelvic CT images

Linear scaling is used to convert raw computed tomography (CT) pixel values into Hounsfield units corresponding to different tissue values. Analysis of a benchmarking study, presented here, where the same CT scan was imported into and then exported from multiple radiotherapy treatment planning systems, found inconsistencies in HU scaling parameter values exported along with the images, particularly when images were exported using the Radiation Therapy Oncology Group format. Several methods of estimating conversion parameters, based on estimating pixel values corresponding to air and water within the image, for pelvic CT images from a large multi-centre trial were compared against original Digital Imaging and Communications in Medicine export parameters. In general using the mean of a sample region at the centroid of the bladder to estimate the value of water was more accurate than using the minimum or maximum or a single value at the centroid. Accuracy of methods of air estimation tested were dependent in part on features of the CT scanners and treatment planning systems, making it difficult to pick one method as superior that was independent of scanner and treatment planning system type. Based on the above analysis, methods for estimating air and water were selected for use in performing linear scaling of a set of pelvic CT images prior to their use in an interpatient image registration application. The selected methods were validated against a more recent and homogeneous dataset. Estimation error was found to be much lower within the validation set.