A Novel Phantom and a Dedicated Developed Software for Image Quality Controls in X-Ray Intraoral Devices.

BACKGROUND: Periodic quality control (QC) procedures are important in order to guarantee the image quality of radiological equipment and are also conducted using phantoms simulating human body.
OBJECTIVE: To perform (QC) measurements in intraoral imaging devices, a new and simple phantom was manufactured. Besides, to simplify QC procedures, computerized LabView-based software has been devised, enabling determination of image quantitative parameters in real time or during post processing.
MATERIAL AND METHODS: In this experimental study, the novel developed phantom consists of a Polymethyl methacrylate (PMMA) circular insert. It is able to perform a complete QC image program of X-ray intraoral equipment and also causes the evaluation of image uniformity, high and low contrast spatial resolution, image linearity and artefacts, with only two exposures.
RESULTS: Three raters analyzed the images using the LabView dedicated software and determined the quantitative and qualitative parameters in an innovative and accurate way. Statistical analysis evaluated the reliability of this study. Good accuracy of the quantitative and qualitative measurements for the different intraoral systems was obtained and no statistical differences were found using the inter-rater analysis.
CONCLUSION: The achieved results and the related statistical analysis showed the validity of this methodology, which could be proposed as an alternative to the commonly adopted procedures, and suggested that the novel phantom, coupled with the LabView based software, could be considered as an effective tool to carry out a QC image program in a reproducible manner. Copyright: © Journal of Biomedical Physics and Engineering.

Introduction

To ensure patients’ safety and maintain acceptable diagnostic equipment performances, national and international guidelines defined criteria and minimum accessibility values for all devices used in diagnostic exams and provided information about the Quality Assurance (QA) procedures during periodic quality controls (QC) [ 1 - 4 ]. Periodic QC procedures are important in order to guarantee consistent image quality of radiological equipment. Then QCs have an important role because they enable a complete evaluation of system status as well as image quality and allow the identification of degradation concerning the image quality itself and sources of any possible equipment malfunction, which may reduce the ability to detect and correctly interpret abnormal findings [ 5 , 6 ].

Walker et al., support the need for these procedures. In fact, since the exposure to different types of radiation must be as low as possible [ 7 ], these QCs cannot be done directly on patients. Moreover, if the device is not properly working, patients may be subjected to several risks.

Dental radiographic assessment represents one of the most frequently undertaken radiological examinations in the European Union [ 8 ]. As a matter of fact, it provides a useful aid in the diagnosis and treatment of oral diseases [ 9 ]. In this context, several professional groups, such as NCRP 145 and AAPM, proposed QA recommendations. Initially, the QC protocols were intended for film-based radiography: in fact, NCRP Report 145 refers to digital radiography by acknowledging that “the required standards, apparatus and software for dental system do not currently exist” [ 10 ]. In 2002, the American Association of Physicists in Medicine (AAPM) task group n. 175 report, concerning quality controls of dental image equipment was published [ 11 ].

However, with the development in computer technology and the introduction of digital systems, many specialists have converted from analog (film-based radiography) to digital imaging sensors. The latter has several advantages over film radiography (e.g.: immediate image production, interactive display and the possibility to store DICOM file format) [ 12 ]. Based on the image acquisition process, digital radiographic devices can be categorized as direct and indirect systems. The first ones acquire images using solid state detectors that are connected to a computer in order to display the image after exposure. The most used solid state detectors are charge couple devices (CCD) and complementary metal oxide semiconductors active pixel sensors (CMOS-APS). CCD is composed of an array of X-ray sensitive elements or wells on a silicon chip arranged in a rectangular matrix. In the readout process, resulting current is amplified, digitized, stored and displayed as digital images [ 13 ]. CMOS-APS are similar to CCDs, except they use an active pixel technology in which the pixels are isolated from their neighbors and directly accessed individually [ 14 ]. Indirect digital systems use photostimulable phosphor (PSP) plates. These imaging plates are coated with a radiation sensitive phosphor, storing a latent image after X-ray exposure [ 15 ]. The plate is scanned using a high-speed laser scanner, and the resulting light emitted by the stimulated phosphor is digitized and converted into displayable electronic information [ 16 ]. Today, a wide variety of intraoral digital systems is available in the dental market; in each case, the optimization of image reproducibility motivates the use of acquisition protocols that use dedicated phantoms to guarantee the maintenance of consistent image quality of the diagnostic equipment over lifetime [ 17 ].

Several useful phantoms have been built using different materials, such as Polymethyl methacrylate (PMMA), water, aluminum, etc., to simulate human body. These phantoms allow the measurement of different imaging parameters. For example, Mah et al., manufactured a phantom, in aluminum and PMMA to measure dynamic range, contrast perceptibility and spatial resolution at a source-to-detector distance that was representative of the one used in the clinical environment [ 18 ].

To study any possible source of error in intraoral radiographic exams, Yakoumakis et al., developed a phantom containing a real human tooth at its center, three air steps of different lengths and a background area for contrast and optical density measurements [ 19 ].

Phantoms are also useful to make comparison tests between different receptor systems; Yoshiura et al., manufactured a phantom in aluminum to compare film-based instrumentation performances with digital ones, obtaining a similar perceptibility curve in both cases and ensuring the validity of digital instruments as an alternative to classical diagnostic systems [ 20 ]. Dashpuntsag et al., performed an imaging quality comparison between the images acquired with photostimulable phosphor plates (PSPs) and the charge-coupled devices (CCDs) ones, showing that CCDs had better performances with respect to the PSPs, moreover a high variability in Gray value-dose relation was evaluated [ 21 ]. Furthermore, Pittayapat et al., tested different portable X-ray diagnostic systems using real human teeth as phantoms [ 22 ].

In this context, we developed a novel PMMA phantom to conduct image quality control of X-ray intraoral devices, regardless the X-ray source. This novel dedicated phantom was tested on different, direct and indirect, digital X-ray intraoral devices. In addition, in order to univocally perform the measurements in real time and/or during post processing, we also developed a dedicated LabView based software procedure. In fact, LabView allows easy interaction between acquisition instrumentation and operators [ 23 ]. It realizes virtual instruments in form of graphical programs running on personal computers or work- stations, which simplify the procedure for data acquisition, instrument control and/or data post-processing [ 24 ]. LabView implements a data flow paradigm in which the code is not written, but rather drawn or represented graphically, similarly to a flowchart diagram [ 25 , 26 ].

All the images were analyzed by three raters (two experienced and one novice) and the reliability of this study was evaluated by a detailed statistical analysis, better detailed in the following section.

Material and Methods

The focus of this experimental research regards the development of a novel phantom to conduct image quality controls on intraoral equipments.

Phantom description

The novel dedicated phantom consists of a PMMA circular insert, 3.5 cm in diameter. The phantom is 1.5 cm high and a PMMA grid with 0.2 × 0.2 mm2 square holes is manufactured on the top surface. In the middle of this grid, corresponding to the center of the phantom, a circular hole is present. This hole is dedicated to the measurement of High Contrast Spatial Resolution (HCSR), when a high-density object (i.e. an aluminum wire), 1.0 mm in diameter, is inserted. It is depicted in Figure 1 (a) and (b).

Figure 1

The novel dedicated phantom used for X-ray intraoral image quality controls. Visible view: Front (1 (a)) and lateral (1 (b)) view. X-Ray view Front (1 (c)) and lateral (1 (d)).

Digital imaging systems sensors

The study was conducted on five different intraoral devices, all using digital sensors. In particular, three direct digital systems (one Carestream CS 2100 and two Vatech EzRay Premium Vex S100W) and two indirect digital devices (one Orix 70 and one Kodak 2100) with PSP technology were analyzed.

Test procedure

The novel phantom is able to perform a complete QC image program of X-ray intraoral equipment and allows the evaluation of the following image parameters:

– Image Uniformity (IU)

– High Contrast Spatial Resolution (HCSR)

– Low Contrast Spatial Resolution (LCSR)

– Image Linearity (IL)

– Artefacts.

In order to determine all the necessary imaging parameters, the phantom needs to be exposed only twice. The two different exposures are shown in Figure 1. The first exposure (Figure 1(c)) allows the evaluation of HCSR and LCSR, IL and the possible presence of artefacts; while the second one (Figure 1(d)) allows us to calculate the IU.

LabView developed software

The method for data processing uses a software for quantitative and qualitative image measurements developed for a LabView environment. The software is compatible with nonstandard and standard image formats (DICOM, BMP, TIFF, JPEG, JPEG2000, PNG, and AIPD). In order to perform image quality measurements, the procedure begins by opening the image under investigation and, using the developed tool, quantitative information was obtained.

The developed software consists of the Front Panel, which represents the visual interface, and the block diagram, that represents the main core of the software.

Image measurements

Three raters, two medical physicists, P1 and P2, and one dentistry, D, performed the QC on the X-ray acquired images. The digital images were acquired with the manufacturers’ proprietary software and they were stored in TIFF format in order to avoid information loss. In post processing, the raters, using the developed LabView software, analyzed the stored images. The radiographic parameters set in order to image the phantom is reported in Table 1.

Table 1

Radiographic parameters set on different intraoral devices to acquire the images.

Device	Tube Voltage (kV)	Tube Current (mA)	Exposure Time (ms)
Kodak 210	60	7	200
Orix 70	70	8	100
Carestream CS 2100	60	7	100
Vatech EzRay Premium Vex S100W (1)	60	7	120
Vatech EzRay Premium Vex S100W (2)	60	7	120

Image Uniformity (IU) refers to the ability of the X-ray system to produce a constant signal response throughout the acquired object, imaged homogenous [ 27 ].

The IU measurement is performed using the dedicated section of the novel phantom depicted in Figure 1 (d). Three circular regions of interest (ROI) of diameter about 10% of the linear width of the phantom are drawn: the first one at the center of the phantom image, and the others on its left and right, respectively. For each ROI, a measurement of mean values of gray number is performed. The measurements are conducted using the LabView based software purposely developed. It opens the radiographic image, previously acquired, and measures the mean value of the gray numbers inside the ROI and its standard deviation (Figure 2 (a)). At acceptance, the difference between the central Gray number (G,c) and the peripheral Gray numbers (G,p) must not exceed 20 Gray Unit (GU).

Figure 2

The Front Panel of the developed software during Image Uniformity (IU) (2 (a)) and High Contrast Spatial Resolution (HCSR) (2 (b)) measurements.

|G,c-G,p|≤20 (1)

We chose the 20 GU difference value because it corresponds to about the 10% of the whole gray scale (i.e. 256 GU).

High Contrast Spatial Resolution (HCSR) measures the system’s ability to resolve high contrast objects. It is often determined having a large signal to noise ratio [ 28 ]. The HCSR measurement is performed using the dedicated section of the novel phantom depicted in Figure 1 (c) and the developed LabView based software procedure is used. A similar procedure was developed and validated in a previous study to determine the slice thickness in Computed Tomography (CT) and Magnetic Resonance (MR) devices [ 29 ]. For the evaluation of the HCSR, the procedure uses the “Line Profile” tool, which draws a line, along which gets the pixel values and intensities, on the opened image. Figure 2 (b) shows the “Front Panel” of the developed software, which allows opening the image, previously acquired, and drawing the “Line Profile”. In the same panel, the result of the line profile appears in real time as intensity vs. pixel number plot and the obtained curve is similar to an edge function.

The HCSR is obtained using the equation 2:

(2)

Where D is the phantom diameter in mm, L represents the distance, expressed in pixels, between the minimum and the maximum of the path of rapid change in image intensity, and P is the total pixels number calculated by the LabView software.

Low Contrast Spatial Resolution (LCSR) measures the system’s ability to resolve low contrast objects of small sizes. It measures the capability of an imaging system to distinguish objects when there is no significant noise contribution [ 27 ]. The LCSR is performed using the dedicated section of the novel phantom depicted in Figure 1 (c), consisting of a grid with 0.2 × 0.2 mm2 square holes. This PMMA grid has minimal differences (if converted to gray scale) from background. Different quantitative methods have been proposed. However, the most widely accepted [ 27 ] requires an observer that subjectively detects the objects (i.e. in this study the grid holes), as distinct. The grid with 0.2 × 0.2 mm2 square holes was chosen to augment the differences (in the gray scale representation) between the phantom and the background.

Image Linearity (IL) (or geometrical distortion) refers to the ability of the system to reproduce the dimensions of an object [ 30 ]. The IL measurement is performed using the dedicated section of the novel phantom depicted in Figure 1 (c) using the LabView software. The LabView software opens the previously acquired X-ray image and four diameter values are measured. The IL is defined by the equation 3:

(3)

Where D is the physical diameter of the novel phantom and is the mean value of the distances measured in the image, using the LabView’s ruler tool. At the acceptance, the maximum IL (%) must not differ from 5%.

Artefacts are features that appear in an image but not in the original imaged object. The identification of imaging artefacts depends on the rater’s ability, being often based on subjective criteria [ 31 ].

Statistical analysis

A further validation of the proposed methodology was the statistical analysis conducted [using the MATLAB ‘lillietest’ function, MATLAB® (MathWorks Inc., Natick, MA) ver. R2018b] on the resulting data set. A Lilliefors test was conducted on every set and each rater to ensure that the collected datasets follow a normal distribution [ 32 ].

In order to confirm the validity of the obtained results, three different statistical methods were performed: a one-way analysis of variance (ANOVA) for quantitative, normal distributed parameters, a post-hoc Tukey’s test for multiple comparison of different datasets, and an analysis with the Fleiss’ Kappa coefficient for the qualitative LCSR parameter, always using MATLAB® (MathWorks Inc., Natick, MA) ver. R2018b.

The one-way ANOVA is a statistical model, which allows comparing different and independent datasets of the same variable to assume whether they are generated by the same distribution; this is the ”null hypothesis” on which the test is based [ 33 ].

It can be seen as a special case of linear model, described by the equation 4:

(4)

Where y is the ith measure of the jth group, y is the population mean of the j group, and ε is the random error, independent and normally distributed, with zero mean and constant variance. Therefore, the “anova1” MATLAB function is used; if the p-value returned by the function is higher than 0.05, the test fails to reject the null hypothesis that indicates the two datasets refer to the same distribution, with a probability proportional to the p-value [ 34 ].

The post-hoc Tukey’s test is a single-step multiple comparison test used to compare all possible pairs of means in a group of different datasets and to find which ones of them statistically differ from the others [ 35 ]. The ‘multcompare’ MATLAB function is used [ 36 ].

The Fleiss’ Kappa coefficient is an inter-rater statistical coefficient, which measures the reliability of different sets of measures by different raters. It is an “extension” of Cohen’s Kappa coefficient that is valid for two raters only. It can be calculated from a matrix and Fleiss’ Kappa can be computed from this matrix to show the degree of agreement between the psychiatrists above the level of agreement expected by chance [ 37 ]. The expression for the Fleiss’ Kappa is (eq. 5):

(5)

Where is the relative observed agreement among raters, given by the equation 6:

(6)

With P rater-rater pair agreement for subject i, while P is the hypothetical probability of chance agreement, given by the equation 7:

(7)

With p proportion of judgements for class k.

In this case, a self-made MATLAB script is used to calculate the Fleiss’ Kappa coefficient (reported in supplementary material section).

Results

Quantitative and qualitative results are following. In Table 2, IU and HCSR measurements conducted on the Vatech EzRay Premium Vex S100W (1) are reported. In particular, for IU measurements, the differences of the peripheral ROIs compared to the central one are highlighted. Two different measurements were performed by each rater and, in all cases, both the differences are less than 20 GU. For each evaluated difference, the error is calculated by using the method of quadratic sum and its value is also reported. In Table 2, the results of the HCSR measurements, carried out by the three raters, using the LabView dedicated software, are shown, and an error of 5%, obtained from software development parameters, for each measurement, was considered. Each rater performed four measurements and the obtained mean value ± standard deviation (SD) is also calculated.

Table 2

Image Uniformity (IU) and High Contrast Spatial Resolution (HCSR) measurements conducted on Vatech EzRay Premium Vex S100W (1) device.

		P1	P2	D
IU measurements	|G,c-G,r| (GU) (1st)	|6.53 ± 4.68|	|8.86 ± 4.86|	|11.10 ± 4.90|
	|G,c-G,r| (GU) (2nd)	|6.79 ± 4.72|	|7.96 ± 4.86|	|10.42 ± 4.89|
	|G,c-G,l| (GU) (1st)	|3.24 ± 4.56|	|1.24 ± 4.58|	|1.45 ± 4.53|
	|G,c-G,l| (GU) (2nd)	|2.66 ± 4.57|	|2.10 ± 4.65|	|1.94 ± 4.58|
HCSR measurements	1st measure (10-1 mm)	1.30 ± 0.07	1.30 ± 0.07	1.20 ± 0.06
	2nd measure (10-1 mm)	1.20 ± 0.06	1.20 ± 0.06	1.20 ± 0.06
	3rd measure (10-1 mm)	1.20 ± 0.06	1.30 ± 0.07	1.10 ± 0.06
	4th measure (10-1 mm)	1.30 ± 0.07	1.20 ± 0.06	1.40 ± 0.07
	Mean ± SD (10-1 mm)	1.25 ± 0.05	1.25 ± 0.02	1.23 ± 0.11

IU: Image Uniformity, HCSR: High Contrast Spatial Resolution

The Lilliefors test, performed on HCSR results, assured that all the datasets had a normal distribution (p>0.05 in all datasets, 0.19 for P1, 0.19 for P2, and 0.13 for D). The results of the ANOVA test showed no statistical difference across datasets (p-value: 0.89) (Figure 3 (a)). In Table 3, the results of the DM measurements, conducted on the Vatech EzRay Premium Vex S100W(1) and carried out by the three raters using the LabView software, are shown. In this case, an error of 5%, obtained from software development parameters, for each measurement was considered. Each rater performed four measurements and, in Table 3, the IL value, obtained using equation. 3, is also calculated. The IL error is calculated using the standard propagation of errors and reported on Table 3.

Figure 3

Box Plot of High Contrast Spatial Resolution (HCSR) (3 (a)) and Image Linearity (IL) (3(b)) values. The central mark represents the median of the observables and the edges of the box are the 25th and 75th percentiles. The whiskers extend to the most extreme data points that are not considered outliers.

Table 3

Image Linearity (IL) measurements conducted on Vatech EzRay Premium Vex S100W(1) equipment.

	P1	P2	D
1st measure of DM (mm)	7.99 ± 0.40	8.00 ± 0.40	7.99 ± 0.40
2nd measure of DM (mm)	8.00 ± 0.40	7.99 ± 0.40	7.99 ± 0.40
3rd measure of DM (mm)	8.00 ± 0.40	8.00 ± 0.40	8.00 ± 0.40
4th measure of DM (mm)	7.99 ± 0.40	7.99 ± 0.40	7.98 ± 0.40
Image Linearity (%)	0.06 ± 0.01	0.06 ± 0.01	0.13 ± 0.01

The Lilliefors test identified that all the datasets had a normal distribution (p>0.05 in all datasets, 0.19 for P1, 0.19 for P2, and 0.50 for D). The ANOVA test provided, for the matrix containing the three datasets, a p-value of 0.50 (p > 0.05). In Figure 3 (b), the plot of the ANOVA test is reported. In this case, it is highlighted that no statistical differences are present.

In Table 4, the results of the image QC conducted on all the considered X-ray devices and the evaluation of each rater are reported. For the IU, the maximum difference values are considered for each device. In addition, the ROI mean values strictly depend on the Window Gray Level, set on the X-ray dental acquired images.

Table 4

Quantitative quality control (QC) measurements conducted on all considered devices.

Device	Rater	IU (GU)	HCSR (10-1 mm)	IL (%)
Kodak 210	P1	|17.79 ± 7.41|	0.83 ± 0.05	0.07 ± 0.01
	P2	|14.62 ± 6.27|	0.83 ± 0.06	0.07 ± 0.01
	D	|15.38 ± 6.25|	0.90 ± 0.05	0.08 ± 0.01
Orix 70	P1	|18.07 ± 7.62|	0.85 ± 0.10	0.05 ± 0.01
	P2	|18.52 ± 7.83|	0.70 ± 0.08	0.05 ± 0.01
	D	|16.35 ± 6.84|	0.86 ± 0.07	0.05 ± 0.01
Carestream CS 2100	P1	|5.58 ± 4.23|	1.40 ± 0.06	0.06 ± 0.01
	P2	|4.94 ± 3.54|	1.43 ± 0.06	0.06 ± 0.01
	D	|5.24 ± 3.76|	1.30 ± 0.05	0.08 ± 0.01
Vatech EzRay Premium Vex S100W (1)	P1	|6.79 ± 4.72|	1.25 ± 0.05	0.06 ± 0.01
	P2	|8.86 ± 4.86|	1.25 ± 0.02	0.06 ± 0.01
	D	|11.10 ± 4.90|	1.23 ± 0.11	0.13 ± 0.01
Vatech EzRay Premium Vex S100W (2)	P1	|7.42 ± 4.07|	1.20 ± 0.05	0.07 ± 0.01
	P2	|8.96 ± 4.64|	1.30 ± 0.05	0.09 ± 0.01
	D	|8.64 ± 4.18|	1.30 ± 0.03	0.06 ± 0.01

IU: Image Uniformity, HCSR: High Contrast Spatial Resolution, IL: Image Linearity

For each device, in order to compare the dataset of every physical parameter, the ANOVA test was performed. The obtained p-values are < 0.001 (IU and HCSR) and 0.39 (IL): it is possible to notice that the values for IU and HCSR do not follow the same distribution. Then we performed a post-hoc multiple comparison using the Tukey’s test, and we found statistical differences only between indirect digital devices and direct digital devices (p < 0.001 for all the comparisons), while we did not find statistical differences between devices of the same type.

In Figure 4, the box plots of IU (4 (a)) and HCSR (4(b)) are reported. In the Figure 4 , it is possible to notice how devices with the same receptor technology are statistically correlated.

Figure 4

Tukey’s test results that put in evidence how only devices of the same type (direct or indirect systems) are statistical correlated. (a) is referred to Image Uniformity (IU), while (b) is about High Contrast Spatial Resolution (HCSR).

Statistical analysis of the LCSR qualitative parameter was conducted using the Fleiss’ Kappa coefficient. The three raters were able to distinguish the visible square grid, except for the Orix 70 X-ray images, in which one rater was not able to clearly distinguish the grid. Therefore, we used a matrix with two columns (the two classes “grid” and “not grid”) and ten rows (two rows for every image obtained with a specific device, one row for the grid region and the other for the non-grid region) in order to calculate the Fleiss’ Kappa coefficient. The Fleiss’ Kappa resulting from this matrix is 0.86: this means that there is a high grade of agreement among raters. The three raters did not detect artefacts among the acquired images; for this reason no statistical analysis was conducted on this parameter.

Discussion

In this study, a novel phantom was developed to conduct a complete image QC program on intraoral devices. The reliability of this novel phantom was evaluated on five intraoral devices different in receptor technology. The developed phantom allows us to obtain information about the following image parameters: Image Uniformity, High Contrast Spatial Resolution, Low Contrast Spatial Resolution, Image Linearity and Artifacts. In this way, a complete QC program could be carried out with only two exposures. Moreover, the novel phantom is coupled with a dedicated LabView software, which provides an innovative, accurate and automated determination of the quantitative parameters, in real time and/or during post processing, with an intuitive visual programming language. The great advantage of the LabView software is that it is adaptable and modifiable according to the operator’s needs together with providing a huge free documentation. Thus, no dedicated software and phantom, distributed by manufacturers, are necessary.

Conclusion

QCs represent an important instrument to monitor the image quality and verify the maintenance of several physical characteristics of the medical imaging devices. Due to the impossibility of performing QCs directly on patients for safety issues, a sort of replacement is needed: for this reason, phantoms are realized in order to simulate human body or parts of it, observing the standard criteria defined by the various protocols for any application. Usually commercial QC phantoms are expensive and often suitable for evaluating just a few parameters during image quality test. Therefore, several dedicated phantoms are necessary in order to perform a complete image QC program.

With this research, we tried successfully to realize a phantom able to unify all the necessary characteristics for a complete image QC program. We obtained a good accuracy of both quantitative and qualitative measurements for all of the different systems, and no statistical differences were found between the raters with the inter-rater analysis, suggesting that a single phantom can be used to evaluate the performance of the X-ray intraoral devices.

The obtained results and the related statistical analysis showed the validity of this methodology and suggested that the novel phantom could be considered as an effective tool that allows an immediate determination of the image parameters interested in the image QC program; furthermore, it is easy to be used and could represent a reproducible instrument for routinely QC programs.