CN105759252B - A kind of insect dimension measurement method based on multifrequency scattering modeling - Google Patents
A kind of insect dimension measurement method based on multifrequency scattering modeling Download PDFInfo
- Publication number
- CN105759252B CN105759252B CN201610171604.1A CN201610171604A CN105759252B CN 105759252 B CN105759252 B CN 105759252B CN 201610171604 A CN201610171604 A CN 201610171604A CN 105759252 B CN105759252 B CN 105759252B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- dbf
- mtd
- frequency
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Catching Or Destruction (AREA)
Abstract
The invention discloses a kind of insect dimension measurement method based on multifrequency scattering modeling, this method can measure simultaneously insect body length and body it is wide, can the further quality of inverting insect and the wide ratio of the long body of body based on this, the insect size range that this method can cover greatly increases compared with traditional entomological radar dimension measurement method, therefore also enhances the ability that entomological radar identifies to caste.
Description
Technical field
The invention belongs to entomological radar technical field, and in particular to a kind of insect dimensional measurement based on multifrequency scattering modeling
Method.
Background technology
Radar is arising at the historic moment during the Second World War, is target to be detected using electromagnetic wave and the electronic equipment of ranging.Army
There are many both echo points of non-military objective nor atmosphere elements with being found that very early with weather radar in the air, confirm these later
Echo is migrating property insect.It is good electromagnetic-wave reflector because the water content in insect bodies is more than 80%, this is also elder brother
Worm radar is able to the important prerequisite successfully developed and applied.Hereafter, entomologists all over the world using entomological radar to moving
Fly property insect and carry out a series of researchs.The U.S., Australia and Britain have been owned by the entomological radar system of oneself at present, I
State is also configured with entomological radar system in many places such as Shandong, Henan, Liaoning.
Traditional entomological radar is limited by factors such as working system, systemic-function, algorithm and indexs, can not accurately obtain elder brother
Worm biological parameter.Traditional entomological radar is usually operated at Rayleigh region, and only single frequency band, therefore can only obtain radar surveying
To entomological radar scattering resonance state compared with the empirical value of Rayleigh region, so as to extrapolate the quality of insect, and further root
According to quality inverting build or classified.Relatively simple, particularly larger to the build insect of information content that this method obtains,
Its scattering can enter resonance region, and this can largely effect on the precision of traditional build inversion method.Therefore the insect of traditional entomological radar
Dimension measurement method can not often cover wider insect size range, limit the validity of entomological radar.Actually insect
Shape parameter also include the wide ratio of the long body of body, this between different types of insect there is also notable difference, if can measure simultaneously
The wide ratio of the long body of body and quality of insect, then can greatly improve the precision that entomological radar identifies to caste.
The content of the invention
In view of this, can it is an object of the invention to provide a kind of insect dimension measurement method based on multifrequency scattering modeling
So that resonance region scattering modeling is carried out to insect using rotational symmetry dielectric ellipsoid, and it is wide by the body length and body of radar inverting insect
It is significant to entomological radar target recognition and classification Deng biological parameter.
A kind of insect dimension measurement method based on multifrequency scattering modeling, comprises the following steps:
Step 1, multiband scattering resonance state RCS data acquisitions:
Insect to be measured is detected using multiband radar system, obtains the scattered signal of N number of frequency altogether, wherein, each wave band
Frequency number be not less than 10;Log-frequency corresponding to each frequency is dBfn=10log10(fn), fnFor the frequency of n-th of frequency
Value, n=1~N;For n-th of frequency, Complete polarimetry matrix corresponding to insect to be measured is obtained:
Wherein,For horizontal polarized components,For perpendicular polarisation components,WithTo intersect
Polarization components;
By radar Complete polarimetry matrix Sn, substitute into following formula:
Obtain the measurement data σ of the scattering resonance state RCS maximums of n-th of frequencymax(dBfn);N number of Frequency point is traveled through,
Obtain the measurement data of scattering resonance state RCS maximums corresponding to each Frequency point;
Step 2, the search based on least square meaning obtain parameter dBf respectively0With K estimationWithSpecifically
For:
1), difference setup parameter dBf0The region of search and scouting interval with K;Wherein, dBf0Represent scattering resonance state RCS
Log-frequency during maximum, K are the variables for the absolute intensity for influenceing scattering, represent target in dBf0When scattering resonance state
RCS;Using electromagnetic simulation software simulating insects in the scattering properties of different-waveband, obtained σmaxOn rising for log-frequency dBf
Volt property is by following piecewise function approximate expression:
2), for log-frequency dBf corresponding to each Frequency point, Selecting All Parameters dBf0With K in the respective region of search
One group of parameter, find σ under current logarithmic frequency from formula (5)maxCorresponding function expression, by current logarithmic frequency, choosing
The parameter dBf taken0The σ found is substituted into K valuesmaxIn function expression, the σ under current logarithmic frequency is obtainedmaxCalculated value, by this
The σ for the respective frequencies point that calculated value obtains with step 1maxMeasured value subtracts each other, and obtains error amount, then error amount is squared;Time
All log-frequency points are gone through, obtain square of error amount corresponding to each Frequency point, finally, to square asking for all error amounts
With;
3), using method 2), according to the scouting interval of setting, dBf is constantly changed in the respective region of search0And K,
Error amount square corresponding to obtaining and value, choose corresponding one group of parameter dBf when making and being worth minimum0With K estimationWith
Step 3, insect Size calculation:
By estimated resultSubstitute into following formula:
Obtain the length 2a of insect;
WillIt is updated in following formula:
Obtain insect volume V;
Insect width 2b is calculated according to the cubature formula of ellipsoid simultaneously:
The multiband radar system covering X-band, Ku wave bands and Ka wave bands, each band are not less than 2GHz.
The present invention has the advantages that:
A kind of insect dimension measurement method based on multifrequency scattering modeling of the present invention, this method can measure insect simultaneously
Body length and body it is wide, based on this can the further quality of inverting insect and the wide ratio of the long body of body, the insect that this method can cover
Size range greatly increases compared with traditional entomological radar dimension measurement method, therefore also enhances entomological radar to insect kind
The ability of class identification.
Brief description of the drawings
Fig. 1 is single portion's entomological radar Vertical Observation geometric representation.
Fig. 2 is that insect resonance region scatters fluctuation characteristic and fitting.
Fig. 3 is the emulation of different size insect broadband radar signals product and matched curve.
Embodiment
The present invention will now be described in detail with reference to the accompanying drawings and examples.
For the insect that build is larger, such as loxostege sticticalis (body grows 1~1.2cm), in X, Ku, Ka wave band, it is normally at humorous
Shake area, can not carry out Analytic modeling, therefore to insect with rotational symmetry dielectric ellipsoid model modeling, and it is based on rotational symmetry medium
Actual scattering properties of the ellipsoid in resonance region carries out dimensional measurement.
Insect is observed vertically upward using single portion's radar, as shown in Figure 1.Generally it is considered that Migrating Insects are all
Flat winged motion is being done, i.e., velocity is parallel to the ground, if insect bodies axle is in radar antenna polarization reception planeInterior
Projection and radar horizon polarised direction vectorAngle be α, radar vertical polarization directions vector representation is
Generally, polarization radar observes that the Complete polarimetry matrix of insect can be expressed as:
Wherein ShhFor horizontal polarized components, SvvFor perpendicular polarisation components, ShvAnd SvhFor cross polar component.
Assuming that polarization response matrix such as (1) formula, rotational symmetry body insect that radar Polarimetry obtains have Shv=
Svh, then the maximum polarization intensity that receives of radar can be calculated as:
Because of the maximum σ of the Radar Cross Section (RCS) of insectmaxIt is expressed as:
σmax=| Smax|2 (3)
, can approximation power function piecewise fitting σ based on the WB-RCS emulation to rotational symmetry dielectric ellipsoid targetmax
The fluctuation characteristic of resonance region, represent for convenience, define decibel frequency dBf first:
DBf=10 × log10(f/109) (4)
Wherein f is frequency, unit Hz.As shown in Figure 2.Utilize electromagnetic simulation software simulating insects dissipating in different-waveband
Characteristic is penetrated, is found based on a large amount of emulation data analyses, σmaxFluctuation characteristic on dBf can use following piecewise function expression formula
It is approximate:
σ in above formulamaxUnit be dBsm, dBf0Decibel frequency during RCS maximums is represented, it is with ellipsoid semi-major axis a's
Relation can using approximate representation as:
In other words, when body shaft length changes, σmaxFluctuation characteristic in resonance region is present along the flat of frequency axis
Move.K is the variable for the absolute intensity for influenceing scattering, and unit dBsm, it represents target in dBf0When RCS.It is relative in insect
In the case that dielectric constant determines, K and insect volume V (unit mm3) between relation can using approximate representation as:
V=(58510 × 10K/10-2.89)4/3 (7)
In other words, when insect volume changes, σmaxThe translation along RCS axles be present in the fluctuation characteristic in resonance region.
Therefore, the invention provides a kind of insect dimension measurement method based on multifrequency scattering modeling, comprise the following steps.
Step 1, multifrequency RCS data acquisitions:
Assuming that multiband radar system launches the signal of N number of frequency simultaneously, it is desirable to frequency covers tri- wave bands of X, Ku and Ka,
Each band is not less than 2G, and each wave band frequency number is not less than 10.Corresponding decibel frequency is dBfj, j=1~N.It is right
Each frequency, we can obtain Complete polarimetry matrix corresponding to single insect:
Utilize polarization radar collision matrix measured value Sn, the rcs measurement of N number of frequency can be obtained by substituting into (2) formula and (3) formula
Data σmax(dBfj), j=1~N.
Step 2, utilizes (5) formula, and the search based on least square meaning can be obtained to parameter dBf0With K estimation
With
1) dBf, is set0It is spaced with K traversal section and traversal, dBf0Can be respectively according to the frequency model for obtaining data with K
Enclose and determine traversal section with amplitude range;
2), for log-frequency dBf corresponding to each Frequency point, one group of parameter dBf being chosen in interval range0,p
And Kq, σ under current logarithmic frequency is found from formula (5)maxCorresponding function expression, by current logarithmic frequency, the ginseng chosen
Number dBf0The σ found is substituted into K valuesmaxIn function expression, the σ under current logarithmic frequency is obtainedmaxCalculated value, by the calculated value
The measured value of the respective frequencies point obtained with step 1 subtracts each other, and obtains error amount, then error amount is squared;Travel through all logarithms
Frequency point, square of error amount corresponding to each Frequency point is obtained, finally, square summation to all error amounts:
3), using method 2), according to the scouting interval of setting, dBf is constantly changed in the respective region of search0And K,
Error amount square corresponding to obtaining and value, choose corresponding one group of parameter dBf when making and being worth minimum0With K estimationWith
Step 3, insect Size calculation:
By estimated resultSubstitute into following formula:
Obtain the length 2a of insect;
WillIt is updated in following formula:
Obtain insect volume V;
Insect width 2b is calculated according to the cubature formula of ellipsoid simultaneously:
Embodiment:
To verify aforementioned dimensions measuring method, to four kinds of typical Migrating Insects (loxostege sticticalis, bollworm, mythimna separata and oryza chinensis)
Be modeled using rotational symmetry dielectric ellipsoid, if be sized respectively 12mm × 2mm, 20mm × 4mm, 24mm × 4mm and
30mm × 6mm (major axis × short axle).The WB-RCS of four different size medium ellipsoids is emulated using CST.And provide corresponding point
Section Function Fitting result, as shown in Figure 3.It can be seen that the piecewise function of given form can meet to different size insect broadbands
RCS emulates the fitting of data.
Point frequency RCS emulation data based on tri- wave bands of X, Ku and Ka, more frequency dispersions are based on using one kind of the present invention
The insect dimension measurement method of modeling is penetrated, completes the insect dimensional measurement emulation under the parameter, idiographic flow is as follows:
1. assume that insect body axis direction is parallel with horizontal polarization directions, can be in the hope of insect RCS according to (2) formula and (3) formula
The maximum that polarizes (unit dBsm) is as shown in table 1.
The insect multifrequency RCS data of table 1
2. the RCS data according to given by table 1, substitute into formula (5) formula and carry out two-dimensional search, parameter can be obtained
Estimated result, as shown in table 2:
The parameter of curve estimated result of table 2
3. the parameter estimation result in table 2 is substituted into (6) (7) (11) formula, can be in the hope of the shape parameter of insect, such as table 3
It is shown.
The insect dimension measurement result of table 3
It can be obtained based on above-mentioned simulation result to draw a conclusion:
To various sizes of big insect, making curve matching using resonance region multifrequency point RCS data can effective inverting insect
Shape parameter.
By simulation result it can be seen that utilizing this high efficiency and validity that modeling method is scattered based on multifrequency.Utilize
This method can realize the target size measurement of entomological radar.
In summary, presently preferred embodiments of the present invention is these are only, is not intended to limit the scope of the present invention.
Within the spirit and principles of the invention, any modification, equivalent substitution and improvements made etc., it should be included in the present invention's
Within protection domain.
Claims (2)
1. a kind of insect dimension measurement method based on multifrequency scattering modeling, it is characterised in that comprise the following steps:
Step 1, multiband scattering resonance state RCS data acquisitions:
Insect to be measured is detected using multiband radar system, obtains the scattered signal of N number of frequency altogether, wherein, the frequency of each wave band
Points are not less than 10;Log-frequency corresponding to each frequency is dBfn=10log10(fn), fnFor the frequency values of n-th of frequency, n
=1~N;For n-th of frequency, Complete polarimetry matrix corresponding to insect to be measured is obtained:
Wherein,For horizontal polarized components,For perpendicular polarisation components,WithFor cross polarization
Component;
By radar Complete polarimetry matrix Sn, substitute into following formula:
<mrow>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>dBf</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>|</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mo>|</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>|</mo>
<mo>+</mo>
<msqrt>
<mrow>
<mo>|</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
<mi>n</mi>
</msubsup>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>4</mn>
<mo>|</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>h</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>|</mo>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
Obtain the measurement data σ of the scattering resonance state RCS maximums of n-th of frequencymax(dBfn);N number of Frequency point is traveled through, is obtained each
The measurement data of scattering resonance state RCS maximums corresponding to Frequency point;
Step 2, the search based on least square meaning obtain parameter dBf respectively0With K estimationWithSpecially:
1), difference setup parameter dBf0The region of search and scouting interval with K;Wherein, dBf0Represent that scattering resonance state RCS is maximum
Log-frequency during value, K are the variables for the absolute intensity for influenceing scattering, represent target in dBf0When scattering resonance state RCS;Profit
With electromagnetic simulation software simulating insects in the scattering properties of different-waveband, Radar Cross Section maximum σ is obtainedmaxOn right
Number frequency dBf fluctuation characteristic is by following piecewise function approximate expression:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>max</mi>
</msub>
<mo>=</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>d</mi>
<mi>B</mi>
<mi>f</mi>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<mrow>
<msub>
<mi>dBf</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<mi>K</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>3.3</mn>
<mo>&times;</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>d</mi>
<mi>B</mi>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>dBf</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>1.5</mn>
</msup>
<mo>+</mo>
<mi>K</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>d</mi>
<mi>B</mi>
<mi>f</mi>
<mo>&le;</mo>
<msub>
<mi>dBf</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>|</mo>
<mi>d</mi>
<mi>B</mi>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>dBf</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mn>3.3</mn>
<msup>
<mo>|</mo>
<mn>1.6</mn>
</msup>
</mrow>
<mn>0.9623</mn>
</mfrac>
<mo>+</mo>
<mi>K</mi>
<mo>-</mo>
<mn>3.3</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>dBf</mi>
<mn>0</mn>
</msub>
<mo><</mo>
<mi>d</mi>
<mi>B</mi>
<mi>f</mi>
<mo>&le;</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>dBf</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mn>3.3</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>|</mo>
<mi>d</mi>
<mi>B</mi>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>dBf</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mn>3.3</mn>
<msup>
<mo>|</mo>
<mn>1.2</mn>
</msup>
</mrow>
<mn>1.5513</mn>
</mfrac>
<mo>+</mo>
<mi>K</mi>
<mo>-</mo>
<mn>3.3</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>dBf</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mn>3.3</mn>
</mrow>
<mo>)</mo>
</mrow>
<mo><</mo>
<mi>d</mi>
<mi>B</mi>
<mi>f</mi>
<mo>&le;</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>dBf</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mn>4.4</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.35</mn>
<mo>&times;</mo>
<mi>d</mi>
<mi>B</mi>
<mi>f</mi>
<mo>+</mo>
<mi>K</mi>
<mo>-</mo>
<mn>4.7607</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>0.35</mn>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>dBf</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mn>4.4</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>dBf</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mn>4.4</mn>
</mrow>
<mo>)</mo>
</mrow>
<mo><</mo>
<mi>d</mi>
<mi>B</mi>
<mi>f</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
2), for log-frequency dBf corresponding to each Frequency point, Selecting All Parameters dBf0With one group of K in the respective region of search
Parameter, σ under current logarithmic frequency is found from formula (5)maxCorresponding function expression, by current logarithmic frequency, the ginseng chosen
Number dBf0The σ found is substituted into K valuesmaxIn function expression, the σ under current logarithmic frequency is obtainedmaxCalculated value, by the calculated value
The σ of the respective frequencies point obtained with step 1maxMeasurement data σmax(dBfn) subtract each other, error amount is obtained, then error amount is asked
Square;All log-frequency points are traveled through, obtain square of error amount corresponding to each Frequency point, finally, to all error amounts
Square summation;
3), using method 2), according to the scouting interval of setting, dBf is constantly changed in the respective region of search0And K, obtain
Corresponding error amount square and value, selection corresponding one group of parameter dBf when making and being worth minimum0With K estimationWith
Step 3, insect Size calculation:
By estimated resultSubstitute into following formula:
<mrow>
<mi>a</mi>
<mo>=</mo>
<mn>70.246</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mi>d</mi>
<mover>
<mi>B</mi>
<mo>^</mo>
</mover>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>/</mo>
<mn>10</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>0.482</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
Obtain the length 2a of insect;
WillIt is updated in following formula:
<mrow>
<mi>V</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>58510</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mover>
<mi>K</mi>
<mo>^</mo>
</mover>
<mo>/</mo>
<mn>10</mn>
</mrow>
</msup>
<mo>-</mo>
<mn>2.89</mn>
<mo>)</mo>
</mrow>
<mrow>
<mn>4</mn>
<mo>/</mo>
<mn>3</mn>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
Obtain insect volume V;
Insect width 2b is calculated according to the cubature formula of ellipsoid simultaneously:
<mrow>
<mi>b</mi>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<mn>3</mn>
<mi>V</mi>
</mrow>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
<mi>a</mi>
</mrow>
</mfrac>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
2. a kind of insect dimension measurement method based on multifrequency scattering modeling as claimed in claim 1, it is characterised in that described
Multiband radar system covering X-band, Ku wave bands and Ka wave bands, each band are not less than 2GHz.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610171604.1A CN105759252B (en) | 2016-03-24 | 2016-03-24 | A kind of insect dimension measurement method based on multifrequency scattering modeling |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610171604.1A CN105759252B (en) | 2016-03-24 | 2016-03-24 | A kind of insect dimension measurement method based on multifrequency scattering modeling |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105759252A CN105759252A (en) | 2016-07-13 |
CN105759252B true CN105759252B (en) | 2018-02-09 |
Family
ID=56346248
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610171604.1A Active CN105759252B (en) | 2016-03-24 | 2016-03-24 | A kind of insect dimension measurement method based on multifrequency scattering modeling |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105759252B (en) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106501802B (en) * | 2016-04-18 | 2019-01-18 | 北京理工大学 | High-resolution multidimensional synergistic insect is migrated Radar Measurement Instrument |
CN106908786B (en) * | 2017-03-09 | 2019-05-03 | 北京理工大学 | A kind of high-precision insect bodies axis based on polarization scattering matrix estimation is towards extracting method |
CN107688169B (en) * | 2017-08-08 | 2020-09-29 | 北京理工大学 | Insect characteristic parameter inversion method based on polarization power matrix determinant |
CN108646240A (en) * | 2018-04-25 | 2018-10-12 | 北京理工大学 | A kind of high-resolution complete polarization entomological radar detection system and its detection method |
CN109164445B (en) * | 2018-07-12 | 2023-09-12 | 北京理工大学 | Frequency modulation step frequency waveform design method and use method of distributed insect radar |
CN109034356B (en) * | 2018-07-23 | 2021-06-01 | 北京理工大学 | Insect density statistical method based on nearest neighbor method correlation and Gaussian beam volume |
CN111007499B (en) * | 2019-11-29 | 2021-09-17 | 北京理工大学 | Insect body width inversion method based on full polarization radar power matrix eigenvalue |
CN111007498B (en) * | 2019-11-29 | 2021-09-17 | 北京理工大学 | Insect body width inversion method based on full polarization radar power matrix determinant |
CN111007500B (en) * | 2019-12-04 | 2021-11-26 | 北京理工大学 | Radon transform-based insect ascending and descending rate rapid measurement method |
CN112327295B (en) * | 2020-10-23 | 2024-05-03 | 北京理工大学 | Method for acquiring and correcting radar extracted insect orientation error |
CN112327294B (en) * | 2020-10-23 | 2023-10-20 | 北京理工大学 | Insect body type parameter inversion method and device based on multi-frequency point polarization characteristics |
CN114396891B (en) * | 2022-02-15 | 2023-09-22 | 集美大学 | Bernoulli adsorption ratchet wheel self-adjustment furniture product contour acquisition device |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102008048719A1 (en) * | 2007-09-21 | 2009-04-02 | Technische Universität Ilmenau | Insects/larvae detecting method for wooden construction of house, involves detecting movement of insects, suppressing interfering signals by using filters, and connecting device with object, where radiation of device penetrates object |
CN104142499A (en) * | 2014-02-26 | 2014-11-12 | 侯加林 | Cryptic insect detecting system based on Doppler effect |
CN104597443A (en) * | 2015-01-22 | 2015-05-06 | 成都锦江电子系统工程有限公司 | Millimeter-wave radar networking based insect detection system |
CN105277929A (en) * | 2015-10-22 | 2016-01-27 | 北京理工大学 | Flapping frequency measurement method based on carrier phase doppler |
-
2016
- 2016-03-24 CN CN201610171604.1A patent/CN105759252B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102008048719A1 (en) * | 2007-09-21 | 2009-04-02 | Technische Universität Ilmenau | Insects/larvae detecting method for wooden construction of house, involves detecting movement of insects, suppressing interfering signals by using filters, and connecting device with object, where radiation of device penetrates object |
CN104142499A (en) * | 2014-02-26 | 2014-11-12 | 侯加林 | Cryptic insect detecting system based on Doppler effect |
CN104597443A (en) * | 2015-01-22 | 2015-05-06 | 成都锦江电子系统工程有限公司 | Millimeter-wave radar networking based insect detection system |
CN105277929A (en) * | 2015-10-22 | 2016-01-27 | 北京理工大学 | Flapping frequency measurement method based on carrier phase doppler |
Non-Patent Citations (2)
Title |
---|
THE FEASIBILITY OF USING VERTICAL-LOOKING RADAR TO MONITOR THE MIGRATION OF BROWN PLANTHOPPER AND OTHER INSECT PESTS OF RICE IN CHINA;Joe R.Riley et al.;《ENTOMOLOGIA SINICA》;20030331;第10卷(第1期);第1-19页 * |
垂直监测昆虫雷达研究进展;张智等;《昆虫学报》;20120731;第55卷(第7期);第849-859页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105759252A (en) | 2016-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105759252B (en) | A kind of insect dimension measurement method based on multifrequency scattering modeling | |
Kim et al. | Experimental analysis of small drone polarimetry based on micro-Doppler signature | |
CN102854504B (en) | Method for sparse synthetic aperture radars imaging on basis of echo simulation operators | |
Ash et al. | Wide-angle synthetic aperture radar imaging: Models and algorithms for anisotropic scattering | |
CN105044693B (en) | Microwave relevance imaging radar amplitude and phase error correction method based on auxiliary array element | |
CN105866770B (en) | A kind of long quality union measuring method of insect bodies based on multiband scatterometry | |
CN105928448B (en) | A kind of insect dimension measurement method that scattering modeling is parsed based on Rayleigh region | |
Harman | Analysis of the radar return of micro-UAVs in flight | |
Yeary et al. | A brief overview of weather radar technologies and instrumentation | |
Park et al. | Simulation and analysis of polarimetric radar signatures of human gaits | |
Glushko et al. | Analysis of differential Doppler velocity for remote sensing of clouds and precipitation with dual-polarization S-band radar | |
Brown | FM airborne passive radar | |
Faia et al. | Detection and location of nonlinear scatterers using DORT applied with pulse inversion | |
McDonald et al. | Statistical and spectral modelling of medium grazing angle coherent radar data | |
Rogers et al. | The synthetic aperture radar trans-ionospheric radio propagation simulator (SAR-TIRPS) | |
CN108490450B (en) | Airborne forward-looking scanning radar imaging method based on broadband autocorrelation fitting | |
Ammar et al. | Performance enhancement of convolution noise jamming against SAR | |
Wu et al. | A mini radar system for low altitude targets detection | |
CN111796277A (en) | Through-wall radar rapid imaging method based on unmanned aerial vehicle platform | |
Sotirelis et al. | Frequency-dependent reflectivity image reconstruction | |
Xie et al. | Back-projection algorithm based on elliptical polar coordinate for low frequency ultra wide band one-stationary bistatic SAR imaging | |
Ma et al. | Bistatic high frequency radar ocean surface cross section for an FMCW source with an antenna on a floating platform | |
Byrd | Simulation of Polarimetric Phased Array Weather Radars | |
Huang et al. | Optimal scattering polarization characteristic for cylinder target in rain at millimeter wave band | |
Xiong et al. | Micro-motion parameters estimation in non-Gaussian noise via mutual correntropy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Hu Cheng Inventor after: Long Teng Inventor after: Liu Changjiang Inventor after: Wang Rui Inventor after: Zeng Tao Inventor before: Zeng Tao Inventor before: Long Teng Inventor before: Liu Changjiang Inventor before: Wang Rui Inventor before: Hu Cheng |
|
COR | Change of bibliographic data | ||
GR01 | Patent grant | ||
GR01 | Patent grant |