CN104783836B - Ultrasonic system Doppler signal gap interpolation method and system - Google Patents
Ultrasonic system Doppler signal gap interpolation method and system Download PDFInfo
- Publication number
- CN104783836B CN104783836B CN201510140766.4A CN201510140766A CN104783836B CN 104783836 B CN104783836 B CN 104783836B CN 201510140766 A CN201510140766 A CN 201510140766A CN 104783836 B CN104783836 B CN 104783836B
- Authority
- CN
- China
- Prior art keywords
- gap
- frequency
- complementary series
- interpolation
- 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.)
- Active
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The present invention provides a kind of ultrasonic system Doppler signal gap interpolation method and system；Methods described includes：Obtain two sections of adjacent Doppler's scanning signals；Gap between two sections of Doppler's scanning signals carries out interpolation calculation, obtains the interpolation sequence in the insertion gap；Frequency domain correction is carried out to the interpolation sequence, the time-domain signal in the insertion gap is obtained.Ultrasonic system Doppler signal gap interpolation method of the invention and device, the interpolation sequence in insertion gap is obtained using real-time adaptive Doppler interpolation scheme；Further, frequency domain correction is carried out to the interpolation sequence, the time-domain signal in the insertion gap is obtained；Make the high frequency pseudo of the spectrogram of acquisition less, and it is more continuous in the fast-changing place of frequency spectrum, improve the interpolation quality of Doppler frequency spectrum.
Description
Technical field
The invention belongs to ultrasonic diagnostic imaging technical field, it is related to a kind of ultrasonic system Doppler signal gap interpolation method
And system.
Background technology
Diasonograph（B ultrasonic machine）With exclusive characteristic, such as real-time；Wherein, the dynamic commonly used in diasonograph
Blood flow and doppler imaging, as one of indispensable means of auxiliary diagnosis of modern medicine；Doppler imaging can be conveniently fast
The specific flow velocity of blood flow is measured promptly, the criterion as some illnesss in clinical diagnosis.
Doppler（Abbreviation PW）In checking process, generally using both of which：I.e. single synchronous mode and many synchronous modes；Institute
State that single synchronous mode is simpler, whole scanning process is only according to certain speed（Abbreviation PRF in Doppler's scanning）How general update
Data are strangled, and other patterns, such as B-mode, color mode are in and stop frozen state；When needing to be swept using other patterns
When retouching, doctor's active activation other patterns are not described in detail herein.However, there are some defects in single synchronous mode, for example：
, it is necessary to keep probe in stationary state during Doppler's scanning.
Many synchronization scanning patterns solve the problems, such as appeal, during whole scanning, B-mode, color mode and Doppler
The successive scanning of pattern, doctor can simultaneously observe the image of multiple patterns, however, the Doppler's scanning under multi-mode there is also one
A little problems, wherein, topmost problem is：After multiple successive scannings of pattern, cause have larger space between doppler data,
Time which results in doppler data is upper discontinuous, Doppler frequency spectrum Quality Down.
Further, in order to solve the Doppler's scanning produced problem under multi-mode.In the prior art, employ various
Scheme solves the above problems.
As shown in figure 1, Chinese patent application, publication number：CN101336830, discloses a kind of for ultrasonic diagnostic imaging
Orthogonal Doppler signal gap filling method and device.
The method fills gap by the signal positive sequence before and after gap, and takes before and after gap each one section as intermediate zone,
Filling signal and original signal weighted superposition in intermediate zone（Weight coefficient sum is 1）.Filled out by auto-correlation function method
The frequency parameter of front and rear end signal, signal averaging frequency, signal power and signal bandwidth are filled, when certain threshold value is met, to inserting
Value signal is multiplied by a complex signal for frequency, with compensate signal before and after gap it is larger when interpolation caused by frequency spectrum it is discontinuous.However,
The method there is also drawback, and it can not be solved the problems, such as in sideband signal phase hit；Meanwhile, the method is simply to inserting
Complement signal is multiplied by the multiple sinusoidal signal of front and rear difference on the frequency, does not account for frequency change transition in time, and gap when
Between width influence, cause frequency to have a mutation, form high frequency pseudomorphism, causing has obvious vertical line on the spectrogram shown in Fig. 1.
In addition, the background section of the Chinese patent application also refer to U.S. Patent number：US5016641、
US5476097, US5891036, US5642732, US4559952, and U.S. Patent Application No. US2007/0049823 etc. are existing
There is technology.However, in the prior art of above-mentioned reference, in carrying out Interpolation Process between two Doppler's scanning signals, can go out
Existing SPA sudden phase anomalies problem, frequency discontinuity problem, at least one in amplitude hit problem and real time problems.
As shown in Fig. 2 U.S. Patent number：US5016641, it is entitled：《Spectral interpolation of
ultrasound doppler signal》Technical scheme, its utilize the signal before gap power spectrum and generation one with
Machine phase, time-domain signal is obtained by IFFT, and the time-domain signal obtains filling out for no saltus step after windowing process, the addition that partly overlaps
Fill signal.The generation method of random phase is that, plus Gaussian distributed on the basis of original signal phase, variance is（PI/10,
PI/2）Random phase.By on the spectrogram that the technical scheme is obtained, because its interpolation method is in Doppler signal and gap
Adjacent place, it is understood that there may be mutation, so have obvious vertical line, while sound has " Dth Dth " sound.
As shown in figure 3, U.S. Patent number：US5476097, it is entitled：《Simultaneous Ultrasonic
Imaging and Doppler Display System》Technical scheme, it uses frequency spectrum harmony cent to open filling, wherein sound
Message number is directly spliced with the signal inverted sequence at gap two ends, it is ensured that gap two ends it is continuous, but there is saltus step in the middle of splicing, pass through
Carrying out LPF after digital-to-analogue conversion again can eliminate the discontinuous noise for causing.The spectrogram obtained by the technical scheme
On, due to merely with spectrogram interpolation, therefore the spectrogram for obtaining is very fuzzy.
As shown in figure 4, U.S. Patent Application No.：US2007/0049823, it is entitled：《Method for Processing
Doppler Signal Gaps》Technical scheme, it uses and data before and after gap is carried out with high-pass filtering, inverted sequence conjugation weighting
Addition obtains signal in gap, and this method can ensure the continuity in stitching portion signal, not be mutated acoustically, but such as
Shown in Fig. 4, there are some thinner lines on the frequency spectrum for obtaining by the method, therefore do not ensure that equally and be not mutated.
As shown in figure 5, U.S. Patent number：US5891036, it is entitled：《Ultrasonic wave Doppler
diagnosing apparatus》Technical scheme, it uses the data inverted sequence before gap to read, take conjugation or real part with it is empty
Portion is stored in interstitial site after exchanging.Produce phase discontinuous due to taking conjugation post gap starting point, solve the phase deviation amount, it is right
Compensation data phase difference in gap.Second method carries out same operation to the data before and after gap, after then processing
Two segment data weighted superpositions, it is ensured that all continuously smooth transition of gap two ends.As shown in figure 5, the frequency spectrum obtained using the method
There is no obvious vertical line, simply relatively former frequency spectrum somewhat " void ", the place for particularly being changed greatly in frequency, in figure
Frequency spectrum has the phenomenon of " fracture " at frequency rapid increase.
Further, U.S. Patent number：US5642732, it is entitled：《Apparatus and Method for
Estimating Missing Doppler Signals and Spectra》Technical scheme, it is based on autoregression（AR）Mould
The Doppler signal gap filling method that type is estimated, in AR models signal extend be Power estimation natural result, and this spectrum
Method of estimation does not need any adding window, and preferable effect can be obtained in theory；But its exponent number and model parameter estimation method are all
Can have an impact to performance.
Further, U.S. Patent number：US4559952, it is entitled：《Method of Ultrasonically
Measuring Blood Flow Velocity 》Technical scheme, its spectrum distribution for trying to achieve gap front signal first, by width
The coefficient of wave filter is controlled to produce the spectral characteristic for needing, directly to produce signal to be filled by a wave filter with noise
The back-page Doppler signal of reading is connect, between the two window function windowed overlapping.The method and U.S. Patent number
US5016641, U.S. Patent Application No. US2007/0049823 are compared, and even more without clear superiority part, it is according to frequency spectrum in addition
Generation filter coefficient, realizes more complicated in real time.
The content of the invention
Mesh of the invention is to provide a kind of ultrasonic system Doppler signal gap interpolation method and device.
One of for achieving the above object, between an embodiment of the present invention provides a kind of ultrasonic system Doppler signal
Gap interpolation method, the described method comprises the following steps：
M1, adjacent two sections of Doppler's scanning signals are obtained, Doppler behind Doppler signal and gap respectively before gap
Signal；
M2, the gap between two sections of Doppler's scanning signals carry out interpolation calculation, obtain the interpolation in the insertion gap
Sequence；
M3, frequency domain correction is carried out to the interpolation sequence, obtain the time-domain signal in the insertion gap.
It is as a further improvement on the present invention, described that " row interpolation meter is entered in the gap between two sections of Doppler's scanning signals
Calculate, obtain the interpolation sequence in the insertion gap " specifically include：
S1, two sections of Doppler's scanning signals are carried out with value, obtain preceding complementary series to be inserted and rear complementary series to be checked；
S2, the preceding complementary series to be inserted and rear complementary series to be checked are carried out respectively phase place is corrected, frequency correction, bandwidth rectify
Just；
S3, step S2 is corrected after preceding complementary series to be inserted be added after rear complementary series to be checked difference adding window, obtain interpolation
Sequence.
As a further improvement on the present invention, the frequency correction is specifically included：To Doppler signal before gap and gap
General Le signal carries out FFT afterwards, calculates its FFT result；
The corresponding centre frequency of forward and backward complementary series to be inserted and bandwidth are obtained according to the FFT result；
Complementary series frequency compensation scope to be inserted before being determined according to the centre frequency for obtaining, and forward and backward treat interpolation sequence
Row compensation frequency sequence；
Complementary series forward and backward to be checked after frequency correction is obtained according to the forward and backward complementary series compensation frequency sequence to be inserted.
As a further improvement on the present invention, the bandwidth correction is specifically included：
By setting the filtering that wave filter is carried out to forward and backward complementary series to be inserted respectively；
The wave filter is IIR types or FIR types.
As a further improvement on the present invention, it is described " frequency domain correction to be carried out to the interpolation sequence, between obtaining described in insertion
The time-domain signal of gap " is specifically included：
P1, to the interpolation its superimposition random phase noise；
P2, the interpolation sequence after step P1 is done frequency domain correction；
P3, adding window is done to the interpolation sequence after step P2 fill up, obtain the time-domain signal in the insertion gap.
As a further improvement on the present invention, the frequency domain correction is specifically included：
According to forward and backward complementary series to be inserted, and the interpolation sequence in the gap is inserted, obtain the gap in the gap
Frequency spectrum；
Judge the mould of the gap frequency spectrum whether between the mould of forward and backward complementary series frequency spectrum to be inserted；
If so, continuing step P3；
If it is not, after being adjusted to the amplitude of the gap frequency spectrum, then carry out step P3.
One of for achieving the above object, between an embodiment of the present invention provides a kind of ultrasonic system Doppler signal
Gap interpolating apparatus, described device includes：Signal acquisition module, for obtaining two sections of adjacent Doppler's scanning signals, respectively
Doppler signal behind Doppler signal and gap before gap；
Interpolating module, interpolation calculation is carried out for the gap between two sections of Doppler's scanning signals, obtains insertion described
The interpolation sequence in gap；
Frequency-domain correction module, for carrying out frequency domain correction to the interpolation sequence, obtains the time domain letter in the insertion gap
Number.
As a further improvement on the present invention, the interpolating module specifically for：
Two sections of Doppler's scanning signals are carried out with value, preceding complementary series to be inserted and rear complementary series to be checked is obtained；
Carry out phase place is corrected, frequency correction, bandwidth correction respectively to the preceding complementary series to be inserted and rear complementary series to be checked；
It is added after preceding complementary series to be inserted after above-mentioned correction is distinguished into adding window with rear complementary series to be checked, obtains interpolation sequence.
As a further improvement on the present invention, the interpolating module enter line frequency correct when specifically for：
FFT is carried out to general Le signal behind Doppler signal before gap and gap, its FFT result is calculated；
The corresponding centre frequency of forward and backward complementary series to be inserted and bandwidth are obtained according to the FFT result；
Complementary series frequency compensation scope to be inserted before being determined according to the centre frequency for obtaining, and forward and backward treat interpolation sequence
Row compensation frequency sequence；
Complementary series forward and backward to be checked after frequency correction is obtained according to the forward and backward complementary series compensation frequency sequence to be inserted.
As a further improvement on the present invention, the interpolating module carry out bandwidth correct when specifically for：
By setting the filtering that wave filter is carried out to forward and backward complementary series to be inserted respectively；
The wave filter is IIR types or FIR types.
As a further improvement on the present invention, the frequency-domain correction module specifically for：
To the interpolation its superimposition random phase noise；
Frequency domain correction is done to the interpolation sequence after superposition random phase noise；
Do adding window and fill up to doing the interpolation sequence after frequency domain correction, obtain the time-domain signal in the insertion gap.
As a further improvement on the present invention, the frequency-domain correction module be additionally operable to when frequency domain is corrected：
According to forward and backward complementary series to be inserted, and the interpolation sequence in the gap is inserted, obtain the gap in the gap
Frequency spectrum；
Judge the mould of the gap frequency spectrum whether between the mould of forward and backward complementary series frequency spectrum to be inserted；
If so, doing adding window and filling up to doing the interpolation sequence after frequency domain correction, the time-domain signal in the insertion gap is obtained；
If it is not, after being adjusted to the amplitude of the gap frequency spectrum, then do adding window to doing the interpolation sequence after frequency domain correction
Fill up, obtain the time-domain signal in the insertion gap.
Compared with prior art, the beneficial effects of the invention are as follows：Ultrasonic system Doppler signal gap interpolation of the invention
Method and device, the interpolation sequence in insertion gap is obtained using real-time adaptive Doppler interpolation scheme；Further, to described
Interpolation sequence carries out frequency domain correction, obtains the time-domain signal in the insertion gap；Make the high frequency pseudo of the spectrogram of acquisition compared with
It is few and more continuous in the fast-changing place of frequency spectrum, improve the interpolation quality of Doppler frequency spectrum.
Brief description of the drawings
Fig. 1 is the spectrogram of the ultrasonic system Doppler signal of first embodiment in background technology；
Fig. 2 is the spectrogram of the ultrasonic system Doppler signal of second embodiment in background technology；
Fig. 3 is the spectrogram of the ultrasonic system Doppler signal of the 3rd implementation method in background technology；
Fig. 4 is the spectrogram of the ultrasonic system Doppler signal of the 4th implementation method in background technology；
Fig. 5 is the spectrogram of the ultrasonic system Doppler signal of the 5th implementation method in background technology；
Fig. 6 is ultrasonic imaging flowage structure schematic diagram in ultrasonic system in an embodiment of the present invention；
Fig. 7 is the flow chart of ultrasonic system Doppler signal gap interpolation method in an embodiment of the present invention；
Fig. 8 is a kind of structural representation of ultrasonic system Doppler signal gap interplotation system in an embodiment of the present invention
Figure；
Fig. 9 is Doppler's scanning structural representation of double synchronous modes in an embodiment of the present invention；
Figure 10 is the structural representation for carrying out windowing process in an embodiment of the present invention to Doppler signal；
Figure 11 A are the height of the spectrogram obtained respectively using the implementation method of background technology the 5th and using method shown in Fig. 7
Frequency pseudo noise contrast schematic diagram；
Figure 11 B are the frequency of the spectrogram obtained respectively using the implementation method of background technology the 5th and using method shown in Fig. 7
The contrast schematic diagram of spectrum speed change.
Specific embodiment
Below with reference to each implementation method shown in the drawings, the present invention will be described in detail.But these implementation methods are not
The limitation present invention, structure, method or change functionally that one of ordinary skill in the art is made according to these implementation methods
Change and be all contained in protection scope of the present invention.
As shown in fig. 6, ultrasonic imaging flowage structure schematic diagram in ultrasonic system.
Ultrasonic signal in tissue is changed into electric analoging signal via the different primitives of probe transducers, by preceding amplification module
Amplify, then data signal is converted to by A/D D/A converter modules；The data signal of each different primitive is by Beam synthesis mould
Block, synthesizes radiofrequency signal；Radiofrequency signal is by after quadrature demodulation, quadrature demodulation result i/q signal feeding being processed accordingly
Module.
Under many synchronous modes, ultrasonic system utilizes time-multiplexed method, periodically different in different time scanning
Pattern.Hereinafter, will be with double synchronous modes（B +PW）As a example by do specific introduction.
With reference to shown in Fig. 7, with reference to shown in Fig. 9, the ultrasonic system Doppler signal gap that an embodiment of the present invention is provided
Interpolation method is comprised the following steps：
Two sections of adjacent Doppler's scanning signals of M1, acquisition；
M2, the gap between two sections of Doppler's scanning signals carry out interpolation calculation, obtain the interpolation in the insertion gap
Sequence.
Although it is understood that the application is with double synchronous modes（B +PW）As a example by do specific introduction, but the application sweeps
The fruit that comes to an end is obtained based on two sections of Doppler's scanning signals, therefore, can also be by B between two sections of adjacent Doppler's scanning signals
Patten transformation is that color mode or B-mode variegate the forms such as pattern, and the application is not between two sections of Doppler's scanning signals
Scanning pattern be specifically limited.
M3, frequency domain correction is carried out to the interpolation sequence, obtain the time-domain signal in the insertion gap.
In present embodiment, describe for convenience, before adjacent two ends Doppler's scanning signal is respectively defined as into gap
Doppler signal Pre_Gap, and Doppler signal Post_Gap behind gap.
It should be noted that describing for convenience, in the following embodiments, the explanation to each alphabetical noun is general
Occur once, but the explanation, it is applicable with later description content, it is not described in detail herein.
Further, the interpolation sequence that the gap is obtained in step M2 is specifically included：
S1, two sections of Doppler's scanning signals are carried out with value, obtain preceding complementary series to be inserted and rear complementary series to be checked.
It should be noted that in the following embodiments, for convenience, current complementary series to be inserted treats interpolation sequence with after
Row occur simultaneously when, forward and backward complementary series to be inserted is referred to as.
In present embodiment, the step S1 is specifically included：
Set the repeat its transmission frequency of Doppler's scanning；
Between obtaining between every section of time cycle of Doppler's scanning, and two sections of time cycles of Doppler's scanning of acquisition
The gap scanning time cycle；During scanning, every section of time cycle all same of Doppler's scanning；
According to repeat its transmission frequency, the time cycle of above-mentioned Doppler's scanning, and gap scanning time cycle, obtain every
Interpolation line number in section Doppler's number of scanning line, and the gap；In present embodiment：Calculated using following algorithm.
Assuming that the time cycle of every section of Doppler's scanning is：T_pw, the gap scanning time cycle is T_gap, how general
Strangle scanning repeat its transmission frequency be：prf；
So, every section of Doppler's number of scanning line：N_pw = T_pw*prf；
Interpolation line number in the gap：N_gap = T_gap*prf.
Further, according to obtaining every section of Doppler's number of scanning line, and the interpolation line number in the gap, take it is described between
The conjugation of the Doppler's scanning signal after Doppler's scanning signal and the gap before gap, obtains preceding complementary series to be inserted and is treated with after
Look into complementary series.
In present embodiment：Calculated using following algorithm.
Pre_Sig_Conj(k) =conj( Pre_Gap(N_gap-k));k<= N_inter;
Post _Sig_ Conj (k) =conj(Post _Gap(N_inter - k)); k<= N_inter;
Wherein, Pre_Sig_Conj (k) represents preamble interpolation sequence；Post _ Sig_ Conj (k) represents follow-up interpolation
Sequence；N_inter=min (N_pw, N_gap), i.e. N_inter are both N_gap and N_pw smaller value, conj（.）It is right to represent
Value inside bracket takes conjugation.
Further, the interpolation sequence for obtaining the gap also includes：
S2, the preceding complementary series to be inserted and rear complementary series to be checked are carried out respectively phase place is corrected, frequency correction, bandwidth rectify
Just.
It is described phase place is corrected is carried out to the preceding complementary series to be inserted and rear complementary series to be checked specifically to wrap in present embodiment
Include：
Calculating preamble interpolation sequence and follow-up interpolation sequence needs RA.
In present embodiment：Calculated using following algorithm.
Angle_pre = angle(Pre_Gap(N_gap)) - angle (Pre_Sig_ Conj (1));
Angle_post = angle (Post _Gap(1) ) - angle (Post_Sig_ Conj (N_
inter));
Wherein, Angle_pre represents the angle of preceding complementary series correction to be inserted；Complementary series correction to be inserted after Angle_post
Angle, operator angle () represents the phase angle for taking bracket the inside value..
Further, the unit vector of RA is multiplied by preceding sequence to be inserted and rear sequence to be inserted respectively, enters line phase
Correction, obtains new complementary series forward and backward to be inserted.
In present embodiment：Calculated using following algorithm.
Pre_Sig_ForGap(k) = Pre_Sig_ Conj (k)*exp(Angle_pre *i);
Post _Sig_ForGap(k) = Post _Sig_ Conj (k)*exp(Angle_ post *i)
Wherein, Pre_Sig_ForGap (k) is the complementary series preceding to be inserted after phase place is corrected；Post _Sig_ForGap(k)
Complementary series rear to be inserted after for phase place is corrected；I represents imaginary unit.
In present embodiment, after phase place is corrected is carried out to forward and backward complementary series to be inserted, further, to described preceding to be inserted
Complementary series and rear complementary series to be checked are entered line frequency correction and are specifically included：
FFT is carried out to general Le signal Post_Gap behind Doppler signal Pre_Gap before gap and gap, it is calculated
FFT result, the complementary series frequency spectrum forward and backward to be inserted of acquisition is respectively：Pre_fft and post_fft；And according to pre_fft and
Post_fft obtains the corresponding centre frequency of forward and backward complementary series to be inserted：Pre_f0, post_f0；With a width of：Pre_B and
post_B；
The FFT is fast Fourier transform（The abbreviation of Fast Fourier Transformation）, it is that discrete Fourier becomes
The fast algorithm for changing, it is the characteristic such as odd, even, empty, real according to discrete fourier transform, and the algorithm to DFT enters
Row improves what is obtained.
It is keep frequency smooth transition in present embodiment, forward and backward sequence to be filled up needs in chronological order successively progressively
Change its centre frequency, while in view of the time slot between Doppler's scanning signal and gap, it is determined that preceding complementary series to be inserted
Frequency compensation scope is [0, (post_f0-pre_f0）* T_gap/T_pw], rear complementary series frequency compensation scope to be inserted is：[-
(post_f0-pre_f0）*T_gap/T_pw, 0 ]；
Further, forward and backward complementary series compensation frequency sequence to be inserted is determined according to the centre frequency for obtaining；And root
The complementary series forward and backward to be checked after frequency is corrected is obtained according to the forward and backward complementary series compensation frequency sequence to be inserted.
In present embodiment：Calculated using following algorithm and entered line frequency adjustment to forward and backward complementary series to be inserted respectively：
Pre_fre_Diff(k) = k*(post_f0-pre_f0）*T_gap/T_pw/ N_gap；
Pre_Sig_ForGap(k) = Pre_Sig_ ForGap (k)*exp(k*Pre_fre_Diff(k) *i);
Post_fre_Diff(k) = (N_gap-k)*( pre_f0 - post_f0）*T_gap/T_pw/ N_gap;
Post _Sig_ForGap(k) = Post _Sig_ ForGap (k)*exp((N_gap-k)*Post_fre_
Diff(k) *i)；
Wherein, 0< k<= N_inter；Pre_fre_Diff is preceding complementary series compensation frequency sequence to be inserted, Post_fre_
Diff is rear complementary series compensation frequency sequence to be inserted；Pre_Sig_ForGap (k) is the complementary series preceding to be inserted after frequency correction；
Post _ Sig_ForGap (k) is the complementary series rear to be inserted after frequency correction；
In present embodiment, after line frequency correction is entered to forward and backward complementary series to be inserted, further, to described preceding to be inserted
Complementary series and rear complementary series to be checked carry out bandwidth correction and specifically include：
In order to ensure that complementary series to be inserted keeps the bandwidth uniformity with neighboring Doppler data respectively, treated to forward and backward respectively
Interpolation sequence is filtered accordingly.
In present embodiment, forward and backward complementary series to be inserted can accordingly be filtered respectively by setting wave filter.
The wave filter can be IIR or FIR types.
In present embodiment：Calculated using following algorithm.
It is to preceding sequence filter design requirement to be inserted：
Centre frequency is pre_filter_f0=cof1*pre_f0;
Bandwidth requirement is pre_filter_B=cof2*pre_B;
It is to rear sequence filter design requirement to be inserted：
Centre frequency is post_filter_f0=cof3*post_f0;
Bandwidth requirement is post_filter_B=cof4*post_B;
Wherein, cof1, cof2, cof3, cof4 are coefficient, and its value is defaulted as 1.0；
In present embodiment, after being filtered accordingly to forward and backward complementary series to be inserted, Pre_Sig_ForGap_ is obtained
Filtered and post _ Sig_ForGap_filtered, Pre_Sig_ForGap_filtered be after bandwidth correction before treat
Interpolation sequence；Post _ Sig_ForGap_filtered is the complementary series rear to be inserted after bandwidth correction.
It should be noted that in present embodiment, sequentially might not be corrected according to phase place is corrected, frequency correction, bandwidth
Order forward and backward complementary series to be inserted is corrected, as long as ensure bandwidth correction frequency correction after complete, you can with complete
Method shown in the present invention, is not described in detail herein.
Further, the interpolation sequence for obtaining the gap also includes：
S3, step S2 is corrected after preceding complementary series to be inserted be added after rear complementary series to be checked difference adding window, obtain interpolation
Sequence.
To ensure progressively transition, to Pre_Sig_ForGap_filtered and post _ Sig_ForGap_filtered points
It is added after other adding window and obtains interpolation sequence：Sig_For_Gap.
In present embodiment, if N_inter<N_gap, needs the benefit 0 after Pre_Sig_ForGap_filtered to make it
Length is equal to N_gap, same its length is equal to N_Gap, it is necessary to mend 0 before post _ Sig_ForGap_filtered.
In present embodiment：Calculated using following algorithm.
Sig_For_Gap(k) =
Pre_Sig_ForGap_filtered(k)*win(k) + (1-win(k))* post _Sig_ForGap_
filtered(k);
Wherein, Win is the coefficient of institute's adding window, for example：Can be oblique angle window win (k)=1- k/N_Gap；Certainly,
The shape of other gradually transition is can be designed as, is not described in detail herein.
Further, frequency domain correction is carried out in step M3 to the interpolation sequence, the time domain letter in the insertion gap is obtained
Number specifically include：
P1, to the interpolation its superimposition random phase noise；
In present embodiment, because the more front and rear Doppler signal bandwidth of fir filter interpolator bandwidth may be narrower, to keep frequency
The uniformity of the background white-noise level of spectrum is, it is necessary to the corresponding random phase noise of Signal averaging.
In present embodiment：Calculated using following algorithm.
Sig_For_Gap(k) = Sig_For_Gap(k)*exp(Phase_random(k)*i);
Wherein, Phase_random is random phase, and it can be designed as being evenly distributed, and is also designed to Gauss and waits other
Distribution form, is not described in detail herein.
Further, the step M3 also includes：
P2, the interpolation sequence after step P1 is done frequency domain correction；
In present embodiment：Calculated using following algorithm.
With data grid technology in gap, length is taken for N_pw data are used to calculate FFT, wherein the N_pw data for obtaining are used
Sig_gap_ForFFT is represented.
In this way, record a demerit crossing coincidence length for N_Overlap=(N_pw-N_gap)/2;
Sig_gap_ForFFT ( k) = Pre_Gap(N_pw- N_ Overlap +k);
k <= N_ Overlap;
Sig_gap_ForFFT ( k) = Sig_For_Gap(k - N_ Overlap);
N_ Overlap< k <= N_ Overlap+N_gap;
Sig_gap_ForFFT ( k) = Post _Gap (k- N_ Overlap- N_gap);
N_ Overlap+N_gap < k <= N_pw; Overlap+N_gap
Further, FFT is done to Sig_gap_ForFFT and obtains frequency spectrum gap_FFT, according to pre_FFT and post_FFT
Gap_FFT is adjusted, gap_adjusted_FFT is obtained.
Further, whether the mould of gap_FFT (k) is judged between pre_FFT (k) and the mould of post_FFT (k),
If so, continue step P3, if it is not, then the amplitude to gap_FFT (k) is adjusted；Wait to be adjusted between pre_FFT (k)
After and the mould of post_FFT (k) between, then carry out step P3.
Further, in present embodiment, adjustment ratio can use nonlinear way, also can be using simple linear side
Formula.
In present embodiment, further, set following program segments judge gap_FFT (k) whether between pre_FFT
Between the mould of (k) and post_FFT (k), and it is adjusted.
If( fabs(gap_FFT(k)) < min( fabs(pre_FFT(k)), fabs(post_FFT(k))) ||
fabs(gap_FFT(k)) > max( fabs(pre_FFT(k)), fabs(post_FFT(k))))
gap_adjusted_FFT(k) =
gap_FFT(k)/ fabs(gap_FFT(k))*( fabs(pre_FFT(k))+ fabs(post_FFT(k)));
else
gap_adjusted_FFT(k) = gap_FFT(k);
Wherein, 0<k <fft_num;Fft_num represents fft variation lengths；Operator fabs () is represented and takes the modulus of complex number.
Further, in present embodiment, the step M3 is further comprising the steps of：
P3, adding window is done to the interpolation sequence after step P2 fill up, obtain the time-domain signal in the insertion gap.
In present embodiment, for the smooth transition at two ends, it is necessary to signal windowing process, certainly, the selection of window
Can be with diversity, such as：Shown in Figure 10, in this way, to gap_adjusted_FFT inverse transformations, obtaining being eventually used for interpolation space
Time-domain signal sig_gap_Result, is not described in detail herein.
With reference to shown in Figure 11 A, 11B, the spectrogram on the left of screen is the U.S. Patent No. by background technology
The spectrogram that the technical scheme of US5891036 is obtained.Screen right side is the spectrogram obtained by method shown in the application.Pass through
The contrast of Figure 11 A understands：The high frequency pseudo of the spectrogram obtained by method shown in the application is less.
By the contrast of Figure 11 B：By method shown in the application obtain spectrogram it is fast-changing in frequency spectrum
Place is more continuous.
To sum up, ultrasonic system Doppler signal gap interpolation method of the invention, using real-time adaptive Doppler's interpolation
Scheme obtains the interpolation sequence in insertion gap；Further, frequency domain correction is carried out to the interpolation sequence, between obtaining described in insertion
The time-domain signal of gap；Make the high frequency pseudo of the spectrogram of acquisition less, and it is more continuous in the fast-changing place of frequency spectrum, carry
The interpolation quality of Doppler frequency spectrum high.
A kind of ultrasonic system Doppler signal gap interpolating apparatus are disclosed with reference to shown in Fig. 7, in an embodiment of the present invention,
Described device includes：Signal acquisition module 100, interpolating module 200, frequency-domain correction module 300；Interpolating module 200 includes：Value
Unit 201, phase place is corrected unit 203, frequency correcting unit 205, bandwidth correcting unit 207, and synthesis unit 209；Frequency domain
Correction module 300 includes noise superpositing unit 301, frequency-domain correction unit 303, and windowing unit 305.
Signal acquisition module 100 is used to obtain two sections of adjacent Doppler's scanning signals；
The gap that interpolating module 200 is used between two sections of Doppler's scanning signals carries out interpolation calculation, obtains insertion institute
State the interpolation sequence in gap.
Although it is understood that the application is with double synchronous modes（B +PW）As a example by do specific introduction, but the application sweeps
The fruit that comes to an end is obtained based on two sections of Doppler's scanning signals, therefore, can also be by B between two sections of adjacent Doppler's scanning signals
Patten transformation is that color mode or B-mode variegate the forms such as pattern, and the application is not between two sections of Doppler's scanning signals
Scanning pattern be specifically limited.
Frequency-domain correction module 300 is used to carry out the interpolation sequence frequency domain correction, obtains the time domain in the insertion gap
Signal.
In present embodiment, describe for convenience, before adjacent two ends Doppler's scanning signal is respectively defined as into gap
Doppler signal Pre_Gap, and Doppler signal Post_Gap behind gap.
It should be noted that describing for convenience, in the following embodiments, the explanation to each alphabetical noun is general
Occur once, but the explanation, it is applicable with later description content, it is not described in detail herein.
Further, value unit 201 is used to carry out two sections of Doppler's scanning signals value, complementary series to be inserted before obtaining
With rear complementary series to be checked.
It should be noted that in the following embodiments, for convenience, current complementary series to be inserted treats interpolation sequence with after
Row occur simultaneously when, forward and backward complementary series to be inserted is referred to as.
In present embodiment, value unit 201 is used to set the repeat its transmission frequency of Doppler's scanning；
Between obtaining between every section of time cycle of Doppler's scanning, and two sections of time cycles of Doppler's scanning of acquisition
The gap scanning time cycle；During scanning, every section of time cycle all same of Doppler's scanning；
According to repeat its transmission frequency, the time cycle of above-mentioned Doppler's scanning, and gap scanning time cycle, obtain every
Interpolation line number in section Doppler's number of scanning line, and the gap；In present embodiment：Calculated using following algorithm.
Assuming that the time cycle of every section of Doppler's scanning is：T_pw, the gap scanning time cycle is T_gap, how general
Strangle scanning repeat its transmission frequency be：prf；
So, every section of Doppler's number of scanning line：N_pw = T_pw*prf；
Interpolation line number in the gap：N_gap = T_gap*prf.
Further, value unit 201 is used for according to every section of Doppler's number of scanning line of acquisition, and inserting in the gap
Value line number, takes the conjugation of the Doppler's scanning signal behind the Doppler's scanning signal before the gap and the gap, before obtaining
Complementary series to be inserted and rear complementary series to be checked.
In present embodiment：Calculated using following algorithm.
Pre_Sig_Conj(k) =conj( Pre_Gap(N_gap-k));k<= N_inter;
Post _Sig_ Conj (k) =conj(Post _Gap(N_inter - k)); k<= N_inter;
Wherein, Pre_Sig_Conj (k) represents preamble interpolation sequence；Post _ Sig_ Conj (k) represents follow-up interpolation
Sequence；N_inter=min (N_pw, N_gap), i.e. N_inter are both N_gap and N_pw smaller value, conj（.）It is right to represent
Value inside bracket takes conjugation.
Further, phase place is corrected unit 203 is used to carry out the preceding complementary series to be inserted and rear complementary series to be checked respectively
Phase place is corrected.
In present embodiment, phase place is corrected unit 203 calculates preamble interpolation sequence and follow-up interpolation sequence to be needed to correct angle
Degree.
In present embodiment：Calculated using following algorithm.
Angle_pre = angle(Pre_Gap(N_gap)) - angle (Pre_Sig_ Conj (1));
Angle_post = angle (Post _Gap(1) ) - angle (Post_Sig_ Conj (N_
inter));
Wherein, Angle_pre represents the angle of preceding complementary series correction to be inserted；Complementary series correction to be inserted after Angle_post
Angle, operator angle () represents the phase angle for taking bracket the inside value..
Further, phase place is corrected unit 203 is multiplied by the list of RA to preceding sequence to be inserted and rear sequence to be inserted respectively
Bit vector, carries out phase place is corrected, obtains new complementary series forward and backward to be inserted.
In present embodiment：Calculated using following algorithm.
Pre_Sig_ForGap(k) = Pre_Sig_ Conj (k)*exp(Angle_pre *i);
Post _Sig_ForGap(k) = Post _Sig_ Conj (k)*exp(Angle_ post *i)
Wherein, Pre_Sig_ForGap (k) is the complementary series preceding to be inserted after phase place is corrected；Post _Sig_ForGap(k)
Complementary series rear to be inserted after for phase place is corrected；I represents imaginary unit.
Further, frequency correcting unit 205 is used to enter line frequency to the preceding complementary series to be inserted and rear complementary series to be checked
Correction；
In present embodiment, frequency correcting unit 205 is used for by treating interpolation before described after phase place is corrected module 203
Line frequency correction is entered in sequence and rear complementary series to be checked.
Frequency correcting unit 205 is carried out to general Le signal Post_Gap behind Doppler signal Pre_Gap before gap and gap
FFT, calculates its FFT result, and the complementary series frequency spectrum forward and backward to be inserted of acquisition is respectively：Pre_fft and post_fft；And root
Obtaining the corresponding centre frequency of forward and backward complementary series to be inserted according to pre_fft and post_fft is：Pre_f0, post_f0；With a width of：
Pre_B and post_B；
In present embodiment, frequency correcting unit 205 is keep frequency smooth transition, forward and backward sequence to be filled up need by
Time sequencing is altered in steps its centre frequency successively, while in view of the time slot between Doppler's scanning signal and gap,
It is determined that preceding complementary series frequency compensation scope to be inserted is [0, (post_f0-pre_f0）* T_gap/T_pw], rear complementary series to be inserted
Frequency compensation scope is：[-(post_f0-pre_f0）*T_gap/T_pw, 0 ].
Further, frequency correcting unit 205 determines forward and backward complementary series compensation to be inserted according to the centre frequency for obtaining
Frequency sequence；And the benefit sequence forward and backward to be checked after frequency correction is obtained according to the forward and backward complementary series compensation frequency sequence to be inserted
Row.
In present embodiment：Frequency correcting unit 205 is calculated respectively to forward and backward complementary series to be inserted using following algorithm
Enter line frequency adjustment：
Pre_fre_Diff(k) = k*(post_f0-pre_f0）*T_gap/T_pw/ N_gap；
Pre_Sig_ForGap(k) = Pre_Sig_ ForGap (k)*exp(k*Pre_fre_Diff(k) *i);
Post_fre_Diff(k) =(N_gap-k)*( pre_f0 - post_f0）*T_gap/T_pw/ N_gap;
Post _Sig_ForGap(k)=Post_Sig_ ForGap(k)*exp((N_gap-k)*Post_fre_Diff
(k)*i)；
Wherein, 0< k<= N_inter；Pre_fre_Diff is preceding complementary series compensation frequency sequence to be inserted, Post_fre_
Diff is rear complementary series compensation frequency sequence to be inserted；Pre_Sig_ForGap (k) is the complementary series preceding to be inserted after frequency correction；
Post _ Sig_ForGap (k) is the complementary series rear to be inserted after frequency correction；
Further, bandwidth correcting unit 207 is used to carry out bandwidth to the preceding complementary series to be inserted and rear complementary series to be checked
Correction.
In present embodiment, bandwidth correcting unit 207 is used for by treating interpolation before described after phase frequency module 205
Line frequency correction is entered in sequence and rear complementary series to be checked.
Bandwidth correcting unit 207 keeps consistent with the bandwidth of neighboring Doppler data respectively in order to ensure complementary series to be inserted
Property, forward and backward complementary series to be inserted is filtered accordingly respectively.
In present embodiment, bandwidth correcting unit 207 can be carried out to forward and backward complementary series to be inserted respectively by setting wave filter
Corresponding filtering.
The wave filter can be IIR or FIR types.
In present embodiment：Calculated using following algorithm.
It is to preceding sequence filter design requirement to be inserted：
Centre frequency is pre_filter_f0=cof1*pre_f0;
Bandwidth requirement is pre_filter_B=cof2*pre_B;
It is to rear sequence filter design requirement to be inserted：
Centre frequency is post_filter_f0=cof3*post_f0;
Bandwidth requirement is post_filter_B=cof4*post_B;
Wherein, cof1, cof2, cof3, cof4 are coefficient, and its value is defaulted as 1.0；
In present embodiment, after bandwidth correcting unit 207 is filtered accordingly to forward and backward complementary series to be inserted, obtain
Pre_Sig_ForGap_filtered and post _ Sig_ForGap_filtered, Pre_Sig_ForGap_filtered are
Complementary series preceding to be inserted after bandwidth correction；Post _ Sig_ForGap_filtered is the complementary series rear to be inserted after bandwidth correction.
It should be noted that in present embodiment, phase place is corrected unit 203, frequency correcting unit 205, bandwidth correction 207
Priority execution sequence might not position said sequence, as long as ensure bandwidth correcting unit performed after frequency correcting unit,
The present invention can be completed, is not described in detail herein.
Further, synthesis unit 209 is used to be rectified by phase place is corrected unit 203, frequency correcting unit 205, bandwidth
Positive 207 complementary series preceding to be inserted is added after distinguishing adding window with rear complementary series to be checked, obtains interpolation sequence.
To ensure progressively transition, synthesis unit 209 is to Pre_Sig_ForGap_filtered and post _ Sig_
ForGap_filtered is added after adding window obtains interpolation sequence respectively：Sig_For_Gap.
In present embodiment, if N_inter<N_gap, synthesis unit 209 is needed in Pre_Sig_ForGap_
0 is mended after filtered makes its length be equal to N_gap, same, it is necessary to 0 is mended before post _ Sig_ForGap_filtered makes it
Length is equal to N_Gap.
In present embodiment：Calculated using following algorithm.
Sig_For_Gap(k) =
Pre_Sig_ForGap_filtered(k)*win(k) + (1-win(k))* post _Sig_ForGap_
filtered(k);
Wherein, Win is the coefficient of institute's adding window, for example：Can be oblique angle window win (k)=1- k/N_Gap；Certainly,
The shape of other gradually transition is can be designed as, is not described in detail herein.
Further, the interpolation sequence of frequency-domain correction module 300 pairs carries out frequency domain correction, obtains the insertion gap
Time-domain signal.
In present embodiment, noise superpositing unit 301 is used for the interpolation its superimposition random phase noise.
Because the more front and rear Doppler signal bandwidth of fir filter interpolator bandwidth may be narrower, noise superpositing unit 301 is holding
The uniformity of the background white-noise level of frequency spectrum is, it is necessary to the corresponding random phase noise of Signal averaging.
In present embodiment：Noise superpositing unit 301 is calculated using following algorithm.
Sig_For_Gap(k) = Sig_For_Gap(k)*exp(Phase_random(k)*i);
Wherein, Phase_random is random phase, and it can be designed as being evenly distributed, and is also designed to Gauss and waits other
Distribution form, is not described in detail herein.
Further, frequency-domain correction unit 303 is used to do the interpolation sequence after step P1 frequency domain correction；
In present embodiment：Frequency-domain correction unit 303 is calculated using following algorithm.
With data grid technology in gap, length is taken for N_pw data are used to calculate FFT, wherein the N_pw data for obtaining are used
Sig_gap_ForFFT is represented.
In this way, record a demerit crossing coincidence length for N_Overlap=(N_pw-N_gap)/2;
Sig_gap_ForFFT ( k) = Pre_Gap(N_pw- N_ Overlap +k);k <= N_ Overlap;
Sig_gap_ForFFT ( k) = Sig_For_Gap(k - N_ Overlap);
N_ Overlap< k <= N_ Overlap+N_gap;
Sig_gap_ForFFT ( k) = Post _Gap (k- N_ Overlap- N_gap);
N_ Overlap+N_gap < k <= N_pw; Overlap+N_gap
Further, frequency-domain correction unit 303 does FFT to Sig_gap_ForFFT and obtains frequency spectrum gap_FFT, according to
Pre_FFT and post_FFT are adjusted to gap_FFT, obtain gap_adjusted_FFT.
Further, frequency-domain correction unit 303 be additionally operable to judge gap_FFT (k) mould whether between pre_FFT (k) and
Between the mould of post_FFT (k), if so, do the interpolation sequence after frequency domain is corrected and do adding window by windowing unit 305 pairs fill up,
Obtain the time-domain signal in the insertion gap；If it is not, then the amplitude to gap_FFT (k) is adjusted；Wait to be adjusted between
After between pre_FFT (k) and the mould of post_FFT (k), then the interpolation sequence after frequency domain is corrected is done by windowing unit 305 pairs
Do adding window to fill up, obtain the time-domain signal in the insertion gap.
Further, in present embodiment, frequency-domain correction unit 303 enters line amplitude adjustment ratio to gap_FFT (k) can be with
Using nonlinear way, simple linear mode can be also used.
In present embodiment, further, what frequency-domain correction unit 303 set following program segments judgement gap_FFT (k) is
It is no between pre_FFT (k) and the mould of post_FFT (k), and it is adjusted.
If( fabs(gap_FFT(k)) < min( fabs(pre_FFT(k)), fabs(post_FFT(k))) ||
fabs(gap_FFT(k)) > max( fabs(pre_FFT(k)), fabs(post_FFT(k))))
gap_adjusted_FFT(k) = gap_FFT(k)/ fabs(gap_FFT(k))*( fabs(pre_FFT(k))
+ fabs(post_FFT(k)));
else
gap_adjusted_FFT(k) = gap_FFT(k);
Wherein, 0<k <fft_num;Fft_num represents fft variation lengths；Operator fabs () is represented and takes the modulus of complex number.
Further, in present embodiment, windowing unit 305 is used for by the interpolation sequence after frequency-domain correction unit 303
Row do adding window and fill up, and obtain the time-domain signal in the insertion gap.
In present embodiment, for the smooth transition at two ends, windowing unit 305 is needed to signal windowing process, certainly,
The selection of window can also diversity, for example：Shown in Figure 10, in this way, to gap_adjusted_FFT inverse transformations, being eventually used for
The time-domain signal sig_gap_Result in interpolation space, is not described in detail herein.
With reference to shown in Figure 11 A, 11B, the spectrogram on the left of screen is the U.S. Patent No. by background technology
The spectrogram that the technical scheme of US5891036 is obtained.Screen right side is the spectrogram obtained by method shown in the application.Pass through
The contrast of Figure 11 A understands：The high frequency pseudo of the spectrogram obtained by the application shown device is less.
By the contrast of Figure 11 B：The spectrogram obtained by the application shown device it is fast-changing in frequency spectrum
Place is more continuous.
In sum, ultrasonic system Doppler signal gap interpolation method of the invention and device, using real-time adaptive
Doppler's interpolation scheme obtains the interpolation sequence in insertion gap；Further, frequency domain correction is carried out to the interpolation sequence, is obtained
Insert the time-domain signal in the gap；Make the high frequency pseudo of the spectrogram of acquisition less, and in the fast-changing place of frequency spectrum
It is more continuous, improve the interpolation quality of Doppler frequency spectrum.
For convenience of description, it is divided into various modules with function during description apparatus above to describe respectively.Certainly, this is being implemented
The function of each module can be realized in same or multiple softwares and/or hardware during application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Realized by the mode of software plus required general hardware platform.Based on such understanding, the technical scheme essence of the application
On the part that is contributed to prior art in other words can be embodied in the form of software product, the computer software product
Can be stored in preservation medium, such as ROM/RAM, magnetic disc, CD, including some instructions are used to so that a computer equipment
（Can be personal computer, Information Push Server, or network equipment etc.）Perform each implementation method of the application or reality
Apply the method described in some parts of mode.
Device embodiments described above are only schematical, wherein the module illustrated as separating component
Can be or may not be physically separate, the part shown as module can be or may not be physics mould
Block, you can with positioned at a place, or can also be distributed on multiple mixed-media network modules mixed-medias.It can according to the actual needs be selected
In some or all of module realize the purpose of present embodiment scheme.Those of ordinary skill in the art are not paying creation
Property work in the case of, you can to understand and implement.
The application can be used in numerous general or special purpose computing system environments or configuration.For example：Personal computer, information
Push server computer, handheld device or portable set, laptop device, multi-processing module system, based on microprocessor mould
The system of block, set top box, programmable consumer-elcetronics devices, network PC, minicom, mainframe computer, including taking up an official post
DCE of what system or equipment etc..
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure etc..The application can also be in a distributed computing environment put into practice, in these DCEs, by
Remote processing devices connected by communication network perform task.In a distributed computing environment, program module can be with
In medium being preserved positioned at the local and remote computer including preserving including equipment.
It should be understood that, although the present specification is described in terms of embodiments, but not each implementation method only includes one
Individual independent technical scheme, this narrating mode of specification is only that for clarity, those skilled in the art will should say
Used as an entirety, technical scheme in each implementation method can also be through appropriately combined, and forming those skilled in the art can for bright book
With the other embodiment for understanding.
Those listed above is a series of to be described in detail only for feasibility implementation method of the invention specifically
Bright, they simultaneously are not used to limit the scope of the invention, all equivalent implementations made without departing from skill spirit of the present invention
Or change should be included within the scope of the present invention.
Claims (6)
1. a kind of ultrasonic system Doppler signal gap interpolation method, it is characterised in that the described method comprises the following steps：
M1, adjacent two sections of Doppler's scanning signals are obtained, Doppler signal behind Doppler signal and gap respectively before gap；
S1, two sections of Doppler's scanning signals are carried out with value, obtain preceding complementary series to be inserted and rear complementary series to be inserted；
S2, the preceding complementary series to be inserted and rear complementary series to be inserted are carried out respectively phase place is corrected, frequency correction, bandwidth correction；
S3, step S2 is corrected after preceding complementary series to be inserted be added after rear complementary series to be inserted difference adding window, obtain interpolation sequence；
P1, to the interpolation its superimposition random phase noise；
P2, the interpolation sequence after step P1 is done frequency domain correction；
Wherein, the frequency domain correction in the step P2 is specifically included：
According to forward and backward complementary series to be inserted, and the interpolation sequence in the gap is inserted, obtain the gap in the gap frequently
Spectrum；
Judge the mould of the gap frequency spectrum whether between the mould of forward and backward complementary series frequency spectrum to be inserted；
If so, continuing step P3；
If it is not, be adjusted to the amplitude of the gap frequency spectrum, treat the mould of the gap frequency spectrum between forward and backward complementary series to be inserted frequently
After between the mould of spectrum, then carry out step P3；
P3, adding window is done to the interpolation sequence after step P2 fill up, obtain the time-domain signal in the insertion gap.
2. ultrasonic system Doppler signal gap interpolation method according to claim 1, it is characterised in that
The frequency correction is specifically included：FFT is carried out to general Le signal behind Doppler signal before gap and gap, it is calculated
FFT result；
The corresponding centre frequency of forward and backward complementary series to be inserted and bandwidth are obtained according to the FFT result；
Complementary series frequency compensation scope to be inserted before being determined according to the centre frequency for obtaining, and forward and backward complementary series to be inserted is mended
Repay frequency sequence；
Complementary series forward and backward to be inserted after frequency correction is obtained according to the forward and backward complementary series compensation frequency sequence to be inserted.
3. ultrasonic system Doppler signal gap interpolation method according to claim 2, it is characterised in that
The bandwidth correction is specifically included：
By setting the filtering that wave filter is carried out to forward and backward complementary series to be inserted respectively；
The wave filter is IIR types or FIR types.
4. a kind of ultrasonic system Doppler signal gap interpolating apparatus, it is characterised in that described device includes：
Signal acquisition module, for obtaining two sections of adjacent Doppler's scanning signals, respectively before gap Doppler signal and
Doppler signal after gap；
The interpolating module is used for：Two sections of Doppler's scanning signals are carried out with value, complementary series to be inserted treats interpolation with after before obtaining
Sequence；
Carry out phase place is corrected, frequency correction, bandwidth correction respectively to the preceding complementary series to be inserted and rear complementary series to be inserted；
It is added after preceding complementary series to be inserted after above-mentioned correction is distinguished into adding window with rear complementary series to be inserted, obtains interpolation sequence；
Frequency-domain correction module, for the interpolation its superimposition random phase noise；Inserting after to superposition random phase noise
Frequency domain correction is done in complementary series；Do adding window and fill up to doing the interpolation sequence after frequency domain correction, obtain the time domain letter in the insertion gap
Number；
Wherein, the frequency-domain correction module be additionally operable to when frequency domain is corrected：
According to forward and backward complementary series to be inserted, and the interpolation sequence in the gap is inserted, obtain the gap in the gap frequently
Spectrum；
Judge the mould of the gap frequency spectrum whether between the mould of forward and backward complementary series frequency spectrum to be inserted；
If so, doing adding window and filling up to doing the interpolation sequence after frequency domain correction, the time-domain signal in the insertion gap is obtained；
If it is not, be adjusted to the amplitude of the gap frequency spectrum, treat the mould of the gap frequency spectrum between forward and backward complementary series to be inserted frequently
After between the mould of spectrum, then do adding window and fill up to doing the interpolation sequence after frequency domain correction, obtain the time-domain signal in the insertion gap.
5. ultrasonic system Doppler signal gap interpolating apparatus according to claim 4, it is characterised in that
The interpolating module enter line frequency correct when specifically for：
FFT is carried out to general Le signal behind Doppler signal before gap and gap, its FFT result is calculated；
The corresponding centre frequency of forward and backward complementary series to be inserted and bandwidth are obtained according to the FFT result；
Complementary series frequency compensation scope to be inserted before being determined according to the centre frequency for obtaining, and forward and backward complementary series to be inserted is mended
Repay frequency sequence；
Complementary series forward and backward to be inserted after frequency correction is obtained according to the forward and backward complementary series compensation frequency sequence to be inserted.
6. ultrasonic system Doppler signal gap interpolating apparatus according to claim 5, it is characterised in that
The interpolating module carry out bandwidth correct when specifically for：
By setting the filtering that wave filter is carried out to forward and backward complementary series to be inserted respectively；
The wave filter is IIR types or FIR types.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510140766.4A CN104783836B (en) | 2015-03-27 | 2015-03-27 | Ultrasonic system Doppler signal gap interpolation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510140766.4A CN104783836B (en) | 2015-03-27 | 2015-03-27 | Ultrasonic system Doppler signal gap interpolation method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104783836A CN104783836A (en) | 2015-07-22 |
| CN104783836B true CN104783836B (en) | 2017-06-23 |
Family
ID=53549358
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510140766.4A Active CN104783836B (en) | 2015-03-27 | 2015-03-27 | Ultrasonic system Doppler signal gap interpolation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104783836B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3340887B1 (en) * | 2015-08-27 | 2020-06-17 | Koninklijke Philips N.V. | Spectral doppler processing with adaptive sample window size |
| CN110141270B (en) * | 2019-06-21 | 2022-01-28 | 青岛海信医疗设备股份有限公司 | Beam synthesis method and device |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5016641A (en) * | 1989-11-13 | 1991-05-21 | Advanced Technology Laboratories, Inc. | Spectral interpolation of ultrasound Doppler signal |
| US5891036A (en) * | 1996-09-11 | 1999-04-06 | Hitachi Medical Corporation | Ultrasonic wave Doppler diagnosing apparatus |
| CN1915174A (en) * | 2005-08-16 | 2007-02-21 | 深圳迈瑞生物医疗电子股份有限公司 | Method for processing gaps of Doppler signal |
| CN101336830A (en) * | 2007-07-03 | 2009-01-07 | 深圳迈瑞生物医疗电子股份有限公司 | Orthogonal Doppler signal gap filling method and device for ultrasonic dignosis imaging |
-
2015
- 2015-03-27 CN CN201510140766.4A patent/CN104783836B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5016641A (en) * | 1989-11-13 | 1991-05-21 | Advanced Technology Laboratories, Inc. | Spectral interpolation of ultrasound Doppler signal |
| US5891036A (en) * | 1996-09-11 | 1999-04-06 | Hitachi Medical Corporation | Ultrasonic wave Doppler diagnosing apparatus |
| CN1915174A (en) * | 2005-08-16 | 2007-02-21 | 深圳迈瑞生物医疗电子股份有限公司 | Method for processing gaps of Doppler signal |
| CN101336830A (en) * | 2007-07-03 | 2009-01-07 | 深圳迈瑞生物医疗电子股份有限公司 | Orthogonal Doppler signal gap filling method and device for ultrasonic dignosis imaging |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104783836A (en) | 2015-07-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20240206851A1 (en) | Filtering method and system for ultrasonic imaging | |
| CN104783836B (en) | Ultrasonic system Doppler signal gap interpolation method and system | |
| JP2005081151A (en) | Tracking clutter filter for spectral and audio doppler | |
| JP2001286472A (en) | Ultrasonic diagnostic apparatus | |
| CN108535613B (en) | Voltage flicker parameter detection method based on combined window function | |
| CN103156637B (en) | Ultrasound volume image data processing method and equipment | |
| Alessandrini et al. | Simulation of realistic echocardiographic sequences for ground-truth validation of motion estimation | |
| CN102727257A (en) | Puncture needle visualized method and puncture needle visualized device | |
| CN114241078A (en) | Accelerated magnetic resonance super-resolution imaging method, device, equipment and medium | |
| US8162838B2 (en) | Doppler signal spectrum calculating method and apparatus | |
| JP2004337457A (en) | Ultrasonic image processor | |
| CN110840484B (en) | Ultrasonic imaging method and device for adaptively matching optimal sound velocity and ultrasonic equipment | |
| CN113297987B (en) | Variational modal decomposition signal noise reduction method based on dual-objective function optimization | |
| CN111458122A (en) | Rotary machine fault diagnosis method based on matching enhancement time-frequency representation | |
| CN114791597A (en) | Ultra-short baseline positioning method and device based on three-dimensional array and electronic equipment | |
| CN110680398A (en) | Ultrasonic echo imaging method and system and ultrasonic equipment | |
| WO2016066497A1 (en) | Method of visualizing a sequence of ultrasound images, computer program product and ultrasound system | |
| JPH0616782B2 (en) | Method and apparatus for combining continuous estimated signals in real time | |
| CN105708495B (en) | Ultrasonic color blood flow imaging boundary processing method and system | |
| CN108254607A (en) | Digital oscilloscope method for displaying waveform and device | |
| WO2021036641A1 (en) | Coupling mismatch detection method and related device | |
| CN106725615A (en) | A kind of ivus image pulse interference suppression method and device | |
| US20230366998A1 (en) | Signal processing methods and systems | |
| WO2014004052A1 (en) | Ultrasound diagnostic apparatus and method for generating doppler spectrum signal | |
| CN105212963B (en) | Realize the wall filtering method and system of ultrasonic color doppler blood flow imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province Patentee after: Feiyinuo Technology Co.,Ltd. Address before: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province Patentee before: Feiyinuo Technology (Suzhou) Co.,Ltd. |
|
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province Patentee after: Feiyinuo Technology (Suzhou) Co.,Ltd. Address before: 215123 unit 501, building C8, bio nano Park, No. 218, Xinghu street, industrial park, Suzhou, Jiangsu Province Patentee before: VINNO TECHNOLOGY (SUZHOU) Co.,Ltd. |