CN101923063A - Recognition method of foreign body in grain pile - Google Patents
Recognition method of foreign body in grain pile Download PDFInfo
- Publication number
- CN101923063A CN101923063A CN 201010257470 CN201010257470A CN101923063A CN 101923063 A CN101923063 A CN 101923063A CN 201010257470 CN201010257470 CN 201010257470 CN 201010257470 A CN201010257470 A CN 201010257470A CN 101923063 A CN101923063 A CN 101923063A
- Authority
- CN
- China
- Prior art keywords
- amplitude
- virtual level
- foreign matter
- antenna
- inductive capacity
- 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
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention relates to a recognition method of foreign bodies in a grain pile, which is used for solving the problems of large interference, inaccurate classification and incapability of classifying multiple targets in the prior art. The method comprises the following steps of: firstly positioning a foreign body target in a ground penetrating radar map by utilizing an image processing technology and a Hough transformation method, then carrying out multilayer single-channel inversion on a dielectric constant on the channel where the foreign body target is positioned on the basis of refraction and refection laws, and finally determining the kind of the foreign body according to the dielectric constant subjected to the multilayer single-channel inversion and the corresponding relation of the dielectric constant and the common foreign body to realize the recognition of the foreign body. The method only needs echo amplitude information and does not need echo phase information, so that the algorithm is more easy to realize. Meanwhile, because of the difference of the dielectric constants of different kinds of objects, the method is more accurate for classifying and recognizing the foreign bodies and can simultaneously classify different kinds of foreign bodies in the same area.
Description
Technical field
The present invention relates to a kind of electromagnetic wave detection field that belongs to, particularly a kind of electromagnetic wave detection technology of utilizing is discerned and the technology of classifying the large-scale grain inner a plurality of unusual targets of heap (foreign matter).
Background technology
In the electromagnetic wave detection field, the main processing means that the medium internal object is discerned comprise imaging identification and characteristic variable identification at present.Imaging processing is passed through the processing to the echoed signal of a plurality of sections, obtains the geometric properties of target object, thereby according to geometric properties (mainly being profile) target is differentiated.Promptly differentiate from the two and three dimensions image that is whether objective body is arranged,, determine the position and the shape of objective body from image if there is objective body.Under more satisfactory situation, adopt effective ground penetrating radar imaging algorithm, can obtain imaging effect preferably, can clearly know whether from one-tenth's image to have foreign matter, and can finish the location of foreign matter and the identification of foreign matter shape.But the imaging recognizer depends on accurate system modelling, comprise the form of excitation source signal, the calculating of antenna surface electric current, the calculating of dyadic Green's function, in imaging process, need carry out repeatedly finding the solution of integral equation and calculate with scattered field, calculated amount is huge, is not suitable for engineering and uses.
Identification mainly is to finish according to the extraction that the echoed signal of visiting the ground ground penetrating radar is carried out characteristic variable to the characteristic variable of foreign matter material.Visiting the ground ground penetrating radar detection of a target is to rely on the emission electromagnetic pulse, the echoed signal of receiving target body is finished, be based on the electromagnetic wave propagation theory, extract, and then realize Target Recognition so can carry out clarification of objective by the electromagnetic phase of wave of evaluating objects.At present, the foreign matter method of identification based on characteristic quantity mainly is echo amplitude and direction analytic approach and two kinds of methods of target scattering echo spectrum method for feature analysis.
Echo amplitude and direction analytic approach are according to the reflection in the electromagnetic wave propagation theory, refraction law, the energy of reflection of electromagnetic wave and phase place depend on the difference of reflecting surface levels medium relative dielectric constant, the electromagnetic property difference of both sides, dielectric interface place is big more, and reflection wave is strong more.When the specific inductive capacity of interface top dielectric during greater than the specific inductive capacity of interface layer dielectric, the phase place homophase of reflection wave and incident wave; Otherwise when the specific inductive capacity of interface top dielectric during less than the specific inductive capacity of interface layer dielectric, the phase place of reflection wave and incident wave is anti-phase.This is a foundation of judging interface media of both sides character and attribute, thereby reflection wave amplitude and direction character are the important evidence of differentiation destination media.
The spectrum sigtral response method of target scattering echo is carried out ultra wideband narrow-pulse scanning to the medium internal object, the echoed signal spectral phase of different target has different situations of change, therefore, a present more feasible method is by comparing the phase spectrum of received signal and reference signal, the target material properties being classified.
Echo amplitude and direction analytic approach are simple and practical, but because electromagnetic wave exists chromatic dispersion and relaxation phenomenon in the grain medium, if echo amplitude is not corrected, then analysis result may differ bigger with actual; The echo spectrum method for feature analysis is to the signal to noise ratio (S/N ratio) of echo, and the loss of medium is had relatively high expectations, and when having multiple foreign matter in the medium, can the phase spectral analysis method be lost efficacy.Grain is the dispersive medium that loss is big, signal to noise ratio (S/N ratio) is lower, and that the method phase of echo relative method that can adopt is at present disturbed is big, the classification out of true, can't carry out the plurality of target classification, is difficult to obtain gratifying recognition effect.
Summary of the invention
The purpose of this invention is to provide recognition method of foreign bodies in a kind of grain heap, disturb problem big, that classify out of true, can't carry out the plurality of target classification in order to solve prior art.
For achieving the above object, the solution of the present invention step is as follows:
A) survey with ground penetrating radar, foreign matter in the grain heap is positioned;
B), the grain heap vertically is divided at least two virtual levels from top to bottom successively according to the survey grain heap degree of depth and the foreign matter degree of depth; Based on inversion algorithm and reflected refraction law, draw inversion formula successively:
The ground floor specific inductive capacity:
Layer specific inductive capacity after the second layer reaches:
In the formula (12),
Be the specific inductive capacity of each virtual level, i represents the sequence number of the number of plies from top to bottom,
Be dielectric constant of air, Ai is the reflection echo amplitude of each virtual level, and ki is the linear error calibration factor of each virtual level, and fi is each virtual level amplitude fading factor; In the formula (11), Am is the amplitude of launching electromagnetic wave, and A0 is the first virtual level reflection echo amplitude, and k1 is the linear error calibration factor of first virtual level, and f1 is the first virtual level amplitude fading factor;
C) ground exploring radar antenna is placed on grain heap surface, measures the reflection echo amplitude A i of each virtual level on the corresponding vertical section of ground exploring radar antenna;
D) adopt well logging method to carry out twice measurement and determine parameter f 1 in formula (11), (12), k1, fi, ki; Adopt an emitting antenna and two receiving antennas, when measuring for the first time, emitting antenna and receiving antenna are inserted grain respectively pile first virtual level, according to the electromagnetic wave amplitude fading and the phase shift of measuring, and, determine the first virtual level specific inductive capacity in conjunction with the electromagnetic wave conduction model of proofreading and correct through geometric divergence
, and the amplitude A m of the launching electromagnetic wave that records according to ground penetrating radar, dielectric constant of air
, ground floor reflection echo amplitude A 0 in conjunction with formula (11), is established f1=1, determines k1; When measuring for the second time, emitting antenna and receiving antenna are inserted grain respectively pile second virtual level,, and, determine the second virtual level specific inductive capacity in conjunction with the electromagnetic wave conduction model of proofreading and correct through geonetrical attenuation according to the electromagnetic wave amplitude fading and the phase shift of measuring
, and, establish ki=k1 in conjunction with formula (12), determine amplitude fading calibration factor fi;
E) utilize the specific inductive capacity of each virtual level on formula (12) the iterative computation ground exploring radar antenna correspondence position vertical scan line;
F), judge the character and the kind of foreign matter according to the specific inductive capacity of each virtual level.
This method only needs echo amplitude information, need not phase of echo information, thereby algorithm is easier to realize.Simultaneously, because the difference of variety classes object specific inductive capacity, this method is more accurate to the classification and the identification of foreign matter, and can classify simultaneously to variety classes foreign matter in the same area.Concrete advantage is as follows:
1, this programme adopts image processing techniques and intellectual analysis algorithm, only need to locate and to discern the foreign matter kind by the ground penetrating radar sectional view of scanning, complex analyses processes such as spectrum analysis have been saved, the process that makes grain pile inner foreign matter identification becomes easy relatively, have more intellectuality, thereby improved the efficient of foreign matter identification.
2, traditional pass through the foreign matter identification that phase-comparison method is carried out, can only carry out the coarseness classification to recognition result, as whether being foreign matter, be metal or nonmetal etc., can not effectively discern the kind of foreign matter.This programme distributes by the specific inductive capacity to inverting and carries out numerical analysis, has different specific inductive capacity according to different material, again according to the Hyperbolic Feature of the foreign matter on the ground penetrating radar figure, can accurately estimate the character and the kind of foreign matter.
3, traditional phase-comparison method is when surveying the territory multiple foreign matter is arranged, because influencing each other between the foreign matter echo, that gives scatter echo phase place and reference wave phase place relatively is with very big interference, has had a strong impact on the result of phase bit comparison.This programme at first adopts foreign matter position, Hough transformation location, on the vertical scan line of foreign matter place, successively carry out the specific inductive capacity inverting then, different types of foreign matter is independent of each other in refutation process, if adopt parallel algorithm also can accelerate the complementary operation process, so this programme can adapt to the situation that has multiple foreign matter to exist.
The amplitude A m of launching electromagnetic wave is by ground exploring radar antenna being close to a sheet metal that is used as mirror surface, being measured.The antenna that the described well logging method of step c) adopts, emitting antenna connects sweep generator, and receiving antenna connects oscillograph.In the step a) foreign matter location is comprised: enroll original ground penetrating radar sectional view; Original ground penetrating radar sectional view is carried out pre-service; Adopt the windowing statistic law to extract the foreign matter region; Adopt the level and smooth foreign matter of two-dimensional filtering region image to cut apart; Utilize Hough transformation to carry out target detection and location.
Description of drawings
Fig. 1 is an original ground penetrating radar figure;
Fig. 2 is through pretreated ground penetrating radar figure;
Fig. 3 is the ground penetrating radar figure through Hough transformation;
Fig. 4 is a grain heap layering synoptic diagram;
Fig. 5 is electromagnetic wave is piled a plurality of virtual levels at grain an equivalent propagation model;
Fig. 6 is the well logging method instrumentation plan;
Fig. 7 is that specific inductive capacity is with grain heap depth profile synoptic diagram.
Embodiment
The present invention will be further described in detail below in conjunction with accompanying drawing, and embodiment is as follows:
1, the ground penetrating radar sectional view of enrolling is carried out the image pre-service.Comprise that random noise inhibition, direct wave and Radio frequency interference (RFI) suppress, the straight coupling ripple of antenna suppresses, improves signal to noise ratio (S/N ratio) and figure image intensifying etc.Original ground penetrating radar figure as shown in Figure 1, can see in the middle of the image that what represent near the hyperbolic curve of top is exactly the foreign matter of needs detection, the bright wisp of image top is represented grain face reflection wave (being called direct wave), and the cancellated bright wisp of bottom represents to be embedded in the bar-mat reinforcement in the floor.As seen, except the ground penetrating radar echo of the foreign matter of required detection, also exist the noise that a large amount of powers do not wait among the original ground penetrating radar figure.
2, adopt the windowing statistic law to extract the foreign matter region.This step is carried out Preliminary detection based on the statistical method of energy measuring, extracts the foreign matter region from mass data, is used for apace most nontarget areas being got rid of.To reduce follow-up calculated amount.The window function of institute's windowing can adopt simple rectangular window function.As adopt following window function:
Wherein M is a sampling number in the time window.Window width
Consider calculated amount, need the choose reasonable window width.The model of setting up binary hypothesis test is as follows:
In the formula, N is a number of scans,
Be target echo signal,
Be noise and interference.If
Normal Distribution, average are 0, and variance is
And establish
Each component linear independence.
Get threshold decision and can divide the foreign matter region roughly, can roughly determine the position range of target by the peak value of statistic curve and the window function of choosing, thereby reduce follow-up calculated amount.
3, adopt two-dimensional filtering level and smooth foreign matter region image and carry out image segmentation.The interference and the noise that except that the target hyperbolic curve, also have a large amount of random fluctuations in the foreign matter place 2 dimensional region that the process above-mentioned steps extracts, can influence the effect of Hough transformation, can adopt the two-dimensional filtering smoothed image for this reason, remove unusual bright spot and irrelevant component, improve the signal to noise ratio (S/N ratio) of foreign matter region, so that improve the estimated accuracy of parameter.As figure Fig. 2 is that original ground penetrating radar figure (Fig. 1) is carried out pretreated result, and the result demonstration has filtered out most noise and interference, has only kept the image of the foreign matter that needs detection and has carried out the figure image intensifying.
4, utilize Hough transformation to carry out target detection and location.For taking into account accuracy of detection and calculated amount, improve the traditional Hough transformation method that adopts rim detection, directly Hough transformation is carried out in the foreign matter region that extracts and focus on hyperbolic curve, do not carry out image pre-service such as binaryzation.And according to ground penetrating radar echo characteristics, with the Hough transformation totalizer merely according to positional information change the position into, amplitude information is comprehensively voted.
According to the Hyperbolic Feature of target echo, adopt the Hough transformation formula:
In the formula
In the foreign matter zone
Locational gray-scale value function, promptly corresponding target echo range value.
Expression after image segmentation by unique point
The hyperbolic curve that constitutes, the hyperbolic curve vertex position
The parameter space of forming Hough transformation.
According to
Distribution situation can judge whether that target occurs, when having target to occur
Precipitous peak can occur, mostly be mild fluctuating during driftlessness.Right
The search peak-peak just can obtain hyp summit estimated value
As Fig. 3 is that ground penetrating radar figure (Fig. 2) to treated mistake has carried out the result behind the Hough transformation.Light line in the hyperbolic curve is exactly the result who behind the Hough transformation hyperbolic curve is focused on.If known foreign matter is a layered object, then need not to carry out Hough transformation, only need carry out rim detection, to determine the position and the degree of depth of layered object.
5, determine detecting parameter and divide the iteration number of plies according to the survey grain heap foreign matter degree of depth.After having finished above-mentioned four steps, generally can determine the position (being the degree of depth of foreign matter and the horizontal range of range observation initial point) of foreign matter.According to the degree of depth of foreign matter, determine to need to use the degree of depth of electromagnetic wave detection.As the antenna frequencies used, the time window and whenever sweep sampling number.Foreign matter volume as estimation is bigger, can adopt the less sampling number of per pass, with the quickening iteration time, but can influence the measurement result precision.According to the time window and whenever sweep sampling number and can determine every iteration layer height of dividing, the height of iteration layer is not less than foreign matter size in principle.Window and whenever sweep these two parameters of counting of sampling when the structure level number purpose divide to rely on.Owing to be that depth is measured, adopt emitting antenna and the incorporate 200M antenna of receiving antenna, or the 80M antenna that separates with receiving antenna of emitting antenna, it is directly placed on the grain face, with minimizing propagation attenuation and geonetrical attenuation.Adopt the mode of operation that regularly triggers during measurement, enroll repeatedly the single track echoed signal on the vertical direction.As the window and the acquisition parameter of whenever sweeping 128 sampled points when adopting 40ns, suppose that the average velocity that electromagnetic wave is propagated is that 6cm/ns(notices that electromagnetic wave propagation speed need adopt the known target depth method to demarcate in advance in the grain heap), the then data acquisition degree of depth is 40ns * 6cm/ns=240cm.Then the degree of depth of each virtual level may be defined as 240cm/128 ≈ 1.875cm in the grain heap.General in the large storehouse of standard, according to national food industry standards in 2002 (LS/T 1203-2002), grain identity distance floor level is 6m, therefore can with the time window be adjusted into 600cm/6cm/s=100ns, the degree of depth of each virtual level is 600cm/128 ≈ 5cm.In actual measurement, consider that the resolution of specific inductive capacity is generally lower, can merge facing virtual level mutually, divide to eliminate duplication and accelerate to iterate operation time, but measuring accuracy can corresponding decline.
6, determine the amplitude A m of launching electromagnetic wave.Because reflection and refraction can take place at common dielectric surface in electromagnetic wave, so there is loss in antenna from the reflection echo energy that the ground floor medium receives.In order to determine the amplitude of launching electromagnetic wave, can utilize sheet metal as mirror surface, and make antenna be close to sheet metal.Like this, the electromagnetic wave energy that emitting antenna sends almost all is reflected back toward receiving antenna, and the energy attenuation that does not exist geometric divergence to cause, and received echo amplitude can be similar to thinks the amplitude A m of transmitted wave.
7, ground exploring radar antenna is placed on grain heap surface, measure the reflection echo amplitude A i of each virtual level on the corresponding vertical section of ground exploring radar antenna; And carry out the method for inversion successively, detailed process is as follows, and dielectric stratifying as shown in Figure 4.Because the grain heap is a bulk solid, scattering can take place when propagating in electromagnetic wave in the grain heap.Echo on certain sampled point that receiving antenna receives in the vertical is the summation of the scatter echo on this sampled point corresponding virtual aspect, and this scatter echo summation can be similar to the reflection echo of regarding this aspect as.The scattering wave of no show receiving antenna and refraction wave will synthesize the refraction wave of this aspect in the lump on this aspect, and through after such abstract, the propagation model of electromagnetic wave in the grain heap as shown in Figure 5.Ignored the above repeatedly reflection wave of secondary among the figure,
Represent respectively
The reflection echo amplitude and the time of arrival of individual medium interface are poor.
Antenna is in the position
The place transmits
, received signal can be similar to the stack of regarding each aspect echo as:
Wherein L is the number of plies that medium is divided,
With
Amplitude and time delay for each aspect echo.Because each layer of grain heap medium is approximate even, and does not generally have other targets in the medium, and be to survey in the storehouse, obtained good treatment and inhibition by pre-service noise and interference, so the echo that receives can be used this information model.
Only consider primary reflection, based on the refraction and the principle of reflection of plane wave, the relational expression that can obtain specific inductive capacity and aspect echo is:
Wherein
Represent the interface of air and grain.
Be
The average attenuation of individual aspect echo,
Be
The specific inductive capacity of individual medium layer,
Be the specific inductive capacity of air, get 1 usually.Based on last relation of plane, the recursion formula that can obtain specific inductive capacity is:
Can successively calculate the specific inductive capacity of each layer by iteration, but need initial parameter, and the precision of initial parameter will directly have influence on the recursion precision.Initial parameter obtains by calibration accurately.
Successively inversion method is generally supposed antenna plane of departure ripple, and grain heap is assumed to smooth surface, do not have the layered medium of loss and chromatic dispersion, but the time to this near field of silo target measurement, antenna institute launching electromagnetic wave can not simply be considered as plane wave, and the grain in the grain heap is typical dispersive medium, for improving measuring accuracy, calibration factor is introduced inversion formula successively:
(12)
In the following formula, Am is an inverting echo amplitude calibration value, and Ai is the reflection echo amplitude of each aspect,
Be the specific inductive capacity of each layer medium, i represents the sequence number of medium.Consider that what calculate is one dimension specific inductive capacity sequence, so omitted position x and depth z in the formula, two dimension is specific inductive capacity continuously
Be reduced to
.fi with ki the calibration factor of corresponding different layers, fi is used to reduce the decay of each layer medium, the influence of electromagnetic how much factors such as diffusion, and ki is the linear error calibration factor, be again a total regulatory factor simultaneously, can reduce the influence that time delay and specific inductive capacity evaluated error etc. are brought.
8, this method is the specific inductive capacity by each layer of iterative computation, so the formula of deriving above (11), (12) need initial parameter.The 6th step was determined can transmitted wave amplitude A m, previous step can determine to measure the reflection wave amplitude A 0 of ground floor medium, and other initial parameter need adopt well logging method to carry out twice measurement, adopt an emitting antenna and two receiving antennas, when measuring for the first time, emitting antenna and receiving antenna are inserted grain respectively pile first virtual level, grain heap moisture measurement is because measurement range is little, belong to the near field working method, therefore adopt elder generation with the plane wave model specific inductive capacity of deriving, and then carried out the method that geonetrical attenuation is proofreaied and correct.
Employing adopts the method principle of well logging method derivation specific inductive capacity as follows through the plane wave model of geometry correction:
Regard plane wave as with the propagation of electromagnetic wave in the stratum is approximate, promptly
In the formula, k is a wave number.Then for the single-emission and double-receiving electrode system of being made up of T-R1-R2, as shown in Figure 6, the ratio of the electric field intensity that two receiving antennas are received is:
Then amplitude fading is:
Phase shift is:
Can release thus:
(5)
Then the specific inductive capacity of medium and conductivity are:
Because electromagnetic wave is actually a kind of spherical wave, it is a kind of approximate regarding it as plane wave, and spherical wave has geonetrical attenuation, therefore need carry out geonetrical attenuation and proofread and correct.Definition:
Then the geometric divergence correction factor is As=44.65+tp1 (1.324+0.184tp1).
Amplitude fading after geometric divergence is proofreaied and correct is EATTc=EATT-As, after the correction
=(EATT-As)/8.686.
Will
Substitution (6) formula can be calculated the specific inductive capacity at measuring point place.
When using above method to carry out this pacing amount, notice that antenna did not insert deeply (had better not surpass ground floor) in the well, in order to avoid initial specific inductive capacity
Calculate inaccurate.
Obtained the specific inductive capacity of ground floor utilizing well logging method
, after launching electromagnetic wave amplitude A m and the reflection wave amplitude A 0, just can utilize (11) formula to calculate the specific inductive capacity factor k1 that calibrates for error.Attention: supposed f1=1 here, this be because: antenna has adopted mask pattern, and with grain face tight coupling, both do not had the decay of echo amplitude, does not have the influence of factors such as how much diffusions yet, the echo of scattering almost can all be received the antenna reception.
Adopt well logging method once more, will transmit and receive in the virtual second layer that antenna is inserted into division, Measuring Dielectric Constant reads the amplitude rreturn value of visiting place, ground penetrating radar echo insertion point, the end according to the degree of depth of inserting then.The linear error calibration factor k that obtained according to the last step and visit the echo amplitude value that end ground penetrating radar receives, substitution stepping type (12), ask amplitude fading calibration factor f1, consider it is same medium, each layer decling phase together can be f1 as each layer amplitude fading factor fk, i.e. f1=fk (k=2,, n).The initial value and the parameter of formula (12) have all been determined like this.
9, behind the initial parameter of having determined iterative (17) and iteration coefficient, can be finally inversed by each layer specific inductive capacity (comprising the specific inductive capacity of foreign matter place layer) according to visiting the echo amplitude value that end ground penetrating radar is gathered.
10, the distribution plan when drawing the specific inductive capacity single-channel scanning according to the numerical value at catastrophe point place, is judged the character and the kind of foreign matter, and specific inductive capacity, the conductivity of various common foreign matters are as shown in table 1, can judge the foreign matter kind according to table 1.The specific inductive capacity inversion result distribution plan of drawing as shown in Figure 7.As can be seen from the figure, the specific inductive capacity obviously different (being almost 1) of other positions in the specific inductive capacity at foreign matter place and the medium can be inferred the kind (is that the foreign matter at 1 place is the old soldi of burying underground in advance as specific inductive capacity) of foreign matter according to the feature of foreign matter specific inductive capacity in view of the above.
Claims (5)
1. recognition method of foreign bodies during a grain is piled is characterized in that step is as follows:
A) survey with ground penetrating radar, foreign matter in the grain heap is positioned;
B), the grain heap vertically is divided at least two virtual levels from top to bottom successively according to the survey grain heap degree of depth and the foreign matter degree of depth; Based on inversion algorithm and reflected refraction law, draw inversion formula successively:
The ground floor specific inductive capacity:
Layer specific inductive capacity after the second layer reaches:
In the formula (12),
Be the specific inductive capacity of each virtual level, i represents the sequence number of the number of plies from top to bottom,
Be dielectric constant of air, Ai is the reflection echo amplitude of each virtual level, and ki is the linear error calibration factor of each virtual level, and fi is each virtual level amplitude fading factor; In the formula (11), Am is the amplitude of launching electromagnetic wave, and A0 is the first virtual level reflection echo amplitude, and k1 is the linear error calibration factor of first virtual level, and f1 is the first virtual level amplitude fading factor;
C) ground exploring radar antenna is placed on grain heap surface, measures the reflection echo amplitude A i of each virtual level on the ground exploring radar antenna correspondence position vertical scan line;
D) adopt well logging method to carry out twice measurement and determine parameter f 1 in formula (11), (12), k1, fi, ki; Adopt an emitting antenna and two receiving antennas, when measuring for the first time, emitting antenna and receiving antenna are inserted first virtual level respectively, according to the electromagnetic wave amplitude fading and the phase shift of measuring, and, determine the first virtual level specific inductive capacity in conjunction with the electromagnetic wave conduction model of proofreading and correct through geonetrical attenuation
, and the amplitude A m of the launching electromagnetic wave that records according to ground penetrating radar, dielectric constant of air
, ground floor reflection echo amplitude A 0 in conjunction with formula (11), is established f1=1, determines k1; When measuring for the second time, emitting antenna and receiving antenna are inserted second virtual level respectively,, and, determine the second virtual level specific inductive capacity in conjunction with the electromagnetic wave conduction model of proofreading and correct through geonetrical attenuation according to the electromagnetic wave amplitude fading and the phase shift of measuring
, and, establish ki=k1 in conjunction with formula (12), determine amplitude fading calibration factor fi;
E) utilize the specific inductive capacity of each virtual level on formula (12) the iterative computation ground exploring radar antenna correspondence position vertical scan line;
F), judge the character and the kind of foreign matter according to the specific inductive capacity of each virtual level.
2. recognition method of foreign bodies in a kind of grain heap according to claim 1 is characterized in that the amplitude A m of launching electromagnetic wave is by ground exploring radar antenna being close to a sheet metal that is used as mirror surface, being measured.
3. recognition method of foreign bodies in a kind of grain heap according to claim 1 is characterized in that, the antenna that the described well logging method of step c) adopts, and emitting antenna connects sweep generator, and receiving antenna connects oscillograph.
4. recognition method of foreign bodies in a kind of grain heap according to claim 1 is characterized in that, in the step a) foreign matter location is comprised: enroll original ground penetrating radar sectional view; Original ground penetrating radar sectional view is carried out pre-service; Adopt the windowing statistic law to extract the foreign matter region; Adopt the level and smooth foreign matter of two-dimensional filtering region image to cut apart; Utilize Hough transformation to carry out target detection and location.
5. according to recognition method of foreign bodies in each described a kind of grain heap among the claim 1-4, it is characterized in that layering is grain to be piled every 5cm be divided into one deck totally 128 layers in the step b).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102574708A CN101923063B (en) | 2010-08-19 | 2010-08-19 | Recognition method of foreign body in grain pile |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102574708A CN101923063B (en) | 2010-08-19 | 2010-08-19 | Recognition method of foreign body in grain pile |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101923063A true CN101923063A (en) | 2010-12-22 |
CN101923063B CN101923063B (en) | 2012-11-14 |
Family
ID=43338117
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102574708A Expired - Fee Related CN101923063B (en) | 2010-08-19 | 2010-08-19 | Recognition method of foreign body in grain pile |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101923063B (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103439707A (en) * | 2013-08-27 | 2013-12-11 | 中南大学 | Method for migration imaging of windowing distance of ground penetrating radar |
CN104755956A (en) * | 2012-11-02 | 2015-07-01 | 艾普斯科有限公司 | Method and apparatus for inspection of cooling towers |
CN105005042A (en) * | 2015-07-27 | 2015-10-28 | 河南工业大学 | Ground penetrating radar underground target locating method |
CN105403883A (en) * | 2015-10-29 | 2016-03-16 | 河南工业大学 | Ground penetrating radar underground target position detection method |
CN107621634A (en) * | 2017-10-27 | 2018-01-23 | 河南科技大学 | A kind of potted-seedling transplanting plants deep monitoring system and method |
CN108876079A (en) * | 2017-05-15 | 2018-11-23 | 北京国双科技有限公司 | A kind of method and device of sequence |
CN109444875A (en) * | 2018-09-19 | 2019-03-08 | 厦门精益远达智能科技有限公司 | Silo monitoring method, device, equipment and system based on radar migration imaging |
CN110009000A (en) * | 2019-03-11 | 2019-07-12 | 东北大学 | The grain heap object detection method of sorting algorithm is improved based on ADABOOST+SVM |
CN110308498A (en) * | 2019-06-24 | 2019-10-08 | 天津天航智远科技有限公司 | Meteorological Observation System and method based near space dirigible |
CN112505646A (en) * | 2020-11-18 | 2021-03-16 | 安洁无线科技(苏州)有限公司 | Foreign matter shielding judgment method and system based on millimeter wave radar |
CN113884550A (en) * | 2021-09-09 | 2022-01-04 | 水利部交通运输部国家能源局南京水利科学研究院 | Method for measuring dielectric constant of soil body/concrete in field layering manner and application thereof |
CN116755088A (en) * | 2023-08-09 | 2023-09-15 | 中国科学院空天信息创新研究院 | Granary depth and foreign matter detection and imaging method based on radar |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5814731A (en) * | 1997-01-28 | 1998-09-29 | Alexander; Alton Michel | Ultrasonic scanning apparatus for nondestructive site characterization of structures using a planar based acoustic transmitter and receiver in a rolling pond |
US6986294B2 (en) * | 2000-02-17 | 2006-01-17 | Bintech Lllp | Bulk materials management apparatus and method |
CN101738642A (en) * | 2008-11-10 | 2010-06-16 | 同济大学 | Method for processing ground penetrating radar data |
-
2010
- 2010-08-19 CN CN2010102574708A patent/CN101923063B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5814731A (en) * | 1997-01-28 | 1998-09-29 | Alexander; Alton Michel | Ultrasonic scanning apparatus for nondestructive site characterization of structures using a planar based acoustic transmitter and receiver in a rolling pond |
US6986294B2 (en) * | 2000-02-17 | 2006-01-17 | Bintech Lllp | Bulk materials management apparatus and method |
CN101738642A (en) * | 2008-11-10 | 2010-06-16 | 同济大学 | Method for processing ground penetrating radar data |
Non-Patent Citations (1)
Title |
---|
《河南工业大学学报(自然科学版)》 20080831 付麦霞等 《应用探地雷达技术探测粮仓异物》 第50-55页 1-5 第29卷, 第4期 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104755956A (en) * | 2012-11-02 | 2015-07-01 | 艾普斯科有限公司 | Method and apparatus for inspection of cooling towers |
CN103439707A (en) * | 2013-08-27 | 2013-12-11 | 中南大学 | Method for migration imaging of windowing distance of ground penetrating radar |
CN105005042A (en) * | 2015-07-27 | 2015-10-28 | 河南工业大学 | Ground penetrating radar underground target locating method |
CN105005042B (en) * | 2015-07-27 | 2017-10-10 | 河南工业大学 | A kind of GPR buried target localization method |
CN105403883A (en) * | 2015-10-29 | 2016-03-16 | 河南工业大学 | Ground penetrating radar underground target position detection method |
CN108876079A (en) * | 2017-05-15 | 2018-11-23 | 北京国双科技有限公司 | A kind of method and device of sequence |
CN107621634B (en) * | 2017-10-27 | 2023-10-27 | 河南科技大学 | Pot seedling transplanting and planting depth monitoring system and method |
CN107621634A (en) * | 2017-10-27 | 2018-01-23 | 河南科技大学 | A kind of potted-seedling transplanting plants deep monitoring system and method |
CN109444875A (en) * | 2018-09-19 | 2019-03-08 | 厦门精益远达智能科技有限公司 | Silo monitoring method, device, equipment and system based on radar migration imaging |
CN110009000A (en) * | 2019-03-11 | 2019-07-12 | 东北大学 | The grain heap object detection method of sorting algorithm is improved based on ADABOOST+SVM |
CN110009000B (en) * | 2019-03-11 | 2022-09-02 | 东北大学 | Grain pile target detection method based on ADABOOST + SVM improved classification algorithm |
CN110308498A (en) * | 2019-06-24 | 2019-10-08 | 天津天航智远科技有限公司 | Meteorological Observation System and method based near space dirigible |
CN112505646A (en) * | 2020-11-18 | 2021-03-16 | 安洁无线科技(苏州)有限公司 | Foreign matter shielding judgment method and system based on millimeter wave radar |
CN112505646B (en) * | 2020-11-18 | 2024-04-12 | 安洁无线科技(苏州)有限公司 | Foreign matter shielding judging method and system based on millimeter wave radar |
CN113884550A (en) * | 2021-09-09 | 2022-01-04 | 水利部交通运输部国家能源局南京水利科学研究院 | Method for measuring dielectric constant of soil body/concrete in field layering manner and application thereof |
CN113884550B (en) * | 2021-09-09 | 2022-05-10 | 水利部交通运输部国家能源局南京水利科学研究院 | Method for measuring dielectric constant of soil body/concrete in field layering manner and application thereof |
CN116755088A (en) * | 2023-08-09 | 2023-09-15 | 中国科学院空天信息创新研究院 | Granary depth and foreign matter detection and imaging method based on radar |
CN116755088B (en) * | 2023-08-09 | 2023-11-17 | 中国科学院空天信息创新研究院 | Granary depth and foreign matter detection and imaging method based on radar |
Also Published As
Publication number | Publication date |
---|---|
CN101923063B (en) | 2012-11-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101923063B (en) | Recognition method of foreign body in grain pile | |
CN104020495B (en) | Automatic underground pipeline parameter recognizing method on basis of ground penetrating radar | |
CN101915789B (en) | Electromagnetic wave detection method for grain pile water content | |
CN110866545A (en) | Method and system for automatically identifying pipeline target in ground penetrating radar data | |
CN109581317B (en) | Corner target positioning method based on echo peak matching | |
CN107390213B (en) | A kind of time lag curve extracting method of the Ground Penetrating Radar record section based on sliding window | |
Potin et al. | An abrupt change detection algorithm for buried landmines localization | |
KR101897763B1 (en) | Method and apparatus for distance measurement using radar | |
CN104730148A (en) | Metal material inner inclusion three-dimensional reconstruction method based on ultrasonic testing technology | |
CN111025286A (en) | Ground penetrating radar map self-adaptive selection method for water damage detection | |
Abed et al. | Echo amplitude normalization of full-waveform airborne laser scanning data based on robust incidence angle estimation | |
Chen et al. | Probabilistic conic mixture model and its applications to mining spatial ground penetrating radar data | |
CN115685096B (en) | Secondary radar side lobe suppression method based on logistic regression | |
CN115561753A (en) | Method, device, equipment and storage medium for determining underground target | |
CN111359913A (en) | Method for sorting ores through laser radar | |
CN105738881A (en) | Offshore Envisat satellite height finding waveform resetting method based on waveform classification | |
CN107589420A (en) | A kind of interior of articles component detection method, apparatus and system | |
CN110488273A (en) | A kind of vehicle tracking detection method and device based on radar | |
Gonzalez-Huici | A strategy for landmine detection and recognition using simulated GPR responses | |
Milisavljević et al. | Improving mine recognition through processing and Dempster–Shafer fusion of ground-penetrating radar data | |
Li et al. | Artifact suppression of back-projection algorithm under multiple buried objects | |
CN110412609B (en) | Multi-pulse laser radar target detection method | |
Cao et al. | Evaluation of pavement system based on ground-penetrating radar full-waveform simulation | |
CN101545969B (en) | Method for evaluating point target resolution of side-glance SAR | |
Tantum et al. | Comparison of algorithms for land mine detection and discrimination using ground penetrating radar |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121114 Termination date: 20140819 |
|
EXPY | Termination of patent right or utility model |