CN104783836A - Interval interpolation method and system for Doppler signals of ultrasonic system - Google Patents
Interval interpolation method and system for Doppler signals of ultrasonic system Download PDFInfo
- Publication number
- CN104783836A CN104783836A CN201510140766.4A CN201510140766A CN104783836A CN 104783836 A CN104783836 A CN 104783836A CN 201510140766 A CN201510140766 A CN 201510140766A CN 104783836 A CN104783836 A CN 104783836A
- 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.)
- Granted
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention provides an interval interpolation method and system for Doppler signals of an ultrasonic system. The method comprises the steps that two adjacent Doppler scanning signals are obtained; interpolation calculation is conducted on the interval between the two adjacent Doppler scanning signals, and an interpolation sequence inserted in the interval is obtained; frequency domain correction is conducted on the interpolation sequence, and time domain signals inserted in the interval are obtained. According to the interval interpolation method and device for Doppler signals, the interpolation sequence inserted in the interval is obtained through a real-time self-adaptive Doppler interpolation scheme; furthermore, the frequency domain correction is conducted on the interpolation sequence, and the time domain signals inserted in the interval are obtained; an obtained spectrogram has less high-frequency pseudo noise, the continuity at the places with high frequency spectrum change is better, and the interpolation quality of the Doppler frequency spectrum is improved.
Description
Technical field
The invention belongs to ultrasonic diagnostic imaging technical field, relate to a kind of ultrasonic system Doppler signal gap interpolation method and system.
Background technology
Diasonograph (B ultrasonic machine) has exclusive characteristic, as real-time; Wherein, dynamic blood flow conventional in diasonograph and doppler imaging, become one of means of the indispensable auxiliary diagnosis of modern medicine; Doppler imaging can measure the concrete flow velocity of blood flow quickly and easily, becomes the criterion of some disease in clinical diagnosis.
In Doppler's (being called for short PW) checking process, usually adopt two kinds of patterns: i.e. single synchronous mode and many synchronous modes; Described single synchronous mode is comparatively simple, and whole scanning process only upgrades doppler data according to certain speed (being called for short PRF in Doppler's scanning), and other patterns, such as B-mode, color mode are in stopping frozen state; When needs adopt other pattern to scan, other patterns of doctor's active activation, are not described in detail at this.But single synchronous mode exists some defects, such as: in Doppler's scanning process, need to keep probe to be in stationary state.
Many synchronous scanning patterns solve appeal problem, in whole scanning process, the scanning in succession of B-mode, color mode and doppler mode, doctor can observe the image of multiple pattern simultaneously, but, also there are some problems in the Doppler's scanning under multi-mode, wherein, topmost problem is: multiple pattern, in succession after scanning, causes there is larger space between doppler data, the 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 prior art, have employed kinds of schemes and solve the problem.
As shown in Figure 1, Chinese patent application, publication number: CN101336830, discloses a kind of orthogonal Doppler signal gap filling method for ultrasonic diagnostic imaging and device.
The method by before and after gap signal positive sequence fill gap, and get before and after gap each one section as intermediate zone, filling signal and original signal weighted superposition (weight coefficient sum is 1) in intermediate zone.The frequency parameter of filling front and back end signal is obtained by auto-correlation function method, signal averaging frequency, signal power and signal bandwidth, when meeting certain threshold value, interpolated signal is multiplied by the complex signal of a frequency, discontinuous with the frequency spectrum that interpolation when signal is larger before and after compensate for clearance causes.But the method also exists drawback, it can not solve the problem at sideband signal phase hit; Meanwhile, the method is just multiplied by the multiple sinusoidal signal of front and back difference on the frequency simply to interpolation signal, do not consider frequency change transition in time, and the impact of the time width in gap, cause frequency to have sudden change, form high frequency pseudomorphism, cause on the spectrogram shown in Fig. 1 and have obvious vertical line.
In addition, the background technology part of this Chinese patent application also refer to U.S. Patent number: US5016641, US5476097, US5891036, US5642732, US4559952, and the prior art such as U.S. Patent Application No. US2007/0049823.But the above-mentioned prior art quoted, carries out in Interpolation Process, all there will be SPA sudden phase anomalies problem, frequency discontinuity problem between to two Doppler's scanning signals, at least one in amplitude hit problem and real time problems.
As shown in Figure 2, U.S. Patent number: US5016641, name is called: the technical scheme of " Spectral interpolation of ultrasound doppler signal ", it utilizes the power spectrum of the signal before gap and a random phase of generation, obtain time-domain signal by IFFT, this time-domain signal obtains not having the filling signal of saltus step after windowing process, the addition that partly overlaps.The generation method of random phase adds Gaussian distributed on the basis of original signal phase place, and variance is the random phase of (PI/10, PI/2).On the spectrogram obtained by this technical scheme, because its interpolation method is in the Doppler signal place adjacent with gap, may there is sudden change, so there is obvious vertical line, sound has " Dth Dth " sound simultaneously.
As shown in Figure 3, U.S. Patent number: US5476097, name is called: the technical scheme of " Simultaneous Ultrasonic Imaging and Doppler Display System ", it adopts frequency spectrum harmony cent to open filling, wherein acoustical signal directly splices by the signal inverted sequence at two ends, gap, ensureing the continuous of two ends, gap, but there is saltus step in the middle of splicing, this discontinuous noise caused can be eliminated by carrying out low-pass filtering after digital-to-analogue conversion again.On the spectrogram obtained by this technical scheme, owing to only utilizing spectrogram interpolation, therefore the spectrogram obtained is very fuzzy.
As shown in Figure 4, U.S. Patent Application No.: US2007/0049823, name is called: the technical scheme of " Method for Processing Doppler Signal Gaps ", it adopts and carries out high-pass filtering to data before and after gap, inverted sequence conjugation weighting summation obtains signal in gap, this method can ensure the seriality at stitching portion signal, acoustically do not suddenly change, but as shown in Figure 4, on the frequency spectrum obtained, there are some thinner lines by the method, therefore not sudden change can not be ensured equally.
As shown in Figure 5, U.S. Patent number: US5891036, name be called: the technical scheme of " Ultrasonic wave Doppler diagnosing apparatus ", it adopts the data inverted sequence before gap to read, get conjugation or real part and imaginary part exchange after stored in interstitial site.Discontinuous owing to getting conjugation post gap starting point generation phase place, solve this phase deviation amount, to this phase contrast of compensation data in gap.Second method carries out same operation to the data before and after gap, then by two segment data weighted superposition after process, ensures all continuously smooth transition of two ends, gap.As shown in Figure 5, the frequency spectrum adopting the method to obtain does not have obvious vertical line, and just relatively former frequency spectrum a little some " void ", and particularly in the place that frequency change is larger, rising portion frequency spectrum as quick in the frequency in figure has the phenomenon of " fracture ".
Further, U.S. Patent number: US5642732, name is called: the technical scheme of " Apparatus and Method for Estimating Missing Doppler Signals and Spectra ", its Doppler signal gap filling method estimated based on autoregression (AR) model, in AR model, signal extension is the natural result of Power estimation, and this Power estimation method is without any need for windowing, can obtain good effect in theory; But its exponent number and model parameter estimation method all can have impact to performance.
Further, U.S. Patent number: US4559952, name is called: the technical scheme of " Method of Ultrasonically Measuring Blood Flow Velocity ", first it try to achieve the spectrum distribution of gap front signal, by broadband noise through a wave filter to produce signal to be filled, the coefficient controlling wave filter produces the spectral characteristic needed, and directly reads back-page Doppler signal, between the two window function windowed overlapping.The method is compared with U.S. Patent number US5016641, U.S. Patent Application No. US2007/0049823, and especially without clear superiority part, it generates filter coefficient in real time according to frequency spectrum in addition, realizes more complicated.
Summary of the invention
Order of the present invention is to provide a kind of ultrasonic system Doppler signal gap interpolation method and device.
One of for achieving the above object, an embodiment of the present invention provides a kind of ultrasonic system Doppler signal gap interpolation method, said method comprising the steps of:
M1, obtain two sections of adjacent Doppler's scanning signals, to be respectively before gap Doppler signal after Doppler signal and gap;
Interpolation calculation is carried out in M2, gap between two sections of Doppler's scanning signals, obtains the interpolation sequence inserting described gap;
M3, frequency domain rectification is carried out to described interpolation sequence, obtain the time-domain signal inserting described gap.
As a further improvement on the present invention, described " interpolation calculation is carried out in the gap between two sections of Doppler's scanning signals, obtains the interpolation sequence inserting described gap " specifically comprises:
S1, value is carried out to two sections of Doppler's scanning signals, obtain front complementary series to be inserted and rear complementary series to be checked;
S2, described front complementary series to be inserted and rear complementary series to be checked carried out respectively to phase place rectification, frequency is corrected, bandwidth is corrected;
S3, step S2 is corrected after front complementary series to be inserted and rear complementary series to be checked be added after windowing respectively, obtain interpolation sequence.
As a further improvement on the present invention, described frequency is corrected and is specifically comprised: carry out FFT conversion to signal of strangling general behind Doppler signal before gap and gap, calculate its FFT result;
Mid frequency corresponding to forward and backward complementary series to be inserted and bandwidth is obtained according to described FFT result;
Described mid frequency according to obtaining determines front complementary series frequency compensation scope to be inserted, and forward and backward complementary series compensating frequency sequence to be inserted;
Complementary series forward and backward to be checked after correcting according to described forward and backward complementary series compensating frequency retrieval frequency to be inserted.
As a further improvement on the present invention, described bandwidth is corrected and is specifically comprised:
By arranging wave filter respectively to the filtering that forward and backward complementary series to be inserted is carried out;
Described wave filter is IIR type or FIR type.
As a further improvement on the present invention, described " carry out frequency domain rectification to described interpolation sequence, obtain the time-domain signal inserting described gap " specifically comprises:
P1, to described interpolation its superimposition random phase noise;
P2, the interpolation sequence after step P1 done to frequency domain and correct;
P3, windowing is done to the interpolation sequence after step P2 fill up, obtain the time-domain signal inserting described gap.
As a further improvement on the present invention, described frequency domain is corrected and is specifically comprised:
According to forward and backward complementary series to be inserted, and insert the described interpolation sequence in described gap, obtain the gap frequency spectrum in described gap;
Judge that the mould of described gap frequency spectrum is whether between the mould of forward and backward complementary series frequency spectrum to be inserted;
If so, step P3 is continued;
If not, after the amplitude of described gap frequency spectrum is adjusted, then step P3 is carried out.
One of for achieving the above object, an embodiment of the present invention provides a kind of ultrasonic system Doppler signal gap interpolating apparatus, described device comprises: signal acquisition module, for obtaining two sections of adjacent Doppler's scanning signals, to be respectively before gap Doppler signal after Doppler signal and gap;
Interpolating module, carries out interpolation calculation for the gap between two sections of Doppler's scanning signals, obtains the interpolation sequence inserting described gap;
Frequency-domain correction module, for carrying out frequency domain rectification to described interpolation sequence, obtains the time-domain signal inserting described gap.
As a further improvement on the present invention, described interpolating module specifically for:
Value is carried out to two sections of Doppler's scanning signals, obtains front complementary series to be inserted and rear complementary series to be checked;
Phase place rectification is carried out respectively to described front complementary series to be inserted and rear complementary series to be checked, frequency is corrected, bandwidth is corrected;
Be added after complementary series to be inserted front after above-mentioned rectification and rear complementary series to be checked respectively windowing, obtain interpolation sequence.
As a further improvement on the present invention, described interpolating module carry out frequency correct time specifically for:
FFT conversion is carried out to signal of strangling general behind Doppler signal before gap and gap, calculates its FFT result;
Mid frequency corresponding to forward and backward complementary series to be inserted and bandwidth is obtained according to described FFT result;
Described mid frequency according to obtaining determines front complementary series frequency compensation scope to be inserted, and forward and backward complementary series compensating frequency sequence to be inserted;
Complementary series forward and backward to be checked after correcting according to described forward and backward complementary series compensating frequency retrieval frequency to be inserted.
As a further improvement on the present invention, described interpolating module carry out bandwidth correct time specifically for:
By arranging wave filter respectively to the filtering that forward and backward complementary series to be inserted is carried out;
Described wave filter is IIR type or FIR type.
As a further improvement on the present invention, described frequency-domain correction module specifically for:
To described interpolation its superimposition random phase noise;
Do frequency domain to the interpolation sequence after superposition random phase noise to correct;
To do frequency domain correct after interpolation sequence do windowing and fill up, obtain the time-domain signal inserting described gap.
As a further improvement on the present invention, described frequency-domain correction module carry out frequency domain correct time also for:
According to forward and backward complementary series to be inserted, and insert the described interpolation sequence in described gap, obtain the gap frequency spectrum in described gap;
Judge that the mould of described gap frequency spectrum is whether between the mould of forward and backward complementary series frequency spectrum to be inserted;
If so, to do frequency domain correct after interpolation sequence do windowing and fill up, obtain the time-domain signal inserting described gap;
If not, after the amplitude of described gap frequency spectrum is adjusted, then to do frequency domain correct after interpolation sequence do windowing and fill up, obtain the time-domain signal inserting described gap.
Compared with prior art, the invention has the beneficial effects as follows: ultrasonic system Doppler signal gap interpolation method of the present invention and device, adopt real-time adaptive Doppler interpolation scheme to obtain the interpolation sequence inserting gap; Further, frequency domain rectification is carried out to described interpolation sequence, obtain the time-domain signal inserting described gap; Make the high frequency pseudo of the spectrogram of acquisition less, and more continuous in the fast-changing place of frequency spectrum, improve the interpolation quality of Doppler frequency spectrum.
Accompanying drawing explanation
Fig. 1 is the spectrogram of the ultrasonic system Doppler signal of the first embodiment in background technology;
Fig. 2 is the spectrogram of the ultrasonic system Doppler signal of the second embodiment in background technology;
Fig. 3 is the spectrogram of the ultrasonic system Doppler signal of the 3rd embodiment in background technology;
Fig. 4 is the spectrogram of the ultrasonic system Doppler signal of the 4th embodiment in background technology;
Fig. 5 is the spectrogram of the ultrasonic system Doppler signal of the 5th embodiment in background technology;
Fig. 6 is ultra sonic 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 the structural representation of a kind of ultrasonic system Doppler signal gap interplotation system in an embodiment of the present invention;
Fig. 9 is Doppler's scanning structural representation of two synchronous mode in an embodiment of the present invention;
Figure 10 is the structural representation in an embodiment of the present invention, Doppler signal being carried out to windowing process;
Figure 11 A is for adopting background technology the 5th embodiment and adopting the high frequency pseudo of the spectrogram that method obtains respectively shown in Fig. 7 to contrast schematic diagram;
Figure 11 B is the contrast schematic diagram adopting background technology the 5th embodiment and adopt the frequency spectrum speed of the spectrogram that method obtains respectively shown in Fig. 7 to change.
Detailed description of the invention
Describe the present invention below with reference to each embodiment shown in the drawings.But these embodiments do not limit the present invention, the structure that those of ordinary skill in the art makes according to these embodiments, method or conversion functionally are all included in protection scope of the present invention.
As shown in Figure 6, ultra sonic imaging flowage structure schematic diagram in ultrasonic system.
Ultrasonic signal in tissue changes electric analoging signal into via the different primitives of probe transducers, is amplified by front amplification module, then is converted to digital signal by A/D D/A converter module; The digital signal of each different primitive, through Beam synthesis module, synthesizes radiofrequency signal; Quadrature demodulation result i/q signal, after quadrature demodulation, is sent into corresponding processing module by radiofrequency signal.
Under many synchronous modes, ultrasonic system utilizes time-multiplexed method, periodically at different time scanning different mode.Below, concrete introduction will be done for two synchronous mode (B+PW).
Shown in composition graphs 7, shown in composition graphs 9, the ultrasonic system Doppler signal gap interpolation method that an embodiment of the present invention provides comprises the following steps:
M1, obtain two sections of adjacent Doppler's scanning signals;
Interpolation calculation is carried out in M2, gap between two sections of Doppler's scanning signals, obtains the interpolation sequence inserting described gap.
Be understandable that, although the application does concrete introduction for two synchronous mode (B+PW), but the scanning result of the application is based on two sections of Doppler's scanning signal acquisitions, therefore, between adjacent two sections of Doppler's scanning signals, B-mode can also be converted to color mode or B-mode variegates the forms such as pattern, the application is not specifically limited the scanning pattern between two sections of Doppler's scanning signals.
M3, frequency domain rectification is carried out to described interpolation sequence, obtain the time-domain signal inserting described gap.
In present embodiment, for convenience of description, adjacent two ends Doppler's scanning signal is defined as respectively Doppler signal Pre_Gap before gap, and Doppler signal Post_Gap behind gap.
It should be noted that, for convenience of description, in the following embodiments, the explanation of each alphabetical noun is generally all occurred once, but this explanation, be all suitable for later description content, be not described in detail at this.
Further, the interpolation sequence obtaining described gap in step M2 specifically comprises:
S1, value is carried out to two sections of Doppler's scanning signals, obtain front complementary series to be inserted and rear complementary series to be checked.
It should be noted that, in the following embodiments, for convenience, when current complementary series to be inserted and rear complementary series to be inserted occur, by it referred to as forward and backward complementary series to be inserted simultaneously.
In present embodiment, described step S1 specifically comprises:
The repeat its transmission frequency of setting Doppler scanning;
Obtain the time cycle of every section of Doppler's scanning, and the gap scanning time cycle between the time cycle obtaining two sections of Doppler's scannings; In scanning process, the time cycle of every section of described Doppler's scanning is all identical;
According to repeat its transmission frequency, the time cycle of above-mentioned Doppler's scanning, and the gap scanning time cycle, obtain every section of Doppler's number of scanning line, and the interpolation line number in described gap; In present embodiment: adopt following algorithm to calculate.
Suppose that the time cycle of every section of Doppler's scanning is: T_pw, the described gap scanning time cycle is T_gap, and the repeat its transmission frequency of Doppler's scanning is: prf;
So, every section of Doppler's number of scanning line: N_pw=T_pw*prf;
Interpolation line number in described gap: N_gap=T_gap*prf.
Further, according to acquisition every section of Doppler's number of scanning line, and the interpolation line number in described gap, get the conjugation of the Doppler's scanning signal behind the Doppler's scanning signal before described gap and described gap, obtain front complementary series to be inserted and rear complementary series to be checked.
In present embodiment: adopt following algorithm to calculate.
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 preorder interpolation sequence; Post _ Sig_ Conj (k) represents follow-up interpolation sequence; N_inter=min (N_pw, N_gap), namely N_inter is N_gap and N_pw smaller value, conj(.) represent and conjugation is got to the value inside bracket.
Further, the interpolation sequence obtaining described gap also comprises:
S2, described front complementary series to be inserted and rear complementary series to be checked carried out respectively to phase place rectification, frequency is corrected, bandwidth is corrected.
In present embodiment, described to described front complementary series to be inserted and rear complementary series to be checked carry out phase place correct specifically comprise:
Calculating preorder interpolation sequence and follow-up interpolation sequence need RA.
In present embodiment: adopt following algorithm to calculate.
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 that front complementary series to be inserted corrects; The angle that after Angle_post, complementary series to be inserted corrects, the phase angle be worth inside bracket is got in operator angle (.) expression.。
Further, respectively front sequence to be inserted and rear sequence to be inserted are multiplied by the unit vector of RA, carry out phase place rectification, obtain new complementary series forward and backward to be inserted.
In present embodiment: adopt following algorithm to calculate.
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) for phase place correct after before complementary series to be inserted; Post _ Sig_ForGap (k) is the complementary series rear to be inserted after phase place rectification; I represents imaginary unit.
In present embodiment, after carrying out phase place rectification to forward and backward complementary series to be inserted, further, carry out frequency rectification to described front complementary series to be inserted and rear complementary series to be checked specifically to comprise:
Carry out FFT conversion to the signal Post_Gap that strangles general behind Doppler signal Pre_Gap before gap and gap, calculate its 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 the mid frequency that pre_fft and post_fft obtains forward and backward complementary series to be inserted corresponding be: pre_f0, post_f0; Bandwidth is: pre_B and post_B;
Described FFT is fast Fourier transform (abbreviation of Fast Fourier Transformation), it is the fast algorithm of discrete fourier transform, it is the characteristic such as odd, even, empty, real according to discrete fourier transform, carries out improving obtaining the algorithm of discrete Fourier transform (DFT).
In present embodiment, for holding frequency smooth transition, forward and backward sequence to be filled up needs progressively to change its mid frequency successively in chronological order, consider the time slot between Doppler's scanning signal and gap simultaneously, determine that front complementary series frequency compensation scope to be inserted 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 compensating frequency sequence to be inserted is determined according to the described mid frequency obtained; And the complementary series forward and backward to be checked after correcting according to described forward and backward complementary series compensating frequency retrieval frequency to be inserted.
In present embodiment: adopt following algorithm to carry out calculating and carry out 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 front complementary series compensating frequency sequence to be inserted, and Post_fre_Diff is rear complementary series compensating frequency sequence to be inserted; Pre_Sig_ForGap (k) for frequency correct after before complementary series to be inserted; Post _ Sig_ForGap (k) is the complementary series rear to be inserted after frequency rectification;
In present embodiment, after carrying out frequency rectification to forward and backward complementary series to be inserted, further, carry out bandwidth rectification to described front complementary series to be inserted and rear complementary series to be checked specifically to comprise:
In order to ensure that complementary series to be inserted keeps the bandwidth concordance with neighboring Doppler data respectively, respectively corresponding filtering is carried out to forward and backward complementary series to be inserted.
In present embodiment, respectively corresponding filtering is carried out to forward and backward complementary series to be inserted by arranging wave filter.
Described wave filter can be IIR or FIR type.
In present embodiment: adopt following algorithm to calculate.
To front sequence filter designing requirement to be inserted be:
Mid frequency is pre_filter_f0=cof1*pre_f0;
Bandwidth requirement is pre_filter_B=cof2*pre_B;
To rear sequence filter designing requirement to be inserted be:
Mid 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 all defaulted as 1.0;
In present embodiment, after carrying out corresponding filtering to forward and backward complementary series to be inserted, acquisition Pre_Sig_ForGap_filtered and post _ Sig_ForGap_filtered, Pre_Sig_ForGap_filtered is the complementary series front to be inserted after bandwidth is corrected; Post _ Sig_ForGap_filtered is the complementary series rear to be inserted after bandwidth is corrected.
It should be noted that, in present embodiment, order according to the order of phase place rectification, frequency rectification, bandwidth rectification, forward and backward complementary series to be inserted might not be corrected, as long as ensure that bandwidth is corrected to complete after frequency is corrected, namely can complete illustrated method, not be described in detail at this.
Further, the interpolation sequence obtaining described gap also comprises:
S3, step S2 is corrected after front complementary series to be inserted and rear complementary series to be checked be added after windowing respectively, obtain interpolation sequence.
For ensureing progressively transition, obtain interpolation sequence to being added after Pre_Sig_ForGap_filtered and post _ Sig_ForGap_filtered respectively windowing: Sig_For_Gap.
In present embodiment, if N_inter< is N_gap, needs after Pre_Sig_ForGap_filtered, mend 0 and make its length equal N_gap, equally, need before post _ Sig_ForGap_filtered, mend 0 and make its length equal N_Gap.
In present embodiment: adopt following algorithm to calculate.
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 windowing, such as: can be oblique angle window win (k)=1-k/N_Gap; Certainly, also can be designed as the shape of other transition gradually, be not described in detail at this.
Further, in step M3, frequency domain rectification is carried out to described interpolation sequence, obtains the time-domain signal inserting described gap and specifically comprise:
P1, to described interpolation its superimposition random phase noise;
In present embodiment, because fir filter interpolator bandwidth may be narrower compared with front and back Doppler signal bandwidth, for keeping the concordance of the background white-noise level of frequency spectrum, need the corresponding random phase noise of Signal averaging.
In present embodiment: adopt following algorithm to calculate.
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 and is evenly distributed, and also can be designed to other distribution forms such as Gauss, not be described in detail at this.
Further, described step M3 also comprises:
P2, the interpolation sequence after step P1 done to frequency domain and correct;
In present embodiment: adopt following algorithm to calculate.
With data grid technology in gap, get length be N_pw data for calculating FFT, the N_pw data Sig_gap_ForFFT wherein obtained represents.
So, record a demerit that to cross coincidence length be 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 conversion is done to Sig_gap_ForFFT and obtains frequency spectrum gap_FFT, according to pre_FFT and post_FFT, gap_FFT is adjusted, obtain gap_adjusted_FFT.
Further, judge that the mould of gap_FFT (k) is whether between pre_FFT (k) and the mould of post_FFT (k), if so, continues step P3, if not, then adjusts the amplitude of gap_FFT (k); Wait being adjusted to after between pre_FFT (k) and the mould of post_FFT (k), then carry out step P3.
Further, in present embodiment, resize ratio can adopt nonlinear way, also can adopt simple linear mode.
In present embodiment, further, set following program segment judge gap_FFT (k) whether between pre_FFT (k) and the mould of post_FFT (k), and it to be 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 length; The modulus of complex number is got in operator fabs () expression.
Further, in present embodiment, described step M3 is further comprising the steps of:
P3, windowing is done to the interpolation sequence after step P2 fill up, obtain the time-domain signal inserting described gap.
In present embodiment, in order in the smooth transition of two ends place, need signal windowing process, certainly, the selection of window also can multiformity, such as: shown in Figure 10, so, to gap_adjusted_FFT inverse transformation, obtain, finally for the time-domain signal sig_gap_Result in interpolation space, not being described in detail at this.
Shown in composition graphs 11A, 11B, the spectrogram on the left of screen is by background technology, and U.S. Patent number is the spectrogram that the technical scheme of US5891036 obtains.Spectrogram for being obtained by method shown in the application on the right side of screen.Known by the contrast of Figure 11 A: the high frequency pseudo of the spectrogram obtained by method shown in the application is less.
Known by the contrast of Figure 11 B: spectrogram more continuous in the fast-changing place of frequency spectrum obtained by method shown in the application.
To sum up, ultrasonic system Doppler signal gap interpolation method of the present invention, adopts real-time adaptive Doppler interpolation scheme to obtain the interpolation sequence inserting gap; Further, frequency domain rectification is carried out to described interpolation sequence, obtain the time-domain signal inserting described gap; Make the high frequency pseudo of the spectrogram of acquisition less, and more continuous in the fast-changing place of frequency spectrum, improve the interpolation quality of Doppler frequency spectrum.
Shown in composition graphs 7, disclose a kind of ultrasonic system Doppler signal gap interpolating apparatus in an embodiment of the present invention, described device comprises: signal acquisition module 100, interpolating module 200, frequency-domain correction module 300; Interpolating module 200 comprises: value unit 201, phase place correcting unit 203, frequency correcting unit 205, bandwidth correcting unit 207, and synthesis unit 209; Frequency-domain correction module 300 comprises noise superpositing unit 301, frequency-domain correction unit 303, and windowing unit 305.
Signal acquisition module 100 is for obtaining two sections of adjacent Doppler's scanning signals;
Interpolating module 200 carries out interpolation calculation for the gap between two sections of Doppler's scanning signals, obtains the interpolation sequence inserting described gap.
Be understandable that, although the application does concrete introduction for two synchronous mode (B+PW), but the scanning result of the application is based on two sections of Doppler's scanning signal acquisitions, therefore, between adjacent two sections of Doppler's scanning signals, B-mode can also be converted to color mode or B-mode variegates the forms such as pattern, the application is not specifically limited the scanning pattern between two sections of Doppler's scanning signals.
Frequency-domain correction module 300, for carrying out frequency domain rectification to described interpolation sequence, obtains the time-domain signal inserting described gap.
In present embodiment, for convenience of description, adjacent two ends Doppler's scanning signal is defined as respectively Doppler signal Pre_Gap before gap, and Doppler signal Post_Gap behind gap.
It should be noted that, for convenience of description, in the following embodiments, the explanation of each alphabetical noun is generally all occurred once, but this explanation, be all suitable for later description content, be not described in detail at this.
Further, value unit 201, for carrying out value to two sections of Doppler's scanning signals, obtains front complementary series to be inserted and rear complementary series to be checked.
It should be noted that, in the following embodiments, for convenience, when current complementary series to be inserted and rear complementary series to be inserted occur, by it referred to as forward and backward complementary series to be inserted simultaneously.
In present embodiment, value unit 201 is for the repeat its transmission frequency for setting Doppler's scanning;
Obtain the time cycle of every section of Doppler's scanning, and the gap scanning time cycle between the time cycle obtaining two sections of Doppler's scannings; In scanning process, the time cycle of every section of described Doppler's scanning is all identical;
According to repeat its transmission frequency, the time cycle of above-mentioned Doppler's scanning, and the gap scanning time cycle, obtain every section of Doppler's number of scanning line, and the interpolation line number in described gap; In present embodiment: adopt following algorithm to calculate.
Suppose that the time cycle of every section of Doppler's scanning is: T_pw, the described gap scanning time cycle is T_gap, and the repeat its transmission frequency of Doppler's scanning is: prf;
So, every section of Doppler's number of scanning line: N_pw=T_pw*prf;
Interpolation line number in described gap: N_gap=T_gap*prf.
Further, value unit 201 for according to acquisition every section of Doppler's number of scanning line, and the interpolation line number in described gap, get the conjugation of the Doppler's scanning signal behind the Doppler's scanning signal before described gap and described gap, obtain front complementary series to be inserted and rear complementary series to be checked.
In present embodiment: adopt following algorithm to calculate.
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 preorder interpolation sequence; Post _ Sig_ Conj (k) represents follow-up interpolation sequence; N_inter=min (N_pw, N_gap), namely N_inter is N_gap and N_pw smaller value, conj(.) represent and conjugation is got to the value inside bracket.
Further, phase place correcting unit 203 is for carrying out phase place rectification respectively to described front complementary series to be inserted and rear complementary series to be checked.
In present embodiment, phase place correcting unit 203 calculates preorder interpolation sequence and follow-up interpolation sequence needs RA.
In present embodiment: adopt following algorithm to calculate.
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 that front complementary series to be inserted corrects; The angle that after Angle_post, complementary series to be inserted corrects, the phase angle be worth inside bracket is got in operator angle (.) expression.。
Further, phase place correcting unit 203 is multiplied by the unit vector of RA respectively to front sequence to be inserted and rear sequence to be inserted, carry out phase place rectification, obtain new complementary series forward and backward to be inserted.
In present embodiment: adopt following algorithm to calculate.
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) for phase place correct after before complementary series to be inserted; Post _ Sig_ForGap (k) is the complementary series rear to be inserted after phase place rectification; I represents imaginary unit.
Further, frequency correcting unit 205 is for carrying out frequency rectification to described front complementary series to be inserted and rear complementary series to be checked;
In present embodiment, frequency correcting unit 205 is for carrying out frequency rectification to by the described front complementary series to be inserted after phase place rectification module 203 and rear complementary series to be checked.
Before frequency correcting unit 205 pairs of gaps, behind Doppler signal Pre_Gap and gap, the general signal Post_Gap that strangles carries out FFT conversion, and calculate its 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 the mid frequency that pre_fft and post_fft obtains forward and backward complementary series to be inserted corresponding be: pre_f0, post_f0; Bandwidth is: pre_B and post_B;
In present embodiment, frequency correcting unit 205 is holding frequency smooth transition, forward and backward sequence to be filled up needs progressively to change its mid frequency successively in chronological order, consider the time slot between Doppler's scanning signal and gap simultaneously, determine that front complementary series frequency compensation scope to be inserted 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, frequency correcting unit 205 determines forward and backward complementary series compensating frequency sequence to be inserted according to the described mid frequency obtained; And the complementary series forward and backward to be checked after correcting according to described forward and backward complementary series compensating frequency retrieval frequency to be inserted.
In present embodiment: frequency correcting unit 205 adopts following algorithm to carry out calculating and carry out 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 front complementary series compensating frequency sequence to be inserted, and Post_fre_Diff is rear complementary series compensating frequency sequence to be inserted; Pre_Sig_ForGap (k) for frequency correct after before complementary series to be inserted; Post _ Sig_ForGap (k) is the complementary series rear to be inserted after frequency rectification;
Further, bandwidth correcting unit 207 is for carrying out bandwidth rectification to described front complementary series to be inserted and rear complementary series to be checked.
In present embodiment, bandwidth correcting unit 207 is for carrying out frequency rectification to by the described front complementary series to be inserted after phase frequency module 205 and rear complementary series to be checked.
Bandwidth correcting unit 207, in order to ensure that complementary series to be inserted keeps the bandwidth concordance with neighboring Doppler data respectively, carries out corresponding filtering to forward and backward complementary series to be inserted respectively.
In present embodiment, bandwidth correcting unit 207 carries out corresponding filtering to forward and backward complementary series to be inserted respectively by arranging wave filter.
Described wave filter can be IIR or FIR type.
In present embodiment: adopt following algorithm to calculate.
To front sequence filter designing requirement to be inserted be:
Mid frequency is pre_filter_f0=cof1*pre_f0;
Bandwidth requirement is pre_filter_B=cof2*pre_B;
To rear sequence filter designing requirement to be inserted be:
Mid 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 all defaulted as 1.0;
In present embodiment, after bandwidth correcting unit 207 carries out corresponding filtering to forward and backward complementary series to be inserted, acquisition Pre_Sig_ForGap_filtered and post _ Sig_ForGap_filtered, Pre_Sig_ForGap_filtered is the complementary series front to be inserted after bandwidth is corrected; Post _ Sig_ForGap_filtered is the complementary series rear to be inserted after bandwidth is corrected.
It should be noted that, in present embodiment, phase place correcting unit 203, frequency correcting unit 205, bandwidth correct 207 priority execution sequence might not position said sequence, as long as ensure that bandwidth correcting unit performs after frequency correcting unit, namely can complete the present invention, not be described in detail at this.
Further, synthesis unit 209 for by correct through phase place correcting unit 203, frequency correcting unit 205, bandwidth before 207 complementary series to be inserted and after complementary series to be checked be added after windowing respectively, obtain interpolation sequence.
For ensureing progressively transition, synthesis unit 209 couples of Pre_Sig_ForGap_filtered and post _ Sig_ForGap_filtered are added after windowing respectively and obtain interpolation sequence: Sig_For_Gap.
In present embodiment, if N_inter< is N_gap, synthesis unit 209 needs after Pre_Sig_ForGap_filtered, mend 0 makes its length equal N_gap, equally, needs before post _ Sig_ForGap_filtered, mend 0 and makes its length equal N_Gap.
In present embodiment: adopt following algorithm to calculate.
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 windowing, such as: can be oblique angle window win (k)=1-k/N_Gap; Certainly, also can be designed as the shape of other transition gradually, be not described in detail at this.
Further, frequency-domain correction module 300 carries out frequency domain rectification to described interpolation sequence, obtains the time-domain signal inserting described gap.
In present embodiment, noise superpositing unit 301 is for described interpolation its superimposition random phase noise.
Because fir filter interpolator bandwidth may be narrower compared with front and back Doppler signal bandwidth, noise superpositing unit 301 is the concordance of the background white-noise level keeping frequency spectrum, needs the corresponding random phase noise of Signal averaging.
In present embodiment: noise superpositing unit 301 adopts following algorithm to calculate.
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 and is evenly distributed, and also can be designed to other distribution forms such as Gauss, not be described in detail at this.
Further, frequency-domain correction unit 303 is corrected for doing frequency domain to the interpolation sequence after step P1;
In present embodiment: frequency-domain correction unit 303 adopts following algorithm to calculate.
With data grid technology in gap, get length be N_pw data for calculating FFT, the N_pw data Sig_gap_ForFFT wherein obtained represents.
So, record a demerit that to cross coincidence length be 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 couples of Sig_gap_ForFFT do FFT conversion and obtain frequency spectrum gap_FFT, adjust, obtain gap_adjusted_FFT according to pre_FFT and post_FFT to gap_FFT.
Further, frequency-domain correction unit 303 is also for judging that the mould of gap_FFT (k) is whether between pre_FFT (k) and the mould of post_FFT (k), if, by windowing unit 305 to do frequency domain correct after interpolation sequence do windowing and fill up, obtain the time-domain signal inserting described gap; If not, then the amplitude of gap_FFT (k) is adjusted; Wait being adjusted to after between pre_FFT (k) and the mould of post_FFT (k), then by windowing unit 305 to do frequency domain correct after interpolation sequence do windowing and fill up, obtain the time-domain signal inserting described gap.
Further, in present embodiment, frequency-domain correction unit 303 couples of gap_FFT (k) carry out amplitude resize ratio and can adopt nonlinear way, also can adopt simple linear mode.
In present embodiment, further, frequency-domain correction unit 303 set following program segment judge gap_FFT (k) whether between pre_FFT (k) and the mould of post_FFT (k), and it to be 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 length; The modulus of complex number is got in operator fabs () expression.
Further, in present embodiment, windowing unit 305 is filled up for doing windowing to the interpolation sequence after frequency-domain correction unit 303, obtains the time-domain signal inserting described gap.
In present embodiment, in order in the smooth transition of two ends place, windowing unit 305 needs signal windowing process, certainly, the selection of window also can multiformity, such as: shown in Figure 10, so, to gap_adjusted_FFT inverse transformation, obtain, finally for the time-domain signal sig_gap_Result in interpolation space, not being described in detail at this.
Shown in composition graphs 11A, 11B, the spectrogram on the left of screen is by background technology, and U.S. Patent number is the spectrogram that the technical scheme of US5891036 obtains.Spectrogram for being obtained by method shown in the application on the right side of screen.Known by the contrast of Figure 11 A: the high frequency pseudo of the spectrogram obtained by the application's shown device is less.
Known by the contrast of Figure 11 B: the spectrogram obtained by the application's shown device more continuous in the fast-changing place of frequency spectrum.
In sum, ultrasonic system Doppler signal gap interpolation method of the present invention and device, adopt real-time adaptive Doppler interpolation scheme to obtain the interpolation sequence inserting gap; Further, frequency domain rectification is carried out to described interpolation sequence, obtain the time-domain signal inserting described gap; Make the high frequency pseudo of the spectrogram of acquisition less, and more continuous in the fast-changing place of frequency spectrum, improve the interpolation quality of Doppler frequency spectrum.
For convenience of description, various module is divided into describe respectively with function when describing above device.Certainly, the function of each module can be realized in same or multiple software and/or hardware when implementing the application.
As seen through the above description of the embodiments, those skilled in the art can be well understood to the mode that the application can add required general hardware platform by software and realizes.Based on such understanding, the technical scheme of the application can embody with the form of software product the part that prior art contributes in essence in other words, this computer software product can be kept to be preserved in medium, as ROM/RAM, magnetic disc, CD etc., comprising some instructions in order to make a computer equipment (can be personal computer, Information Push Server, or the network equipment etc.) perform the method described in some part of each embodiment of the application or embodiment.
Device embodiments described above is only schematic, the wherein said module illustrated as separating component can or may not be physically separates, parts as module display can be or may not be physical module, namely can be positioned at a place, or also can be distributed on multiple mixed-media network modules mixed-media.Some or all of module wherein can be selected according to the actual needs to realize the object of present embodiment scheme.Those of ordinary skill in the art, when not paying creative work, are namely appreciated that and implement.
The application can be used in numerous general or special purpose computing system environment or configuration.Such as: personal computer, Information Push Server computer, handheld device or portable set, laptop device, multi-processing module system, system, set top box, programmable consumer-elcetronics devices, network PC, minicomputer, mainframe computer, the distributed computing environment comprising above any system or equipment etc. based on micro treatment module.
The application can describe in the general context of computer executable instructions, such as program module.Usually, program module comprises the routine, program, object, assembly, data structure etc. that perform particular task or realize particular abstract data type.Also can put into practice the application in a distributed computing environment, in these distributed computing environment, be executed the task by the remote processing devices be connected by communication network.In a distributed computing environment, program module can be arranged in the local and remote computer preservation medium comprising preservation equipment.
Be to be understood that, although this description is described according to embodiment, but not each embodiment only comprises an independently technical scheme, this narrating mode of description is only for clarity sake, those skilled in the art should by description integrally, technical scheme in each embodiment also through appropriately combined, can form other embodiments that it will be appreciated by those skilled in the art that.
A series of detailed description listed is above only illustrating for feasibility embodiment of the present invention; they are also not used to limit the scope of the invention, all do not depart from the skill of the present invention equivalent implementations done of spirit or change all should be included within protection scope of the present invention.
Claims (12)
1. a ultrasonic system Doppler signal gap interpolation method, is characterized in that, said method comprising the steps of:
M1, obtain two sections of adjacent Doppler's scanning signals, to be respectively before gap Doppler signal after Doppler signal and gap;
Interpolation calculation is carried out in M2, gap between two sections of Doppler's scanning signals, obtains the interpolation sequence inserting described gap;
M3, frequency domain rectification is carried out to described interpolation sequence, obtain the time-domain signal inserting described gap.
2. ultrasonic system Doppler signal gap interpolation method according to claim 1, is characterized in that, described " interpolation calculation is carried out in the gap between two sections of Doppler's scanning signals, obtains the interpolation sequence inserting described gap " specifically comprises:
S1, value is carried out to two sections of Doppler's scanning signals, obtain front complementary series to be inserted and rear complementary series to be checked;
S2, described front complementary series to be inserted and rear complementary series to be checked carried out respectively to phase place rectification, frequency is corrected, bandwidth is corrected;
S3, step S2 is corrected after front complementary series to be inserted and rear complementary series to be checked be added after windowing respectively, obtain interpolation sequence.
3. ultrasonic system Doppler signal gap interpolation method according to claim 2, is characterized in that,
Described frequency is corrected and is specifically comprised: carry out FFT conversion to signal of strangling general behind Doppler signal before gap and gap, calculate its FFT result;
Mid frequency corresponding to forward and backward complementary series to be inserted and bandwidth is obtained according to described FFT result;
Described mid frequency according to obtaining determines front complementary series frequency compensation scope to be inserted, and forward and backward complementary series compensating frequency sequence to be inserted;
Complementary series forward and backward to be checked after correcting according to described forward and backward complementary series compensating frequency retrieval frequency to be inserted.
4. ultrasonic system Doppler signal gap interpolation method according to claim 3, is characterized in that,
Described bandwidth is corrected and is specifically comprised:
By arranging wave filter respectively to the filtering that forward and backward complementary series to be inserted is carried out;
Described wave filter is IIR type or FIR type.
5. the ultrasonic system Doppler signal gap interpolation method according to any one of Claims 1-4, is characterized in that, described " carry out frequency domain rectification to described interpolation sequence, obtain the time-domain signal inserting described gap " specifically comprises:
P1, to described interpolation its superimposition random phase noise;
P2, the interpolation sequence after step P1 done to frequency domain and correct;
P3, windowing is done to the interpolation sequence after step P2 fill up, obtain the time-domain signal inserting described gap.
6. ultrasonic system Doppler signal gap interpolation method according to claim 5, is characterized in that,
Described frequency domain is corrected and is specifically comprised:
According to forward and backward complementary series to be inserted, and insert the described interpolation sequence in described gap, obtain the gap frequency spectrum in described gap;
Judge that the mould of described gap frequency spectrum is whether between the mould of forward and backward complementary series frequency spectrum to be inserted;
If so, step P3 is continued;
If not, after the amplitude of described gap frequency spectrum is adjusted, then step P3 is carried out.
7. a ultrasonic system Doppler signal gap interpolating apparatus, is characterized in that, described device comprises:
Signal acquisition module, for obtaining two sections of adjacent Doppler's scanning signals, to be respectively before gap Doppler signal after Doppler signal and gap;
Interpolating module, carries out interpolation calculation for the gap between two sections of Doppler's scanning signals, obtains the interpolation sequence inserting described gap;
Frequency-domain correction module, for carrying out frequency domain rectification to described interpolation sequence, obtains the time-domain signal inserting described gap.
8. ultrasonic system Doppler signal gap interpolating apparatus according to claim 7, is characterized in that,
Described interpolating module specifically for:
Value is carried out to two sections of Doppler's scanning signals, obtains front complementary series to be inserted and rear complementary series to be checked;
Phase place rectification is carried out respectively to described front complementary series to be inserted and rear complementary series to be checked, frequency is corrected, bandwidth is corrected;
Be added after complementary series to be inserted front after above-mentioned rectification and rear complementary series to be checked respectively windowing, obtain interpolation sequence.
9. ultrasonic system Doppler signal gap interpolating apparatus according to claim 8, is characterized in that,
Described interpolating module carry out frequency correct time specifically for:
FFT conversion is carried out to signal of strangling general behind Doppler signal before gap and gap, calculates its FFT result;
Mid frequency corresponding to forward and backward complementary series to be inserted and bandwidth is obtained according to described FFT result;
Described mid frequency according to obtaining determines front complementary series frequency compensation scope to be inserted, and forward and backward complementary series compensating frequency sequence to be inserted;
Complementary series forward and backward to be checked after correcting according to described forward and backward complementary series compensating frequency retrieval frequency to be inserted.
10. ultrasonic system Doppler signal gap interpolating apparatus according to claim 9, is characterized in that,
Described interpolating module carry out bandwidth correct time specifically for:
By arranging wave filter respectively to the filtering that forward and backward complementary series to be inserted is carried out;
Described wave filter is IIR type or FIR type.
11. ultrasonic system Doppler signal gap interpolating apparatus according to any one of claim 7 to 10, is characterized in that,
Described frequency-domain correction module specifically for:
To described interpolation its superimposition random phase noise;
Do frequency domain to the interpolation sequence after superposition random phase noise to correct;
To do frequency domain correct after interpolation sequence do windowing and fill up, obtain the time-domain signal inserting described gap.
12. ultrasonic system Doppler signal gap interpolating apparatus according to claim 11, is characterized in that,
Described frequency-domain correction module carry out frequency domain correct time also for:
According to forward and backward complementary series to be inserted, and insert the described interpolation sequence in described gap, obtain the gap frequency spectrum in described gap;
Judge that the mould of described gap frequency spectrum is whether between the mould of forward and backward complementary series frequency spectrum to be inserted;
If so, to do frequency domain correct after interpolation sequence do windowing and fill up, obtain the time-domain signal inserting described gap;
If not, after the amplitude of described gap frequency spectrum is adjusted, then to do frequency domain correct after interpolation sequence do windowing and fill up, obtain the time-domain signal inserting described gap.
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 true CN104783836A (en) | 2015-07-22 |
| CN104783836B 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) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107920802A (en) * | 2015-08-27 | 2018-04-17 | 皇家飞利浦有限公司 | Utilize the spectrum doppler processing of adaptively sampled window size |
| CN110141270A (en) * | 2019-06-21 | 2019-08-20 | 青岛海信医疗设备股份有限公司 | Beam synthesizing method and equipment |
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 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107920802A (en) * | 2015-08-27 | 2018-04-17 | 皇家飞利浦有限公司 | Utilize the spectrum doppler processing of adaptively sampled window size |
| CN110141270A (en) * | 2019-06-21 | 2019-08-20 | 青岛海信医疗设备股份有限公司 | Beam synthesizing method and equipment |
| CN110141270B (en) * | 2019-06-21 | 2022-01-28 | 青岛海信医疗设备股份有限公司 | Beam synthesis method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104783836B (en) | 2017-06-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8206302B2 (en) | Method and apparatus for filling doppler signal gaps in ultrasound diagnostic imaging | |
| US10729407B2 (en) | Ultrasonic diagnostic apparatus, image processing apparatus, and image processing method | |
| US7611464B2 (en) | Method for processing Doppler signal gaps | |
| JP5642159B2 (en) | Ultrasound blood flow Doppler sound with pitch shift | |
| US5016641A (en) | Spectral interpolation of ultrasound Doppler signal | |
| JP2005081151A (en) | Tracking clutter filter for spectral and audio doppler | |
| US8992429B2 (en) | Ultrasonic diagnostic apparatus | |
| US20240206851A1 (en) | Filtering method and system for ultrasonic imaging | |
| JP2011254862A (en) | Ultrasonic diagnostic apparatus | |
| JP2001286472A (en) | Ultrasonic diagnostic apparatus | |
| WO2014156269A1 (en) | Ultrasonic imaging device and ultrasonic image display method | |
| JP4094317B2 (en) | Ultrasonic diagnostic equipment | |
| US8162838B2 (en) | Doppler signal spectrum calculating method and apparatus | |
| JP2007202617A (en) | Ultrasonic diagnostic system | |
| CN104783836A (en) | Interval interpolation method and system for Doppler signals of ultrasonic system | |
| JP2010063732A (en) | Ultrasonic diagnostic equipment | |
| JP4698003B2 (en) | Ultrasonic diagnostic equipment | |
| US8652049B2 (en) | Ultrasonic diagnostic apparatus | |
| JP2015119949A (en) | Ultrasonic diagnostic device and control method | |
| JP5838383B2 (en) | Ultrasonic diagnostic apparatus and control program therefor | |
| CN101683276B (en) | Doppler signal gap filling method, device and ultrasound imaging system | |
| JP2002186615A (en) | Ultrasonic daignostic device | |
| JP3693631B2 (en) | Ultrasonic tissue diagnostic equipment | |
| JP4698073B2 (en) | Ultrasonic diagnostic equipment | |
| CN108042155A (en) | Ultrasound echo signal automatic time frequency demodulation method and processing system |
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. |