CN105662472A - Method and device for generating Doppler frequency spectrogram - Google Patents
Method and device for generating Doppler frequency spectrogram Download PDFInfo
- Publication number
- CN105662472A CN105662472A CN201610022327.8A CN201610022327A CN105662472A CN 105662472 A CN105662472 A CN 105662472A CN 201610022327 A CN201610022327 A CN 201610022327A CN 105662472 A CN105662472 A CN 105662472A
- Authority
- CN
- China
- Prior art keywords
- signal
- spectrum
- gain
- ambient noise
- doppler
- 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
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
Abstract
The invention provides a method and device for generating a Doppler frequency spectrogram. The method comprises the steps that background noise signals are acquired for at least a part of frequency spectrum signals used for generating the Doppler frequency spectrogram, gains are determined according to the intensity of the background noise signals, and the Doppler frequency spectrogram is generated based on the determined gains. By means of the method and device for generating the Doppler frequency spectrogram, the Doppler frequency spectrogram with a more ideal display effect can be automatically generated, user intervention is not needed, the working intensity of a user is greatly lowered, and the working efficiency is improved.
Description
Technical field
The present invention relates to medical instruments field, in particular it relates to a kind of method and apparatus generating Doppler spectrum.
Background technology
Doppler effect refers to that the wavelength of the ripple of object radiation produces change due to the relative motion of this object and observer. Before the wave source of motion, ripple is compressed, and wavelength becomes shorter, and frequency becomes higher; Time after the wave source of motion, it may occur that contrary phenomenon, wavelength becomes longer, and frequency becomes relatively low. The movement velocity of wave source is more high, and produced Doppler effect is more notable. Therefore, the degree changed according to wave frequency, it is possible to calculate wave source and follow the speed of observed direction motion.
Dopplcr ultrasound blood analysis is to utilize Doppler effect, check, by Noninvasive, a kind of method evaluating different physiologic characters. Transcranial Doppler sonography analyser is the ultrasonic device of a kind of customization, is specifically designed to the ultrasonic examination through skull. Transcranial Doppler sonography analyser uses the wafer in external supersonic probe to launch ultrasound wave (abbreviation transmitted wave) through gap or " window " of skull to cerebrovascular, the existence of blood flow will cause the generation of Doppler effect (Doppler frequency shift), last ultrasound wave is reflected back to probe (abbreviation echo), through same wafer receipt, analyser carry out data process and draw corresponding blood flow information. Owing to adopting degree of depth gating (or range gating) technology, fixed point measuring of blood flow can be carried out, thus there is significantly high range resolution, it is possible to certain character putting blood flow be made and analyzes accurately. Specifically, if a certain detection degree of depth is absent from blood flowing, then not producing Doppler effect, compared with transmitted wave, the mid frequency of echo will not change; And if a certain detection degree of depth exists blood flowing, then can produce Doppler effect, compared with transmitted wave, the mid frequency of echo can offset. If using wave filter that mid frequency is filtered, only retain Doppler frequency deviation composition, then for being absent from the degree of depth of blood flowing, by only remaining ambient noise signal (garbage signal) in spectrum signal, and for there is the degree of depth of blood flowing, spectrum signal will include Doppler signal (useful signal) and ambient noise signal. Ambient noise signal therein is produced by the physical restriction of transcranial Doppler sonography analyser itself.
In order to highlight Doppler signal, it is thus achieved that the better Doppler spectrum of display effect, existing transcranial doppler equipment needs user to manually adjust parameter in Doppler spectrum gatherer process. This requires that user carries out operation bidirectional, adds user's specific works amount. Additionally, manually adjust parameter because the display effect after adjusting cannot be expected, often require that user has certain experience, suitable to be adjusted as early as possible.
Summary of the invention
In order to solve problems of the prior art at least in part, according to an aspect of the invention, it is provided a kind of method generating Doppler spectrum, including:
At least some of for what be used for generating in the spectrum signal of Doppler spectrum, obtain ambient noise signal therein;
Intensity according to described ambient noise signal determines gain; And
Based on determined Gain generating Doppler spectrum.
According to a further aspect in the invention, additionally provide a kind of equipment generating Doppler spectrum, determine device and mapping arrangements including background noise acquisition device, gain.
Background noise acquisition device, at least some of for what be used for generating in the spectrum signal of Doppler spectrum, obtains ambient noise signal therein. Gain determines that device is for determining gain according to the intensity of described ambient noise signal. Mapping arrangements is for based on determined Gain generating Doppler spectrum.
The method and apparatus of above-mentioned generation Doppler spectrum can automatically obtain the better Doppler spectrum of display effect, it is not necessary to user intervention, is substantially reduced user job intensity, improves work efficiency.
Introducing the concept of a series of simplification in summary of the invention, these concepts will further describe in detailed description of the invention part. Present invention part is not meant to the key feature and the essential features that attempt to limit technical scheme required for protection, does not more mean that the protection domain attempting to determine technical scheme required for protection.
Below in conjunction with accompanying drawing, describe advantages and features of the invention in detail.
Accompanying drawing explanation
The drawings below of the present invention is used for understanding the present invention in this as the part of the present invention. Shown in the drawings of embodiments of the present invention and description thereof, it is used for explaining principles of the invention. In the accompanying drawings,
Fig. 1 illustrates the schematic block diagram of transcranial Doppler sonography analyser according to an embodiment of the invention;
Fig. 2 illustrates Doppler spectrum according to an embodiment of the invention;
Fig. 3 illustrates the Doppler spectrum that will obtain after the gain reduction of the Doppler spectrum shown in Fig. 2;
Fig. 4 illustrates the indicative flowchart of the method generating Doppler spectrum according to an embodiment of the invention;
Fig. 5 illustrates the indicative flowchart of the method generating Doppler spectrum in accordance with another embodiment of the present invention;
Fig. 6 illustrates the Doppler spectrum according to further embodiment of the present invention;
Fig. 7 illustrates the indicative flowchart of the method generating Doppler spectrum according to another embodiment of the present invention; And
Fig. 8 illustrates the schematic block diagram of the equipment generating Doppler spectrum according to an embodiment of the invention.
Detailed description of the invention
In the following description, it is provided that substantial amounts of details is so as to understand the present invention up hill and dale. But, those skilled in the art it will be seen that, the presently preferred embodiments of the present invention that only relates to described below, and the present invention can be carried out without one or more such details. Additionally, in order to avoid obscuring with the present invention, be not described for technical characteristics more well known in the art.
Fig. 1 illustrates the schematic block diagram of transcranial Doppler sonography analyser 1000 according to an embodiment of the invention. As it is shown in figure 1, this transcranial Doppler sonography analyser can include main frame 1100 and probe 1200. Main frame 1100 can include discharger 1110, probe socket 1120 and process circuit 1130. Probe 1200 can be multiple.
Discharger 1110 is used for providing transmitting signal. Discharger 1110 can include mission controller 1111 and drive circuit 1112. Mission controller 1111 can provide the pulse train of characteristic frequency as launching signal. Drive circuit 1112 may be used for this transmitting signal is converted to high-voltage signal, to drive probe 1200. The signal of every time launching of mission controller 1111 is equivalent to carry out on a timeline once sampling. This transmitting signal can be expressed as time dependent one-dimensional signal.
Probe 1200 can receive, via probe socket 1120, the high-voltage signal that drive circuit 1112 is changed, and is carried out electro-acoustic conversion, to launch ultrasound wave. Ultrasound wave is sent in tissue and skeleton, and part energy can return to probe. Probe 1200 can also receive the ultrasound wave after measured is reflected, and carries out sound-electric conversion, and the echo comprising velocity of sound information is changed into the signal of telecommunication, to generate reception signal.
Process circuit 1130 and can include amplifier 1131, AD sample circuit 1132 and signal processing module 1133. The reception signal generated through the sound-electric conversion of probe 1200 is typically more faint, therefore alternatively, processes circuit 1130 and can include amplifier 1131, so that the faint signal of telecommunication is converted into the stronger signal of telecommunication. Additionally, the signal of telecommunication of probe 1200 generation is analogue signal, therefore alternatively, processing circuit 1130 and can include AD sample circuit 1132, it is by with sample frequency Fs1To analog signal sampling, it is transformed into digital signal. Sample frequency Fs1It is properly termed as systematic sampling rate. Can also including signal processing module 1133 it addition, process in circuit 1130, it is for processing digital signal, to generate Doppler spectrum. Doppler signal is essentially non-stationary signal, changes over its frequency and also can change therewith. Because there are mapping relations one by one in frequency and measured's blood flow rate, so Doppler signal includes the information of time that detects and the information of blood flow rate detected. This digital signal can be carried out Fourier transformation by signal processing module 1133, so that it is carried out spectrum analysis. Short Time Fourier Transform is a kind of conventional signal processing method for processing signal f (t). Its thought is to select analysis window function g (t) of a Time-Frequency Localization, assuming that analysis window function g (t) is steadily (pseudo-steady) in a short time interval, this analysis window function g (t) mobile, make f (t) g (t) is stationary signal in different finite time width. Thus, calculate signal f (t) at each not power spectrum in the same time. Such that obtain the spectrum expression formula of primary signal f (t). Here, primary signal f (t) is that AD sample circuit 1132 carries out the digital signal that analog digital conversion generates.
Alternatively, main frame 1100 can connect host computer 1300, to be shown according to the Doppler spectrum that spectrum signal generates by host computer 1300. Doppler spectrum is 3-D view, and every display line therein is corresponding to a power spectrum, it is possible to be called spectral line. According to Doppler effect principle, the frequency in spectrum signal is proportional to blood flow rate. So, generally, the abscissa representing time of Doppler spectrum, vertical coordinate represents blood flow rate.
The gain of Doppler spectrum is the increment of the intensity of the primary signal for generating this Doppler spectrum. The gray value one_to_one corresponding of the intensity of signal and the pixel of Doppler spectrum. Each pixel value of the Doppler spectrum generated according to primary signal is all deducted the value corresponding with gain, then obtains desired Doppler spectrum. It is appreciated that if pixel value is less than the value corresponding with gain, then this pixel value assignment is 0 by subtraction, say, that this pixel value will be displayed as ater. In other words, the existence of gain so that the luminance-reduction in Doppler spectrum. For Doppler spectrum, it may be desirable to its gain is arranged so that wherein noise does not just show.
When such as ageing equipment, external interference, the noise of Doppler Analyzer itself can raise, and now much noise appears in Doppler spectrum, affects Consumer's Experience. At this moment need gain reduction, to filter out more picture noise.
Fig. 2 illustrates Doppler spectrum according to an embodiment of the invention. Wherein, shade represents energy relative intensity. In image, the color of pixel is more shallow, then the signal energy corresponding to this pixel is more strong. In image, waveform figure that is relatively light, that have periodic regularity to change is corresponding to Doppler signal. Part relatively dark in image includes figure irregular, random distribution, and it corresponds to ambient noise signal. As in figure 2 it is shown, ambient noise signal generally can be uniformly distributed in signal passband. This Doppler spectrum has illustrated more background noise. Fig. 3 illustrates the Doppler spectrum that will obtain after the gain reduction of the Doppler spectrum shown in Fig. 2. As it is shown on figure 3, after by gain reduction, the Doppler spectrum obtained illustrate only lesser amount of background noise. From the Doppler spectrum shown in Fig. 3 it can be seen that the while that background noise therein being repressed, Doppler signal intensity is also affected by certain impact.
Conversely, for the situation of background noise reduction, such as chip and probe body constitution preferably situation, gain can be properly increased, to observe more Doppler signal.
Therefore, suitable gain is to ensure that the key factor of Doppler spectrum optimal display result. Doppler Analyzer in the market needs user to manually adjust Doppler spectrum, and reasonably to show Doppler signal, this brings operational inconvenience to user. For Consumer's Experience better, according to one aspect of the invention, it is provided that a kind of method generating Doppler spectrum. Fig. 4 illustrates the indicative flowchart of the method 400 generating Doppler spectrum according to an embodiment of the invention. As shown in Figure 4, the method 400 includes step S420, step S440 and step S460.
In the step s 420, at least some of for what be used for generating in the spectrum signal of Doppler spectrum, obtain ambient noise signal therein. As it has been described above, except including Doppler signal in spectrum signal, also include ambient noise signal. There is a unique corresponding Doppler signal in each ambient noise signal. Pixel that this ambient noise signal generates and the pixel that corresponding Doppler signal generates are co-located on the same spectral line of Doppler spectrum. Doppler signal includes the information of the blood circumstance of reflection measured. Ambient noise signal is garbage signal. Doppler signal is different from the intensity of ambient noise signal. It is said that in general, the intensity of ambient noise signal is lower than the intensity of Doppler signal.
Alternatively, for being used for generating the spectrum signal of the spectral line of Doppler spectrum, based on signal intensity, Doppler signal therein and ambient noise signal are separated, to obtain ambient noise signal therein. In other words, namely this separation process separates ambient noise signal and corresponding Doppler signal. This separation process can carry out for the spectrum signal corresponding to a spectral line. This separation process can also be respectively directed to the spectrum signal corresponding to each spectral line in a plurality of spectral line and carry out. In other words, it is respectively directed to the spectrum signal for generating a plurality of spectral line in Doppler spectrum, obtains ambient noise signal therein. Thus, the ambient noise signal obtained corresponds respectively to each bar spectral line in a plurality of spectral line of Doppler spectrum. The bar number of spectral line can be any appropriate numerical value, for instance 300. It is appreciated that these spectral lines may be located at the optional position of Doppler spectrum. Such as, these spectral lines can be Time Continuous, and is positioned at a cardiac cycle of measured; These spectral lines can also is that the time is discontinuous, and it may be located in the different cardiac cycles of measured.
Alternatively, K means clustering method is utilized ambient noise signal and corresponding Doppler signal to be separated. This separation process is substantially the categorizing process of ambient noise signal and Doppler signal. In other words, it is divided into ambient noise signal and Doppler signal by spectrum signal, utilizes K means clustering method to realize this process and not only classify accurately, and simple and quick.
Alternatively, it is possible to first calculate the average of the intensity of spectrum signal. Then, will be greater than the signal of average as Doppler signal, by the as background noise signal of the signal less than average.
Can being obtained by priori, Doppler signal and ambient noise signal are Relatively centralized from respective frequency. Doppler signal is positioned at the scope that frequency is relatively low, and ambient noise signal is positioned at the scope that frequency is higher mostly. Alternatively, above-mentioned priori is utilized can to revise above-mentioned initial gross separation result further. Such as, although be divided for Doppler signal based on some ambient noise signal of strength factor, but be based on the reason that frequency is higher, it is possible to be adjusted to ambient noise signal. Otherwise, although it is divided for ambient noise signal based on some Doppler signal of strength factor, but is based on the reason that frequency is relatively low, it is possible to be adjusted to Doppler signal. Thus, it is possible to obtain more structurally sound final separation result.
It is presented above a kind of specific implementation obtaining ambient noise signal, it will be understood that it is merely illustrative and unrestricted. For example, it is also possible to be directly based upon frequency to obtain ambient noise signal. Such as, will be above the signal of characteristic frequency threshold value as background noise signal.
It is alternatively possible to separated ambient noise signal is put into system cache, process for further analysis.
In step S440, determine the gain of Doppler spectrum according to the intensity of acquired ambient noise signal. The brightness of each pixel of Doppler spectrum is mainly determined by the intensity of corresponding signal. As it was previously stated, gain determines the overall brightness of Doppler spectrum. It is generally desirable in Doppler spectrum, show the least possible ambient noise signal. So, determine that gain can effectively shield the display of the ambient noise signal in Doppler spectrum according to the intensity of ambient noise signal.
In one example, it is possible to using the intermediate value of the intensity of this ambient noise signal as gain. It is appreciated that if the ambient noise signal obtained in the step s 420 corresponds only to a spectral line of Doppler spectrum, then will only calculate the intermediate value of the intensity of this ambient noise signal at this. If the ambient noise signal of each bar spectral line in obtaining a plurality of spectral line corresponding respectively to Doppler spectrum in the step s 420, then the intermediate value of intensity of these ambient noise signals can be calculated at this. In another example, it is possible to this intermediate value is added an empirical value as gain.
In step S460, based on determined Gain generating Doppler spectrum. Based on determined gain, it may be determined that the value that each pixel value of the Doppler spectrum generated according to primary signal should be deducted. After each pixel value of the Doppler spectrum generated according to primary signal is performed subtraction with this value, then generate desired Doppler spectrum. This Doppler spectrum can shield all less than the ambient noise signal corresponding to the pixel of this value, and other ambient noise signals are also effectively weakened, thus highlighting Doppler signal better.
In the method 400 of above-mentioned generation Doppler spectrum, determine gain according to the intensity of ambient noise signal, and according to this Gain Automatic generation Doppler spectrum. The Doppler spectrum more desirably showing Doppler signal can be generated when without user intervention, be substantially reduced user job intensity, improve work efficiency.
Additionally, step S420 obtains ambient noise signal determines gain independent of step S460, which ensure that method 400 occurs without positive feedback reforming phenomena. So that it is guaranteed that generate desired Doppler spectrum.
Alternatively, above-mentioned steps S440 can specifically include step S441 and step S445. In step S441, calculate the mean μ of the intensity of above-mentioned ambient noise signal. In step S445, determine gain according to the mean μ of the intensity of ambient noise signal. As a rule, background noise belongs to thermal noise, and its statistical property meets random normal distribution. Intensity is about 50% more than the ambient noise signal of average. The average of intensity embodies the intensity of ambient noise signal better. The mean μ of the intensity according to ambient noise signal determines that gain is more accurate, thus being more beneficial for method 400 to generate more preferably Doppler spectrum.
Alternatively, before step S445, above-mentioned steps S440 can also include step S443, calculates the variances sigma of the intensity of above-mentioned ambient noise signal. In step S445, it is determined that gain except that according to the average of the intensity of ambient noise signal always according to the variance of the intensity of ambient noise signal. As it was previously stated, background noise belongs to thermal noise, its statistical property meets random normal distribution. The variances sigma of intensity substantially embodies the degree of scatter of ambient noise signal. The variances sigma of intensity is more little, then the intensity of ambient noise signal is more concentrated; The variances sigma of intensity is more big, then the more dispersion of the intensity of ambient noise signal. So, when step S445 determines gain, the average and the variance two aspect factor that consider intensity will determine gain more accurately such that it is able to obtain display effect Doppler spectrum more preferably.
Alternatively, determine that gain farther includes to determine described gain according to equation below according to the average of intensity and variance:
G=μ+k* σ,
Wherein, G represents described gain, and μ represents described average, and σ represents described variance, and k represents the coefficient of described variance, and k is greater than or equal to any real number of 0.
Different users is different for the demand of Doppler spectrum. Some users wish that background noise is completely clean; Some users wish slightly there are some background noises; Some users then wish to become totally visible background noise, see Doppler signal with more vivid. One feasible noise rating method is Parameter Estimation Method. As it was previously stated, the statistical property of ambient noise signal meets normal distribution. Thus, the ratio of the intensity ambient noise signal more than μ+σ is about 15.8%; The ratio of the intensity ambient noise signal more than μ+2* σ is about 2.2%; The ratio of the intensity ambient noise signal more than μ+3* σ is about 0.1%. In other words, by providing different coefficient of variation k, it is provided that the gear that noise level is different, and each gear has relatively-stationary noise level. Thus, meet the demand of different user, improve Consumer's Experience.
Fig. 5 illustrates the flow chart of the method 500 generating Doppler spectrum in accordance with another embodiment of the present invention. As it is shown in figure 5, the method 500 includes step S505, step S510, step S515, step S520, step S541, step S542, step S543, step S544, step S545 and step S560. Wherein step S520 and step S560 is similar with the corresponding step in said method 400 respectively, for sake of simplicity, in this not go into detail.
Before step S520, method 500 also includes step S505: remove the noise in Doppler signal in spectrum signal. Directly all spectrum signals gathered in the time period corresponding to noise signal are abandoned, be not used in subsequent calculations analysis.
In image acquisition process, measured and user can not remain stationary for a long time. Therefore, accidentally there will be little body action (such as measured's cough) and image disruption may be caused. In addition, extraneous various physical conditions are likely to and image are produced interference, such as strong electromagnetic. This interference signal source is not belonging to measured, step S520 can be obtained ambient noise signal and interfere, thus affect the accuracy of gain, finally make the Doppler spectrum automatically generated in method 500 undesirable. Fig. 6 illustrates Doppler spectrum according to an embodiment of the invention, and as shown in Figure 6, this Doppler spectrum includes the part that interference signal is corresponding, i.e. left part elongated portion in figure.
In general, the image that interference signal and the Doppler signal that generates because of blood flow generate has notable difference, for instance in the following areas: energy range and distribution, persistent period, form, periodic regularity etc. Such as, interference signal can make to occur in Doppler spectrum the figure that form is high and sharp. Doppler signal includes the information of time that detects and the information of blood flow rate detected. Noise can be the Doppler signal met the following conditions: the persistent period is less than time threshold, and blood flow rate exceedes the particular percentile of current scale. Such as: the persistent period is less than 100ms, and blood flow rate is significantly high, exceedes the 80% of current scale, it is believed that be exist to disturb in short-term. After correctly identifying and removing noise, it is ensured that the data that subsequent step is analyzed are all useful signal, for obtaining correct ambient noise signal and determining therefrom that gain provides strong guarantee, and then ensure that the Doppler spectrum generated is better.
Before step S520, method 500 also includes step S510 and step S515. In step S510, spectrum signal is carried out cardiac cycle analysis, spectrum signal is divided into peri odic spectrum signal corresponding with a cardiac cycle respectively according to cardiac cycle. In step S515, each peri odic spectrum signal in multiple peri odic spectrum signals selects the part being used for generating the spectrum signal of a plurality of spectral line in Doppler spectrum. Such as, in each peri odic spectrum signal, select the spectrum signal for generating 30 spectral lines in Doppler spectrum. In step S520 and its subsequent step, it is analyzed calculating to the spectrum signal from the n*30 bar spectral line for generating Doppler spectrum selected in n peri odic spectrum signal.
In human body, the speed of blood flowing is typically all and is changed according to cardiac cycle. Heart is at systole, and it outwards penetrates blood, and Ink vessel transfusing blood flow rate is accelerated; Heart is at relaxing period, and Ink vessel transfusing blood flow rate slows down. Cardiac cycle is the normal rhythm of human body, is analyzed calculating for the spectrum signal of multiple complete cardiac cycles and ensure that the information obtained is more valuable.
Step S541 to step S545 in method 500 shown in Fig. 5, corresponding to the step S440 in the method 400 shown in Fig. 4, determines gain according to the intensity of ambient noise signal. These steps in method 500 are described more fully below.
Step S541 and step S542 achieves the average of the intensity calculating acquired ambient noise signal jointly. In step S541, for each peri odic spectrum signal, calculate the Periodic Mean of the intensity of ambient noise signal therein. In step S542, calculate the average of the intensity of acquired ambient noise signal according to described Periodic Mean. For example, it is possible to by rhythmic for institute Periodic Mean averaging, using as described average. Again for example, it is possible to select the intermediate value of rhythmic Periodic Mean as described average.
Step S543 and step S544 achieves the variance of the intensity calculating acquired ambient noise signal jointly. In step S543, for each peri odic spectrum signal, calculate the periodic variance of the intensity of ambient noise signal therein. In step S544, calculate the variance of the intensity of acquired ambient noise signal according to described periodic variance. For example, it is possible to by rhythmic for institute periodic variance averaging, using as described variance. Again for example, it is possible to select the intermediate value of rhythmic periodic variance as described variance.
In step S545, determine gain according to the average of the intensity of acquired ambient noise signal and variance. Step S445 in this process such as said method 400, for sake of simplicity, do not repeat them here.
This method 500 can accumulate the parameter of multiple cardiac cycle, considers to determine gain, obtains more reliable result. The analytical calculation carrying out spectrum signal in units of cardiac cycle ensure that the accuracy of signal computing, it is to avoid error.
This area ordinary person is appreciated that in said method 500, is described with the order of step S505, step S510, step S515, step S520, step S541, step S542, step S543, step S544, step S545 and step S560. But this implementation is merely to illustrate the example of embodiments of the invention, and the present invention is not caused restriction by it. Such as, step S505 can perform after step S510, step S515. Again such as, step S542 can perform after step S543. Step S505 and step S510, step S515, step S541 are absent from dependence each other to these 7 steps of step S545, and it can be individually present, it is also possible to coexists.
Fig. 7 illustrates the method 700 generating Doppler spectrum according to further embodiment of this invention. As it is shown in fig. 7, the method 700 includes step S705, step S710, step S715, step S720, step S741, step S743, step S746, step S747 and step S760. Wherein step S705, step S710, step S715, step S720, step S741, step S743 and step S760 are similar with the corresponding step in said method 500 respectively, for sake of simplicity, in this not go into detail.
Step S741 in method 700 shown in Fig. 7, step S743, step S746 and step S747, corresponding to the step S440 in the method 400 shown in Fig. 4, determine gain according to the intensity of ambient noise signal.
As, described in method 500, in step S741 and step S743, for each peri odic spectrum signal, calculated Periodic Mean and the periodic variance of the intensity of ambient noise signal therein respectively.
In step S746, for each peri odic spectrum signal, determine its cycle gain according to the Periodic Mean of the intensity of ambient noise signal therein and periodic variance. This process is similar to the step S445 in said method 400, for sake of simplicity, do not repeat them here.
In step S747, determine gain according to cycle gain. For example, it is possible to by rhythmic for institute cycle gain averaging, using as described gain. Again for example, it is possible to select the intermediate value of rhythmic cycle gain as described gain.
This method 700 also accumulates the parameter of multiple cardiac cycle, considers to determine gain, obtains more reliable result.
This area ordinary person is appreciated that with method 500 similarly, and in said method 600, the execution sequence of step is merely to illustrate the example of embodiments of the invention, and the present invention is not caused restriction by it.
According to a further aspect of the invention, a kind of equipment for generating Doppler spectrum is additionally provided. Fig. 8 illustrates the schematic block diagram of the equipment 800 generating Doppler spectrum according to an embodiment of the invention. As shown in Figure 8, equipment 800 includes background noise acquisition device 820, device 840 is determined in gain and mapping arrangements 860.
At least some of for in the spectrum signal being used for generating Doppler spectrum of background noise acquisition device 820, obtains ambient noise signal therein. Alternatively, background noise acquisition device 820 is further used for, for the spectrum signal for generating a plurality of spectral line in described Doppler spectrum, obtaining ambient noise signal therein.
Background noise acquisition device 820 farther includes separation module, for utilizing K means clustering method described ambient noise signal and corresponding Doppler signal to be separated.
Gain determines that device 840 is for determining gain according to the intensity of described ambient noise signal.
Mapping arrangements 860 is for based on determined Gain generating Doppler spectrum.
In the equipment 800 of above-mentioned generation Doppler spectrum, determine gain according to the intensity of ambient noise signal, and according to this Gain Automatic generation Doppler spectrum. The Doppler spectrum more desirably showing Doppler signal can be generated when without user intervention, be substantially reduced user job intensity, improve work efficiency.
Alternatively, gain determines that device 840 farther includes mean value computation module and module is determined in gain. Mean value computation module is for calculating the average of the intensity of described ambient noise signal. Gain determines that module is for determining described gain according to described average.
Alternatively, gain determines that device 840 farther includes variance computing module, for calculating the variance of the intensity of described ambient noise signal. Described gain determines that module determines that described gain is always according to described variance.
Alternatively, described gain determines that module determines described gain according to equation below: G=μ+k* σ,
Wherein, G represents described gain, and μ represents described average, and σ represents described variance, and k represents the coefficient of described variance, and k is greater than or equal to any real number of 0.
In one example, equipment 800 farther includes cycle analysis device and selects device. Cycle analysis device is for carrying out cardiac cycle analysis to described spectrum signal, described spectrum signal to be divided into peri odic spectrum signal corresponding with a cardiac cycle respectively according to cardiac cycle. Select device for each peri odic spectrum signal in multiple peri odic spectrum signals selects a part for the spectrum signal of the described a plurality of spectral line for generating in described Doppler spectrum.
Described mean value computation module farther includes Periodic Mean computing unit and average calculation unit. Periodic Mean computing unit is for for each peri odic spectrum signal, calculating the Periodic Mean of the intensity of ambient noise signal therein. Average calculation unit is for calculating described average according to described Periodic Mean.
Described variance computing module farther includes periodic variance computing unit and variance computing unit. Periodic variance computing unit is for for each peri odic spectrum signal, calculating the periodic variance of the intensity of ambient noise signal therein. Variance computing unit is for calculating described variance according to described periodic variance.
In another example, equipment 800 also includes above-mentioned cycle analysis device and above-mentioned selection device. Wherein gain determines device 840 except including above-mentioned Periodic Mean computing unit and periodic variance computing unit, also includes cycle gain and determines unit. Cycle gain determines that unit determines its cycle gain for the periodic variance computed for each peri odic spectrum signal, the Periodic Mean computed according to Periodic Mean computing unit and periodic variance computing unit. Gain determines that device 840 also includes gain determination unit, for determining gain according to above-mentioned cycle gain.
Alternatively, equipment 800 farther includes denoising device, for removing the noise in described Doppler signal. Described Doppler signal includes the information of time that detects and the information of blood flow rate detected, and described noise is the Doppler signal met the following conditions:
Persistent period is less than time threshold; And
Described blood flow rate exceedes the particular percentile of current scale.
Those of ordinary skill in the art are by reading the detailed description above for the method generating Doppler spectrum, it is to be understood that the structure of equipment of above-mentioned generation Doppler spectrum, realization and advantage, therefore repeat no more here.
The present invention is illustrated already by above-described embodiment, but it is to be understood that, above-described embodiment is only intended to citing and descriptive purpose, and is not intended to limit the invention in described scope of embodiments. In addition it will be appreciated by persons skilled in the art that and the invention is not limited in above-described embodiment, more kinds of variants and modifications can also be made according to the teachings of the present invention, within these variants and modifications all fall within present invention scope required for protection. Protection scope of the present invention is defined by the appended claims and equivalent scope thereof.
Claims (12)
1. the method generating Doppler spectrum, including:
At least some of for what be used for generating in the spectrum signal of Doppler spectrum, obtain ambient noise signal therein;
Intensity according to described ambient noise signal determines gain; And
Based on determined Gain generating Doppler spectrum.
2. method according to claim 1, it is characterised in that the described intensity according to described ambient noise signal determines that gain farther includes:
Calculate the average of the intensity of described ambient noise signal;
Described gain is determined according to described average.
3. method according to claim 2, it is characterised in that the described intensity according to described ambient noise signal determines that gain farther includes:
Calculate the variance of the intensity of described ambient noise signal;
Described determine that described gain farther includes according to described average:
Described gain is determined according to described average and described variance.
4. method according to claim 3, it is characterised in that described determine that described gain farther includes to determine described gain according to equation below according to described average and described variance:
G=μ+k* σ,
Wherein, G represents described gain, and μ represents described average, and σ represents described variance, and k represents the coefficient of described variance, and k is greater than or equal to any real number of 0.
5. the method according to any one of Claims 1-4, it is characterised in that described at least some of for what be used for generating in the spectrum signal of Doppler spectrum, obtains ambient noise signal therein and farther includes:
K means clustering method is utilized described ambient noise signal and corresponding Doppler signal to be separated.
6. the method according to claim 3 or 4, it is characterised in that described at least some of for what be used for generating in the spectrum signal of Doppler spectrum, obtains ambient noise signal therein and farther includes:
For being used for generating the spectrum signal of a plurality of spectral line in described Doppler spectrum, obtain ambient noise signal therein.
7. method according to claim 6, it is characterised in that for being used for generating the spectrum signal of a plurality of spectral line in described Doppler spectrum, before obtaining ambient noise signal therein, described method farther includes:
Described spectrum signal is carried out cardiac cycle analysis, described spectrum signal is divided into peri odic spectrum signal corresponding with a cardiac cycle respectively according to cardiac cycle; And
Each peri odic spectrum signal in multiple peri odic spectrum signals selects a part for the spectrum signal of the described a plurality of spectral line for generating in described Doppler spectrum;
The average of the intensity of the described ambient noise signal of described calculating farther includes:
For each peri odic spectrum signal, calculate the Periodic Mean of the intensity of ambient noise signal therein; And
Described average is calculated according to described Periodic Mean;
The variance of the intensity of the described ambient noise signal of described calculating farther includes:
For each peri odic spectrum signal, calculate the periodic variance of the intensity of ambient noise signal therein; And
Described variance is calculated according to described periodic variance.
8. the method according to any one of Claims 1-4, it is characterised in that before the described intensity according to described ambient noise signal determines gain, described method farther includes:
Remove the noise in Doppler signal in described spectrum signal.
9. method according to claim 8, it is characterised in that wherein said Doppler signal includes the information of time that detects and the information of blood flow rate detected, and described noise is the Doppler signal met the following conditions:
Persistent period is less than time threshold; And
Described blood flow rate exceedes the particular percentile of current scale.
10. generate an equipment for Doppler spectrum, including:
Background noise acquisition device, at least some of for what be used for generating in the spectrum signal of Doppler spectrum, obtains ambient noise signal therein;
Device is determined in gain, determines gain for the intensity according to described ambient noise signal; And
Mapping arrangements, for based on determined Gain generating Doppler spectrum.
11. equipment according to claim 10, it is characterised in that described gain determines that device farther includes:
Mean value computation module, for calculating the average of the intensity of described ambient noise signal;
Module is determined in gain, for determining described gain according to described average.
12. equipment according to claim 11, it is characterised in that described gain determines that device farther includes:
Variance computing module, for calculating the variance of the intensity of described ambient noise signal;
Described gain determines that module determines that described gain is always according to described variance.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610022327.8A CN105662472A (en) | 2016-01-13 | 2016-01-13 | Method and device for generating Doppler frequency spectrogram |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610022327.8A CN105662472A (en) | 2016-01-13 | 2016-01-13 | Method and device for generating Doppler frequency spectrogram |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105662472A true CN105662472A (en) | 2016-06-15 |
Family
ID=56300484
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610022327.8A Pending CN105662472A (en) | 2016-01-13 | 2016-01-13 | Method and device for generating Doppler frequency spectrogram |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105662472A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018218441A1 (en) * | 2017-05-27 | 2018-12-06 | 北京悦琦创通科技有限公司 | Spectrogram analysis method, apparatus and device, and computer readable storage medium |
| CN109445882A (en) * | 2017-08-21 | 2019-03-08 | 深圳市鼎阳科技有限公司 | A kind of frequency spectrum three-dimensional display apparatus, method and computer readable storage medium |
| CN111699378A (en) * | 2018-01-05 | 2020-09-22 | 统一半导体公司 | Method and system for optically inspecting a substrate |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5800358A (en) * | 1997-03-31 | 1998-09-01 | Hewlett Packard Company | Undersampled omnidirectional ultrasonic flow detector |
| CN1907230A (en) * | 2005-08-02 | 2007-02-07 | 深圳迈瑞生物医疗电子股份有限公司 | Automatic detection method of frequency spectrum Doppler blood flow velocity |
| CN101138507A (en) * | 2006-09-05 | 2008-03-12 | 深圳迈瑞生物医疗电子股份有限公司 | Doppler bloodstream aural signal processing method and device thereof |
| CN101766497A (en) * | 2008-12-31 | 2010-07-07 | 深圳迈瑞生物医疗电子股份有限公司 | Method for processing signal of sound spectrogram image and system therefor |
| CN101879076A (en) * | 2009-05-08 | 2010-11-10 | 深圳迈瑞生物医疗电子股份有限公司 | Method and device for automatically optimizing Doppler ultrasonic imaging |
| CN102743190A (en) * | 2012-08-01 | 2012-10-24 | 深圳市理邦精密仪器股份有限公司 | Method and system for improving noise ratio of transcranial doppler signals |
-
2016
- 2016-01-13 CN CN201610022327.8A patent/CN105662472A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5800358A (en) * | 1997-03-31 | 1998-09-01 | Hewlett Packard Company | Undersampled omnidirectional ultrasonic flow detector |
| CN1907230A (en) * | 2005-08-02 | 2007-02-07 | 深圳迈瑞生物医疗电子股份有限公司 | Automatic detection method of frequency spectrum Doppler blood flow velocity |
| CN101138507A (en) * | 2006-09-05 | 2008-03-12 | 深圳迈瑞生物医疗电子股份有限公司 | Doppler bloodstream aural signal processing method and device thereof |
| CN101766497A (en) * | 2008-12-31 | 2010-07-07 | 深圳迈瑞生物医疗电子股份有限公司 | Method for processing signal of sound spectrogram image and system therefor |
| CN101879076A (en) * | 2009-05-08 | 2010-11-10 | 深圳迈瑞生物医疗电子股份有限公司 | Method and device for automatically optimizing Doppler ultrasonic imaging |
| CN102743190A (en) * | 2012-08-01 | 2012-10-24 | 深圳市理邦精密仪器股份有限公司 | Method and system for improving noise ratio of transcranial doppler signals |
Non-Patent Citations (1)
| Title |
|---|
| 禹晶等: "《数字图像处理》", 30 September 2015, 北京机械工业出版社 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018218441A1 (en) * | 2017-05-27 | 2018-12-06 | 北京悦琦创通科技有限公司 | Spectrogram analysis method, apparatus and device, and computer readable storage medium |
| CN109445882A (en) * | 2017-08-21 | 2019-03-08 | 深圳市鼎阳科技有限公司 | A kind of frequency spectrum three-dimensional display apparatus, method and computer readable storage medium |
| CN109445882B (en) * | 2017-08-21 | 2020-08-11 | 深圳市鼎阳科技股份有限公司 | Frequency spectrum three-dimensional display device and method and computer readable storage medium |
| CN111699378A (en) * | 2018-01-05 | 2020-09-22 | 统一半导体公司 | Method and system for optically inspecting a substrate |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109875529B (en) | Vital sign detection method and system based on ultra-wideband radar | |
| JP4831465B2 (en) | Optimization of ultrasonic collection based on ultrasonic detection index | |
| US8979760B2 (en) | Method and apparatus for automatic optimization of doppler imaging parameters | |
| DE102013001230B4 (en) | Axis-related characterization of shear waves with ultrasound | |
| US10278671B2 (en) | Shear wave detection in medical ultrasound imaging | |
| US20110118606A1 (en) | Adaptively performing clutter filtering in an ultrasound system | |
| JP5333427B2 (en) | HEART RATE DETECTOR, HEART RATE DETECTING METHOD, AND PROGRAM | |
| CN110141218B (en) | Electrocardiosignal classification method, electrocardiosignal classification device, program product and storage medium | |
| KR101117035B1 (en) | Ultrasound system and method of performing surface-rendering on volume data | |
| CN107468250B (en) | Biological tissue terahertz imaging method, system and equipment based on multi-scale entropy | |
| Wang et al. | Noncontact heart rate measurement based on an improved convolutional sparse coding method using IR-UWB radar | |
| CN105662472A (en) | Method and device for generating Doppler frequency spectrogram | |
| US20220202390A1 (en) | Fetal ultrasound processing unit for separating heart rate signals | |
| CN105476661A (en) | Method and equipment for regulating Doppler frequency spectrogram | |
| CN116172539A (en) | Vital sign detection method, system, equipment and medium based on machine learning | |
| US20230086332A1 (en) | High-Sensitivity and Real-Time Ultrasound Blood Flow Imaging Based on Adaptive and Localized Spatiotemporal Clutter Filtering | |
| CN113591769A (en) | Non-contact heart rate detection method based on photoplethysmography | |
| CN112168210B (en) | Medical image processing terminal, ultrasonic diagnostic apparatus, and fetal image processing method | |
| CN105476664B (en) | The method and apparatus for generating Doppler spectrum | |
| CN105476663A (en) | Method and equipment for regulating Doppler frequency spectrogram | |
| KR102447020B1 (en) | Apparatus and method for displaying ultrasound image | |
| CN116115202A (en) | Method for detecting human breath and heartbeat frequency by ultra-wideband radar based on improved EMD | |
| US20220296209A1 (en) | Fetal ultrasound processing unit | |
| EP2853918B1 (en) | Shear wave detection in medical ultrasound imaging | |
| KR102223058B1 (en) | Ultrasound system and clutter filtering method thereof |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160615 |