NL2026449A - Abnormal grain condition detection method based on radio tomographic imaging - Google Patents

Abnormal grain condition detection method based on radio tomographic imaging Download PDF

Info

Publication number
NL2026449A
NL2026449A NL2026449A NL2026449A NL2026449A NL 2026449 A NL2026449 A NL 2026449A NL 2026449 A NL2026449 A NL 2026449A NL 2026449 A NL2026449 A NL 2026449A NL 2026449 A NL2026449 A NL 2026449A
Authority
NL
Netherlands
Prior art keywords
abnormal grain
grain condition
formula
follows
condition
Prior art date
Application number
NL2026449A
Other languages
Dutch (nl)
Other versions
NL2026449B1 (en
Inventor
Zhu Chunhua
Tian Jiake
Deng Miaolei
Yang Tiejun
Yang Weidong
Guo Xinying
Yang Jing
Fu Hongliang
Wang Jiaojiao
Original Assignee
Univ Henan Technology
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 Univ Henan Technology filed Critical Univ Henan Technology
Publication of NL2026449A publication Critical patent/NL2026449A/en
Application granted granted Critical
Publication of NL2026449B1 publication Critical patent/NL2026449B1/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N22/00Investigating or analysing materials by the use of microwaves or radio waves, i.e. electromagnetic waves with a wavelength of one millimetre or more
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/02Food
    • G01N33/10Starch-containing substances, e.g. dough
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Business, Economics & Management (AREA)
  • Primary Health Care (AREA)
  • General Business, Economics & Management (AREA)
  • Mining & Mineral Resources (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Animal Husbandry (AREA)
  • Economics (AREA)
  • Agronomy & Crop Science (AREA)
  • Theoretical Computer Science (AREA)
  • Electromagnetism (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Arrangements For Transmission Of Measured Signals (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

An abnormal grain condition detection method based on radio tomographic imaging includes the following sequential steps: (1) arranging sensors: arranging a plurality of sensor nodes around a grain depot; (2) establishing a non-abnormal grain condition 5 database: transmitting radio signals, and summarizing radio signals received by the sensors to establish the non-abnormal grain condition database; (3) collecting experimental data of an abnormal grain condition: transmitting radio signals, and summarizing radio signals received by the sensors to collect the experimental data of the abnormal grain condition; (4) preprocessing the collected experimental data; (5) 10 establishing a linear model of communication link attenuation on the pre-processed data; (6) solving a vector matrix; (7) solving the linear model of communication link attenuation to obtain an image attenuation value; and (8) performing image reconstruction to find out a position of the abnormal grain condition. The method of the present invention improves the accuracy of the abnormal grain condition detection; 15 which is simple and easy to implement.

Description

-1-
ABNORMAL GRAIN CONDITION DETECTION METHOD BASED ON RADIO TOMOGRAPHIC IMAGING
TECHNICAL FIELD The present invention belongs to the technical field of abnormal grain condition detection, and more particularly, relates to an abnormal grain condition detection method based on radio tomographic imaging.
BACKGROUND Grain security is always a serious problem worldwide. Grain condition detection can find the abnormal grain condition in the stored grain heap, so as to make early warning and policy regulation for grain security ahead of time. At present, the abnormal grain condition detection means mainly adopts manual detection to the grain depot, depot penetrating radar and wired sensor network detection. Manual detection to the grain depot requires lots of manpower, material resources and financial resources, which is low in the speed and accuracy, and it is difficult to find the abnormal grain condition in the grain heap. With respect to wired sensor network detection, various sensors are embedded into and equipped into the grain heap to collect various grain condition monitoring data in the coverage by the sensor nodes in the grain depot in real time.
However, existing wired sensor network detection has problems such as difficult layout, single monitoring parameter, small coverage, difficult upgrading and others. Conventional ground penetrating radar is called as “depot penetrating radar” when used for the abnormal grain condition detection. When detecting the grain depot, the receiving and transmitting antenna is arranged on the grain surface, the abnormal condition is estimated by microwave reflection theory and electromagnetic wave reflection imaging caused by dielectric characteristic change interface of the grain heap. Depot penetrating radar based on reflection theory faces two problems in practical applications. The first problem is the contradiction between the detection depth and the resolution of detection target. Specifically, the imaging resolution of detection target is inversely proportional to the detection frequency, and the low frequency of the radar seriously affects the accuracy of anomaly target recognition. The second problem is the imaging difficulty caused by small difference of relative dielectric constants between the grain and the abnormal grain condition. Since the magnitude of the detected reflection wave of the
-2- anomaly target is mainly determined by difference of relative dielectric constants between the grain and the abnormal grain condition at the interface, when the difference between the dielectric constant of media of the grain and the dielectric constant of media of the abnormal grain condition is not large enough, the reflection wave is weak, while the transmission wave is strong, so that it is difficult to obtain the return wave signal to cause the difficulty of imaging. Therefore, from the aspect of detection resolution, the microwave detection based on transmission theory presents obvious advantages. Recently, a radio tomographic imaging technique based on transmission theory has emerged. The radio tomographic imaging uses distributed wireless sensor nodes to measure power attenuation of the link and obtain the signal attenuation in region by inversion. Radio tomographic imaging has low hardware requirements, and the radio signals have good penetration. Radio tomographic imaging has been widely used in target positioning and tracking, vehicle recognition, through-the-wall detection, medical sanitation and other fields. Compared with conventional ground penetrating radar, the radio tomographic imaging has great advantages. Firstly, radio tomographic imaging can work at a relative high frequency band to obtain a better spatial resolution. Secondly, radio tomographic imaging only requires to measure power of radio signals, and has a low requirement for nodes, which can reduce cost of the system and facilitate miniaturization of the system.
SUMMARY An objective of the present invention is to provide an abnormal grain condition detection method based on radio tomographic imaging with a high detection accuracy.
In order to solve the technical problems mentioned above, the present invention provides the following technical solution. An abnormal grain condition detection method based on radio tomographic imaging includes the following sequential steps: (1) arranging sensors: arranging a plurality of sensor nodes around a grain depot; (2) establishing a non-abnormal grain condition database: transmitting radio signals, and summarizing radio signals received by the sensors to establish the non-abnormal grain condition database;
-3-
(3) collecting experimental data of an abnormal grain condition: transmitting radio signals, and summarizing radio signals received by the sensors to collect the experimental data of the abnormal grain condition;
(4) preprocessing the collected experimental data according to formula (1);
ite NN f(jjx h4 => SA, (1), {0 where, f(i) =exp((i-r) /(2*6°)/N2*S° is a Gaussian mask; Af, An, Ar are signal strength values received by the sensors on an /™ link at (1-0)! moment, 2 moment, (7+ )'® moment, respectively; and @ is a filtering window;
(5) establishing a linear model of communication link attenuation on the pre-
processed data to obtain formula (2);
Ar, = WAx, +n, (2),
where, Ar, represents a measurement matrix consisting of received signals of all communication links at moment; Ax, represents a shadow fading vector of a received signal strength of all grids at f moment; n, represents a noise vector, and W is a vector matrix of a weight of each grid on each link; (6) solving a vector matrix by formula (3); ch í e” d()+d,(2)<d, +4 0 otherwise ’
where, d, is alength of the / link; d, (1) and d,(2) are distances from aj" grid to two sensors, respectively; A is a parameter adjusting an range of an ellipse; and 4 is a distance from a grid inside the elliptic to a current link;
(7) solving the linear model of communication link attenuation by formula (4) to obtain an image attenuation value;
Ak, = (WIW +al' Ty W'Ar, (4),
where, Ax, is a derivative of Ax,, and WTisa transpose of W; and
(8) performing image reconstruction to find out a position of the abnormal grain condition;
(1) comparing each image attenuation value to obtain a grid coordinate Jmax of the maximum image attenuation value;
-4- (11) obtaining the position of the abnormal grain condition by formula (5); P= (5), where, P is the position of the abnormal grain condition. In step (4), the filtering window takes a value of 5. In step (5), a method for establishing the linear model of communication link attenuation on the pre-processed data is as follows: the received signal strength measurement #, of the / link at t moment is expressed as follows: hy = P, 4 Ss, -h, 4%, (6), where, ©, represents a transmitting power of a sensor node, £; represents a path loss of the /? link, S,, represents a shadow fading caused by shielding the link through the abnormal grain condition, /;, represents a multi-path fading noise, and 7%, represents a measured noise of a monitoring area; a received signal strength change caused by the abnormal grain condition is follows: A, =h-h = Si +, +n, (7), the monitoring area of the grain depot is divided evenly into N grids; the shadow fading of the ? link caused by the abnormal grain condition is expressed by a weighted sum of all grid values; and the shadow fading of the link is expressed by using a spatial integral as follows:
N S= D Wij A, (8), j=1 where, Ax ; represents a signal strength change of the monitoring area of the grain depot at the /!! grid, and wy, represents a weight of the / grid on the / link; since a medium of the monitoring area of the grain depot is single, is similar to an outdoor environment and has a small multi-path fading, +, 30, the received signal strength change of the /! link in the monitoring area is as follows:
N Ar, = Dw, Ax a HR, (9), J=1 formula (9) only considers the received signal strength change of the / link, and when all links are considered, formula (2) is obtained as follows:
-5- Ar =WAx +n, (2), where, Ar, represents the measurement matrix consisting of the received signals of all communication links at ‘moment; Ax, represents the shadow fading vector of the received signal strength of all grids at # moment; n, represents the noise vector; and W is the vector matrix of the weight of each grid on each link.
In step (6), a method for obtaining the formula (3) is as follows: modeling an area influenced by the abnormal grain condition as the ellipse, wherein a focus of the ellipse is a receiving and transmitting node; since grids inside the ellipse are considered to cause influence on a result, giving a weight to the grids inside the ellipse; while since grids outside the ellipse are considered to cause no influence on the result, giving a weight of 0 to the grids outside the ellipse.
In step (7), a method for obtaining the formula (4) is as follows: firstly, minimizing an objective function as follows: min Ar, - WA, | + a UAx | (10), 4x, where, T is a Tikhonov matrix and represents a prior information of a model solution; |TAx,| represents a penalty term for Tikhonov regularization; and @ is an adjustable regularization parameter, a value of & determines whether a final solution is biased to measured data or biased to the prior information; then, an approximate estimation of AX, is obtained by taking a derivative of formula (9) to obtain formula (4) as follows: Ak, = (WW +al"T)" WT Ar, (4) By adopting the above technical solution, the present invention has the following advantages. Compared with traditional methods, simulation experiment shows that the method of the present invention has more advantages. In order to solve those difficult problems such as large workload of grain detection, small difference of dielectric constants of abnormal grain condition and low accuracy, a new detection method is proposed by the present invention to provide new technical solution and research direction for the abnormal grain condition detection. This new method can effectively find out the position of the abnormal grain condition, which has a high accuracy and a fast detection speed.
-6-
BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a flow chart of the present invention; FIG. 2 shows an experimental result when water anomaly exists in the grain condition in case of four sensor nodes; FIG. 3 shows an experimental result when water anomaly exists in the grain condition in case of six sensor nodes; FIG. 4 shows an experimental result when water anomaly exists in the grain condition in case of eight sensor nodes; FIG. 5 is a comparative chart of median errors when air anomaly, metal anomaly and plastic anomaly exist in the grain condition; FIG. 6 shows an experimental result when air anomaly exists in the grain condition; FIG. 7 shows an experimental result when metal anomaly exists in the grain condition; FIG. 8 shows an experimental result when plastic anomaly exists in the grain condition; FIG. 9 is a schematic diagram showing positions of the abnormal grain condition; FIG. 10 shows median errors of different positions; and FIG. 11 is a comparative chart of median errors of the abnormal grain condition with image preprocessing and without image preprocessing.
DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS As shown in FIG. 1, an abnormal grain condition detection method based on radio tomographic imaging includes the following sequential steps. (1) Sensors are arranged. A plurality of sensor nodes are arranged around a grain depot. In the present embodiment, a resin-based transparent plastic box with a size of
-7- 200 mm * 300 m * 500 mm is adopted to simulate the grain depot. four sensor nodes, six sensor nodes, and eight sensor nodes are successively arranged to perform grouped experiments, and the sensor nodes are required to be located at a same horizontal plane.
(2) A non-abnormal grain condition database is established. Radio signals are transmitted, and radio signals received by the sensors are summarized to establish the non-abnormal grain condition database. The non-abnormal grain is poured into the grain depot. Then, a sensor network transmits and receives the radio signals automatically for ten minutes, summarizes the received signal strength data, extracts the received signal strength feature of each sensor node, and establishes the non-abnormal grain condition database.
(3) Experimental data of an abnormal grain condition are collected. Radio signals are transmitted, and radio signals received by the sensors are summarized to collect the experimental data of the abnormal grain condition. The abnormal grain is poured into the grain depot. Then, the sensor network continuously receives and transmits the radio frequency signals, and extracts the received signal strength of each sensor node periodically as the experimental data of the abnormal grain condition.
(4) Data collected in step (3) are processed. Grain belongs to granular body. When the radio signal is transmitted in the grain depot through transmission, the grain surface easily reflects the radio signal. Since the slight change of the received signals caused by the grain reflection belongs to fast fading, it is usually modeled as the noise. When the system processes the signals in time domain, a low-pass filter is used to reduce the influence of the fast fading. The low-pass filter employs a one-dimensional Gaussian filter. Assuming that the size of the Gaussian filtering window is © and the standard deviation of the Gaussian function is & , the output signal can be expressed as follows: +o ne=2 SOA, to where, f(D) =exp((i-r) /(Q*6°)/ 2s? ‚which is a Gaussian mask; Af, Ar, Ar. are the signal strength values received by the sensors on the /" link at (t-0)% moment, i moment, (7+ )™ moment, respectively; and © is the filtering window. In the present embodiment, © takes a value of 5. Since the excessive large value of the filtering window © increases the time delay and the excessive small value of the filtering window © causes the poor denoising effect, after weighing the time
-8- delay and the validity, the present invention takes a value of 5 as the size of the filtering window.
In this step, the data is pre-processed to improve the accuracy of the detection result.
From the above, grain is granular body.
If the data is not pre-processed, the existing excessive large noise decreases the accuracy greatly, while after the data being pre- processed, a large part of noise is removed, and the accuracy is improved significantly and is far higher than the accuracy obtained in case that the data is not pre-processed.
As shown in FIG. 11, after pre-processing, the simulation experiment is closer to the actual location and has a relatively small median error. (5) A linear model of communication link attenuation is established on the pre- processed data.
When there is no abnormal grain condition in the grain depot, the received signal strength measurement of the /" link is recorded as # . When an abnormal grain condition exists in the grain depot, the abnormal grain condition reflects and shields the radio signals to cause shadow fading, so that the received signal strength measurement of the link changes.
The received signal strength measurement 4%, of the / link at t moment is expressed as follows: he = P -f Ss, Pi, TM, (6), where, P, represents a transmitting power of a sensor node, I} represents a path loss of the MM link, 5;, represents the shadow fading caused by shielding the link through the abnormal grain condition, +, represents the multi-path fading noise, and #;, represents the measured noise of the monitoring area.
The received signal strength change caused by the abnormal grain condition is follows: Ar, hhh; Tr 8 +k, +n, (7). The monitoring area of the grain depot is divided evenly into N grids.
The shadow fading of the /* link caused by the abnormal grain condition can be expressed by a weighted sum of all grid values . The shadow fading of the /M link is expressed by using a spatial integral as follows:
-9-
N Si = 2 wax, (8),
JT where, AX, represents the signal strength change of the monitoring area of the grain depot at aj" grid, and wy, represents a weight of the / grid on the / link. Since the medium of the monitoring area of the grain depot is single, is similar to the outdoor environment and has a small multi-path fading, assuming that £,, 0 the received signal strength change of the / link in the monitoring area is as follows:
N Ar, = 2 yA, +1, (9). I= The above formula (9) only considers the received signal strength (RSS) change of the /* link, and the monitoring area of the grain depot includes L links totally, so that a system of linear equations containing L RSS changes is obtained and expressed in a matrix form as follows: Ar, = WAx, +n, (2), where, Ar, represents a measurement matrix consisting of the received signals of all communication links at # moment; Ax, represents a shadow fading vector of the received signal strength of all grids at f moment; n, represents a noise vector; and W 1s a vector matrix of a weight of each grid on each link.
(6) The vector matrix is solved by formula (3) as follows: Je” d,()+dy(2)<d, +4 hi | 0 otherwise (3), where, d, is a length of the / link; d,(1) and d,(2) are distances from the j* grid to two sensors, respectively; A is a parameter adjusting an range of an ellipse; and / is a distance from a grid inside the elliptic to a current link.
Formula (3) is established depending on the following basis: the area influenced by the abnormal grain condition is modeled as the ellipse whose focus is the receiving and transmitting node. Grids inside the ellipse are considered to cause influence on the result and thus are given a weight, while grids outside the ellipse are considered to cause no influence on the result and thus have a weight of 0.
-10- (7) The linear model of communication link attenuation is solved by formula (4) to obtain an image attenuation value as follows: Ak, = (WIW +a") WTAr, (4), where, AX, is a derivative of Ax,, and WT is a transpose of W.
After the weight matrix W of the ellipse is determined, the formula (2) is solved to obtain Ax, so that image reconstruction is performed to find out the position of the abnormal grain condition.
When solving Ax,, since the number of unknowns in the system of equations is larger than the number of equations, regularization method is required. Tikhonov regularization is an effective means to solve the ill-conditioned equations. Firstly, the objective function is minimized as follows: min Jar, - Wax, [+ ocfrax (10), where, F is a Tikhonov matrix, which represents the prior information of the model solution; [Ax | represents a penalty term for Tikhonov regularization; and « is an adjustable regularization parameter whose value determines whether the final solution is biased to the measured data or biased to the prior information. Then, an approximate estimation of AX, can be obtained by taking a derivative of the minimized objective function as follows: Ak, = (WIW +a Ty WTAr, (4) where, AX, is derivative of Ax;, and WT is the transpose of W. (8) Image reconstruction is performed to find out the position of the abnormal grain condition. (1) Each image attenuation value is compared to obtain the grid coordinate Jmax of the maximum image attenuation value. (11) The position of the abnormal grain condition is obtained by formula (5) as follows: P=Limax (5), where, P is the position of the abnormal grain condition. Experimental analysis:
-11- Simulation analysis is performed on the method of the present invention. Firstly, water anomaly experiments are successively performed on four nodes, six nodes and eight nodes, and the experimental results are shown in FIGS. 2-4. The experimental results show that the experimental median error of four nodes 1s 0.0162m, the experimental median error of six nodes is 0.0135m, and the experimental median error of eight nodes is 0.0128m. Compared with four nodes, the size of imaging area of foreign body of six nodes and eight nodes is more condensed and has a smaller range. Meanwhile, considering the complexity of the equipment, data processing operation and other factors, it is considered that six nodes are suitable for this experiment.
Then, experiments are performed to explore the influence of different grain conditions on experimental results. Experiments are successively performed on void anomaly, metal anomaly and plastic anomaly, and the experimental results are shown in FIGS. 6-8. The experimental results show that when the abnormal grain condition is the metal, in the radio tomographic imaging results, the surface is relatively smooth, and the position of the foreign body even has a basic shape only because the large dielectric constant of the metal absorbs the radio frequency signals more obviously, and the changes in communication link attenuation becomes more obvious.
Finally, experiments are performed to explore the radio tomographic imaging effects of the metal anomaly at different positions in case of six nodes to further research the influence of the position on the imaging results. The positions of the metal during experiment are shown in FIG. 9, and the experimental median errors are shown in FIG.
10. The experimental results show that the median error of the position 3 is minimal while the median error of other positions is larger than that of the position 3 because the position 3 has a larger number of links than other positions and has a more accurate imaging result.
The method of the present invention can detect the abnormal grain condition, which has a high-accuracy detection result and a fast detection speed. Compared with traditional detection methods, the method of the present invention is easier to achieve online measurement, and can fast find the position of the abnormal grain condition without moving grain.

Claims (5)

-12 - NL2026449 Conclusies-12 - NL2026449 Conclusions 1. Detectiewerkwijze voor een abnormale graanomstandigheid op basis van radiotomografische beeldvorming, waarbij de detectiewerkwijze de volgende opeenvolgende stappen omvat: (I) het opstellen van sensoren: het opstellen van een veelvoud aan sensorknopen rond een graandepot; (2) het vaststellen van een database van niet-abnormale graanomstandigheid: het plaatsen van graan met een niet-abnormale graanomstandigheid in het gaandepot, het uitzenden van radiosignalen en het samenvatten van radiosignalen die door de sensoren ontvangen zijn om de database van niet-abnormale graanomstandigheid vast te stellen; (3) het verzamelen van experimentele gegevens van een abnormale graanomstandigheid: het plaatsen van graan in het graandepot, het uitzenden van radiosignalen en het samenvatten van radiosignalen die door de sensoren ontvangen zijn om de experimentele gegevens van de abnormale graanomstandigheid te verzamelen; (4) het voorbewerken van verzamelde experimentele gegevens volgens formule (1); fy = Sr “AR, 1. eo waar, fli) = exp(-(i-r)? / (2*8%)) / N(2* 6°) een Gaussiaans masker is; Ar rm, Ary, Art signaalsterktewaarden zijn die ontvangen worden door de sensoren op een 1° verbinding op respectievelijk (#-»)® moment, i© moment, (#+®)° moment; en @ een filterraam is; (5) het vaststellen van een lineair model van communicatieverbindingverzwakking op de voorbewerkte gegevens om formule (2) te verkrijgen; Art = WAX + nt (2), waar Art een metingenmatrix voorstelt die bestaat uit ontvangen signalen van alle communicatieverbindingen op t moment; Ax¢ een schaduwvervagingsvector van een ontvangen signaalsterkte van alle roosters op f moment voorstelt, nt eenA detection method for an abnormal grain condition based on radiotomographic imaging, the detection method comprising the following sequential steps: (I) arranging sensors: arranging a plurality of sensor nodes around a grain depot; (2) establishing a non-abnormal grain condition database: placing grain with a non-abnormal grain condition in the gantry, transmitting radio signals and summarizing radio signals received by the sensors to create the database of non-abnormal grain condition determine grain condition; (3) collecting experimental data of an abnormal grain condition: placing grain in the grain depot, transmitting radio signals and summarizing radio signals received by the sensors to collect the experimental data of the abnormal grain condition; (4) preprocessing collected experimental data according to formula (1); fy = Sr “AR, 1. eo where, fli) = exp(-(i-r)? / (2*8%)) / N(2* 6°) is a Gaussian mask; Ar rm, Ary, Art are signal strength values received by the sensors on a 1° connection at (#-»)® moment, 1© moment, (#+®)° moment, respectively; and @ is a filter frame; (5) establishing a linear model of communication link attenuation on the preprocessed data to obtain formula (2); Art = WAX + nt(2), where Art represents a measurement matrix consisting of received signals from all communication links at t time; Ax¢ represents a shadow blur vector of a received signal strength of all grids at f moment, nt a -13 - NL2026449 achtergrondvector voorstelt; en W een vectormatrix is van een gewicht van elk rooster op elke verbinding; (6) het oplossen van een vectormatrix middels formule (3); nn jet died 2 ad, +4 i. ie anders or waar d; een lengte van de /° verbinding is; d;(1) en d;(2) respectievelijk afstanden zijn van een j° rooster tot twee sensoren; 4 een parameter is die een bereik van een ellips aanpast; en / een afstand is van een rooster binnen de ellipsvormige tot een huidige verbinding; (7) het oplossen van het lineaire model van communicatieverbindingverzwakking middels formule (4) om een beeldverzwakkingswaarde te verkrijgen; Ak, = (WW an War, ES waar “een afgeleide is van Ax, en WT een transpositie is van W; en (8) het uitvoeren van beeldreconstructie om de positie van de abnormale graanomstandigheid uit te vinden; (1) het vergelijken van elke beeldverzwakkingswaarde om een roostercoördinaat Jmax te verkrijgen van een maximale beeldverzwakkingswaarde; (ii) _ het verkrijgen van de positie van de abnormale graanomstandigheid middels formule (5); P=Ilmax (5), waar P de positie van de abnormale graanomstandigheid is.-13 - NL2026449 represents background vector; and W is a vector matrix of a weight of each grid on each link; (6) solving a vector matrix by formula (3); nn jet died 2 ad, +4 i. ie else or true d; is a length of the /° link; d;(1) and d;(2) are distances from a j° grid to two sensors, respectively; 4 is a parameter that adjusts a range of an ellipse; and / is a distance from a grid within the ellipsoidal to a current connection; (7) solving the linear model of communication link attenuation by formula (4) to obtain an image attenuation value; Ak, = (WW an War, ES where “is a derivative of Ax, and WT is a transposition of W; and (8) performing image reconstruction to find out the position of the abnormal grain condition; (1) comparing each image attenuation value to obtain a grid coordinate Jmax of a maximum image attenuation value, (ii) obtaining the position of the abnormal grain condition by formula (5): P=IImax (5), where P is the position of the abnormal grain condition. 2. Detectiewerkwijze voor een abnormale graanomstandigheid op basis van radiotomografische beeldvorming volgens conclusie 1, met het kenmerk, dat in stap (4) het filterraam een waarde van 5 aanneemt.A grain abnormality detection method based on radiotomographic imaging according to claim 1, characterized in that in step (4) the filter window assumes a value of 5. 3. Detectiewerkwijze voor een abnormale graanomstandigheid op basis van radiotomografische beeldvorming volgens conclusie 1, met het kenmerk, dat in stap (5) een werkwijze voor het vaststellen van het lineaire model van communicatieverbindingverzwakking op de voorbewerkte gegevens als volgt is:A grain abnormality detection method based on radiotomographic imaging according to claim 1, characterized in that in step (5), a method for determining the linear model of communication link attenuation on the preprocessed data is as follows: - 14 - NL2026449 een ontvangen signaalsterktemeting r;, van de /° verbinding op moment wordt als volgt uitgedrukt: Fir = Ps Pi Ste Fit Mut (6), waar P een uitzendvermogen van een sensorknoop voorstelt, P; een padverlies van de FF verbinding voorstelt, Si; een schaduwvervaging voorstelt die veroorzaakt wordt door het afschermen van de verbinding door de abnormale graanomstandigheid, #;; een multipad vervagingsachtergrond voorstelt en +1: een gemeten achtergrond van een geobserveerd gebied voorstelt; een verandering in ontvangen signaalsterkte die veroorzaakt wordt door de abnormale graanomstandigheid is als volgt: A, = 5 “= Se +F + 4, {73 het geobserveerde gebied van het graandepot wordt gelijkmatig in N roosters verdeeld; de schaduwvervaging van de “verbinding die veroorzaakt wordt door de abnormale graanomstandigheid wordt uitgedrukt door een gewogen som van alle roosterwaarden; en de schaduwvervaging van de A verbinding wordt uitgedrukt met behulp van een ruimtelijke integraal als volgt: £ Sie = DM, (8), waar Ax;, een verandering in signaalsterkte van het observeringsgebied van het graandepot bij het /© rooster voorstelt, en w;, een gewicht van het j° rooster op de verbinding voorstelt, en F; =0; de verandering in ontvangen signaalsterkte van de A verbinding in het observeringsgebied is als volgt: x An, = GA, + {8 = formule (9) beschouwt slechts de verandering in ontvangen signaalsterkte van de Ff verbinding, en wanneer alle verbindingen beschouwd zijn, wordt formule (2) als volgt verkregen: Art = WAX: + Nt (2), waar Art de metingenmatrix voorstelt die bestaat uit ontvangen signalen van alle communicatieverbindingen op # moment; Ax: een schaduwvervagingsvector van de ontvangen signaalsterkte van alle roosters op + moment voorstelt, mt dea received signal strength measurement r i , of the /° connection at the moment is expressed as follows: Fir = Ps Pi Ste Fit Mut (6), where P represents a transmit power of a sensor node, P; represents a path loss of the FF compound, Si; represents a shadow blur caused by the shielding of the joint by the abnormal grain condition, #;; represents a multipath blur background and +1 represents a measured background of an observed area; a change in received signal strength caused by the abnormal grain condition is as follows: A, = 5 “ = Se +F + 4, {73 the observed area of the grain depot is evenly divided into N grids; the shadow blur of the “connection caused by the abnormal grain condition is expressed by a weighted sum of all grid values; and the shadow blur of the A compound is expressed using a spatial integral as follows: Sie = DM, (8), where Ax;, represents a change in signal strength of the grain deposit observation area at the /© grid, and w ;, represents a weight of the j° grid on the joint, and F; =0; the change in received signal strength of the A link in the observation area is as follows: x An, = GA, + {8 = formula (9) considers only the change in received signal strength of the Ff link, and when all links are considered, it becomes formula (2) obtained as follows: Art = WAX: + Nt (2), where Art represents the measurement matrix consisting of received signals from all communication links at # moment; Ax: represents a shadow blur vector of the received signal strength of all grids at + moment, mt de -15- NL2026449 achtergrondvector voorstelt; en W de vectormatrix is van het gewicht van elk rooster op elke verbinding.-15- NL2026449 represents background vector; and W is the vector matrix of the weight of each grid on each connection. 4. Detectiewerkwijze voor een abnormale graanomstandigheid op basis van radiotomografische beeldvorming volgens conclusie 1, met het kenmerk, dat in stap (6) een werkwijze voor het verkrijgen van de formule (3) als volgt is: het modelleren van een gebied dat door de abnormale graanomstandigheid beïnvloed wordt als de ellips, waarbij een focus van de ellips een ontvangende en uitzendende knoop is; aangezien roosters binnen de ellips verondersteld worden invloed op een resultaat te hebben, het geven van een gewicht aan de roosters binnen de ellips; terwijl, aangezien de roosters buiten de ellips verondersteld worden geen invloed op het resultaat te hebben, het geven van een gewicht van 0 aan de roosters buiten de ellips.The abnormal grain condition detection method based on radiotomographic imaging according to claim 1, characterized in that in step (6), a method for obtaining the formula (3) is as follows: modeling a region covered by the abnormal grain condition is affected as the ellipse, a focus of the ellipse being a receiving and transmitting node; since lattices within the ellipse are supposed to affect a result, giving a weight to the lattices within the ellipse; whereas, since the lattices outside the ellipse are supposed to have no influence on the result, giving a weight of 0 to the lattices outside the ellipse. 5. Detectiewerkwijze voor een abnormale graanomstandigheid op basis van radiotomografische beeldvorming volgens conclusie 1, met het kenmerk, dat in stap (7) een werkwijze voor het verkrijgen van formule (4) als volgt is: eerst, het als volgt minimaliseren van een objectieve functie min/| Art - WAX + of|TAxy||? (10), waar T een Tikhonov-matrix is en een voorgaande informatie van een modeloplossing voorstelt, |[TAx¢|| een strafterm voorstelt voor Tikhonov- regularisatie; en a een aanpasbare regularisatieparameter is, waarbij een waarde van a bepaalt of een uiteindelijke oplossing losgelaten is op gemeten gegevens of losgelaten is op de voorgaande informatie; vervolgens wordt een benaderende schatting van Ax, verkregen middels het nemen van een afgeleide van formule (9) om formule (4) als volgt te verkrijgen: at, = (WIW rar ry War, fa,The radiotomographic imaging abnormal grain condition detecting method according to claim 1, characterized in that in step (7), a method for obtaining formula (4) is as follows: first, minimizing an objective function as follows min/| Art - WAX + or|TAxy||? (10), where T is a Tikhonov matrix and represents a previous information of a model solution, |[TAx¢|| propose a penalty term for Tikhonov regularization; and a is an adjustable regularization parameter, wherein a value of a determines whether a final solution has been released on measured data or has been released on the preceding information; then an approximate estimate of Ax is obtained by taking a derivative of formula (9) to obtain formula (4) as follows: at, = (WIW rar ry War, fa,
NL2026449A 2019-09-11 2020-09-11 Abnormal grain condition detection method based on radio tomographic imaging NL2026449B1 (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910858421.0A CN110766561A (en) 2019-09-11 2019-09-11 Abnormal grain condition detection method based on wireless tomography technology

Publications (2)

Publication Number Publication Date
NL2026449A true NL2026449A (en) 2021-05-12
NL2026449B1 NL2026449B1 (en) 2021-10-14

Family

ID=69329739

Family Applications (1)

Application Number Title Priority Date Filing Date
NL2026449A NL2026449B1 (en) 2019-09-11 2020-09-11 Abnormal grain condition detection method based on radio tomographic imaging

Country Status (2)

Country Link
CN (1) CN110766561A (en)
NL (1) NL2026449B1 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112818514B (en) * 2021-01-05 2023-02-07 河南工业大学 Wireless tomography method and ellipse weight model combining horizontal distance attenuation
CN114428083A (en) * 2021-12-17 2022-05-03 中国电子科技集团公司第四十一研究所 On-line raw material cigarette packet mildew composite detection method

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101087218A (en) * 2006-06-07 2007-12-12 中国科学院电子学研究所 Grain information monitoring system based on wireless sensor network
CN101441805B (en) * 2007-11-21 2010-09-22 中国科学院电子学研究所 Grain condition monitoring system of sensor
US20130226660A1 (en) * 2010-03-04 2013-08-29 Lusine Yepremyan Cyclicality-Based Rules for Data Anomaly Detection
CN105374016B (en) * 2015-10-30 2018-10-26 沈阳工业大学 A kind of image rebuilding method of three-dimensional capacitance tomography
CN107422296B (en) * 2017-07-31 2018-07-13 河南工业大学 Wireless tomography determines method with grid pixel pad value
CN108490387B (en) * 2018-03-08 2021-11-30 河南工业大学 Grid pixel attenuation value self-adaptive determination method for wireless tomography

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ABDULLAH M S M ET AL: "A Review on Moisture Measurement Technique in Agricultural Silos", IOP CONF. SER.: MATER. SCI. ENG, 27 August 2019 (2019-08-27), XP055839325, Retrieved from the Internet <URL:https://iopscience.iop.org/article/10.1088/1757-899X/705/1/012001/pdf> [retrieved on 20210908] *
ASEFI MOHAMMAD ET AL: "Grain bin monitoring via electromagnetic imaging", COMPUTERS AND ELECTRONICS IN AGRICULTURE, ELSEVIER, AMSTERDAM, NL, vol. 119, 5 November 2015 (2015-11-05), pages 133 - 141, XP029320918, ISSN: 0168-1699, DOI: 10.1016/J.COMPAG.2015.10.016 *
HÜSEYIN YIGITLER ET AL: "Detector Based Radio Tomographic Imaging", ARXIV.ORG, CORNELL UNIVERSITY LIBRARY, 201 OLIN LIBRARY CORNELL UNIVERSITY ITHACA, NY 14853, 11 April 2016 (2016-04-11), XP080694737 *
LIM ET AL: "Rice Moisture Imaging Using Electromagnetic Measurement Technique", FOOD AND BIOPRODUCTS PROCESSING, INSTITUTION OF CHEMICAL ENGINEERS, RUGBY, GB, vol. 81, no. 3, 1 September 2003 (2003-09-01), pages 159 - 169, XP022525085, ISSN: 0960-3085, DOI: 10.1205/096030803322437929 *

Also Published As

Publication number Publication date
CN110766561A (en) 2020-02-07
NL2026449B1 (en) 2021-10-14

Similar Documents

Publication Publication Date Title
NL2026449B1 (en) Abnormal grain condition detection method based on radio tomographic imaging
CN102419437B (en) Track-before-detect method based on flight path inspection
CN105445701B (en) The pulse angle estimating method of DDMA MIMO radar targets
CN111999716B (en) Clutter prior information-based target adaptive fusion detection method
CN108957399A (en) A method of motion detection is realized using LTE signaling channel status information
CN107576963B (en) Estimation method of dual-polarization radar differential propagation phase shift based on particle filtering
CN108614245B (en) Radio frequency stealth-based interference method while tracking
CN109239683B (en) Clarmero boundary analysis method of broadband passive MIMO radar
WO2022151708A1 (en) Signal proxy-based multi-target constant false alarm rate measuring method
CN111954153B (en) Intelligent taxi calling navigation positioning method based on UWB positioning
CN106383339B (en) A kind of mirror target suppressing method of multi-site radar signal grade joint-detection
CN109521412B (en) Radar networking airspace target detection method based on local statistic fusion
CN108089166A (en) A kind of adaptive frequency domain detection method based on millimeter wave cloud detection radar
CN109655827B (en) Hidden target positioning method for single building corner
CN101881822A (en) Method for inhibiting same frequency interference of shared-spectrum radars
CN110007299A (en) A kind of dim target detection tracking based on hybrid coordinate puppet spectral technology
CN115113155A (en) Airborne distributed aperture coherent synthetic radar testing and evaluating method
CN108152796B (en) Main lobe moving interference elimination method based on gray Kalman filtering
Jagadesh et al. Modeling Target Detection and Performance Analysis of Electronic Countermeasures for Phased Radar.
CN111999715A (en) Target knowledge auxiliary self-adaptive fusion detection method under heterogeneous clutter
CN113093174B (en) PHD filter radar fluctuation weak multi-target-based pre-detection tracking method
CN108508413B (en) Target detection method based on probability statistics under low signal-to-noise ratio condition
CN112986975B (en) Distance weighting-based passive radar network centralized detection method
CN112083388B (en) Trace point quality evaluation method based on secondary radar response signal
CN115201863A (en) GEO satellite orbit maneuver detection method based on interferometric data