WO2009015563A1 - Procédé et système d'identification de matériau par le biais d'images binoculaires de transmission d'énergie multiple - Google Patents

Procédé et système d'identification de matériau par le biais d'images binoculaires de transmission d'énergie multiple Download PDF

Info

Publication number
WO2009015563A1
WO2009015563A1 PCT/CN2008/001418 CN2008001418W WO2009015563A1 WO 2009015563 A1 WO2009015563 A1 WO 2009015563A1 CN 2008001418 W CN2008001418 W CN 2008001418W WO 2009015563 A1 WO2009015563 A1 WO 2009015563A1
Authority
WO
WIPO (PCT)
Prior art keywords
energy
image
grayscale
fault
tomogram
Prior art date
Application number
PCT/CN2008/001418
Other languages
English (en)
French (fr)
Chinese (zh)
Inventor
Yali Xie
Qitian Miao
Hua Peng
Kejun Kang
Haifeng Hu
Zhiqiang Chen
Xueguang Cao
Chuanxiang Tang
Jianping Gu
Xuewu Wang
Hongsheng Wen
Bei He
Yaohong Liu
Shangmin Sun
Quanwei Song
Jin Lin
Xianli Ding
Original Assignee
Nuctech Company Limited
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Priority claimed from CN2008100813251A external-priority patent/CN101358936B/zh
Application filed by Nuctech Company Limited, Tsinghua University filed Critical Nuctech Company Limited
Publication of WO2009015563A1 publication Critical patent/WO2009015563A1/zh

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/55Depth or shape recovery from multiple images
    • G06T7/593Depth or shape recovery from multiple images from stereo images
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N13/00Stereoscopic video systems; Multi-view video systems; Details thereof
    • H04N13/20Image signal generators
    • H04N13/204Image signal generators using stereoscopic image cameras
    • H04N13/239Image signal generators using stereoscopic image cameras using two 2D image sensors having a relative position equal to or related to the interocular distance
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating 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/02Investigating 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/04Investigating 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 forming images of the material
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating 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/02Investigating 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/06Investigating 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
    • G01N23/083Investigating 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 the radiation being X-rays
    • G01N23/087Investigating 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 the radiation being X-rays using polyenergetic X-rays
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V5/00Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
    • G01V5/20Detecting prohibited goods, e.g. weapons, explosives, hazardous substances, contraband or smuggled objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N13/00Stereoscopic video systems; Multi-view video systems; Details thereof
    • H04N13/20Image signal generators
    • H04N13/204Image signal generators using stereoscopic image cameras
    • H04N13/254Image signal generators using stereoscopic image cameras in combination with electromagnetic radiation sources for illuminating objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection

Definitions

  • the present invention relates to the field of radiation imaging technology, and more particularly to a method for material identification using dual-view multiple-energy transmission images. Background technique
  • Radiation imaging technology enables the transmission of the object's transmission image in a non-contact state by means of the high-energy X-ray penetration capability.
  • the working principle of scanning radiation imaging uses X-rays emitted by the radiation source, and the X-rays are received by the detector after being received by the detector, and then converted into an electrical signal input image acquisition system, by image The acquisition system inputs an image signal into a computer display to display the detected image.
  • the transmission image formed by radiation imaging is a projection of all objects penetrated by the X-ray beam stream, and does not include information in the depth direction.
  • the scanned image is an image formed by superimposing all objects along the direction of the scanning beam, which is not conducive to checking the items hidden behind other objects, and also cannot identify the material of the transmitted image.
  • a more mature main body reconstruction technology uses computed tomography, but the disadvantages of this technology are: not only complicated equipment, high cost, but also unable to quickly check large objects. The pass rate is low, and again, the material of the transmissive object cannot be identified.
  • the dual-view radiation transmission image processing technique is a radiation imaging technique that can detach an object located at a different depth of the slice in the detection space from the image, thereby removing the occlusion. This technique is used to strip certain overlapping objects in the transmitted image so that the obscured object is more visible in the image, but the material properties cannot be identified.
  • the multi-energy transmission image recognition technology is a radiation transmission image processing technique that identifies properties of materials such as organic substances, mixtures, metals, etc., by utilizing the characteristics that different materials have different ray attenuation capabilities for different energies.
  • this technique can only identify the properties of the object that accounts for the main part of the attenuation. If the object to be identified absorbs only a small fraction of the total attenuation dose, it is not possible to use this technique to identify the properties of such an object.
  • the present invention provides a scanning radiation identification-imaging method for a large object radiation inspection system having a simple structure.
  • the method is a material identification method for combining a dual-view technique and a multi-energy transmission image technique for performing a transmission image.
  • the method firstly detects the space by the dual-view technique.
  • the tomographic template of the intermediate object along the depth direction is then reconstructed from the transmission image, and finally the multi-energy technique is used to identify the object with successful gray reconstruction in the fault.
  • a method for material identification using a dual-view multi-energy transmission image includes the steps of: passing two X-rays at a certain angle through an object to be inspected, thereby obtaining left and right transmission image data, and transmitting the left and right directions The image is segmented, and the results of the segmentation are matched;
  • Material identification is performed on objects with successful grayscale reconstruction in the fault.
  • the invention provides a method for material identification by using dual-view multi-energy transmission images, by which an object covering the main component of the ray absorption can be removed for an object overlapping in the direction of the beam, so that the original Objects that are not apparent due to the secondary component of the ray absorption become apparent and can be identified by their material properties: such as organic matter, mixtures, metals, etc.).
  • the method according to the invention makes it possible to identify the non-primary components in the direction of the radiation, which lays the foundation for the automatic identification of explosives, drugs and other harmful objects contained in the container.
  • FIG. 1 shows a schematic diagram of detecting a spatial fault according to the dual view technique of the present invention
  • FIG. 2 shows a specific flow chart of image segmentation in accordance with the present invention
  • Figure 3 shows an example of a template map of a matched fault according to the present invention.
  • Figure 5 shows an example of a left view obtained by the method according to the invention and its greyscale reconstruction
  • Figure 6 shows a flow of grayscale grayscale reconstruction including different energies in accordance with the present invention
  • Figure 7 shows a detailed flow of grayscale reconstruction for a single energy in accordance with the present invention.
  • Figure 8 is a flow chart showing the material identification of an object in any layer according to the present invention.
  • Figure 9 is a multi-energy recognition effect diagram for a plurality of materials without overlap;
  • Figure 10 shows a non-peeling recognition effect diagram
  • Figure 11 is a diagram showing the multi-energy material recognition effect of the double-view occlusion peeling process.
  • Figure 12 is a block diagram of a dual-view multi-energy scanning radiation imaging system in accordance with the present invention.
  • FIGS. 13A and 13B are top views of a schematic layout of a dual-view multi-energy scanning radiation imaging system in accordance with the present invention.
  • Figure 14 is a side elevational view showing the arrangement of a dual angle scanning multi-energy radiation imaging system in accordance with the present invention. detailed description
  • a specific implementation of the method of dual-view multi-energy transmission image material identification according to the present invention is described in three parts hereinafter. I. Using a dual-view processing technique for each energy double-view image to obtain a tomogram template corresponding to the transmission image of the energy, and combining the fault templates of different energies into a template of a set of tomographic images.
  • FIG. 1 is a schematic illustration of the detection of spatial faults in accordance with the dual viewing angle technique of the present invention.
  • the detection space is a three-dimensional space formed by a region scanned by a sector of a radiation source (coordinate origin 0) to a vertical array of detectors and ?: where 0 ⁇ Coordinate origin; S: source position; left detector array; p: right detector array; 0L left beam; OR: right beam; ⁇ : left and right beam angle.
  • the scan direction (vertical upward) is the positive direction of the f-axis
  • the coordinate value is the scan ordinal number
  • the direction of the detector arrangement (perpendicular to the outside of the paper) is the positive axis
  • the coordinate value is the detector ordinal number to the horizontal To the right is the Z-axis positive direction
  • the coordinate value is the fault ordinal number.
  • the space rectangular coordinate system is established with the position 0 of the ray source as the coordinate origin, and we call the space where the angular coordinate system is located as the detection space.
  • the fault is a series of spatial planes parallel to the X- 0-y plane, in the figure
  • the dotted line shows the projection of the plane of the fault on the plane, and the depth of each layer is the distance between the plane of the layer and the plane.
  • A is the projection of the main beam direction of the left beam on the - ⁇ z plane, /? is right
  • 0 is the angle formed by the projection of the left and right beams on the - plane.
  • a template map of the detected space object can be obtained by the image edge extraction technique. That is, several edges are obtained by detecting local discontinuities, and then they are connected.
  • the edge extraction method is reliable in X-ray transmission image segmentation due to the inherent characteristics of the X-ray transmission image when the objects overlap.
  • the present invention simultaneously uses Sobel and Canny edge detection operators to extract edges, and then combines the edges detected by the two detection operators.
  • the edge obtained by the edge is connected to form a closed area, thereby realizing the division of the left and right views.
  • Figure 2 shows a detailed flow chart of image segmentation in accordance with the present invention.
  • edges are extracted.
  • Sobel and Canny edge detection operators are simultaneously used to extract edges.
  • represents the pixel gray value of the i-th column and the j-th row in the image
  • ⁇ ⁇ ' ⁇ is a set of all pixel points of the image
  • the present invention uses the Sobel edge detection operator for the number
  • Each pixel of the image [inspects the weighted difference of the gray levels of the upper, lower, left and right neighbors, the neighboring point weight is significant, and the next nearest neighboring point weight is small. Its definition is as follows:
  • , ⁇ ⁇ ⁇ 1 are the convolution operator
  • is the convolution of the i-th column and the j-th row
  • the matrix definition of the convolution operator is Then select the threshold 73 ⁇ 4, and the point ( ) that meets the condition s(/, )>73 ⁇ 4 is the step edge point, which is the edge image.
  • the general steps of the Canny edge detection algorithm are: first smoothing the image with a Gaussian filter; then using the finite difference of the first-order partial derivative to calculate the amplitude and direction of the gradient; non-maximum suppression of the gradient amplitude; and finally using the double threshold
  • the algorithm detects and joins edges.
  • the Canny operator uses a double threshold algorithm to reduce false edges. That is, the non-maximum suppression image is binarized with two thresholds A and 7, ⁇ 2 ⁇ , to obtain two threshold edge images N, and N 2 (i,j)o N 2 (, ) is obtained using a high threshold. Contains few false edges, but with discontinuities. Then connect the edges into outlines in N 2 ( , _ ).
  • the algorithm finds the edges that can be joined at the endpoints of the edges in 7 '), and the algorithm continuously finds the edges in the corresponding 8 neighborhoods, so that the algorithm continuously collects the edges until ⁇ ('',') is connected to form the contour. .
  • a closed edge is obtained.
  • the edges detected by the Sobel and Canny edge detection operators are combined and edge-joined to close them.
  • the preliminary edge map of the present invention is the result of finding a logical "or" for the binary edge image obtained by the above two operators. Due to the influence of noise and the like, the edges obtained by the above methods are generally still discontinuous, so they need to be connected.
  • the invention connects the edge pixels according to their continuity in the gradient amplitude or gradient direction. For example, if the pixel ( ⁇ t) is in the neighborhood of the pixel (x, and their gradient magnitude and gradient direction are satisfied at a given threshold:
  • step 03 the segmentation of the view is performed according to the obtained closed edge. Since the closed edge divides the image into two areas, the morphological expansion-corrosion operation can be used to find points that belong to the interior of the area. Then, the region growth method can be used to fill the pixels in the region with the value "1", and the pixels outside the region are filled with the value "0", and the binary template of the inner region of each closed edge is obtained.
  • the size of the map is equal to the projection of the detection space on the plane, ie the number of scans (width) of the number of X detectors (height). This completes the image segmentation and gets the object template.
  • objects on the two template maps are matched according to a certain rule by a dual-view technique. That is, the connected area with each internal value of 1 in the left template picture is compared with each template in the right template picture one by one, and the corresponding template of the same object in the right view is found. Thus, each object that matches successfully corresponds to a template in the left and right views.
  • the difference in position between the two templates in the horizontal direction is the parallax ⁇ ".
  • ⁇ , and / ⁇ are the position of the center of gravity of the matching object in the tomographic template in the left and right views, and the parallax is proportional to the depth of each layer.
  • Figure 3 is a diagram showing an example of a template map of a matched fault in accordance with the present invention. See Figure 3, which shows the left and right view segmentation results.
  • Figure 3 (a) is a left view object template
  • Figure 3 (b) is a right view object template, the template of which can be seen as a rectangle.
  • the fault template of the transmission image obtained according to the dual-view tomography technique reflects the front-back position of the object in the depth direction in the detected space, and the geometry of the fault template reflects the shape contour of the object.
  • FIG. 4 is a flow chart showing a template diagram of each energy layer obtained in accordance with the present invention.
  • a set of transmission images containing the respective radiant energy is first established.
  • a double loop operation with nesting inside and outside is performed, and the inner and outer dashed boxes respectively represent the inner and outer double loops.
  • the inner loop is a template tomogram generation loop.
  • the system establishes a matching object set to match the object number as a loop variable, according to steps 01 to 03 described in (1), to create, match, and object the object in the dual-view transmission image of an energy. Parallax calculation.
  • a template of parallax approximation is merged into the same tomogram in a logical or logical manner to obtain a template tomogram set of matched objects in each fault under a certain radiant energy.
  • the outer loop is a loop for different radiant energies.
  • an inner loop is generated, a tomogram containing a template matching the good objects at each fault depth, and a transmission image generated under the next energy are selected.
  • the image repeats the above steps until the transmitted image of all energy is processed.
  • the template tomogram sets of different energies are layered in a logical or logical manner and are identified by a fault depth. Each layer contains a set of template tomograms of several object templates. 2. According to the fault template, the process of gray reconstruction of multiple energies is performed separately.
  • the fault template map obtained above only reflects the geometry of the object and its spatial position in the container. If material identification is to be performed, gray value reconstruction of multiple energies must also be performed. After gray value reconstruction, we can obtain the gray values of various energies of each segmented object. Thereby material identification of the object can be achieved.
  • grayscale reconstruction for each energy, we reconstruct the objects in each fault based on the two-view grayscale reconstruction method that is stripped from the outside to the inside. That is: first, the gray scale of the matched object in the X- 0-y plane is directly reconstructed (directly adjacent to the background area), and the background is still the original background value, and the object is curved in the outline of the object. The reconstructed grayscale image of the gray value to be obtained is scanned, and the object is stripped from the original image using the reconstructed grayscale. Then do the same for the outer layer of the object, repeat the above procedure until all the matched objects have been reconstructed.
  • Figure 5 (a) Taking the left view shown in Figure 5 (a) as an example, we perform grayscale reconstruction in combination with a template map similar to that shown in Figure 4.
  • three objects overlap each other, from the outside to the inside: a larger rectangle, a smaller rectangle, and a smaller ellipse.
  • Figures 5(b), 5(c) and 5(d) show the effect of grayscale reconstruction.
  • Figure 5 (b) is the outermost side
  • Figure 5 (c> is the middle
  • Figure 5 (d) is the innermost side.
  • Figure 5 (b) is the result of grayscale reconstruction of the outermost object, where the gray value of the light area is the background value in the original picture, and the gray value of the dark area is the gray area of the light area minus the outermost object reconstruction.
  • Gray value, dark area wheel The contour is the same as the contour in the template image of the object, and is a larger rectangle
  • Figure 5 (c) is the grayscale reconstruction result of the intermediate object, wherein the gray value of the light color region is the background value in the original image, and the dark region is The gray value is the gray area minus the intermediate object reconstructed gray value, and the dark area contour is the same as the contour of the object template, which is a smaller rectangle
  • Figure 5 (d) is the innermost object Gray-scale reconstruction results, wherein the gray value of the light-colored area is the background value in the original image, and the gray value of the dark-colored area is the gray-scale area minus the gray-scale area, and the dark area contour and the The template of the object has the same outline and is an elli
  • Figure 6 shows the flow of grayscale grayscale reconstruction including different energies.
  • the first weight is a fault set, which contains all template tomograms with the fault depth as the element distinguishing mark; the second weight is the object set, and the object serial number is used to distinguish the mark, including Each matched object at a specific fault depth; the third weight is an energy set, with different radiant energy as a distinguishing mark, and a transmission gradation reconstructed map at a specific energy of a specific matched object.
  • the above method of reconstructing the grayscale image of the object and completing the reconstruction of the grayscale image of an energy is to perform the grayscale reconstruction of the object by the method of peeling off from the outer to the inner and the grayscale.
  • the objects have all been completed.
  • the flow chart is shown in Figure 7.
  • FIG. 7 shows a detailed flow chart of grayscale reconstruction in accordance with the present invention. Referring to Figure 7, a detailed implementation process of the grayscale reconstruction according to the present invention will be described in detail.
  • the gray scale reconstruction of the object is performed by the method of sequentially peeling off from the outside to the inside and the gray scale.
  • step 01 a set of grayscale reconstruction candidate objects is created using the object obtained by the image segmentation process; in step 02, an attribute of an object is obtained.
  • step 03 it is confirmed whether the acquired object has an edge adjacent to the background.
  • step 04 the acquired object has a grayscale reconstructed object if it has an edge adjacent to the background; if the object has an object occlusion, the grayscale of the occlusion is reconstructed.
  • the object is erased in the image.
  • the gray scale of the reconstructed object be equal to the outer gray level of the edge minus the gray level of the edge side, that is,
  • gray values of various energies of the object represented by each template in each layer can be obtained. These gray values differ depending on the energy. By analyzing these differences, the material in any layer can be identified.
  • FIG. 8 shows the flow of the above identification process.
  • a "fault-object-energy" triple set obtained by the second part of the implementation process is introduced.
  • a double cycle of "fault-object” is then carried out.
  • the outer dotted line frame is a cycle for each fault depth.
  • the present invention creates a "object-energy" double set containing the matched objects of the fault and containing the grayscale reconstruction map for each energy of each object.
  • the object loop in the inner dashed box is performed, that is, the loop in the object with the object number is marked for all the objects in the matching graph in a certain fault.
  • Painted into a color tomogram in such a color tomogram: the identified object
  • the outline of the body is determined by the outline of the template, and the color of the fill in the outline is determined by the combination of the gray scale reconstruction result and the material recognition result. Determination of Color
  • the present invention will be described below in the related description of FIG. Then, create a new "object-energy" set for the next fault, and perform material identification and color map reconstruction until all faults have been processed the same.
  • all the color material recognition effect tomograms obtained above are integrated into a set of color material recognition effect tomograms of objects on different faults, which is the final result obtained by the dual-view and multi-energy image processing technology of the present invention. .
  • the multi-energy material identification method can identify the difference of the gray levels of different energy transmission maps for objects without overlap.
  • Orange is the identification color of organic or light materials
  • green is the identification color of light metals or mixtures
  • blue is the metal identification color.
  • the recognition effect diagram is as shown in Fig. 9.
  • FIG. 9 there is shown a multi-energy recognition effect diagram of graphite, aluminum, iron, lead, and polyethylene rectangular at a mass thickness of 30 g.
  • the left-to-right targets are arranged in the order of graphite, aluminum, iron, lead, and polyethylene.
  • the recognition order is: orange, green, blue, blue, and orange. The material is identified correctly.
  • Fig. 10 shows a non-peeling recognition effect diagram.
  • FIG. 10 there is a large rectangular steel plate obstruction, in the middle is a blocked cylinder containing petroleum liquefied gas, the left is a carton-mounted disc, and the right is a carton containing cigarettes. It can be seen that the identification of the petroleum liquefied gas in the cylinder and the cigarette in the right side of the carton basically recognizes the error, and a partial recognition error occurs in the carton loaded on the left side of the carton. If the double-view layering method is used to peel off the obstruction and then recognize it, the effect is as shown in Fig. 11.
  • Figure 11 for the dual-view occluded multi-energy material recognition effect diagram: where Figure 11 (a) is the occlusion; Figure 11 (b) is the identified object.
  • Figure 11 shows that after using the dual-view occlusion stripping multi-energy material identification technology, the steel sheet covering in Figure 11 (a) is recognized as blue, which is metal, and the carton-mounted disc in Figure 11 (b) is liquefied. Both LPG and carton cigarettes are identified as orange, organic, and the recognition results are correct.
  • Figure 12 is a block diagram of a dual viewing angle multi-energy scanning radiation imaging system in accordance with the present invention. As shown in Fig. 12, the dual-view multi-energy scanning radiation imaging system of the present invention includes the following devices:
  • Multi-energy radiation source 1 X-ray generator, capable of generating X-ray beams of different energies.
  • the beam current controller receives the X-rays emitted by the radiation source 1 and emits two symmetrical or asymmetrical X-rays at an angle.
  • the left detector array 4 receives X-rays of different energies emitted by the multi-energy radiation source and converts them into electrical signals for input to the left image acquisition system 6.
  • the right detector array 5 receives X-rays of different energies emitted by the multi-energy radiation source and converts them into electrical signals for input to the right image acquisition system 7.
  • the left image acquisition system 6 receives an electrical signal from the left detector array 4 and acquires left image data therefrom.
  • the right image acquisition system 7 receives the electrical signal from the right detector array 5 and obtains the right image data therefrom.
  • the computer processing system 8 receives the left and right image data from the left image acquisition system 6 and the right image acquisition system 7 and processes them, respectively displaying the image of the measured object in the computer display, and also displaying the different depths reconstructed by the two. Tomographic image of the fault.
  • the radiation source 1 emits two symmetrical or asymmetrical X-rays at an angle through the beam controller 2, and each X-ray beam passes through the object to be inspected 3 by the left detector array 4 and
  • the right detector array 5 receives; then converts the electrical signals into the left image acquisition system 6 and the right image acquisition system 7, respectively, and the image data in the left and right image acquisition systems 6 and 7 can be processed by the computer processing system 8 respectively.
  • the image of the measured object is displayed on the computer display, and the tomographic image of the depth at different depths is also displayed.
  • the dual-view multi-energy scanning radiation imaging system can respectively obtain a tomographic template map of a transmission image corresponding to the energy of the double-view image of each energy by using a dual-view processing technique, and merge the tomographic templates of different energies.
  • a template of a set of tomographic images according to the tomographic template, respectively, a gray reconstruction process of multiple energies is performed, and the reconstructed tomographic image is subjected to material identification of any layer.
  • DETAILED DESCRIPTION OF THE INVENTION The description of the dual-view multi-energy scanning radiation imaging method of the present invention is not repeated here.
  • a preferred embodiment of the invention utilizes a double-slit collimator as a beam current controller for beam current control of the radiation emitted by the source.
  • FIG. 13 and 14 are a plan view and a side view, respectively, showing a schematic view of the arrangement of the apparatus embodying the present invention.
  • 13A is a case where the beam ray is symmetrical
  • FIG. 13B is a case where the beam ray is asymmetrical.
  • the beam controller is provided with two collimating slits to form a symmetric or asymmetric X-ray of the radiation source with a certain angle of the beam, left and right detector arrays.
  • 4 and 5 respectively scan the object to be measured at a symmetrical angle to the beam fan defined by the collimating slit of the double-slit collimator, and transmit the respective electrical signals to the corresponding left and right images.
  • the system then performs image processing by computer processing system 8 to derive a tomographic image containing depth information.
  • this method can not only enter the main components of the transmitted image in the ray direction.
  • the traditional multi-energy material identification method can only identify the main components in the ray direction. For example, in the direction of the ray, a thicker steel plate overlaps with a bag of smaller drugs. If the traditional multi-energy material identification method is used, in the ray direction, only the steel plate can be recognized, and the drug cannot be identified. .
  • the material identification can be performed separately in each fault, and not only can It is also recognized that steel maple (the main component in the direction of the ray attenuates the ray is large), and it is also possible to identify the drug (the non-primary component in the ray direction is less attenuated by the ray).
  • This method is particularly useful for material identification of container transmission images. For the transmission scanning image of the container, because the container has a large thickness and a long radiation penetration distance, harmful objects such as explosives and drugs are often non-main components in the ray direction. Therefore, the method lays a foundation for automatically identifying harmful substances such as explosives and drugs in a container transmission scanned image.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Quality & Reliability (AREA)
  • Toxicology (AREA)
  • Electromagnetism (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
PCT/CN2008/001418 2007-08-02 2008-08-04 Procédé et système d'identification de matériau par le biais d'images binoculaires de transmission d'énergie multiple WO2009015563A1 (fr)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
CN200710119871.5 2007-08-02
CN200710119871 2007-08-02
CN200810081325.1 2008-02-25
CN2008100813251A CN101358936B (zh) 2007-08-02 2008-02-25 一种利用双视角多能量透射图像进行材料识别的方法及系统

Publications (1)

Publication Number Publication Date
WO2009015563A1 true WO2009015563A1 (fr) 2009-02-05

Family

ID=39767414

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2008/001418 WO2009015563A1 (fr) 2007-08-02 2008-08-04 Procédé et système d'identification de matériau par le biais d'images binoculaires de transmission d'énergie multiple

Country Status (4)

Country Link
DE (1) DE102008040952A1 (de)
FR (1) FR2919780B1 (de)
GB (1) GB2454047B (de)
WO (1) WO2009015563A1 (de)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104567758B (zh) 2013-10-29 2017-11-17 同方威视技术股份有限公司 立体成像系统及其方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0610084A2 (de) * 1993-02-05 1994-08-10 Max Robinson Von 3D-Bildsystem abgeleitete visuelle Darstellung von Information
US20030081720A1 (en) * 2001-10-31 2003-05-01 Swift David C. 3D stereoscopic X-ray system
GB2390005A (en) * 2002-06-17 2003-12-24 Royal Holloway University Of L Screening Apparatus
JP2007017304A (ja) * 2005-07-08 2007-01-25 Matsushita Electric Ind Co Ltd X線検査装置およびx線検査方法
CN101210895A (zh) * 2006-12-28 2008-07-02 清华同方威视技术股份有限公司 一种双视角扫描辐射成像的方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
GB9720658D0 (en) * 1997-09-29 1997-11-26 Univ Nottingham Trent Detecting, improving and charecterising material in 3-D space
US7031425B2 (en) * 2002-11-27 2006-04-18 Ge Medical Systems Global Technology Company, Llc Methods and apparatus for generating CT scout images
CA2574402A1 (en) * 2004-07-20 2006-01-26 William Awad System and method for detecting the presence of a threat in a package
CN1995993B (zh) 2005-12-31 2010-07-14 清华大学 一种利用多种能量辐射扫描物质的方法及其装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0610084A2 (de) * 1993-02-05 1994-08-10 Max Robinson Von 3D-Bildsystem abgeleitete visuelle Darstellung von Information
US20030081720A1 (en) * 2001-10-31 2003-05-01 Swift David C. 3D stereoscopic X-ray system
GB2390005A (en) * 2002-06-17 2003-12-24 Royal Holloway University Of L Screening Apparatus
JP2007017304A (ja) * 2005-07-08 2007-01-25 Matsushita Electric Ind Co Ltd X線検査装置およびx線検査方法
CN101210895A (zh) * 2006-12-28 2008-07-02 清华同方威视技术股份有限公司 一种双视角扫描辐射成像的方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
YUE W. ET AL.: "Automated 2-D Cephalometric Analysis on X-ray Images by a Model Based Approach", IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, vol. 53, no. 8, August 2006 (2006-08-01), pages 1615 - 1623 *

Also Published As

Publication number Publication date
FR2919780A1 (fr) 2009-02-06
FR2919780B1 (fr) 2017-09-08
GB2454047A (en) 2009-04-29
GB0814150D0 (en) 2008-09-10
DE102008040952A1 (de) 2009-02-19
GB2454047B (en) 2011-02-23

Similar Documents

Publication Publication Date Title
US8194953B2 (en) Method and system of material identification using binocular steroscopic and multi-energy transmission images
US7545906B2 (en) Method and system for binocular steroscopic scanning radiographic imaging
AU2002307053B2 (en) X-ray inspection system
EP2856430B1 (de) Bestimmung eines z-effektiven wertes für voxelset unter verwendung eines ct-dichtebildes und zerstreuter mehrfachenergiedaten
US8917927B2 (en) Portable backscatter advanced imaging technology scanner with automated target recognition
US20090245463A1 (en) Automatic material discrimination by using computer tomography
US8180139B2 (en) Method and system for inspection of containers
US8670522B2 (en) Stereo X-ray inspection apparatus and method for forming three-dimensional image through volume reconstruction of image acquired from the same
JP6487703B2 (ja) X線検査装置及びx線検査方法
US9414788B2 (en) X-ray imaging apparatus and method of generating stereoscopic image
JP2005530153A (ja) 三次元座標を取得するための立体x線撮像装置
US20160335756A1 (en) Inspection system
EP2869094A1 (de) Stereoskopische Abbildungssysteme und Verfahren
EP3460531B1 (de) Scanning-bildgebungssystem zur sicheren inspektion eines objekts und bildgebungsverfahren dafür
US8542893B2 (en) System and method of pulmonary emboli detection for computed tomography
CN207263664U (zh) 辐射检查系统
EP1918738A2 (de) Verfahren zur Bestimmung der Lage und Form eines Gepäckstücks in einem Gepäckumschlagcontainer durch Röntgenbildanalyse
WO2009015563A1 (fr) Procédé et système d'identification de matériau par le biais d'images binoculaires de transmission d'énergie multiple
Frosio et al. Optimized acquisition geometry for X-ray inspection
JP2015524061A (ja) 放射線撮影装置で生じた映像で板型物体を分離する方法と装置
CN107292888A (zh) 待测物中屏蔽区域的获取方法、装置以及检测系统
CN106910165A (zh) 修复原始ct投影数据的方法及装置、ct成像系统

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 08783607

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 08783607

Country of ref document: EP

Kind code of ref document: A1