Home Prediction of Perceptible Artifacts in JPEG2000 Compressed Abdomen CT Images Using a Perceptual Image Quality Metric
Post
Cancel

Prediction of Perceptible Artifacts in JPEG2000 Compressed Abdomen CT Images Using a Perceptual Image Quality Metric

Rationale and Objectives

To test a perceptual quality metric (high-dynamic range visual difference predictor, HDR-VDP) in predicting perceptible artifacts in Joint Photographic Experts Group 2000 compressed thin- and thick-section abdomen computed tomography images.

Materials and Methods

A total of 120 thin (0.67 mm) and corresponding thick (5 mm) sections were compressed to ratios from 4:1 to 15:1. Peak signal-to-noise ratio (PSNR), HDR-VDP results (paired t -tests), and five radiologists’ pooled responses for the presence of artifacts (exact tests for paired proportions) were compared between the thin and thick sections. For three subsets of 120 thin- (subset A), 120 thick- (subset B), and 60 thin- and 60 thick-section compressed images (subset C), receiver operating curve analysis was performed to compare PSNR and HDR-VDP in predicting the radiologists’ responses. Using the cutoff values where the sum of sensitivity and specificity was the maximum in subset C, visually lossless thresholds (VLTs) were estimated for the 240 original images and the estimation accuracy was compared (McNemar test).

Results

Thin sections showed more artifacts in terms of PSNR, HDR-VDP, and radiologists’ responses ( p < .0001). HDR-VDP outperformed PSNR for subset C (area under the curve: 0.97 versus 0.93, p = 0.03), whereas they did not differ significantly for subset A or B. Using the cutoff values, PSNR and HDR-VDP predicted the VLT accurately for 124 (51.7%) and 183 (76.3%) images, respectively ( p < .0001).

Conclusions

HDR-VDP can predict the perceptible compression artifacts, and therefore can be potentially used to estimate the VLT for such compressions.

Irreversible image compression appears to be an immediate and effective means to reduce enormous data ( ) generated by modern computed tomography (CT) scanners ( ). Previous studies ( ) have reported acceptable compression levels for CT images as 8:1 to 20:1. Most of these investigated a compression threshold that does not cause a loss of diagnostic information. However, it is difficult to generalize such study results, because the diagnostically lossless threshold varies with diagnostic tasks ( ). For instance, it has been reported that detection performance of CT for hepatic nodules is preserved with up to 10:1 compression ( ); however, it is uncertain whether this compression level is acceptable for the characterization of the nodules and for the detection of any coincidental findings which might be clinically important in the same CT dataset. Furthermore, compression artifacts are affected by image content itself ( ) and scanning parameters such as section thickness ( ). Therefore it is very unlikely that a reported acceptable compression level for images of a certain type (eg, abdomen CT) can be an optimal compression guideline for all images of the same type.

If a compressed image is indistinguishable from its original, there is no basis for arguing the compression hinders any diagnostic accuracy ( ). Although this visually lossless threshold (VLT) allows relatively low-level compressions, it has been gaining support as a practicable compression level ( ). To estimate the VLT, human readers need to determine whether a compressed image is distinguishable from its original or not at various compression levels, which seems impractical in clinical practice.

Get Radiology Tree app to read full this article<

Materials and methods

Get Radiology Tree app to read full this article<

CT Scanning

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Table 1

Lesions in 120 Original Thick Sections

Lesion Number of Images Normal 70 Cystic focal lesion in the solid organ 6 Arterial luminal irregularity or wall calcification 5 Peritoneal carcinomatosis or peritonitis 5 Tube, anastomosis, or other postoperative finding 5 Bowel wall thickening 4 Solid focal lesion in the solid organ 4 Acute appendicitis 3 Ascites or pleural effusion 3 Acute pancreatitis or retroperitoneal edema 2 Colonic diverticulosis or diverticulitis 2 Focal calcification or stone 2 Gall bladder distension or wall thickening 2 Lymph node enlargement 2 Organ enlargement or atrophy 2 Severe bowel distension 2 Biliary dilatation 1

Get Radiology Tree app to read full this article<

Image Compression

Get Radiology Tree app to read full this article<

Table 2

Actual Compression Levels, PSNR and HDR-VDP Results, and Pooled Readers’ Responses

Nominal Compression Level Actual Compression Level ⁎ PSNR (dB) HDR-VDP Readers’ Responses † Thin Thick Thin Thick Thin Thick Thin Thick 4:1 4.0 ± 0.1 — 52.9 ± 2.4 — 2.8 ± 2.4 — 0% (0/120) — 6:1 6.0 ± 0.1 6.2 ± 0.4 46.3 ± 3.0 53.7 ± 1.7 21.0 ± 13.4 3.4 ± 2.1 35.0% (42/120) 0% (0/120) 8:1 8.0 ± 0.0 8.0 ± 0.0 42.3 ± 3.1 50.7 ± 2.5 50.6 ± 21.1 12.4 ± 7.5 95.0% (114/120) 21.7% (26/120) 10:1 10.0 ± 0.1 10.0 ± 0.1 39.7 ± 3.1 48.2 ± 2.5 82.5 ± 26.2 29.9 ± 13.9 100% (120/120) 81.7% (98/120) 15:1 — 15.0 ± 0.1 — 44.0 ± 2.5 — 79.1 ± 22.5 — 100% (120/120)

PSNR: peak signal-to-noise ratio; HDR-VDP: high-dynamic range visual difference predictor.

Data are mean ± standard deviations.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Human Observer Analysis

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

VLT

Get Radiology Tree app to read full this article<

Peak Signal-to-Noise Ratio

Get Radiology Tree app to read full this article<

PSNR=20log10(255RMSE), P

S

N

R

=

20

log

10

(

255

R

M

S

E

)

,

where

RMSE(root - mean - square error)=∑512x=1∑512y=1(f(x,y)−g(x,y))25122−−−−−−−−−−−−−−−−−√, R

M

S

E

(

root - mean - square error

)

=

x

=

1

512

y

=

1

512

(

f

(

x

,

y

)

g

(

x

,

y

)

)

2

512

2

,

where f(x, y) and g(x, y) are the pixel values in the original and compressed images, respectively.

Get Radiology Tree app to read full this article<

Perceptual Model

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Statistical Analysis

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Results

Get Radiology Tree app to read full this article<

Figure 1, Individual readers’ responses at each compression level for the thin (a) and thick (b) sections. Each bar indicates the percentage of positive responses (compressed images being rated as distinguishable from their originals). Error bars indicate 95% confidence intervals. Each gray shade represents a different reader.

Figure 2, JPEG2000 compression artifacts in (a) thin- and (b) thick-section abdomen computed tomography images in a 38-year-old male with colonic diverticulitis (arrows). According to the pooled readers’ responses, the visually lossless threshold range was 4:1–6:1 and 8:1–10:1 for the thin and thick sections, respectively. The compression artifacts are best demonstrated if the original and compressed images are downloaded (supplementary materials) and displayed alternately on the same monitor. Subtraction images (second columns) and high-dynamic range visual difference predictor (HDR-VDP) maps (right columns) represent the mathematical and predicted perceptual differences, respectively, between the original and compressed images at each compression level. The regions of interest for the original and compressed images are smaller than those for the subtraction images and the HDR-VDP maps. For the original and compressed images, window width and level are 400 and 20 Hounsfield units, respectively.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Figure 3, Correlation between the metric results and the number of positively responding readers. Symbols □, ○, ▵, +, and × represent 4:1, 6:1, 8:1, 10:1, and 15:1 compressions, respectively. The solid and dotted horizontal lines represent the cutoff values yielding the maximum sum of sensitivity and specificity and yielding 100% sensitivity, respectively, in the receiver operating characteristic analyses. (a) Peak signal-to-noise ratio (PSNR) ( r = −0.853) and (b) high-dynamic range visual difference predictor (HDR-VDP) results ( r = 0.909) for subset A (120 compressed thin sections). (c) PSNR ( r = −0.874) and (d) HDR-VDP results ( r = 0.904) for subset B (120 compressed thick sections). (e) PSNR ( r = −0.843) and (f) HDR-VDP ( r = 0.902) results for subset C (60 compressed thin and 60 compressed thick sections).

Table 3

Results of ROC Analysis

PSNR HDR-VDP_P_ Value Subset A (120 thin sections) AUC 0.98 (0.94–0.99) 0.99 (0.95, 0.99) .35 Cutoff value yielding the maximum sum of sensitivity and specificity 45.6 dB 24.8 Sensitivity (%) 95.6 (87.6–99.0) 94.1 (85.6–98.3) Specificity (%) 92.3 (81.4–97.8) 100 (93.1–100) Cutoff value yielding 100% sensitivity 47.0 dB 14.8 Specificity (%) 76.9 (63.2–87.5) 80.8 (67.5–90.4) Subset B (120 thick sections) AUC 0.97 (0.93–0.99) 0.98 (0.93–0.99) .53 Cutoff value yielding the maximum sum of sensitivity and specificity 49.7 dB 18.0 Sensitivity (%) 90.2 (79.8–96.3) 91.8 (81.9–97.3) Specificity (%) 96.6 (88.3–99.5) 100 (93.9–100) Cutoff value yielding 100% sensitivity 52.3 dB 6.5 Specificity (%) 57.6 (44.1–70.4) 54.2 (40.8–67.3) Subset C (60 thin and 60 thick sections) AUC 0.93 (0.87–0.97) 0.97 (0.92–0.99) .03 Cutoff value yielding the maximum sum of sensitivity and specificity 48.4 dB 20.3 Sensitivity (%) 89.4 (79.4–95.6) 89.3 (79.4–95.6) Specificity (%) 83.3 (70.7–92.1) 100 (93.3–100) Cutoff value yielding 100% sensitivity 52.3 dB 6.5 Specificity (%) 53.7 (39.6–67.4) 59.3 (45.0–72.4)

PSNR: peak signal-to-noise ratio; HDR-VDP: high-dynamic range visual difference predictor; ROC: receiver operating characteristic; AUC: area under the curve.

Data in parentheses are the 95% confidence intervals.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Figure 4, VLT estimation results by peak signal-to-noise ratio (PSNR) and high-dynamic range visual difference predictor (HDR-VDP) using the cutoff values that yielded the maximum sum of sensitivity and specificity in subset C. Solid line bubbles indicate that the visually lossless threshold (VLT) predicted by the metric ( y coordinate) matched the pooled readers’ decision ( x coordinate), which was regarded as a reference standard. Dashed line bubbles indicate over- or underestimations of the VLT. The bubble area is proportional to the number of superimposed data points. (a) PSNR and (b) HDR-VDP results for the 120 thin sections. (c) PSNR and (d) HDR-VDP results for the 120 thick sections.

Figure 5, Visually lossless threshold (VLT) estimation results by peak signal-to-noise ratio (PSNR) and high-dynamic range visual difference predictor (HDR-VDP) using the cutoff values that yielded 100% sensitivity in subset C. Solid line bubbles indicate that the VLT predicted by the metric ( y coordinate) matched the pooled readers’ decision ( x coordinate), which was regarded as a reference standard. Dashed line bubbles indicate over- or underestimations of the VLT. The bubble area is proportional to the number of superimposed data points. (a) PSNR and (b) HDR-VDP results for the 120 thin sections. (c) PSNR, and (d) HDR-VDP results for the 120 thick sections.

Table 4

VLT Range Prediction for 120 Thin and 120 Thick Sections Using the Cutoff Values Determined from ROC Analyses for Subset C

Concordant Underestimated Overestimated Cutoff value yielding the maximum sum of sensitivity and specificity PSNR 124 (57, 67) 65 (61, 4) 51 (2, 49) HDR-VDP 183 (94, 89) 20 (16, 4) 37 (10, 27) Cutoff value yielding 100% sensitivity PSNR 40 (11, 29) 199 (109, 90) 1 (0, 1) HDR-VDP 60 (36, 24) 180 (84, 96) 0 (0, 0)

See Table 3 for abbreviations.

Data are the number of images. Data in parentheses (separated by a comma) are the numbers of thin and thick sections.

Get Radiology Tree app to read full this article<

Discussion

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Acknowledgments

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

References

  • 1. Erickson B.J., Manduca A., Palisson P., et. al.: Wavelet compression of medical images. Radiology 1998; 206: pp. 599-607.

  • 2. Schomer D.F., Elekes A.A., Hazle J.D., et. al.: Introduction to wavelet-based compression of medical images. Radiographics 1998; 18: pp. 469-481.

  • 3. Bak P.R.G.: Will the use of irreversible compression become a standard of practice?. http://www.scarnet.org/pdf/SCARNews_Winter06.pdf Accessed July 28, 2006

  • 4. Lee K.H., Lee H.J., Kim J.H., et. al.: Managing the CT data explosion: initial experiences of archiving volumetric datasets in a mini-PACS. J Digit Imaging 2005; 18: pp. 188-195.

  • 5. Rubin G.D.: Data explosion: the challenge of multidetector-row CT. Eur J Radiol 2000; 36: pp. 74-80.

  • 6. Tamm E.P., Thompson S., Venable S.L., McEnery K.: Impact of multislice CT on PACS resources. J Digit Imaging 2002; 15: pp. 96-101.

  • 7. Cosman P.C., Davidson H.C., Bergin C.J., et. al.: Thoracic CT images: effect of lossy image compression on diagnostic accuracy. Radiology 1994; 190: pp. 517-524.

  • 8. Goldberg M.A., Gazelle G.S., Boland G.W., et. al.: Focal hepatic lesions: effect of three-dimensional wavelet compression on detection at CT. Radiology 1997; 202: pp. 159-165.

  • 9. Ko J.P., Chang J., Bomsztyk E., et. al.: Effect of CT image compression on computer-assisted lung nodule volume measurement. Radiology 2005; 237: pp. 83-88.

  • 10. Ko J.P., Rusinek H., Naidich D.P., et. al.: Wavelet compression of low-dose chest CT data: effect on lung nodule detection. Radiology 2003; 228: pp. 70-75.

  • 11. Li F., Sone S., Takashima S., et. al.: Effects of JPEG and wavelet compression of spiral low-dose CT images on detection of small lung cancers. Acta Radiol 2001; 42: pp. 156-160.

  • 12. Megibow A.J., Rusinek H., Lisi V., et. al.: Computed tomography diagnosis utilizing compressed image data: an ROC analysis using acute appendicitis as a model. J Digit Imaging 2002; 15: pp. 84-90.

  • 13. Ohgiya Y., Gokan T., Nobusawa H., et. al.: Acute cerebral infarction: effect of JPEG compression on detection at CT. Radiology 2003; 227: pp. 124-127.

  • 14. Zalis M.E., Hahn P.F., Arellano R.S., et. al.: CT colonography with teleradiology: effect of lossy wavelet compression on polyp detection-initial observations. Radiology 2001; 220: pp. 387-392.

  • 15. Zheng L.M., Sone S., Itani Y., et. al.: Effect of CT digital image compression on detection of coronary artery calcification. Acta Radiol 2000; 41: pp. 116-121.

  • 16. Kim TJ, Lee KW, Kim B, et al. Regional variance of visually lossless threshold in compressed chest CT Images: lung versus mediastinum and chest Wall. Eur J Radiol. In press.

  • 17. Woo H.S., Kim K.J., Kim T.J., et. al.: JPEG 2000 compression of abdominal CT: difference in compression tolerance between thin- and thick-section images. Am J Roentgenol 2007; 189: pp. 535-541.

  • 18. Slone R.M., Foos D.H., Whiting B.R., et. al.: Assessment of visually lossless irreversible image compression: comparison of three methods by using an image-comparison workstation. Radiology 2000; 215: pp. 543-553.

  • 19. Ringl H., Schernthaner R.E., Bankier A.A., et. al.: JPEG2000 compression of thin-section CT images of the lung: effect of compression ratio on image quality. Radiology 2006; 240: pp. 869-877.

  • 20. Slone R.M., Muka E., Pilgram T.K.: Irreversible JPEG compression of digital chest radiographs for primary interpretation: assessment of visually lossless threshold. Radiology 2003; 228: pp. 425-429.

  • 21. Lee K.H., Kim Y.H., Kim B.H., et. al.: Irreversible JPEG 2000 compression of abdominal CT for primary interpretation: assessment of visually lossless threshold. Eur Radiol 2007; 17: pp. 1529-1534.

  • 22. Kim B, Lee KH, Kim KJ, et al. Prediction of perceptible artifacts in JPEG2000 compressed chest CT images using mathematical and perceptual quality metrics. Am J Roentgenol. In press.

  • 23. Krupinski E.A., Kallergi M.: Choosing a radiology workstation: technical and clinical considerations. Radiology 2007; 242: pp. 671-682.

  • 24. Mantiuk R., Daly S., Myszkowski K., et. al.: Predicting visible differences in high dynamic range images-model and its calibration. Proc Human Vision and Electronic Imaging X, IS&T/SPIE’s 17th Annual Symposium on Electronic Imaging 2005; pp. 204-214.

  • 25. Mahesh M., Scatarige J.C., Cooper J., et. al.: Dose and pitch relationship for a particular multislice CT scanner. Am J Roentgenol 2001; 177: pp. 1273-1275.

  • 26. Kim KJ, Kim B, Choi SW, et al. Definition of compression ratio: difference between two commercial JPEG2000 program libraries. Telemed J E Health. In press.

  • 27. Digital Imaging and Communications in Medicine (DICOM). Part 14: gray scale standard display function. NEMA standards publication no. PS 3.14-2006.2006.National Electrical Manufacturers AssociationVa

  • 28. Wang Z., Sheikh H.R., Bovik A.C.: Objective video quality assessment.Furht B.Marqure O.The handbook of video databases: design and applications.2003.CRC PressBoca Raton, FL:pp. 1041-1078.

  • 29. Quick R.F.: A vector magnitude model of contrast detection. Kybernetik 1974; 16: pp. 65-67.

  • 30. Eckert M.P., Bradley A.P.: Perceptual quality metrics applied to still image compression. Signal Processing 1998; 70: pp. 177-200.

  • 31. Fleiss J.L., Cuzick J.: The reliability of dichotomous judgements: unequal numbers of judges per subject. Appl Psychol Measurement 1979; 3: pp. 537-542.

  • 32. Persons K., Palisson P., Manduca A., et. al.: An analytical look at the effects of compression on medical images. J Digit Imaging 1997; 10: pp. 60-66.

  • 33. Daly S.: The visible differences predictor: an algorithm for the assessment of image fidelity.Watson A.B.Digital images and human vision.1993.MIT PressCambridge, MA:pp. 179-206.

  • 34. Lubin J.: The use of psychophysical data and models in the analysis of display system performance.Watson A.B.Digital images and human vision.1993.MIT PressCambridge, MA:pp. 163-178.

  • 35. Siddiqui K.M., Siegel E.L., Reiner B.I., et. al.: Correlation of radiologists’ image quality perception with quantitative assessment parameters: just-noticeable difference vs. peak signal-to-noise ratios. Proc SPIE 2005; 5748: pp. 58-64.

  • 36. Krupinski E.A., Johnson J., Roehrig H., et. al.: Using a human visual system model to optimize soft-copy mammography display. Acad Radiol 2003; 10: pp. 1030-1035.

  • 37. Li B., Meyer G.W., Klassen R.V.: A comparison of two image quality models. Proc SPIE 1998; 3299: pp. 98-109.

  • 38. Lubin J., Brill M.H., Crane R.L.: Vision model-based assessment of distortion magnitudes in digital video. http://www.mpeg.org/MPEG/JND/ Accessed July 4, 2006

  • 39. Daly S.: A visual model for optimizing the design of image processing algorithms. IEEE Int Conf Image Proc 1994; pp. 16-20.

  • 40. Siddiqui K.M., Johnson J.P., Reiner B.I., et. al.: Discrete cosine transform JPEG compression vs. 2D JPEG2000 compression: JNDmetrix visual discrimination model image quality analysis. Proc SPIE 2005; 5748: pp. 202-207.

  • 41. Mantiuk R.: HDR visual difference predictor. http://sourceforge.net/projects/hdrvdp Accessed July 4, 2006

  • 42. Eskicioglu A.M., Fisher P.S.: Image quality measures and their performance. IEEE Trans Comm 1995; 43: pp. 2959-2965.

  • 43. Final report from the video quality experts group on the validation of objective models of video quality assessment, phase II. http://www.its.bldrdoc.gov/vqeg/projects/frtv_phaseII/downloads/VQEGII_Final_Report.pdf Accessed July 4, 2006

  • 44. Wang Z., Bovik A.C., Sheikh H.R., et. al.: Image quality assessment: from error visibility to structural similarity. IEEE Trans Image Proc 2004; 13: pp. 600-612.

  • 45. Lee K.H., Kim Y.H., Hahn S., et. al.: Computed tomography diagnosis of acute appendicitis: advantages of reviewing thin-section datasets using sliding slab average intensity projection technique. Invest Radiol 2006; 41: pp. 579-585.

  • 46. Prokop M.: Multislice CT: technical principles and future trends. Eur Radiol 2003; 13: pp. M3-M13.

This post is licensed under CC BY 4.0 by the author.