The present invention is directed to the field of X-ray detection systems.
There exists a need for improved systems and methods of screening baggage for explosives, weapons, and other contraband. Some existing systems employ X-ray scanners, computed tomography (CT) scanners, or other imaging devices to detect concealed objects. In some such systems, a CT scanner is preceded by an X-ray scanner, which performs a “prescanning” function to determine initial information on the contents of an article of baggage. Existing X-ray based systems provide differing degrees of sophistication in terms of their ability to analyze baggage based on the X-ray data obtained. Some, for example, balance the speed of the baggage screening with the accuracy and reliability with which contraband is detected. While the prescanning function discussed above may increase the accuracy and reliability with which contraband is detected, there exists a need for improved systems and methods of screening baggage.
One embodiment of the invention is directed to a method or apparatus for analyzing an object in which a dual energy X-ray prescanner performs a prescan of the object to determine prescan information about the object. Then, a CT scanner performs a CT scan on at least one plane of the object based on the prescan information. If the CT scan of the object includes or is in the vicinity of metal, then metal artifact correction of a reconstructed image from the CT scan may be performed using the prescan and CT scan information.
Another embodiment of the invention is directed to a method or apparatus for analyzing an object in which a prescanner, which need not be a dual energy prescanner, performs a prescan of the object to determine prescan information. Then, a CT scanner performs a CT scan of the object to determine CT information. A processor analyzes the CT information and the prescan information to determine whether to update the prescan information based on the CT information.
While the description and claims herein recite use of a CT scanner, such term is intended to cover any device that measures at least density of an object scanned by the device.
The present invention relates to a system or method in which a prescanner X-ray device and a downstream (of the prescanner) computed tomography (CT) device scan an object. The object may be located within a piece of baggage, a manufactured product, the human body, or some other item penetrable by X-rays. Information collected on the object may be transmitted from the prescanner to the CT scanner and/or from the CT scanner to the prescanner.
One embodiment of the present invention, illustrated in
In accordance with one illustrative embodiment, information from prescanner device 1 is transmitted from prescanner device 1 to a processor 5 via a data link 7. Data link 7, and any other data link described herein, is not limited to any particular type of link and may be implemented using any suitable means for transmitting information, such as an Ethernet link.
Processor 5 may process the information transmitted from the prescanner device, and transmit the processed information, or a control signal with instructions based on the processed information, to CT scanner device 3 via a data link 9. Processor 5 may be located external or internal to CT scanner device.
It should be appreciated that while
Prescanner device 1 may be any of numerous multiple energy X-ray devices. For example, prescanner device 1 may be a single or multi-view dual energy line scanning X-ray device, a dual energy CT scanner device, or any other device capable of measuring effective atomic number characteristics of an object, the significance of which will be appreciated from the forthcoming discussion. U.S. Pat. No. 5,838,758 (Krug), which is hereby incorporated by reference, teaches dual energy X-ray inspection systems, any of which may be employed as the prescanner device according to an embodiment of the invention.
CT scanner device 3 may be any of numerous devices for performing computed tomography or, more generally, may be any device capable of measuring density characteristics of an object. Prescanner device 1 and CT scanner device 3 may be implemented as separate units, as shown in
The screening systems described herein may be used in a variety of applications to recognize and detect target objects of interest. Target objects may include, but are not limited to, concealed objects (e.g., explosive devices or other weapons) inside a container (e.g., baggage), defects (e.g., cracks, air bubbles, or impurities) in articles of manufacture (e.g., commercial products), and areas of interest (e.g., tumors or other masses, including masses located near bone, metal, or another high-density material from which artifacts may result) within the body. Thus, the invention described herein may be used, for example, in settings such as airports, manufacturing plants, and hospitals, and other settings in the travel, commercial, and medical industries.
Certain characteristics of target objects discussed above can be determined mathematically based on the absorption of X-ray radiation by the object. The absorption of X-ray radiation by a material in an item is proportional to the degree of X-ray attenuation and is dependent on the energy of the X-ray radiation and the following material parameters: thickness, density, and atomic number. The relationship between these values can be described by Equation 1:
Ix=I0 exp [−(μ/ρ)x] (1)
where, Ix is the intensity of the X-ray radiation after passing through a material, I0 is the intensity of the X-ray radiation before passing through a material, μ/ρ is the mass attenuation coefficient; and x is obtained by multiplying the thickness of the material by its density. It should be appreciated that since X-ray absorption by a material is dependent on the thickness, density, and atomic number of the material, absorption and attenuation may be most accurately determined when all three parameters of a material are known. The scanning devices described herein can accurately determine the thickness, density, and/or atomic number of an object, and these parameters may be used to determine whether an object is a target object.
In the embodiment of
This method reduces the number of slices necessary to be taken by the CT scanner, including the number of slices taken through metal, to detect a target object and increases the accuracy with which target objects are detected. A CT scanner device employed alone to scan an item performs CT scans of planes (or “slices”) of the item and provides information on the three dimensional spatial configurations of objects therein. While this technique is useful in identifying target objects within the scanned item, each CT scan is time consuming and has a limited image quality. Numerous of these time-consuming scans are required to ensure no target area is missed. By employing prescanner device 1 upstream of the CT scanner, according to one embodiment of the present invention, possible target objects and their two-dimensional locations are determined in a quick (relative to a CT scan) prescan. A significant advantage lies in reducing the number of slices, and thereby reducing the scan time, for an item.
In addition to reducing the scan time of the CT scanner device, the feeding forward of information from prescanner device 1 to CT scanner device 3 may increase the accuracy of the CT scan images. For example, as will be described in greater detail below, for those slices that are in the vicinity of metal, the fedforward information can be used to perform metal artifact correction, thereby increasing the accuracy of any reconstructed image from the CT scan and ability to detect target objects.
Another embodiment of the present invention, illustrated in
According to one embodiment of the invention (
Another embodiment of the present invention, illustrated in
In the embodiment of
Beginning with step 20, an item (e.g., an article of baggage) to be screened is loaded into a machine of the invention. In step 21, the item is scanned and analyzed using the prescanner device 1. The prescanner device 1 may be a line scanner, such as one of the VIS series offered by PerkinElmer Detection Systems, the assignee herein. The item is initially loaded into the prescanner device 1 for scanning. For example, a human operator may place the item on a conveyor which, with the aid of a motion controller, moves the item through prescanner device 1. In one embodiment, prescanner device 1 has at least two X-ray sources for generating X-ray beams and may have one or more X-ray detectors for receiving X-ray beams. The X-ray image resulting from the scan consists of a two-dimensional array of pixels representing a view of the three-dimensional item from one angle. A processor, either internal or external to prescanner device 1, calculates the attenuation of the generated X-rays penetrating the item for each pixel. According to one embodiment of the invention, alternate pulses of high energy X-rays (e.g., 150 kV) and low energy X-rays (e.g., 75 kV) are respectively generated by dual X-ray sources, and the processor calculates the attenuation for each pixel of the image resulting from the respective high energy and low energy beams.
In a step 23, a table (Table A) is generated containing atomic number and mass characteristics for each object. Table A may be stored electronically by a memory (not shown) coupled to a processor. Both the processor and the memory may be either internal or external to prescanner device 1. An object may be defined as any region having similar atomic number and mass characteristics. The calculated attenuation of the high energy and low energy beam pulses for each pixel of the scanned item are used to determine the effective atomic number of all objects. To derive the effective atomic number of each object based on the attenuation, the attenuation of X-rays at each different energy level is analyzed. One method for doing so is described in U.S. Pat. No. 5,838,758 (Krug), incorporated by reference herein. It is known that materials with a high effective atomic number (e.g., metals) absorb low energy X-ray radiation more strongly, whereas materials with a low effective atomic number (e.g., organic materials) absorb high energy X-ray radiation more strongly. Thus, the effective atomic number of each object may be determined by analyzing the attenuation of low and high energy X-rays by each pixel. To determine the effective atomic number for a particular object, all pixels within the object are compared to pixels surrounding the object and a histogram is created, where the mode (peak of the histogram) represents the effective atomic number.
In addition to effective atomic number information, Table A may also contain mass information for each object. The mass for each pixel may also be determined based on the X-ray attenuation of both the high and low energy X-rays. The relationship between X-ray attenuation and material mass (i.e., thickness) is logarithmic; X-ray radiation decreases logarithmically as the material thickness increases. Thus, mass may be estimated by analyzing the attenuation of X-rays of all energies by materials within an item. To determine the mass for a particular object, mass values for all pixels within an object are added.
In an embodiment, Table A also contains confidence values for the effective atomic number and mass values for each object. Confidence values for the effective atomic number and mass values represent a probability or range of probabilities that the atomic number and mass data are correct. To determine a confidence level for the effective atomic number value or mass value of a particular object, a feature vector denoting properties such as compactness, connectiveness, gradients, histogram spread and other features may be used.
Numerous known procedures are available for determining the confidence level. One such procedure uses machine vision technology for object classification. Machine vision technology includes: (1) segmenting a group of picture elements from their background, (2) describing that group of picture elements by a set of features, and (3) using the resulting feature vector to classify the picture elements.
One software tool available for such object classification is Image Process and Analysis Software offered by Data Translation, Inc. as SP0550. Other software packages that provide similar tools for algorithm development include: Checkpoint® by Cognex Corporation of Natick, Mass., Framework® by DVT of Woodcliff Lake, N.J., and the Powervision® family of products of RVSI of Canton, Mass. The invention need not be limited to the features found in the exemplary software packages mentioned. There are numerous other approaches as described, for examples, in the following textbooks:
A target object, such as an explosive, has a typical effective atomic number and mass value. Further, for a particular range of atomic number values, a particular range of mass values will be characteristic of a target object. Thus, it is useful to consider both atomic number and mass values in determining whether a target object is present.
In a step 25, a list of objects warranting further study (i.e., objects of interest), including the locations for the objects, is generated. The atomic number characteristics of Table A can be used to differentiate potential target objects from the background, since different objects will generally have different effective atomic numbers. A potential target object may comprise a collection of pixels in close proximity having atomic number values that fall within a certain range. For example, a weapon or explosive may comprise a collection of pixels having high effective atomic number values that fall within a particular range. Thus, it is possible to determine two-dimensional coordinates (e.g., X1-X2, z1-z2 in
While the list of objects warranting further study and two-dimensional coordinates associated with each object may be generated automatically, it is also possible that a human operator may manually determine the information. For example, an operator may view an X-ray image to determine objects of interest and their respective locations in two dimensions. Thus, the prescan analysis may be performed automatically or manually, and the invention is not limited to either method of analysis.
Once a location of an object of interest, or a region thereof, has been determined, a CT scan of the object or region of interest may be performed. Locations of slices (i.e., two-dimensional planes) in the item to be scanned are chosen to coincide with a potential target object. Some target objects, such as explosives, are typically found near metal objects (e.g., wires, batteries). Metal, due in part to its high density, may cause artifacts in an image in the region surrounding the metal. Thus, if a potential target object is located near metal, it is preferable to choose a slice that includes the target object, but that is not in the vicinity of the metal. However, if a slice near metal is chosen, according to one aspect of the invention, a metal artifact correction is performed to correct for the image artifacts, as will be described in step 35.
If, after step 25, there are no objects warranting further study, a decision may be made as to an appropriate course of action, based on the prescan information (
In step 31, CT images are generated for the item cross-sections identified in step 25, if any. To form a CT image of a cross-section (i.e., slice) of an item, a finely collimated beam of radiation is passed through the item in the desired slice plane, and the attenuation is measured. The process is repeated and a set of projections is acquired as the X-ray beam is passed through the object at different angles. A reconstructed image of the two-dimensional distribution of the linear attenuation coefficient, μ(x,y), may be obtained from these projections. If the projections could be acquired with an infinitely narrow X-ray beam, and the angular increment at which the X-ray beam is passed was negligible, the result would be a continuous set of projections. Displayed as a two-dimensional function, the continuous set of projections is referred to as the sinogram. An image may be reconstructed from the sinogram by implementing any of a number of well-known reconstruction techniques including, but not limited to, back projection, iteration, Fourier transform, and filtered back projection.
As discussed above, a CT image of a slice results in a two-dimensional image of a cross-sectional plane of the scanned item. The image consists of an array of pixels (e.g., 900 pixels×512 pixels). According to one illustrative embodiment shown in
In step 33, it is determined whether any imaged object of interest is in the vicinity of a metal object. Additionally, it may be determined whether the image of the object of interest is likely to be distorted by metal artifacts caused by the metal object. For example, although a metal object is in close proximity to the object of interest, it may be determined that the size of the metal object relative to the object of interest renders it unlikely that the metal object will have a significant negative effect on the image of the object of interest (e.g., if the metal object is much smaller than the object of interest). If a potential target object is in the vicinity of a metal object, such that the image of the object is likely to be distorted by metal artifacts, a metal artifact correction is performed on the slice containing the metal artifacts, according to one aspect (feedforward mode) of the invention described herein.
If it is determined in step 33 that a potential target object is in the vicinity of a metal object, information fed forward from the prescanner device is used to predict the type and shape of metal responsible for the metal artifacts in step 35. In particular, the mass information and effective atomic number information from Table A are used to identify the metal type and perform a metal artifact correction specific to the type and shape of the metal. The metal artifact correction algorithm is described in detail below in connection with
In step 37, the scanned CT images are analyzed. According to one embodiment, the density, area, and three-dimensional coordinates are determined for each target object, for example using image processing algorithms (e.g., region growing). The area of each target object is specified by a range of two-dimensional (e.g., x1-x2, y1-y2 in
In step 39, a table (Table B) is generated containing the density, area, and three-dimensional coordinates for each target object, and a confidence level for each characteristic of each target object. Table B may be stored electronically by a memory (not shown) coupled to a processor and may, along with or separate from the processor, be either internal or external to CT scanner device 3. The three-dimensional coordinates for each target object are transmitted (“fed back”) to prescanner device 1 in step 41. According to one aspect (feedback mode) of the invention, this information from Table B may be used to augment Table A. The processor, coupled to the memory that stores Table A, considers the fedback information and the information in Table A in determining whether to update any of the information in Table A.
Since prescanner device 1 images the item from only one view, the prescanner device may not be able to discern whether an identified object is a single object or a plurality of objects, as objects that overlap when imaged from a particular perspective may appear as a single merged object. If the prescanner device cannot differentiate a plurality of overlapping objects, it may determine a mass value for an object that is actually the mass values of two or more objects combined. The three-dimensional coordinate information provided in Table B can be used to differentiate objects, and thereby correct erroneous effective atomic number values and mass values of Table A. If the mass of an object changes, the object may no longer be of interest or, conversely, may become interesting. For example, if an original mass determination is based on two merged non-target objects, the mass value will be erroneously high, and may fall within the range corresponding to a target object. When the two merged objects are differentiated and their masses are determined separately, the individual objects may no longer be of interest if the mass value falls below a minimum mass associated with potential target objects.
In sum, the information fedback from Table B by the CT scanner device allows for more accurate determinations of the effective atomic number and mass of each object, as listed in Table A, by the prescanner device. Hence, superior detection by the prescanner device and a lower false alarm rate may be achieved by feeding back information from the CT scanner device to the prescanner device. It should be appreciated that multiple feedforward/feedbackwards loops are possible, whereby information generated by the prescanner device 1 and CT scanner device 3 is alternately transmitted between the two devices. It should be appreciated that the information from Table B need not be transmitted to the prescanner device. Rather, the CT scanner device or an external computer may implement an algorithm, similar to that which may be implemented by the prescanner, to augment Table A based on the Table B information.
In step 27, a decision is made based on the information in Tables A and B as to an appropriate course of action. As discussed above, possible actions include returning the baggage to the passenger, searching the baggage by hand, or calling the bomb squad. An algorithm may be used to synthesize the information of the two tables to determine an appropriate action. For each potential target object, the algorithm may consider the effective atomic number, density, and associated confidence levels for each, as well as the thickness of the potential target object and the proximity of the potential target object to metal. Based on the information, a likelihood is determined that an identified object is a target object. The likelihood is derived from a histogram representing, for example, the probability that an object having a given effective atomic number, density, thickness, mass, and proximity to metal is a target object, and may be represented as a probability that the object is a target object or as an absolute indication that the object is/is not a target object. It should be appreciated that any of the automated decisions or actions described above may alternatively be performed by a human operator.
In step 43, a CT image is generated. Uncorrected CT images may contain metal artifacts when a scan is performed within a certain proximity to metal, which may result in inaccuracies. For example, beam hardening artifacts cause inaccuracies in the estimation of attenuation coefficients for pixels associated with x-rays that traverse highly attenuating structures. Streaky shadows or star patterns of streaks may result near high density objects in regions of pixels where essentially no attenuation information exists. Scatter artifacts may result from the dispersion of X-ray photons by the atoms within the item, and may cause noise in the CT image.
In step 45, the image is clipped so that the image contains only the metal that accounts for the artifacts of the image. The region to be clipped is identified by considering the effective atomic number information of Table A. Each pixel in the image of the metal will have an effective atomic number that falls within a range corresponding to the effective atomic number of the metal. The clipped image contains only the image of the metal, and does not contain the object of interest or artifacts.
The use of dual energy levels in the prescanner device makes it possible to determine the characteristics of the metal in the image. In step 47, the type of metal and thickness of the metal in the image are identified based on the information in Table A. In particular, the effective atomic number information of Table A is used to identify the type of metal and the mass information of Table A is used to determine the thickness of the metal.
A sinogram of the clipped image is generated in step 49. As discussed above, a single sinogram contains the information about a particular slice from all angles, with the information from each angle in its own row.
In step 51, a table (Table C) is generated that contains beam hardening, noise, and scatter correction parameters. The correction parameters are determined according to algorithms well-known in the art for compensating for beam hardening, noise, and scatter, based on the type and thickness of metal responsible for the artifacts.
In step 53, artifacts are introduced into the sinogram of the clipped image using the table (Table C) generated in step 51. In particular, the sinogram is corrupted using beam hardening and scatter effects based on the shape and type of the metal responsible for the artifacts, determined in step 47. The sinogram of the image of the metal and artifacts is reconstructed in step 55.
In step 57, the reconstructed artifact image generated in step 55 is subtracted from the sum of the original CT image generated in step 43 and the clipped image generated in step 45. The result of the image subtraction is a metal artifact corrected image 59. The image will result in a more accurate determination as to whether the object of interest represents a target object.
In an embodiment, the artifact image may also used as a map for determining whether the CT values read in the image are accurate.
Having described several embodiments of the invention in detail, various modifications and improvements will readily occur to those skilled in the art. Such modifications and improvements are intended to be within the spirit and scope of the invention. Accordingly, the foregoing description is by way of example only, and is not intended as limiting. The invention is limited only as defined by the following claims and equivalents thereto.
This application is a continuation of application Ser. No. 10/702,814, filed on Nov. 6, 2003, now U.S. Pat. No. 6,944,264 entitled “Method and Apparatus For Transmitting Information About A Target Object Between A Prescanner and A CT Scanner”; which, in turn, is a continuation of application Ser. No. 10/068,459, filed on Feb. 6, 2002, entitled “Method and Apparatus For Transmitting Information About A Target Object Between A Prescanner and A CT Scanner,” now U.S. Pat. No. 6,816,571, issued Nov. 9, 2004.
Number | Name | Date | Kind |
---|---|---|---|
4020346 | Dennis | Apr 1977 | A |
4029963 | Alvarez et al. | Jun 1977 | A |
4064440 | Roder | Dec 1977 | A |
4217641 | Naparstek | Aug 1980 | A |
4247774 | Brooks | Jan 1981 | A |
4539648 | Schatzki | Sep 1985 | A |
4580219 | Pelc et al. | Apr 1986 | A |
4590558 | Glover et al. | May 1986 | A |
4709333 | Crawford | Nov 1987 | A |
4759047 | Donges et al. | Jul 1988 | A |
4788704 | Donges et al. | Nov 1988 | A |
4941162 | Vartsky et al. | Jul 1990 | A |
4957250 | Hararat-Tehrani | Sep 1990 | A |
5070519 | Stein et al. | Dec 1991 | A |
5109691 | Corrigan et al. | May 1992 | A |
5125015 | Shimoni et al. | Jun 1992 | A |
5162652 | Cohen et al. | Nov 1992 | A |
5182764 | Peschmann et al. | Jan 1993 | A |
5243664 | Tuy | Sep 1993 | A |
5319547 | Krug et al. | Jun 1994 | A |
5323004 | Ettinger et al. | Jun 1994 | A |
5367552 | Peschmann | Nov 1994 | A |
5490218 | Krug et al. | Feb 1996 | A |
5600303 | Husseiny et al. | Feb 1997 | A |
5600700 | Krug et al. | Feb 1997 | A |
5642393 | Krug et al. | Jun 1997 | A |
5661774 | Gordon et al. | Aug 1997 | A |
5666391 | Ohnesorge et al. | Sep 1997 | A |
5796802 | Gordon | Aug 1998 | A |
5805660 | Perion et al. | Sep 1998 | A |
5838758 | Krug et al. | Nov 1998 | A |
5905809 | Timmer | May 1999 | A |
5933471 | Kalvin | Aug 1999 | A |
5953444 | Joseph et al. | Sep 1999 | A |
6018562 | Willson | Jan 2000 | A |
6026143 | Simanovsky et al. | Feb 2000 | A |
6076400 | Bechwati et al. | Jun 2000 | A |
6088423 | Krug et al. | Jul 2000 | A |
6094467 | Gayer et al. | Jul 2000 | A |
6094472 | Smith | Jul 2000 | A |
6118850 | Mayo et al. | Sep 2000 | A |
6125193 | Han | Sep 2000 | A |
6163591 | Benjamin | Dec 2000 | A |
6198795 | Naumann et al. | Mar 2001 | B1 |
6218943 | Ellenbogen | Apr 2001 | B1 |
6256404 | Gordon et al. | Jul 2001 | B1 |
6272230 | Hiraoglu et al. | Aug 2001 | B1 |
6298112 | Acharya et al. | Oct 2001 | B1 |
6345113 | Crawford et al. | Feb 2002 | B1 |
6359961 | Aufrichtig et al. | Mar 2002 | B1 |
6418189 | Schafer | Jul 2002 | B1 |
6430255 | Fenkart et al. | Aug 2002 | B2 |
6437656 | Guynn et al. | Aug 2002 | B1 |
6600801 | Raupach | Jul 2003 | B2 |
6707879 | McClelland et al. | Mar 2004 | B2 |
6721391 | McClelland et al. | Apr 2004 | B2 |
6735272 | Sorenson | May 2004 | B1 |
6788761 | Bijjani et al. | Sep 2004 | B2 |
6816571 | Bijjani et al. | Nov 2004 | B2 |
6944264 | Bijjani et al. | Sep 2005 | B2 |
7023957 | Bijjani et al. | Apr 2006 | B2 |
20020172324 | Ellenbogen et al. | Nov 2002 | A1 |
20020176531 | McClelland et al. | Nov 2002 | A1 |
20020186862 | McClelland et al. | Dec 2002 | A1 |
20030085163 | Chan et al. | May 2003 | A1 |
20040076262 | Shao et al. | Apr 2004 | A1 |
20040120456 | Ellenbogen et al. | Jun 2004 | A1 |
20050008119 | McClelland et al. | Jan 2005 | A1 |
20050031076 | McClelland et al. | Feb 2005 | A1 |
Number | Date | Country | |
---|---|---|---|
20050111619 A1 | May 2005 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 10702814 | Nov 2003 | US |
Child | 11018078 | US | |
Parent | 10068459 | Feb 2002 | US |
Child | 10702814 | US |