CN106097407A - Image processing method and image processing apparatus - Google Patents

Image processing method and image processing apparatus Download PDF

Info

Publication number
CN106097407A
CN106097407A CN201610371890.6A CN201610371890A CN106097407A CN 106097407 A CN106097407 A CN 106097407A CN 201610371890 A CN201610371890 A CN 201610371890A CN 106097407 A CN106097407 A CN 106097407A
Authority
CN
China
Prior art keywords
determinand
grey scale
scale curve
under test
projection
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.)
Pending
Application number
CN201610371890.6A
Other languages
Chinese (zh)
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.)
Tsinghua University
Original Assignee
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
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN201610371890.6A priority Critical patent/CN106097407A/en
Priority to KR1020160121190A priority patent/KR101847730B1/en
Publication of CN106097407A publication Critical patent/CN106097407A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • 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
    • G01V5/22Active interrogation, i.e. by irradiating objects or goods using external radiation sources, e.g. using gamma rays or cosmic rays
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/32Determination of transform parameters for the alignment of images, i.e. image registration using correlation-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • G06T7/337Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods involving reference images or patches
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • G06T7/77Determining position or orientation of objects or cameras using statistical methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Probability & Statistics with Applications (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Measurement Of Radiation (AREA)

Abstract

The present invention relates to a kind of image processing method and image processing apparatus.The method can include that the gradation data arrived according to diverse location sensor acquisition after projecting determinand obtains the grey scale curve step of projection picture;The step of the sensing station that flex point is corresponding in grey scale curve is obtained according to this grey scale curve;The position of the sensing station corresponding according to described flex point and radiation source obtains may combining of the vertex position of determinand;And in the possible combination of the vertex position of described object under test, the step of the vertex position transitivity of determinand is determined according to this grey scale curve.The invention also discloses a kind of image processing apparatus.

Description

Image processing method and image processing apparatus
Technical field
The invention belongs to technical field of image processing, relate more specifically to be mainly used in the image processing method of safety inspection And image processing apparatus.
Background technology
At present, the criminal behavior such as smuggle, steal into another country and frequently occur, threatened the collection of customs duty of national customs, to national economy Bring massive losses.Behavior also serious harm state-owned enterprise of China and national industry are stolen into another country in the smuggling in the field such as chemical industry, food Development, very disruptive China market order.It addition, part lawless person illegally carries drugs to seek exorbitant profit, to the people Health also result in greatly harm.
In recent years, the malice implemented by carrying the dangerous goods such as gun, ammunition, cutter is done violence, the attack of terrorism etc. is serious Harm social public security event also happens occasionally.Therefore national interests people's safety and property safety are had by Security Inspection Equipments There is important meaning.
Existing Security Inspection Equipments mainly includes that Handheld vehicle bottom checks equipment and X-ray rays safety detection apparatus.Handheld vehicle The end, checks that equipment is by staff's handheld device, and vehicle chassis carries out optical scanning, and whether observation vehicle chassis has changes The vestige made or whether carry illegal article etc. on vehicle chassis secretly.But, offender pretends meticulously frequently with all kinds of methods Improved vehicle chassis, the article kept illegally after only cannot detecting such transformation with optical scanning and observation;X-ray safety check sets Standby, use radiographic viewing techniques, by being scanned article, and image output judges whether to carry secretly violated for inspection personnel Product.
Summary of the invention
Some safety detection devices and the processing method of imaging also cannot be obtained the structure within object under test and image, Cannot detect whether vehicle interior carries, mails illegal article quickly and efficiently.For the detection of device internal structure can be detected Mode, currently mainly coordinates corresponding image processing method for radiant image class detection device.But, such device and method are also Cannot be by data messages such as the physical property of the more specifically determinand such as detection the acquisition position of object under test, size, density.And And some determinands being also possible that, some active materials, general detection equipment are difficult to detect that well shielded puts Penetrating property material.
This application provides a kind of image at least two non-parallel direction, determinand projected in least one set Processing method, the method includes that the gradation data arrived according to diverse location sensor acquisition after projecting determinand obtains projection The grey scale curve step of picture;The step of the sensing station that flex point is corresponding in grey scale curve is obtained according to this grey scale curve;According to Sensing station and the possible of vertex position of the position acquisition determinand of radiation source that described flex point is corresponding combine;And root In the possible combination of the vertex position of described object under test, the vertex position transitivity of determinand is determined according to this grey scale curve Step.
On the other hand, this application provides one determinand to be carried out in the non-parallel direction of at least two in least one set The image processing apparatus of projection, this device includes the gray scale for arriving according to diverse location sensor acquisition after projecting determinand Data obtain the module of the grey scale curve of projection picture;For obtaining, according to this grey scale curve, the biography that in grey scale curve, flex point is corresponding The module of sensor position;The top of determinand is obtained for the position of the sensing station corresponding according to described flex point and radiation source The module that may combine of some position;And for according to this grey scale curve at possible group of vertex position of described object under test Conjunction determines the module of the vertex position transitivity of determinand.
According to disclosed embodiments, by analyze determinand gained grey scale curve after transmission projection feature and Its relation with object under test, finds out the corner position acquired a special sense on curve, and combines the structure of optical projection system, will be many Many groups grey scale curve of Angles Projections gained is interrelated to be calculated, thus analyze the size of object under test, position and Some physical property characteristics such as density, the gamma ray absorption coefficient etc. of object.The method can be answered in the case of radioactive source selects suitably Use all kinds of safe, medical etc. check.Owing to active material stores frequently with the container that density is bigger, in this kind of inspection In survey, procurator can check high-density objects according to the physical property characteristic emphasis of the determinand that the method is measured, thus realizes Detection to active material.
Accompanying drawing explanation
By being more clearly understood from the features and advantages of the present invention with reference to accompanying drawing, accompanying drawing is schematic and should not manage Solve as the present invention is carried out any restriction, in the accompanying drawings:
Fig. 1 illustrates that determinand is thrown by one according to an embodiment of the invention in the non-parallel direction of at least two The schematic diagram of the system of shadow;
Fig. 2 shows in the system shown in Fig. 1 according to an embodiment of the invention by a radioactive source to be measured The grey scale curve figure that thing is obtained by detector in the case of projecting;
Fig. 3 illustrates the flow chart of a kind of image processing method according to an embodiment of the invention;
Fig. 4 illustrates the schematic diagram of determinand in a kind of image processing method according to an embodiment of the invention;
Fig. 5 illustrates the vertex position obtaining determinand in a kind of image processing method according to an embodiment of the invention May the schematic diagram of method of combination;
Fig. 6 illustrates the schematic diagram emulating a kind of image processing method according to an embodiment of the invention;
Fig. 7 a illustrates and according to an embodiment of the invention a kind of image processing method emulates obtain to be measured The schematic diagram of thing grey scale curve image;
Fig. 7 b illustrates and according to an embodiment of the invention a kind of image processing method emulates obtain to be measured The schematic diagram of the first derivative of thing grey scale curve image;
Fig. 7 c illustrates and according to an embodiment of the invention a kind of image processing method emulates obtain to be measured The schematic diagram of the second dervative of thing grey scale curve image;
Fig. 8 illustrates that according to an embodiment of the invention carry out a kind of image processing method emulates acquisition determinand The schematic diagram of the method that may combine of vertex position;
Fig. 9 illustrates object under test in emulating a kind of image processing method according to an embodiment of the invention The comparison diagram that may combine corresponding grey scale curve and determinand grey scale curve of vertex position;
Figure 10 illustrate a kind of image processing method according to an embodiment of the invention emulate in object under test The comparison diagram that may combine corresponding grey scale curve and determinand grey scale curve of vertex position;
Figure 11 illustrate a kind of image processing method according to an embodiment of the invention test in the reality of determinand Test design sketch;
Figure 12 illustrates the structural representation of a kind of image processing apparatus according to an embodiment of the invention.
Detailed description of the invention
Feature and the exemplary embodiment of various aspects of the invention are described more fully below.In following detailed description In, it is proposed that many details, in order to complete understanding of the present invention is provided.But, to those skilled in the art It will be apparent that the present invention can implement in the case of some details in need not these details.Below to enforcement The description of example is only used to be provided by the example illustrating the present invention be better understood from the present invention.The present invention never limits In any concrete configuration set forth below and method, but cover under the premise of without departing from the spirit of the present invention element, Parts and any amendment of algorithm, replace and improve.In the the accompanying drawings and the following description, it is shown without known structure and skill Art, in order to avoid the present invention causes unnecessary obscuring.
Below in conjunction with the accompanying drawings, image processing method according to embodiments of the present invention and image processing apparatus are described in detail.
Fig. 1 illustrates that determinand is thrown by one according to an embodiment of the invention in the non-parallel direction of at least two The schematic diagram of shadow system.As it is shown in figure 1, the system projecting determinand 3 can be included in the non-parallel direction of at least two and treat Surveying the radioactive source 1 that thing 3 carries out projecting, can use one or more radioactive source, can use two radioactive sources here, one can To be arranged on the bottom of determinand, another can be arranged on the side of determinand.Can also only use a radioactive source, extremely Determinand is projected and obtains two identical projection pictures of record INFORMATION OF INCOMPLETE by few two non-parallel directions, and radioactive source is Radiation source including radioactive substance.Wherein, radioactive substance such as can include following one or more: tritium, ferrum 55, cobalt 60, nickel 63, selenium 75, antimony 124, ytterbium 169, thulium 170, Ir192, thallium 204, polonium 210, plutonium 238, Krypton 85, strontium 90, Ce 137, promethium 147 With the radium in some actinides such as plutonium 239, americium 241, californium 252 etc., sodium 22, cobalt 57, yttrium 88, cadmium 109, bismuth 207, uranium-radium series 226 etc..In a kind of exemplary embodiment, radioactive source uses Co 60 as radioactive substance.The system of measuring targets projection Can also include that be correspondingly arranged with radioactive source receives the sensor 2 launching particle, sensor 2 such as can include observing, remembering The detector of record number of particles, such as ionization chamber, proportional counter, Geiger Müller counter, scintillation counter, quasiconductor are visited Survey device etc..This system also includes the image processing apparatus 4 processing projection picture.
Fig. 2 shows in the case of being projected determinand 3 by a radioactive source in the system shown in figure 1, passes through The grey scale curve figure that detector obtains.From Figure 2 it can be seen that grey scale curve be the ray with projection by path object density and Length and the curve that changes, its rule is followed as shown in formula (1):
I=I0·e-μt (1)
Wherein, μ represents the object absorptance to ray, and t represents that ray passes the path of object.When ray is not worn When crossing object under test, gray value is maximum I0;And when ray is through object under test, grey scale curve can not pass through less than ray Maximum I during object under test0.It addition, when ray is through polygonal a certain summit, grey scale curve can correspondingly go out The flex point of existing slope sudden change;And during ray other positions inswept, grey scale curve is a smooth curve, do not have The situation of slope sudden change.If polygon is n limit shape, then on curve, the relation of number p of flex point and n is 0 < p≤n.Therefore, ash Write music the corner position of line, the more specific location information of polygon vertex can be reflected;And by the gradation data of grey scale curve, The much informations such as the density of object under test, shape can be reflected again.
Fig. 3 illustrates the flow chart of a kind of image processing method according to an embodiment of the invention, and the method such as may be used Perform with the image processing apparatus 4 in system as shown in Figure 1.The beam formed due to radiographic source is sector, penetrates in judgement Line source, whether when determinand is projected by non-parallel direction, can compare whether the centrage in radiogenic fan-ray beam puts down OK.After determinand is projected by the non-parallel direction of at least two, can obtain at least two width projection pictures, this projection picture can Thinking gray-scale map, its gray scale represents the ray that in the system to determinand projection, detector is corresponding and passes respective substance pair on path Adding up of ray attenuation effect, can be obtained one according to gray-scale map and be changed by density and the length of object on path with ray Grey scale curve.As it is shown on figure 3, in step S31, can arrive according to diverse location sensor acquisition after determinand is projected Gradation data obtains the grey scale curve of projection picture.In step s 32, can turn according in this grey scale curve acquisition grey scale curve The sensing station that point is corresponding.In step S33, the summit of determinand can be obtained according to the sensing station that this flex point is corresponding May combining of position.In step S34, can be according to this grey scale curve at possible group of vertex position of described object under test Conjunction determines the vertex position transitivity of determinand.Here physical property can include that such as density, gamma ray absorption coefficient etc. are to be measured The physical property of thing.
In one example, can be filtered processing to the grey scale curve that determinand projection is obtained, after making filtering, bent Line does not have the higher noise of frequency and fluctuation, the smooth disappearance causing flex point that again will not be excessive or write music with determinand ash Line compare have bigger delayed.Such as to after filtering high frequency noise, grey scale signal is being carried out compensation of delay.In one example, The grey scale curve that projection obtains can be carried out 2 derivations, in second dervative image, corner position can be with impulse form body Existing, find this few group pulse positions, the position of flex point in grey scale curve can be drawn.
Fig. 4 illustrates the schematic diagram of determinand in a kind of image processing method according to an embodiment of the invention.One A little safety checks, import and export in the application such as detection, such as, when carrying out the inspection of core security, in chest often occurs, have less miscellaneous of density Thing and the bigger radioactive source of density or the situation of shielding material.When radioactive source is shielded by high density materials such as using lead Time, passive type radiological measuring often cannot detect that object under test has radioactivity.Can be inspection in the inspection of core security Look into the high-density objects by all kinds of low-density foreign material cinctures focused on as shown in Figure 4.Owing to ray is through anti-after object Reflect is the integration of material on path, and high-density objects is relative to low-density object, and obtain after passing object for ray penetrates Line strength often has large effect, therefore, according to the method shown in Fig. 3, can be write music by the ash after projecting determinand Corner position on line obtains position and the size shape of high-density objects interested, and by its half-tone information, obtains Average data in determinand density.Although Fig. 4 is illustrated that the convex polygon that projection picture is even density of determinand Situation, but the method can also be applied in the middle of the situation that projection picture is concave polygon in determinand cross section, can regard as by height Concave polygon and the polygon being made up of air of density form a new convex polygon, and the concrete data calculated are the two Average effect, although the result possible error obtained in this case through this image processing method can be slightly larger than determinand The situation that projection picture is convex polygon, but remain able to preferably to reduce the physical property of object under test.The method can also Apply in the determinand image of uneven density, determinand can be regarded as the combination of multiple uniform density determinand, pass through The method obtains the meansigma methods of density that determinand density can regard the combination of multiple uniform density determinand as, but this density Meansigma methods remain able to embody the density variation of determinand and surrounding, i.e. in background environment, preferably reduction is to be measured The physical property of object.
Fig. 5 illustrates the vertex position obtaining determinand in a kind of image processing method according to an embodiment of the invention May the schematic diagram of method of combination.As shown in fig. 5, it is assumed that the image obtaining determinand projection is convex four shown in shade Limit shape, can set up two-dimensional coordinate system in plane, the system coordinates that will project determinand, according to the parameter of system structure, Several flex points can be reduced in system, it is possible to by location position in corresponding two-dimensional spatial location.Assume that spy Survey device array and be arranged on three directions of object under test, then detector array can be distributed in the positive axis both sides of y-axis with coordinatograph, puts The source b of penetrating can be arranged on the negative semiaxis of y-axis, and radiation source a can be arranged on the detector array opposite side relative to determinand, The grey scale curve that determinand projection is obtained is found out the position of the flex point of slope sudden change, by its correspondence to corresponding coordinate system Detector array on.Due to the polygonal summit in all corresponding object under test cross section of the flex point in grey scale curve, then radioactive source And the line between flex point is bound to by polygonal summit to be measured.Can according to the position of radiation source and with radiation source pair What the line of the detector position that the flex point answered is corresponding obtained the vertex position of object under test may combination.In one example, Owing to the information relevance in the projection picture that obtains when determinand is carried out rectangular projection is minimum, the valid data obtained are Many, therefore can choose orthogonal projection picture and carry out image procossing.Here assume that and project on a direction wherein To grey scale curve in have p flex point, project in the other directions and the grey scale curve obtained have q flex point, and p≤q.? When plane sets up two-dimensional coordinate system, the detector position coordinate that during twice projection, the point coordinates at radioactive source place is corresponding with flex point Line can get p q different intersection point, then the group of n point during polygonal summit, n limit one is set to q intersection point of this p Close.May be then the group number K on polygonal summit to be measured1For:
K 1 = Σ n = q n = p q C p q n - - - ( 2 )
Be can be seen that by formula (2), it is assumed that work as p=3, q=4, K1=3797, amount of calculation is the biggest.Further, along with p's Yu q Increasing, amount of calculation will increase by a fairly big margin.The straight-line pass radioactive source being linked to be due to object under test any two summit be little generally Rate event, therefore can regard as and have on each ray and the summit of an only object under test.It may be then polygon to be measured The group number K on summit2For:
K2=n! (3)
Now be can be seen that by formula (3), it is assumed that as n=4, K2=24.Amount of calculation is by possible group of determinand summit Close number K1=3797 reduce to 24, and this will increase substantially calculating speed, decrease the complexity that the method calculates.Can be by May being connected two-by-two on the possible summit in combination of these determinand summits, take the outline of all lines as border, by this Image matrixing, it is also possible to converted with physical length according to pixel size in image outline inside and fill, such as, needing Survey object is sentenced numeral 1 and is filled with, and sentences numeral 0 without object and is filled with.Just available institute likely with object under test one The position of the object caused and the image array of expression determinand shape.They are reduced in system from two-dimensional coordinate.Also It is exactly can be with the possible imaginary determinand constituted that combines of the vertex position of measuring targets at radiation source positions and detector Carrying out imaginary determinand projection in the case of invariant position, available a series of new only reflection rays pass the length of object Grey scale curve.Due to grey scale curve here do not comprise determinand with ray by the letter of the variable density of object on path Breath, so grey scale curve here is different from the grey scale curve that the projection of former determinand is obtained by radioactive source.
In one embodiment, the gradation data that can collect according to diverse location detector after imagination determinand projection Obtain the imagination grey scale curve of determinand projection picture and compare with the grey scale curve of object under test projection picture, can choose with The imaginary determinand vertex position that the immediate grey scale curve of grey scale curve of object under test projection picture is corresponding is object under test Vertex position.In one example, the data can chosen on detector between minimum flex point and maximum flex point calculate, The datum mark of determinand grey scale curve with imagination determinand grey scale curve is alignd, then compares the difference of two curves.If to imagination A series of grey scale curve reflection ray that determinand projection obtains is S through the value of object lengthm,x, determinand is projected by system Obtained former gradation data is Im,x, not placing the gradation data obtained by object system is I0If an alignment coefficient is μm,x, M is the m article curve in a series of arc curve, and x is x-th point in a series of grey scale curve.Then there is a following relational expression:
Im,x=I0·exp(-μm,xSm,x) (4)
μ m , x = l n ( I 0 / I m , x ) S m , x - - - ( 5 )
Calculate each μ corresponding to point of two curvesm,xAnd take its average, can obtain should the μ of curvemValue, it After do following computing:
Pm,x=I0·exp(-μmSm,x) (6)
An available new curve, seeks Pm,xWith Im,xArticle two, deviation ν of curvemFor:
v m = Σ ( I m , x - P m , x ) 2 n - - - ( 7 )
Wherein, n is Pm,xThe number of upper corresponding point.Try to achieve the ν of minimummImaginary determinand polygon corresponding to value and Coordinate, is shape and the position of required object under test.
In one embodiment, can be according to detector gained original gradation data when not placing determinand, to determinand Gradation data that after projection, diverse location detector collects and radiation source ray after imagination determinand projection is treated through imagination The length value surveying thing obtains the determinand absorptance to ray, and then obtains the physical property of determinand.In one example, obtaining To with original gradation curve closest to imaginary determinand grey scale curve after, with the μ in its corresponding formula (5) of this grey scale curve Value is the object under test absorptance to system ray.
In one embodiment, in the case of the ray particle energy of the radioactive source selected is metastable, the method can With according to object under test to the absorptance of system ray and the density of the ray particle energy balane determinand of radioactive source, real Differentiation referring now to object under test physical property.Such as, when the radioactive source of system uses Co 60 as radioactive substance, its ray grain Sub-energy is 1.25MeV, primarily serves the purpose of Compton scattering between particle and object under test.Under this energy, particle passes μ/ρ value of object can approximate regards constant value as, thus obtains the density of object under test accordingly.
Fig. 6 illustrates the schematic diagram emulating a kind of image processing method according to an embodiment of the invention.As Shown in Fig. 6, in one embodiment, two radioactive sources can be chosen, radioactive source can be individually positioned in system bottom and side Lower section, can circumferentially arrange detector to guarantee the total data that can obtain ray through object under test in three faces.To system Parameter be set, by system parameter settings be such as: wherein, a is system for a=350cm, d=350cm, d=50cm Width, b is the height of system, and d is the radioactive source 1 distance away from ground.It is assumed here that object under test is an arbitrary convex pentagon, to be measured Thing place height is h, such as, assume that determinand, on large truck, so, places object ground height h onboard For example, it can be set to be 150cm, then the theoretical value on each summit of determinand is given in Table 1, the most each coordinate unit is centimetre.
Table 1
Assume that this object under test can be set as μ=0.9, visiting time unloaded for the absorptance of ray in simulations Survey the initial gray data obtained by device and can be set as I0=10;A detector is arranged every 1cm on detector array.Figure 7a illustrates the determinand grey scale curve emulating a kind of image processing method and obtaining according to an embodiment of the invention The schematic diagram of image.As shown in Figure 7a, send ray measuring targets by radioactive source and carry out the projection in both direction, can obtain The grey scale curve of the projection picture in both direction.Fig. 7 b illustrate according to an embodiment of the invention to a kind of image at Reason method carries out emulating the schematic diagram of the first derivative of the determinand grey scale curve image obtained;Fig. 7 c illustrates according to the present invention's Showing of the second dervative that a kind of image processing method is emulated the determinand grey scale curve image obtained of one embodiment It is intended to.Two grey scale curve can be done as shown in Fig. 7 b and 7c derivation process.In the curve of second dervative, in grey scale curve Corner position have the projection of corresponding pulsed, the corresponding coordinate x in detectordIt is given in Table 2 respectively.
Table 2
These points are placed under two-dimensional coordinate system, according to system structure, transform it into coordinate in projection systems and exist Table 3 is given.
Table 3
Can obtain according to optical projection system, the coordinate of two groups of radioactive sources is given in Table 4.
Table 4
Fig. 8 illustrates that according to an embodiment of the invention carry out a kind of image processing method emulates acquisition determinand The schematic diagram of the method that may combine of vertex position.As shown in Figure 8, can be by flex point coordinate corresponding for radioactive source coordinate Being connected, radiographic source 1 can obtain 10 line segments and 25 different intersection points with radiographic source 2 after being connected from corresponding flex point.At these Intersection point is selected corresponding combination, and the polygon that all possible combination draws is projected, and choose specified point and make It is corresponding with former grey scale curve data that the polygon that all possible combination draws carries out projecting gained grey scale curve data.Fig. 9 Illustrate according to an embodiment of the invention a kind of image processing method emulated in the vertex position of object under test The comparison diagram of corresponding grey scale curve and determinand grey scale curve may be combined;Figure 10 illustrates according to one embodiment of present invention A kind of image processing method emulate in object under test vertex position may grey scale curve corresponding to combination with to be measured The comparison diagram of thing grey scale curve.As shown in Figure 9 and Figure 10, the position of correspondence detector after abscissa represents determinand projection in figure Unit is centimetre to represent, vertical coordinate represents that ray passes determinand length value and the determinand product to gamma ray absorption coefficient.Will Comparing imagination determinand projection the data obtained with determinand projection the data obtained, the combination obtaining deviation value minimum is The coordinate on required determinand summit, the determinand apex coordinate finally given and error amount are given in Table 5.
Table 5
Simulation result according to this image processing method is to the maximum for the determination of polygon vertex position, absolute error 0.28cm, the relative error that attenuation quotient calculates is 0.33%, and every error is the least.
In testing at one, aluminium block can be projected the projection picture obtained and use at this image processing method Reason, verifies the effect of this image processing method.Determinand object be an aluminium material, solid, cross section be foursquare length Cube.Can choose a cross section of determinand, the position at pair cross-section place carries out dimensional measurement, measures gained cross section to be measured Dimensioning modest ability and width are 175mm.Measuring targets carries out twice projection measurement.Convenient for processing, can make at radioactive source At 180 ° with 240 °.In the measurements, same cross section measuring targets can be chosen and take multiple measurements, such as, carry out 4000 times Measure, take and measure the average of numerical value to improve the quality of data.The grey scale curve image that can obtain both direction projection enters Row secondary derivation obtains corner position, obtains the possible combination on determinand summit and to the combination of these summits by corner position The determinand become projects again, the grey scale curve obtained is compared with former grey scale curve, takes closest with virgin curve The summit that summit is determinand corresponding to grey scale curve.Show that this object under test each point coordinate figure provides in table 6, try to achieve and treat Surveying object is 0.0146 to the absorptance of ray.
Table 6
As can be seen from Table 6, this result substantially conforms to the shape of object under test, and its four length of side is respectively 172mm, 178mm, 173mm, 168mm, differ 7mm with its standard size 175mm maximum, meets the range of error analyzed when emulation.
It is pentagonal cylinder that aluminium block in a upper experiment is replaced by solid, cross section, the effect to this image processing method Fruit is verified again.Choosing a cross section of determinand, the position at pair cross-section place carries out dimensional measurement, measures gained and treats The length on five limits surveying the physical dimension in cross section is respectively limit 1=125cm, limit 2=70.71cm, limit 3=125cm, limit 4= 175cm, limit 5=175cm.Choose radioactive source to be at 270 ° and 210 °.In the measurements, same cross section can be chosen to determinand Body carries out 4000 times and measures, and takes and measures the average of numerical value to improve the quality of data.Both direction can be projected the ash obtained Line image of writing music carries out secondary derivation and obtains corner position, obtains may combining and to this of determinand summit by corner position The determinand that a little summits combine projects again, the grey scale curve obtained is compared with former grey scale curve, take and The summit that summit is determinand that the immediate grey scale curve of virgin curve is corresponding.Draw this object under test each point coordinate figure such as table 7 In be given, trying to achieve object under test to the absorptance of ray is 0.0143.
Table 7
As can be seen from Table 7, this result substantially conforms to the shape of object under test, and its 5 length of side is respectively as follows: limit 1= 124.45mm, limit 2=72.31mm, limit 3=129.84mm, limit 4=170.00mm, limit 5=172.80mm, with its standard size 175mm maximum difference 7mm, meets the range of error analyzed when emulation, has higher reduction degree.
Tabling look-up, for aluminum, when photon energy is 1MeV, absorptance can approximate than μ/ρ with determinand density Take constant value 6.146 × 10-3, when photon energy is 1.5MeV, what μ/ρ can approximate takes constant value 5.007 × 10-3, false If this partial value is linear distribution, then, when photon energy is 1.25MeV, what μ/ρ can approximate takes constant value 5.5765 × 10-3.The density then trying to achieve aluminum is ρAl'=2.618g/cm3, and the theoretical value of aluminium material density: ρAl=2.702g/cm3, contrast Measurement result, this image processing method preferably reduces the physical property of this object.
Bucket used by water dispenser can also project the projection picture obtained to use at this image processing method Reason, verifies the effect of this processing method.Bucket used by water dispenser can be filled water be positioned in container this bucket Dress water measures.Radioactive source can be made to be placed in 0 ° of position, and container enters system, and radioactive source and detector are scanned, To the radioactive source radiation image 0 ° of position;To 90 ° of positions, container can be made reversely to remove system, put by movable radiation source again Penetrate source and detector scanning, obtain the radiation image of 90 ° of positions.The a series of objects obtained through this image processing method continuous Cross section, these cross sections can rebuild the 3D model image of object under test.Figure 11 illustrates according to an embodiment of the invention one Plant the experiment effect figure of determinand during image processing method is tested.Every by the bucket that obtained by this image processing method The data in one cross section all carry out identical process, and available a series of length are unequal, high close rectangle.Use display soft Part measuring targets carries out reconstructing three-dimensional model, these rectangles is arranged in order in order, the graphics of available bucket Shape and three-view diagram are as shown in figure 11.
The gamma ray absorption coefficient that can be tried to achieve determinand bottled water by this image processing method is 7.2805 × 10-3.Table look-up Can obtain, for water, when photon energy is 1.25MeV, can approximate thinks μ/ρ=6.2 × 10-3.Then can try to achieve bottled water Density pw'=1.174g/cm3.And the actual value for the density of water isThen from experimental result, Its density measurement ρw' and actual valueComparing, its relative error is 17.4%.By experimental result it can be seen that test base This reflects the position at the shape of bucket, size and place exactly.But, the aluminium block experiment done before contrast, bucket is real For the experiment of the error tested relatively aluminium block, error is bigger.
Can analysis further to experimental result.For reducing error, projection picture can be carried out subtract shadow and process.Due to packaging Case is relatively-stationary background, therefore can the impact that image is caused by container be deducted accordingly.For convenience of processing, can be by image The middle region without object, as the background of container, takes respective column data and takes the mean as " the unloaded data " in experiment, this In " unloaded data " detector gained original gradation data when being equivalent to not place determinand.Bucket position is subtracted Shadow.Measurement result after being optimized digital image processing method is: when photon energy is 1.25MeV, penetrating of bottled water What the absorptance of line can approximate takes constant value 7.0794 × 10-3, then the density p of bottled water is tried to achievew'=1.142g/cm3, The relative error of density measure is down to 14.2% by 17.4% before.Preferably reduce the physical property of this object.
In one embodiment, a kind of image procossing dress determinand projected at least two non-parallel direction Put, can include for projection can be obtained according to the gradation data that diverse location sensor acquisition after projecting determinand arrives The module of the grey scale curve of picture.For the mould of the sensing station that flex point is corresponding in grey scale curve can be obtained according to grey scale curve Block.For can obtain the possibility of the vertex position of determinand according to the position of sensing station corresponding to flex point and radiation source The module of combination.For can determine the top of determinand in the possible combination of the vertex position of object under test according to grey scale curve The module of some position transitivity.In one example, it is also possible to include filtering high frequency for grey scale curve is filtered process The module of noise signal.In one example, it is also possible to include the module for grey scale curve being sought second dervative.Sensor example As the detector observing, recording number of particles can be included.
In one embodiment, this device could be included for can according to the position of radiation source and with described radiation What the line of the detector position that flex point corresponding to source is corresponding obtained the vertex position of object under test may the module of combination.
In one embodiment, this device could be included for combine with the possible of the vertex position of measuring targets The imaginary determinand constituted carries out imagination determinand projection in the case of described radiation source positions and detector position are constant Module.In one example, this device could be included for detect according to diverse location after imagination determinand projection Gradation data that device collects obtain imagination determinand projection picture grey scale curve and with the gray scale of object under test projection picture Curve ratio relatively, chooses the imaginary determinand summit corresponding with the immediate grey scale curve of grey scale curve of object under test projection picture Position is the module of the vertex position of object under test.
In one embodiment, could be included for can be former according to detector gained when not placing determinand for this device Beginning gradation data, determinand is projected after the gradation data that collects of diverse location detector and to imagination determinand projection after Radiation source ray obtains the determinand absorptance to ray through the length value of imagination determinand, and then obtains the thing of determinand The module of property.Figure 12 illustrates the structural representation of an embodiment of a kind of image processing apparatus according to the present invention.Process dress Putting and can use general-purpose computing system structure, computer system can specifically based on processor computer.Described process sets Standby entity includes input and output I/O interface 121, memorizer 122, at least one processor 123 and at least one communication interface 124.Wherein, input and output I/O interface 121, memorizer 122, at least one processor 123 and at least one communication interface 124 Between connected by communication bus 125.Described I/O interface 121, for receiving the text data from subscriber equipment, and by institute State transmission of textual data give described processor 123, wherein, described text data use SQL SQL form represent. Processor 123 can be a general central processor (CPU), microprocessor, ASIC (application- Specific integrated circuit, ASIC), or one or more for controlling the collection that the present invention program program performs Become circuit.Wherein, described communication bus 125 can include a path, transmits information between said modules.Described communication interface 124, use the device of any transceiver one class, be used for and other equipment or communication, such as Ethernet, wireless access network (RAN), WLAN (Wireless Local Area Networks, WLAN) etc..Computer system includes one or many Individual memorizer 122, can be read only memory (read-only memory, ROM) maybe can store static information and instruction its The static storage device of his type, random access memory (random access memory, RAM) or can store information and The other kinds of dynamic memory of instruction, it is also possible to be EEPROM (Electrically Erasable Programmable Read-Only Memory, EEPROM), read-only optical disc (Compact Disc Read- Only Memory, CD-ROM) or other optical disc storage, laser disc storage (include compress laser disc, laser dish, laser disc, digital universal Laser disc, Blu-ray Disc etc.), magnetic disk storage medium or other magnetic storage apparatus or can be used in carrying or storing that there is finger Order or the desired program code of data structure form by any other medium of computer access, but can be not limited to this. These memorizeies 122 are connected with processor 123 by communication bus 125.Wherein, described memorizer 122 is used for storing execution The application code of the present invention program, the application code performing the present invention program preserves in memory, and by processing Device 123 controls to perform.Described processor 123 is for performing the application program of storage in described memorizer 122.
It should be noted that, in the claims, word " comprises " or " including " is not precluded from existing and does not arranges in the claims Element or assembly.Article "a" or "an" before being positioned at element or assembly is also not precluded from there is multiple such element Or the situation of assembly.
Furthermore, it should also be noted that the language used in this specification primarily to the readable and purpose of teaching and select That select rather than select to explain or limit subject of the present invention.Therefore, without departing from appended claims In the case of scope and spirit, many modifications and changes are all apparent for those skilled in the art 's.About the scope of the present invention, the description done in description is all illustrative and not restrictive, the scope of the present invention It is defined by the appended claims.

Claims (14)

1. the image processing method at least two non-parallel direction, determinand projected, including:
The gradation data arrived according to diverse location sensor acquisition after projecting determinand obtains the grey scale curve of projection picture;
The sensing station that in grey scale curve, flex point is corresponding is obtained according to described grey scale curve;
The possible combination of the vertex position of determinand is obtained according to the sensing station that described flex point is corresponding;And
In the possible combination of the vertex position of described object under test, the vertex position of determinand is determined according to described grey scale curve Transitivity.
2. the method for claim 1, also includes being filtered described grey scale curve processing, filters high-frequency noise letter Number.
3. the method for claim 1, also includes described grey scale curve is asked second dervative.
The most the method for claim 1, wherein according to position and the institute corresponding with described radiation source of described radiation source State the possible combination of the vertex position of the line acquisition object under test of sensing station corresponding to flex point.
5. the method for claim 1, also includes what the possible combination of the vertex position to described object under test was constituted Imagination determinand carries out imagination determinand projection in the case of described radiation source positions and sensing station are constant.
6. method as claimed in claim 5, wherein, the ash arrived according to diverse location sensor acquisition after imagination determinand projection Degrees of data obtains the grey scale curve of imagination determinand projection picture and compares with the grey scale curve of object under test projection picture, chooses The imaginary determinand vertex position corresponding with the immediate grey scale curve of grey scale curve of object under test projection picture is determinand The vertex position of body.
7. method as claimed in claim 5, wherein, according to sensor gained original gradation data when not placing determinand, right Determinand projection after diverse location sensor acquisition to gradation data and to imagination determinand projection after radiation source ray pass The length value of imagination determinand obtains the determinand absorptance to ray, and then obtains the physical property of determinand.
8. the image processing apparatus at least two non-parallel direction, determinand projected, including:
For according to determinand project after diverse location sensor acquisition to gradation data obtain projection picture ash write music The module of line;
For obtaining the module of the sensing station that flex point is corresponding in grey scale curve according to described grey scale curve;
For the position of the sensing station corresponding according to described flex point and radiation source obtain determinand vertex position can The module that can combine;And
For determine the summit of determinand in the possible combination of the vertex position of described object under test according to described grey scale curve The module of position transitivity.
9. device as claimed in claim 8, also includes filtering high-frequency noise for described grey scale curve is filtered process The module of signal.
10. device as claimed in claim 8, also includes the module for described grey scale curve is asked second dervative.
11. devices as claimed in claim 8, also include for the position according to described radiation source and with described radiation source pair What the line of the sensing station that the described flex point answered is corresponding obtained the vertex position of object under test may the module of combination.
12. devices as claimed in claim 8, also include the possible combination institute structure for the vertex position to described object under test The imaginary determinand become carries out the mould of imagination determinand projection in the case of described radiation source positions and sensing station are constant Block.
13. devices as claimed in claim 12, also include for adopting according to diverse location sensor after imagination determinand projection The gradation data that arrives of collection obtain imagination determinand projection picture grey scale curve and with the grey scale curve of object under test projection picture Relatively, the imaginary determinand vertex position corresponding with the immediate grey scale curve of grey scale curve of object under test projection picture is chosen Module for the vertex position of object under test.
14. devices as claimed in claim 12, also include for according to sensor gained original gradation when not placing determinand Data, determinand is projected after the gradation data that arrives of diverse location sensor acquisition and to radiation source after imagination determinand projection Ray obtains the determinand absorptance to ray through the length value of imagination determinand, and then obtains the mould of the physical property of determinand Block.
CN201610371890.6A 2016-05-30 2016-05-30 Image processing method and image processing apparatus Pending CN106097407A (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201610371890.6A CN106097407A (en) 2016-05-30 2016-05-30 Image processing method and image processing apparatus
KR1020160121190A KR101847730B1 (en) 2016-05-30 2016-09-22 Image processing method and image processing apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610371890.6A CN106097407A (en) 2016-05-30 2016-05-30 Image processing method and image processing apparatus

Publications (1)

Publication Number Publication Date
CN106097407A true CN106097407A (en) 2016-11-09

Family

ID=57230840

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610371890.6A Pending CN106097407A (en) 2016-05-30 2016-05-30 Image processing method and image processing apparatus

Country Status (2)

Country Link
KR (1) KR101847730B1 (en)
CN (1) CN106097407A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108427940A (en) * 2018-04-04 2018-08-21 浙江安精智能科技有限公司 Water fountain effluent intelligent controlling device based on depth camera and its control method

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110456783A (en) * 2018-05-08 2019-11-15 深圳市优必选科技有限公司 A kind of intelligent terminal and its gray scale tracking sense mould group, gray-scale sensor
CN116109540B (en) * 2023-03-22 2023-07-18 智洋创新科技股份有限公司 Image registration fusion method and system based on particle swarm optimization gray curve matching

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201242531Y (en) * 2008-04-17 2009-05-20 清华大学 Scanning imagery system for straight-line track
CN101470082A (en) * 2007-12-27 2009-07-01 同方威视技术股份有限公司 Article detection apparatus and detection method thereof
US20120128225A1 (en) * 2010-11-22 2012-05-24 Canon Kabushiki Kaisha Image processing apparatus, image processing method, and non-transitory storage medium
CN104749648A (en) * 2013-12-27 2015-07-01 清华大学 Multi-energy spectrum static CT apparatus

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6763083B2 (en) 2002-08-30 2004-07-13 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Article screening system
DE102005023376A1 (en) 2005-05-17 2006-11-23 Carl Zeiss Industrielle Messtechnik Gmbh Method and device for determining material boundaries of a test object
JP5613501B2 (en) 2010-08-26 2014-10-22 富士フイルム株式会社 Pipe thickness measuring apparatus and method
JPWO2012108278A1 (en) 2011-02-09 2014-07-03 株式会社Ihi検査計測 X-ray inspection apparatus and method
JP6391365B2 (en) 2014-08-26 2018-09-19 名古屋電機工業株式会社 X-ray inspection apparatus, X-ray inspection method, and X-ray inspection program

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101470082A (en) * 2007-12-27 2009-07-01 同方威视技术股份有限公司 Article detection apparatus and detection method thereof
CN201242531Y (en) * 2008-04-17 2009-05-20 清华大学 Scanning imagery system for straight-line track
US20120128225A1 (en) * 2010-11-22 2012-05-24 Canon Kabushiki Kaisha Image processing apparatus, image processing method, and non-transitory storage medium
CN104749648A (en) * 2013-12-27 2015-07-01 清华大学 Multi-energy spectrum static CT apparatus

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘建敬 等: "一种地磁传感器对电磁舵机的抗干扰方法", 《测试技术学报》 *
孔超 等: "超声相控阵***中高精度相控发射和相控信号并行处理的实现", 《仪器研制》 *
张宇明: "摆式列车倾摆控制信号的实时检测与处理", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》 *
赵曦 等: "双投影法识别客体形状及其物性", 《核电子学与探测技术》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108427940A (en) * 2018-04-04 2018-08-21 浙江安精智能科技有限公司 Water fountain effluent intelligent controlling device based on depth camera and its control method
CN108427940B (en) * 2018-04-04 2023-11-17 徐育 Intelligent control device and control method for water dispenser water outlet based on depth camera

Also Published As

Publication number Publication date
KR20170135628A (en) 2017-12-08
KR101847730B1 (en) 2018-04-10

Similar Documents

Publication Publication Date Title
US8915425B2 (en) Automatic identification system for randomly oriented objects
US10042079B2 (en) Image-based object detection and feature extraction from a reconstructed charged particle image of a volume of interest
US8923587B2 (en) Image reconstruction from limited or incomplete data
US7573039B2 (en) Compton camera configuration and imaging method
CN106353828A (en) Method and device for estimating weight of inspected object in security inspection system
US20080315091A1 (en) Imaging and sensing based on muon tomography
Backholm et al. Simultaneous reconstruction of emission and attenuation in passive gamma emission tomography of spent nuclear fuel
CN106097407A (en) Image processing method and image processing apparatus
CN107688194A (en) Data processing method and data processing equipment
CN103908278A (en) Image reconstruction method and device and medical imaging system
Bücherl et al. A Bayesian method for the evaluation of segmented gamma scanning measurements–Description of the principle
Dekker et al. Evaluation of an iterative reconstruction algorithm for optical CT radiation dosimetry
Ukaegbu et al. Ground Penetrating Radar as a Contextual Sensor for Multi-Sensor Radiological Characterisation
Olson et al. Passive detection of small low-emission sources: two-dimensional numerical case studies
CN113576504B (en) Mu sub-imaging method for medium-low atomic number substance
Marleau et al. An Implementation of Zero Knowledge Confirmation using a Two-dimensional Time-Encoded Imaging System.
Miller et al. Combining radiography and passive measurements for radiological threat localization in cargo
CN111340127B (en) Energy spectrum CT iterative material decomposition method and device based on material clustering
Robinson et al. Imaging for dismantlement verification: Information management and analysis algorithms
D’Alessandro et al. Muography applied to nuclear waste storage sites
Beaumont et al. Imaging of fast neutrons and gamma rays from 252Cf in a heavily shielded environment
Pastena et al. A novel approach to the localization and estimation of radioactivity in contaminated waste packages via imaging techniques
Zeng et al. Image reconstruction algorithm for a spinning strip CZT SPECT camera with a parallel slat collimator and small pixels: CZT SPECT camera using small pixels
AU2014202327B2 (en) Statistical tomographic reconstruction based on measurements of charged particles
US20230228904A1 (en) Muon tomography method and apparatus

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20161109

RJ01 Rejection of invention patent application after publication