CN104615877A - Method for conducting signal denoising based on wavelet packet - Google Patents

Method for conducting signal denoising based on wavelet packet Download PDF

Info

Publication number
CN104615877A
CN104615877A CN201510042993.3A CN201510042993A CN104615877A CN 104615877 A CN104615877 A CN 104615877A CN 201510042993 A CN201510042993 A CN 201510042993A CN 104615877 A CN104615877 A CN 104615877A
Authority
CN
China
Prior art keywords
wavelet packet
signal
threshold
node
wavelet
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.)
Pending
Application number
CN201510042993.3A
Other languages
Chinese (zh)
Inventor
赵瑞山
杨帆
吴作启
徐佳
张羽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Liaoning Technical University
Original Assignee
Liaoning Technical University
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 Liaoning Technical University filed Critical Liaoning Technical University
Priority to CN201510042993.3A priority Critical patent/CN104615877A/en
Publication of CN104615877A publication Critical patent/CN104615877A/en
Pending legal-status Critical Current

Links

Landscapes

  • Data Exchanges In Wide-Area Networks (AREA)

Abstract

The invention relates to a method for conducting signal denoising based on a wavelet packet. The method is mainly applied to signal processing in engineering surveying. Low-frequency parts and high-frequency parts of signals are respectively processed through wavelet packet denoising, and the resolution ratios of high frequencies and low frequencies are improved. When a signal processing method used at present is used for conducting denoising processing on the signals, negative phenomena of low speed and low efficiency will be caused by limitation of data; meanwhile, as hard threshold functions and soft threshold functions are numerously adopted in data analysis so as to conduct quantization processing on decomposition coefficients, the denoising signals are lack of fidelity, and the relatively-ideal denoising signals cannot be obtained.

Description

A kind of method of carrying out signal denoising based on wavelet packet
Technical field
This method relates to a kind of method of carrying out signal denoising based on wavelet packet, is mainly used in the signal transacting in engineering survey.
Background technology
Wavelet packet denoising theory is amplified by Wavelet Denoising Method, and it processes respectively to the low frequency of signal and HFS, improves the resolution of low-and high-frequency.When the signal processing method of current use carries out denoising to signal, due to the limitation of data, slow, the inefficient negative phenomena of speed can be brought, simultaneously because majority in data analysis all adopts hard threshold function or soft-threshold function to carry out quantification treatment to coefficient of dissociation, it can be denoised signal distortion, can not obtain the denoised signal of relative ideal.
Summary of the invention
The object of the invention is to the defect overcoming prior art, a kind of method of carrying out signal denoising based on wavelet packet is provided.
To achieve these goals, present invention employs following technical scheme:
Wavelet packet analysis has carried out meticulousr analysis on the basis of wavelet analysis, frequency band is carried out multilayer division by it, the high frequency series part that multi-frequency analysis is not segmented is decomposed further, and frequency band can be chosen according to the feature of analyzed signal is adaptive, make it to match with signal spectrum, thus when improving-frequency resolution, there is the prospect more wide than wavelet analysis.Effective wavelet packet analysis directly perceived is exactly carry out a threshold value quantizing to each node coefficient that WAVELET PACKET DECOMPOSITION obtains, and then to the coefficient reconstruct after quantification, thus reaches the object of denoising.
Signal antinoise method based on wavelet packet analysis will experience four basic steps: 1. the WAVELET PACKET DECOMPOSITION of signal, selects a small echo, then carries out WAVELET PACKET DECOMPOSITION to signals and associated noises; 2. choose cost function, conventional cost function is information entropy and logarithmic entropy; 3. the choosing of best wavelet packet basis, for a given entropy standard, adopts top-down searching algorithm to calculate optimal tree, chooses best wavelet packet basis; 4. selected threshold, adopts different choosing methods to obtain different threshold values according to different decomposition frequency domains; 5. the threshold value quantizing of WAVELET PACKET DECOMPOSITION coefficient, for the coefficient of each WAVELET PACKET DECOMPOSITION, adopts the soft-threshold function of the improvement threshold value corresponding with this coefficient to carry out quantification treatment; 6. the wavelet package reconstruction of signal, according to the WAVELET PACKET DECOMPOSITION coefficient of lowermost layer and the coefficient through quantification treatment, carries out wavelet package reconstruction.
In above-mentioned each step, it is crucial that how selected threshold and how to carry out threshold value quantizing, to a certain extent, he is directly connected to the quality of signal being carried out to noise reduction process, and accompanying drawing 1 carries out the concrete steps of signal denoising for the method.
Wavelet packet basis storehouse is made up of many wavelet packet basiss, and different wavelet packet basiss has different character, reacts different signal characteristics, and we wish to choose optimum wavelet packet basis to specific signal, is used for expressing the feature of this signal.From the viewpoint of signal transacting, the search procedure of optimal base is in fact with as far as possible few coefficient, reacts information as much as possible, to reach the object of characteristics extraction.Gap between coefficient is the bigger the better, if only have a few coefficients very large, so these coefficients just can the feature of representation signal, and contrary coefficients comparison is average, and gap not quite, is so just difficult to the feature finding out signal.
When searching for optimal base, first will provide the cost function of a burst, cost function can be defined as any function about burst, and find the base making cost function minimum in all wavelet packet basiss then in wavelet library, this base is just optimal base.In order to reduce calculated amount, improve program operation speed, reduce unknown parameter simultaneously, this method adopts the top-down Binomial Trees based on " beta pruning " to be that cost function carries out choosing of optimal base with logarithmic entropy, that is logarithmic entropy that is that only have the logarithm entropy of lower one deck and that be less than this layer is just decomposed, otherwise with regard to beta pruning, detailed algorithm flow is as follows: to signal node U n jcarry out WAVELET PACKET DECOMPOSITION, calculate each node logarithmic entropy; comparison node U n jentropy and its two child node U decomposed 2n j+1, U 2n+1 j+1entropy and, and mark the less node of entropy; decision node U n jwhether be labeled, if be labeled, then deletion of node U 2n j+1, U 2n+1 j+1decompose downwards no longer further, if be not labeled, then by node U 2n j+1, U 2n+1 j+1as node U n jrepeat above step.
Threshold denoising is that a kind of realization is simple, the good Wavelet Package Denoising Method of effect.Its basic thought be to wavelet decomposition after each layer coefficients in the mould coefficient that is greater than or less than certain threshold value process respectively, then carry out inverse transformation reconstruct denoising after signal.
The selection of threshold value is the key point of denoising, and threshold value is selected large, and the feature of useful signal will be filtered, if contrary Threshold selection is too small, and still noisy existence in the signal after denoising.Because the selection of threshold value directly affects de-noising effect, people propose an a lot of empirical model: soft-threshold, hard-threshold, minimax threshold value, fixed threshold etc., but all models are not general, have its scope of application.
In Practical Project signal, useful information mainly concentrates on low frequency part, and HFS majority is noise, only containing a small amount of useful information.If choose unified threshold value to process coefficient of dissociation, be difficult to reflect the different characteristic of signal at different frequency bands.Therefore need to choose multiple threshold value to carry out quantification treatment for different frequency bands.Natural order and the frequency band order of WAVELET PACKET DECOMPOSITION tree node there are differences, and want to choose different threshold value for different frequency bands, first will understand fully the frequency band order of WAVELET PACKET DECOMPOSITION tree node.When wavelet packet often decomposes one deck, low frequency decomposition part arranges from small to large by frequency, and HFS is by arranging from big to small.Accompanying drawing 2 is the frequency domain distribution rule of three layers of WAVELET PACKET DECOMPOSITION, visible, in WAVELET PACKET DECOMPOSITION, each layer frequency domain arranges and naturally arranges and there is inconsistent phenomenon, to rearrange frequency domain sequence when selected threshold, then adopt different Research on threshold selection selected threshold for different high, normal, basic frequency domains.For low frequency sequence this method thr=σ * [log (n)] with the formula 0.5/ (2n) 0.5for threshold model carries out choosing of threshold value, wherein n is signal length, and σ is the standard deviation of original signal WAVELET PACKET DECOMPOSITION coefficient.The WAVELET PACKET DECOMPOSITION coefficient of medium-high frequency sequence then adopts Stein without partial likelihood estimation principle self-adaptation selected threshold.
After trying to achieve threshold value, wavelet packet analysis majority all adopt hard threshold function (as | ω | during >T, η (ω)=ω, when | ω | during≤T, η (ω)=0) or soft-threshold function (when | ω | during >T, η (ω)=ω-sign (ω) * T, when | ω | during≤T, η (ω)=0) each decomposition node coefficient is processed.Hard thresholding method can be good at retaining the local features such as image border, but image there will be the vision distortion such as ring, pseudo-Gibbs' effect, and relative smooth is wanted in soft-threshold process, but may cause the distortion phenomenons such as edge fog.It shows in one-dimensional signal it is then that hard-threshold can make signal become coarse, and some point there will be interruption, and discontinuous point shrinks by soft-threshold, also can cause the contraction of signal simultaneously.And the soft threshold method of the improvement that this method proposes (when | ω | during >T, η (ω)=ω-sign (ω) * T+ sign (ω) * T/ (2*k+1), when | ω | during≤T, η (ω)=ω 2k+1/ [(2*k+1) * T 2k]) be compromise process to first two method, the large threshold function method of change along with Decomposition order is gradually near soft-threshold disposal route, then closer to hard-threshold disposal route when calculating the little node of Decomposition order, desirable denoising effect can be gathered in the crops like this while stick signal is complete.Accompanying drawing 3 is the function curve of three kinds of methods, and from figure, also can find out that hard thresholding method occurs being interrupted at threshold value T, soft threshold method makes curve produce contraction on η (ω) direction.
The invention has the beneficial effects as follows, it, for traditional denoising method, decreases calculated amount, improves the efficiency of program; The limitation of the obvious distortion of signal brought when overcoming hard-threshold or soft-threshold quantization parameter, while stick signal integrality, can obtain desirable denoised signal.
Accompanying drawing explanation
Fig. 1 is the signal denoising process flow diagram based on wavelet packet analysis.
Fig. 2 is the frequency domain distribution schematic diagram of three layers of WAVELET PACKET DECOMPOSITION.
Fig. 3 is the curve map of three kinds of threshold function functions, and wherein 3-1 is hard threshold function, and 3-2 is soft-threshold function, and 3-3 is the soft-threshold function of the improvement that this granting proposes.
Embodiment
Wavelet packet analysis has carried out meticulousr analysis on the basis of wavelet analysis, frequency band is carried out multilayer division by it, the high frequency series part that multi-frequency analysis is not segmented is decomposed further, and frequency band can be chosen according to the feature of analyzed signal is adaptive, make it to match with signal spectrum, thus when improving-frequency resolution, there is the prospect more wide than wavelet analysis.Effective wavelet packet analysis directly perceived is exactly carry out a threshold value quantizing to each node coefficient that WAVELET PACKET DECOMPOSITION obtains, and then to the coefficient reconstruct after quantification, thus reaches the object of denoising.
Signal antinoise method based on wavelet packet analysis will experience four basic steps: 1. the WAVELET PACKET DECOMPOSITION of signal, selects a small echo, then carries out WAVELET PACKET DECOMPOSITION to signals and associated noises; 2. choose cost function, conventional cost function is information entropy and logarithmic entropy; 3. the choosing of best wavelet packet basis, for a given entropy standard, adopts top-down searching algorithm to calculate optimal tree, chooses best wavelet packet basis; 4. selected threshold, adopts different choosing methods to obtain different threshold values according to different decomposition frequency domains; 5. the threshold value quantizing of WAVELET PACKET DECOMPOSITION coefficient, for the coefficient of each WAVELET PACKET DECOMPOSITION, adopts the soft-threshold function of the improvement threshold value corresponding with this coefficient to carry out quantification treatment; 6. the wavelet package reconstruction of signal, according to the WAVELET PACKET DECOMPOSITION coefficient of lowermost layer and the coefficient through quantification treatment, carries out wavelet package reconstruction.
In above-mentioned each step, it is crucial that how selected threshold and how to carry out threshold value quantizing, to a certain extent, he is directly connected to the quality of signal being carried out to noise reduction process, and accompanying drawing 1 carries out the concrete steps of signal denoising for the method.
Wavelet packet basis storehouse is made up of many wavelet packet basiss, and different wavelet packet basiss has different character, reacts different signal characteristics, and we wish to choose optimum wavelet packet basis to specific signal, is used for expressing the feature of this signal.From the viewpoint of signal transacting, the search procedure of optimal base is in fact with as far as possible few coefficient, reacts information as much as possible, to reach the object of characteristics extraction.Gap between coefficient is the bigger the better, if only have a few coefficients very large, so these coefficients just can the feature of representation signal, and contrary coefficients comparison is average, and gap not quite, is so just difficult to the feature finding out signal.
When searching for optimal base, first will provide the cost function of a burst, cost function can be defined as any function about burst, and find the base making cost function minimum in all wavelet packet basiss then in wavelet library, this base is just optimal base.In order to reduce calculated amount, improve program operation speed, reduce unknown parameter simultaneously, this method adopts the top-down Binomial Trees based on " beta pruning " to be that cost function carries out choosing of optimal base with logarithmic entropy, that is logarithmic entropy that is that only have the logarithm entropy of lower one deck and that be less than this layer is just decomposed, otherwise with regard to beta pruning, detailed algorithm flow is as follows: to signal node U n jcarry out WAVELET PACKET DECOMPOSITION, calculate each node logarithmic entropy; comparison node U n jentropy and its two child node U decomposed 2n j+1, U 2n+1 j+1entropy and, and mark the less node of entropy; decision node U n jwhether be labeled, if be labeled, then deletion of node U 2n j+1, U 2n+1 j+1decompose downwards no longer further, if be not labeled, then by node U 2n j+1, U 2n+1 j+1as node U n jrepeat above step.
Threshold denoising is that a kind of realization is simple, the good Wavelet Package Denoising Method of effect.Its basic thought be to wavelet decomposition after each layer coefficients in the mould coefficient that is greater than or less than certain threshold value process respectively, then carry out inverse transformation reconstruct denoising after signal.
The selection of threshold value is the key point of denoising, and threshold value is selected large, and the feature of useful signal will be filtered, if contrary Threshold selection is too small, and still noisy existence in the signal after denoising.Because the selection of threshold value directly affects de-noising effect, people propose an a lot of empirical model: soft-threshold, hard-threshold, minimax threshold value, fixed threshold etc., but all models are not general, have its scope of application.
In Practical Project signal, useful information mainly concentrates on low frequency part, and HFS majority is noise, only containing a small amount of useful information.If choose unified threshold value to process coefficient of dissociation, be difficult to reflect the different characteristic of signal at different frequency bands.Therefore need to choose multiple threshold value to carry out quantification treatment for different frequency bands.Natural order and the frequency band order of WAVELET PACKET DECOMPOSITION tree node there are differences, and want to choose different threshold value for different frequency bands, first will understand fully the frequency band order of WAVELET PACKET DECOMPOSITION tree node.When wavelet packet often decomposes one deck, low frequency decomposition part arranges from small to large by frequency, and HFS is by arranging from big to small.Accompanying drawing 2 is the frequency domain distribution rule of three layers of WAVELET PACKET DECOMPOSITION, visible, in WAVELET PACKET DECOMPOSITION, each layer frequency domain arranges and naturally arranges and there is inconsistent phenomenon, to rearrange frequency domain sequence when selected threshold, then adopt different Research on threshold selection selected threshold for different high, normal, basic frequency domains.For low frequency sequence this method thr=σ * [log (n)] with the formula 0.5/ (2n) 0.5for threshold model carries out choosing of threshold value, wherein n is signal length, and σ is the standard deviation of original signal WAVELET PACKET DECOMPOSITION coefficient.The WAVELET PACKET DECOMPOSITION coefficient of medium-high frequency sequence then adopts Stein without partial likelihood estimation principle self-adaptation selected threshold.
After trying to achieve threshold value, wavelet packet analysis majority all adopt hard threshold function (as | ω | during >T, η (ω)=ω, when | ω | during≤T, η (ω)=0) or soft-threshold function (when | ω | during >T, η (ω)=ω-sign (ω) * T, when | ω | during≤T, η (ω)=0) each decomposition node coefficient is processed.Hard thresholding method can be good at retaining the local features such as image border, but image there will be the vision distortion such as ring, pseudo-Gibbs' effect, and relative smooth is wanted in soft-threshold process, but may cause the distortion phenomenons such as edge fog.It shows in one-dimensional signal it is then that hard-threshold can make signal become coarse, and some point there will be interruption, and discontinuous point shrinks by soft-threshold, also can cause the contraction of signal simultaneously.And the soft threshold method of the improvement that this method proposes (when | ω | during >T, η (ω)=ω-sign (ω) * T+ sign (ω) * T/ (2*k+1), when | ω | during≤T, η (ω)=ω 2k+1/ [(2*k+1) * T 2k]) be compromise process to first two method, the large threshold function method of change along with Decomposition order is gradually near soft-threshold disposal route, then closer to hard-threshold disposal route when calculating the little node of Decomposition order, desirable denoising effect can be gathered in the crops like this while stick signal is complete.Accompanying drawing 3 is the function curve of three kinds of methods, and from figure, also can find out that hard thresholding method occurs being interrupted at threshold value T, soft threshold method makes curve produce contraction on η (ω) direction.

Claims (4)

1. one kind is carried out the method for signal denoising based on wavelet packet, wavelet packet analysis has carried out meticulousr analysis on the basis of wavelet analysis, frequency band is carried out multilayer division by it, the high frequency series part that multi-frequency analysis is not segmented is decomposed further, and frequency band can be chosen according to the feature of analyzed signal is adaptive, make it to match with signal spectrum, thus when improving-frequency resolution, it is characterized in that, concrete steps are:
1) WAVELET PACKET DECOMPOSITION of signal, selects a small echo, then carries out WAVELET PACKET DECOMPOSITION to signals and associated noises;
2) information entropy or logarithmic entropy is chosen as cost function;
3) the choosing of best wavelet packet basis, for a given entropy standard, adopts top-down searching algorithm to calculate optimal tree, chooses best wavelet packet basis;
4) selected threshold, adopts different choosing methods to obtain different threshold values according to different decomposition frequency domains;
5) threshold value quantizing of WAVELET PACKET DECOMPOSITION coefficient, for the coefficient of each WAVELET PACKET DECOMPOSITION, adopts the soft-threshold function of the improvement threshold value corresponding with this coefficient to carry out quantification treatment;
6) wavelet package reconstruction of signal, according to the WAVELET PACKET DECOMPOSITION coefficient of lowermost layer and the coefficient through quantification treatment, carries out wavelet package reconstruction.
2. method of carrying out signal denoising based on wavelet packet according to claim 1, it is characterized in that, the concrete steps chosen of best wavelet packet basis described in step 3) are: the cost function first providing a burst, then find the base making cost function minimum in all wavelet packet basiss in wavelet library, this base is just optimal base; The top-down Binomial Trees based on " beta pruning " is adopted to be that cost function carries out choosing of optimal base with logarithmic entropy, that is logarithmic entropy that is that only have the logarithm entropy of lower one deck and that be less than this layer is just decomposed, otherwise with regard to beta pruning, detailed algorithm flow is as follows: to signal node U n jcarry out WAVELET PACKET DECOMPOSITION, calculate each node logarithmic entropy; comparison node U n jentropy and its two child node U decomposed 2n j+1, U 2n+1 j+1entropy and, and mark the less node of entropy; decision node U n jwhether be labeled, if be labeled, then deletion of node U 2n j+1, U 2n+1 j+1decompose downwards no longer further, if be not labeled, then by node U 2n j+1, U 2n+1 j+1as node U n jrepeat above step.
3. method of carrying out signal denoising based on wavelet packet according to claim 1, it is characterized in that, step 4) is specially: rearrange frequency domain sequence, different Research on threshold selection selected threshold is adopted, for low frequency sequence thr=σ * [log (n)] with the formula for different high, normal, basic frequency domains 0.5/ (2n) 0.5for threshold model carries out choosing of threshold value, wherein n is signal length, and σ is the standard deviation of original signal WAVELET PACKET DECOMPOSITION coefficient, and the WAVELET PACKET DECOMPOSITION coefficient of medium-high frequency sequence then adopts Stein without partial likelihood estimation principle self-adaptation selected threshold.
4. method of carrying out signal denoising based on wavelet packet according to claim 1, it is characterized in that, the soft-threshold function of the improvement described in step 5) is: when | ω | during >T, η (ω)=ω-sign (ω) * T+ sign (ω) * T/ (2*k+1), when | ω | during≤T, η (ω)=ω 2k+1/ [(2*k+1) * T 2k]; The large threshold function method of change along with Decomposition order, gradually near soft-threshold disposal route, then closer to hard-threshold disposal route when calculating the little node of Decomposition order, can gather in the crops desirable denoising effect like this while stick signal is complete.
CN201510042993.3A 2015-01-28 2015-01-28 Method for conducting signal denoising based on wavelet packet Pending CN104615877A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510042993.3A CN104615877A (en) 2015-01-28 2015-01-28 Method for conducting signal denoising based on wavelet packet

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510042993.3A CN104615877A (en) 2015-01-28 2015-01-28 Method for conducting signal denoising based on wavelet packet

Publications (1)

Publication Number Publication Date
CN104615877A true CN104615877A (en) 2015-05-13

Family

ID=53150318

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510042993.3A Pending CN104615877A (en) 2015-01-28 2015-01-28 Method for conducting signal denoising based on wavelet packet

Country Status (1)

Country Link
CN (1) CN104615877A (en)

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105719658A (en) * 2016-03-04 2016-06-29 安徽师范大学 Wavelet packet speech denoising method based on new threshold function and self-adaptive threshold
CN106228520A (en) * 2016-07-22 2016-12-14 上海联影医疗科技有限公司 Image enchancing method and device
CN106330343A (en) * 2015-07-07 2017-01-11 大连测控技术研究所 Acoustic beacon signal extraction method based on wavelet packet and fuzzy theory
CN106437687A (en) * 2016-08-31 2017-02-22 中国船舶重工集团公司第七八研究所 While-drilling mud impulse signal noise reduction method based on optimal wavelet packet decomposition tree
CN106782588A (en) * 2016-11-15 2017-05-31 中国人民解放军第四军医大学 A kind of millimeter wave bioradar sound enhancement method based on sense of hearing wavelet packet
CN107273585A (en) * 2017-05-25 2017-10-20 国网山东省电力公司青岛供电公司 A kind of load ratio bridging switch fault detection method and device
CN107659314A (en) * 2017-09-19 2018-02-02 电子科技大学 The rarefaction expression of distributing optical fiber sensing space-time two-dimension signal and compression method
CN107688120A (en) * 2016-11-17 2018-02-13 丹阳华神电器有限公司 Signals and associated noises processing method and iteration singular spectrum Soft-threshold Denoising Method based on fuzzy entropy
CN108491355A (en) * 2018-02-05 2018-09-04 南京邮电大学 A kind of ultrasonic signal noise-reduction method based on CEEMD and wavelet packet
CN108647614A (en) * 2018-04-28 2018-10-12 吉林大学 The recognition methods of electrocardiogram beat classification and system
CN108873813A (en) * 2018-06-25 2018-11-23 山东大学 Tool wear degree detection method based on main shaft of numerical control machine tool servo motor current signal
CN108985234A (en) * 2018-07-19 2018-12-11 沈阳建筑大学 A kind of Bayes's wavelet-packet noise reduction method suitable for non-Gaussian signal
CN109561437A (en) * 2017-09-26 2019-04-02 大唐移动通信设备有限公司 The method and device that a kind of pair of signal is handled
US10290108B2 (en) 2015-12-31 2019-05-14 Shanghai United Imaging Healthcare Co., Ltd. Methods and systems for image processing
CN110096998A (en) * 2019-04-29 2019-08-06 天津大学 A kind of decomposition method extracted for machining surface shape characteristic
CN110691229A (en) * 2019-08-23 2020-01-14 昆明理工大学 Hologram compression method, encoder and reproduced image output system
CN111079591A (en) * 2019-12-04 2020-04-28 国网天津市电力公司电力科学研究院 Bad data restoration method and system based on improved multi-scale principal component analysis
CN112613705A (en) * 2020-12-14 2021-04-06 中广核研究院有限公司 Component quality acquisition method, component quality acquisition device, computer equipment and storage medium
CN115389888A (en) * 2022-10-28 2022-11-25 山东科华电力技术有限公司 Partial discharge real-time monitoring system based on high-voltage cable
WO2024099006A1 (en) * 2023-07-12 2024-05-16 广西科学院 Beidou offshore multipath interference elimination method and apparatus based on wavelet packet analysis

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005067640A2 (en) * 2004-01-07 2005-07-28 University Of Hawaii Methods for enhanced detection & analysis of differentially expressed genes using gene chip microarrays
CN102156963A (en) * 2011-01-20 2011-08-17 中山大学 Denoising method for image with mixed noises
CN102508206A (en) * 2011-10-31 2012-06-20 哈尔滨工程大学 Linear frequency modulation (FM) signal parameter estimation method based on small-wave-packet denoising and power spectral entropy
CN103901474A (en) * 2014-04-18 2014-07-02 成都新核中创信息科技有限公司 Microseism weak signal extracting method based on middle course wavelet threshold analysis

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005067640A2 (en) * 2004-01-07 2005-07-28 University Of Hawaii Methods for enhanced detection & analysis of differentially expressed genes using gene chip microarrays
CN102156963A (en) * 2011-01-20 2011-08-17 中山大学 Denoising method for image with mixed noises
CN102508206A (en) * 2011-10-31 2012-06-20 哈尔滨工程大学 Linear frequency modulation (FM) signal parameter estimation method based on small-wave-packet denoising and power spectral entropy
CN103901474A (en) * 2014-04-18 2014-07-02 成都新核中创信息科技有限公司 Microseism weak signal extracting method based on middle course wavelet threshold analysis

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
侯新国 等: "最佳小波包基改进软阈值的消噪方法及应用", 《振动、测试与诊断》 *
章浙涛 等: "小波包多阈值去噪法及其在形变分析中的应用", 《测绘学报》 *
罗永 等: "小波包变换及代价函数设计综述", 《数学理论与应用》 *
郭晓霞 等: "小波去噪中软硬阈值的一种改良折衷法", 《智能系统学报》 *

Cited By (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106330343A (en) * 2015-07-07 2017-01-11 大连测控技术研究所 Acoustic beacon signal extraction method based on wavelet packet and fuzzy theory
CN106330343B (en) * 2015-07-07 2020-04-17 大连测控技术研究所 Underwater acoustic beacon signal extraction method based on wavelet packet and fuzzy theory
US10290108B2 (en) 2015-12-31 2019-05-14 Shanghai United Imaging Healthcare Co., Ltd. Methods and systems for image processing
US11880978B2 (en) 2015-12-31 2024-01-23 Shanghai United Imaging Healthcare Co., Ltd. Methods and systems for image processing
US11049254B2 (en) 2015-12-31 2021-06-29 Shanghai United Imaging Healthcare Co., Ltd. Methods and systems for image processing
CN105719658B (en) * 2016-03-04 2019-11-26 安徽师范大学 Wavelet packet voice de-noising method based on new threshold function table and adaptive threshold
CN105719658A (en) * 2016-03-04 2016-06-29 安徽师范大学 Wavelet packet speech denoising method based on new threshold function and self-adaptive threshold
CN106228520B (en) * 2016-07-22 2019-10-22 上海联影医疗科技有限公司 Image enchancing method and device
CN106228520A (en) * 2016-07-22 2016-12-14 上海联影医疗科技有限公司 Image enchancing method and device
CN106437687A (en) * 2016-08-31 2017-02-22 中国船舶重工集团公司第七八研究所 While-drilling mud impulse signal noise reduction method based on optimal wavelet packet decomposition tree
CN106782588A (en) * 2016-11-15 2017-05-31 中国人民解放军第四军医大学 A kind of millimeter wave bioradar sound enhancement method based on sense of hearing wavelet packet
CN107688120A (en) * 2016-11-17 2018-02-13 丹阳华神电器有限公司 Signals and associated noises processing method and iteration singular spectrum Soft-threshold Denoising Method based on fuzzy entropy
WO2018090467A1 (en) * 2016-11-17 2018-05-24 丹阳华神电器有限公司 Fuzzy entropy based noisy signal processing method and iterative singular spectrum soft thresholding denoising method
CN107273585A (en) * 2017-05-25 2017-10-20 国网山东省电力公司青岛供电公司 A kind of load ratio bridging switch fault detection method and device
CN107273585B (en) * 2017-05-25 2021-04-13 国网山东省电力公司青岛供电公司 On-load tap-changer fault detection method and device
CN107659314A (en) * 2017-09-19 2018-02-02 电子科技大学 The rarefaction expression of distributing optical fiber sensing space-time two-dimension signal and compression method
CN107659314B (en) * 2017-09-19 2021-02-19 电子科技大学 Sparse representation and compression method of distributed optical fiber sensing space-time two-dimensional signal
CN109561437B (en) * 2017-09-26 2020-10-16 大唐移动通信设备有限公司 Method and device for processing signal
CN109561437A (en) * 2017-09-26 2019-04-02 大唐移动通信设备有限公司 The method and device that a kind of pair of signal is handled
WO2019062406A1 (en) * 2017-09-26 2019-04-04 大唐移动通信设备有限公司 Method and apparatus for signal processing
US11283661B2 (en) 2017-09-26 2022-03-22 Datang Mobile Communications Equipment Co., Ltd. Method and apparatus for signal processing
CN108491355A (en) * 2018-02-05 2018-09-04 南京邮电大学 A kind of ultrasonic signal noise-reduction method based on CEEMD and wavelet packet
CN108647614A (en) * 2018-04-28 2018-10-12 吉林大学 The recognition methods of electrocardiogram beat classification and system
CN108873813B (en) * 2018-06-25 2020-04-28 山东大学 Cutter abrasion degree detection method based on numerical control machine tool spindle servo motor current signal
CN108873813A (en) * 2018-06-25 2018-11-23 山东大学 Tool wear degree detection method based on main shaft of numerical control machine tool servo motor current signal
CN108985234A (en) * 2018-07-19 2018-12-11 沈阳建筑大学 A kind of Bayes's wavelet-packet noise reduction method suitable for non-Gaussian signal
CN108985234B (en) * 2018-07-19 2021-08-31 沈阳建筑大学 Bayes wavelet packet noise reduction method suitable for non-Gaussian signals
CN110096998B (en) * 2019-04-29 2023-05-12 天津大学 Decomposition method for extracting machined surface morphology features
CN110096998A (en) * 2019-04-29 2019-08-06 天津大学 A kind of decomposition method extracted for machining surface shape characteristic
CN110691229B (en) * 2019-08-23 2021-10-22 昆明理工大学 Hologram compression method, encoder and reproduced image output system
CN110691229A (en) * 2019-08-23 2020-01-14 昆明理工大学 Hologram compression method, encoder and reproduced image output system
CN111079591B (en) * 2019-12-04 2024-01-02 国网天津市电力公司电力科学研究院 Bad data restoration method and system based on improved multi-scale principal component analysis
CN111079591A (en) * 2019-12-04 2020-04-28 国网天津市电力公司电力科学研究院 Bad data restoration method and system based on improved multi-scale principal component analysis
CN112613705A (en) * 2020-12-14 2021-04-06 中广核研究院有限公司 Component quality acquisition method, component quality acquisition device, computer equipment and storage medium
CN115389888A (en) * 2022-10-28 2022-11-25 山东科华电力技术有限公司 Partial discharge real-time monitoring system based on high-voltage cable
CN115389888B (en) * 2022-10-28 2023-01-31 山东科华电力技术有限公司 Partial discharge real-time monitoring system based on high-voltage cable
WO2024099006A1 (en) * 2023-07-12 2024-05-16 广西科学院 Beidou offshore multipath interference elimination method and apparatus based on wavelet packet analysis

Similar Documents

Publication Publication Date Title
CN104615877A (en) Method for conducting signal denoising based on wavelet packet
CN102156963A (en) Denoising method for image with mixed noises
CN103761719A (en) Self-adaptive wavelet threshold de-noising method based on neighborhood correlation
CN103295204B (en) A kind of image self-adapting enhancement method based on non-down sampling contourlet transform
CN106570843A (en) Adaptive wavelet threshold function image noise suppression method
CN111583152A (en) Image artifact detection and automatic removal method based on U-net structure
CN107909558A (en) A kind of non-local mean image de-noising method based on unsupervised learning
CN110060220A (en) Based on the image de-noising method and system for improving BM3D algorithm
CN109978775B (en) Color denoising method and device
CN104182944A (en) Optical image denoising method based on serial connection of curvelet transform and wavelet transform
CN108269246B (en) Image equalization enhancement method for low-frequency wavelet coefficient interpolation
CN103440870A (en) Method and device for voice frequency noise reduction
Qin et al. An improved method of image denoising based on wavelet transform
CN105719658A (en) Wavelet packet speech denoising method based on new threshold function and self-adaptive threshold
Huang et al. Synthetic aperture radar image compression using tree-structured edge-directed orthogonal wavelet packet transform
Liu et al. Restoration algorithm for noisy complex illumination
CN105259539B (en) A kind of branch's formula noise suppressing method that ice radar data is visited for polar region
Thote et al. Improved denoising technique for natural and synthetic images
CN108664934B (en) Deep space communication image optimal wavelet packet basis searching method
Radhika et al. Image denoising techniques preserving edge
Ahmad et al. Rain Streaks Elimination Using Hybrid Median Filter and Contrast Stretching
CN104422961A (en) Random noise attenuation method for earthquakes
CN115866151B (en) Image communication method
Chen et al. An improved wavelet-based method for SAR images denoising using data fusion technique
CN105631822B (en) A method of improving small echo variable threshold values contracting model reduces picture noise

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20150513

WD01 Invention patent application deemed withdrawn after publication