CN105997147B - A kind of ultrasonic pulse Doppler imaging method and device - Google Patents
A kind of ultrasonic pulse Doppler imaging method and device Download PDFInfo
- Publication number
- CN105997147B CN105997147B CN201610605821.7A CN201610605821A CN105997147B CN 105997147 B CN105997147 B CN 105997147B CN 201610605821 A CN201610605821 A CN 201610605821A CN 105997147 B CN105997147 B CN 105997147B
- Authority
- CN
- China
- Prior art keywords
- energy
- matrix
- signal
- value
- sequence
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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 a kind of ultrasonic pulse Doppler imaging method and devices, wherein the described method includes: establishing energy matrix to the echo-signal for receiving transmitting sequence；The velocity series of scanning target are established according to pulse recurrence frequency and center probe frequency；According to the corresponding energy value of velocity element each in the energy matrix acquisition speed sequence；According to the energy spectrum of the energy value formation speed sequence；Ultrasonography is generated according to the energy spectrum.The order that coding compression filter can be effectively reduced, can be realized doppler imaging without biggish calculation amount.
Description
Technical field
The present invention relates to ultrasonic pulse technical field of imaging more particularly to a kind of ultrasonic pulse Doppler imaging methods
And device.
Background technique
Diasonograph (B ultrasound machine) is with its exclusive characteristic, such as real-time, dynamic blood flow and doppler imaging, becomes existing
For indispensable one of the means of auxiliary diagnosis of medicine.Wherein, impulse Doppler imaging can conveniently and efficiently measure blood flow
Specific flow velocity becomes the judgment criteria of certain illnesss in clinical diagnosis.
Traditional coded excitation impulse Doppler imaging technology, blood flow velocity, which can compress pulse, generates large effect,
When blood flow is in high flow rate, since Doppler effect causes echo deformation.Lead to the pulse designed according to original coding sequence
Compression filter can cause the snr gain of pulse compression result to decline the rising with side lobe levels, while encode compression filtering
The order of device is often very long, and biggish calculation amount is needed just to be able to achieve 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, it is existing to solve
The technical issues of when thering is in technology blood flow be in high flow rate.
In a first aspect, the embodiment of the invention provides a kind of ultrasonic pulse Doppler imaging methods, comprising:
The velocity series of scanning signal are established according to pulse recurrence frequency and center probe frequency；
The corresponding energy value of each velocity element in acquisition speed sequence；
According to the energy spectrum of the energy value formation speed sequence；
The energy spectrum is compressed, and ultrasonography is formed according to the compressed energy spectrum.
Second aspect, the embodiment of the invention also provides a kind of ultrasonic pulse doppler imaging devices, comprising:
Energy matrix establishes module, for establishing energy matrix to the echo-signal for receiving transmitting sequence；
Velocity series establish module, for establishing the speed of scanning target according to pulse recurrence frequency and center probe frequency
Sequence；
Energy value obtains module, for according to the corresponding energy of velocity element each in the energy matrix acquisition speed sequence
Magnitude；
Energy spectrum generation module, for the energy spectrum according to the energy value formation speed sequence；
Image generation module, for generating ultrasonography according to the energy spectrum.
Ultrasonic pulse Doppler imaging method and device provided in an embodiment of the present invention, pass through the speed for scanning signal
Sequence, and in calculating speed sequence each velocity element energy value, generate energy spectrum, and according to energy spectrum formed ultrasonogram
Picture.The order that coding compression filter can be effectively reduced, can be realized doppler imaging without biggish calculation amount.
Detailed description of the invention
By reading a detailed description of non-restrictive embodiments in the light of the attached drawings below, of the invention other
Feature, objects and advantages will become more apparent upon:
Fig. 1 is the flow diagram for the ultrasonic pulse Doppler imaging method that the embodiment of the present invention one provides；
Fig. 2 is the flow diagram of ultrasonic pulse Doppler imaging method provided by Embodiment 2 of the present invention；
Fig. 3 is the flow diagram for the ultrasonic pulse Doppler imaging method that the embodiment of the present invention three provides；
Fig. 4 is the flow diagram for the ultrasonic pulse Doppler imaging method that the embodiment of the present invention four provides；
Fig. 5 is the structural schematic diagram of ultrasonic pulse doppler imaging device provided by Embodiment 2 of the present invention.
Specific embodiment
The present invention is described in further detail with reference to the accompanying drawings and examples.It is understood that this place is retouched
The specific embodiment stated is used only for explaining the present invention rather than limiting the invention.It also should be noted that in order to just
In description, only some but not all contents related to the present invention are shown in the drawings.
Embodiment one
Fig. 1 is the flow diagram for the ultrasonic pulse Doppler imaging method that the embodiment of the present invention one provides, this implementation
Example is applicable to the case where being imaged according to the echo-signal of ultrasonic pulse, and this method can be by ultrasonic pulse doppler imaging
Device executes, which can be realized by software/hardware mode, and can be integrated in corresponding ultrasonic wave detecting system.
Referring to Fig. 1, the ultrasonic pulse Doppler imaging method, comprising:
S110 establishes energy matrix to the echo-signal for receiving transmitting sequence.
Ultrasonic unit encodes the scanning signal of transmitting, illustratively, according to the size of sampling gate, determines transmitting
The time span of pulse, i.e. TimeOfFM=2*SV/c, wherein SV is the size of sampling gate, and C is the velocity of sound, usually desirable
1540m/s.The frequency of each scanning signal in scanning sequence, i.e. fstart=f0-are calculated according to the centre frequency of probe
BandWidthOfFM/2；Wherein f0 indicates the centre frequency of impulse ejection, and BandWidthOfFM is the modulation band of transmitted waveform
Width is setting value.The waveform of each scanning signal in scanning sequence is generated according to the frequency of scanning signal, it can be using such as lower section
Formula calculates basic waveform: 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 are the sample frequency of system transmitting.
In addition, the consistency in order to inhibit frequency component in the side lobe levels and bandwidth that emit, it generally can be to coded excitation
Basic waveform adding window, to obtain final transmitted waveform: under type such as can be used and realize adding window:
TxWave=TxWave_base.*Win, wherein * indicates dot product, and Win indicates window function.
After through the scanning signal in probe transmitting scanning sequence, scanning letter is received by the receiver of ultrasonic unit
Number echo, due to the frequency of sound wave that, there are relative motion, receiver receives between scanning target and receiver and transmitter
Difference can be generated.The frequency energy matrix of sequences of echo signals is established based on the multiple echo-signals received.
S120 establishes the velocity series of scanning target according to pulse recurrence frequency and center probe frequency.
When carrying out scanning using ultrasonic pulse Doppler, scanning frequency and scanning depth have following corresponding relationship Prf
=C/ (2*D), correspondingly, the maximum speed of the scanning target in the multi-strip scanning line emitted in a cycle are as follows: V_max=
Prf*c/(2*f0).Wherein, c is the velocity of sound, and f0 is the centre frequency of probe scanning.Illustratively, the velocity series of scanning target
It can be blood flowing speed sequence.
It is arranged according to maximum speed and system, can establish the velocity series of scanning target.In the present embodiment, the system
System setting can for ultrasonic system according to the pulse recurrence frequency PRF scan line emitted quantity, generally can be 128 or
256 etc..Illustratively, the velocity series of scanning target can be established in the following way:
V_Dis (i)=i*V_max/K-V_max/2, i=(0,1 ... K-1).Wherein V_Dis is velocity series.
S130, according to the corresponding energy value of velocity element each in the energy matrix acquisition speed sequence.
Since the frequency of the corresponding scanning signal of each velocity element is different, scanning signal is determined based on the velocity element
Frequency, and corresponding or similar energy value is searched in the energy matrix for the scanning echo-signal established in advance according to frequency.
S140, according to the energy spectrum of the energy value formation speed sequence.
It adds up to the corresponding energy value of velocity element each in velocity series, obtains the corresponding energy value of velocity series
Power_V(i)；And to the counting of the corresponding energy value of each velocity element points carry out cumulative V_Num (i)=V_Num (i)+
1.And according to Power_V (i)=Power_V (i)/V_Num (i)；The energy spectrum of formation speed sequence.
S150 generates ultrasonography according to the energy spectrum.
It can determine that the corresponding brightness of each velocity element of scanning signal is shown according to the energy spectrum that step S140 is generated.Show
Example property, two-dimensional coordinate system is established, the horizontal axis of the coordinate system is for showing the time, and the longitudinal axis is for showing speed, and root
It is raw according to the corresponding brightness of the point in energy spectrum adjustment coordinate system, it generates ultrasonography and is checked for user.
Ultrasonic pulse Doppler imaging method and device provided in this embodiment, by establishing energy square for echo-signal
Battle array, and establishes the velocity series of scanning signal, and in calculating speed sequence each velocity element energy value, generate energy spectrum,
And ultrasonography is formed according to energy spectrum.The order of coding compression filter can be effectively reduced, biggish calculation amount is not necessarily to
Doppler imaging can be realized.
Embodiment two
Fig. 2 is the flow diagram of ultrasonic pulse Doppler imaging method provided by Embodiment 2 of the present invention.This implementation
Example is optimized based on above-described embodiment, in the present embodiment, will establish energy to the echo-signal for receiving transmitting sequence
Moment matrix specifically optimizes are as follows: docks the echo-signal received and carries out quadrature demodulation according to the centre frequency of the impulse ejection；It is right
Signal sampling after the quadrature demodulation；Energy matrix is generated according to sampled signal.
Correspondingly, ultrasonic pulse Doppler imaging method provided by the present embodiment, specifically includes:
S210 carries out quadrature demodulation according to the centre frequency of the impulse ejection to the echo-signal received.
For the echo-signal of the scanning line of transmission, quadrature demodulation is carried out according to the centre frequency of ultrasonic equipment probe.
It usually realizes analog to digital conversion and sampling, is taken out positive and negative blood flow signal with Orthogonal phase detection method, and will with Spectral Analysis Method
Positive and negative blood flow distinguishes.Signal after sampling is multiplied with the same frequency of numeric field, so that it may obtain I, the Q component of base band.
S220, to the signal sampling after the quadrature demodulation.
The signal after quadrature demodulation is intercepted by sampling gate SV.Illustratively, can sample generation IQ (1, i), IQ (2,
I), (3, i) IQ ... .IQ (N, i) } sequence, the number of signals in sample sequence has the frequency range of sampling gate, and c is that ultrasonic wave exists
Spread speed in tissue, fs are the signal sampling frequencies of system.
S230 generates energy matrix according to sampled signal.
M I/Q signal sequence can be stored using first-in first-out, form the I/Q signal matrix of a M*N.According to adopting
The IQ component of sample signal generates energy matrix.And real-time update is carried out to the matrix, when new sequence is sent into, oldest data will
It is removed.Frequency domain is transferred to by Fourier transformation according to the IQ component of sampled signal to be handled.And calculate Fourier transformation
Signal I/Q signal afterwards calculates energy.Obtain energy matrix.
S240 establishes the velocity series of scanning target according to pulse recurrence frequency and center probe frequency.
S250, according to the corresponding energy value of velocity element each in the energy matrix acquisition speed sequence.
S260, according to the energy spectrum of the energy value formation speed sequence.
S270 generates ultrasonography according to the energy spectrum.
The present embodiment is specifically optimized by will establish energy matrix to the echo-signal for receiving transmitting sequence are as follows: docking is received
The echo-signal arrived carries out quadrature demodulation according to the centre frequency of the impulse ejection；Signal after the quadrature demodulation is adopted
Sample；Energy matrix is generated according to sampled signal.The energy matrix of echo-signal can be accurately generated according to I/Q signal, and can be right
Energy matrix carries out real-time update.
Embodiment three
Fig. 3 is the flow diagram for the ultrasonic pulse Doppler imaging method that the embodiment of the present invention three provides.This implementation
Example is optimized based on above-described embodiment, in the present embodiment, generates energy matrix according to sampled signal for described, specifically
Optimization are as follows: generated according to sample sequence and update sampled signal matrix；Two-dimensional Fourier transform done to the matrix, described in calculating
Sampled signal generates energy matrix in the energy of frequency domain；The mould square of each element in calculating matrix, obtains transformed energy
Moment matrix.
Correspondingly, ultrasonic pulse Doppler imaging method provided by the present embodiment, specifically includes:
S310 carries out quadrature demodulation according to the centre frequency of the impulse ejection to the echo-signal received.
S320, to the signal sampling after the quadrature demodulation.
S330 is generated according to sample sequence and is updated sampled signal matrix.
M I/Q signal sequence can be stored using first-in first-out, form the I/Q signal matrix of a M*N.And to this
Matrix carries out real-time update, and when new sequence is sent into, oldest data will be removed.
S340 does two-dimensional Fourier transform to the matrix, calculates the sampled signal in the energy of frequency domain, generates energy
Matrix.
By carrying out two-dimensional Fourier transform to sampled signal matrix, sampled signal can be gone into frequency domain and handled.
Simultaneously as the frequency of the sampled signal of aliasing has biggish difference, the aliasing in sampled signal can be effectively removed.
S350 calculates the mould square of each element in transformation matrix, obtains transformed energy matrix.
It illustratively, can be by the part I of the I/Q signal in transformation matrix as plural real part；The conduct of the part Q
The imaginary part of plural number, calculates the mould square of each I/Q signal, the energy value as the I/Q signal.Successively calculate each element pair
The mould square for the I/Q signal answered obtains transformed energy matrix.
S360 establishes the velocity series of scanning target according to pulse recurrence frequency and center probe frequency.
S370, according to the corresponding energy value of velocity element each in the energy matrix acquisition speed sequence.
S380, according to the energy spectrum of the energy value formation speed sequence.
S390 generates ultrasonography according to the energy spectrum.
The present embodiment is specific to optimize are as follows: raw according to sample sequence by generating energy matrix according to sampled signal for described
At and update sampled signal matrix；Two-dimensional Fourier transform is done to the matrix, calculates the sampled signal in the energy of frequency domain,
Generate energy matrix；The mould square of each element in calculating matrix, obtains transformed energy matrix.It can be by echo-signal
It goes to frequency domain to be handled, effectively removes the aliasing in echo-signal, improve the accuracy for calculating energy matrix.
In a preferred embodiment of the present embodiment, before doing two-dimensional Fourier transform to the matrix, also wrap
It includes following steps: windowing process is carried out to the matrix；Correspondingly, doing two-dimensional Fourier transform to the matrix for described, have
Body optimization are as follows: two-dimensional Fourier transform is done to the matrix after adding window.Mutation conversion phenomena can be eliminated or be reduced to window function, and increase
The susceptibility for adding two-dimentional Fourier's variation, effectively improves the accuracy of Fourier transform.
In addition, in another preferred embodiment of the present embodiment, after the transformed energy matrix of acquisition,
Increase following steps: transposition is carried out to the matrix of energy after the transformation.Since a line every in energy matrix can represent one
The echo of multiple transmitting signals of transmit cycle.And it is the same tranmitting frequency that each velocity element in velocity series is corresponding
Transmitting signal.Matrix is subjected to transposition accordingly, it can velocity element in more convenient degree series of constant speed really and energy
Corresponding relationship.
Example IV
Fig. 4 is the flow diagram for the ultrasonic pulse Doppler imaging method that the embodiment of the present invention four provides.This implementation
Example is optimized based on above-described embodiment, in the present embodiment, it is corresponding to will acquire each velocity element in velocity series
Energy value, it is specific to optimize are as follows: to calculate the corresponding frequency values of every a line in transposed matrix；Judge whether the respective frequencies value is in
In the bandwidth range of demodulator filter, when in the bandwidth range in demodulator filter, the velocity element is calculated described
Index value and adjacent integers index value in row；Interpolation coefficient is calculated according to the adjacent integers index value；According to energy ratio
The corresponding energy value of each velocity element in coefficient calculating speed sequence.
Correspondingly, ultrasonic pulse Doppler imaging method provided by the present embodiment, specifically includes:
S410 establishes energy matrix to the echo-signal for receiving transmitting sequence.
S420 establishes the velocity series of scanning target according to pulse recurrence frequency and center probe frequency.
S430 calculates the corresponding frequency values of every a line in transposed matrix, judges whether the respective frequencies value is in demodulation
In the bandwidth range of filter.
Since the corresponding I/Q signal of every a line of transposed matrix is different the transmitting of same sequence in transmitting signal sequence
The echo-signal of signal, therefore the frequency of each echo-signal is essentially identical, accordingly, can calculate every a line in transposed matrix
Frequency values.Illustratively, under type such as can be used for i-th _ fft row of transposed matrix to calculate:
F_IQ (i_fft)=i_fft/N_fft*fs-fs/2；, wherein fs is the signal sampling frequencies of system.
Judge whether it is in demodulator filter bandwidth according to the corresponding frequency values of every a line in transposed matrix are calculated
In range.Specifically, whether determination frequency meets following condition:
F_IQ (i_fft)>-Biq and f_IQ (i_fft)<Biq, wherein Biq is demodulator filter bandwidth.
S440 calculates index of the velocity element in the row when in the bandwidth range in demodulator filter
Value and adjacent integers index value.
If the corresponding frequency values of every a line are in demodulator filter bandwidth range in transposed matrix, according to such as lower section
The corresponding index value of each velocity element in formula calculating speed sequence:
Fd_Indx=mod (2*V_Dis (i)/c* (f_IQ (i_fft)+f0)/prf*M_fft, M_fft)；
Mod () is modulo operation.
Since the fd_Indx being calculated is floating number, its adjacent integers index value can be calculated according to such as under type:
Fd_preIndx=ceil (fd_Indx)；
Fd_postIndx=floor (fd_Indx)；Wherein ceil () expression rounds up, and floor () expression takes downwards
It is whole.
S450 calculates interpolation coefficient according to the adjacent integers index value.
Corresponding function can be established for adjacent integers index value, and is calculated by the function in adjacent integers index value area
Between in approximation so that determine frequency it is more accurate.Illustratively, interpolation system can be calculated by square distance function
Number, calculation are 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, interpolation coefficient can also be calculated by absolute value or cube function etc..
S460, according to the corresponding energy value of velocity element each in interpolation coefficient calculating speed sequence.
Energy and count value are added up.Illustratively, it can be calculated 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 the energy value formation speed sequence.
In the energy value that all speed have been calculated and after counting points, illustratively, energy can be calculated in the following way
Amount spectrum:
Power_V (i)=Power_V (i)/V_Num (i).
S480 generates ultrasonography according to the energy spectrum.
The present embodiment is specific to optimize are as follows: to calculate by will acquire the corresponding energy value of each velocity element in velocity series
The corresponding frequency values of every a line in transposed matrix；Judge whether the respective frequencies value is in the bandwidth range of demodulator filter
It is interior, when in the bandwidth range in demodulator filter, calculate index value of the velocity element in the row and adjacent whole
Number index value；Interpolation coefficient is calculated according to the adjacent integers index value；According to each speed in interpolation coefficient calculating speed sequence
Spend the corresponding energy value of element.The energy spectrum of velocity series can be accurately calculated.
Embodiment five
Fig. 5 is the structural schematic diagram for the ultrasonic pulse doppler imaging device that the embodiment of the present invention seven provides, such as Fig. 5 institute
Show, described device includes:
Energy matrix establishes module 510, for establishing energy matrix to the echo-signal for receiving transmitting sequence；
Velocity series establish module 520, for establishing scanning target according to pulse recurrence frequency and center probe frequency
Velocity series；
Energy value obtains module 530, for corresponding according to velocity element each in the energy matrix acquisition speed sequence
Energy value；
Energy spectrum generation module 540, for the energy spectrum according to the energy value formation speed sequence；
Image generation module 550, for generating ultrasonography according to the energy spectrum.
Ultrasonic pulse doppler imaging device provided in this embodiment, by the velocity series for scanning signal, and is counted
The energy value of each velocity element in velocity series is calculated, generates energy spectrum, and ultrasonography is formed according to energy spectrum.Can have
Effect reduces the order of coding compression filter, and doppler imaging can be realized without biggish calculation amount.
On the basis of the various embodiments described above, the energy matrix establishes module and includes:
Quadrature demodulation unit, it is orthogonal for being carried out to the echo-signal received according to the centre frequency of the impulse ejection
Demodulation；
Sampling unit, for the signal sampling after the quadrature demodulation；
Energy matrix generation unit, for generating energy matrix according to sampled signal.
On the basis of the various embodiments described above, the energy matrix generation unit is used for:
Subelement is updated, for generating according to sample sequence and updating sampled signal matrix；
Fourier transformation subelement generates Fourier transform matrix for doing two-dimensional Fourier transform to the matrix；
Mould square computation subunit obtains transformed energy for calculating the mould square of each element in transformation matrix
Moment matrix.
On the basis of the various embodiments described above, described device further include:
Transposition unit, for carrying out transposition to the transformed energy matrix.
On the basis of the various embodiments described above, described device further include:
Windowing process unit, for carrying out windowing process to the matrix；
The Fourier transformation subelement is used for:
Two-dimensional Fourier transform is done to the matrix after adding window.
On the basis of the above embodiments, the capacitation magnitude obtains module and is used for:
Calculate the corresponding frequency values of every a line in transposed matrix；
Judge whether the respective frequencies value is in the bandwidth range of demodulator filter, in the band for being in demodulator filter
When in wide scope, index value and adjacent integers index value of the velocity element in the row are calculated；
Interpolation coefficient is calculated according to the adjacent integers index value；
According to the corresponding energy value of velocity element each in interpolation coefficient calculating speed sequence.
Ultrasonic pulse doppler imaging device provided by the present embodiment can be used for executing any embodiment of that present invention and mention
The ultrasonic pulse Doppler imaging method of confession, has corresponding functional module, realizes identical beneficial effect.
Obviously, it will be understood by those skilled in the art that each module of the above invention or each operation can be by as above
The terminal device is implemented.Optionally, the embodiment of the present invention can be realized with the program that computer installation can be performed, thus
It can be stored in storage device and be executed by a processor, the program can store computer-readable deposits in a kind of
In storage media, storage medium mentioned above can be read-only memory, disk or CD etc.；Or they are fabricated to
Each integrated circuit modules, or by them multiple modules or operation be fabricated to single integrated circuit module and realize.This
Sample, the present invention is not limited to the combinations of any specific hardware and software.
Note that the above is only a better embodiment of the present invention and the applied technical principle.It will be appreciated by those skilled in the art that
The invention is not limited to the specific embodiments described herein, be able to carry out for a person skilled in the art it is various it is apparent variation,
It readjusts and substitutes without departing from protection scope of the present invention.Therefore, although being carried out by above embodiments to the present invention
It is described in further detail, but the present invention is not limited to the above embodiments only, without departing from the inventive concept, also
It may include more other equivalent embodiments, and the scope of the invention is determined by the scope of the appended claims.
Claims (6)
1. a kind of ultrasonic pulse Doppler imaging method characterized by comprising
Energy matrix is established to the echo-signal for receiving transmitting sequence, includes: to the echo-signal received according to the arteries and veins
The centre frequency of punching transmitting carries out quadrature demodulation, to the signal sampling after the quadrature demodulation, generates energy according to sampled signal
Matrix, wherein it is described according to sampled signal generate energy matrix include: that sampled signal square is generated and updated according to sample sequence
Battle array does two-dimensional Fourier transform to the matrix, generates Fourier transform matrix, calculates the mould of each element in transformation matrix
Square, obtain transformed energy matrix；
Transposition is carried out to the transformed energy matrix；
The velocity series of scanning target are established according to pulse recurrence frequency and center probe frequency；
According to the corresponding energy value of velocity element each in the energy matrix acquisition speed sequence after the transposition；
According to the energy spectrum of the energy value formation speed sequence；
Ultrasonography is generated according to the energy spectrum.
2. the method according to claim 1, wherein being gone back before doing two-dimensional Fourier transform to the matrix
Include:
Windowing process is carried out to the matrix；
It is described that two-dimensional Fourier transform is done to the matrix, comprising:
Two-dimensional Fourier transform is done to the matrix after adding window.
3. the method according to claim 1, wherein each velocity element is corresponding in the acquisition speed sequence
Energy value, comprising:
The corresponding frequency values of every a line in energy matrix after calculating the transposition；
Judge whether the respective frequencies value is in the bandwidth range of demodulator filter, in the bandwidth model for being in demodulator filter
When enclosing interior, index value and adjacent integers index value of the velocity element in corresponding row are calculated；
Interpolation coefficient is calculated according to the adjacent integers index value；
According to the corresponding energy value of velocity element each in interpolation coefficient calculating speed sequence.
4. a kind of ultrasonic pulse doppler imaging device characterized by comprising
Energy matrix establishes module, for establishing energy matrix to the echo-signal for receiving transmitting sequence, specifically includes: orthogonal
Demodulating unit, for carrying out quadrature demodulation according to the centre frequency of the impulse ejection to the echo-signal received, sampling is single
Member, for the signal sampling after the quadrature demodulation, energy matrix generation unit, for generating energy square according to sampled signal
Battle array, wherein energy matrix generation unit includes: update subelement, for being generated according to sample sequence and updating sampled signal square
Battle array, Fourier transformation subelement generate Fourier transform matrix, mould square for doing two-dimensional Fourier transform to the matrix
Computation subunit obtains transformed energy matrix for calculating the mould square of each element in transformation matrix；
Transposition unit, for carrying out transposition to the transformed energy matrix；
Velocity series establish module, for establishing the speed sequence of scanning target according to pulse recurrence frequency and center probe frequency
Column；
Energy value obtains module, for corresponding according to velocity element each in the energy matrix acquisition speed sequence after the transposition
Energy value；
Energy spectrum generation module, for the energy spectrum according to the energy value formation speed sequence；
Image generation module, for generating ultrasonography according to the energy spectrum.
5. device according to claim 4, which is characterized in that described device further include:
Windowing process unit, for carrying out windowing process to the matrix；
The Fourier transformation subelement is used for:
Two-dimensional Fourier transform is done to the matrix after adding window.
6. device according to claim 4, which is characterized in that the energy value obtains module and is used for:
The corresponding frequency values of every a line in energy matrix after calculating the transposition；
Judge whether the respective frequencies value is in the bandwidth range of demodulator filter, in the bandwidth model for being in demodulator filter
When enclosing interior, index value and adjacent integers index value of the velocity element in corresponding row are calculated；
Interpolation coefficient is calculated according to the adjacent integers index value；
According to the corresponding energy value of velocity element each in interpolation coefficient calculating speed sequence.
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 CN105997147A (en) | 2016-10-12 |
| CN105997147B true 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) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106955098B (en) * | 2017-05-05 | 2020-10-16 | 飞依诺科技(苏州)有限公司 | Blood vessel flow velocity calculation method and device |
| CN108042157A (en) * | 2017-12-21 | 2018-05-18 | 飞依诺科技(苏州)有限公司 | A kind of ultrasonic imaging method and device for ultrasound scanning device |
| CN108491793A (en) * | 2018-03-22 | 2018-09-04 | 无锡睿勤科技有限公司 | A kind of control method and mobile terminal of fingerprint recognition |
| CN110349256B (en) * | 2019-07-16 | 2023-05-23 | 深圳大学 | Vascular reconstruction method and device and computer terminal |
| CN110840488B (en) * | 2019-11-27 | 2021-09-28 | 深圳开立生物医疗科技股份有限公司 | Imaging method, system and device based on shear wave |
Family Cites Families (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 |
| US5895358A (en) * | 1997-05-07 | 1999-04-20 | General Electric Company | Method and apparatus for mapping color flow velocity data into display intensities |
| WO2005037060A2 (en) * | 2003-10-03 | 2005-04-28 | University Of Washington | Transcutaneous localization of arterial bleeding by ultrasonic imaging |
| CN100496409C (en) * | 2005-08-02 | 2009-06-10 | 深圳迈瑞生物医疗电子股份有限公司 | Automatic detection method of frequency spectrum Doppler blood flow velocity |
| JP5642770B2 (en) * | 2009-05-01 | 2014-12-17 | パルティ、ヨーラム | Flow measurement based on Doppler |
| CN105433986B (en) * | 2016-02-01 | 2018-07-13 | 飞依诺科技(苏州)有限公司 | Automatic calibration method and system of ultrasonic equipment |
| CN105997148B (en) * | 2016-05-26 | 2019-01-29 | 飞依诺科技(苏州)有限公司 | Pulse Doppler ultrahigh spectral resolution imaging processing method and processing system |
-
2016
- 2016-07-28 CN CN201610605821.7A patent/CN105997147B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105997147A (en) | 2016-10-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105997147B (en) | A kind of ultrasonic pulse Doppler imaging method and device | |
| CN102551791B (en) | A kind of ultrasonic imaging method and device | |
| US6687625B2 (en) | Method and apparatus for feature tracking strain estimation for elastography | |
| US8684934B2 (en) | Adaptively performing clutter filtering in an ultrasound system | |
| CN105997148B (en) | Pulse Doppler ultrahigh spectral resolution imaging processing method and processing system | |
| CN106991708A (en) | The processing method and processing system of ultrasonic Doppler blood flow imaging | |
| JP6179940B2 (en) | Doppler imaging signal transmitter, Doppler imaging signal receiver, Doppler imaging system and method | |
| US9474510B2 (en) | Ultrasound and system for forming an ultrasound image | |
| JP2007292668A (en) | Doppler measuring instrument and tidal current meter | |
| JP6441735B2 (en) | Underwater detector | |
| EP3534177A1 (en) | Method of compressing sonar data | |
| CN114167423A (en) | Radar sea wave parameter measuring method based on depth regression network | |
| CN105212964B (en) | Based on RF data ultrasonic imaging processing method and system | |
| CN102928619B (en) | Signal processing method of broadband beam-control acoustic Doppler flow testing system | |
| CN101576618A (en) | Acoustic positioning measurement method based on wavelet transformation and measurement system thereof | |
| CN117740090A (en) | Ultrasonic clutter elimination method, system and flowmeter | |
| US7803114B2 (en) | Ultrasonic diagnostic apparatus and data processing method therefor | |
| CN108652666B (en) | Doppler blood flow imaging generation method and device | |
| JP4698003B2 (en) | Ultrasonic diagnostic equipment | |
| CN114545353B (en) | Pulse Doppler signal undersampling and parameter estimation method based on FRI sampling and Chinese remainder theorem | |
| CN105919624A (en) | Encoding and imaging methods and devices of high-pulse repeated scanning frequency signals | |
| CN106955098B (en) | Blood vessel flow velocity calculation method and device | |
| JP4827330B2 (en) | Radar wave measuring method and apparatus | |
| CN114569152A (en) | Blood vessel flow velocity estimation method based on sparse pulse | |
| JP5936857B2 (en) | Ultrasonic diagnostic apparatus and image processing program thereof |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| 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. |