Converting a digital radiograph to an absolute thickness map

• US 20060058974A1
• Filed: 04/18/2005
• Published: 03/16/2006
• Est. Priority Date: 09/15/2004
• Status: Active Grant

• Current Assignee: BL Technologies Inc. (Ashland Global Holdings, Inc.)
• Original Assignee: GE Betz, Inc. (General Electric Co. (Italy))
• Inventors: Venugopal, Manoharan, Mishra, Debasish, Lasiuk, Brian W., Gopalika, Nityanand, Venkatachalam, Rajashekar, Ramaswamy, Sivaramanivas, Hopple, Michael R., Batzinger, Thomas J., Dixon, Elizabeth L., Reddy Kommareddy, Vamshi Krishna
• Granted Patent: US 7,480,363 B2
• Time in Patent Office: Days
• US Class Current: 702/97
• CPC Class Codes: G01B 15/025 by measuring absorption

Abstract

A digital radiography imaging system for acquiring digital images of an object, and a method for transforming digital images into an absolute thickness map characterizing the object under inspection. The system includes a radiation source for directing radiation through a desired region of the object, and a radiation detector having a plurality of sensing elements for detecting radiation passing through the object. Numerical data generated from each sensing element is calibrated, for example by correcting for variations in radiation paths between the source and detector, by correcting for variations in the spatial frequency response (MTF) of the detector, by correcting for variations in the geometric profile of the object under inspection, and by correcting for material contained in and/or around the object. The calibrated data is processed in order to generate and display an absolute thickness map of the object. The calibration procedures are adapted for extracting a thickness map from both isotope sources and X-ray tube sources.

Claims

1. A method for calibrating digital data representing a digital image of an object, said method comprising the steps of:

a. providing a radiation source for directing radiation through a region of said object;

b. providing a radiation detector having a plurality of sensing elements for detecting said radiation;

c. acquiring a flat field image to equalize a response from each said sensing element;

d. recording numerical data indicative of an intensity of said radiation incident on each said sensing element;

e. correlating said numerical data to a measurable property of said object;

f. calibrating said numerical data, comprising the steps of;

i. correcting for variations in a source-detector-distance (SDD) between said source and said detector;

ii. correcting for variations in a source-object-distance (SOD) between said source and said object;

iii. correcting for variations in a spatial frequency response (MTF) of said detector; and

iv. correcting for variations in a geometric profile defining said object; and

g. utilizing said calibrated data to generate a map representing said measurable property.

2. The method of claim 1, wherein said calibrating step (f) further comprises the step of correcting for material contained within said object.

3. The method of claim 2, wherein said calibrating step (f) further comprises the step of correcting for material attached to an outside surface of said object.

4. The method of claim 3, wherein said object is a pipeline and said measurable parameter is a wall thickness of said pipeline.

5. The method of claim 1, wherein said calibrating step (f) further comprises the step of correcting for variations in EM scattering and undercutting of said radiation.

6. The method of claim 1, wherein said detector is a linear array detector.

7. The method of claim 6, wherein said detector is a flat panel digital detector.

8. The method of claim 7, further comprising the step of parameterizing absorption measurements with a multi-component exponential function.

9. A system for calibrating digital images, comprising:

a. a radiation source for directing radiation through a region of an object;

b. a radiation detector having a plurality of sensing elements for detecting said radiation;

c. means for acquiring a flat field image to equalize a response from each said sensing element;

d. means for recording numerical data indicative of an intensity of said radiation incident on each said sensing element;

e. means for correlating said numerical data to a measurable property of said object;

f. means for calibrating said numerical data, comprising;

i. means for correcting for variations in a source-detector-distance (SDD) between said source and said detector;

ii. means for correcting for variations in a source-object-distance (SOD) between said source and said object;

iii. means for correcting for variations in a spatial frequency response (MTF) of said detector;

iv. means for correcting for variations in a geometric profile defining said object; and

g. means for generating a map representing said measurable property from said calibrated data.

10. The imaging system of claim 9, wherein said means for calibrating further comprises means for correcting for variations in EM scattering and undercutting of said radiation.

11. The imaging system of claim 10, wherein said means for calibrating further comprises means for correcting for material contained within said object.

12. The imaging system of claim 11, wherein said means for calibrating further comprises means for correcting for material attached to an outside surface of said object.

13. The imaging system of claim 12, wherein said object is a pipeline, and said measurable parameter is a wall thickness of said pipeline.

14. A method for calibrating digital data, comprising the steps of:

a. providing a calibration block having a plurality of thickness regions disposed along a length of said block;

b. acquiring a digital image representing each said thickness region;

c. detecting the boundaries of each said thickness region;

d. sampling an inner area of said image at each said thickness region so as to generate a mean pixel value and an uncertainty representing each said thickness region;

e. mapping each said mean pixel value to a thickness dimension defining each said thickness region; and

f. generating a look-up table to store results from said mapping step (e).

15. The method of claim 14, wherein said detecting step (c) further comprises the steps of:

calculating a center-line-path proximate a mid-width of said image, said center-line-path being substantially parallel to a length of said image;

extracting pixel values from said image along said center-line-path to obtain a first reference distribution characterizing said thickness regions;

differentiating said first reference distribution to identify boundaries between each said thickness region; and

determining the number of said thickness regions and the outer dimensions defining said thickness regions.

16. The method of claim 15, wherein said detecting step (c) further comprises the steps of:

i. placing an attenuation marker proximate a midsection of one of said thickness regions;

ii. calculating an edge-line-path proximate an edge of said image, said edge-line-path being substantially parallel to said center-line-path;

iii. extracting pixel values from said image along said edge-line-path to obtain a second reference distribution characterizing said thickness regions;

iv. differentiating said second reference distribution;

v. comparing said first and second differentiated reference distributions to identify presence or non-presence of said attenuating marker in said second reference distribution;

vi. re-calculating said edge-line-path incrementally closer to said center-line-path if said comparing step (v) fails to identify presence of said attenuation marker; and

vii. repeating steps (iii) through (vi) until said comparing step (v) identifies presence of said attenuation marker.

17. The method of claim 16, wherein said detecting step (c) further comprises the use of edge filters and non-linear transforms.

18. The method of claim 17, wherein said non-linear transforms are Hough transforms.

19. The method of claim 18, wherein said sampling step (d) further comprises the step of filtering said reference data with a high order non-linear filter.

20. The method of claim 19, further comprising the step of masking the outer edges of said block with radiation attenuating material to reduce the amount of radiation scatter and undercutting around said outer edges.

21. The method of claim 20 wherein said inner area comprises an inner central about 25% of the surface area of said thickness regions.

22. The method of claim 21, wherein said mapping step (e) further comprises the step of fitting said mean pixel values with a chi-square or least-squares fitting algorithm.

23. The method of claim 22, further comprising the step of scaling said digital data by a linear factor calculated from an exposure ratio.

24. The method of claim 23, further comprising the step of scaling said data as a function of the energy of said X-ray source or as a function of the detector response at predetermined voltage levels.