CN105997147A - Ultrasonic pulse Doppler imaging method and device - Google Patents
Ultrasonic pulse Doppler imaging method and device Download PDFInfo
- Publication number
- CN105997147A CN105997147A CN201610605821.7A CN201610605821A CN105997147A CN 105997147 A CN105997147 A CN 105997147A CN 201610605821 A CN201610605821 A CN 201610605821A CN 105997147 A CN105997147 A CN 105997147A
- Authority
- CN
- China
- Prior art keywords
- matrix
- energy
- signal
- value
- velocity
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Hematology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention discloses an ultrasonic pulse Doppler imaging method and device. The method comprises the steps that an energy matrix is established for echo signals of a received transmitting sequence; a speed sequence for scanning a target is established according to pulse repetition frequency and the central frequency of a probe; energy values corresponding to all speed elements in the speed sequence are obtained according to the energy matrix; an energy spectrum of the speed sequence is generated according to the energy values; an ultrasonic image is generated according to the energy spectrum. The order of a coding compression filter can be effectively reduced, and Doppler imaging can be achieved without a large calculation amount.
Description
Technical field
The present invention relates to ultrasonic pulse technical field of imaging, particularly relate to a kind of ultrasonic pulse Doppler and become
As method and device.
Background technology
Diasonograph (B ultrasonic machine) is with its exclusive characteristic, such as real-time, and dynamic blood flow and doppler imaging,
Become one of indispensable means of auxiliary diagnosis of modern medicine.Wherein, impulse Doppler imaging can be convenient
Measure the concrete flow velocity of blood flow quickly, clinical diagnosis becomes the criterion of some disease.
Traditional coded excitation impulse Doppler imaging technology, pulse compression can be produced bigger by blood flow velocity
Impact, when blood flow is in high flow rate, owing to Doppler effect causes echo deformation.Cause according to original volume
The pulse compression filter of code sequential design can cause the snr gain of pulse compression result to decline and secondary lobe water
Flat rising, the exponent number of compression coding wave filter is the longest simultaneously, needs bigger amount of calculation to realize
Doppler imaging.
Summary of the invention
In view of this, the embodiment of the present invention provides a kind of ultrasonic pulse Doppler imaging method and device, with
Solve technical problem when blood flow in prior art is in high flow rate.
First aspect, embodiments provides a kind of ultrasonic pulse Doppler imaging method, including:
The velocity series of scanning signal is set up according to pulse recurrence frequency and center probe frequency；
The energy value that in acquisition speed sequence, each velocity element is corresponding；
Energy spectrum according to described energy value formation speed sequence；
Compress described energy spectrum, and form ultrasonography according to the energy spectrum after described compression.
Second aspect, the embodiment of the present invention additionally provides a kind of ultrasonic pulse doppler imaging device, including:
Energy matrix sets up module, for the echo-signal receiving transmitting sequence is set up energy matrix；
Velocity series sets up module, for setting up scanning target according to pulse recurrence frequency and center probe frequency
Velocity series；
Energy value acquisition module, for according to each velocity element pair in described energy matrix acquisition speed sequence
The energy value answered；
Energy spectrum generation module, for the energy spectrum according to described energy value formation speed sequence；
Image generation module, for generating ultrasonography according to described energy spectrum.
Ultrasonic pulse Doppler imaging method that the embodiment of the present invention provides and device, by for scanning signal
Velocity series, and calculate the energy value of each velocity element in velocity series, generate energy spectrum, and according to
Energy spectrum forms ultrasonography.Can effectively reduce the exponent number of compression coding wave filter, it is not necessary to bigger meter
Calculation amount can realize doppler imaging.
Accompanying drawing explanation
The detailed description that non-limiting example is made made with reference to the following drawings by reading, the present invention
Other features, objects and advantages will become more apparent upon:
Fig. 1 is the schematic flow sheet of the ultrasonic pulse Doppler imaging method that the embodiment of the present invention one provides；
Fig. 2 is the schematic flow sheet of the ultrasonic pulse Doppler imaging method that the embodiment of the present invention two provides；
Fig. 3 is the schematic flow sheet of the ultrasonic pulse Doppler imaging method that the embodiment of the present invention three provides；
Fig. 4 is the schematic flow sheet of the ultrasonic pulse Doppler imaging method that the embodiment of the present invention four provides；
Fig. 5 is the structural representation of the ultrasonic pulse doppler imaging device that the embodiment of the present invention two provides.
Detailed description of the invention
The present invention is described in further detail with embodiment below in conjunction with the accompanying drawings.It is understood that this
Specific embodiment described by place is used only for explaining the present invention, rather than limitation of the invention.The most also need
Be noted that for the ease of describing, accompanying drawing illustrate only part related to the present invention and not all in
Hold.
Embodiment one
The schematic flow sheet of the ultrasonic pulse Doppler imaging method that Fig. 1 provides for the embodiment of the present invention one,
The present embodiment is applicable to the situation of the echo-signal imaging according to ultrasonic pulse, and the method can be by ultrasonic
Wave impulse doppler imaging device performs, and this device can be realized by software/hardware mode, and can be integrated in phase
In the ultrasonic wave detecting system answered.
See Fig. 1, described ultrasonic pulse Doppler imaging method, including:
S110, sets up energy matrix to the echo-signal receiving transmitting sequence.
The scanning signal launched is encoded by ultrasonic unit, exemplary, according to the size of sampling gate,
Determining exomonental time span, i.e. TimeOfFM=2*SV/c, wherein, SV is the size of sampling gate,
C is the velocity of sound, usual desirable 1540m/s.Mid frequency according to probe calculates each scanning letter in scanning sequence
Number frequency, i.e. fstart=f0 BandWidthOfFM/2；During wherein f0 represents impulse ejection
Frequency of heart, BandWidthOfFM is the modulation bandwidth of transmitted waveform, for arranging value.Frequency according to scanning signal
Rate generates the waveform of each scanning signal in scanning sequence, can calculate basic waveform in the following way:
TxWave_base=cos (2*pi* (fstart*t+BandWidthOfFM/2/TimeOfFM*t.^2)),
Wherein t=(0,1/Tx_fs, 2/Tx_fs ..., TimeOfFM 1/Tx_fs), Tx_fs be system launch
Sample frequency.
Additionally, in order to suppress the concordance of frequency component in the side lobe levels launched and bandwidth, typically can be to volume
Code launches basic waveform windowing, thus obtains final transmitted waveform: can realize windowing in the following way:
TxWave=TxWave_base.*Win, wherein * represents that dot product, Win represent window function.
After launching the scanning signal in scanning sequence by probe, received by the receptor of ultrasonic unit
, owing to there is relative motion, receptor between scanning target and receptor and emitter in the echo of scanning signal
The frequency of sound wave received can produce difference.Sequences of echo signals is set up based on the multiple echo-signals received
Frequency energy matrix.
S120, sets up the velocity series of scanning target according to pulse recurrence frequency and center probe frequency.
When using ultrasonic pulse Doppler to carry out scanning, scanning frequency and the scanning degree of depth have following corresponding pass
Being Prf=C/ (2*D), accordingly, the scanning target in the multi-strip scanning line launched in cycle is
Big speed is: V_max=Prf*c/ (2*f0).Wherein, c is the velocity of sound, and f0 is the center frequency of probe scanning
Rate.Exemplary, the velocity series of scanning target can be blood flowing speed sequence.
Arrange according to maximal rate and system, the velocity series of scanning target can be set up.In the present embodiment,
It can be the quantity of the scan line that ultrasonic system is launched according to pulse recurrence frequency PRF that described system is arranged,
Generally can be 128 or 256 etc..Exemplary, the speed of scanning target can be set up in the following way
Sequence:
V_Dis (i)=i*V_max/K V_max/2, i=(0,1 ... K-1).Wherein V_Dis is speed
Degree series.
S130, according to the energy value that each velocity element in described energy matrix acquisition speed sequence is corresponding.
Owing to the frequency of scanning signal corresponding to each velocity element is different, determine scanning based on this velocity element
The frequency of signal, and according to frequency the energy matrix of the scanning echo-signal set up in advance search corresponding or
Close energy value.
S140, according to the energy spectrum of described energy value formation speed sequence.
The energy value that velocity element each in velocity series is corresponding is added up, obtains velocity series corresponding
Energy value Power_V (i)；And the counting of the energy value corresponding to each velocity element counts and adds up
V_Num (i)=V_Num (i)+1.And according to Power_V (i)=Power_V (i)/V_Num (i)；Raw
Become the energy spectrum of velocity series.
S150, generates ultrasonography according to described energy spectrum.
The energy spectrum generated according to step S140 may determine that the brightness that each velocity element of scanning signal is corresponding shows
Show.Exemplary, set up two-dimensional coordinate system, the transverse axis of described coordinate system is used for showing time, the described longitudinal axis
For display speed, and it is raw to adjust, according to energy spectrum, the brightness that in coordinate system, this point is corresponding, generates ultrasonogram
As checking for user.
The ultrasonic pulse Doppler imaging method of the present embodiment offer and device, by setting up for echo-signal
Energy matrix, and set up the velocity series of scanning signal, and calculate the energy of each velocity element in velocity series
Value, generates energy spectrum, and forms ultrasonography according to energy spectrum.Can effectively reduce compression coding filter
The exponent number of ripple device, it is not necessary to bigger amount of calculation can realize doppler imaging.
Embodiment two
The schematic flow sheet of the ultrasonic pulse Doppler imaging method that Fig. 2 provides for the embodiment of the present invention two.
The present embodiment is optimized based on above-described embodiment, in the present embodiment, and will be to receiving transmitting sequence
Echo-signal set up energy matrix and be specifically optimized for: to the echo-signal received according to described impulse ejection
Mid frequency carry out quadrature demodulation；To the signal sampling after described quadrature demodulation；Generate according to sampled signal
Energy matrix.
Accordingly, the ultrasonic pulse Doppler imaging method that the present embodiment is provided, specifically include:
S210, carries out quadrature demodulation to the echo-signal received according to the mid frequency of described impulse ejection.
The echo-signal of the scanning line for sending, carries out orthogonal according to the mid frequency of ultrasonic equipment probe
Demodulation.It is typically implemented analog to digital conversion and sampling, by Orthogonal phase detection method, positive and negative blood flow signal is taken out,
And by Spectral Analysis Method, positive and negative blood flow is made a distinction.Signal after sampling is multiplied, just with the same frequency of numeric field
The I of base band, Q component can be obtained.
S220, to the signal sampling after described quadrature demodulation.
The signal after quadrature demodulation is intercepted by sampling gate SV.Exemplary, generation of can sampling IQ (1, i),
IQ (2, i), IQ (3, i) ... .IQ (N, i) } sequence, the number of signals in sample sequence has the frequency of sampling gate
Scope, c be ultrasound wave in in-house spread speed, fs is the signal sampling frequency of system.
S230, generates energy matrix according to sampled signal.
First-in first-out can be used to store M I/Q signal sequence, form the I/Q signal square of a M*N
Battle array.IQ component according to sampled signal generates energy matrix.And this matrix is carried out real-time update, when new
Sequence is sent into, and the oldest data will be eliminated.IQ component according to sampled signal is shifted by Fourier transformation
Process to frequency domain.And calculate the calculating energy of the signal I/Q signal after Fourier transformation.Obtain energy matrix.
S240, sets up the velocity series of scanning target according to pulse recurrence frequency and center probe frequency.
S250, according to the energy value that each velocity element in described energy matrix acquisition speed sequence is corresponding.
S260, according to the energy spectrum of described energy value formation speed sequence.
S270, generates ultrasonography according to described energy spectrum.
The present embodiment is by being specifically optimized for the echo-signal receiving transmitting sequence is set up energy matrix:
The echo-signal received is carried out quadrature demodulation according to the mid frequency of described impulse ejection；To described orthogonal
Signal sampling after demodulation；Energy matrix is generated according to sampled signal.Can generate accurately according to I/Q signal
The energy matrix of echo-signal, and energy matrix can be carried out real-time update.
Embodiment three
The schematic flow sheet of the ultrasonic pulse Doppler imaging method that Fig. 3 provides for the embodiment of the present invention three.
The present embodiment is optimized based on above-described embodiment, in the present embodiment, by described according to sampled signal
Generate energy matrix, be specifically optimized for: generate according to sample sequence and update sampled signal matrix；To described
Matrix does two-dimensional Fourier transform, calculates the described sampled signal energy at frequency domain, generates energy matrix；Meter
Calculate the mould square of each element in matrix, it is thus achieved that the energy matrix after conversion.
Accordingly, the ultrasonic pulse Doppler imaging method that the present embodiment is provided, specifically include:
S310, carries out quadrature demodulation to the echo-signal received according to the mid frequency of described impulse ejection.
S320, to the signal sampling after described quadrature demodulation.
S330, generates according to sample sequence and updates sampled signal matrix.
First-in first-out can be used to store M I/Q signal sequence, form the I/Q signal square of a M*N
Battle array.And this matrix is carried out real-time update, when new sequence is sent into, and the oldest data will be eliminated.
S340, does two-dimensional Fourier transform to described matrix, calculates the described sampled signal energy at frequency domain,
Generate energy matrix.
By sampled signal matrix is carried out two-dimensional Fourier transform, sampled signal can be forwarded frequency domain to and carry out
Process.Simultaneously as the frequency of the sampled signal of aliasing has bigger difference, can effectively remove sampling letter
Aliasing in number.
S350, calculates the mould square of each element in transformation matrix, it is thus achieved that the energy matrix after conversion.
Exemplary, can be using the I part of the I/Q signal in transformation matrix as plural real part；Q
Part, as the imaginary part of plural number, calculates the mould square of each I/Q signal, as the energy of this I/Q signal
Value.Calculate the mould square of I/Q signal corresponding to each element successively, obtain the energy matrix after conversion.
S360, sets up the velocity series of scanning target according to pulse recurrence frequency and center probe frequency.
S370, according to the energy value that each velocity element in described energy matrix acquisition speed sequence is corresponding.
S380, according to the energy spectrum of described energy value formation speed sequence.
S390, generates ultrasonography according to described energy spectrum.
The present embodiment, by generating energy matrix by described according to sampled signal, is specifically optimized for: according to sampling
Sequence generates and updates sampled signal matrix；Described matrix is done two-dimensional Fourier transform, calculates described sampling
Signal, at the energy of frequency domain, generates energy matrix；Calculate the mould square of each element in matrix, it is thus achieved that become
Energy matrix after changing.Echo-signal can be forwarded to frequency domain process, effectively remove in echo-signal
Aliasing, improves the accuracy calculating energy matrix.
In a preferred implementation of the present embodiment, before described matrix is done two-dimensional Fourier transform,
Also comprise the steps: described matrix is carried out windowing process；Accordingly, by described, described matrix is done two
Dimension Fourier transformation, is specifically optimized for: the matrix after windowing is done two-dimensional Fourier transform.Window function can
Eliminate or reduce sudden change conversion phenomena, and increase the sensitivity that two dimension Fourier changes, be effectively improved Fourier
The accuracy of conversion.
Additionally, in another preferred implementation of the present embodiment, at the described energy square obtained after converting
After Zhen, increase following steps: the matrix of energy after described conversion is carried out transposition.Due in energy matrix
Every a line can represent the echo of the multiple transmitting signals launching the cycle.And each speed in velocity series
What degree element was corresponding is the transmitting signal of same tranmitting frequency.Accordingly matrix is carried out transposition, can be more
Determine the corresponding relation of the velocity element in velocity series and energy easily.
Embodiment four
The schematic flow sheet of the ultrasonic pulse Doppler imaging method that Fig. 4 provides for the embodiment of the present invention four.
The present embodiment is optimized based on above-described embodiment, in the present embodiment, by every in acquisition speed sequence
The energy value that individual velocity element is corresponding, is specifically optimized for: calculate the frequency values that in transposed matrix, every a line is corresponding；
Judge whether described respective frequencies value is in the bandwidth range of demodulator filter, be in demodulator filter
Time in bandwidth range, calculate described velocity element index value in described row and adjacent integers index value；Root
Interpolation coefficient is calculated according to described adjacent integers index value；According to each in energy ratio coefficient calculations velocity series
The energy value that velocity element is corresponding.
Accordingly, the ultrasonic pulse Doppler imaging method that the present embodiment is provided, specifically include:
S410, sets up energy matrix to the echo-signal receiving transmitting sequence.
S420, sets up the velocity series of scanning target according to pulse recurrence frequency and center probe frequency.
S430, calculates the frequency values that in transposed matrix, every a line is corresponding, it is judged that whether described respective frequencies value is located
In the bandwidth range of demodulator filter.
Owing to the I/Q signal that every a line of transposed matrix is corresponding is same order in different transmitting signal sequences
Launching the echo-signal of signal, the frequency of the most each echo-signal is essentially identical, accordingly, can calculate and turn
Put the frequency values of every a line in matrix.Exemplary, the i-th _ fft row for transposed matrix can use as follows
Mode calculates:
F_IQ (i_fft)=i_fft/N_fft*fs-fs/2；, wherein fs is the signal sampling frequency of system
Rate.
Frequency values according to being calculated in transposed matrix every a line corresponding judges whether it is in demodulator filter
In bandwidth range.Concrete, it is judged that whether frequency meets following condition:
F_IQ (i_fft)>-Biq and f_IQ (i_fft)<Biq, wherein, Biq is demodulator filter bandwidth.
S440, time in the bandwidth range being in demodulator filter, calculates described velocity element in described row
Index value and adjacent integers index value.
If the frequency values that in transposed matrix, every a line is corresponding is in demodulator filter bandwidth range, then according to
The index value that in following manner calculating velocity series, each velocity element is corresponding:
Fd_Indx=mod (2*V_Dis (i)/c* (f_IQ (i_fft)+f0)/prf*M_fft, M_fft)；
Mod () is modulo operation.
Owing to calculated fd_Indx is floating number, its adjacent integers rope can be calculated according to following manner
Draw value:
Fd_preIndx=ceil (fd_Indx)；
Fd_postIndx=floor (fd_Indx)；Wherein ceil () expression rounds up, floor () table
Show and round downwards.
S450, calculates interpolation coefficient according to described adjacent integers index value.
Can be that adjacent integers index value sets up corresponding function, and be calculated at adjacent integers rope by this function
Draw the approximation in value interval, so that the frequency determined is more accurate.Exemplary, distance can be passed through
Chi square function calculates interpolation coefficient, and calculation is as follows:
Dis_a=abs (fd_Indx-fd_preIndx)；
Dis_b=abs (fd_postIndx-fd_Indx)；
C1=Dis_b^2/ (Dis_a^2+Dis_b^2)；
C2=1-c1；
Wherein abs () expression takes absolute value.
In addition, it is possible to calculate interpolation coefficient by absolute value or cube function etc..
S460, calculates, according to interpolation coefficient, the energy value that in velocity series, each velocity element is corresponding.
Energy and count value are added up.Exemplary, can calculate in the following way:
Power_V (i+1)=Power_V (i)+c1*P_2DFFT (i_fft, fd_preIndx)+c2*
P_2DFFT(i_fft,fd_postIndx)
V_Num (i+1)=V_Num (i)+1；
S470, according to the energy spectrum of described energy value formation speed sequence.
After the energy value and counting having calculated all speed is counted, exemplary, can be in the following way
Calculating energy spectrum:
Power_V (i)=Power_V (i)/V_Num (i).
S480, generates ultrasonography according to described energy spectrum.
The present embodiment, by by energy value corresponding for velocity element each in acquisition speed sequence, is specifically optimized for:
Calculate the frequency values that in transposed matrix, every a line is corresponding；Judge whether described respective frequencies value is in demodulation filtering
In the bandwidth range of device, time in the bandwidth range being in demodulator filter, calculate described velocity element in institute
State the index value in row and adjacent integers index value；Interpolation coefficient is calculated according to described adjacent integers index value；
The energy value that in velocity series, each velocity element is corresponding is calculated according to interpolation coefficient.Can calculate exactly
The energy spectrum of velocity series.
Embodiment five
Fig. 5 is the structural representation of the ultrasonic pulse doppler imaging device that the embodiment of the present invention seven provides,
As it is shown in figure 5, described device includes:
Energy matrix sets up module 510, for the echo-signal receiving transmitting sequence is set up energy matrix；
Velocity series sets up module 520, for setting up scanning mesh according to pulse recurrence frequency and center probe frequency
Target velocity series；
Energy value acquisition module 530, for according to each velocity element in described energy matrix acquisition speed sequence
Corresponding energy value；
Energy spectrum generation module 540, for the energy spectrum according to described energy value formation speed sequence；
Image generation module 550, for generating ultrasonography according to described energy spectrum.
The ultrasonic pulse doppler imaging device that the present embodiment provides, by the velocity series for scanning signal,
And calculate the energy value of each velocity element in velocity series, generate energy spectrum, and formed according to energy spectrum super
Sonogram picture.Can effectively reduce the exponent number of compression coding wave filter, it is not necessary to bigger amount of calculation can realize
Doppler imaging.
On the basis of the various embodiments described above, described energy matrix is set up module and is included:
Quadrature demodulation unit, for entering according to the mid frequency of described impulse ejection the echo-signal received
Row quadrature demodulation；
Sampling unit, for the signal sampling after described quadrature demodulation；
Energy matrix signal generating unit, for generating energy matrix according to sampled signal.
On the basis of the various embodiments described above, described energy matrix signal generating unit is used for:
Update subelement, for generating according to sample sequence and updating sampled signal matrix；
Fourier transformation subelement, for described matrix does two-dimensional Fourier transform, generates Fourier transformation
Matrix；
Mould square computation subunit, for calculating the mould square of each element in transformation matrix, it is thus achieved that conversion
After energy matrix.
On the basis of the various embodiments described above, described device also includes:
Transposition unit, for carrying out transposition to the energy matrix after described conversion.
On the basis of the various embodiments described above, described device also includes:
Windowing process unit, for carrying out windowing process to described matrix；
Described Fourier transformation subelement is used for:
Matrix after windowing is done two-dimensional Fourier transform.
On the basis of above-described embodiment, described capacitation value acquisition module is used for:
Calculate the frequency values that in transposed matrix, every a line is corresponding；
Judge whether described respective frequencies value is in the bandwidth range of demodulator filter, be in demodulation filtering
Time in the bandwidth range of device, calculate described velocity element index value in described row and adjacent integers index value；
Interpolation coefficient is calculated according to described adjacent integers index value；
The energy value that in velocity series, each velocity element is corresponding is calculated according to interpolation coefficient.
The ultrasonic pulse doppler imaging device that the present embodiment is provided can be used for performing the present invention and arbitrarily implements
The ultrasonic pulse Doppler imaging method that example provides, possesses corresponding functional module, it is achieved identical is useful
Effect.
Obviously, it will be understood by those skilled in the art that each module or each operation of the above-mentioned present invention can be led to
Cross terminal unit as above to implement.Alternatively, the embodiment of the present invention can be able to perform with computer installation
Program realize, perform such that it is able to be stored in storing in device by processor, described journey
Sequence can be stored in a kind of computer-readable recording medium, and storage medium mentioned above can be read-only depositing
Reservoir, disk or CD etc.；Or they to be fabricated to respectively each integrated circuit modules, or by them
In multiple modules or operation be fabricated to single integrated circuit module and realize.So, the present invention is not restricted to
The combination of any specific hardware and software.
Note, above are only presently preferred embodiments of the present invention and institute's application technology principle.Those skilled in the art
It will be appreciated that the invention is not restricted to specific embodiment described here, can enter for a person skilled in the art
Row various obvious changes, readjust and substitute without departing from protection scope of the present invention.Therefore, though
So by above example, the present invention is described in further detail, but the present invention be not limited only to
Upper embodiment, without departing from the inventive concept, it is also possible to include other Equivalent embodiments more,
And the scope of the present invention is determined by scope of the appended claims.
Claims (12)
1. a ultrasonic pulse Doppler imaging method, it is characterised in that including:
The echo-signal receiving transmitting sequence is set up energy matrix；
The velocity series of scanning target is set up according to pulse recurrence frequency and center probe frequency；
According to the energy value that each velocity element in described energy matrix acquisition speed sequence is corresponding；
Energy spectrum according to described energy value formation speed sequence；
Ultrasonography is generated according to described energy spectrum.
Method the most according to claim 1, it is characterised in that the described echo-signal to receiving is built
Vertical energy matrix, including:
The echo-signal received is carried out quadrature demodulation according to the mid frequency of described impulse ejection；
To the signal sampling after described quadrature demodulation；
Energy matrix is generated according to sampled signal.
Method the most according to claim 2, it is characterised in that described according to sampled signal generation energy
Matrix, including:
Generate according to sample sequence and update sampled signal matrix；
Described matrix is done two-dimensional Fourier transform, generates Fourier transform matrix；
Calculate the mould square of each element in transformation matrix, it is thus achieved that the energy matrix after conversion.
Method the most according to claim 3, it is characterised in that at the described energy square obtained after converting
After Zhen, including:
Energy matrix after described conversion is carried out transposition.
Method the most according to claim 3, it is characterised in that described matrix is being two dimension Fourier
Before conversion, also include:
Described matrix is carried out windowing process；
Described described matrix is done two-dimensional Fourier transform, including:
Matrix after windowing is done two-dimensional Fourier transform.
Method the most according to claim 3, it is characterised in that each speed in described acquisition speed sequence
The energy value that degree element is corresponding, including:
Calculate the frequency values that in transposed matrix, every a line is corresponding；
Judge whether described respective frequencies value is in the bandwidth range of demodulator filter, be in demodulation filtering
Time in the bandwidth range of device, calculate described velocity element index value in described row and adjacent integers index value；
Interpolation coefficient is calculated according to described adjacent integers index value；
The energy value that in velocity series, each velocity element is corresponding is calculated according to interpolation coefficient.
7. a ultrasonic pulse doppler imaging device, it is characterised in that including:
Energy matrix sets up module, for the echo-signal receiving transmitting sequence is set up energy matrix；
Velocity series sets up module, for setting up scanning target according to pulse recurrence frequency and center probe frequency
Velocity series；
Energy value acquisition module, for according to each velocity element pair in described energy matrix acquisition speed sequence
The energy value answered；
Energy spectrum generation module, for the energy spectrum according to described energy value formation speed sequence；
Image generation module, for generating ultrasonography according to described energy spectrum.
Device the most according to claim 7, it is characterised in that described energy matrix is set up module and included:
Quadrature demodulation unit, for entering according to the mid frequency of described impulse ejection the echo-signal received
Row quadrature demodulation；
Sampling unit, for the signal sampling after described quadrature demodulation；
Energy matrix signal generating unit, for generating energy matrix according to sampled signal.
Device the most according to claim 8, it is characterised in that described energy matrix signal generating unit is used for:
Update subelement, for generating according to sample sequence and updating sampled signal matrix；
Fourier transformation subelement, for described matrix does two-dimensional Fourier transform, generates Fourier transformation
Matrix；
Mould square computation subunit, for calculating the mould square of each element in transformation matrix, it is thus achieved that conversion
After energy matrix.
Device the most according to claim 9, it is characterised in that described device also includes:
Transposition unit, for carrying out transposition to the energy matrix after described conversion.
11. devices according to claim 9, it is characterised in that described device also includes:
Windowing process unit, for carrying out windowing process to described matrix；
Described Fourier transformation subelement is used for:
Matrix after windowing is done two-dimensional Fourier transform.
12. devices according to claim 9, it is characterised in that described capacitation value acquisition module is used
In:
Calculate the frequency values that in transposed matrix, every a line is corresponding；
Judge whether described respective frequencies value is in the bandwidth range of demodulator filter, be in demodulation filtering
Time in the bandwidth range of device, calculate described velocity element index value in described row and adjacent integers index value；
Interpolation coefficient is calculated according to described adjacent integers index value；
The energy value that in velocity series, each velocity element is corresponding is calculated according to interpolation coefficient.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610605821.7A CN105997147B (en) | 2016-07-28 | 2016-07-28 | A kind of ultrasonic pulse Doppler imaging method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610605821.7A CN105997147B (en) | 2016-07-28 | 2016-07-28 | A kind of ultrasonic pulse Doppler imaging method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105997147A true CN105997147A (en) | 2016-10-12 |
| CN105997147B CN105997147B (en) | 2019-04-16 |
Family
ID=57114861
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610605821.7A Active CN105997147B (en) | 2016-07-28 | 2016-07-28 | A kind of ultrasonic pulse Doppler imaging method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105997147B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106955098A (en) * | 2017-05-05 | 2017-07-18 | 飞依诺科技(苏州)有限公司 | A kind of vascular flow rate computational methods and device |
| CN108491793A (en) * | 2018-03-22 | 2018-09-04 | 无锡睿勤科技有限公司 | A kind of control method and mobile terminal of fingerprint recognition |
| CN108652666A (en) * | 2017-12-21 | 2018-10-16 | 飞依诺科技(苏州)有限公司 | A kind of generation method and device of doppler blood flow imaging |
| CN110349256A (en) * | 2019-07-16 | 2019-10-18 | 深圳大学 | Reconstructing blood vessel method, apparatus and terminal |
| WO2021103493A1 (en) * | 2019-11-27 | 2021-06-03 | 深圳开立生物医疗科技股份有限公司 | Shear wave-based imaging method, system and apparatus |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5433206A (en) * | 1994-06-29 | 1995-07-18 | Elscint, Ltd. | System and a method for simultaneous, real time ultrasound imaging of biological tissue and measuring of blood flow velocity |
| CN1207889A (en) * | 1997-05-07 | 1999-02-17 | 通用电气公司 | Method and apparatus for optimizing ultrasound color imaging of flow |
| WO2005037060A3 (en) * | 2003-10-03 | 2005-08-25 | Univ Washington | Transcutaneous localization of arterial bleeding by ultrasonic imaging |
| CN1907230A (en) * | 2005-08-02 | 2007-02-07 | 深圳迈瑞生物医疗电子股份有限公司 | Automatic detection method of frequency spectrum Doppler blood flow velocity |
| CN102481141A (en) * | 2009-05-01 | 2012-05-30 | 尤伦·帕提 | Doppler based flow measurements |
| CN105433986A (en) * | 2016-02-01 | 2016-03-30 | 飞依诺科技(苏州)有限公司 | Automatic calibration method and system of ultrasonic equipment |
| CN105997148A (en) * | 2016-05-26 | 2016-10-12 | 飞依诺科技(苏州)有限公司 | Pulse Doppler ultrahigh spectral resolution imaging processing method and processing system |
-
2016
- 2016-07-28 CN CN201610605821.7A patent/CN105997147B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5433206A (en) * | 1994-06-29 | 1995-07-18 | Elscint, Ltd. | System and a method for simultaneous, real time ultrasound imaging of biological tissue and measuring of blood flow velocity |
| CN1207889A (en) * | 1997-05-07 | 1999-02-17 | 通用电气公司 | Method and apparatus for optimizing ultrasound color imaging of flow |
| WO2005037060A3 (en) * | 2003-10-03 | 2005-08-25 | Univ Washington | Transcutaneous localization of arterial bleeding by ultrasonic imaging |
| CN1907230A (en) * | 2005-08-02 | 2007-02-07 | 深圳迈瑞生物医疗电子股份有限公司 | Automatic detection method of frequency spectrum Doppler blood flow velocity |
| CN102481141A (en) * | 2009-05-01 | 2012-05-30 | 尤伦·帕提 | Doppler based flow measurements |
| CN105433986A (en) * | 2016-02-01 | 2016-03-30 | 飞依诺科技(苏州)有限公司 | Automatic calibration method and system of ultrasonic equipment |
| CN105997148A (en) * | 2016-05-26 | 2016-10-12 | 飞依诺科技(苏州)有限公司 | Pulse Doppler ultrahigh spectral resolution imaging processing method and processing system |
Non-Patent Citations (3)
| Title |
|---|
| 李伟锋: "超声多普勒内窥成像***及信号处理的研究", 《万方学位论文数据库》 * |
| 董海艳: "新型血管超声图像处理关键技术研究", 《中国博士学位论文全文数据库 医药卫生科技辑》 * |
| 贾兆航: "便携式超声多普勒血流频谱分析***设计与实现", 《中国优秀硕士学位论文全文数据 医药卫生科技辑》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106955098A (en) * | 2017-05-05 | 2017-07-18 | 飞依诺科技(苏州)有限公司 | A kind of vascular flow rate computational methods and device |
| CN108652666A (en) * | 2017-12-21 | 2018-10-16 | 飞依诺科技(苏州)有限公司 | A kind of generation method and device of doppler blood flow imaging |
| CN108491793A (en) * | 2018-03-22 | 2018-09-04 | 无锡睿勤科技有限公司 | A kind of control method and mobile terminal of fingerprint recognition |
| CN110349256A (en) * | 2019-07-16 | 2019-10-18 | 深圳大学 | Reconstructing blood vessel method, apparatus and terminal |
| CN110349256B (en) * | 2019-07-16 | 2023-05-23 | 深圳大学 | Vascular reconstruction method and device and computer terminal |
| WO2021103493A1 (en) * | 2019-11-27 | 2021-06-03 | 深圳开立生物医疗科技股份有限公司 | Shear wave-based imaging method, system and apparatus |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105997147B (en) | 2019-04-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105997147A (en) | Ultrasonic pulse Doppler imaging method and device | |
| CN102551791B (en) | A kind of ultrasonic imaging method and device | |
| CN106405541B (en) | Full-coherent continuous wave Doppler radar and distance and speed measuring method thereof | |
| JP4928801B2 (en) | Ultrasonic diagnostic equipment | |
| EP3466343B1 (en) | Pulse doppler ultrahigh spectrum resolution imaging processing method and processing system | |
| US6335701B1 (en) | Radar system and coherent integrating method thereof | |
| US20110118606A1 (en) | Adaptively performing clutter filtering in an ultrasound system | |
| US8864671B2 (en) | Methods and systems for color flow imaging | |
| CN101184443B (en) | Velocity measuring method and velocity measuring device using the same | |
| CN103156647A (en) | Ultrasonic diagnosis apparatus, image processing apparatus, and image processing method | |
| US9474510B2 (en) | Ultrasound and system for forming an ultrasound image | |
| CN106991708A (en) | The processing method and processing system of ultrasonic Doppler blood flow imaging | |
| CN104168232A (en) | Method for determining multipath time delay and Doppler frequency shift in underwater acoustic channel | |
| CN115508820A (en) | Target detection method of linear frequency modulation pulse radar | |
| CN105212964B (en) | Based on RF data ultrasonic imaging processing method and system | |
| CN101961251B (en) | Method and device for computing apodization curve in real time in medical ultrasonic diagnosis system | |
| CN108652666B (en) | Doppler blood flow imaging generation method and device | |
| US7803114B2 (en) | Ultrasonic diagnostic apparatus and data processing method therefor | |
| CN106955098B (en) | Blood vessel flow velocity calculation method and device | |
| CN105919624A (en) | Encoding and imaging methods and devices of high-pulse repeated scanning frequency signals | |
| CN107870034B (en) | Underwater acoustic velocity measurement method based on phase difference | |
| CN114545353A (en) | Pulse Doppler signal undersampling and parameter estimation method based on FRI sampling and Chinese remainder theorem | |
| Chi et al. | Improving broadband acoustic Doppler current profiler with orthogonal coprime pulse pairs and robust Chinese remainder theorem | |
| CN114569152A (en) | Blood vessel flow velocity estimation method based on sparse pulse | |
| JP3881078B2 (en) | Frequency estimation method, frequency estimation device, Doppler sonar and tidal meter |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into 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: 215000 5th floor, building a, 4th floor, building C, 27 Xinfa Road, Suzhou Industrial Park, Suzhou City, Jiangsu Province Patentee before: VINNO TECHNOLOGY (SUZHOU) Co.,Ltd. |