The statistical method (SM) applied in this work is described in detail. It is generally acknowledged that LCD is influenced by both, the amount of noise and the noise frequency spectrum. Hence the standard deviation σ of the pixel values in a homogenous image is not by itself useful in determining LCD. However, considering several identical regions of interest (ROI) placed on a uniform background, each with a mean pixel value µ, the standard deviation of these means (σµ) is relevant when determining the LCD, since the ROI limits noise to the spatial frequencies of the object.
Since, according to the central limit theorem, is normally distributed, a prediction on the detectability of an insert, identical in size to the considered ROI, can be made. The LCth characterised by 95% confidence level is given by
LCth=3.29*σµ,
where σµ is the standard deviation of all means from ROIs of the same size as the insert under consideration.
The statistical method is applicable to any homogenous test object. The simple test object used in this work consisted of an acetate sheet with the following details fixed on it with thin tape on the detail borders:
a central 99.99% pure aluminium platelet (3 x 3 cm2, 0.49 mm thick) to determine LCth in terms of %contrast;
a 99.99% pure aluminium step wedge (from 0.20 to 1.00 mm in 0.20 mm steps, each step is 0.6 x 0.8 cm2) for correlation of LCth to the absolute physical quantity “aluminium thickness” (mm Al);
a lead platelet (2 x 2 cm2, 0.5 mm thick) in lateral position for automatic pattern recognition.
Radiografic image of the simple homemade test object for the SM method
A Matlab code for automated analysis of LCD test object images was developed. A square ROI (120 x 120 pixels) was used to evaluate LCth according to the statistical method. To decrease uncertainty, the evaluation was carried out for 9 non overlapping ROIs arranged in the central uniform Al region and the average value LCth considered for each detail size.
In order to express LCth in terms of a more physical quantity, namely threshold Al thickness, the correlation between image contrast and aluminium thickness has been determined. For that purpose, square ROIs of 50 pixels were placed on each step of the Al step wedge and linear regression between mean pixel value and Al thickness calculated for each image.
The Statistical LCD calculation algorithm consists in a pixel-binning process repeated for different pixel binning matrix sizes (2x2, 4x4, etc).
Finally, the so determined data points (LCth) were interpolated to calculate LCth values for detail sizes corresponding to those in the CDMAM phantom or CDRAD phantom.
To further quantify IQ, for each contrast detail curve a summarising figure of merit was used, the Inverse Image Quality Figure IQFinv, calculated as:
IQFinv=1/(∑di*LCth),
where diameter and contrast threshold refer to the contrast detail curve.
Villa R, Oberhofer N, Paruccini N, 2020. "Validation of a statistical method for low contrast detectability as a simple tool for QC in digital mammography," Proc. SPIE 11513, 15th International Workshop on Breast Imaging (IWBI2020), 115131K (22 May 2020); https://doi.org/10.1117/12.2563987
Villa R., Paruccini N., Baglivi A., Signoriello M., Montezuma Velasquez R.A., Morzenti S., De Ponti E., Crespi, A., 2019. "Model observers for Low Contrast Detectability evaluation in dynamic angiography: A feasible approach". Physica Medica 64, 89–97. doi:10.1016/j.ejmp.2019.06.015
Spadavecchia C., Villa R., Pasquali C., Paruccini N., Oberhofer N., Crespi, A., 2016. "A statistical method for low contrast detectability assessment in digital mammography", in: Lecture Notes in Computer Science (Including Subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics). Springer Verlag, pp. 532–539. doi:10.1007/978-3-319-41546-8_67