Converting a digital radiograph to an absolute thickness map

• US 7,480,363 B2
• Filed: 04/18/2005
• Issued: 01/20/2009
• 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: Gopalika, Nityanand, Ramaswamy, Sivaramanivas, Hopple, Michael R., Batzinger, Thomas J., Dixon, Elizabeth L., Kommareddy, Vamishi Krishna Reddy, Venugopal, Manoharan, Mishra, Debasish, Lasiuk, Brian W., Venkatachalam, Rajashekar
• Primary Examiner(s): Ho; Allen C.
• Application Number:
  US11/108,498
• Publication Number:
  US 20060058974A1
• Time in Patent Office: 1,373 Days
• Field of Search: 378/51, 378/54, 378/56, 378/57, 378/58, 378/59, 378/207, 378/55
• US Class Current: 378/54
• 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 generating a wall thickness map of an object by converting digital data representing a digital radiographic image into a display map representing the absolute thickness of the object, said method comprising the steps of:

a. providing an X-ray or gamma ray 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 measurement of the thickness 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 by correcting the flux of X-rays or gamma rays to account for scattering;

ii. correcting for variations in a source-object-distance (SOD) between said source and said object by correcting the flux of X-rays or gamma rays to account for scattering;

iii. correcting for variations in a spatial frequency response (MTF) of said detector to determine the spatial extent and penetration depth of features within the object; and

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

g. generating an absolute thickness map for a region of said object from said calibrated data.

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 thickness 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 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. A system for generating a wall thickness map of an object by converting digital data representing digital radiographic images into a display map representing the absolute thickness of a region of the object, the system comprising:

a. a radiation source for directing X-ray or gamma ray 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 measurement of the thickness 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 an absolute thickness of the region of the object from said calibrated data.

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

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

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

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

13. A method for calibrating digital data and using the data to generate a thickness map of an object, comprising the steps of:

providing an X-ray gamma ray radiation source and radiation detector;

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

acquiring a digital image with said radiation detector representing each of said thickness regions;

detecting the boundaries of each said thickness region;

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;

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

generating a look-up table to store results from said mapping step; and

generating a thickness map using the look up table.

14. The method of claim 13, wherein said detecting step 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.

15. The method of claim 14, wherein said detecting step 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.

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

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

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

19. The method of claim 18, 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.

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

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

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

23. The method of claim 22, 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.