Development of a neonate X-ray phantom for 2D imaging applications using single-tone inkjet printing

Purpose Inkjet printers can be used to fabricate anthropomorphic phantoms by the use of iodine-doped ink. However, challenges persist in implementing this technique. The calibration from grayscale to ink density is complex and time-consuming. The purpose of this work is to develop a printing methodology that requires a simpler calibration and is less dependent on printer characteristics to produce the desired range of x-ray attenuation values. Methods Conventional grayscale printing was substituted by single-tone printing; that is, the superposition of pure black layers of iodinated ink. Printing was performed with a consumer-grade inkjet printer using ink made of potassium-iodide (KI) dissolved in water at 1 g/ml. A calibration for the attenuation of ink was measured using a commercial x-ray system at 70 kVp. A neonate radiograph obtained at 70 kVp served as an anatomical model. The attenuation map of the neonate radiograph was processed into a series of single-tone images. Single-tone images were printed, stacked, and imaged at 70 kVp. The phantom was evaluated by comparing attenuation values between the printed phantom and the original radiograph; attenuation maps were compared using the structural similarity index measure (SSIM), while attenuation histograms were compared using the Kullback-Leibler (KL) divergence. A region of interest (ROI)-based analysis was also performed, where the attenuation distribution within given ROIs was compared between phantom and patient. The phantom sharpness was evaluated in terms of modulation transfer function (MTF) estimates and signal spread profiles of high spatial resolution features in the image. Results The printed phantom required 36 pages. The printing queue was automated and it took about 2 h to print the phantom. The radiograph of the printed phantom demonstrated a close resemblance to the original neonate radiograph. The SSIM of the phantom with respect to that of the patient was 0.53. Both patient and phantom attenuation histograms followed similar distributions, and the KL divergence between such histograms was 0.20. The ROI-based analysis showed that the largest deviations from patient attenuation values were observed at the higher and lower ends of the attenuation range. The limiting resolution of the proposed methodology was about 1 mm. Conclusion A methodology to generate a neonate phantom for 2D imaging applications, using single-tone printing, was developed. This method only requires a single-value calibration and required less than 2 h to print a complete phantom.