CN110430340A - A kind of noise-reduction method and system of pulse array signals - Google Patents

A kind of noise-reduction method and system of pulse array signals Download PDF

Info

Publication number
CN110430340A
CN110430340A CN201910577862.3A CN201910577862A CN110430340A CN 110430340 A CN110430340 A CN 110430340A CN 201910577862 A CN201910577862 A CN 201910577862A CN 110430340 A CN110430340 A CN 110430340A
Authority
CN
China
Prior art keywords
pulse
search
array signals
region
neighborhood
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910577862.3A
Other languages
Chinese (zh)
Other versions
CN110430340B (en
Inventor
田永鸿
朱林
李家宁
付溢华
董思维
黄铁军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Pulse vision (Beijing) Technology Co.,Ltd.
Original Assignee
Peking 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 Peking University filed Critical Peking University
Priority to CN201910577862.3A priority Critical patent/CN110430340B/en
Publication of CN110430340A publication Critical patent/CN110430340A/en
Application granted granted Critical
Publication of CN110430340B publication Critical patent/CN110430340B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N5/00Details of television systems
    • H04N5/14Picture signal circuitry for video frequency region
    • H04N5/21Circuitry for suppressing or minimising disturbance, e.g. moiré or halo
    • H04N5/213Circuitry for suppressing or minimising impulsive noise

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Image Processing (AREA)

Abstract

The invention discloses a kind of noise-reduction methods of pulse array signals, it include: that the spatial neighborhood for treating denoising pulse and the granting characteristics of the pulse array signals on timestamp neighborhood are modeled, the granting characteristic model of pulse array signals is established, pulse array signals spatial-temporal distribution characteristic is extracted;According to the noise profile characteristic in pulse array signals, space time filter is constructed, the feature of extraction is filtered;According to the granting characteristic model, filtered pulse spatial-temporal distribution characteristic is reduced to pulse array signals.The present invention has fully considered spatial-temporal distribution characteristic of the noise in pulse array signals, randomness on its transitivity in the time domain, airspace is modeled, in conjunction with the physical significance that pulse itself represents, space time filter is constructed and pulse signal is efficiently denoised, it is high-efficient, excellent noise reduction effect.

Description

A kind of noise-reduction method and system of pulse array signals
Technical field
The present invention relates to signal processing technology fields, and in particular to a kind of noise-reduction method and system of pulse array signals.
Background technique
Traditional visual sensor carries out complete sample to scene usually as unit of frame, according to preset fixed frequency. This sampling based on fixed frame per second, cannot reflect the dynamic change of scene, be easy to lead current scene over-sampling or lack sampling Cause the problems such as video data redundancy is big, time resolution is low, easily fuzzy under high-speed motion.It is inspired in the vision of biological retina Sampling mechanism, the Novel camera of acquisition pulse array signal enters the people visual field, including provides arteries and veins based on intensity of illumination variation The sensor for rushing signal, such as dynamic visual sensor (Dynamic Vision Sensor, DVS), the image based on asynchronous time Sensor (Asynchronous Time-based image Sensor, ATIS), dynamic active pixel visual sensor (Dynamic and Active Pixel Vision Sensor, DAVIS) etc. provides signal based on intensity of illumination integrated intensity Sensor, such as light intensity accumulate sensor.The sensor of this video camera acquires certain time, optical signal in certain area Information has many advantages, such as high dynamic range, high time resolution.
Noise is an important problem in field of signal processing.Pulse array signals have high temporal resolution, Temporal correlation between signal is stronger, and influence of noise has the characteristics that the transitivity in time domain, the randomness on airspace.With biography System image and signal denoising difference, denoise pulse array signals, need the characteristics of fully considering pulse signal and physics Meaning efficiently handles a large amount of pulse data.
Summary of the invention
One purpose of the disclosure is to provide a kind of new technical solution of the noise reduction of pulse array signals.In order to disclosure The some aspects of embodiment have a basic understanding, simple summary is shown below.The summarized section is not generally to comment It states, nor to determine key/critical component or describe the protection scope of these embodiments.Its sole purpose is with simple Form some concepts are presented, in this, as the preamble of following detailed description.
According to a first aspect of the embodiments of the present invention, a kind of noise-reduction method of pulse array signals is provided, comprising: treat It makes an uproar the spatial neighborhood of pulse and the granting characteristic of the pulse array signals on timestamp neighborhood is modeled, establish pulse array letter Number granting characteristic model, extract pulse array signals spatial-temporal distribution characteristic;
According to the noise profile characteristic in pulse array signals, space time filter is constructed, the feature of extraction is filtered;
According to the granting characteristic model, filtered pulse spatial-temporal distribution characteristic is reduced to pulse array signals.
Further, the pulse to be denoised is specific pulse on time-space domain, is in a certain specific pixel, from a certain Impulse response sending starts, one section of pulse signal before occurring to next pulse response.
Further, the noise profile characteristic according in pulse array signals constructs space time filter, to extraction Feature be filtered, comprising:
Pulse array signals of the pulse to be denoised on spatial neighborhood and timestamp neighborhood are set as its first field of search Domain;
Pulse array signals in the first region of search are analyzed, the pulse release information of several pulse trains are obtained, according to arteries and veins Release information is rushed to obtain each pulse expression feature or convert transform domain extraction feature for pulse release information;
The pulse array signals that first region of search includes are converted to pulse strength, on pulse strength domain, are calculated special Determine the similarity of other pulses in pulse and the first region of search, set first threshold, retains similarity in the first region of search Greater than the pulse of first threshold, it is denoted as the second region of search;
For the second region of search, second threshold is set, judges whether the number of pulses in the second region of search is greater than the Two threshold values;If so, being estimated using the pulse pair pulse to be denoised in the second region of search, and it is replaced using estimated value Former intensity;If it is not, then subtracting the second region of search using the first region of search, third region of search is obtained, is searched for using third Pulse pair pulse to be denoised in region is estimated, and replaces its former intensity using estimated value.
Further, obtaining each pulse according to pulse release information indicates feature, are as follows: according to pulse Time Of Release interval It determines pulse strength, or pulse strength is determined according to the inverse at pulse Time Of Release interval, or according to pulse Time Of Release interval Logarithm determine pulse strength.
Further, the similarity of the calculating certain pulses and other pulses in the first region of search, comprising: calculate The variance of pulse strength in region of search determines a confidence interval according to the intensity size of variance and pulse to be denoised, according to Distribution situation of the pulse in confidence interval determines similarity.
Further, the first threshold, second threshold are dynamically set according to neighborhood pulse relationship, to adapt to difference Pulse data.
Further, the first threshold is different from the second threshold.
Further, the pulse pair pulse to be denoised in second region of search of use is estimated, including but not It is limited to: weighted mean method, Neighborhood Filtering method.
Further, the denoising pulse for the treatment of is estimated as being estimated by filter, space time filter Weight or convolution kernel assign bigger value to the pulse closer to pulse to be denoised.
Further, according to the granting characteristic model, filtered pulse spatial-temporal distribution characteristic is reduced to pulse battle array Column signal, comprising: the intensity value of certain pulses adjusted is reverted into pulse array signals according to pulse sequence relationship, Method is the inverse operation that each pulse expression feature this operation step is obtained according to pulse release information.
According to a second aspect of the embodiments of the present invention, a kind of noise reduction system of pulse array signals is provided, which is characterized in that Include:
Extraction module is modeled, the pulse array signals on spatial neighborhood and timestamp neighborhood for treating denoising pulse It provides characteristic to be modeled, establishes the granting characteristic model of pulse array signals, extract pulse array signals spatial-temporal distribution characteristic;
Filter module, for space time filter being constructed, to extraction according to the noise profile characteristic in pulse array signals Feature is filtered;
Recovery module, for according to the granting characteristic model, filtered pulse spatial-temporal distribution characteristic to be reduced to arteries and veins Rush array signal.
According to a third aspect of the embodiments of the present invention, a kind of noise reduction system of pulse array signals is provided, for realizing upper The method stated, the system include: sequentially connected neighborhood search module, pulse analysis and conversion module, similarity estimation module With pulse denoising and recombination module;
The neighborhood search module obtains target search region for searching for the neighborhood pulse array data of target pulse;
The pulse analysis and conversion module obtain each arteries and veins for corresponding to each pixel search region for analyzing pulse data The pulse release information of sequence is rushed, and is converted into pulse strength;
The similarity estimation module, for calculating the intensity similarity of each pulse and target pulse in region of search;
The pulse denoising and recombination module, for being denoised to target pulse, and the pulse data weight after denoising Group is pulse array signals.
According to a fourth aspect of the embodiments of the present invention, a kind of electronic equipment is provided, including memory, processor and is stored in On the memory and the computer program that can run on the processor, the processor execute described program, to realize Above-mentioned method.
According to a fifth aspect of the embodiments of the present invention, a kind of non-transitorycomputer readable storage medium is provided, is deposited thereon Computer program is contained, which is executed by processor, to realize above-mentioned method.
Technical solution provided in an embodiment of the present invention can include the following benefits:
The noise-reduction method and system of pulse array signals provided by the invention, for specific arteries and veins to be denoised on time-space domain A certain range of pulse array signals on its spatial neighborhood and timestamp neighborhood, are converted into pulse strength information by punching;Meter The similarity for calculating pulse to be denoised Yu neighborhood pulse is filtered using the intensity of neighborhood pulse pair pulse to be denoised;For adjusting The intensity value of certain pulses after whole reverts to pulse array signals according to its pulse sequence relationship;The present invention fully considers Spatial-temporal distribution characteristic of the noise in pulse array signals builds the randomness on its transitivity in the time domain, airspace Mould constructs space time filter and is efficiently denoised to pulse signal in conjunction with the physical significance that pulse itself represents, high-efficient, drop Effect of making an uproar is good.
Other features and advantages of the present invention will illustrate in the following description, also, partly become from specification It is clear that alternatively, Partial Feature and advantage can deduce from specification or unambiguously determine, or pass through implementation The embodiment of the present invention understands.The objectives and other advantages of the invention can be by written specification, claims and attached Specifically noted structure is achieved and obtained in figure.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this The some embodiments recorded in invention, for those of ordinary skill in the art, without creative efforts, It is also possible to obtain other drawings based on these drawings.
Fig. 1 is the flow chart of the noise-reduction method of the pulse array signals of one embodiment of the disclosure;
Fig. 2 is a kind of step schematic diagram of the noise-reduction method for pulse array signals that another embodiment of the disclosure provides;
Fig. 3 is that a kind of space-time neighborhood selection of the noise-reduction method for pulse array signals that another embodiment of the disclosure provides is shown It is intended to;
Fig. 4 is pulse signal and the pulse of a kind of noise-reduction method for pulse array signals that another embodiment of the disclosure provides Strength relationship exemplary diagram;
Fig. 5 is a kind of structural block diagram of the noise reduction system for pulse array signals that another embodiment of the disclosure provides.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawing and specific implementation The present invention will be further described for example.It should be appreciated that described herein, specific examples are only used to explain the present invention, and does not have to It is of the invention in limiting.Based on the embodiments of the present invention, those of ordinary skill in the art are not making creative work premise Under every other embodiment obtained, shall fall within the protection scope of the present invention.
Those skilled in the art of the present technique are appreciated that unless otherwise defined, all terms used herein (including technology art Language and scientific term), there is meaning identical with the general understanding of those of ordinary skill in fields of the present invention.Should also Understand, those terms such as defined in the general dictionary, it should be understood that have in the context of the prior art The consistent meaning of meaning, and unless idealization or meaning too formal otherwise will not be used by specific definitions as here To explain.
As shown in Figure 1, one embodiment of the disclosure, provides a kind of noise-reduction method of pulse array signals, comprising:
Step S1: treat denoising pulse spatial neighborhood and the pulse array signals on timestamp neighborhood granting characteristic into Row modeling, establishes the granting characteristic model of pulse array signals, extracts pulse array signals spatial-temporal distribution characteristic.
The pulse to be denoised is specific pulse on time-space domain;Specific pulse, is defined as on the time-space domain, One section of pulse signal in a certain specific pixel, since issuing a certain impulse response, before occurring to next pulse response.
The spatial neighborhood and timestamp neighborhood, spatial neighborhood are four neighborhoods, eight neighborhood or other relationships, and timestamp is adjacent Domain is that can use the bigger radius of neighbourhood of space neighborhood as unit of pulse.
Step S2: according to the noise profile characteristic in pulse array signals, constructing space time filter, to the feature of extraction into Row filtering;Step S2 is specifically included:
Step S21: for specifically pulse to be denoised on time-space domain, it is set on spatial neighborhood and timestamp neighborhood A certain range of pulse array signals are as its first region of search;The side of the desirable smaller such as 2*2 or 4*4 of its spatial neighborhood The circle that block or radius are 1.5 or 2, is also possible to biggish neighborhood;Its timestamp neighborhood is generally bigger than spatial neighborhood range, with Adapt to the high time resolution of pulse array signals;
Step S22: letter is provided in pulse array signals in the first region of search of analysis, the pulse for obtaining several pulse trains Breath obtains each pulse expression feature or converts transform domain for pulse release information and extracts feature according to pulse release information;
The method that each pulse indicates feature is obtained according to pulse release information, it can be with are as follows: according between pulse Time Of Release Pulse strength is determined every determining pulse strength, or according to the inverse at pulse Time Of Release interval, or according between pulse Time Of Release Every logarithm determine pulse strength;Timestamp neighborhood is that can use the bigger radius of neighbourhood of space neighborhood as unit of pulse;
Step S23: for the first region of search, it includes pulse array signals be converted to pulse strength, in pulse In intensity domain, the similarity of other pulses in certain pulses and the first region of search is calculated, sets first threshold, retains first and searches Similarity is greater than the pulse of first threshold in rope region, is denoted as the second region of search;
Step S24: for the second region of search, second threshold is set, judges that the number of pulses in the second region of search is It is no to be greater than second threshold;
If the number of pulses in the second region of search is greater than second threshold, the pulse pair in the second region of search is used Pulse to be denoised is estimated, and replaces its former intensity using estimated value;
If the number of pulses in the second region of search is less than second threshold, second is subtracted using the first region of search and is searched Rope region obtains third region of search, is estimated using the pulse pair pulse to be denoised in third region of search, and use is estimated Evaluation replaces its former intensity.
The similarity of other pulses in the calculating certain pulses and the first region of search, comprising: calculate region of search The variance of interior pulse strength determines a confidence interval according to the intensity size of variance and pulse to be denoised, is being set according to pulse Distribution situation in letter section determines similarity.
The first threshold, second threshold are dynamically set according to neighborhood pulse relationship, to adapt to different pulse number evidence. The first threshold is different from the second threshold.
Pulse pair pulse to be denoised in second region of search of use is estimated that including but not limited to weighting is flat Equal method, Neighborhood Filtering method etc..
It is described to treat denoising pulse and be estimated as being estimated by filter, the weight or convolution kernel of filter Bigger value is assigned to the pulse closer to pulse to be denoised.
Step S3: according to the granting characteristic model for the pulse array signals established in step S1, when filtered pulse Empty distribution characteristics is reduced to pulse array signals;That is: by the intensity value of certain pulses adjusted, according to pulse sequence relationship, Revert to pulse array signals.The intensity value by certain pulses adjusted reverts to arteries and veins according to pulse sequence relationship Array signal is rushed, method is to obtain each pulse according to pulse release information described in step S22 to indicate feature this operation The inverse operation of step.
A kind of noise reduction system of pulse array signals, comprising:
Extraction module is modeled, the pulse array signals on spatial neighborhood and timestamp neighborhood for treating denoising pulse It provides characteristic to be modeled, establishes the granting characteristic model of pulse array signals, extract pulse array signals spatial-temporal distribution characteristic;
Filter module, for space time filter being constructed, to extraction according to the noise profile characteristic in pulse array signals Feature is filtered;
Recovery module, for according to the granting characteristic model, filtered pulse spatial-temporal distribution characteristic to be reduced to arteries and veins Rush array signal.
A kind of electronic equipment, including memory, processor and be stored on the memory and can be on the processor The computer program of operation, the processor executes described program, to realize noise-reduction method described above.
A kind of non-transitorycomputer readable storage medium, is stored thereon with computer program, which is held by processor Row, to realize noise-reduction method described above.
According to presently filed embodiment, a kind of noise-reduction method of pulse array signals is proposed, as shown in Figure 2, comprising:
A certain range of pulse array letter according to pulse position to be denoised, on its spatial neighborhood and timestamp neighborhood Number it is used as its first region of search, as shown in figure 3, selecting eight neighborhood or radius for 1.5 circle is spatial neighborhood SpatialRange;Assuming that it is t that pulse pair to be denoised, which answers time interval, select 5*t as timestamp neighborhood TimeRange;Really Fixed first region of search.
Pulse array signals in first region of search are converted to pulse strength, as shown in Figure 4;
Calculate all pulse strengties in the first region of search and and quadratic sum, and calculate mean value;
Calculate the whole variance and standard deviation of pulse in the first region of search;
Assuming that number of pulses is M, the standard deviation sigma of pulse strength information, scale coefficient k, setting the in the first region of search One threshold value T1=k* σ.Assuming that the intensity of pulse to be denoised is P, confidence interval bound is calculated according to average value, confidence area is set Between be [p-k* σ, p+k* σ].
The pulse in the first region of search is traversed again, counts the number of pulses in confidence interval.
Be arranged second threshold T2=M/10, if meet condition umber of pulse be less than T2, by current PRF and it is all Pulse in confidence interval excludes, and takes other interior pulse strength average values of the first region of search S1 as the result P ' after denoising:
Wherein, num (S) is the quantity of pulse in the S of region.The intersection of a ∩ b expression region a and region b.It indicates not Belong to the region of S.
If the umber of pulse for meeting condition is greater than T2, all pulses for meeting condition in the first region of search are taken, are remembered For the second region of search, the average intensity value of pulse in the second region of search S2 is calculated, the denoising result P ' as target pulse:
It should be noted that the mode for solving final result can also be uneven by setting according to the distance of distance objective pulse Even weight can obtain better effect in some cases.
After the completion of current PRF denoising, the pulse of next position is denoised using the above method, until traversing pulse battle array All positions in column signal.
As shown in figure 5, a kind of noise reduction system of pulse array signals, for realizing noise-reduction method described in the present embodiment, The system includes: sequentially connected neighborhood search module, pulse analysis and conversion module, similarity estimation module and pulse denoising And recombination module;
The neighborhood search module obtains target search region for searching for the neighborhood pulse array data of target pulse;
The pulse analysis and conversion module obtain each arteries and veins for corresponding to each pixel search region for analyzing pulse data The pulse release information of sequence is rushed, and is converted into pulse strength;
The similarity estimation module, for calculating the intensity similarity of each pulse and target pulse in region of search;
The pulse denoising and recombination module, for being denoised to target pulse, and the pulse data weight after denoising Group is pulse array signals.
Term " module " is not intended to be limited to specific physical form.Depending on concrete application, module be can be implemented as firmly Part, firmware, software and/or combination thereof.In addition, different modules can share common component or even be realized by same components. May exist between disparate modules or there is no clear boundaries.
It should be understood that although each step in the flow chart of attached drawing is successively shown according to the instruction of arrow, These steps are not that the inevitable sequence according to arrow instruction successively executes.Unless expressly stating otherwise herein, these steps Execution there is no stringent sequences to limit, can execute in the other order.Moreover, at least one in the flow chart of attached drawing Part steps may include that perhaps these sub-steps of multiple stages or stage are not necessarily in synchronization to multiple sub-steps Completion is executed, but can be executed at different times, execution sequence, which is also not necessarily, successively to be carried out, but can be with other At least part of the sub-step or stage of step or other steps executes in turn or alternately.
Algorithm and display be not inherently related to any certain computer, virtual bench or other equipment provided herein. Various fexible units can also be used together with teachings based herein.As described above, it constructs required by this kind of device Structure be obvious.In addition, the present invention is also not directed to any particular programming language.It should be understood that can use various Programming language realizes summary of the invention described herein, and the description done above to language-specific is to disclose this hair Bright preferred forms.
In the instructions provided here, numerous specific details are set forth.It is to be appreciated, however, that implementation of the invention Example can be practiced without these specific details.In some instances, well known method, structure is not been shown in detail And technology, so as not to obscure the understanding of this specification.
Similarly, it should be understood that in order to simplify the disclosure and help to understand one or more of the various inventive aspects, In Above in the description of exemplary embodiment of the present invention, each feature of the invention is grouped together into single implementation sometimes In example, figure or descriptions thereof.However, the disclosed method should not be interpreted as reflecting the following intention: i.e. required to protect Shield the present invention claims features more more than feature expressly recited in each claim.More precisely, as following Claims reflect as, inventive aspect is all features less than single embodiment disclosed above.Therefore, Thus the claims for following specific embodiment are expressly incorporated in the specific embodiment, wherein each claim itself All as a separate embodiment of the present invention.
Those skilled in the art will understand that can be carried out adaptively to the module in the equipment in embodiment Change and they are arranged in one or more devices different from this embodiment.It can be the module or list in embodiment Member or component are combined into a module or unit or component, and furthermore they can be divided into multiple submodule or subelement or Sub-component.Other than such feature and/or at least some of process or unit exclude each other, it can use any Combination is to all features disclosed in this specification (including adjoint claim, abstract and attached drawing) and so disclosed All process or units of what method or apparatus are combined.Unless expressly stated otherwise, this specification is (including adjoint power Benefit require, abstract and attached drawing) disclosed in each feature can carry out generation with an alternative feature that provides the same, equivalent, or similar purpose It replaces.
In addition, it will be appreciated by those of skill in the art that although some embodiments described herein include other embodiments In included certain features rather than other feature, but the combination of the feature of different embodiments mean it is of the invention Within the scope of and form different embodiments.For example, in the following claims, embodiment claimed is appointed Meaning one of can in any combination mode come using.
Various component embodiments of the invention can be implemented in hardware, or to run on one or more processors Software module realize, or be implemented in a combination thereof.It will be understood by those of skill in the art that can be used in practice One in the creating device of microprocessor or digital signal processor (DSP) to realize virtual machine according to an embodiment of the present invention The some or all functions of a little or whole components.The present invention is also implemented as executing method as described herein Some or all device or device programs (for example, computer program and computer program product).Such realization Program of the invention can store on a computer-readable medium, or may be in the form of one or more signals.This The signal of sample can be downloaded from an internet website to obtain, and is perhaps provided on the carrier signal or mentions in any other forms For.
Embodiments of the present invention above described embodiment only expresses, the description thereof is more specific and detailed, but can not Therefore limitations on the scope of the patent of the present invention are interpreted as.It should be pointed out that for those of ordinary skill in the art, Without departing from the inventive concept of the premise, various modifications and improvements can be made, these belong to protection model of the invention It encloses.Therefore, the scope of protection of the patent of the invention shall be subject to the appended claims.

Claims (14)

1. a kind of noise-reduction method of pulse array signals characterized by comprising
The granting characteristic of the spatial neighborhood and the pulse array signals on timestamp neighborhood for the treatment of denoising pulse is modeled, and is established The granting characteristic model of pulse array signals extracts pulse array signals spatial-temporal distribution characteristic;
According to the noise profile characteristic in pulse array signals, space time filter is constructed, the feature of extraction is filtered;
According to the granting characteristic model, filtered pulse spatial-temporal distribution characteristic is reduced to pulse array signals.
2. the method according to claim 1, wherein the pulse to be denoised be time-space domain on specific pulse, It is one section of pulse letter in a certain specific pixel, since issuing a certain impulse response, before occurring to next pulse response Number.
3. the method according to claim 1, wherein described special according to the noise profile in pulse array signals Property, space time filter is constructed, the feature of extraction is filtered, comprising:
Pulse array signals of the pulse to be denoised on spatial neighborhood and timestamp neighborhood are set as its first region of search;
Pulse array signals in the first region of search are analyzed, the pulse release information of several pulse trains is obtained, is sent out according to pulse Putting each pulse of information acquisition indicates feature or converts transform domain extraction feature for pulse release information;
The pulse array signals that first region of search includes are converted to pulse strength, on pulse strength domain, calculate specific arteries and veins The similarity of punching and other pulses in the first region of search sets first threshold, retains similarity in the first region of search and is greater than The pulse of first threshold is denoted as the second region of search;
For the second region of search, second threshold is set, judges whether the number of pulses in the second region of search is greater than the second threshold Value;If so, being estimated using the pulse pair pulse to be denoised in the second region of search, and replace it former strong using estimated value Degree;If it is not, then subtracting the second region of search using the first region of search, third region of search is obtained, uses third region of search In pulse pair pulse to be denoised estimated, and replace its former intensity using estimated value.
4. according to the method described in claim 3, it is characterized in that, obtaining each pulse according to pulse release information indicates special Sign, are as follows: pulse strength is determined according to pulse Time Of Release interval, or determines that pulse is strong according to the inverse at pulse Time Of Release interval Degree, or pulse strength is determined according to the logarithm at pulse Time Of Release interval.
5. according to the method described in claim 3, it is characterized in that, its in the calculating certain pulses and the first region of search The similarity of his pulse, comprising: calculate the variance of pulse strength in region of search, the intensity according to variance and pulse to be denoised is big One confidence interval of small determination determines similarity according to distribution situation of the pulse in confidence interval.
6. according to the method described in claim 3, it is characterized in that, the first threshold, second threshold are according to neighborhood pulse Relationship is dynamically set, to adapt to different pulse number evidence.
7. according to the method described in claim 3, it is characterized in that, the first threshold is different from the second threshold.
8. according to the method described in claim 3, it is characterized in that, the pulse pair in second region of search of use waits for Pulse of making an uproar estimated, including but not limited to: weighted mean method, Neighborhood Filtering method.
9. according to the method described in claim 3, it is characterized in that, the denoising pulse for the treatment of carries out being estimated as passing through filtering Device is estimated that the weight or convolution kernel of filter assign bigger value to the pulse closer to pulse to be denoised.
10. according to the method described in claim 4, it is characterized in that, according to the granting characteristic model, filtered pulse Spatial-temporal distribution characteristic is reduced to pulse array signals, comprising: by the intensity value of certain pulses adjusted, is closed according to pulse sequence System, reverts to pulse array signals, and method is described according to each pulse expression this behaviour of feature of pulse release information acquisition Make the inverse operation of step.
11. a kind of noise reduction system of pulse array signals characterized by comprising
Extraction module is modeled, the granting of the pulse array signals on spatial neighborhood and timestamp neighborhood for treating denoising pulse Characteristic is modeled, and the granting characteristic model of pulse array signals is established, and extracts pulse array signals spatial-temporal distribution characteristic;
Filter module, for space time filter being constructed, to the feature of extraction according to the noise profile characteristic in pulse array signals It is filtered;
Recovery module, for according to the granting characteristic model, filtered pulse spatial-temporal distribution characteristic to be reduced to pulse battle array Column signal.
12. a kind of noise reduction system of pulse array signals, which is characterized in that for realizing described in any one of claim 1-10 Method, the system include: sequentially connected neighborhood search module, pulse analysis and conversion module, similarity estimation module and Pulse denoising and recombination module;
The neighborhood search module obtains target search region for searching for the neighborhood pulse array data of target pulse;
The pulse analysis and conversion module obtain each pulse sequence for corresponding to each pixel search region for analyzing pulse data The pulse release information of column, and it is converted into pulse strength;
The similarity estimation module, for calculating the intensity similarity of each pulse and target pulse in region of search;
The pulse denoising and recombination module are reassembled as denoising to target pulse, and the pulse data after denoising Pulse array signals.
13. a kind of electronic equipment, which is characterized in that including memory, processor and be stored on the memory and can be in institute The computer program run on processor is stated, the processor executes described program, to realize as any in claim 1-10 Method described in.
14. a kind of non-transitorycomputer readable storage medium, is stored thereon with computer program, which is characterized in that the program It is executed by processor, to realize such as method of any of claims 1-10.
CN201910577862.3A 2019-06-28 2019-06-28 Method and system for reducing noise of pulse array signal Active CN110430340B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910577862.3A CN110430340B (en) 2019-06-28 2019-06-28 Method and system for reducing noise of pulse array signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910577862.3A CN110430340B (en) 2019-06-28 2019-06-28 Method and system for reducing noise of pulse array signal

Publications (2)

Publication Number Publication Date
CN110430340A true CN110430340A (en) 2019-11-08
CN110430340B CN110430340B (en) 2020-08-25

Family

ID=68408851

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910577862.3A Active CN110430340B (en) 2019-06-28 2019-06-28 Method and system for reducing noise of pulse array signal

Country Status (1)

Country Link
CN (1) CN110430340B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113034634A (en) * 2021-03-04 2021-06-25 北京大学 Adaptive imaging method, system and computer medium based on pulse signal

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107004255A (en) * 2014-09-23 2017-08-01 英特尔公司 The non-local mean image denoising retained with details of mixing is driven using self-similarity
CN107147856A (en) * 2017-03-30 2017-09-08 深圳大学 A kind of pixel cell and its denoising method, dynamic visual sensor, imaging device
US20180253830A1 (en) * 2017-02-24 2018-09-06 Sunnybrook Research Institute Systems and methods for noise reduction in imaging
CN109068024A (en) * 2018-06-29 2018-12-21 北京大学 The method that a kind of pair of space-time signal is filtered
CN109101884A (en) * 2018-07-10 2018-12-28 北京大学 A kind of pulse array prediction technique
CN109902594A (en) * 2019-01-31 2019-06-18 国网江苏省电力有限公司泰州供电分公司 A kind of pulse signal denoising method based on small echo

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107004255A (en) * 2014-09-23 2017-08-01 英特尔公司 The non-local mean image denoising retained with details of mixing is driven using self-similarity
US20180253830A1 (en) * 2017-02-24 2018-09-06 Sunnybrook Research Institute Systems and methods for noise reduction in imaging
CN107147856A (en) * 2017-03-30 2017-09-08 深圳大学 A kind of pixel cell and its denoising method, dynamic visual sensor, imaging device
CN109068024A (en) * 2018-06-29 2018-12-21 北京大学 The method that a kind of pair of space-time signal is filtered
CN109101884A (en) * 2018-07-10 2018-12-28 北京大学 A kind of pulse array prediction technique
CN109902594A (en) * 2019-01-31 2019-06-18 国网江苏省电力有限公司泰州供电分公司 A kind of pulse signal denoising method based on small echo

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113034634A (en) * 2021-03-04 2021-06-25 北京大学 Adaptive imaging method, system and computer medium based on pulse signal
CN113034634B (en) * 2021-03-04 2022-05-06 脉冲视觉(北京)科技有限公司 Adaptive imaging method, system and computer medium based on pulse signal

Also Published As

Publication number Publication date
CN110430340B (en) 2020-08-25

Similar Documents

Publication Publication Date Title
CN109247923B (en) Non-contact type pulse real-time estimation method and device based on video
CN104205169A (en) Method of estimating optical flow on the basis of an asynchronous light sensor
EP3848852B1 (en) Method and apparatus for detecting temporal action of video, electronic device, storage medium, and computer program product
CN112270240B (en) Signal processing method, device, electronic equipment and storage medium
CN112419203B (en) Diffusion weighted image compressed sensing recovery method and device based on countermeasure network
CN116152591B (en) Model training method, infrared small target detection method and device and electronic equipment
KR20090123757A (en) Systems and method for estimating the centers of moving objects in a video sequence
CN113143295A (en) Equipment control method and terminal based on motor imagery electroencephalogram signals
WO2019144455A1 (en) Wavelet analysis-based brain electrical signal denoising method, storage medium and device
CN109068024A (en) The method that a kind of pair of space-time signal is filtered
CN110430340A (en) A kind of noise-reduction method and system of pulse array signals
CN115578262A (en) Polarization image super-resolution reconstruction method based on AFAN model
CN116524402A (en) Multi-time span context modeling action recognition method based on multi-head self-attention
CN115311186B (en) Cross-scale attention confrontation fusion method and terminal for infrared and visible light images
CN112315440A (en) Heart rate detection method, wearable device and readable storage medium
CN116452812A (en) Camouflage object identification and semantic segmentation method
CN116524380A (en) Target detection method based on brain-computer signal fusion
CN115374831B (en) Dynamic and static combination velocity imagery classification method for multi-modal registration and space-time feature attention
CN111479109A (en) Video quality evaluation method, system and terminal based on audio-visual combined attention
CN116725486A (en) Dynamic electrocardiosignal atrial fibrillation detection method and device based on double phases
KR102290243B1 (en) Classification method for Electroencephalography motor imagery
CN115760773A (en) Noise suppression method and system for hyperspectral medical image acquisition
CN112336326B (en) Volume pulse wave signal processing method, blood pressure measuring device, and storage medium
CN114170239A (en) Large-view-field self-adaptive wall filtering method and system based on hierarchical search
CN113827253A (en) Computing device and method for removing noise from electroencephalogram signals

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
TR01 Transfer of patent right

Effective date of registration: 20220127

Address after: 100084 b501b-2, 5th floor, building 8, yard 1, Zhongguancun East Road, Haidian District, Beijing

Patentee after: Pulse vision (Beijing) Technology Co.,Ltd.

Address before: 100871 No. 5, the Summer Palace Road, Beijing, Haidian District

Patentee before: Peking University

TR01 Transfer of patent right