CN101592579A - The method and the device that utilize multi-view X ray that explosive substances in luggage is surveyed automatically - Google Patents

The method and the device that utilize multi-view X ray that explosive substances in luggage is surveyed automatically Download PDF

Info

Publication number
CN101592579A
CN101592579A CNA2009100884957A CN200910088495A CN101592579A CN 101592579 A CN101592579 A CN 101592579A CN A2009100884957 A CNA2009100884957 A CN A2009100884957A CN 200910088495 A CN200910088495 A CN 200910088495A CN 101592579 A CN101592579 A CN 101592579A
Authority
CN
China
Prior art keywords
ray
explosive
luggage
density
visual angle
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.)
Granted
Application number
CNA2009100884957A
Other languages
Chinese (zh)
Other versions
CN101592579B (en
Inventor
杨立瑞
李宏伟
孔维武
杨桂文
李永清
王勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
First Research Institute of Ministry of Public Security
Beijing Zhongdun Anmin Analysis Technology Co Ltd
Original Assignee
First Research Institute of Ministry of Public Security
Beijing Zhongdun Anmin Analysis Technology Co Ltd
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
Application filed by First Research Institute of Ministry of Public Security, Beijing Zhongdun Anmin Analysis Technology Co Ltd filed Critical First Research Institute of Ministry of Public Security
Priority to CN2009100884957A priority Critical patent/CN101592579B/en
Publication of CN101592579A publication Critical patent/CN101592579A/en
Application granted granted Critical
Publication of CN101592579B publication Critical patent/CN101592579B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

A kind of method of utilizing multi-view X ray that explosive in the luggage is surveyed automatically, it at first gathers the picture signal at each visual angle of tested luggage and articles; Secondly each visual angle picture signal of obtaining is carried out material automatic detection, promptly the double energy X-ray image is cut apart, and the high energy after cutting apart and low energy image are carried out material characteristics identification according to double energy X-ray material discrimination principle, obtain meeting the hazardous location of explosive material characteristics, and mark the position that wherein is fit to participate in density calculation; The hazardous location segmentation result is carried out density automatic detection, promptly carry out the hazardous location coupling, find out each hazardous location of the same object of representative, rebuild by various visual angles then and calculate, obtain the density value of hazardous location according to matched rule; At last, according to the result of material automatic detection and density detection, whether decision reports to the police to these luggage and articles.The invention has the advantages that speed of detection is fast, characteristic matching is easy, detection accuracy height, detection efficiency height.

Description

The method and the device that utilize multi-view X ray that explosive substances in luggage is surveyed automatically
Technical field
The present invention relates to a kind of safety detection method and device, particularly relate to a kind of method and apparatus that adopts the X ray digital imaging technology that explosive in the luggage is surveyed.
Background technology
For the luggage safety inspection, X ray safety check technology is the contactless safety inspection technology of a class that is most widely used at present, and in recent years, X ray safety inspection technology has obtained fast development under the ordering about of demands such as aviation safety.
The early stage x-ray inspection equipment that uses mostly is the X ray technology that adopts the single-view monoergic, it can show and represent in the luggage diverse location variety classes material to the perspective view of the absorbing state of ray, thereby and then observes by the safety inspection personnel that each regional shape facility goes to judge whether to have known violated thing in this perspective view.
Continuous development along with the X ray technology, up to the present the double energy X-ray technology of material effective atomic number has appearred estimating, this kind technology utilizes certain material that the difference on the X ray absorptive character of different-energy is differentiated this material, makes can to distinguish by material behavior by the material that shape is distinguished.Two types double energy X-ray rays safety detection apparatus is arranged at present, a kind of is standard or conventional double energy X-ray baggage inspection apparatus, and it adopts a source two spies, accurate double energy X-ray, low pass is crossed speed and hand inspection, is mainly used in conventional hand-baggage and ships in the baggage check; Another kind of then be advanced double energy X-ray rays safety detection apparatus (AT machine), its double source is two to be visited, the real dual-energy X ray, and high pass is crossed speed and is surveyed automatically, is mainly used in the multistage baggage screening system of shipping automatically.
In order to satisfy the safety check demand of high detection accuracy, the X-CT technology has been introduced in the safety check field, and the CT type explosive rays safety detection apparatus of X-CT technology has appearred adopting, its groundwork process is: at first scan by a certain cross section to three-dimensional body, gather the data acquisition relevant with this three-dimensional body cross section structure, and then to these data acquisitions and carry out inverse operation according to certain mathematical principle and obtain and described three-dimensional body cross section structure parameter value one to one, finally recover the cross-sectional image of three-dimensional body at last from parameter value by display technique.The density of equipment resolution characteristic of this employing X-CT type of skill is strong, has higher explosive detection ability, but the data volume that it need be handled is very huge, the time that obtains image is also a lot of than other system head, manufacturing cost is also very high, and these have all restricted popularizing and using of CT type equipment.
CT type explosive rays safety detection apparatus at above-mentioned employing X-CT technology, the rays safety detection apparatus of multi-view X ray explosive auto-Detection Technology has appearred adopting at present, it mainly adopts a plurality of fixed viewpoint double energy X-rays to go material is shone, go to obtain the material characteristics of object then by the double energy X-ray picture signal, and the projected image that passes through various visual angles is rebuild the checking matter cross section, obtain the material density feature, checking matter is carried out more accurate rerum natura detection thereby reach.
For the multi-view X ray detection technique, it has all had application on medical science and industrial detection, but adopt the quick luggage safety inspection of the channel-type multi-view X ray technology of this kind multi-view X ray scanning technique but just just to grow up later in 2000, in this luggage safety inspection technical field, U.S. VIVID company is the producer that develops the automatic detecting devices of multi-view X ray explosive (MVT) the earliest, the three visual angle mode of operations that this MVT adopts 3 sources 3 to visit, on detection method, it adopts the coupling of directly extracting checking matter two-dimensional geometry feature and carrying out the three-dimensional geometry feature from each visual angle projection, after definite checking matter three-dimensional geometry feature, carry out density detection and materials detection (calculating effective atomic number) again, its weak point is that three-dimensional geometry Feature Extraction and coupling are difficult, and the more needs of object that meet three-dimensional geometry characteristic matching condition are investigated one by one, and on the various visual angles layout, at the bottom of adopting two, this MVT shines visual angle and a horizontal view angle, though the correlativity at each visual angle is less in passage, the three-dimensional reconstruction that helps checking matter, but this various visual angles layout type causes in the luggage of horizontal positioned article overlapping serious easily, thereby extracts checking matter feature difficulty in each visual angle projected image.
Though the five visual angle mode of operations that the automatic detecting devices EDtS of multi-view X ray explosive of Germany Heimann company adopts 3 sources 5 to visit, and irradiation line source and 1 top irradiation line source at the bottom of adopting 2 on the various visual angles layout, each end irradiation line source is corresponding 2 groups of detectors respectively, the top corresponding 1 group of detector of irradiation line source forms 5 visual angles, for the MVT of U.S. VIVID company, its various visual angles layout type can not cause the luggage and articles of horizontal positioned to overlap seriously, and carry out the mode of three-dimensional reconstruction by the projected image at 5 visual angles, its reconstruction redundance is bigger, blur level is less, but the mode of operation at 5 visual angles but causes characteristic matching difficulty in the detection process easily.
Summary of the invention
In view of the foregoing, one of the present invention purpose is to provide a kind of method of utilizing multi-view X ray that explosive in the luggage is surveyed automatically, this method is mainly surveyed the explosive in the luggage automatically by material and density two dimensional features, and may further comprise the steps:
1) picture signal at each visual angle of the tested luggage and articles of collection;
2) each visual angle picture signal of obtaining is carried out material automatic detection: at first picture signal is cut apart, and the high energy after cutting apart and low energy image are carried out material characteristics identification according to double energy X-ray material discrimination principle, obtain meeting the hazardous location of explosive material characteristics, and in this hazardous location, mark the position that wherein is fit to participate in density calculation;
3) above-mentioned hazardous location segmentation result is carried out density automatic detection: at first carry out the hazardous location coupling according to the matched rule of setting in advance, find out each hazardous location of the same object of representative in the multi-view image, rebuild by various visual angles then and calculate, obtain the density value of described hazardous location;
4) according to material automatic detection result and hazardous location density detection result, whether decision reports to the police to these luggage and articles.
Wherein, picture signal described in the step 1) is the double energy X-ray image that comprises high energy image signal and low energy image signal;
Also comprise in the step 3) corresponding explosive database is set up at each visual angle; And
Comprise in the step 4) that a plurality of visual angles density detection value is asked on average, if this mean value falls into explosive database middle-high density explosive district and low-density explosive district, then Dui Ying article are judged to explosive, and are reported to the police; And to falling into other zones of explosive database, then corresponding article are judged to safe product, will not report to the police.
Described explosive database is that horizontal ordinate is that checking matter atomic number, the longitudinal axis are the differentiation plane of the two dimension " effective atomic number-density " that density measurement surrounded, and comprises high density explosive thing and two districts of low-density explosive in this explosive database.
In addition, image segmentation specifically comprises step 2):
A) utilize edge detection operator that the double energy X-ray image is carried out Edge Gradient Feature
B) in the neighborhood of extraction edge feature, search background pixel and the foreground pixel that comprises background, and adopt background to reject formula and ask the foreground image that obtains through after the background rejecting, wherein,
Described background is rejected C T = G T + B × G A G B , G BBe background pixel, G T+BFor comprising the foreground pixel of background, G ABe the gray-scale value of double energy X-ray image under the air, and background pixel G BGray-scale value greater than the foreground pixel G that comprises background T+BGray-scale value, G TBe the foreground image after rejecting through background.
In addition, another object of the present invention also is to provide the device that a kind of characteristic matching is easy, detection accuracy is high and utilize multi-view X ray that explosive in the luggage is surveyed automatically.
This described device that utilizes multi-view X ray that explosive in the luggage is surveyed automatically comprise at the bottom of overall treatment computing machine, conveyor, transfer passage, top, left side are according to visual angle module, right side according at the bottom of visual angle module and the left side according to the visual angle module; Wherein, described overall treatment computing machine for be connected to top, described left side according to visual angle module, right side at the bottom of according at the bottom of visual angle module and the left side according to the visual angle module, and top, described left side according to visual angle module, right side at the bottom of according at the bottom of visual angle module and the left side according to the visual angle module for being installed in the described transfer passage, and described conveyor is for being controlled by described overall treatment computing machine; Top, described left side according to the visual angle module comprise the top, left side according to x-ray source and to should the top, left side according to L type first detector of x-ray source; At the bottom of the described right side according to the visual angle module comprise at the bottom of the right side according to x-ray source and at the bottom of should the right side according to L type second detector of x-ray source; At the bottom of the described left side according to the low energy X ray source, with at the bottom of this left side according to photograph sigmatron source at the bottom of corresponding L type the 3rd detector in low energy X ray source, the left side and with at the bottom of the described left side according to corresponding L type the 4th detector in sigmatron source.
Be compared to prior art, method and the device that utilizes multi-view X ray that explosive substances in luggage is surveyed automatically of the present invention has the following advantages:
(1) speed of detection is fast, that is, per hour detectable 1800 parcels are suitable for the luggage safety inspection field of high percent of pass;
(2) characteristic matching is easy, for the luggage of article situation more complicated, and its detection accuracy height, detection efficiency height.
Description of drawings
Fig. 1 is the method flow diagram that utilizes multi-view X ray that explosive in the luggage is surveyed automatically of the present invention;
Fig. 2 is that step 320 is calculated Y direction reconstruction polygon intersection point synoptic diagram among Fig. 1;
Fig. 3 rebuilds the coordinate system synoptic diagram of being set up for step 321 various visual angles among Fig. 1;
Fig. 4 is the various visual angles algebraic reconstruction model synoptic diagram that step 321 is determined among Fig. 1;
Fig. 5 is the designed various visual angles explosive database synoptic diagram of step 321 among Fig. 1;
Fig. 6 is the one group of original high energy image that utilizes one of method that multi-view X ray surveys automatically to explosive in luggage embodiment to gather of the present invention;
Fig. 7 is that result schematic diagram is extracted in the material automatic detection hazardous location of Fig. 6 correspondence;
Fig. 8 is the hazardous location matching result synoptic diagram of Fig. 7 correspondence;
Fig. 9 is the various visual angles reconstructed results synoptic diagram in two bodies cross sections among Fig. 7;
Figure 10 is the equivalent atom ordinal number and the density detection result schematic diagram of two object correspondences among Fig. 6;
Figure 11 is the final alarming result of two objects among Fig. 7;
Figure 12 is the structural representation that utilizes the device that multi-view X ray surveys automatically to explosive in the luggage of the present invention;
Figure 13 is the front elevational schematic of utilizing transfer passage and x-ray source in the device that multi-view X ray surveys automatically to explosive in the luggage, detector position of the present invention.
Embodiment
Come the device and method that utilizes multi-view X ray that explosive in the luggage is surveyed automatically of the present invention is described in further detail below in conjunction with accompanying drawing and specific embodiment.
The method overview flow chart of Fig. 1 for utilizing multi-view X ray that explosive in the luggage is surveyed automatically described in the present invention, it mainly comprises following 4 steps:
Step 30: the picture signal of gathering three visual angles;
Step 31: material automatic detection;
Step 32: density automatic detection;
Step 33: whether report to the police according to material characteristics and density feature decision.
Wherein, and in conjunction with shown in Figure 12 and Fig. 3, described three visual angles comprise the V1 visual angle that is made of according to x-ray source 1 and first detector 2 top, described left side respectively; Shine the V2 visual angle that the x-ray source 4 and second detector 3 constitute at the bottom of the described right side; At the bottom of the described left side according to low energy X ray source 9, with at the bottom of this left side according at the bottom of the 3rd detector 5 in low energy X ray source 9, the left side according to sigmatron source 10 and with shining sigmatron source 10 at the bottom of the described left side V3 visual angle of corresponding the 4th detector 6 formations; Described picture signal is the double energy X-ray image that comprises high energy image signal and low energy image signal.
In above-mentioned steps 30, luggage 12 is sent into passage 7 by conveyor 11, and in regular turn through shining sigmatron source 10 and the 4th detector 6 at the bottom of the left side, low energy X ray source, left side 9 and the 3rd detector 5, shine the x-ray source 4 and second detector 3 at the bottom of the right side, the top, left side is according to the x-ray source 1 and first detector 2, produce V3 respectively, V2, the dual energy images signal at V1 visual angle, and be sent to the operation of carrying out subsequent step 31 to 33 in the overall treatment computing machine 8, and then realize demonstration to automatic detection of explosive and result of detection.
Material automatic detection further comprises described in the step 31:
Step 310: the double energy X-ray image is carried out image segmentation, and high energy after cutting apart and low energy image are carried out material characteristics identification, and then obtain meeting the hazardous location of explosive material characteristics according to double energy X-ray material discrimination principle.
Step 311: set reliability criterion, the position of determining to be fit in the hazardous location do density calculation.
Wherein, for described step 310, at first, utilize edge detection operator to carry out Edge Gradient Feature the double energy X-ray image that collects, and in the neighborhood of edge feature, search background pixel G BWith the foreground pixel G that comprises background T+B, background pixel G here BGray-scale value greater than the foreground pixel value G that comprises background T+B, use following formula comprising the foreground pixel value G of background T+BSubtracting background pixel value G B, obtain cutting apart foreground image G through the quilt after the background rejecting T
G T = G T + B × G A G B , Wherein, G AIt is the gray-scale value under the air double energy X-ray image.
Secondly, according to double energy X-ray material discrimination principle high energy after cutting apart and low energy image are carried out material characteristics identification on this basis, obtain meeting the hazardous location of explosive material characteristics.
Described double energy X-ray material discrimination principle is: double energy X-ray shines on certain material, and is different with the ratio of low energy X ray uptake to the sigmatron uptake according to this material, differentiates different materials.
And for the hazardous location that extracts through step 310, under ideal state, the hazardous location section projection that is fit to do density calculation should have uniform gray level, do not blocked, have certain width by high z number material and ray can be penetrated the essential characteristic of wearing danger fully.
Therefore, in step 311, and determine to be fit in the hazardous location do the slice position of density calculation according to following reliability criterion: the section interior intensity changes mild, and gray scale is normal near the edge of section, the width of section is suitable, and the edge or the inside of section are not blocked by high z number material.In the practical operation, also can utilize the rule opposite with above-mentioned reliability criterion to remove to determine the slice position that is not suitable for doing density calculation in the hazardous location, so remaining section is exactly the section that is fit to participate in density calculation.
In addition, owing to whether checked article are reported to the police and need be decided in conjunction with the materials detection result and the density detection result of these article, therefore in that being carried out, each visual angle high energy image also needs density automatic detection is carried out in the hazardous location after material automatic detection obtains wherein hazardous location.
Density automatic detection further comprises described in the step 32:
Step 320: various visual angles hazardous location coupling, set matched rule, find out the hazardous location of the same object of representative in the different visual angles image.
Step 321: rebuild by various visual angles, obtain the density detection result of hazardous location.
Wherein, for step 320, after multi-view image is passing through step 31, all may comprise many groups hazardous location information in each visual angle image, in order to calculate the density result of hazardous location, step 32 at first needs one group of hazardous location of the same object of representative is correctly found out from all hazardous locations, and the present invention finishes the coupling work of hazardous location on the same group of three visual angles in step 320.
The coupling flow process is: at first travel through the one group of hazardous location data that obtains in the hazardous location, three visual angles by the hazardous location.Judging whether according to various array modes can the match is successful, finishes the matching operation of carrying out next group hazardous location after the judgement again, all finishes judgement up to all traversal hazardous locations.
In order to realize the various combinations of hazardous location in the different visual angles, need carry out the coupling work of 7 kinds of permutation and combination.Need to prove that because initial visual angle difference, the result of its coupling also can be different, need in actual applications therefore to consider that being is that mate at initial visual angle with which visual angle.
Suppose to represent one group of hazardous location information with the order of V3-V2-V1, then the method needs that carry out hazardous location traversal according to all combinations of hazardous location, three visual angles are at first respectively taken out a hazardous location from V3, V2, V1 visual angle, definition warning hazardous location is expressed as a, obtain one group of permutation and combination, mate as initial visual angle respectively according to V3, V2, V1 order all positions that are labeled as a in will making up then.
Suppose that V1, V2, V3 represent the hazardous location at three visual angles, with V1 as initial visual angle, V1 SBe V1 explosive area directions X origin coordinates, V1 EFor V1 explosive area directions X stops coordinate, V2 SBe V2 explosive area directions X origin coordinates, V2 EFor V2 explosive area directions X stops coordinate, V3 SBe V3 explosive area directions X origin coordinates, V3 EFor V3 explosive area directions X stops coordinate, T1 is the empirical constant value that a large amount of experiments of a process obtain, and T2 is another constant value greater than T1, judges whether that the condition of mating has three, is divided into combination 1, combination 2, combination 3, and their implication is as follows:
1) combination 1 principle:
Condition 1:abs (V1 S-V2 S)≤T1 and abs (V1 S-V3 S)≤T1
Condition 2:abs (V1 E-V2 E)≤T1 and abs (V1 E-V3 E)≤T1
Condition 3: after condition 1 and condition 2 are all set up, adopt the Y direction to rebuild polygon intersection point determination methods and judge whether satisfy condition in the ray scanning direction.If also satisfy, then the hazardous location of mating is just thought in this group hazardous location.
2) combination 2 principles:
Condition 1:abs (V1 S-V2 S)＞T1 and abs (V1 S-V2 S)≤T2 and abs (V1 S-V3 S)≤T1
Perhaps: abs (V1 S-V3 S)＞T1 and abs (V1 S-V3 S)≤T2 and abs (V1 S-V2 S)≤T1
Condition 2:abs (V1 E-V2 E)≤T1 and abs (V1 E-V3 E)≤T1
Condition 3: after condition 1 and condition 2 are all set up, adopt the Y direction to rebuild polygon intersection point determination methods and judge whether satisfy condition in the ray scanning direction.If also satisfy, then the hazardous location of mating is just thought in this group hazardous location.
3) combination 3 principles:
Condition 1:abs (V1 E-V2 E)＞T1 and abs (V1 E-V2 E)≤T2 and abs (V1 E-V3 E)≤T1
Perhaps: abs (V1 E-V3 E)＞T1 and abs (V1 E-V3 E)≤T2 and abs (V1 E-V2 E)≤T1 condition 2:
Abs (V1 S-V2 S)≤T1 and abs (V1 S-V3 S)≤T1
Condition 3: after condition 1 and condition 2 are all set up, adopt the Y direction to rebuild polygon intersection point determination methods and judge whether satisfy condition in the ray scanning direction.If also satisfy, then the hazardous location of mating is just thought in this group hazardous location.
The situation of combination 3 is similar to combination 2 principles, and the limit of just not satisfying less than T1 is to stop the limit.
It is to intersect according to three visual angle initial sums termination rays and checking matter that Y direction in the coupling combination condition of hazardous location is rebuild polygon intersection point determination methods, utilize the polygonal intersection point number that is surrounded to judge, if the intersection point number, is just thought Satisfying Matching Conditions more than or equal to 4.Referring to Fig. 2, each radiographic source and luggage 12 (checking matter) all have initial sum to stop two intersections (following corresponding to going up of hazardous location), six intersections at three visual angles can surround a polygon ABCDEF, if this polygon vertex number is more than or equal to 4, and all coordinate system effective range (being the passage effective range), Y direction matching condition is just satisfied in this group hazardous location so.
Step 321 is rebuild the density detection result who calculates the hazardous location by carrying out checking matter cross section multi-view angle three-dimensional.
In the step 321, the checking matter cross section is adopted the various visual angles algebraic reconstruction of the present invention's design to calculate and is obtained.The various visual angles algebraic reconstruction calculates to be finished according to the following steps:
1) the multi-view X ray sniffer that proposes according to the present invention is set up various visual angles and is rebuild coordinate system, channel cross-section is divided into some little squares with the width of a pixel, and be each cell initialize 0, set up coordinate system in mode shown in Figure 3 again, (V1x, V1y), (V2x, V2y), (V3x, V3y) respectively in the representative graph 4 top, left side according to x-ray source 1, right side at the bottom of according at the bottom of x-ray source 4, the left side according to the central coordinate of circle in sigmatron source 10.
2) determine various visual angles algebraic reconstruction model, various visual angles algebraic reconstruction model is seen Fig. 4, and the width of its ray is zero in this model, and the spacing between the ray is a pixel wide, weight factor wij is defined as i bar ray and j the length that pixel intersects, and pi is the projection value of i bar ray.
3) determine the zone that the various visual angles algebraic reconstruction calculates, obtain three visual angles with the polygonal region ABCDEF that the object cross section surrounds, see Fig. 2 according to the projected image at three visual angles.
4) according to the zone of various visual angles algebraic reconstruction model of determining and the calculating of various visual angles algebraic reconstruction, the various visual angles algebraic reconstruction calculates and carries out as follows:
f ( k + 1 ) ij = f ( k ) ij + &lambda; &times; p - &Sigma; n = 1 v ( w ijn &times; f ( k ) ijn ) &Sigma; n = 1 v w ijn &times; w ijn &times; w ijn k < t k f ( k ) ij = 0 , 1 N &Sigma; f ^ i , j > &beta; k f ( k ) i , j , k &GreaterEqual; t k
Wherein, k is an iterations, f (k) IjBe f after the k time iteration IjThe value of cell, f (k+1) IjBe f after the k+1 time iteration IjThe value of cell, λ is a coefficient of relaxation, P is the projection value of ray L correspondence, w IjnExpression ray L passes cell f IjLength, v is that ray L passes the sum that is included in the cell in the polygon ABCDEF, w IjBe that ray L passes cell f IjLength, t kFor carrying out the threshold value of quickening interative computation, α kBe the upper limit segmentation threshold of the k time iteration, β kBe the lower limit segmentation threshold of the k time iteration, N is for satisfying f (k) jThe pixel number of＞β (k) condition.
According to above formula, the various visual angles algebraic reconstruction calculates and comprises conventional iterative computation and quicken two iterative computation processes of iterative computation, utilizes and quickens interative computation, can shorten the used time of iterative computation, improves the efficient of iterative computation.
It is to be noted, it is that the density of hypothesis checking matter inside is to carry out under the prerequisite of even substantially matter that above various visual angles algebraic reconstruction calculates, be that the value on each pixel is all identical on the checking matter cross section, and other zone is not because exist checking matter, can think that these regional pixel values are 0, that is to say, the pixel value of last minute checking matter region of image be remarked pixel size equate on the occasion of, other parts are 0.
In step 321, a kind of multi-view X ray projection imaging explosive database design principle that the present invention relates to is: if X ray passes the thickness of checking matter is t, and resulting projection gray level value is G t, so can be according to G tDetermine the electron density of checking matter with t.Like this, it is that checking matter atomic number, the longitudinal axis are the differentiation plane of the two dimension " effective atomic number-density " that density measurement surrounded that multi-view X ray transmission imaging explosive database design becomes horizontal ordinate, as shown in Figure 5.According to the needs of explosive detection, two districts of high density explosive thing and low-density explosive are arranged in the database.Need to prove that in order to improve the accuracy of density detection, each visual angle all needs to set up corresponding explosive database, final result of detection is got the average of three visual angle probe values.
In step 33, the result according to step 32 hazardous location is surveyed is judged to explosive with the article that fall into explosive database middle-high density explosive district and low-density explosive district, is reported to the police; Be judged to safe product to falling into other regional corresponding article of explosive database, will not report to the police.
Shown in Figure 12 and Figure 13, of the present inventionly utilize device that multi-view X ray surveys automatically to explosive in the luggage for explosive in the luggage being surveyed automatically by material and density two dimensional features and multi-view X ray, it comprise at the bottom of overall treatment computing machine 8, conveyor 11, transfer passage 7, top, left side are according to visual angle module, right side according at the bottom of visual angle module and the left side according to the visual angle module; Wherein, top, described left side comprises that according to the visual angle module top, left side is according to the x-ray source 1 and first detector 2; Shine at the bottom of the described right side and shine the x-ray source 4 and second detector 3 at the bottom of the visual angle module comprises the right side; At the bottom of the described left side according to the visual angle module comprise at the bottom of the described left side according to low energy X ray source 9, with at the bottom of this left side according to shining sigmatron source 10 at the bottom of the 3rd corresponding detector 5 of low energy X ray source 9, the left side and corresponding the 4th detector 6 with shining sigmatron source 10 at the bottom of the described left side.
See Figure 13, top, above-mentioned left side is according to x-ray source 1, shine x-ray source 4 at the bottom of the right side, the end, left side, is according to sigmatron source 10, the end, left side is the different azimuth that lays respectively at transfer passage 7 according to low energy X ray source 9, promptly, on the positive apparent direction of described transfer passage 7, top, described left side is the left side top position that is arranged at transfer passage 7 according to x-ray source 1, shine x-ray source 4 at the bottom of the described right side for being arranged at the right side, bottom of transfer passage 7, at the bottom of the described left side according at the bottom of sigmatron source 10 and the left side according to low energy X ray source 9 for being arranged at the bottom left of transfer passage 7, thereby in the left side top of transfer passage 7, bottom, right side and three angles of left bottom constitute the layout pattern at one three visual angle.
Equally, push up according to shining low energy X ray source 9 at the bottom of x-ray source 1, the right side at the bottom of shining sigmatron source 10, left side at the bottom of x-ray source 4, the left side corresponding to above-mentioned left side, first detector 2, second detector 3, four detector 6, three detector 5 corresponding with it also is attached to the diverse location of described transfer passage 7 respectively.
Wherein, top, described left side constitutes the V1 visual angle according to the x-ray source 1 and first detector 2; Constitute the V2 visual angle according to the x-ray source 4 and second detector 3 at the bottom of the described right side; At the bottom of the described left side according to low energy X ray source 9, with at the bottom of this left side according at the bottom of the 3rd detector 5 in low energy X ray source 9, the left side according to sigmatron source 10 and with at the bottom of the described left side according to corresponding the 4th detector 6 in sigmatron source 10 formation V3 visual angles jointly.
Described first detector 2, second detector 3, the 4th detector 6 and the 3rd detector 5 are L type detector.
In order to embody effect of the present invention better, and in conjunction with shown in Figure 12 and Figure 13, now a kind of multi-view X ray that utilizes that relates among the present invention is described further the method that explosive in the luggage carries out fast automatic detection with one group of actual treatment data.
Fig. 6 has provided the three visual angle high energy view data that obtain through step 30, and wherein A and B represent two images of items respectively, and the hazardous location that process step 310 obtains as shown in Figure 7.Matching result for hazardous location in the step 320 is referring to shown in Fig. 8, then enter step 321, Fig. 8 is that a certain slice position of two groups of article is rebuild the schematic cross-section that calculates through various visual angles, utilize explosive database illustrated in Figure 7 that the atomic number and the density of two groups of article are compared again, and referring again to Fig. 8, the result who obtains is: second group of article B is after process step 31, though on material character near explosive, but through after the step 32, confirm its density characteristics and do not meet the requirement of explosive, be confirmed to be safe product; First group of article A meets the requirement of explosive on material character and density characteristics, in the explosive database and fall into high density explosive thing district, be confirmed to be explosive.
What Figure 11 provided is final alarming result, in the reality, and the article of not reporting to the police, promptly second group of article B is the pure water that one bottle of 380ml adorns, and the warning article, promptly first group of article A is a column type simulation explosive, alarming result meets the actual conditions of inspected object.

Claims (9)

1. method of utilizing multi-view X ray that explosive in the luggage is surveyed automatically, this method is mainly surveyed the explosive in the luggage automatically by material and density two dimensional features, it is characterized in that, may further comprise the steps:
1) picture signal at each visual angle of the tested luggage and articles of collection;
2) picture signal of obtaining is carried out material automatic detection: at first picture signal is carried out image segmentation, and the high energy after cutting apart and low energy image are carried out material characteristics identification according to double energy X-ray material discrimination principle, the hazardous location that meets the explosive material characteristics is marked the position that wherein is fit to participate in density calculation;
3) above-mentioned hazardous location segmentation result is carried out density automatic detection: at first carry out the hazardous location coupling according to the matched rule of setting in advance, find out each hazardous location of the same object of representative in the multi-view image, rebuild by various visual angles then and calculate, obtain the density value of hazardous location;
4) according to material automatic detection result and hazardous location density detection result, whether decision reports to the police to these luggage and articles;
Wherein, picture signal described in the step 1) is the double energy X-ray image that comprises high energy image signal and low energy image signal.
2. according to the described method of utilizing multi-view X ray that explosive in the luggage is surveyed automatically of claim 1, it is characterized in that, in step 2) in, described image segmentation specifically comprises:
A) utilize edge detection operator that the double energy X-ray image is carried out Edge Gradient Feature
B) in the neighborhood of extraction edge feature, search background pixel and the foreground pixel that comprises background, and adopt background to reject formula and ask the foreground image that obtains through after the background rejecting, wherein,
Described background is rejected G T = G T + B &times; G A G B , G BBe background pixel, G T+BFor comprising the foreground pixel of background, G ABe the gray-scale value of double energy X-ray image under the air, and background pixel G BGray-scale value greater than the foreground pixel G that comprises background T+BGray-scale value, G TBe the foreground image after rejecting through background.
3. according to the described method of utilizing multi-view X ray that explosive in the luggage is surveyed automatically of claim 1, it is characterized in that, also comprise in the step 3) corresponding explosive database is set up at each visual angle.
4. according to the described method of utilizing multi-view X ray that explosive in the luggage is surveyed automatically of claim 1, it is characterized in that step 4) further comprises, a plurality of visual angles density detection value is asked average.
5. according to the described method of utilizing multi-view X ray that explosive in the luggage is surveyed automatically of claim 1, it is characterized in that, calculate at the algebraic reconstruction of various visual angles described in the step 3) and finish according to the following steps:
A) set up various visual angles and rebuild coordinate system, soon the transfer passage cross section is divided into some little squares with the width of a pixel, and is each cell initialize 0, and sets up coordinate with the center of circle of each visual angle x-ray source again;
B) determine various visual angles algebraic reconstruction model;
C) determine the zone that the various visual angles algebraic reconstruction calculates, obtain the polygonal region that a plurality of visual angles surround the object cross section according to the projected image at each visual angle;
D) according to the zone of the various visual angles algebraic reconstruction model of determining with the calculating of various visual angles algebraic reconstruction, and by various visual angles algebraic reconstruction computing formula density value is calculated, described computing formula is as follows:
f ( k + 1 ) ij = f ( k ) ij + &lambda; &times; p - &Sigma; n = 1 v ( w ijn &times; f ( k ) ijn ) &Sigma; n = 1 v w ijn &times; w ijn &times; w ijn k < t k f ( k ) ij = 0 , 1 N &Sigma; f ^ ij > &beta; k f ( k ) i , j , k &GreaterEqual; t k
Wherein, k is an iterations, f (k) IjBe f after the k time iteration IjThe value of cell, f (k+1) IjBe f after the k+1 time iteration IjThe value of cell, λ is a coefficient of relaxation, P is the projection value of ray L correspondence, w IjnExpression ray L passes cell f IjLength, v is that ray L passes the sum that is included in the cell in the polygon ABCDEF, w IjBe that ray L passes cell f IjLength, t kFor carrying out the threshold value of quickening interative computation, α kBe the upper limit segmentation threshold of the k time iteration, β kBe the lower limit segmentation threshold of the k time iteration, N is for satisfying f (k) jThe pixel number of＞β (k) condition.
6. according to the described method of utilizing multi-view X ray that explosive in the luggage is surveyed automatically of claim 3, it is characterized in that, described explosive database is that horizontal ordinate is that checking matter atomic number, the longitudinal axis are the differentiation plane of the two dimension " effective atomic number-density " that density measurement surrounded, and comprises high density explosive thing and two districts of low-density explosive in this explosive database.
7. according to claim 1 or the 6 described methods of utilizing multi-view X ray that explosive in the luggage is surveyed automatically, it is characterized in that, step 4) further comprises, if this mean value falls into explosive database middle-high density explosive district and low-density explosive district, then the article with correspondence are judged to explosive, and are reported to the police.
8. according to claim 1 or the 6 described methods of utilizing multi-view X ray that explosive in the luggage is surveyed automatically, it is characterized in that step 4) further comprises, if this mean value falls into other zones of explosive database, then the corresponding article of institute are judged to safe product, and will not report to the police.
9. a device that utilizes multi-view X ray that explosive in the luggage is surveyed automatically is characterized in that, comprise at the bottom of overall treatment computing machine, conveyor, transfer passage, top, left side are according to visual angle module, right side according at the bottom of visual angle module and the left side according to the visual angle module; Wherein, described overall treatment computing machine for be connected to top, described left side according to visual angle module, right side at the bottom of according at the bottom of visual angle module and the left side according to the visual angle module, and top, described left side according to visual angle module, right side at the bottom of according at the bottom of visual angle module and the left side according to the visual angle module for being installed in the described transfer passage, and described conveyor is for being controlled by described overall treatment computing machine; Top, described left side according to the visual angle module comprise the top, left side according to x-ray source and to should the top, left side according to L type first detector of x-ray source; At the bottom of the described right side according to the visual angle module comprise at the bottom of the right side according to x-ray source and at the bottom of should the right side according to L type second detector of x-ray source; At the bottom of the described left side according to the low energy X ray source, with at the bottom of this left side according to photograph sigmatron source at the bottom of corresponding L type the 3rd detector in low energy X ray source, the left side and with at the bottom of the described left side according to corresponding L type the 4th detector in sigmatron source.
CN2009100884957A 2009-07-03 2009-07-03 Method and device for automatically detecting explosive substances in luggage by using multi-view X rays Active CN101592579B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100884957A CN101592579B (en) 2009-07-03 2009-07-03 Method and device for automatically detecting explosive substances in luggage by using multi-view X rays

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100884957A CN101592579B (en) 2009-07-03 2009-07-03 Method and device for automatically detecting explosive substances in luggage by using multi-view X rays

Publications (2)

Publication Number Publication Date
CN101592579A true CN101592579A (en) 2009-12-02
CN101592579B CN101592579B (en) 2011-07-06

Family

ID=41407349

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100884957A Active CN101592579B (en) 2009-07-03 2009-07-03 Method and device for automatically detecting explosive substances in luggage by using multi-view X rays

Country Status (1)

Country Link
CN (1) CN101592579B (en)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102455305A (en) * 2010-11-01 2012-05-16 北京中盾安民分析技术有限公司 Device for detecting hidden matters on human body by using dual-energy X-rays and image processing method used by device
CN102901740A (en) * 2012-07-30 2013-01-30 公安部第一研究所 Channel type four visual angles X ray liquid state article safety inspection system
CN102928448A (en) * 2012-10-30 2013-02-13 公安部第一研究所 Channel type four-perspective X-ray liquid goods safety inspection method and inspection device
CN103630563A (en) * 2013-12-04 2014-03-12 公安部第一研究所 X-ray helical scanning and navigating method
CN103903297A (en) * 2012-12-27 2014-07-02 同方威视技术股份有限公司 Three-dimensional data processing and identification method
CN105068134A (en) * 2015-07-29 2015-11-18 公安部第一研究所 Method for detecting dangerous goods hidden in shoe by using X-ray multi-view image
WO2016095798A1 (en) * 2014-12-18 2016-06-23 同方威视技术股份有限公司 Method for positioning target in three-dimensional ct image and security check system
CN107202808A (en) * 2017-06-14 2017-09-26 上海英曼尼安全装备有限公司 A kind of double source double-visual angle safety inspection method and its system
WO2018121325A1 (en) * 2016-12-29 2018-07-05 同方威视技术股份有限公司 Multi-angle imaging data processing method and apparatus
CN108345040A (en) * 2018-03-02 2018-07-31 上海瑞示电子科技有限公司 Dual-beam source detection apparatus and its control device go out Shu Fangfa
CN109410226A (en) * 2018-09-20 2019-03-01 公安部第研究所 A kind of package image cutting method and corresponding rays safety detection apparatus for safety check
CN110133739A (en) * 2019-04-04 2019-08-16 南京全设智能科技有限公司 A kind of X-ray rays safety detection apparatus and its drawing method is known automatically
CN111290040A (en) * 2020-03-12 2020-06-16 安徽启新明智科技有限公司 Active double-view-angle correlation method based on image recognition
CN112542007A (en) * 2020-11-30 2021-03-23 福州外语外贸学院 Dangerous target detection method and system between financial withdrawals
WO2024043190A1 (en) * 2022-08-22 2024-02-29 株式会社 東芝 Inspection device, inspection system, and inspection method

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102455305A (en) * 2010-11-01 2012-05-16 北京中盾安民分析技术有限公司 Device for detecting hidden matters on human body by using dual-energy X-rays and image processing method used by device
CN102901740A (en) * 2012-07-30 2013-01-30 公安部第一研究所 Channel type four visual angles X ray liquid state article safety inspection system
CN102901740B (en) * 2012-07-30 2014-10-22 公安部第一研究所 Channel type four visual angles X ray liquid state article safety inspection system
CN102928448A (en) * 2012-10-30 2013-02-13 公安部第一研究所 Channel type four-perspective X-ray liquid goods safety inspection method and inspection device
CN103903297A (en) * 2012-12-27 2014-07-02 同方威视技术股份有限公司 Three-dimensional data processing and identification method
CN103903297B (en) * 2012-12-27 2016-12-28 同方威视技术股份有限公司 Three-dimensional data processes and recognition methods
CN103630563A (en) * 2013-12-04 2014-03-12 公安部第一研究所 X-ray helical scanning and navigating method
US10145977B2 (en) 2014-12-18 2018-12-04 Nuctech Company Limited Method for positioning target in three-dimensional CT image and security check system
WO2016095798A1 (en) * 2014-12-18 2016-06-23 同方威视技术股份有限公司 Method for positioning target in three-dimensional ct image and security check system
CN105068134A (en) * 2015-07-29 2015-11-18 公安部第一研究所 Method for detecting dangerous goods hidden in shoe by using X-ray multi-view image
CN105068134B (en) * 2015-07-29 2017-09-15 公安部第一研究所 A kind of method that dangerous material are concealed in utilization X-ray multi-view image detection footwear
WO2018121325A1 (en) * 2016-12-29 2018-07-05 同方威视技术股份有限公司 Multi-angle imaging data processing method and apparatus
CN108267465A (en) * 2016-12-29 2018-07-10 同方威视技术股份有限公司 Various visual angles imaging data processing method and equipment
CN107202808A (en) * 2017-06-14 2017-09-26 上海英曼尼安全装备有限公司 A kind of double source double-visual angle safety inspection method and its system
CN108345040A (en) * 2018-03-02 2018-07-31 上海瑞示电子科技有限公司 Dual-beam source detection apparatus and its control device go out Shu Fangfa
CN109410226A (en) * 2018-09-20 2019-03-01 公安部第研究所 A kind of package image cutting method and corresponding rays safety detection apparatus for safety check
CN109410226B (en) * 2018-09-20 2021-03-30 公安部第一研究所 Package image segmentation method for security inspection and corresponding security inspection equipment
CN110133739A (en) * 2019-04-04 2019-08-16 南京全设智能科技有限公司 A kind of X-ray rays safety detection apparatus and its drawing method is known automatically
CN111290040A (en) * 2020-03-12 2020-06-16 安徽启新明智科技有限公司 Active double-view-angle correlation method based on image recognition
CN112542007A (en) * 2020-11-30 2021-03-23 福州外语外贸学院 Dangerous target detection method and system between financial withdrawals
WO2024043190A1 (en) * 2022-08-22 2024-02-29 株式会社 東芝 Inspection device, inspection system, and inspection method

Also Published As

Publication number Publication date
CN101592579B (en) 2011-07-06

Similar Documents

Publication Publication Date Title
CN101592579B (en) Method and device for automatically detecting explosive substances in luggage by using multi-view X rays
CN201514387U (en) System for utilizing multi-view X ray to automatically detect explosive in luggage
Liu et al. Automatic recognition of pavement cracks from combined GPR B-scan and C-scan images using multiscale feature fusion deep neural networks
US10042079B2 (en) Image-based object detection and feature extraction from a reconstructed charged particle image of a volume of interest
CN102928448B (en) Channel-type four visual angle X ray liquid article safety detection method and testing fixture
CN101231254B (en) Double source three-dimensional image forming method and system
CN102460067B (en) Systems and methods for automated, rapid detection of high-atomic-number materials
CN103675931B (en) CT system and the detection device for CT system
CN101403710A (en) Liquid article examination method and equipment
CN102920459B (en) Human body circumference parameter measuring method based on three-dimensional point cloud
WO2006053279A2 (en) Non-intrusive container inspection system using forward-scattered radiation
CN104165896A (en) Liquid goods safety inspection method and device
CN101592622A (en) Multi-view X ray luggage explosive substance automatic detection device with real dual-energy
CN202948145U (en) CT system and detecting device applied to the CT system
Zhao et al. Testing gravity theories via transverse Doppler and gravitational redshifts in galaxy clusters
CN102901740B (en) Channel type four visual angles X ray liquid state article safety inspection system
CN101900694B (en) Dual-energy under-sampling substance recognition system and method based on straight path scanning
CN104181178A (en) Channel type two-viewing-angle X-ray liquid substance safety check system
CN106164707A (en) It is used in particular for method and the X ray detection system of the Non-Destructive Testing of object
CN104567758A (en) Three-dimensional imaging system and method thereof
CN103903298B (en) Three-dimensional data is processed and recognition methodss
JP2021516760A (en) Methods and Devices for Detecting and / or Identifying Materials and Articles Using Charged Particles
CN102162798B (en) Method and equipment for inspecting liquid article
CN105784741A (en) Multi-visual-angle X ray safety check device
CN202903699U (en) Channel type four-visual angle X-ray liquid article safety check system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant