US20090052762A1 - Multi-energy radiographic system for estimating effective atomic number using multiple ratios - Google Patents
Multi-energy radiographic system for estimating effective atomic number using multiple ratios Download PDFInfo
- Publication number
- US20090052762A1 US20090052762A1 US12/129,371 US12937108A US2009052762A1 US 20090052762 A1 US20090052762 A1 US 20090052762A1 US 12937108 A US12937108 A US 12937108A US 2009052762 A1 US2009052762 A1 US 2009052762A1
- Authority
- US
- United States
- Prior art keywords
- eff
- attenuation
- values
- radiographic images
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/06—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/217—Validation; Performance evaluation; Active pattern learning techniques
- G06F18/2178—Validation; Performance evaluation; Active pattern learning techniques based on feedback of a supervisor
Definitions
- the present invention relates generally to radiographic imaging, and, more particularly, to determination of the effective atomic number, Z eff , of a material using multi-energy radiographic imaging.
- Standard techniques using high-energy radiographic systems exist to determine or estimate atomic composition of a material. For example, two X-ray energy levels may be used to image an object of interest. The gray level intensity values measured for each of the two energy levels are used to compute a corresponding ratio. The ratio of the intensity values for an unknown material is compared against known materials. The known material with the closest ratio to that measured is used to estimate the unknown material's effective atomic number (Z eff ) of the unknown material.
- Z eff effective atomic number
- these radiographic systems may be prone to noise and other non-linear effects that can cause errors in Z eff determination, especially in high-Z materials.
- Elements with high-Z include special nuclear materials, such as plutonium and highly enriched uranium, as well as elements that would be extremely effective in shielding special nuclear materials from passive radiation detection techniques.
- insufficient penetration may occur when high-Z or high density materials do not allow enough energy to penetrate through the object of interest to a detector.
- low-Z or low density materials may induce over-saturation. In the case of over-saturation, little to no attenuation may occur at a particular X-ray energy level.
- higher energies may be used in cases of insufficient penetration. Higher energies may solve insufficient penetration issues but exacerbate over-saturation issues and vice versa.
- the solution for overcoming over-saturation issues may be at odds with the solution for overcoming insufficient penetration.
- Embodiments of the present invention may address the above-mentioned problems and limitations, among other things.
- An embodiment of the present invention can include a method of determining an effective atomic number, Z eff , of an object including obtaining a plurality of radiographic images of the object.
- Each radiographic image can be obtained using a different independent X-ray energy level.
- Each image can also include a plurality of pixels.
- the images can be registered and normalized. Noise may be removed from the images.
- An intensity value for each pixel in the region of interest in each radiographic image can be determined.
- a plurality of measured ratios, R, for each pixel in the region of interest can be formed. Each measured ratio can be formed using intensity values from corresponding pixels in a pair of different radiographic images.
- the method can include calculating at least one adjusted measured ratio, R m , for each pixel in the region of interest based on the plurality of measured ratios, R, and at least one corresponding estimation coefficient, ⁇ , and comparing the at least one adjusted measured ratio, R m , for each pixel in the region of interest to a material attenuation database.
- the method can include assigning a Z eff value to each pixel in the region of interest based on the comparison and outputting the assigned Z eff values.
- Another embodiment may include a computer program product including a computer readable medium encoded with software instructions for causing a computer to perform the steps of receiving more than two radiographic images of an object and determining an intensity value for each pixel in each of the plurality of radiographic images.
- Each image can be obtained using a different independent X-ray energy level and can have a plurality of pixels.
- the steps can also include forming a plurality of measured ratios, R.
- Each measured ratio, R can be formed using intensity values from corresponding pixels in a pair of different radiographic images.
- the steps can also include calculating at least one adjusted measured ratio, R m , based on the plurality of measured ratios, R, and at least one corresponding estimation coefficient, ⁇ , and comparing the at least one adjusted measured ratio, R m , to a material attenuation database.
- the steps can include assigning a Z eff value based on the comparison and outputting the assigned Z eff value.
- Another embodiment may include a system for determining an effective atomic number, Z eff , of an object.
- the system may include at least one processor.
- the at least one processor can form a plurality of measured ratios, R. Each measured ratio may be formed using intensity values taken from a pair of different radiographic images.
- the at least one processor can calculate at least one adjusted measured ratio, R m , based on the plurality of measured ratios, R, and at least one corresponding estimation coefficient, ⁇ .
- the at least one processor can output a Z eff value based on a comparison of the at least one adjusted measured ratio, R m , to a material attenuation database.
- FIG. 1 is a block diagram showing an overview of an embodiment of a system for estimating Z eff .
- FIG. 2 is a flowchart showing a process overview of an embodiment of a process for estimating Z eff .
- FIG. 3 is a block diagram showing an embodiment of a dynamic material attenuation processor for use in estimating Z eff .
- FIG. 4 is a flowchart showing a process overview of an embodiment of Z assignment arbitration for use in estimating Z eff .
- the present invention is directed to identifying material compositions by estimating their effective atomic number or Z eff value.
- the present invention further seeks to overcome the limitations of the above-mentioned ratio method by using a multi-energy approach. Multiple energies are used to determine intensity ratios that are proportional to attenuation coefficient ratios. These intensity ratios can then be used to build a series of measured ratios. These measured ratios can then be used to build a series of adjusted measured ratios for use in a Z eff determination that may be more accurate than conventional methods.
- Beer's law equation can be stated using an arbitrary energy level k as,
- I ( k ) I o ( k ) e ⁇ k t , (1)
- I(k) is the measured intensity of the radiation at energy level k
- I O (k) is the input intensity of the radiation at energy level k
- ⁇ k is the linear attenuation coefficient of the material at energy level k
- t is the energy independent material thickness
- the linear attenuation coefficient depends on the material cross-section coefficient, ⁇ . Since attenuation is energy dependent, each material has the basic regions of energy scatter, including photoelectric absorption ⁇ pe , coherent scatter ⁇ cs , incoherent scatter ⁇ is , and pair-wise production ⁇ pp . Each material has a material combined cross-sectional coefficient which depends on the sum of the individual cross-sections. For example,
- the cross-sections and linear attenuation coefficients are related by:
- ⁇ is the material linear attenuation coefficient
- ⁇ is the material combined cross-sectional coefficient
- N is Avogadro's number
- A is atomic weight
- ⁇ density
- R j is the ratio of the logarithm of the respective intensity ratios and x and y represent the upper and lower energy levels. This expression eliminates the dependency on mass thickness. Thus, the R j value may be compared against attenuation curves at the upper and lower energy levels for various materials to determine the Z eff of an unknown material.
- Intensity values can be determined from measured radiographic images and used in the determination of a plurality of ratios for different pairs of radiographic source energy levels. Any number of energy levels in the X-ray regime may be used. For example, at least two energy levels and corresponding measured ratios can be used. In another example, four different energy levels and corresponding measured ratios can be used.
- Table 1 shows an example of an application employing four different source energy levels.
- Four different source energy levels (E 1 -E 4 ) are used to generate four different radiographic images of an object.
- E 1 corresponds to the lowest X-ray source energy level.
- E 1 may be 1MeV.
- E 4 corresponds to the highest X-ray source energy level.
- E 4 may be 10 MeV.
- Measured intensity values I(x) for each energy level, x are extracted from the corresponding radiographic images and used in the determination of each measured ratio, R j .
- Source intensity values I o (x) for each energy level are known from measurement conditions of the radiographic images and are also used in the determination of the measured ratios, R j .
- a family of attenuation curves for a group of materials with respect to X-ray energy is typically non-linear, with various ranges of polynomial order.
- the attenuation curves themselves are non-linear, the ratios from energy level to energy level can be modeled using linear methods. At least two energy solutions form a linear attenuation fit based on the ratio. Determining the atomic structure of a material can be based on comparing the slopes of the intensity values from radiographic images with those values of known materials residing in a database.
- the attenuation coefficient is proportional to measured X-ray intensities.
- Radiographic imagers can use measured X-ray intensity to create a corresponding gray level picture, otherwise known as a radiographic image.
- the conversion from intensity to gray level can be defined by an inverse linear ratio, but the process of imaging various materials poses several non-linear settings. Despite the nonlinearity, ratios between the measured intensity values are approximately the same as ratios between the true material linear attenuation coefficients when adjusted for gray level conversion. In other words, an adjusted ratio of measured intensity values, as defined in equation (4), may be correlated with a similar ratio of known linear attenuation coefficients in determining the Z eff of an unknown material.
- R j is the measured ratio of attenuation coefficients at two different source energy levels
- ⁇ j is an estimation coefficient corresponding to the particular two energy levels used in determining R j
- j is a subscript referring to the particular pair of source energy levels.
- the estimation coefficient, ⁇ j can be determined from the mapping of the photon intensity values to a gray level value.
- Separate adjusted measured ratios, R m may be determined for each measured ratio, R, used. Alternately, the multiple separate measured ratios may be combined into a single vector. For example, using the set of ratios defined by Table 1, equation (5) would become:
- î j represents a vector in n-space corresponding to attenuation coefficient R j .
- the n-space may be a multi-dimensional space based on the number of energy levels employed or the number of measured ratios formed.
- the measured energy ratios, R using the intensity values for that particular energy level may be excluded from the calculation of the adjusted measured ratio, R m,j .
- the ratio R m,j corresponding to those conditions of over-saturation or non-penetration may be excluded from use in the determination of Z eff .
- Gray level mapping of the estimation coefficient is nonlinear across material and energy levels.
- the estimation coefficients, ⁇ may be determined empirically through experimental evaluation of known materials. Alternately, an algorithm can be used for minimum error tuning so as to optimize the estimation coefficients, ⁇ . It should be appreciated that other methods may be used to determine the estimation coefficients.
- FIG. 1 is a block diagram showing an overview of an exemplary embodiment of a system for estimating Z eff .
- a material domain imaging processor 100 may receive multiple radiographic images and other data as input and may output a multi-energy high Z-mapping for identification of threats.
- An image processor 108 may receive a set of radiographic images of an object of interest from radiographic imaging system 102 or another processor, system, or database. Each of the radiographic images may be taken at a different X-ray energy level. Each image may be processed by the image processor 108 to register the images, normalize the images, and remove any noise in the images.
- the processed images may be sent to the Z eff processor 116 .
- the Z eff processor 116 may extract corresponding intensity values from each pixel in each image.
- the Z eff processor 116 may only determine an intensity value for each pixel in each image that is within a region of interest.
- the images may be sent to an external analysis system, such as segmentation processor 114 .
- Segmentation processor 114 may use, for example, edge boundary detection, texture based detection, or other region processing in order to determine boundaries within the images and to select a common region of interest for evaluation by the material domain imaging processor 100 .
- the region of interest in each image may then be communicated to the Z eff processor 116 for determination of intensity values.
- the relevant pixels for subsequent processing may be limited to those pixels within the region of interest instead of the entirety of pixels in each image.
- Relevant pixels may also be selected based on penetration conditions (e.g., oversaturation, non-penetration). For example, like regions across the various registered images may be compared to determine penetration conditions.
- the Z eff processor 116 may use the intensities to form a set of measured ratios, R j , for the relevant pixels.
- the measured ratios, R j may be formed using the relation set forth in equation (4) above. Normalized intensity values from corresponding pixels in a pair of different radiographic images may also be used to form each measured ratio.
- Corresponding pixels may refer to pixels in different images which correspond to the same point or location on an imaged object.
- the set of ratios for each pixel may take the form shown in Table 1.
- the Z eff processor 116 may use estimation coefficients, ⁇ j , with the set of measured ratios, R j , to determine at least one adjusted measured ratio, R m,j , for each relevant pixel. Separate adjusted measured ratios, R m,j , may be determined for each energy ratio in the set or only for certain energy ratios within the set. Alternately, the multiple energy ratios may also be combined into a single measured ratio vector, R m .
- the estimation coefficients may be provided by tuning module 110 .
- the relevant pixels may include, for example, a single pixel, all the pixels in the entire image, or just the pixels within a region of interest. In an alternative embodiment, relevant pixels may be selected based on penetration conditions (e.g., over-saturation, non-penetration, etc.).
- the material domain imaging processor 100 or a separate system may compare regions across various images to determine penetration conditions.
- the material domain imaging processor 100 or a separate system i.e., segmentation processor 114 or context analysis system 106
- regions with pixel values at an upper extreme of the photon intensity scale may be correlated as regions of over-saturation.
- the Z eff processor 116 may be configured to exclude these regions from the list of relevant pixels for further processing.
- tuning module 110 can employ an algorithm for minimum error tuning of the estimation coefficients.
- Accepted attenuation database 104 may be provided to the tuning module 110 .
- accepted attenuation data may be provided by an integrated database, a memory device, or a separate system or processor.
- the tuning module 110 may use ratios from radiographic images of a known material to optimize the value of the estimation coefficients such that adjusted measured ratios correspond to ratios of accepted attenuation data.
- the accepted attenuation data may be data from public attenuation sources, such as the NIST public data source.
- the tuning module 110 may include a database of previously determined estimation coefficients for use by the Z eff processor 116 . The estimation coefficients may be determined via experimental evaluation or other means.
- the resulting adjusted measured ratios, R m,j can be output from the Z eff processor 116 to the material assignment module 118 .
- the material assignment module 118 can compare the adjusted measured ratios to a material attenuation database from dynamic material attenuation (DMA) module 112 .
- the Z eff processor 116 may send the measured ratios, R j , to the material assignment module 118 .
- material assignment module 118 would use the estimation coefficients and the measured ratios to form the adjusted measured ratios, R m,j , for comparison to the material attenuation database.
- the material assignment module 118 may assign a Z eff value to each relevant pixel or to a region of pixels based on the comparison.
- the DMA module 112 may be a processor, such as dynamic material attenuation processor 300 , that uses configurable settings to dynamically create a densely populated attenuation ratio lookup table with variable resolution in both the energy scale and the effective atomic number (Z eff ) scale.
- the DMA module 112 can be used for the two or more energies used in the radiographic images for determining Z eff values.
- Sparse attenuation data from public sources for any material (liquid, solid, gas) can be stored on a disk using a standard format.
- accepted attenuation data i.e., attenuation coefficients
- Accepted attenuation data may also be provided by an integrated database, a memory device, or a separate system or processor.
- the accepted attenuation data may be data from public attenuation sources, such as the NIST public data source.
- the attenuation data can be created using a variety of scattering approaches, such as coherent, incoherent, photoelectric, pairwise production nuclear field, pairwise production electric field, total scatter with coherent, and total scatter without coherent.
- a user may select the materials used to create the attenuation lookup table, such as water, peroxide, lead, carbon, etc.
- the attenuation lookup table may be created as a function of the X-ray energy level and the Z eff of the material.
- the DMA module 112 may include a database of previously determined attenuation lookup tables, such as material attenuation table 322 in FIG. 3 , for use by the material assignment module 118 .
- Material assignment module 118 may include a plurality of independent material assignment algorithms. Each material assignment algorithm may employ a different methodology or use different measured ratios for comparison to the same attenuation database to generate a set of candidate Z eff values for each relevant pixel. Each material assignment algorithm can also assign a confidence value to the candidate Z eff values.
- the set of candidate Z eff values may be sent to the Z eff processor 116 .
- the Z eff processor 116 may compare the confidence values for each Z eff value in the set and may select the Z eff value for each relevant pixel with the highest associated confidence value.
- the result can be a Z eff image with each pixel having an assigned Z eff value with the highest confidence.
- the Z eff processor 116 may output the result to external analysis system, such as context analysis system 106 , for evaluation of regions of non-penetration or over-saturation for inclusion in the final processor output.
- context analysis system 106 may use a-priori information, in the way of configuration data, to assist in identifying non-penetrable and false alarm cases.
- the result from the Z eff processor 116 can be sent to a color image assignment module 120 .
- the color image assignment module 120 can map a color scale to a range of corresponding Z eff values.
- each relevant pixel may be assigned a color based on the assigned Z eff value, thereby creating a color image of the object.
- This color image can be output from the material domain imaging processor 100 to threat decision processor 124 .
- the color image and Z eff values may be further processed by a threat assignment module 122 .
- the Z eff value for each pixel can be compared with a threat threshold by the threat assignment module 122 to determine regions where the threshold is exceeded, i.e., those regions where a threat exists. These regions and associated confidence values may be output from the material domain imaging processor 100 to the threat decision processor 124 for further processing or integrated decision making.
- FIG. 2 shows a process flow 200 of an exemplary embodiment of a method of material domain image processing.
- the process begins at step 201 and proceeds to step 202 .
- radiographic images may be obtained of an object.
- Each radiographic image of an object may be obtained at a different X-ray energy level.
- at least two independent X-ray energy levels can be used.
- four X-ray energy levels can be used to generate four independent radiographic images of an object.
- the radiographic images may be subject to image processing.
- the image processing may be any of a number of processing steps known in the art.
- each radiographic image may be registered with the other radiographic images such that regions of interest in the object are aligned.
- the radiographic images may be normalized to a gray scale.
- noise filtering may be employed to reduce noise artifacts in the image.
- an intensity value may be determined for each relevant pixel in each processed image.
- the relevant pixels may include a single pixel, pixels within a designated region of interest, or all pixels within each image.
- the intensity values may be used to create a set of measured ratios, R j , for each relevant pixel based on the different images, with each image corresponding to a different X-ray energy level.
- the set of measured ratios, R j for each pixel may be formed as shown in Table 1.
- at least one adjusted measured ratio, R m,j for each relevant pixel can be determined.
- the measured ratio can be based on the product of at least one estimation coefficient, ⁇ j , and the set of measured ratios, R j , from step 208 . Separate adjusted measured ratios, R m,j , may be determined for each measured ratio, R j , in the set or only for certain energy ratios within the set.
- the resulting set of adjusted measured ratios, R m,j may be compared to a material attenuation database in step 212 using at least one algorithm.
- the material attenuation database can be a densely populated attenuation ratio lookup table based on data from public sources of attenuation data. For example, this comparison may be a determination of an error, ⁇ , for each measured ratio with respect to the corresponding accepted attenuation ratio, as given by:
- ⁇ j is the estimation coefficient corresponding to the measured ratio
- R j and R T,j is an accepted attenuation ratio derived from the densely populated attenuation ratio lookup table.
- Results of the comparison in step 212 can be used in step 214 to determine the Z eff for each relevant pixel.
- step 212 can provide a comparison using a set of algorithms. Each algorithm may employ a different methodology, different adjusted measured ratios, and/or different estimation coefficients for comparison to the same database. The results of this comparison may be used to generate a set of candidate Z eff values and associated confidence values for each relevant pixel in step 214 .
- An arbitration step may be included in step 214 to select the Z eff value that has the highest associated confidence value for the given measurement conditions.
- the Z eff value for each relevant pixel may be used in the construction of a color image.
- colors may be assigned to correspond with different Z eff values.
- the assignment of colors may be predetermined or set dynamically to correspond with the range of measured Z eff values.
- colors may be assigned to each pixel according to the Z eff value of the pixel so as to generate a Z eff color image of the object. The method can end at step 220 .
- FIG. 3 illustrates a dynamic material attenuation processor 300 , which may be used as DMA module 112 in the embodiment of FIG. 1 .
- the dynamic material attenuation processor 300 may be used to dynamically create a densely populated attenuation ratio lookup table for use by the material assignment module 118 in the determination of Z eff values.
- the processor 300 can use configurable settings to create, for example, a table of accepted attenuation coefficient ratios, R T,j , with variable resolution in both the energy scale and Z eff . Similar to equation (4) above, the accepted attenuation coefficient ratio, R T,j , can be expressed as:
- Inputs into the dynamic material attenuation processor 300 may be derived from public attenuation data source 302 , such as periodic attenuations, liquid attenuations, or other forms available through the NIST public data source, for example.
- the public attenuation data source 302 may include a sparsely populated database.
- the data from the public attenuation data source 302 may provide attenuation data of various material compositions at different radiographic energy levels under a variety of scattering conditions, including, but not limited to, coherent, incoherent, photoelectric, pair-wise production nuclear field, pair-wise production electric field, total coherent, and total non-coherent.
- the processor 300 may build a fully connected table of Z eff values versus ratio of accepted attenuation coefficients at different energy levels.
- the table may have rows with different Z eff values.
- Each column in the table may have an accepted attenuation coefficient ratio, R T,j , that corresponds with a different pair of source energy levels used in the radiographic images by Z eff processor 116 .
- the densely populated attenuation ratio lookup table may be dynamically configurable based on user inputs or configuration data.
- the dynamic material attenuation processor 300 may receive a subsystem request 304 .
- This request 304 may include, for example, a particular source energy range and a particular Z eff range. This information may also be provided by configuration data.
- Configuration data can control the energy density needed for the resultant attenuation data table as well as the Z eff density needed.
- the requested energy values and energy resolution may be obtained from the configuration data.
- the requested Z eff values and Z eff resolution may be obtained from the configuration data.
- the configuration data may be integrated and stored with the processor 300 .
- the configuration data may be predetermined based on the radiographic system or object under test. Alternately, the configuration data may be input by a user.
- the dynamic material attenuation processor 300 can receive the request 304 , the energy configuration data 306 , and/or the Z eff configuration data 312 and may obtain the necessary attenuation data from public attenuation data source 302 .
- the attenuation data may be normalized for both energy and Z eff values.
- the resultant attenuation data may then be formed into accepted attenuation coefficient ratios, R T,j , and organized into tabular form at 324 .
- Publicly available attenuation data source 302 may lack some of the Z eff values or the source energy values requested in subsystem request 304 . Accordingly, the dynamic material attenuation processor 300 may create a high-resolution table so as to account for this missing information.
- the tabular form of accepted attenuation coefficient ratios, R T,j , from 324 may be expanded based on the subsystem request 304 , the energy configuration data 306 , and/or the Z eff configuration data 312 .
- Any data missing in the table generated by 326 may be interpolated using a polynomial spline fit.
- Poly-spline fit 310 may interface with the high-resolution table created in 326 to interpolate missing information with respect to source energy values.
- Poly-spline fit 316 may interface with the high-resolution table created in 326 to interpolated missing information with respect to Z eff values (i.e., material composition).
- the result can be an attenuation ratio lookup table based on source energy values and Z eff values.
- sparse input data and configuration information can be used to create a densely populated configurable attenuation ratio lookup table.
- a color map to the Z eff values may also be created to coincide with the Z eff value settings.
- the attenuation ratio lookup table may be output to 318 for assigning colors to the various Z eff values.
- the resultant attenuation ratio lookup table and the assigned colors may then be output to 320 for formatting to a given data standard.
- the attenuation ratio lookup table may be sent directly from 326 to 320 for formatting to the standard.
- the format may be any available data standard.
- the standard can be an N42 standard, such as ANSI N42.42.
- the attenuation ratio table may be output as a material attenuation ratio table 322 for use by material assignment module 118 .
- FIG. 4 illustrates a Z algorithm arbitration process 400 for use in estimating Z eff .
- the process starts at step 401 and continues to step 402 .
- a pixel may be selected in a region of interest for which at least one adjusted measured ratio, R m , has been determined.
- an algorithm may be selected from a set of algorithms. Using the selected algorithm, the at least one adjusted measured ratio, R m , may be compared to a densely-populated attenuation lookup table in step 406 . Based on the comparison, the selected algorithm can determine a Z eff value for the relevant pixel at step 408 .
- the selected algorithm may also determine a confidence value for the Z eff value.
- This confidence value may be based on measurement conditions. Alternately, the confidence value may be based on a particular set of estimation coefficients employed and the resulting range of Z eff determined.
- the arbitration process can check to see if all algorithms within the set have been evaluated. The process can repeat steps 404 - 410 until all algorithms have been evaluated. Once all algorithms have been evaluated, thereby generating a set of candidate Z eff values with associated confidence values, the process may proceed to step 414 .
- the Z eff value with the highest confidence value can be selected and assigned to the relevant pixel.
- the arbitration process can check to see if all relevant pixels have been evaluated. The process can repeat steps 402 - 414 until all relevant pixels have been evaluated.
- the process may proceed to step 418 , wherein the Z eff values and associated confidence values for each relevant pixel may be output for further processing.
- the Z eff values may be output to color image assignment module 120 for subsequent color imaging of the Z eff values.
- the process may continue to step 420 where the process may end.
- steps of the present invention may be repeated in whole or in part in order to perform the contemplated Z eff estimation. Further, it should be appreciated that the steps mentioned above may be performed on a single or distributed processor. Also, the processes, modules, and units described in the various figures of the embodiments above may be distributed across multiple computers or systems or may be co-located in a single processor or system.
- Embodiments of the method, system, and computer program product for determining an effective atomic number, Z eff , of an object may be implemented on a general-purpose computer, a special-purpose computer, a programmed microprocessor or microcontroller and peripheral integrated circuit element, an ASIC or other integrated circuit, a digital signal processor, a hardwired electronic or logic circuit such as a discrete element circuit, a programmed logic circuit such as a PLD, PLA, FPGA, PAL, or the like.
- any process capable of implementing the functions or steps described herein can be used to implement embodiments of the method, system, or computer program product for determining an effective atomic number, Z eff , of an object.
- embodiments of the disclosed method, system, and computer program product for determining an effective atomic number, Z eff , of an object may be readily implemented, fully or partially, in software using, for example, object or object-oriented software development environments that provide portable source code that can be used on a variety of computer platforms.
- embodiments of the disclosed method, system, and computer program product for determining an effective atomic number, Z eff , of an object can be implemented partially or fully in hardware using, for example, standard logic circuits or a VLSI design.
- Other hardware or software can be used to implement embodiments depending on the speed and/or efficiency requirements of the systems, the particular function, and/or particular software or hardware system, microprocessor, or microcomputer being utilized.
- Embodiments of the method, system, and computer program product for determining an effective atomic number, Z eff , of an object can be implemented in hardware and/or software using any known or later developed systems or structures, devices and/or software by those of ordinary skill in the applicable art from the function description provided herein and with a general basic knowledge of the computer, radiographic, and image processing arts.
- embodiments of the disclosed method, system, and computer program product for determining an effective atomic number, Z eff , of an object can be implemented in software executed on a programmed general purpose computer, a special purpose computer, a microprocessor, or the like.
- the Z eff determination method of this invention can be implemented as a program embedded on a personal computer such as a JAVA® or CGI script, as a resource residing on a server or image processing workstation, as a routine embedded in a dedicated processing system, or the like.
- the method and system can also be implemented by physically incorporating the method for determining Z eff of an object into a software and/or hardware system, such as the hardware and software systems of multi-energy X-ray inspection systems.
Abstract
A method of determining an effective atomic number, Zeff, of an object may include obtaining a plurality of radiographic images of the object. Each radiographic image can be obtained using a different independent X-ray energy level. An intensity value for each pixel in a region of interest in each radiographic image can be determined. A plurality of measured ratios, R, can be formed using attenuation coefficients from a pair of different radiographic images. At least one adjusted measured ratio, Rm, can be calculated based on the plurality of measured ratios, R, and at least one corresponding estimation coefficient, α. Zeff values can be assigned based on a comparison of the at least one adjusted measured ratio, Rm to a material attenuation database.
Description
- The present application claims the benefit of U.S. Provisional Patent Application No. 60/940,632, entitled “Threat Detection System”, filed May 29, 2007, which is hereby incorporated by reference in its entirety.
- The present invention relates generally to radiographic imaging, and, more particularly, to determination of the effective atomic number, Zeff, of a material using multi-energy radiographic imaging.
- Standard techniques using high-energy radiographic systems exist to determine or estimate atomic composition of a material. For example, two X-ray energy levels may be used to image an object of interest. The gray level intensity values measured for each of the two energy levels are used to compute a corresponding ratio. The ratio of the intensity values for an unknown material is compared against known materials. The known material with the closest ratio to that measured is used to estimate the unknown material's effective atomic number (Zeff) of the unknown material. However, these radiographic systems may be prone to noise and other non-linear effects that can cause errors in Zeff determination, especially in high-Z materials. Elements with high-Z include special nuclear materials, such as plutonium and highly enriched uranium, as well as elements that would be extremely effective in shielding special nuclear materials from passive radiation detection techniques.
- Using a radiographic system with only two X-ray energy levels, two common problems can occur. First, insufficient penetration may occur when high-Z or high density materials do not allow enough energy to penetrate through the object of interest to a detector. Second, low-Z or low density materials may induce over-saturation. In the case of over-saturation, little to no attenuation may occur at a particular X-ray energy level. To overcome the problem of over-saturation, lower energies may be used, whereas higher energies may be used in cases of insufficient penetration. Higher energies may solve insufficient penetration issues but exacerbate over-saturation issues and vice versa. Thus, in dual energy systems, the solution for overcoming over-saturation issues may be at odds with the solution for overcoming insufficient penetration. Embodiments of the present invention may address the above-mentioned problems and limitations, among other things.
- An embodiment of the present invention can include a method of determining an effective atomic number, Zeff, of an object including obtaining a plurality of radiographic images of the object. Each radiographic image can be obtained using a different independent X-ray energy level. Each image can also include a plurality of pixels. The images can be registered and normalized. Noise may be removed from the images. An intensity value for each pixel in the region of interest in each radiographic image can be determined. A plurality of measured ratios, R, for each pixel in the region of interest can be formed. Each measured ratio can be formed using intensity values from corresponding pixels in a pair of different radiographic images. The method can include calculating at least one adjusted measured ratio, Rm, for each pixel in the region of interest based on the plurality of measured ratios, R, and at least one corresponding estimation coefficient, α, and comparing the at least one adjusted measured ratio, Rm, for each pixel in the region of interest to a material attenuation database. The method can include assigning a Zeff value to each pixel in the region of interest based on the comparison and outputting the assigned Zeff values.
- Another embodiment may include a computer program product including a computer readable medium encoded with software instructions for causing a computer to perform the steps of receiving more than two radiographic images of an object and determining an intensity value for each pixel in each of the plurality of radiographic images. Each image can be obtained using a different independent X-ray energy level and can have a plurality of pixels. The steps can also include forming a plurality of measured ratios, R. Each measured ratio, R, can be formed using intensity values from corresponding pixels in a pair of different radiographic images. The steps can also include calculating at least one adjusted measured ratio, Rm, based on the plurality of measured ratios, R, and at least one corresponding estimation coefficient, α, and comparing the at least one adjusted measured ratio, Rm, to a material attenuation database. The steps can include assigning a Zeff value based on the comparison and outputting the assigned Zeff value.
- Another embodiment may include a system for determining an effective atomic number, Zeff, of an object. The system may include at least one processor. The at least one processor can form a plurality of measured ratios, R. Each measured ratio may be formed using intensity values taken from a pair of different radiographic images. The at least one processor can calculate at least one adjusted measured ratio, Rm, based on the plurality of measured ratios, R, and at least one corresponding estimation coefficient, α. The at least one processor can output a Zeff value based on a comparison of the at least one adjusted measured ratio, Rm, to a material attenuation database.
-
FIG. 1 is a block diagram showing an overview of an embodiment of a system for estimating Zeff. -
FIG. 2 is a flowchart showing a process overview of an embodiment of a process for estimating Zeff. -
FIG. 3 is a block diagram showing an embodiment of a dynamic material attenuation processor for use in estimating Zeff. -
FIG. 4 is a flowchart showing a process overview of an embodiment of Z assignment arbitration for use in estimating Zeff. - The present invention is directed to identifying material compositions by estimating their effective atomic number or Zeff value. The present invention further seeks to overcome the limitations of the above-mentioned ratio method by using a multi-energy approach. Multiple energies are used to determine intensity ratios that are proportional to attenuation coefficient ratios. These intensity ratios can then be used to build a series of measured ratios. These measured ratios can then be used to build a series of adjusted measured ratios for use in a Zeff determination that may be more accurate than conventional methods.
- Imaging of a material relies on the application of Beer's law. Beer's law equation can be stated using an arbitrary energy level k as,
-
I(k)=I o(k)e −μk t, (1) - where I(k) is the measured intensity of the radiation at energy level k, IO(k) is the input intensity of the radiation at energy level k, μk is the linear attenuation coefficient of the material at energy level k, and t is the energy independent material thickness.
- Experimentally, the linear attenuation coefficient depends on the material cross-section coefficient, σ. Since attenuation is energy dependent, each material has the basic regions of energy scatter, including photoelectric absorption σpe, coherent scatter σcs, incoherent scatter σis, and pair-wise production σpp. Each material has a material combined cross-sectional coefficient which depends on the sum of the individual cross-sections. For example,
-
σ=σpe+σcs+σis+σpp. (2) - The cross-sections and linear attenuation coefficients are related by:
-
- where μ is the material linear attenuation coefficient, σ is the material combined cross-sectional coefficient, N is Avogadro's number, A is atomic weight, and ρ is density.
- Using any two energies, denoted as x and y, the following intensity ratio may apply:
-
- where Rj is the ratio of the logarithm of the respective intensity ratios and x and y represent the upper and lower energy levels. This expression eliminates the dependency on mass thickness. Thus, the Rj value may be compared against attenuation curves at the upper and lower energy levels for various materials to determine the Zeff of an unknown material.
- Intensity values can be determined from measured radiographic images and used in the determination of a plurality of ratios for different pairs of radiographic source energy levels. Any number of energy levels in the X-ray regime may be used. For example, at least two energy levels and corresponding measured ratios can be used. In another example, four different energy levels and corresponding measured ratios can be used.
-
TABLE 1 Multi-Energy Attenuation Coefficient and Ratio Table Source EnergyLevels(1 < 2 < 3 < 4) CorrespondingAttenuationCoefficients(μx, μy) EnergyRegime E1, E2 μ1, μ2 Low Ri E1, E3 μ1, μ3 Low- Rii Medium E1, E4 μ1, μ4 Low-High Riii E2, E3 μ2, μ3 Medium Riv E2, E4 μ2, μ4 Medium- Rv High E3, E4 μ3, μ4 High Rvi - Table 1 shows an example of an application employing four different source energy levels. Four different source energy levels (E1-E4) are used to generate four different radiographic images of an object. E1 corresponds to the lowest X-ray source energy level. For example, E1 may be 1MeV. E4 corresponds to the highest X-ray source energy level. For example, E4 may be 10 MeV. Measured intensity values I(x) for each energy level, x, are extracted from the corresponding radiographic images and used in the determination of each measured ratio, Rj. Source intensity values Io(x) for each energy level are known from measurement conditions of the radiographic images and are also used in the determination of the measured ratios, Rj. Although only four source energy levels are discussed in detail with respect to Table 1, it is of course contemplated that more than four or less than four energy levels may be employed.
- Ratios exist for various combinations of intensity values based on the different source energy levels. Therefore, it may be feasible that one could select between low, medium, or high energy ratios for estimating the effective atomic number of the material of the object to minimize errors due to saturation, noise, or insufficient penetration.
- A family of attenuation curves for a group of materials with respect to X-ray energy is typically non-linear, with various ranges of polynomial order. Although the attenuation curves themselves are non-linear, the ratios from energy level to energy level can be modeled using linear methods. At least two energy solutions form a linear attenuation fit based on the ratio. Determining the atomic structure of a material can be based on comparing the slopes of the intensity values from radiographic images with those values of known materials residing in a database. The attenuation coefficient is proportional to measured X-ray intensities.
- Radiographic imagers can use measured X-ray intensity to create a corresponding gray level picture, otherwise known as a radiographic image. The conversion from intensity to gray level can be defined by an inverse linear ratio, but the process of imaging various materials poses several non-linear settings. Despite the nonlinearity, ratios between the measured intensity values are approximately the same as ratios between the true material linear attenuation coefficients when adjusted for gray level conversion. In other words, an adjusted ratio of measured intensity values, as defined in equation (4), may be correlated with a similar ratio of known linear attenuation coefficients in determining the Zeff of an unknown material.
- An adjusted measured ratio, Rm,j, can be determined from the equation:
-
Rm,j=αjRj - where Rj is the measured ratio of attenuation coefficients at two different source energy levels, αj is an estimation coefficient corresponding to the particular two energy levels used in determining Rj, and j is a subscript referring to the particular pair of source energy levels. The estimation coefficient, αj, can be determined from the mapping of the photon intensity values to a gray level value. Separate adjusted measured ratios, Rm, may be determined for each measured ratio, R, used. Alternately, the multiple separate measured ratios may be combined into a single vector. For example, using the set of ratios defined by Table 1, equation (5) would become:
-
R m=αi R iî i+αii R iiî ii=αiii R iiiî iii+αiv R ivî iv+αv R vî v+αiv R ivî iv - wherein îj represents a vector in n-space corresponding to attenuation coefficient Rj. The n-space may be a multi-dimensional space based on the number of energy levels employed or the number of measured ratios formed.
- While it is contemplated that all radiographic images corresponding to different energy levels can be used in the determination of the Zeff, it is also possible that only specific measured ratios of intensity values may be selected depending on the conditions associated with a particular energy level. For example, for a given energy level, an intensity value may not be applicable due to an over-saturation or non-penetration condition. Therefore, the measured energy ratios, R, using the intensity values for that particular energy level may be excluded from the calculation of the adjusted measured ratio, Rm,j. Alternately, the ratio Rm,j corresponding to those conditions of over-saturation or non-penetration may be excluded from use in the determination of Zeff.
- Gray level mapping of the estimation coefficient is nonlinear across material and energy levels. To compensate for the nonlinearity, the estimation coefficients, α, may be determined empirically through experimental evaluation of known materials. Alternately, an algorithm can be used for minimum error tuning so as to optimize the estimation coefficients, α. It should be appreciated that other methods may be used to determine the estimation coefficients.
- In view of the foregoing features described above, structures and methodologies in accordance with various aspects of the present invention will be better appreciated with reference to
FIGS. 1-5 . -
FIG. 1 is a block diagram showing an overview of an exemplary embodiment of a system for estimating Zeff. A materialdomain imaging processor 100 may receive multiple radiographic images and other data as input and may output a multi-energy high Z-mapping for identification of threats. Animage processor 108 may receive a set of radiographic images of an object of interest fromradiographic imaging system 102 or another processor, system, or database. Each of the radiographic images may be taken at a different X-ray energy level. Each image may be processed by theimage processor 108 to register the images, normalize the images, and remove any noise in the images. - The processed images may be sent to the Zeff processor 116. The Zeff processor 116 may extract corresponding intensity values from each pixel in each image. In an alternate embodiment, the Zeff processor 116 may only determine an intensity value for each pixel in each image that is within a region of interest. For determination of a region of interest, the images may be sent to an external analysis system, such as
segmentation processor 114.Segmentation processor 114 may use, for example, edge boundary detection, texture based detection, or other region processing in order to determine boundaries within the images and to select a common region of interest for evaluation by the materialdomain imaging processor 100. The region of interest in each image may then be communicated to the Zeff processor 116 for determination of intensity values. The relevant pixels for subsequent processing may be limited to those pixels within the region of interest instead of the entirety of pixels in each image. Relevant pixels may also be selected based on penetration conditions (e.g., oversaturation, non-penetration). For example, like regions across the various registered images may be compared to determine penetration conditions. - After determination of the intensity values, the Zeff processor 116 may use the intensities to form a set of measured ratios, Rj, for the relevant pixels. The measured ratios, Rj, may be formed using the relation set forth in equation (4) above. Normalized intensity values from corresponding pixels in a pair of different radiographic images may also be used to form each measured ratio. Corresponding pixels may refer to pixels in different images which correspond to the same point or location on an imaged object. For example, the set of ratios for each pixel may take the form shown in Table 1.
- The Zeff processor 116 may use estimation coefficients, αj, with the set of measured ratios, Rj, to determine at least one adjusted measured ratio, Rm,j, for each relevant pixel. Separate adjusted measured ratios, Rm,j, may be determined for each energy ratio in the set or only for certain energy ratios within the set. Alternately, the multiple energy ratios may also be combined into a single measured ratio vector, Rm. The estimation coefficients may be provided by tuning
module 110. - The relevant pixels may include, for example, a single pixel, all the pixels in the entire image, or just the pixels within a region of interest. In an alternative embodiment, relevant pixels may be selected based on penetration conditions (e.g., over-saturation, non-penetration, etc.). The material
domain imaging processor 100 or a separate system may compare regions across various images to determine penetration conditions. The materialdomain imaging processor 100 or a separate system (i.e.,segmentation processor 114 or context analysis system 106) may correlate regions with pixel values at a lower extreme of a photon intensity scale as regions of non-penetration, while regions with pixel values at an upper extreme of the photon intensity scale may be correlated as regions of over-saturation. The Zeff processor 116 may be configured to exclude these regions from the list of relevant pixels for further processing. - In an exemplary embodiment,
tuning module 110 can employ an algorithm for minimum error tuning of the estimation coefficients.Accepted attenuation database 104 may be provided to thetuning module 110. Alternately, accepted attenuation data may be provided by an integrated database, a memory device, or a separate system or processor. Thetuning module 110 may use ratios from radiographic images of a known material to optimize the value of the estimation coefficients such that adjusted measured ratios correspond to ratios of accepted attenuation data. The accepted attenuation data may be data from public attenuation sources, such as the NIST public data source. In an alternate embodiment, thetuning module 110 may include a database of previously determined estimation coefficients for use by the Zeff processor 116. The estimation coefficients may be determined via experimental evaluation or other means. - The resulting adjusted measured ratios, Rm,j, can be output from the Zeff processor 116 to the
material assignment module 118. Thematerial assignment module 118 can compare the adjusted measured ratios to a material attenuation database from dynamic material attenuation (DMA)module 112. Alternately, the Zeff processor 116 may send the measured ratios, Rj, to thematerial assignment module 118. In such a scenario,material assignment module 118 would use the estimation coefficients and the measured ratios to form the adjusted measured ratios, Rm,j, for comparison to the material attenuation database. In either embodiment, thematerial assignment module 118 may assign a Zeff value to each relevant pixel or to a region of pixels based on the comparison. - In an exemplary embodiment, the
DMA module 112 may be a processor, such as dynamicmaterial attenuation processor 300, that uses configurable settings to dynamically create a densely populated attenuation ratio lookup table with variable resolution in both the energy scale and the effective atomic number (Zeff) scale. TheDMA module 112 can be used for the two or more energies used in the radiographic images for determining Zeff values. Sparse attenuation data from public sources for any material (liquid, solid, gas) can be stored on a disk using a standard format. Alternately, accepted attenuation data (i.e., attenuation coefficients) may be input to theDMA module 112 from acceptedattenuation database 104. Accepted attenuation data may also be provided by an integrated database, a memory device, or a separate system or processor. The accepted attenuation data may be data from public attenuation sources, such as the NIST public data source. The attenuation data can be created using a variety of scattering approaches, such as coherent, incoherent, photoelectric, pairwise production nuclear field, pairwise production electric field, total scatter with coherent, and total scatter without coherent. A user may select the materials used to create the attenuation lookup table, such as water, peroxide, lead, carbon, etc. Using the attenuation data, the attenuation lookup table may be created as a function of the X-ray energy level and the Zeff of the material. Alternately, theDMA module 112 may include a database of previously determined attenuation lookup tables, such as material attenuation table 322 inFIG. 3 , for use by thematerial assignment module 118. -
Material assignment module 118 may include a plurality of independent material assignment algorithms. Each material assignment algorithm may employ a different methodology or use different measured ratios for comparison to the same attenuation database to generate a set of candidate Zeff values for each relevant pixel. Each material assignment algorithm can also assign a confidence value to the candidate Zeff values. - The set of candidate Zeff values may be sent to the Zeff processor 116. The Zeff processor 116 may compare the confidence values for each Zeff value in the set and may select the Zeff value for each relevant pixel with the highest associated confidence value. The result can be a Zeff image with each pixel having an assigned Zeff value with the highest confidence. In an exemplary embodiment, the Zeff processor 116 may output the result to external analysis system, such as
context analysis system 106, for evaluation of regions of non-penetration or over-saturation for inclusion in the final processor output. For example,context analysis system 106 may use a-priori information, in the way of configuration data, to assist in identifying non-penetrable and false alarm cases. - The result from the Zeff processor 116 can be sent to a color
image assignment module 120. The colorimage assignment module 120 can map a color scale to a range of corresponding Zeff values. Thus, each relevant pixel may be assigned a color based on the assigned Zeff value, thereby creating a color image of the object. This color image can be output from the materialdomain imaging processor 100 tothreat decision processor 124. In addition, the color image and Zeff values may be further processed by athreat assignment module 122. The Zeff value for each pixel can be compared with a threat threshold by thethreat assignment module 122 to determine regions where the threshold is exceeded, i.e., those regions where a threat exists. These regions and associated confidence values may be output from the materialdomain imaging processor 100 to thethreat decision processor 124 for further processing or integrated decision making. -
FIG. 2 shows aprocess flow 200 of an exemplary embodiment of a method of material domain image processing. The process begins atstep 201 and proceeds to step 202. Instep 202, radiographic images may be obtained of an object. Each radiographic image of an object may be obtained at a different X-ray energy level. For example, at least two independent X-ray energy levels can be used. In another example, four X-ray energy levels can be used to generate four independent radiographic images of an object. Instep 204, the radiographic images may be subject to image processing. The image processing may be any of a number of processing steps known in the art. For example, each radiographic image may be registered with the other radiographic images such that regions of interest in the object are aligned. Further, the radiographic images may be normalized to a gray scale. In addition, noise filtering may be employed to reduce noise artifacts in the image. - In
step 206, an intensity value may be determined for each relevant pixel in each processed image. The relevant pixels may include a single pixel, pixels within a designated region of interest, or all pixels within each image. Atstep 208, the intensity values may be used to create a set of measured ratios, Rj, for each relevant pixel based on the different images, with each image corresponding to a different X-ray energy level. For example, the set of measured ratios, Rj, for each pixel may be formed as shown in Table 1. Atstep 210, at least one adjusted measured ratio, Rm,j, for each relevant pixel can be determined. The measured ratio can be based on the product of at least one estimation coefficient, αj, and the set of measured ratios, Rj, fromstep 208. Separate adjusted measured ratios, Rm,j, may be determined for each measured ratio, Rj, in the set or only for certain energy ratios within the set. - The resulting set of adjusted measured ratios, Rm,j, may be compared to a material attenuation database in
step 212 using at least one algorithm. The material attenuation database can be a densely populated attenuation ratio lookup table based on data from public sources of attenuation data. For example, this comparison may be a determination of an error, ε, for each measured ratio with respect to the corresponding accepted attenuation ratio, as given by: -
ε=|R m,j −R T,j|=|(αj *R j)−R T,j| (7) - where αj is the estimation coefficient corresponding to the measured ratio, Rj, and RT,j is an accepted attenuation ratio derived from the densely populated attenuation ratio lookup table.
- Results of the comparison in
step 212 can be used instep 214 to determine the Zeff for each relevant pixel. In an exemplary embodiment, step 212 can provide a comparison using a set of algorithms. Each algorithm may employ a different methodology, different adjusted measured ratios, and/or different estimation coefficients for comparison to the same database. The results of this comparison may be used to generate a set of candidate Zeff values and associated confidence values for each relevant pixel instep 214. An arbitration step may be included instep 214 to select the Zeff value that has the highest associated confidence value for the given measurement conditions. - In
steps step 216, colors may be assigned to correspond with different Zeff values. For example, the assignment of colors may be predetermined or set dynamically to correspond with the range of measured Zeff values. Instep 218, colors may be assigned to each pixel according to the Zeff value of the pixel so as to generate a Zeff color image of the object. The method can end atstep 220. -
FIG. 3 illustrates a dynamicmaterial attenuation processor 300, which may be used asDMA module 112 in the embodiment ofFIG. 1 . The dynamicmaterial attenuation processor 300 may be used to dynamically create a densely populated attenuation ratio lookup table for use by thematerial assignment module 118 in the determination of Zeff values. Theprocessor 300 can use configurable settings to create, for example, a table of accepted attenuation coefficient ratios, RT,j, with variable resolution in both the energy scale and Zeff. Similar to equation (4) above, the accepted attenuation coefficient ratio, RT,j, can be expressed as: -
- Inputs into the dynamic
material attenuation processor 300 may be derived from publicattenuation data source 302, such as periodic attenuations, liquid attenuations, or other forms available through the NIST public data source, for example. The publicattenuation data source 302 may include a sparsely populated database. The data from the publicattenuation data source 302 may provide attenuation data of various material compositions at different radiographic energy levels under a variety of scattering conditions, including, but not limited to, coherent, incoherent, photoelectric, pair-wise production nuclear field, pair-wise production electric field, total coherent, and total non-coherent. Using this data, theprocessor 300 may build a fully connected table of Zeff values versus ratio of accepted attenuation coefficients at different energy levels. For example, the table may have rows with different Zeff values. Each column in the table may have an accepted attenuation coefficient ratio, RT,j, that corresponds with a different pair of source energy levels used in the radiographic images by Zeff processor 116. - The densely populated attenuation ratio lookup table may be dynamically configurable based on user inputs or configuration data. For example, the dynamic
material attenuation processor 300 may receive asubsystem request 304. Thisrequest 304 may include, for example, a particular source energy range and a particular Zeff range. This information may also be provided by configuration data. Configuration data can control the energy density needed for the resultant attenuation data table as well as the Zeff density needed. At 306, the requested energy values and energy resolution may be obtained from the configuration data. At 312, the requested Zeff values and Zeff resolution may be obtained from the configuration data. Although shown inFIG. 3 as being external toprocessor 300, the configuration data may be integrated and stored with theprocessor 300. The configuration data may be predetermined based on the radiographic system or object under test. Alternately, the configuration data may be input by a user. - At 308, the dynamic
material attenuation processor 300 can receive therequest 304, theenergy configuration data 306, and/or the Zeff configuration data 312 and may obtain the necessary attenuation data from publicattenuation data source 302. At 314, the attenuation data may be normalized for both energy and Zeff values. The resultant attenuation data may then be formed into accepted attenuation coefficient ratios, RT,j, and organized into tabular form at 324. Publicly availableattenuation data source 302 may lack some of the Zeff values or the source energy values requested insubsystem request 304. Accordingly, the dynamicmaterial attenuation processor 300 may create a high-resolution table so as to account for this missing information. At 326, the tabular form of accepted attenuation coefficient ratios, RT,j, from 324 may be expanded based on thesubsystem request 304, theenergy configuration data 306, and/or the Zeff configuration data 312. Any data missing in the table generated by 326 may be interpolated using a polynomial spline fit. Poly-spline fit 310 may interface with the high-resolution table created in 326 to interpolate missing information with respect to source energy values. Poly-spline fit 316 may interface with the high-resolution table created in 326 to interpolated missing information with respect to Zeff values (i.e., material composition). The result can be an attenuation ratio lookup table based on source energy values and Zeff values. Thus, sparse input data and configuration information can be used to create a densely populated configurable attenuation ratio lookup table. - A color map to the Zeff values may also be created to coincide with the Zeff value settings. The attenuation ratio lookup table may be output to 318 for assigning colors to the various Zeff values. The resultant attenuation ratio lookup table and the assigned colors may then be output to 320 for formatting to a given data standard. Alternately, the attenuation ratio lookup table may be sent directly from 326 to 320 for formatting to the standard. The format may be any available data standard. For example, the standard can be an N42 standard, such as ANSI N42.42. After formatting, the attenuation ratio table may be output as a material attenuation ratio table 322 for use by
material assignment module 118. -
FIG. 4 illustrates a Zalgorithm arbitration process 400 for use in estimating Zeff. The process starts atstep 401 and continues to step 402. Instep 402, a pixel may be selected in a region of interest for which at least one adjusted measured ratio, Rm, has been determined. Atstep 404, an algorithm may be selected from a set of algorithms. Using the selected algorithm, the at least one adjusted measured ratio, Rm, may be compared to a densely-populated attenuation lookup table instep 406. Based on the comparison, the selected algorithm can determine a Zeff value for the relevant pixel atstep 408. Atstep 410, the selected algorithm may also determine a confidence value for the Zeff value. This confidence value may be based on measurement conditions. Alternately, the confidence value may be based on a particular set of estimation coefficients employed and the resulting range of Zeff determined. Atstep 412, the arbitration process can check to see if all algorithms within the set have been evaluated. The process can repeat steps 404-410 until all algorithms have been evaluated. Once all algorithms have been evaluated, thereby generating a set of candidate Zeff values with associated confidence values, the process may proceed to step 414. Atstep 414, the Zeff value with the highest confidence value can be selected and assigned to the relevant pixel. Atstep 416, the arbitration process can check to see if all relevant pixels have been evaluated. The process can repeat steps 402-414 until all relevant pixels have been evaluated. If all pixels have been evaluated, the process may proceed to step 418, wherein the Zeff values and associated confidence values for each relevant pixel may be output for further processing. For example, the Zeff values may be output to colorimage assignment module 120 for subsequent color imaging of the Zeff values. The process may continue to step 420 where the process may end. - It should be appreciated that the steps of the present invention may be repeated in whole or in part in order to perform the contemplated Zeff estimation. Further, it should be appreciated that the steps mentioned above may be performed on a single or distributed processor. Also, the processes, modules, and units described in the various figures of the embodiments above may be distributed across multiple computers or systems or may be co-located in a single processor or system.
- Embodiments of the method, system, and computer program product for determining an effective atomic number, Zeff, of an object, may be implemented on a general-purpose computer, a special-purpose computer, a programmed microprocessor or microcontroller and peripheral integrated circuit element, an ASIC or other integrated circuit, a digital signal processor, a hardwired electronic or logic circuit such as a discrete element circuit, a programmed logic circuit such as a PLD, PLA, FPGA, PAL, or the like. In general, any process capable of implementing the functions or steps described herein can be used to implement embodiments of the method, system, or computer program product for determining an effective atomic number, Zeff, of an object.
- Furthermore, embodiments of the disclosed method, system, and computer program product for determining an effective atomic number, Zeff, of an object may be readily implemented, fully or partially, in software using, for example, object or object-oriented software development environments that provide portable source code that can be used on a variety of computer platforms. Alternatively, embodiments of the disclosed method, system, and computer program product for determining an effective atomic number, Zeff, of an object can be implemented partially or fully in hardware using, for example, standard logic circuits or a VLSI design. Other hardware or software can be used to implement embodiments depending on the speed and/or efficiency requirements of the systems, the particular function, and/or particular software or hardware system, microprocessor, or microcomputer being utilized. Embodiments of the method, system, and computer program product for determining an effective atomic number, Zeff, of an object can be implemented in hardware and/or software using any known or later developed systems or structures, devices and/or software by those of ordinary skill in the applicable art from the function description provided herein and with a general basic knowledge of the computer, radiographic, and image processing arts.
- Moreover, embodiments of the disclosed method, system, and computer program product for determining an effective atomic number, Zeff, of an object can be implemented in software executed on a programmed general purpose computer, a special purpose computer, a microprocessor, or the like. Also, the Zeff determination method of this invention can be implemented as a program embedded on a personal computer such as a JAVA® or CGI script, as a resource residing on a server or image processing workstation, as a routine embedded in a dedicated processing system, or the like. The method and system can also be implemented by physically incorporating the method for determining Zeff of an object into a software and/or hardware system, such as the hardware and software systems of multi-energy X-ray inspection systems.
- It is, therefore, apparent that there is provided, in accordance with the present invention, a method, system, and computer program product for determining Zeff of an object. While this invention has been described in conjunction with a number of embodiments, it is evident that many alternatives, modifications and variations would be or are apparent to those of ordinary skill in the applicable arts. Accordingly, Applicants intend to embrace all such alternatives, modifications, equivalents and variations that are within the spirit and scope of this invention.
Claims (20)
1. A method of determining an effective atomic number, Zeff, of an object, the method comprising the steps of:
obtaining a plurality of radiographic images of the object, each radiographic image obtained using a different independent X-ray energy level and each image including a plurality of pixels;
registering the plurality of radiographic images with each other;
normalizing the plurality of radiographic images;
removing noise from the plurality of radiographic images;
determining an intensity value for each pixel in a region of interest in each radiographic image;
forming a plurality of measured ratios for each pixel in the region of interest, each ratio formed using intensity values from corresponding pixels in a pair of different radiographic images;
calculating at least one adjusted measured ratio for each pixel in the region of interest based on the plurality of measured ratios and at least one corresponding estimation coefficient;
comparing the at least one adjusted measured ratio for each pixel in the region of interest to a material attenuation database;
assigning a Zeff value to each pixel in the region of interest based on the comparing; and
outputting the assigned Zeff values.
2. The method of claim 1 , wherein the step of obtaining a plurality of radiographic images of the object comprises obtaining more than two radiographic images of the object.
3. The method of claim 1 , wherein said at least one adjusted measured ratio includes a set of adjusted measured ratios and each adjusted measured ratio of the set is based on a corresponding one of the plurality of measured ratios and a corresponding one of a plurality of estimation coefficients.
4. The method of claim 3 , wherein comparing the at least one adjusted measured ratio to a database includes:
providing a plurality of independent material assignment systems, each material assignment system configured to determine a Zeff value based on at least one different adjusted measured ratio from the set of adjusted measured ratios and the material attenuation database;
generating a plurality of candidate Zeff values and associated confidence values using the plurality of material assignment systems; and
selecting the Zeff value from the plurality of candidate Zeff values with the highest associated confidence value.
5. The method of claim 1 , further comprising:
mapping a color scale to a range of Zeff values; and
generating a color image of the object using the color mapping and the assigned Zeff values.
6. The method of claim 1 , further comprising:
comparing the Zeff values to a threat threshold; and
generating an output of a threat region of interest and a confidence value for the threat region of interest based on the comparison of the Zeff values to the threat threshold.
7. The method of claim 1 , wherein the material assignment database comprises a densely populated attenuation ratio lookup table.
8. The method of claim 7 , further comprising:
receiving accepted attenuation data of various materials, the accepted attenuation data providing information regarding attenuation of X-rays at various energy levels based on the corresponding Zeff values of said materials;
receiving configuration data including energy scale resolution and Zeff resolution for said attenuation lookup table; and
generating said attenuation lookup table based on said accepted attenuation data and said configuration data.
9. The method of claim 1 , further comprising:
determining said region of interest through image processing and analysis of the radiographic images.
10. A computer program product comprising:
a computer readable medium encoded with software instructions that, when executed by a computer, cause the computer to perform the steps of:
receiving more than two radiographic images of an object, each image obtained using a different independent X-ray energy level and each image having a plurality of pixels;
determining a normalized intensity value for each pixel in each of the radiographic images of the object;
forming a plurality of measured ratios, each measured ratio formed using the normalized intensity values from corresponding pixels in a pair of different radiographic images;
calculating at least one adjusted measured ratio based on the plurality of measured ratios and at least one corresponding estimation coefficient;
comparing the at least one adjusted measured ratio to a material attenuation database;
assigning a Zeff value based on the comparing; and
outputting the assigned Zeff value.
11. The computer program product of claim 10 , wherein the determining a normalized intensity value includes:
registering the radiographic images with each other;
normalizing the radiographic images;
processing the radiographic images to remove noise; and
determining a plurality of intensity values, each of the plurality of intensity values being determined for a different pixel of the plurality of pixels in one of registered, normalized, and processed radiographic images.
12. The computer program product of claim 10 , wherein the comparing the at least one adjusted measured ratio, includes:
providing a plurality of independent material assignment algorithms for determining Zeff based on the at least one adjusted measured ratio and the material attenuation database;
generating a plurality of candidate Zeff values and associated confidence values using the plurality of independent material assignment algorithms; and
selecting the Zeff value from the plurality of candidate Zeff values with the highest associated confidence value.
13. The computer program product of claim 10 , further comprising the steps of:
mapping a color scale to a range of Zeff values;
generating a color image of the object using the color mapping and the assigned Zeff values;
comparing the Zeff values to a threat threshold to determine a threat region of interest; and
generating an output of the color image, a threat region of interest, and a confidence value for the threat region of interest based on the comparison of the Zeff values to the threat threshold.
14. The computer program product of claim 10 , wherein the material attenuation database comprises a densely populated attenuation lookup table based on accepted attenuation data, the accepted attenuation data providing information regarding attenuation of X-rays at various energy levels based on Zeff values.
15. A system for determining an effective atomic number, Zeff, of an object, the system comprising:
at least one processor configured to perform the steps of:
forming a plurality of measured ratios, each measured ratio formed using intensity values from a pair of different radiographic images of an object;
calculating at least one adjusted measured ratio based on the plurality of measured ratios and at least one corresponding estimation coefficient; and
outputting a Zeff value based on a comparison of the at least one adjusted measured ratio to a material attenuation database.
16. The system of claim 15 , further comprising:
an image processing module configured to register, normalize, and remove noise from each of the radiographic images.
17. The system of claim 15 , further comprising:
a plurality of independent material assignment algorithms for determining a plurality of Zeff values and corresponding confidence values, each material assignment algorithm configured to compare the at least one adjusted measured ratio to the material attenuation database to determine a Zeff value and an associated confidence value,
wherein the at least one processor is configured to assign a Zeff value from the plurality of determined Zeff values based on the associated confidence value and to generate an output at least based on the assigned Zeff value.
18. The system of claim 15 , further comprising:
a color image assignment module configured to map a color scale to a range of Zeff values and to generate a color image of the object using the color mapping and assigned Zeff values.
19. The system of claim 15 , further comprising:
a threat assignment module configured to compare the Zeff values to a threat threshold and to generate an output of at least one threat region of interest and a corresponding confidence value for each threat region of interest.
20. The system of claim 15 , further comprising:
a material attenuation module configured to dynamically generate a densely populated attenuation ratio lookup table based on accepted attenuation data and configuration data.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/129,371 US20090052762A1 (en) | 2007-05-29 | 2008-05-29 | Multi-energy radiographic system for estimating effective atomic number using multiple ratios |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US94063207P | 2007-05-29 | 2007-05-29 | |
US12/129,371 US20090052762A1 (en) | 2007-05-29 | 2008-05-29 | Multi-energy radiographic system for estimating effective atomic number using multiple ratios |
Publications (1)
Publication Number | Publication Date |
---|---|
US20090052762A1 true US20090052762A1 (en) | 2009-02-26 |
Family
ID=40088192
Family Applications (7)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/129,036 Abandoned US20090003651A1 (en) | 2007-05-29 | 2008-05-29 | Object segmentation recognition |
US12/129,055 Abandoned US20090052622A1 (en) | 2007-05-29 | 2008-05-29 | Nuclear material detection system |
US12/129,371 Abandoned US20090052762A1 (en) | 2007-05-29 | 2008-05-29 | Multi-energy radiographic system for estimating effective atomic number using multiple ratios |
US12/129,383 Expired - Fee Related US8094874B2 (en) | 2007-05-29 | 2008-05-29 | Material context analysis |
US12/129,393 Abandoned US20090055344A1 (en) | 2007-05-29 | 2008-05-29 | System and method for arbitrating outputs from a plurality of threat analysis systems |
US12/129,439 Abandoned US20080298544A1 (en) | 2007-05-29 | 2008-05-29 | Genetic tuning of coefficients in a threat detection system |
US12/129,410 Abandoned US20090003699A1 (en) | 2007-05-29 | 2008-05-29 | User guided object segmentation recognition |
Family Applications Before (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/129,036 Abandoned US20090003651A1 (en) | 2007-05-29 | 2008-05-29 | Object segmentation recognition |
US12/129,055 Abandoned US20090052622A1 (en) | 2007-05-29 | 2008-05-29 | Nuclear material detection system |
Family Applications After (4)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/129,383 Expired - Fee Related US8094874B2 (en) | 2007-05-29 | 2008-05-29 | Material context analysis |
US12/129,393 Abandoned US20090055344A1 (en) | 2007-05-29 | 2008-05-29 | System and method for arbitrating outputs from a plurality of threat analysis systems |
US12/129,439 Abandoned US20080298544A1 (en) | 2007-05-29 | 2008-05-29 | Genetic tuning of coefficients in a threat detection system |
US12/129,410 Abandoned US20090003699A1 (en) | 2007-05-29 | 2008-05-29 | User guided object segmentation recognition |
Country Status (1)
Country | Link |
---|---|
US (7) | US20090003651A1 (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120069963A1 (en) * | 2010-08-06 | 2012-03-22 | Telesecurity Sciences, Inc. | Dual energy backscatter x-ray shoe scanning device |
US20120170826A1 (en) * | 2009-09-22 | 2012-07-05 | Kwang Eun Jang | Multi-energy x-ray system, image processing apparatus for discriminating multi-energy x-ray material, and image processing method for material discrimination of multi-energy x-ray system |
US20140140601A1 (en) * | 2011-06-30 | 2014-05-22 | Andrew Litvin | Iterative image reconstruction |
US9476923B2 (en) | 2011-06-30 | 2016-10-25 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Method and device for identifying a material by the spectral analysis of electromagnetic radiation passing through said material |
US10078150B2 (en) | 2015-04-14 | 2018-09-18 | Board Of Regents, The University Of Texas System | Detecting and quantifying materials in containers utilizing an inverse algorithm with adaptive regularization |
US10216866B2 (en) * | 2011-02-25 | 2019-02-26 | Smiths Heimann Gmbh | Image reconstruction based on parametric models |
US10379252B2 (en) | 2015-07-22 | 2019-08-13 | Tsinghua University | Methods and devices for estimating weight of an object to be inspected in an inspection system |
Families Citing this family (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090003651A1 (en) * | 2007-05-29 | 2009-01-01 | Peter Dugan | Object segmentation recognition |
DE102007028895B4 (en) * | 2007-06-22 | 2010-07-15 | Siemens Ag | Method for segmenting structures in 3D image data sets |
US8200015B2 (en) * | 2007-06-22 | 2012-06-12 | Siemens Aktiengesellschaft | Method for interactively segmenting structures in image data records and image processing unit for carrying out the method |
KR20100038046A (en) * | 2008-10-02 | 2010-04-12 | 가부시키가이샤 한도오따이 에네루기 켄큐쇼 | Touch panel and method for driving the same |
JP5740132B2 (en) | 2009-10-26 | 2015-06-24 | 株式会社半導体エネルギー研究所 | Display device and semiconductor device |
US20120113146A1 (en) * | 2010-11-10 | 2012-05-10 | Patrick Michael Virtue | Methods, apparatus and articles of manufacture to combine segmentations of medical diagnostic images |
US8924325B1 (en) * | 2011-02-08 | 2014-12-30 | Lockheed Martin Corporation | Computerized target hostility determination and countermeasure |
WO2013052549A1 (en) * | 2011-10-03 | 2013-04-11 | Cornell University | System and methods of acoustic monitoring |
JP5895624B2 (en) * | 2012-03-14 | 2016-03-30 | オムロン株式会社 | Image processing apparatus, image processing method, control program, and recording medium |
US9589188B2 (en) * | 2012-11-14 | 2017-03-07 | Varian Medical Systems, Inc. | Method and apparatus pertaining to identifying objects of interest in a high-energy image |
GB2508841A (en) * | 2012-12-12 | 2014-06-18 | Ibm | Computing prioritised general arbitration rules for conflicting rules |
US9697467B2 (en) | 2014-05-21 | 2017-07-04 | International Business Machines Corporation | Goal-driven composition with preferences method and system |
US9785755B2 (en) | 2014-05-21 | 2017-10-10 | International Business Machines Corporation | Predictive hypothesis exploration using planning |
US9118714B1 (en) * | 2014-07-23 | 2015-08-25 | Lookingglass Cyber Solutions, Inc. | Apparatuses, methods and systems for a cyber threat visualization and editing user interface |
GB2530252B (en) * | 2014-09-10 | 2020-04-01 | Smiths Heimann Sas | Determination of a degree of homogeneity in images |
CN104482996B (en) * | 2014-12-24 | 2019-03-15 | 胡桂标 | The material kind of passive nuclear level sensing device corrects measuring system |
CN104778444B (en) * | 2015-03-10 | 2018-01-16 | 公安部交通管理科学研究所 | The appearance features analysis method of vehicle image under road scene |
US9687207B2 (en) * | 2015-04-01 | 2017-06-27 | Toshiba Medical Systems Corporation | Pre-reconstruction calibration, data correction, and material decomposition method and apparatus for photon-counting spectrally-resolving X-ray detectors and X-ray imaging |
US9760801B2 (en) | 2015-05-12 | 2017-09-12 | Lawrence Livermore National Security, Llc | Identification of uncommon objects in containers |
IL239191A0 (en) * | 2015-06-03 | 2015-11-30 | Amir B Geva | Image classification system |
US11841789B2 (en) | 2016-01-27 | 2023-12-12 | Microsoft Technology Licensing, Llc | Visual aids for debugging |
US11868896B2 (en) | 2016-01-27 | 2024-01-09 | Microsoft Technology Licensing, Llc | Interface for working with simulations on premises |
US11836650B2 (en) | 2016-01-27 | 2023-12-05 | Microsoft Technology Licensing, Llc | Artificial intelligence engine for mixing and enhancing features from one or more trained pre-existing machine-learning models |
US10671938B2 (en) | 2016-01-27 | 2020-06-02 | Bonsai AI, Inc. | Artificial intelligence engine configured to work with a pedagogical programming language to train one or more trained artificial intelligence models |
US20180357543A1 (en) * | 2016-01-27 | 2018-12-13 | Bonsai AI, Inc. | Artificial intelligence system configured to measure performance of artificial intelligence over time |
US11775850B2 (en) | 2016-01-27 | 2023-10-03 | Microsoft Technology Licensing, Llc | Artificial intelligence engine having various algorithms to build different concepts contained within a same AI model |
US10204226B2 (en) | 2016-12-07 | 2019-02-12 | General Electric Company | Feature and boundary tuning for threat detection in industrial asset control system |
US11120297B2 (en) * | 2018-11-30 | 2021-09-14 | International Business Machines Corporation | Segmentation of target areas in images |
US10939044B1 (en) * | 2019-08-27 | 2021-03-02 | Adobe Inc. | Automatically setting zoom level for image capture |
Citations (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US538758A (en) * | 1895-05-07 | Richard watkins | ||
US4794257A (en) * | 1984-05-14 | 1988-12-27 | Matsushita Electric Industrial Co., Ltd. | Quantum-counting radiography |
US5132998A (en) * | 1989-03-03 | 1992-07-21 | Matsushita Electric Industrial Co., Ltd. | Radiographic image processing method and photographic imaging apparatus therefor |
US5319547A (en) * | 1990-08-10 | 1994-06-07 | Vivid Technologies, Inc. | Device and method for inspection of baggage and other objects |
US5600700A (en) * | 1995-09-25 | 1997-02-04 | Vivid Technologies, Inc. | Detecting explosives or other contraband by employing transmitted and scattered X-rays |
US5600303A (en) * | 1993-01-15 | 1997-02-04 | Technology International Incorporated | Detection of concealed explosives and contraband |
US5642393A (en) * | 1995-09-26 | 1997-06-24 | Vivid Technologies, Inc. | Detecting contraband by employing interactive multiprobe tomography |
US6018562A (en) * | 1995-11-13 | 2000-01-25 | The United States Of America As Represented By The Secretary Of The Army | Apparatus and method for automatic recognition of concealed objects using multiple energy computed tomography |
US6026171A (en) * | 1998-02-11 | 2000-02-15 | Analogic Corporation | Apparatus and method for detection of liquids in computed tomography data |
US6236709B1 (en) * | 1998-05-04 | 2001-05-22 | Ensco, Inc. | Continuous high speed tomographic imaging system and method |
US6282305B1 (en) * | 1998-06-05 | 2001-08-28 | Arch Development Corporation | Method and system for the computerized assessment of breast cancer risk |
US20010033636A1 (en) * | 1999-11-13 | 2001-10-25 | Martin Hartick | Method and apparatus for determining a material of a detected item |
US6567496B1 (en) * | 1999-10-14 | 2003-05-20 | Sychev Boris S | Cargo inspection apparatus and process |
US20040247075A1 (en) * | 2003-06-06 | 2004-12-09 | Johnson James H. | Vehicle mounted inspection systems and methods |
US20050002550A1 (en) * | 2003-07-03 | 2005-01-06 | Ge Medical Systems Global Technology Company, Llc | Imaging chain for digital tomosynthesis on a flat panel detector |
US20050025280A1 (en) * | 2002-12-10 | 2005-02-03 | Robert Schulte | Volumetric 3D x-ray imaging system for baggage inspection including the detection of explosives |
US20050031075A1 (en) * | 2003-08-07 | 2005-02-10 | Hopkins Forrest Frank | System and method for detecting an object |
US20050058242A1 (en) * | 2003-09-15 | 2005-03-17 | Peschmann Kristian R. | Methods and systems for the rapid detection of concealed objects |
US20050111619A1 (en) * | 2002-02-06 | 2005-05-26 | L-3 Communications Security And Detection Systems Corporation Delaware | Method and apparatus for target transmitting information about a target object between a prescanner and a CT scanner |
US20050180542A1 (en) * | 2004-02-17 | 2005-08-18 | General Electric Company | CT-Guided system and method for analyzing regions of interest for contraband detection |
US20050256820A1 (en) * | 2004-05-14 | 2005-11-17 | Lockheed Martin Corporation | Cognitive arbitration system |
US6969861B2 (en) * | 2001-10-02 | 2005-11-29 | Konica Corporation | Cassette for radiographic imaging, radiographic image reading apparatus and radiographic image reading method |
US20060098773A1 (en) * | 2003-09-15 | 2006-05-11 | Peschmann Kristian R | Methods and systems for rapid detection of concealed objects using fluorescence |
US7092485B2 (en) * | 2003-05-27 | 2006-08-15 | Control Screening, Llc | X-ray inspection system for detecting explosives and other contraband |
US7103137B2 (en) * | 2002-07-24 | 2006-09-05 | Varian Medical Systems Technology, Inc. | Radiation scanning of objects for contraband |
US20060204107A1 (en) * | 2005-03-04 | 2006-09-14 | Lockheed Martin Corporation | Object recognition system using dynamic length genetic training |
US7123762B2 (en) * | 2002-02-08 | 2006-10-17 | University Of Chicago | Method and system for risk-modulated diagnosis of disease |
US20060233302A1 (en) * | 2004-10-22 | 2006-10-19 | Might Matthew B | Angled-beam detection system for container inspection |
US7130371B2 (en) * | 2002-09-27 | 2006-10-31 | Scantech Holdings, Llc | System for alternately pulsing energy of accelerated electrons bombarding a conversion target |
US20060256914A1 (en) * | 2004-11-12 | 2006-11-16 | Might Matthew B | Non-intrusive container inspection system using forward-scattered radiation |
US20060269114A1 (en) * | 2003-07-03 | 2006-11-30 | General Electric Company | Methods and systems for prescribing parameters for tomosynthesis |
US7162007B2 (en) * | 2004-02-06 | 2007-01-09 | Elyan Vladimir V | Non-intrusive inspection systems for large container screening and inspection |
US7162005B2 (en) * | 2002-07-19 | 2007-01-09 | Varian Medical Systems Technologies, Inc. | Radiation sources and compact radiation scanning systems |
US20070009084A1 (en) * | 2005-06-01 | 2007-01-11 | Endicott Interconnect Technologies, Inc. | Imaging inspection apparatus with directional cooling |
US7190757B2 (en) * | 2004-05-21 | 2007-03-13 | Analogic Corporation | Method of and system for computing effective atomic number images in multi-energy computed tomography |
US20070248212A1 (en) * | 2004-10-22 | 2007-10-25 | Might Matthew B | Cryptographic container security system |
US7399976B2 (en) * | 2005-11-03 | 2008-07-15 | Tsinghua University | Method and equipment for discriminating materials by employing fast neutron and continuous spectral x-ray |
US20080298544A1 (en) * | 2007-05-29 | 2008-12-04 | Peter Dugan | Genetic tuning of coefficients in a threat detection system |
US7491958B2 (en) * | 2003-08-27 | 2009-02-17 | Scantech Holdings, Llc | Radiographic inspection system for inspecting the contents of a container having dual injector and dual accelerating section |
US7817833B2 (en) * | 2004-05-26 | 2010-10-19 | Guardian Technologies International, Inc. | System and method for identifying feature of interest in hyperspectral data |
Family Cites Families (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US73518A (en) * | 1868-01-21 | Luke fitzpatrick and jacob schinneller | ||
US7394363B1 (en) * | 1998-05-12 | 2008-07-01 | Bahador Ghahramani | Intelligent multi purpose early warning system for shipping containers, components therefor and methods of making the same |
CA2348150C (en) * | 2000-05-25 | 2007-03-13 | Esam M.A. Hussein | Non-rotating x-ray system for three-dimensional, three-parameter imaging |
US20020186875A1 (en) * | 2001-04-09 | 2002-12-12 | Burmer Glenna C. | Computer methods for image pattern recognition in organic material |
US7444309B2 (en) * | 2001-10-31 | 2008-10-28 | Icosystem Corporation | Method and system for implementing evolutionary algorithms |
US7356115B2 (en) * | 2002-12-04 | 2008-04-08 | Varian Medical Systems Technology, Inc. | Radiation scanning units including a movable platform |
WO2005024845A2 (en) * | 2003-04-08 | 2005-03-17 | Lawrence Berkeley National Laboratory | Detecting special nuclear materials in containers using high-energy gamma rays emitted by fission products |
US7856081B2 (en) * | 2003-09-15 | 2010-12-21 | Rapiscan Systems, Inc. | Methods and systems for rapid detection of concealed objects using fluorescence |
WO2005084352A2 (en) * | 2004-03-01 | 2005-09-15 | Varian Medical Systems Technologies, Inc. | Dual energy radiation scanning of objects |
US7847260B2 (en) * | 2005-02-04 | 2010-12-07 | Dan Inbar | Nuclear threat detection |
US7336767B1 (en) * | 2005-03-08 | 2008-02-26 | Khai Minh Le | Back-scattered X-ray radiation attenuation method and apparatus |
CA2608119A1 (en) * | 2005-05-11 | 2006-11-16 | Optosecurity Inc. | Method and system for screening luggage items, cargo containers or persons |
US7545907B2 (en) * | 2005-11-09 | 2009-06-09 | Dexela Limited | Methods and apparatus for obtaining low-dose imaging |
US7536365B2 (en) * | 2005-12-08 | 2009-05-19 | Northrop Grumman Corporation | Hybrid architecture for acquisition, recognition, and fusion |
US20070211248A1 (en) * | 2006-01-17 | 2007-09-13 | Innovative American Technology, Inc. | Advanced pattern recognition systems for spectral analysis |
US7483511B2 (en) * | 2006-06-06 | 2009-01-27 | Ge Homeland Protection, Inc. | Inspection system and method |
US8015127B2 (en) * | 2006-09-12 | 2011-09-06 | New York University | System, method, and computer-accessible medium for providing a multi-objective evolutionary optimization of agent-based models |
EP2074413B1 (en) * | 2006-10-25 | 2017-12-20 | Soreq Nuclear Research Center | Method for detecting nitrogenous materials via gamma-resonance absorption (gra) |
US7492862B2 (en) * | 2007-01-17 | 2009-02-17 | Ge Homeland Protection, Inc. | Computed tomography cargo inspection system and method |
US7706502B2 (en) * | 2007-05-31 | 2010-04-27 | Morpho Detection, Inc. | Cargo container inspection system and apparatus |
-
2008
- 2008-05-29 US US12/129,036 patent/US20090003651A1/en not_active Abandoned
- 2008-05-29 US US12/129,055 patent/US20090052622A1/en not_active Abandoned
- 2008-05-29 US US12/129,371 patent/US20090052762A1/en not_active Abandoned
- 2008-05-29 US US12/129,383 patent/US8094874B2/en not_active Expired - Fee Related
- 2008-05-29 US US12/129,393 patent/US20090055344A1/en not_active Abandoned
- 2008-05-29 US US12/129,439 patent/US20080298544A1/en not_active Abandoned
- 2008-05-29 US US12/129,410 patent/US20090003699A1/en not_active Abandoned
Patent Citations (46)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US538758A (en) * | 1895-05-07 | Richard watkins | ||
US4794257A (en) * | 1984-05-14 | 1988-12-27 | Matsushita Electric Industrial Co., Ltd. | Quantum-counting radiography |
US5132998A (en) * | 1989-03-03 | 1992-07-21 | Matsushita Electric Industrial Co., Ltd. | Radiographic image processing method and photographic imaging apparatus therefor |
US5319547A (en) * | 1990-08-10 | 1994-06-07 | Vivid Technologies, Inc. | Device and method for inspection of baggage and other objects |
US5490218A (en) * | 1990-08-10 | 1996-02-06 | Vivid Technologies, Inc. | Device and method for inspection of baggage and other objects |
US5600303A (en) * | 1993-01-15 | 1997-02-04 | Technology International Incorporated | Detection of concealed explosives and contraband |
US5600700A (en) * | 1995-09-25 | 1997-02-04 | Vivid Technologies, Inc. | Detecting explosives or other contraband by employing transmitted and scattered X-rays |
US5642393A (en) * | 1995-09-26 | 1997-06-24 | Vivid Technologies, Inc. | Detecting contraband by employing interactive multiprobe tomography |
US6018562A (en) * | 1995-11-13 | 2000-01-25 | The United States Of America As Represented By The Secretary Of The Army | Apparatus and method for automatic recognition of concealed objects using multiple energy computed tomography |
US6026171A (en) * | 1998-02-11 | 2000-02-15 | Analogic Corporation | Apparatus and method for detection of liquids in computed tomography data |
US6236709B1 (en) * | 1998-05-04 | 2001-05-22 | Ensco, Inc. | Continuous high speed tomographic imaging system and method |
US6282305B1 (en) * | 1998-06-05 | 2001-08-28 | Arch Development Corporation | Method and system for the computerized assessment of breast cancer risk |
US6567496B1 (en) * | 1999-10-14 | 2003-05-20 | Sychev Boris S | Cargo inspection apparatus and process |
US20010033636A1 (en) * | 1999-11-13 | 2001-10-25 | Martin Hartick | Method and apparatus for determining a material of a detected item |
US6969861B2 (en) * | 2001-10-02 | 2005-11-29 | Konica Corporation | Cassette for radiographic imaging, radiographic image reading apparatus and radiographic image reading method |
US20050111619A1 (en) * | 2002-02-06 | 2005-05-26 | L-3 Communications Security And Detection Systems Corporation Delaware | Method and apparatus for target transmitting information about a target object between a prescanner and a CT scanner |
US7023957B2 (en) * | 2002-02-06 | 2006-04-04 | L-3 Communications Security And Detection Systems, Inc. | Method and apparatus for transmitting information about a target object between a prescanner and a CT scanner |
US7123762B2 (en) * | 2002-02-08 | 2006-10-17 | University Of Chicago | Method and system for risk-modulated diagnosis of disease |
US7162005B2 (en) * | 2002-07-19 | 2007-01-09 | Varian Medical Systems Technologies, Inc. | Radiation sources and compact radiation scanning systems |
US7103137B2 (en) * | 2002-07-24 | 2006-09-05 | Varian Medical Systems Technology, Inc. | Radiation scanning of objects for contraband |
US7130371B2 (en) * | 2002-09-27 | 2006-10-31 | Scantech Holdings, Llc | System for alternately pulsing energy of accelerated electrons bombarding a conversion target |
US20050025280A1 (en) * | 2002-12-10 | 2005-02-03 | Robert Schulte | Volumetric 3D x-ray imaging system for baggage inspection including the detection of explosives |
US7092485B2 (en) * | 2003-05-27 | 2006-08-15 | Control Screening, Llc | X-ray inspection system for detecting explosives and other contraband |
US20040247075A1 (en) * | 2003-06-06 | 2004-12-09 | Johnson James H. | Vehicle mounted inspection systems and methods |
US6937692B2 (en) * | 2003-06-06 | 2005-08-30 | Varian Medical Systems Technologies, Inc. | Vehicle mounted inspection systems and methods |
US20060269114A1 (en) * | 2003-07-03 | 2006-11-30 | General Electric Company | Methods and systems for prescribing parameters for tomosynthesis |
US20050002550A1 (en) * | 2003-07-03 | 2005-01-06 | Ge Medical Systems Global Technology Company, Llc | Imaging chain for digital tomosynthesis on a flat panel detector |
US20050031075A1 (en) * | 2003-08-07 | 2005-02-10 | Hopkins Forrest Frank | System and method for detecting an object |
US7491958B2 (en) * | 2003-08-27 | 2009-02-17 | Scantech Holdings, Llc | Radiographic inspection system for inspecting the contents of a container having dual injector and dual accelerating section |
US20060098773A1 (en) * | 2003-09-15 | 2006-05-11 | Peschmann Kristian R | Methods and systems for rapid detection of concealed objects using fluorescence |
US20050058242A1 (en) * | 2003-09-15 | 2005-03-17 | Peschmann Kristian R. | Methods and systems for the rapid detection of concealed objects |
US7162007B2 (en) * | 2004-02-06 | 2007-01-09 | Elyan Vladimir V | Non-intrusive inspection systems for large container screening and inspection |
US20050180542A1 (en) * | 2004-02-17 | 2005-08-18 | General Electric Company | CT-Guided system and method for analyzing regions of interest for contraband detection |
US20050256820A1 (en) * | 2004-05-14 | 2005-11-17 | Lockheed Martin Corporation | Cognitive arbitration system |
US7190757B2 (en) * | 2004-05-21 | 2007-03-13 | Analogic Corporation | Method of and system for computing effective atomic number images in multi-energy computed tomography |
US7817833B2 (en) * | 2004-05-26 | 2010-10-19 | Guardian Technologies International, Inc. | System and method for identifying feature of interest in hyperspectral data |
US20060233302A1 (en) * | 2004-10-22 | 2006-10-19 | Might Matthew B | Angled-beam detection system for container inspection |
US20070248212A1 (en) * | 2004-10-22 | 2007-10-25 | Might Matthew B | Cryptographic container security system |
US7356118B2 (en) * | 2004-10-22 | 2008-04-08 | Scantech Holdings, Llc | Angled-beam detection system for container inspection |
US20060256914A1 (en) * | 2004-11-12 | 2006-11-16 | Might Matthew B | Non-intrusive container inspection system using forward-scattered radiation |
US20060204107A1 (en) * | 2005-03-04 | 2006-09-14 | Lockheed Martin Corporation | Object recognition system using dynamic length genetic training |
US20070009084A1 (en) * | 2005-06-01 | 2007-01-11 | Endicott Interconnect Technologies, Inc. | Imaging inspection apparatus with directional cooling |
US7399976B2 (en) * | 2005-11-03 | 2008-07-15 | Tsinghua University | Method and equipment for discriminating materials by employing fast neutron and continuous spectral x-ray |
US20090003699A1 (en) * | 2007-05-29 | 2009-01-01 | Peter Dugan | User guided object segmentation recognition |
US20080298544A1 (en) * | 2007-05-29 | 2008-12-04 | Peter Dugan | Genetic tuning of coefficients in a threat detection system |
US20090055344A1 (en) * | 2007-05-29 | 2009-02-26 | Peter Dugan | System and method for arbitrating outputs from a plurality of threat analysis systems |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120170826A1 (en) * | 2009-09-22 | 2012-07-05 | Kwang Eun Jang | Multi-energy x-ray system, image processing apparatus for discriminating multi-energy x-ray material, and image processing method for material discrimination of multi-energy x-ray system |
US8837801B2 (en) * | 2009-09-22 | 2014-09-16 | Samsung Electronics Co., Ltd. | Multi-energy X-ray system, image processing apparatus for discriminating multi-energy X-ray material, and image processing method for material discrimination of multi-energy X-ray system |
US20120069963A1 (en) * | 2010-08-06 | 2012-03-22 | Telesecurity Sciences, Inc. | Dual energy backscatter x-ray shoe scanning device |
US9036782B2 (en) * | 2010-08-06 | 2015-05-19 | Telesecurity Sciences, Inc. | Dual energy backscatter X-ray shoe scanning device |
US10216866B2 (en) * | 2011-02-25 | 2019-02-26 | Smiths Heimann Gmbh | Image reconstruction based on parametric models |
US20140140601A1 (en) * | 2011-06-30 | 2014-05-22 | Andrew Litvin | Iterative image reconstruction |
US9476923B2 (en) | 2011-06-30 | 2016-10-25 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Method and device for identifying a material by the spectral analysis of electromagnetic radiation passing through said material |
US9824467B2 (en) * | 2011-06-30 | 2017-11-21 | Analogic Corporation | Iterative image reconstruction |
US10078150B2 (en) | 2015-04-14 | 2018-09-18 | Board Of Regents, The University Of Texas System | Detecting and quantifying materials in containers utilizing an inverse algorithm with adaptive regularization |
US10379252B2 (en) | 2015-07-22 | 2019-08-13 | Tsinghua University | Methods and devices for estimating weight of an object to be inspected in an inspection system |
Also Published As
Publication number | Publication date |
---|---|
US20090003699A1 (en) | 2009-01-01 |
US20080298544A1 (en) | 2008-12-04 |
US20090055344A1 (en) | 2009-02-26 |
US8094874B2 (en) | 2012-01-10 |
US20090003651A1 (en) | 2009-01-01 |
US20090052622A1 (en) | 2009-02-26 |
US20090052732A1 (en) | 2009-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20090052762A1 (en) | Multi-energy radiographic system for estimating effective atomic number using multiple ratios | |
DK1970731T3 (en) | Method for calibrating a dual-energy computed tomography system and image reconstruction method | |
Rivetti et al. | Comparison of different commercial FFDM units by means of physical characterization and contrast‐detail analysis | |
JP6456699B2 (en) | X-ray photon counting computed tomography apparatus, spectral correction method, and spectral correction program | |
Mackenzie et al. | Conversion of mammographic images to appear with the noise and sharpness characteristics of a different detector and x‐ray system | |
US10379252B2 (en) | Methods and devices for estimating weight of an object to be inspected in an inspection system | |
US20090304249A1 (en) | Material composition detection from effective atomic number computation | |
US9589373B2 (en) | Monte carlo modeling of field angle-dependent spectra for radiographic imaging systems | |
Kappadath et al. | Dual‐energy digital mammography for calcification imaging: Scatter and nonuniformity corrections | |
JP2015525360A (en) | Dynamic modeling of imperfections for photon counting detectors | |
WO2019235087A1 (en) | Radiography device, radiography method, and program | |
Osipov et al. | Parameter identification method for dual-energy X-ray imaging | |
Howard et al. | Sampling-based uncertainty quantification in deconvolution of X-ray radiographs | |
US20180164231A1 (en) | Dual-energy microfocus radiographic imaging method for meat inspection | |
Kuhls‐Gilcrist et al. | Accurate MTF measurement in digital radiography using noise response | |
Rink et al. | Investigating the feasibility of photon‐counting K‐edge imaging at high x‐ray fluxes using nonlinearity corrections | |
JP2014176565A (en) | Image processor, radiographic apparatus, image processing method, computer program and recording medium | |
Wang et al. | Locally linear transform based three‐dimensional gradient‐norm minimization for spectral CT reconstruction | |
Choi et al. | A unified statistical framework for material decomposition using multienergy photon counting x‐ray detectors | |
Ullherr et al. | SNR spectra as a quantitative model for image quality in polychromatic X-ray imaging | |
Anam et al. | Impact of ROI size on the accuracy of noise measurement in CT on computational and ACR phantoms | |
Ghammraoui et al. | Inclusion of a GaAs detector model in the Photon Counting Toolkit software for the study of breast imaging systems | |
Marchal | Extension of x‐ray imaging linear systems analysis to detectors with energy discrimination capability | |
Lee et al. | Nonlinear lag correction based on the autoregressive model for dynamic flat-panel detectors | |
Goodman et al. | Angular Detection and Shielding Characterization of Simulated 235U Using Time Encoded Imaging and 3D Position Sensitive CdZnTe Detectors |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: LOCKHEED MARTIN CORPORATION, MARYLAND Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:DUGAN, PETER;FINCH, ROBERT L.;MUNLEY, JOHN M.;AND OTHERS;REEL/FRAME:021809/0851;SIGNING DATES FROM 20080627 TO 20081105 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |