CN107480619A - The noise-reduction method and system of GPR B-scan image based on EEMD and arrangement entropy - Google Patents
The noise-reduction method and system of GPR B-scan image based on EEMD and arrangement entropy Download PDFInfo
- Publication number
- CN107480619A CN107480619A CN201710656970.0A CN201710656970A CN107480619A CN 107480619 A CN107480619 A CN 107480619A CN 201710656970 A CN201710656970 A CN 201710656970A CN 107480619 A CN107480619 A CN 107480619A
- Authority
- CN
- China
- Prior art keywords
- mtd
- mrow
- msub
- mtr
- noise
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 31
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000001174 ascending effect Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 238000000205 computational method Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000004083 survival effect Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005764 inhibitory process Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- -1 radius 0.2m Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10044—Radar image
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Image Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a kind of noise-reduction method and system of the GPR B-scan image based on EEMD and arrangement entropy, GPR two dimension B-scan picture signal is obtained first, then denoising is carried out to every one of B-scan picture signal of acquisition, to obtain the signal after the road B-scan picture signal suppresses noise, finally the signal after each road denoising is reconfigured, the two-dimentional B-scan picture signal after the noise that is inhibited;Wherein, the method for denoising is carried out for any one of B-scan picture signal to be included:EEMD decomposition is carried out for the road B-scan picture signal, the K IMF component arranged from high frequency to low frequency is obtained, calculates the arrangement entropy of each IMF components, choose arrangement entropy no more than the IMF components of preset value to reconstruct, obtain the signal after denoising.The present invention solves signal mode Aliasing Problem present in EMD decomposition, can effectively reduce noise.
Description
Technical field
The present invention relates to digital processing field is belonged to, the B-scan image procossing of GPR is related specifically to, specifically
The noise-reduction method and system of GPR B-scan image based on EEMD and arrangement entropy.
Background technology
GPR is a kind of discontinuity using underground medium come effective detection method of Underground target, extensively
Be applied to the field such as geology, resource, environment, engineering and military affairs.During actual detection, because underground medium structure is answered
Miscellaneous, physical parameter difference and the presence of noise jamming, the echo of GPR is unavoidably by noise pollution so that echo is believed
Number have non-stationary, nonlinear feature.Noise pollution causes certain difficulty to follow-up object detection and recognition, therefore visits
The noise suppressed of ground radar is particularly important for follow-up data processing.
Hilbert-Huang transform (Hilbert-Huang transform, HHT) is that one kind analyzes non-linear, non-stationary letter
Number new method, be suitable for handling the B-scan image of GPR.B-scan is to include the echo-signal of scanning
In two dimensional surface, abscissa represents the horizontal displacement direction of scanning, and ordinate represents the depth direction of scanning.This method handles B
During scan image, by empirical mode decomposition (Empirical mode decomposition, EMD) by signal decomposition to be intrinsic
Mode function (Intrinsic mode function, IMF), it can be obtained preferably when handling non-stationary, nonlinear signal
Effect.Although EMD methods have many advantages, but modal overlap in decomposition be present, cause to include in some IMF component
The signal of different scale, or similar magnitude signal are present in different IMF components, have impact on noise to a certain extent
Inhibition, cause to cause inhibition poor.
The content of the invention
The technical problem to be solved in the present invention is, presses down for above-mentioned EMD conversion when handling GPR B-scan image
The technological deficiency of effect processed difference, there is provided a kind of noise-reduction method of GPR B-scan image based on EEMD and arrangement entropy and
System.
According to the wherein one side of the present invention, the present invention is its technical problem of solution, there is provided one kind is based on EEMD and row
The noise-reduction method of the GPR B-scan image of row entropy, comprises the following steps:
S1, obtain GPR two dimension B-scan picture signal X ∈ RM×N, wherein M is road number, and N is the sampling per track data
Points;
S2, every one of B-scan picture signal to acquisition carry out denoising, are suppressed with obtaining the road B-scan picture signal
Signal after noise;Wherein, the method for denoising is carried out for any one of B-scan picture signal to be included:
S21, for the road B-scan picture signal { x (t), t=1,2 ..., N } carry out EEMD decomposition, obtain from high frequency to
K IMF components c of low frequency arrangement1(t)、c2(t)、…、cK(t);
The arrangement entropy of each IMF components in S22, calculation procedure S21;
S23, selection arrangement entropy are reconstructed no more than the IMF components of preset value, obtain the signal x ' (t) after denoising,
X ' (t) is the IMF component sums that arrangement entropy is not more than preset value;
S3, the signal after each road denoising is reconfigured, the two-dimentional B-scan picture signal X' after the noise that is inhibited
∈RM×N。
Further, comprised the following steps that in the S21 of the noise-reduction method of the present invention:
S211, default T noise signal is separately added into x (t), obtains T signals and associated noises, wherein
xi(t)=x (t)+ni(t),
In formula, i=1,2 ..., T, xi(t) it is i-th of signals and associated noises, ni(t) made an uproar for what is added in i-th of signals and associated noises
Sound, ni(t) it is that average is 0, standard deviation is the white Gaussian noise of constant;
S212, EMD decomposition is carried out respectively to each signals and associated noises, obtain the IMF components of each signals and associated noises;
S213, the corresponding IMF components for decomposing to obtain to T signals and associated noises using following formula are averaged, and obtain EEMD
The IMF components of decomposition;
In formula, cp(t) p-th of IMF component for decomposing to obtain by EEMD for x (t), ci,p(t) i-th of signals and associated noises warp
Cross p-th of IMF component that EMD decomposes to obtain, p=1,2,3 ..., K.
Further, in the step S22 of the noise-reduction method of the present invention, for any one x (t), the row of each IMF components
Row entropy computational methods are as follows:
S221, each IMF components arrived for decomposition carry out phase space reconfiguration, obtain matrix:
Wherein, i=1,2 ... Q, m and λ are respectively Embedded dimensions and time delay, and Q is that vector is reconstructed in phase space reconstruction
Number, Q=N- (m-1) λ;
S222, component (c will be reconstructedp(i),cp(i+λ),…,cp[i+ (m-1) λ]) arranged according to ascending order, it is as follows:
cp[i+(j1-1)λ]≤cp[i+(j2-1)λ]≤…≤cp[i+(jm-1)λ]
In formula, j1,j2,…,jmRepresent the index of each element column in reconstruct component;
For each vectorial c of reconstruct in phase space reconstructionp(i) one group of symbol for reflecting its element size order, is obtained respectively
Number sequence S (g)=(j1,j2,…,jm), wherein g=1,2 ..., q, q≤m!;
S223, the Probability p for calculating the appearance of each symbol sebolic addressing1,p2,…,pq, and according to the following formula sequence of calculation { cp
(i), i=1,2 ..., K arrangement entropy Hp(m):
Further, also include after the step S223 of the noise-reduction method of the present invention:
S224, the arrangement entropy H that will be calculated in S223p(m) new arrangement entropy is updated to after being normalized, new
Arrangement entropy normalized formula be:
Hp=Hp(m)/ln(m!).
Further, in the noise-reduction method of the present invention, c in step S211(t)、c2(t)、…、cK(t) it is according to from height
Frequency arranges to low frequency, is specifically comprised the following steps in step S23:
By the preset value successively compared with the arrangement entropy of each IMF components, when the arrangement entropy of k-th of IMF component
When value is less than or equal to the preset value, it is not more than preset value using each IMF components after k and k as the arrangement entropy
IMF components.
In order to solve the above technical problems, present invention also offers a kind of GPR B-scan based on EEMD and arrangement entropy
The noise reduction system of image, including:
Signal acquisition module, for obtaining GPR two dimension B-scan picture signal X ∈ RM×N, wherein M is road number, and N is
Sampling number per track data;
Signal denoising module, for carrying out denoising to every one of B-scan picture signal of acquisition, swept with obtaining road B
Retouch picture signal and suppress the signal after noise;
Signals revivification module, for being reconfigured to the signal after each road denoising, the two-dimentional B after the noise that is inhibited
Scan image signal X' ∈ RM×N;
Wherein, signal denoising module is handled any one of B-scan picture signal using following submodules:
EEMD decomposes submodule, for carrying out EEMD points for one of B-scan picture signal { x (t), t=1,2 ..., N }
Solution, obtain the K IMF components c arranged from high frequency to low frequency1(t)、c2(t)、…、cK(t);
Entropy calculating sub module is arranged, the arrangement entropy of each IMF components in submodule is decomposed for calculation procedure EEMD;
Signal reconstruction submodule, reconstructed for choosing arrangement entropy no more than the IMF components of preset value, after obtaining denoising
Signal x ' (t), x ' (t) be arrangement entropy be not more than preset value IMF component sums.
Further, in the noise reduction system of the present invention, EEMD decomposes submodule and specifically included:
Noise adding device, for being separately added into default T noise signal into x (t), T signals and associated noises are obtained, its
In
xi(t)=x (t)+ni(t),
In formula, i=1,2 ..., T, xi(t) it is i-th of signals and associated noises, ni(t) made an uproar for what is added in i-th of signals and associated noises
Sound, ni(t) it is that average is 0, standard deviation is the white Gaussian noise of constant;
EMD resolving cells, for carrying out EMD decomposition respectively to each signals and associated noises, obtain the IMF of each signals and associated noises
Component;
IMF converter units, the corresponding IMF components for decomposing to obtain to T signals and associated noises using following formula are put down
, the IMF components of EEMD decomposition are obtained;
In formula, cp(t) p-th of IMF component for decomposing to obtain by EEMD for x (t), ci,p(t) i-th of signals and associated noises warp
Cross p-th of IMF component that EMD decomposes to obtain, p=1,2,3 ..., K.
Further, in the noise reduction system of the present invention, for any one x (t), under arrangement entropy calculating sub module utilizes
State the arrangement entropy that unit calculates each IMF components:
Phase space reconfiguration unit, phase space reconfiguration is carried out for each IMF components arrived for decomposition, obtains matrix:
Wherein, i=1,2 ... Q, m and λ are respectively Embedded dimensions and time delay, and Q is that vector is reconstructed in phase space reconstruction
Number, Q=N- (m-1) λ;
Symbol sebolic addressing acquiring unit, for component (c will to be reconstructedp(i),cp(i+λ),…,cp[i+ (m-1) λ]) according to ascending order
Arrangement, it is as follows:
cp[i+(j1-1)λ]≤cp[i+(j2-1)λ]≤…≤cp[i+(jm-1)λ]
In formula, j1,j2,…,jmRepresent the index of each element column in reconstruct component;
For each vectorial c of reconstruct in phase space reconstructionp(i) one group of symbol for reflecting its element size order, is obtained respectively
Number sequence S (g)=(j1,j2,…,jm), wherein g=1,2 ..., q, q≤m!;
Entropy computing unit is arranged, for calculating the Probability p of each symbol sebolic addressing appearance1,p2,…,pq, and according to following
The formula sequence of calculation { cp(i), i=1,2 ..., K arrangement entropy Hp(m):
Further, in the noise reduction system of the present invention, arrangement entropy calculating sub module also includes:
Entropy updating block is arranged, for the arrangement entropy H calculated in entropy computing unit will to be arrangedp(m) place is normalized
New arrangement entropy is updated to after reason, the normalized formula of new arrangement entropy is:
Hp=Hp(m)/ln(m!).
Further, in the noise reduction system of the present invention, specifically following module is included in signal reconstruction submodule:
By the preset value successively compared with the arrangement entropy of each IMF components, when the arrangement entropy of k-th of IMF component
When value is less than or equal to the preset value, it is not more than preset value using each IMF components after k and k as the arrangement entropy
IMF components.
Implement the noise-reduction method and system of the GPR B-scan image based on EEMD and arrangement entropy of the present invention, have
Following beneficial effect:The present invention is decomposed using EEMD methods to Gpr Signal, obtains what is arranged from high frequency to low frequency
IMF, solve signal mode Aliasing Problem present in EMD decomposition, can effectively reduce noise.And it is further, calculating
During the arrangement entropy of each IMF components, determine that noise signal corresponds to IMF IMFs corresponding with echo signal separation by preset value, enter
And IMF corresponding to echo signal is reconstructed, it can effectively suppress noise, more retain target information.
Brief description of the drawings
Below in conjunction with drawings and Examples, the invention will be further described, in accompanying drawing:
Fig. 1 is an embodiment of the noise-reduction method of the GPR B-scan image based on EEMD and arrangement entropy of the present invention
Flow chart;
Fig. 2 is the original B-scan picture signal of GPR;
Fig. 3 is to add the image after making an uproar to original B-scan image;
Fig. 4 is to add the 40th signal and its EEMD exploded views after making an uproar;
Fig. 5 is the 40th signals and associated noises and suppresses the signal after noise;
Fig. 6 is the image to being reconfigured after all road signals progress denoising in Fig. 3;
Fig. 7 is an embodiment of the noise reduction system of the GPR B-scan image based on EEMD and arrangement entropy of the present invention
Theory diagram.
Embodiment
In order to which technical characteristic, purpose and the effect of the present invention is more clearly understood, now compares accompanying drawing and describe in detail
The embodiment of the present invention.
As shown in figure 1, it is the noise-reduction method based on EEMD and the GPR B-scan image for arranging entropy of the invention
The flow chart of one embodiment.In the noise-reduction method of the present embodiment, following steps is included:
S1, obtain GPR two dimension B-scan picture signal X ∈ RM×N, wherein M is road number, and N is the sampling per track data
Points.
S2, every one of B-scan picture signal to acquisition carry out denoising, are suppressed with obtaining the road B-scan picture signal
Signal after noise.Wherein, the method being reconstructed for any one of B-scan picture signal includes:
S21, carry out EEMD decomposition for one B-scan picture signal { x (t), t=1,2 ..., N }, obtain from high frequency to
K IMF components c of low frequency arrangement1(t)、c2(t)、…、cK(t).Wherein,
Wherein K be decompose number, namely EEMD decompose IMF components number, c1(t)、c2(t)、…、cK(t) it is point
The IMF components arranged from high frequency to low frequency that solution obtains, r (t) are to decompose the survival function finally obtained.
The arrangement entropy of each IMF components in S22, calculation procedure 2.Because frequency is higher, arrangement entropy is bigger, therefore arranges
The descending arrangement of entropy.
S23, the arrangement entropy according to each IMF components, selected threshold Th determine that noise signal corresponds to IMF and echo signal pair
IMF separation k, 1≤k≤K, the separation k obtained according to step are answered, is entered using k-th and k-th later IMF component
Row reconstruct, obtains the signal x ' (t) after denoising.K-th and k-th later IMF component be and echo signal (non-noise letter
Number) corresponding to IMF components, the IMF components before k-th are IMF components corresponding with noise signal.Wherein,
Wherein x ' (t) is suppression noise Hou Gai roads B-scan picture signal.
Specifically, can be by the preset value successively compared with the arrangement entropy of each IMF components, when k-th of IMF component
Arrangement entropy when being less than or equal to the preset value, each IMF components after k and k is little as the arrangement entropy
In the IMF components of preset value.
S3, the signal after each road denoising is reconfigured, the two-dimentional B-scan picture signal X' after the noise that is inhibited
∈RM×N。
In the step S21 of the present embodiment, one of Gpr Signal EEMD decomposition method is as follows:
S211, for one of primary signal x (t), be separately added into default T noise signal, obtain T signals and associated noises,
Wherein:
xi(t)=x (t)+ni(t) i=1,2 ..., T
Wherein xi(t) it is i-th of signals and associated noises, ni(t) it is the noise added in i-th of signals and associated noises, it is that average is 0,
Standard deviation is the white Gaussian noise of constant.
S212, to xi(t) EMD decomposition is carried out, obtains a series of IMF components, wherein:
Wherein K be IMF components quantity, ci,p(t) p-th of IMF component of i-th of signals and associated noises, ri(t) contain for i-th
The survival function of noise cancellation signal.
S213, the corresponding IMF components for decomposing to obtain to T signals and associated noises are averaged, and obtain the IMF components of EEMD decomposition
For:
Wherein cp(t) p-th of IMF component for decomposing to obtain by EEMD for primary signal.
Primary signal x (t) may finally be expressed as by what EEMD was decomposed:
Wherein r (t) is final survival function, represents the average tendency of primary signal.
In the step S22 of the present embodiment, the arrangement entropy computational methods of each IMF components are as follows:
S221, the IMF components { c arrived for decompositionp(i), i=1,2 ..., K } phase space reconfiguration is carried out, obtain square
Battle array:
Wherein, i=1,2 ... Q, m and λ are respectively Embedded dimensions and time delay, and Q is that vector is reconstructed in phase space reconstruction
Number, Q=N- (m-1) λ;
S222, component (c will be reconstructedp(i),cp(i+λ),…,cp[i+ (m-1) λ]) arranged according to ascending order, it is as follows:
cp[i+(j1-1)λ]≤cp[i+(j2-1)λ]≤…≤cp[i+(jm-1)λ]
In formula, j1,j2,…,jmRepresent the index of each element column in reconstruct component.
Accordingly, for any vectorial c of reconstruct in phase space reconstructionp(i), can obtain reflecting the one of its element size order
Group code sequence S (g)=(j1,j2,…,jm), wherein g=1,2 ..., q, q≤m!.M different symbol j1,j2,…,jmAltogether
There is m!The different symbol sebolic addressing of kind, S (g) is m!One kind in the different symbol sebolic addressing of kind;
S223, the Probability p for calculating the appearance of each symbol sebolic addressing1,p2,…,pq, and according to the following formula sequence of calculation { cp
(i), i=1,2 ..., K arrangement entropy Hp(m):
In actual treatment, those skilled in the art are accustomed to Hp(m) normalized is done, using normalized arrangement entropy as row
The end value of row entropy.Wherein, p is worked asg=1/m!When, Hp(m) maximum ln (m are reached!), Hp(m) use is normalized
Following formula are carried out:
Hp=Hp(m)/ln(m!).
The B-scan image for producing GPR is emulated using FDTD methods, as shown in Figure 2.Simulation parameters are as follows:
(1) underground medium is concrete, and its relative dielectric constant is 6.0, and center of antenna frequency is 900MHz;
(2) simulating area width is 3m, depth 2m;Three targets are ideal cylindrical conductor, radius 0.2m, water
Prosposition puts respectively 0.9m, 1.5m and 2.1m, and depth is 0.6m;There are an a length of 1m, a width of 0.4m air at depth 1.5m
The rectangular area of filling;
(3) B-scan image shares 80, and per pass has 2036 sampled points;
Processing of making an uproar is added to the GPR B-scan image of emulation, obtains noisy image, its signal to noise ratio is 0.985dB, is such as schemed
Shown in 3.
40th signals and associated noises are carried out with EEMD decomposition, primary signal is with decomposing obtained each IMF components such as Fig. 4 institutes
Show.
The arrangement entropy of each IMF components is calculated, determines that noise corresponds to IMF IMFs corresponding with echo signal separation, step
It is as follows:
(1) calculating EEMD decomposes the arrangement entropy of obtained each IMF components;
(2) it is 0.4 to set threshold value Th, and the arrangement entropy of each IMF components is obtained into noise signal compared with threshold value successively
IMF separations k with echo signal is 5;
(3) the separation k=5 obtained according to step (2), IMF components to the 5th and behind are reconstructed, obtained
Signal after denoising, as shown in Figure 5.
Pair plus all road signals after making an uproar handle, the ground penetrating radar image after the noise that is inhibited, its signal to noise ratio is
15.385dB as shown in Figure 6.
With reference to figure 7, it is the one of the noise reduction system based on EEMD and the GPR B-scan image for arranging entropy of the invention
The theory diagram of embodiment, the noise reduction system of the present embodiment include signal acquisition module 71, signal denoising module 72 and signal also
Former module 73.
Signal acquisition module 71 obtains GPR two dimension B-scan picture signal X ∈ RM×N, signal denoising module 72, to obtaining
Every one of B-scan picture signal taken carries out denoising, to obtain the signal after the road B-scan picture signal suppresses noise, letter
Number recovery module 73, is reconfigured to the signal after each road denoising, the two-dimentional B-scan image letter after the noise that is inhibited
Number.
Wherein, signal denoising module 72 decomposes submodule, arrangement entropy calculating sub module and signal weight using following EEMD
Every one of B-scan picture signal is reconstructed structure submodule.EEMD decomposes submodule to any one of B-scan picture signal { x
(t), t=1,2 ..., N } EEMD decomposition is carried out, obtain the K IMF components c arranged from high frequency to low frequency1(t)、c2(t)、…、cK
(t) the arrangement entropy that entropy calculating sub module calculation procedure EEMD decomposes each IMF components in submodule, signal reconstruction submodule, are arranged
Block chooses arrangement entropy no more than the IMF components of preset value to reconstruct, and obtains the signal x ' (t) after denoising, and x ' (t) is arrangement entropy
IMF component sum of the value no more than preset value.
In the present embodiment, EEMD decomposes submodule specifically comprising noise adding device, EMD resolving cells and IMF conversion
Unit.Noise adding device is separately added into default T noise signal into x (t), obtains T signals and associated noises, and EMD decomposes single
Member carries out EMD decomposition respectively to each signals and associated noises, obtains the IMF components of each signals and associated noises, IMF converter units are to T
The corresponding IMF components that signals and associated noises decompose to obtain are averaged, and obtain the IMF components of EEMD decomposition.
In the present embodiment, calculated using following units each IMF points for any one x (t), arrangement entropy calculating sub module
The arrangement entropy of amount:
Phase space reconfiguration unit, phase space reconfiguration is carried out for each IMF components arrived for decomposition, obtains matrix:
Wherein, i=1,2 ... Q, m and λ are respectively Embedded dimensions and time delay, and Q is that vector is reconstructed in phase space reconstruction
Number, Q=N- (m-1) λ;
Symbol sebolic addressing acquiring unit, for component (c will to be reconstructedp(i),cp(i+λ),…,cp[i+ (m-1) λ]) according to ascending order
Arrangement, it is as follows:
cp[i+(j1-1)λ]≤cp[i+(j2-1)λ]≤…≤cp[i+(jm-1)λ]
In formula, j1,j2,…,jmRepresent the index of each element column in reconstruct component;
For each vectorial c of reconstruct in phase space reconstructionp(i) one group of symbol for reflecting its element size order, is obtained respectively
Number sequence S (g)=(j1,j2,…,jm), wherein g=1,2 ..., q, q≤m!;
Entropy computing unit is arranged, for calculating the Probability p of each symbol sebolic addressing appearance1,p2,…,pq, and according to following
The formula sequence of calculation { cp(i), i=1,2 ..., K arrangement entropy Hp(m):
Entropy updating block is arranged, for the arrangement entropy H calculated in entropy computing unit will to be arrangedp(m) place is normalized
New arrangement entropy is updated to after reason, the normalized formula of new arrangement entropy is:
Hp=Hp(m)/ln(m!).
Embodiments of the invention are described above in conjunction with accompanying drawing, but the invention is not limited in above-mentioned specific
Embodiment, above-mentioned embodiment is only schematical, rather than restricted, one of ordinary skill in the art
Under the enlightenment of the present invention, in the case of present inventive concept and scope of the claimed protection is not departed from, it can also make a lot
Form, these are belonged within the protection of the present invention.
Claims (10)
1. a kind of noise-reduction method of the GPR B-scan image based on EEMD and arrangement entropy, it is characterised in that including following step
Suddenly:
S1, obtain GPR two dimension B-scan picture signal X ∈ RM×N, wherein M is road number, and N is the sampling number per track data;
S2, every one of B-scan picture signal to acquisition carry out denoising, suppress noise to obtain the road B-scan picture signal
Signal afterwards;Wherein, the method for denoising is carried out for any one of B-scan picture signal to be included:
S21, EEMD decomposition is carried out for the road B-scan picture signal { x (t), t=1,2 ..., N }, obtained from high frequency to low frequency
K IMF components c of arrangement1(t)、c2(t)、…、cK(t);
The arrangement entropy of each IMF components in S22, calculation procedure S21;
S23, selection arrangement entropy are reconstructed no more than the IMF components of preset value, obtain the signal x ' (t) after denoising, x ' (t)
It is not more than the IMF component sums of preset value for arrangement entropy;
S3, the signal after each road denoising is reconfigured, the two-dimentional B-scan picture signal X' ∈ R after the noise that is inhibitedM ×N。
2. noise-reduction method according to claim 1, it is characterised in that the S21's comprises the following steps that:
S211, default T noise signal is separately added into x (t), obtains T signals and associated noises, wherein
xi(t)=x (t)+ni(t),
In formula, i=1,2 ..., T, xi(t) it is i-th of signals and associated noises, ni(t) it is the noise added in i-th of signals and associated noises, ni
(t) it is that average is 0, standard deviation is the white Gaussian noise of constant;
S212, EMD decomposition is carried out respectively to each signals and associated noises, obtain the IMF components of each signals and associated noises;
S213, the corresponding IMF components for decomposing to obtain to T signals and associated noises using following formula are averaged, and obtain EEMD decomposition
IMF components;
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<msub>
<mi>c</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
In formula, cp(t) p-th of IMF component for decomposing to obtain by EEMD for x (t), ci,p(t) i-th of signals and associated noises passes through EMD
Decompose p-th obtained of IMF components, p=1,2,3 ..., K.
3. noise-reduction method according to claim 1, it is characterised in that in the step S22, for any one x (t),
The arrangement entropy computational methods of each IMF components are as follows:
S221, each IMF components arrived for decomposition carry out phase space reconfiguration, obtain matrix:
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>+</mo>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>+</mo>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
Wherein, i=1,2 ... Q, m and λ are respectively Embedded dimensions and time delay, and Q is to reconstruct vectorial number in phase space reconstruction,
Q=N- (m-1) λ;
S222, component (c will be reconstructedp(i),cp(i+λ),…,cp[i+ (m-1) λ]) arranged according to ascending order, it is as follows:
cp[i+(j1-1)λ]≤cp[i+(j2-1)λ]≤…≤cp[i+(jm-1)λ]
In formula, j1,j2,…,jmRepresent the index of each element column in reconstruct component;
For each vectorial c of reconstruct in phase space reconstructionp(i) the group code sequence for reflecting its element size order, is obtained respectively
S (g)=(j1,j2,…,jm), wherein g=1,2 ..., q, q≤m!;
S223, the Probability p for calculating the appearance of each symbol sebolic addressing1,p2,…,pq, and according to the following formula sequence of calculation { cp(i),i
=1,2 ..., K arrangement entropy Hp(m):
<mrow>
<msub>
<mi>H</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>g</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>q</mi>
</munderover>
<msub>
<mi>p</mi>
<mi>g</mi>
</msub>
<mi>ln</mi>
<mi> </mi>
<msub>
<mi>p</mi>
<mi>g</mi>
</msub>
<mo>.</mo>
</mrow>
4. noise-reduction method according to claim 3, it is characterised in that also include after the step S223:
S224, the arrangement entropy H that will be calculated in S223p(m) new arrangement entropy, new arrangement are updated to after being normalized
The normalized formula of entropy is:
Hp=Hp(m)/ln(m!).
5. noise-reduction method according to claim 1, it is characterised in that c in step S211(t)、c2(t)、…、cK(t) be by
Arrange according to from high frequency to low frequency, specifically comprised the following steps in step S23:
By the preset value successively compared with the arrangement entropy of each IMF components, when the arrangement entropy of k-th of IMF component is small
When the preset value, the IMF of preset value is not more than using each IMF components after k and k as the arrangement entropy
Component.
A kind of 6. noise reduction system of the GPR B-scan image based on EEMD and arrangement entropy, it is characterised in that including:
Signal acquisition module, for obtaining GPR two dimension B-scan picture signal X ∈ RM×N, wherein M is road number, and N is per pass
The sampling number of data;
Signal denoising module, for carrying out denoising to every one of B-scan picture signal of acquisition, to obtain the road B-scan figure
Suppress the signal after noise as signal;
Signals revivification module, for being reconfigured to the signal after each road denoising, the two-dimentional B-scan after the noise that is inhibited
Picture signal X' ∈ RM×N;
Wherein, signal denoising module is handled any one of B-scan picture signal using following submodules:
EEMD decomposes submodule, for carrying out EEMD decomposition for one of B-scan picture signal { x (t), t=1,2 ..., N }, obtains
To the K IMF components c arranged from high frequency to low frequency1(t)、c2(t)、…、cK(t);
Entropy calculating sub module is arranged, the arrangement entropy of each IMF components in submodule is decomposed for calculation procedure EEMD;
Signal reconstruction submodule, reconstructed for choosing arrangement entropy no more than the IMF components of preset value, obtain the letter after denoising
Number x ' (t), x ' (t) are the IMF component sums that arrangement entropy is not more than preset value.
7. noise reduction system according to claim 6, it is characterised in that the EEMD decomposes submodule and specifically included:
Noise adding device, for being separately added into default T noise signal into x (t), T signals and associated noises are obtained, wherein
xi(t)=x (t)+ni(t),
In formula, i=1,2 ..., T, xi(t) it is i-th of signals and associated noises, ni(t) it is the noise added in i-th of signals and associated noises, ni
(t) it is that average is 0, standard deviation is the white Gaussian noise of constant;
EMD resolving cells, for carrying out EMD decomposition respectively to each signals and associated noises, obtain IMF points of each signals and associated noises
Amount;
IMF converter units, the corresponding IMF components for decomposing to obtain to T signals and associated noises using following formula are averaged, obtained
The IMF components decomposed to EEMD;
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<msub>
<mi>c</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
In formula, cp(t) p-th of IMF component for decomposing to obtain by EEMD for x (t), ci,p(t) i-th of signals and associated noises passes through EMD
Decompose p-th obtained of IMF component, p=1,2,3 ..., K.
8. noise reduction system according to claim 6, it is characterised in that for any one x (t), arrangement entropy calculates submodule
Block calculates the arrangement entropy of each IMF components using following units:
Phase space reconfiguration unit, phase space reconfiguration is carried out for each IMF components arrived for decomposition, obtains matrix:
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>+</mo>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>+</mo>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
Wherein, i=1,2 ... Q, m and λ are respectively Embedded dimensions and time delay, and Q is to reconstruct vectorial number in phase space reconstruction,
Q=N- (m-1) λ;
Symbol sebolic addressing acquiring unit, for component (c will to be reconstructedp(i),cp(i+λ),…,cp[i+ (m-1) λ]) arranged according to ascending order
Row, it is as follows:
cp[i+(j1-1)λ]≤cp[i+(j2-1)λ]≤…≤cp[i+(jm-1)λ]
In formula, j1,j2,…,jmRepresent the index of each element column in reconstruct component;
For each vectorial c of reconstruct in phase space reconstructionp(i) the group code sequence for reflecting its element size order, is obtained respectively
S (g)=(j1,j2,…,jm), wherein g=1,2 ..., q, q≤m!;
Entropy computing unit is arranged, for calculating the Probability p of each symbol sebolic addressing appearance1,p2,…,pq, and according to following formula
The sequence of calculation { cp(i), i=1,2 ..., K arrangement entropy Hp(m):
<mrow>
<msub>
<mi>H</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>g</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>q</mi>
</munderover>
<msub>
<mi>p</mi>
<mi>g</mi>
</msub>
<mi>ln</mi>
<mi> </mi>
<msub>
<mi>p</mi>
<mi>g</mi>
</msub>
<mo>.</mo>
</mrow>
9. noise reduction system according to claim 8, it is characterised in that arrangement entropy calculating sub module also includes:
Entropy updating block is arranged, for the arrangement entropy H calculated in entropy computing unit will to be arrangedp(m) after being normalized more
It is newly new arrangement entropy, the normalized formula of new arrangement entropy is:
Hp=Hp(m)/ln(m!).
10. noise reduction system according to claim 6, it is characterised in that following mould is specifically included in signal reconstruction submodule
Block:
By the preset value successively compared with the arrangement entropy of each IMF components, when the arrangement entropy of k-th of IMF component is small
When the preset value, the IMF of preset value is not more than using each IMF components after k and k as the arrangement entropy
Component.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710656970.0A CN107480619B (en) | 2017-08-03 | 2017-08-03 | Noise-reduction method and system based on EEMD and the Ground Penetrating Radar B-scan image for arranging entropy |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710656970.0A CN107480619B (en) | 2017-08-03 | 2017-08-03 | Noise-reduction method and system based on EEMD and the Ground Penetrating Radar B-scan image for arranging entropy |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107480619A true CN107480619A (en) | 2017-12-15 |
CN107480619B CN107480619B (en) | 2019-11-15 |
Family
ID=60597606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710656970.0A Active CN107480619B (en) | 2017-08-03 | 2017-08-03 | Noise-reduction method and system based on EEMD and the Ground Penetrating Radar B-scan image for arranging entropy |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107480619B (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108267784A (en) * | 2018-01-18 | 2018-07-10 | 国网江西省电力有限公司电力科学研究院 | A kind of seismic signal random noise compression process method |
CN108594161A (en) * | 2018-05-03 | 2018-09-28 | 国网重庆市电力公司电力科学研究院 | Foreign matter voice signal noise-reduction method, system in a kind of electric energy meter |
CN110096956A (en) * | 2019-03-25 | 2019-08-06 | 中国地质大学(武汉) | Signal antinoise method and device based on EEMD and arrangement entropy second differnce |
CN110112757A (en) * | 2019-05-17 | 2019-08-09 | 福州大学 | Low-frequency oscillation analysis method based on SURE wavelet noise and improvement HHT |
CN110659621A (en) * | 2019-09-27 | 2020-01-07 | 山东科技大学 | Combined noise reduction method based on variational modal decomposition and permutation entropy |
CN111380680A (en) * | 2020-03-10 | 2020-07-07 | 昆明理工大学 | Check valve fault feature extraction method based on improved permutation entropy |
CN113297932A (en) * | 2021-05-11 | 2021-08-24 | 中铁第四勘察设计院集团有限公司 | Satellite data denoising method, device, equipment and storage medium |
CN113702968A (en) * | 2021-08-03 | 2021-11-26 | 煤炭科学研究总院 | Life detection method and device |
CN113822235A (en) * | 2021-11-22 | 2021-12-21 | 中南大学 | CEEMD denoising and signal reconstructing method and device based on permutation entropy standard measurement |
CN116340731A (en) * | 2023-05-29 | 2023-06-27 | 北京科技大学 | Noise reduction method and system for rock acoustic emission signals |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104198997A (en) * | 2014-09-29 | 2014-12-10 | 中国矿业大学(北京) | Data registering method for ground penetrating radar based on mutual information |
US20150164413A1 (en) * | 2013-12-13 | 2015-06-18 | National Chung Shan Institute Of Science And Technology | Method of creating anesthetic consciousness index with artificial neural network |
CN104748961A (en) * | 2015-03-30 | 2015-07-01 | 中国矿业大学 | Gear fault diagnosis method based on SVD decomposition and noise reduction and correlation EEMD entropy features |
CN106251368A (en) * | 2016-06-12 | 2016-12-21 | 中国科学院遥感与数字地球研究所 | SAR image based on BEMD and the fusion method of multispectral image |
CN106644484A (en) * | 2016-09-14 | 2017-05-10 | 西安工业大学 | Turboprop Engine rotor system fault diagnosis method through combination of EEMD and neighborhood rough set |
-
2017
- 2017-08-03 CN CN201710656970.0A patent/CN107480619B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150164413A1 (en) * | 2013-12-13 | 2015-06-18 | National Chung Shan Institute Of Science And Technology | Method of creating anesthetic consciousness index with artificial neural network |
CN104198997A (en) * | 2014-09-29 | 2014-12-10 | 中国矿业大学(北京) | Data registering method for ground penetrating radar based on mutual information |
CN104748961A (en) * | 2015-03-30 | 2015-07-01 | 中国矿业大学 | Gear fault diagnosis method based on SVD decomposition and noise reduction and correlation EEMD entropy features |
CN106251368A (en) * | 2016-06-12 | 2016-12-21 | 中国科学院遥感与数字地球研究所 | SAR image based on BEMD and the fusion method of multispectral image |
CN106644484A (en) * | 2016-09-14 | 2017-05-10 | 西安工业大学 | Turboprop Engine rotor system fault diagnosis method through combination of EEMD and neighborhood rough set |
Non-Patent Citations (1)
Title |
---|
周涛涛: "《基于CEEMD和排列熵的故障数据小波阈值降噪方法》", 《振动与冲击》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108267784A (en) * | 2018-01-18 | 2018-07-10 | 国网江西省电力有限公司电力科学研究院 | A kind of seismic signal random noise compression process method |
CN108594161A (en) * | 2018-05-03 | 2018-09-28 | 国网重庆市电力公司电力科学研究院 | Foreign matter voice signal noise-reduction method, system in a kind of electric energy meter |
CN108594161B (en) * | 2018-05-03 | 2020-06-19 | 国网重庆市电力公司电力科学研究院 | Noise reduction method and system for foreign matter sound signals in electric energy meter |
CN110096956A (en) * | 2019-03-25 | 2019-08-06 | 中国地质大学(武汉) | Signal antinoise method and device based on EEMD and arrangement entropy second differnce |
CN110096956B (en) * | 2019-03-25 | 2023-02-10 | 中国地质大学(武汉) | Signal denoising method and device based on EEMD and permutation entropy second-order difference |
CN110112757B (en) * | 2019-05-17 | 2022-04-12 | 福州大学 | Low-frequency oscillation analysis method based on SURE wavelet denoising and improved HHT |
CN110112757A (en) * | 2019-05-17 | 2019-08-09 | 福州大学 | Low-frequency oscillation analysis method based on SURE wavelet noise and improvement HHT |
CN110659621A (en) * | 2019-09-27 | 2020-01-07 | 山东科技大学 | Combined noise reduction method based on variational modal decomposition and permutation entropy |
CN111380680A (en) * | 2020-03-10 | 2020-07-07 | 昆明理工大学 | Check valve fault feature extraction method based on improved permutation entropy |
CN113297932A (en) * | 2021-05-11 | 2021-08-24 | 中铁第四勘察设计院集团有限公司 | Satellite data denoising method, device, equipment and storage medium |
CN113702968A (en) * | 2021-08-03 | 2021-11-26 | 煤炭科学研究总院 | Life detection method and device |
CN113822235A (en) * | 2021-11-22 | 2021-12-21 | 中南大学 | CEEMD denoising and signal reconstructing method and device based on permutation entropy standard measurement |
CN113822235B (en) * | 2021-11-22 | 2022-02-18 | 中南大学 | CEEMD denoising and signal reconstructing method and device based on permutation entropy standard measurement |
CN116340731A (en) * | 2023-05-29 | 2023-06-27 | 北京科技大学 | Noise reduction method and system for rock acoustic emission signals |
CN116340731B (en) * | 2023-05-29 | 2023-08-22 | 北京科技大学 | Noise reduction method and system for rock acoustic emission signals |
Also Published As
Publication number | Publication date |
---|---|
CN107480619B (en) | 2019-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107480619A (en) | The noise-reduction method and system of GPR B-scan image based on EEMD and arrangement entropy | |
US6292592B1 (en) | Efficient multi-resolution space-time adaptive processor | |
CN106408530A (en) | Sparse and low-rank matrix approximation-based hyperspectral image restoration method | |
CN102509263B (en) | K-SVD (K-means singular value decomposition) speckle inhibiting method based on SAR (synthetic aperture radar) image local statistic characteristic | |
CN101482617A (en) | Synthetic aperture radar image denoising method based on non-down sampling profile wave | |
CN113361367B (en) | Underground target electromagnetic inversion method and system based on deep learning | |
CN110830043B (en) | Image compressed sensing reconstruction method based on mixed weighted total variation and non-local low rank | |
CN112578471B (en) | Clutter noise removing method for ground penetrating radar | |
CN108230280A (en) | Image speckle noise minimizing technology based on tensor model and compressive sensing theory | |
CN110780274A (en) | Improved L for scanning radar 1Regularization azimuth super-resolution imaging method | |
Fukuda et al. | Smoothing effect of wavelet-based speckle filtering: The Haar basis case | |
CN108428221A (en) | A kind of neighborhood bivariate shrinkage function denoising method based on shearlet transformation | |
Yüksel et al. | A simple generalized neuro-fuzzy operator for efficient removal of impulse noise from highly corrupted digital images | |
CN111461999B (en) | SAR image speckle suppression method based on super-pixel similarity measurement | |
CN113935246A (en) | Signal robust sparse time-frequency analysis method, terminal equipment and storage medium | |
Caracena | Analytic approximation of discrete field samples with weighted sums and the gridless computation of field derivatives | |
CN112906716A (en) | Noisy SAR image target identification method based on wavelet de-noising threshold self-learning | |
CN113064165A (en) | Scanning radar pitch-azimuth two-dimensional super-resolution method | |
CN116577749A (en) | Scanning radar super-resolution method under unknown broadening of antenna pattern | |
CN116027293A (en) | Rapid sparse angle super-resolution method for scanning radar | |
CN114998137A (en) | Ground penetrating radar image clutter suppression method based on generation countermeasure network | |
Maarel | On the interpretability of ordination diagrams | |
Parkes et al. | The compression of raw SAR and SAR image data | |
CN108564098B (en) | Polarization SAR classification method based on scattering complete convolution model | |
CN116068500B (en) | Novel space-time interference suppression method for multichannel synthetic aperture radar |
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 |