CN111007498A - Insect body width inversion method based on full polarization radar power matrix determinant - Google Patents

Insect body width inversion method based on full polarization radar power matrix determinant Download PDF

Info

Publication number
CN111007498A
CN111007498A CN201911203556.XA CN201911203556A CN111007498A CN 111007498 A CN111007498 A CN 111007498A CN 201911203556 A CN201911203556 A CN 201911203556A CN 111007498 A CN111007498 A CN 111007498A
Authority
CN
China
Prior art keywords
insect
body width
determinant
detg
insects
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911203556.XA
Other languages
Chinese (zh)
Other versions
CN111007498B (en
Inventor
胡程
李卫东
王锐
龙腾
曾涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201911203556.XA priority Critical patent/CN111007498B/en
Publication of CN111007498A publication Critical patent/CN111007498A/en
Application granted granted Critical
Publication of CN111007498B publication Critical patent/CN111007498B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
    • G01B7/02—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness
    • 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/024—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using polarisation effects
    • 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
    • G01S7/418—Theoretical aspects

Abstract

The invention discloses an insect body width inversion method based on a complete polarization radar power matrix determinant, which comprises the steps of firstly, measuring different insect samples by utilizing a complete polarization radar to obtain a polarization scattering matrix of each sample; obtaining a Graves power matrix of each insect sample according to the scattering matrix, and then obtaining a determinant value; taking the determinant value of each insect sample as an independent variable, taking the body width of the insect sample as a dependent variable, and obtaining a fitting formula between the determinant value and the body width of the insect through curve fitting; and aiming at the insects to be detected, acquiring determinant values corresponding to the insects to be detected, and respectively and correspondingly acquiring the body widths of the insects to be detected according to a fitting formula. The invention realizes the measurement of the body width parameters of the insects for the first time, obviously improves the measurement capability of the VLR on the biological parameters of the insects, greatly improves the species identification capability of the migratory flying insects, and has important significance for the research of an aerial ecosystem and the control of pests.

Description

Insect body width inversion method based on full polarization radar power matrix determinant
Technical Field
The invention belongs to the technical field of insect radars, and particularly relates to an insect body width inversion method based on a full-polarization radar power matrix determinant.
Background
The migratory flying of insects is an important reason for the outbreak of insect pests in different places and the pandemic of diseases, and seriously threatens the food production safety of China. If the effective monitoring of the migratory flying insects can be realized, a new early accurate early warning system for insect pests is established, necessary decision information and information support can be provided for realizing 'source' management, and the safety of food production in China is fundamentally guaranteed. The radar is one of the most effective means for monitoring the air migration of insects, and since the first special insect radar appears in the sixties of the 20 th century, the migration phenomena of common orientation, aggregation and stratification of insects and the like are discovered in sequence, so that the research of the migration entomology is greatly promoted.
Currently, a mainstream X-band Vertical insect Radar (VLR) adopts a non-coherent rotating linear polarization system, and can extract characteristic parameters such as body axis orientation, weight, and wing vibration frequency of an insect, but cannot realize species identification of the insect. To realize species identification, morphological parameters such as body length and body width of the insect need to be measured. An insect characteristic parameter inversion method based on a polarization power matrix eigenvalue (patent application number is 201710671461.5) and an insect characteristic parameter inversion method based on a polarization power matrix determinant (patent application number is 201710671464.9) respectively provide insect body length parameter inversion methods based on the polarization power matrix eigenvalue and the polarization power matrix determinant, but the insect body length parameter inversion methods cannot be used for insect body width inversion. Currently, no insect body width inversion method exists.
Disclosure of Invention
In view of the above, the invention provides an insect body width inversion method based on a full polarization radar power matrix determinant, which realizes the inversion of the insect body width and is beneficial to distinguishing flying insects, thereby researching an aerial ecosystem and predicting and controlling insect pests.
The insect body width inversion method based on the complete polarization radar power matrix determinant comprises the following steps:
measuring different insect samples by using a full-polarization radar to obtain a polarization scattering matrix S of each sample;
step two, obtaining a Graves power matrix G ═ S of each insect sample according to the scattering matrix SHS；
Step three, obtaining determinant values DetG of Graves power matrixes of all insect samples;
step four, taking the determinant value DetG of each insect sample as an independent variable, taking the body width of the insect sample as a dependent variable, and obtaining a fitting formula between the determinant value DetG and the body width of the insect through curve fitting;
and step five, aiming at the insects to be detected, obtaining determinant values DetG corresponding to the insects to be detected, and respectively and correspondingly obtaining the body widths of the insects to be detected according to the fitting formula of the step four.
Optionally, the step four-way fitting is performed through a curve, and a fitting formula between the insect determinant value DetG and the body width is obtained in a mode that:
and fitting by using a third-order polynomial model and adopting a linear least square method to obtain a fitting formula between the determinant value DetG and the body width of the insect.
Optionally, in step four, the fitting formula of the body width W of the insect to be tested is:
W＝0.052[lg(DetG)]3+1.8763[lg(DetG)]2+22.7013·lg(DetG)+95.3938。
the invention has the following beneficial effects:
the invention relates to an insect body width inversion method based on a complete polarization radar power matrix determinant, which provides an effective means for accurately inverting the body width of an insect from data observed by a complete polarization radar; according to the invention, through the research on a Polarization Scattering Matrix (PSM) measured by a full polarization radar, the strong correlation relationship between the determinant value DetG of the Graves power matrix and the body width of the insect is found, so that a method for inverting the body width of the insect by the determinant value DetG of the full polarization Graves power matrix is provided; the method realizes the measurement of the body width parameters of the insects for the first time, remarkably improves the measurement capability of the insect radar on the biological parameters of the insects, greatly improves the type identification capability of the migratory flying insects, and has important significance for the research of an aerial ecosystem and the treatment of pests.
Drawings
FIG. 1 is a plot of insect power matrix determinant versus insect body width scatter based on fully polarimetric data; the circle is the measured insect data point; the solid line is a third-order polynomial fitting curve of the determinant and the body width of the insect.
FIG. 2 is a flow chart of the method of the present invention.
Fig. 3 is an insect full polarization data measuring device.
FIG. 4 is a plot of insect power matrix determinant versus insect body width scatter based on rotational polarization measurement data.
Detailed Description
The invention is described in detail below by way of example with reference to the accompanying drawings.
The full polarization radar directly measures full polarization echoes of a target, i.e., HH, HV, VH, and VV polarization echoes (H denotes horizontal, V denotes vertical; the first letter denotes the polarization of a received signal, and the second letter denotes the polarization of a transmitted signal) by simultaneously or time-divisionally transmitting the horizontal and vertical polarization signals and simultaneously receiving the horizontal and vertical polarization signals. The 4 polarization echoes are directly corresponding to 4 elements of the target PSM after signal post-processing, thereby obtaining the target PSM. Compared with PSM measured by the rotary linear polarization radar in the patent of an insect characteristic parameter inversion method based on a polarization power matrix eigenvalue (patent application number 201710671461.5) and the patent of an insect characteristic parameter inversion method based on a polarization power matrix determinant (patent application number 201710671464.9), the PSM can be directly measured by the full polarization radar, and the target PSM can be obtained by the rotary linear polarization radar only by performing parameter estimation on an echo signal.
The insect PSM measured at the fully polarized radar is:
Figure BDA0002296446320000041
wherein s is11、s12、s21And s22β, β' and γ are the corresponding relative phases, respectively, as the square root of the target HH, HV, VH and VV polarization RCS (Radar cross-section).
The Graves power matrix of the target is defined as
Figure BDA0002296446320000042
Wherein the superscript H denotes the conjugate transpose.
The determinant value of the power matrix may be expressed as
Figure BDA0002296446320000043
We will next study the relationship between DetG and insect body width based on PSM of 80 insects measured by X-band fully-polarized radar.
First, a scatter plot of the correspondence between lg (DetG) and the body width (unit: mm) of the insect is shown in FIG. 1, in which the circle represents the measured insect data (80) of the full-polarization radar, and the body width, body length, and body weight of the corresponding insect range from 2.3mm to 13.9mm, 10.8mm to 46.5mm, and 25.6mg to 935.3 mg. It can be seen that the DetG has strong correlation with the body width, and a fitting formula between the insect determinant value DetG and the body width can be obtained through curve fitting. This example uses a third order polynomial to fit the relationship between DetG and body width. First, assume that the body width W can be expressed in the form of a third order polynomial with DetG as an argument:
W＝a3[lg(DetG)]3+a2[lg(DetG)]2+a1[lg(DetG)]+a0(4)
in the above formula, the coefficient a of the polynomialiI ═ 1,2,3 are unknown and can be estimated by linear least squares from the measured body widths W and DetG:
A＝(XTX)-1XTW (5)
wherein A is a vector consisting of polynomial coefficients:
A＝[a3,a2,a1,a0]T(6)
x is a matrix consisting of insects DetG:
Figure BDA0002296446320000051
wherein, subscript is insect number, N is 80; w is a vector consisting of the body width of the insect:
W＝[w1w2… wN]T(8)
the fitting result of the polynomial coefficient can be obtained by taking the actually measured data of 80 insects into formula (5), and then the fitting polynomial coefficient into formula (4) can be obtained as the fitting formula of the body width W of the insects:
w (unit: mm) ═ 0.052[ lg (DetG)]3+1.8763[lg(DetG)]2+22.7013·lg(DetG)+95.3938 (9)
The fitted curve is shown as the curve in fig. 1. The correlation coefficient of the insect body width to the fitted curve was 0.89(P < 0.001). The formula (9) is an empirical formula for inverting the body width of the insect.
Therefore, the invention provides an insect body width inversion method based on a complete polarization radar power matrix determinant, and the implementation steps are described in the following specific embodiment:
in order to verify the insect body width inversion method, 80 insects are subjected to body width inversion by the insect body width inversion method based on the power matrix determinant of the fully-polarized radar based on X-band fully-polarized radar actually-measured PSM data, as shown in FIG. 2:
the method comprises the following steps of firstly, measuring insect PSM in a complete polarization mode;
a fully polarized radar system is built based on a 4-port vector network analyzer and a device shown in figure 3, and insect PSM is measured. 2 dual-polarized antennas are arranged at the lower end port of the hollow horn with the inner wall fully adhered with the wave-absorbing material and vertically irradiate upwards; the PE thread is adhered to the back of the insect to fix the insect on the upper port of the loudspeaker, and the head direction of the insect can be adjusted by rotating the metal ring. And 1 of the 2 dual-polarized antennas is used as a transmitting antenna, the other one is used as a receiving antenna, each antenna is provided with a horizontal (H) polarized port and a vertical (H) polarized port, and 4 ports of the vector network analyzer are respectively connected with 4 ports of the two antennas. Full polarization measurement of the target PSM can be achieved by transmitting both H and V polarized signals, while receiving both H and V polarized signals.
Step two, extracting DetG from PSM of insects by using formula (2) -formula (3);
thirdly, inverting the body width of the insect by using a formula (9) based on the characteristic parameter DetG;
the average relative error of the insect body width based on the characteristic parameter DetG inversion is 15.78%. The body width inversion formula given here is a unified formula that can be used directly for all insects.
The method is only suitable for the full polarization radar, and the inversion of the body width of the insect is realized. Fig. 4 shows a scatter diagram of DetG and body width of 168 insect data measured by a rotating polarization radar in a patent "a polar power matrix eigenvalue-based insect characteristic parameter inversion method" (patent application No. 201710671461.5) and a patent "a polar power matrix determinant-based insect characteristic parameter inversion method" (patent application No. 201710671464.9), and it can be seen that although DetG is positively correlated with insect body width, the correlation is weak and insufficient for body width inversion. Therefore, the method is not suitable for a rotary polarization radar.
In summary, the above is only an embodiment of the present invention, and is not intended to limit the scope of the present invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (3)

1. An insect body width inversion method based on a full polarization radar power matrix determinant is characterized by comprising the following steps:
measuring different insect samples by using a full-polarization radar to obtain a polarization scattering matrix S of each sample;
step two, obtaining Graves work of each insect sample according to the scattering matrix SRate matrix G ═ SHS；
Step three, obtaining determinant values DetG of Graves power matrixes of all insect samples;
step four, taking the determinant value DetG of each insect sample as an independent variable, taking the body width of the insect sample as a dependent variable, and obtaining a fitting formula between the determinant value DetG and the body width of the insect through curve fitting;
and step five, aiming at the insects to be detected, obtaining determinant values DetG corresponding to the insects to be detected, and respectively and correspondingly obtaining the body widths of the insects to be detected according to the fitting formula of the step four.
2. The insect body width inversion method based on the full polarization radar power matrix determinant as claimed in claim 1, wherein the step four-way fitting is through a curve fitting, and a fitting formula between the insect determinant value DetG and the body width is obtained by:
and fitting by using a third-order polynomial model and adopting a linear least square method to obtain a fitting formula between the determinant value DetG and the body width of the insect.
3. The insect body width inversion method based on the complete polarization radar power matrix determinant as claimed in claim 1,
in the fourth step, the fitting formula of the body width W of the insect to be detected is as follows:
W＝0.052[lg(DetG)]3+1.8763[lg(DetG)]2+22.7013·lg(DetG)+95.3938。
CN201911203556.XA 2019-11-29 2019-11-29 Insect body width inversion method based on full polarization radar power matrix determinant Active CN111007498B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911203556.XA CN111007498B (en) 2019-11-29 2019-11-29 Insect body width inversion method based on full polarization radar power matrix determinant

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911203556.XA CN111007498B (en) 2019-11-29 2019-11-29 Insect body width inversion method based on full polarization radar power matrix determinant

Publications (2)

Publication Number Publication Date
CN111007498A true CN111007498A (en) 2020-04-14
CN111007498B CN111007498B (en) 2021-09-17

Family

ID=70113654

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911203556.XA Active CN111007498B (en) 2019-11-29 2019-11-29 Insect body width inversion method based on full polarization radar power matrix determinant

Country Status (1)

Country Link
CN (1) CN111007498B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016022637A1 (en) * 2014-08-08 2016-02-11 Urthecast Corp. Apparatus and methods for quad-polarized synthetic aperture radar
CN105759252A (en) * 2016-03-24 2016-07-13 北京理工大学 Insect dimension measurement method based on multi-frequency scattering modeling
CN105928448A (en) * 2016-04-18 2016-09-07 北京理工大学 Insect dimension measurement method based on Rayleigh region analysis scattering modeling
CN107589412A (en) * 2017-08-08 2018-01-16 北京理工大学 A kind of insect characteristic parameter inversion method based on polar power matrix exgenvalue
CN107688169A (en) * 2017-08-08 2018-02-13 北京理工大学 A kind of insect characteristic parameter inversion method based on polar power matrix determinant

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016022637A1 (en) * 2014-08-08 2016-02-11 Urthecast Corp. Apparatus and methods for quad-polarized synthetic aperture radar
CN105759252A (en) * 2016-03-24 2016-07-13 北京理工大学 Insect dimension measurement method based on multi-frequency scattering modeling
CN105928448A (en) * 2016-04-18 2016-09-07 北京理工大学 Insect dimension measurement method based on Rayleigh region analysis scattering modeling
CN107589412A (en) * 2017-08-08 2018-01-16 北京理工大学 A kind of insect characteristic parameter inversion method based on polar power matrix exgenvalue
CN107688169A (en) * 2017-08-08 2018-02-13 北京理工大学 A kind of insect characteristic parameter inversion method based on polar power matrix determinant

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHENG HU等: "Insect Biological Parameter Estimation Based on the Invariant Target Parameters of the Scattering Matrix", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 *
刘召庆等: "一种介质体宽角域散射的快速分析方法", 《火控雷达技术》 *

Also Published As

Publication number Publication date
CN111007498B (en) 2021-09-17

Similar Documents

Publication Publication Date Title
CN107589412B (en) Insect characteristic parameter inversion method based on polarization power matrix characteristic value
CN107688169B (en) Insect characteristic parameter inversion method based on polarization power matrix determinant
Hernández-Mesa et al. Current applications and perspectives of ion mobility spectrometry to answer chemical food safety issues
Mirkovic et al. Electromagnetic model reliably predicts radar scattering characteristics of airborne organisms
CN106443611B (en) RCS (radar cross section) measuring method of weak scattering target
CN106646344B (en) A kind of Wave arrival direction estimating method using relatively prime battle array
Zito et al. A 90nm CMOS SoC UWB pulse radar for respiratory rate monitoring
CN105866770B (en) A kind of long quality union measuring method of insect bodies based on multiband scatterometry
CN108957415A (en) The real-time Calibration Method and system of Analysis of Differential Reflectivity Factor Measured
CN106908786B (en) A kind of high-precision insect bodies axis based on polarization scattering matrix estimation is towards extracting method
CN111007498B (en) Insect body width inversion method based on full polarization radar power matrix determinant
CN109696711A (en) A kind of convolutional neural networks target detection recognition methods based on polarization scattering characteristics
CN111007499B (en) Insect body width inversion method based on full polarization radar power matrix eigenvalue
CN103955008B (en) A kind of amplitude calibration method for Multi probe near-field scattering imaging
Chen et al. Time-domain ultra-wideband microwave imaging radar system
CN111007497B (en) Parallel and vertical insect distinguishing method based on characteristic phase
Naghibi et al. Skin artifact removal algorithms for radar-based microwave breast cancer detection
CN109212483A (en) It is a kind of to solve 90 degree of the big figure insect body axis towards mistake towards extracting method
CN103983495A (en) Tetrodotoxin positive quality control sample as well as preparation method and application thereof
CN106646458B (en) Two-dimentional through-wall detection target radar system
de Moustier OS-CFAR detection of targets in the water column and on the seafloor with a multibeam echosounder
CN104320205B (en) Sparse DOA algorithm for estimating in Spatial Doppler domain
Yang et al. Microwave imaging for breast cancer detection using Vivaldi antenna array
Chufeng et al. Signal processing techniques in high-resolution RCS measurement system
McKenzie et al. Identifying Little Sandy Desert bat species from their echolocation calls.

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant