US20110178406A1 - Methods and apparatus for hprf doppler ultrasonic imaging - Google Patents
Methods and apparatus for hprf doppler ultrasonic imaging Download PDFInfo
- Publication number
- US20110178406A1 US20110178406A1 US13/007,437 US201113007437A US2011178406A1 US 20110178406 A1 US20110178406 A1 US 20110178406A1 US 201113007437 A US201113007437 A US 201113007437A US 2011178406 A1 US2011178406 A1 US 2011178406A1
- Authority
- US
- United States
- Prior art keywords
- prf
- level
- selectable
- sample volume
- sampling gate
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/5205—Means for monitoring or calibrating
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
- G01S15/8986—Combined Doppler and pulse-echo imaging systems with measures taken for suppressing velocity ambiguities, i.e. anti-aliasing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52019—Details of transmitters
Definitions
- This disclosure relates to pulsed wave spectral Doppler imaging and, in particular, to High Pulse Repeat Frequency (HPRF) Doppler ultrasonic imaging.
- HPRF High Pulse Repeat Frequency
- FIG. 1 is a flowchart of a method for HPRF Doppler ultrasonic imaging
- FIG. 2 is a schematic view of a sampling gate position
- FIG. 3 is a timing diagram showing a period of dead time
- FIG. 4 is a flowchart of method for HPRF Doppler ultrasonic imaging
- FIG. 5 shows Doppler spectrograms before and after automatic optimization of PRF emitted under HPRF mode
- FIG. 6 is a schematic block diagram of an apparatus for HPRF Doppler ultrasonic imaging
- FIG. 7 is a schematic block diagram of an apparatus for HPRF Doppler ultrasonic imaging
- FIG. 8 is a schematic block diagram of a Doppler ultrasonic imaging system including an apparatus for HPRF Doppler ultrasonic imaging.
- PW imaging pulsed wave spectral Doppler imaging
- the ultrasonic front-end emits an ultrasonic pulsed signal at regular intervals into human tissue.
- the Pulse Repeat Interval (“PRI”) will be limited by the SVD. That is, the neighbor PRI should be larger than the time required for the ultrasonic wave to go from the probe to the SVD and back to the probe. If the next ultrasonic wave is sent before the previous ultrasonic wave returns back to the probe, the information received will be mixed with speed information from another target location.
- the pulse repeat frequency (“PRF”) i.e., the reciprocal of PRI, determines the maximum speed range recognized by the current Doppler frequency.
- PRF pulse repeat frequency
- the velocity of blood flow is higher and needs a higher PRF.
- the time for the ultrasonic signal to be sent and received is longer, so larger PRI or lower PRF is needed, and conventional PW imaging fails to meet this requirement.
- high PRF (HPRF) pulsed wave spectral Doppler imaging is used. Higher PRF is utilized so that the echo signal of the nth emission pulse comes back to the probe at the same time with those of the (n ⁇ 1)th, (n ⁇ 2)th, . . . , (n ⁇ m+1)th emission pulse.
- the received echo signal contains the echo information of the emission pulse of the former m times. Therefore, the received spectral information may come from m depth positions.
- the other m ⁇ 1 depth positions are referred to as virtual sample positions, and the original SVD is referred to as the real sample position.
- the user determines which depth position the received velocity information comes from according based on experience. For example, if a certain sample position has no blood flow or tissue movement, the received large velocity won't be from the sample position.
- the depth relation between m ⁇ 1 virtual sample position and the real sample position is as follow:
- c is the velocity of ultrasonic wave in human tissue.
- the PRF is fixed level preset in the system.
- the Doppler spectrogram associated with each PRF is acquired, after which the system calculates one or more corresponding virtual sample position(s) using the above formula so that user can determine where the received spectral information comes from.
- the aforesaid PRF is a preset fixed level and only one value is used in emission and reception in actual use.
- the position of the real sample is determined by the target blood vessel or tissue of the target patient, and the position is variable.
- the PRF may lead to the time of the target signal to the probe being too near to the pulse emission time, such that the received target spectral SNR will be reduced, and the judgment of acquiring target velocity information will be effected. That is, to some real sample positions, which are given arbitrarily, the preset PRFs might be not optimal or even not advisable.
- the present disclosure provides a method and apparatus for improving the performance of HPRF Doppler ultrasonic imaging by means of transmitting with the best PRF, wherein the best PRF is calculated by a PRF level selected by the user and a sampling gate position.
- a method for improving HPRF performance in Doppler ultrasonic imaging includes collecting parameters selected by a user through an interface, the parameters including a selected PRF level, a real sample volume depth, and a sampling gate size; setting a selectable PRF range according to the standard PRF value corresponding to the selected PRF level; evaluating each PRF successively selected from the selectable PRF range with a preset step; and selecting the PRF with best evaluation as the best PRF of the selected PRF level, which is then utilized in transmitting ultrasonic pulses.
- an apparatus for improving HPRF performance in Doppler ultrasonic imaging includes: a collecting unit configured for collecting parameters selected by a user through an interface, the parameters including the selected PRF level, real sample volume depth and sampling gate size; a setting unit configured for setting an selectable PRF range according to the standard PRF value corresponding to the selected PRF level; and an evaluating unit configured for evaluating each PRF successively selected from the selectable PRF range with a preset step; wherein the PRF with best evaluation is taken as the best PRF of the selected PRF level and utilized in transmitting ultrasonic pulses.
- the ultrasonic pulse with the best PRF is emitted in real time, wherein the best PRF is calculated with the PRF level and the SVD chosen by user, and the SNR of the Doppler spectrogram is enhanced under HPRF mode, thus the imaging effect of the Doppler spectrum is improved.
- a method for HPRF Doppler ultrasonic imaging includes a collecting step 100 , a setting step 102 , and an evaluating step 106 .
- the method may optionally include a judging step of sampling gate size 104 and/or a reinforcing judgment step 108 .
- step 100 parameters selected by user through a user interface are collected, wherein the parameters include information relating to the sampling gate position and selected PRF level.
- the information relating to the sampling gate position includes a sample volume depth (SVD), i.e., GateDepth in FIG. 2 , and a sampling gate size, i.e., GateSize in FIG. 2 .
- the sampling gate herein is the real sampling gate and the corresponding SVD is the real SVD.
- the standard PRF value corresponding to the selected PRF level selected through the user interface is marked as usePRF.
- the selectable PRF range can be directly set by user through a user interface, or automatically generated by an ultrasonic system.
- the selectable PRF range can be set by multiplying the standard PRF value corresponding to the selected PRF level by a preset upper limit parameter and a preset lower limit parameter respectively.
- a preset upper limit parameter and a preset lower limit parameter respectively.
- the lower limit and upper limit parameters mainly affect the optimized adjustment range of each PRF level. They can be determined based on user's demand and the effect of image optimization.
- the selectable PRF range can be set by subtracting and adding the standard PRF value corresponding to the selected PRF level by preset adjustment parameters respectively. For example, in one embodiment,
- Prf_Low usePRF ⁇ a
- Prf_High usePRF+ b
- the parameters mainly affect the optimized adjustment range of each PRF level. They can be determined based on user's demand and the effect of image optimization.
- the parameters ⁇ , ⁇ , a and b may be directly set by user through a user interface, or automatically set by default by an ultrasonic system.
- a plurality of methods for setting a selectable PRF range is described in the aforesaid embodiments. It will be understood by those skilled in the art that other methods can be utilized based on actual conditions, user demands, and the effect of image optimization.
- the methods for setting the PRF range mentioned in the above embodiments are exemplary and should not be construed as limiting.
- FIG. 3 shows a time series and the positions of dead time of ultrasonic pulse transmission and reception, wherein Rate means the beginning of each transmission and reception cycle, and nTxgate is a transmission time series.
- the effective time T 3 of transmission in each cycle is a dead time.
- the sampling line of HPRF is determined, the emission order of each transducer in the probe and the relative time delay can be determined.
- the transmitting effective time T 3 in the time series can be determined.
- the time length of dead time may be the effective time of transmission added by an additional moment, that is:
- the AdditionalDeadTime is the additional dead time which is a fixed value preset in the system.
- the fixed value is related to the method for designing the time series in the system and the performance of hardware used for emission.
- the system calculates the times of the target signal returns to the probe at the beginning and end sampling gate positions, for example:
- GateStart 2 ⁇ ( GateDepth - GateSize / 2 ) c + DeadTime
- GateEnd 2 ⁇ ( GateDepth + GateSize / 2 ) c + DeadTime
- c is the velocity of sound in the tissue.
- An upper limit parameter of the number of dead time set in the system is MaxGateNumber.
- the parameter which corresponds to the maximum number of sampling gates under HPRF mode, can be set based on the user's demand or the desired image effect.
- GateNumber and usePRF may be adjusted in a forcible manner as follow:
- ⁇ is a fixed parameter set in the system.
- the parameter is usually set as a small value and can be adjusted with image effect.
- the time of each depth of the sampling gate is GateStart ⁇ GateEnd corresponding from emission to reception.
- the GateStart cannot be fallen in the previous adjacent dead time and the GateEnd can not be in the next adjacent dead time.
- GatePrf_Low fix ⁇ ( GateNumber )
- GatePrf_High fix ⁇ ( GateNumber ) + 1 GateEnd
- fix means rounding a number to the nearest integer.
- the real emission PRF range should be in the selectable PRF range restricted by the sampling gate, such that:
- Temp_Prf_Low max(GatePrf_Low,Prf_Low)
- Temp_Prf_High min(GatePrf_High,Prf_High)
- the system evaluates (step 106 ) each PRF successively selected from the selectable PRF range with a preset step.
- the step Prf_Step may be set flexibly according to the user's demand, the effect of image optimization, or the time and resource spent on evaluation.
- the system then takes the PRF evaluated currently as Prf_temp and calculates the corresponding number of dead time with, for example, the following formula:
- the scoring may contain three scoring items: the scoring distance to the previous adjacent dead time, the scoring distance to the next adjacent dead time, and the scoring similarity with the standard PRF value corresponding to current PRF level. The three steps will be further specifically described herein.
- the time distance from a beginning time to the previous dead time may be as follow, where the beginning time is the one at which the signal at the beginning position of sampling gate returns to the probe
- the score may be function of the time distance. The larger the time distance, the better the score.
- the time distance from an end time to the next dead time may be as follows, wherein the end time is the one at which the signal at the end position of sampling gate returns to the probe:
- DisEnd fix ⁇ ( GateNumber ) + 1 Prf_temp - GateEnd
- the score may be function of the time distance. The larger the time distance, the better the score.
- the final score may be function of one of the three items (including the case where the score is one item itself, which is also called as the function of the item in the disclosure), or any two of the three items (i.e., the product of any two items), or the function of the three items (i.e., the product of the three scoring items); wherein the final score function can be set flexibly, for example, it can be each item itself, or can be other functions.
- the PRF with best total score is taken as the best PRF of the current PRF level.
- the size of the sampling gate after adjustment is output at the same time.
- the scoring of PRF is then finished.
- the definition of “the best” or “the worst” may be defined flexibly and is not limited in the disclosure.
- “the best” can be defined as the higher score and “the worst” as the lower score; alternatively, “the best” can be defined as the lower score and “the worst” as the higher score.
- the score items may be other which can measure the performance of current Prf_temp as understood by those skilled in the art.
- the scoring methods are not limited in the three score items. They may be any score items which can be used to evaluate the performance of current Prf_temp.
- the number of virtual sampling gates and their corresponding SVD can be calculated with the current real SVD as follows:
- the current PRF level is not suitable.
- the next PRF level is should be scored to search for the best PRF corresponding to next PRF level (step 108 ).
- the threshold DepthThre is preset and may be flexibly adjusted based on actual image effect.
- the reinforcing judgment step 108 and the evaluating step of PRF 106 may be interchanged. That is, the reinforcing judgment step may be performed by utilizing the standard PRF value of the current level so as to determine whether to adjust the level. After searching the suitable level, the method may proceed to the evaluating step to search for the best PRF value corresponding to the level, and set the best PRF to be the PRF utilized in transmitting ultrasonic pulses in the system. In this case, the selected PRF is adopted in the aforesaid formulas about the depth of each virtual sampling gate.
- each of the aforesaid embodiments there are a plurality of preset parameters. These parameters may be set flexibly according to actual situations. They may be directly set by a user through a user interface, or may be set by default by an ultrasonic system.
- FIG. 5 shows the image before and after optimization, wherein FIG. 5 a is the Doppler spectrum before optimization and FIG. 5 b is the optimized Doppler spectrum. From FIG. 5 , after HPRF is optimized, the Doppler signals are markedly enhanced, the SNR is improved, and the contour of the velocity information is more clearly shown.
- FIG. 6 illustrates an apparatus 600 for HPRF Doppler ultrasonic imaging according to one embodiment.
- the apparatus 600 may include a collecting unit 620 , a setting unit 630 , and an evaluating unit 650 .
- the apparatus 600 optionally further includes a judging unit of sampling gate size 640 and/or a reinforcing judgment unit 660 .
- the collecting unit 620 is configured to implement step 100 of FIG. 1 ;
- the setting unit 630 is configured to implement step 102 ;
- the judging unit of sampling gate size 640 is configured to implement step 104 ;
- the evaluating unit 650 is configured to implement step 104 ;
- the reinforcing judgment unit 660 is configured to implement step 108 .
- the apparatus 600 for HPRF Doppler ultrasonic imaging optionally further includes an interface 680 and an output unit 690 .
- the interface 680 may be configured to collect the information selected by the user, i.e., PRF level, SVD, and the sampling gate size.
- the output unit 690 may be configured to output the best PRF value.
- the apparatus 600 for HPRF Doppler ultrasonic imaging comprises a processing unit 613 , i.e., DSP or CPU.
- the processing unit 613 may be a single unit or a plurality of units for implementing the above-mentioned steps.
- the apparatus optionally further includes an interface 680 and an output unit 690 .
- the apparatus 600 further comprises at least one non-transitory computer-readable medium 610 which may be in form of non-volatile memory, i.e., EEPROM, flash or disk drive.
- the computer-readable medium 610 may include a computer program 611 including program codes, which, when executed by the processing unit 613 , implements the steps shown in FIG. 1 and/or FIG. 4 .
- the computer program 611 in the apparatus 600 may include a collecting module 611 a for implementing step 100 , a setting module 611 b for implementing step 102 , a judging module of sampling gate size 611 c for implementing step 104 , an evaluating module 611 d for implementing step 106 , and a reinforcing judgment module 611 e for implementing step 108 .
- the modules 611 a - 611 e when running the modules 611 a - 611 e on the processing unit 613 , the modules respectively correspond to the units 620 , 630 , 640 , 650 , 660 shown in FIG. 6 .
- the apparatus 600 for HPRF Doppler ultrasonic imaging can be realized in a Doppler ultrasonic imaging system with software, hardware, firmware, or their combination, as shown in FIG. 8 .
Abstract
Description
- This application claims the benefit of Chinese Patent Application No. 201010044456.X, filed on Jan. 21, 2010, which is incorporated herein by reference in its entirety.
- This disclosure relates to pulsed wave spectral Doppler imaging and, in particular, to High Pulse Repeat Frequency (HPRF) Doppler ultrasonic imaging.
-
FIG. 1 is a flowchart of a method for HPRF Doppler ultrasonic imaging; -
FIG. 2 is a schematic view of a sampling gate position; -
FIG. 3 is a timing diagram showing a period of dead time; -
FIG. 4 is a flowchart of method for HPRF Doppler ultrasonic imaging; -
FIG. 5 shows Doppler spectrograms before and after automatic optimization of PRF emitted under HPRF mode; -
FIG. 6 is a schematic block diagram of an apparatus for HPRF Doppler ultrasonic imaging; -
FIG. 7 is a schematic block diagram of an apparatus for HPRF Doppler ultrasonic imaging; -
FIG. 8 is a schematic block diagram of a Doppler ultrasonic imaging system including an apparatus for HPRF Doppler ultrasonic imaging. - During the process of pulsed wave spectral Doppler imaging (hereinafter “PW imaging”) in a medical ultrasonic imaging system, to detect the Doppler frequency shift, the ultrasonic front-end emits an ultrasonic pulsed signal at regular intervals into human tissue. In conventional PW imaging, to ensure that speed information is accurately obtained at the target location, i.e., Sample Volume Depth (“SVD”), the Pulse Repeat Interval (“PRI”) will be limited by the SVD. That is, the neighbor PRI should be larger than the time required for the ultrasonic wave to go from the probe to the SVD and back to the probe. If the next ultrasonic wave is sent before the previous ultrasonic wave returns back to the probe, the information received will be mixed with speed information from another target location. The pulse repeat frequency (“PRF”), i.e., the reciprocal of PRI, determines the maximum speed range recognized by the current Doppler frequency. In some cases, in a large arterial vessel, for example, the velocity of blood flow is higher and needs a higher PRF. However, if the location of the vessel is deeper, the time for the ultrasonic signal to be sent and received is longer, so larger PRI or lower PRF is needed, and conventional PW imaging fails to meet this requirement.
- To solve these problems, high PRF (HPRF) pulsed wave spectral Doppler imaging is used. Higher PRF is utilized so that the echo signal of the nth emission pulse comes back to the probe at the same time with those of the (n−1)th, (n−2)th, . . . , (n−m+1)th emission pulse. Thus, the received echo signal contains the echo information of the emission pulse of the former m times. Therefore, the received spectral information may come from m depth positions. Except for the real target position SVD, the other m−1 depth positions are referred to as virtual sample positions, and the original SVD is referred to as the real sample position. In the actual use, the user determines which depth position the received velocity information comes from according based on experience. For example, if a certain sample position has no blood flow or tissue movement, the received large velocity won't be from the sample position. The depth relation between m−1 virtual sample position and the real sample position is as follow:
-
- where c is the velocity of ultrasonic wave in human tissue.
- According to the HPRF technology, the PRF is fixed level preset in the system. The Doppler spectrogram associated with each PRF is acquired, after which the system calculates one or more corresponding virtual sample position(s) using the above formula so that user can determine where the received spectral information comes from.
- The shortcoming of existing HPRF technology is that the aforesaid PRF is a preset fixed level and only one value is used in emission and reception in actual use. However, the position of the real sample is determined by the target blood vessel or tissue of the target patient, and the position is variable. In some real sample positions, the PRF may lead to the time of the target signal to the probe being too near to the pulse emission time, such that the received target spectral SNR will be reduced, and the judgment of acquiring target velocity information will be effected. That is, to some real sample positions, which are given arbitrarily, the preset PRFs might be not optimal or even not advisable.
- The present disclosure provides a method and apparatus for improving the performance of HPRF Doppler ultrasonic imaging by means of transmitting with the best PRF, wherein the best PRF is calculated by a PRF level selected by the user and a sampling gate position.
- According to the one aspect of the disclosure, a method for improving HPRF performance in Doppler ultrasonic imaging includes collecting parameters selected by a user through an interface, the parameters including a selected PRF level, a real sample volume depth, and a sampling gate size; setting a selectable PRF range according to the standard PRF value corresponding to the selected PRF level; evaluating each PRF successively selected from the selectable PRF range with a preset step; and selecting the PRF with best evaluation as the best PRF of the selected PRF level, which is then utilized in transmitting ultrasonic pulses.
- According to another aspect, an apparatus for improving HPRF performance in Doppler ultrasonic imaging includes: a collecting unit configured for collecting parameters selected by a user through an interface, the parameters including the selected PRF level, real sample volume depth and sampling gate size; a setting unit configured for setting an selectable PRF range according to the standard PRF value corresponding to the selected PRF level; and an evaluating unit configured for evaluating each PRF successively selected from the selectable PRF range with a preset step; wherein the PRF with best evaluation is taken as the best PRF of the selected PRF level and utilized in transmitting ultrasonic pulses.
- In one embodiment, the ultrasonic pulse with the best PRF is emitted in real time, wherein the best PRF is calculated with the PRF level and the SVD chosen by user, and the SNR of the Doppler spectrogram is enhanced under HPRF mode, thus the imaging effect of the Doppler spectrum is improved.
- Referring to
FIG. 1 , a method for HPRF Doppler ultrasonic imaging includes a collectingstep 100, asetting step 102, and an evaluatingstep 106. In other embodiments, the method may optionally include a judging step ofsampling gate size 104 and/or a reinforcingjudgment step 108. - In
step 100, parameters selected by user through a user interface are collected, wherein the parameters include information relating to the sampling gate position and selected PRF level. The information relating to the sampling gate position includes a sample volume depth (SVD), i.e., GateDepth inFIG. 2 , and a sampling gate size, i.e., GateSize inFIG. 2 . The sampling gate herein is the real sampling gate and the corresponding SVD is the real SVD. Additionally, the standard PRF value corresponding to the selected PRF level selected through the user interface is marked as usePRF. - In
step 102, a selectable PRF range Prf_Low˜Prf_High for optimization is set, wherein Prf_Low<=usePRF and Prf_High>=usePRF. The selectable PRF range can be directly set by user through a user interface, or automatically generated by an ultrasonic system. - For example, in one embodiment, the selectable PRF range can be set by multiplying the standard PRF value corresponding to the selected PRF level by a preset upper limit parameter and a preset lower limit parameter respectively. For example, in one embodiment,
-
Prf_Low=α*usePRF -
Prf_High=β*usePRF - where α and β are the lower limit and upper limit parameters respectively preset by the system, and α<=1, β>=1. The lower limit and upper limit parameters mainly affect the optimized adjustment range of each PRF level. They can be determined based on user's demand and the effect of image optimization.
- In one embodiment, the selectable PRF range can be set by subtracting and adding the standard PRF value corresponding to the selected PRF level by preset adjustment parameters respectively. For example, in one embodiment,
-
Prf_Low=usePRF−a -
Prf_High=usePRF+b - where a and b are the PRF adjustment parameters respectively preset by the system, and a>=0, b>=0. The parameters mainly affect the optimized adjustment range of each PRF level. They can be determined based on user's demand and the effect of image optimization.
- In the aforesaid embodiments, according to user's demand, the parameters α, β, a and b may be directly set by user through a user interface, or automatically set by default by an ultrasonic system.
- A plurality of methods for setting a selectable PRF range is described in the aforesaid embodiments. It will be understood by those skilled in the art that other methods can be utilized based on actual conditions, user demands, and the effect of image optimization. The methods for setting the PRF range mentioned in the above embodiments are exemplary and should not be construed as limiting.
-
FIG. 3 shows a time series and the positions of dead time of ultrasonic pulse transmission and reception, wherein Rate means the beginning of each transmission and reception cycle, and nTxgate is a transmission time series. The effective time T3 of transmission in each cycle is a dead time. The time interval between two adjacent Rate is a repeat time interval PRI=1/PRF of transmission corresponding to the current transmitting PRF. Under the HPRF mode, if the time of the target signal at the current SVD returns to the probe is just falling in the dead time or very close to the dead time, the receiving of the target signal may meet obstacles or interferences, which reduces the SNR of Doppler spectral information. - If the sampling line of HPRF is determined, the emission order of each transducer in the probe and the relative time delay can be determined. Thus the transmitting effective time T3 in the time series can be determined. In actual use, the time length of dead time may be the effective time of transmission added by an additional moment, that is:
-
DeadTime=T 3+AdditionalDeadTime - where the AdditionalDeadTime is the additional dead time which is a fixed value preset in the system. The fixed value is related to the method for designing the time series in the system and the performance of hardware used for emission.
- The system calculates the times of the target signal returns to the probe at the beginning and end sampling gate positions, for example:
-
- wherein c is the velocity of sound in the tissue.
- The number of dead time then is:
-
GateNumber=GateEnd*Prf_Low - An upper limit parameter of the number of dead time set in the system is MaxGateNumber. The parameter, which corresponds to the maximum number of sampling gates under HPRF mode, can be set based on the user's demand or the desired image effect. When the calculated GateNumber is equal to or larger than the upper limit parameter, GateNumber and usePRF may be adjusted in a forcible manner as follow:
-
- where ε is a fixed parameter set in the system. The parameter is usually set as a small value and can be adjusted with image effect.
- Due to the size of sampling gate, the time of each depth of the sampling gate is GateStart˜GateEnd corresponding from emission to reception. To eliminate the interference from dead time, the GateStart cannot be fallen in the previous adjacent dead time and the GateEnd can not be in the next adjacent dead time. Accordingly, in one embodiment, the selectable PRF range is restricted by the current sampling gate may be:
-
- wherein fix means rounding a number to the nearest integer.
- The real emission PRF range should be in the selectable PRF range restricted by the sampling gate, such that:
-
Temp_Prf_Low=max(GatePrf_Low,Prf_Low) -
Temp_Prf_High=min(GatePrf_High,Prf_High) - If Temp_Prf_Low<=Temp_Prf_High, an effective PRF can be searched in the adjustment range by the system; otherwise, the system fails to search the effective PRF and the sampling gate size must be reduced (step 104).
- Additionally, after reducing the sampling gate size, a new judgment is needed based on the reduced sampling gate until a suitable size of the sampling gate is searched or the size is adjusted to the minimum sampling gate size restricted by the system.
- The system evaluates (step 106) each PRF successively selected from the selectable PRF range with a preset step. The step Prf_Step may be set flexibly according to the user's demand, the effect of image optimization, or the time and resource spent on evaluation.
- The system then takes the PRF evaluated currently as Prf_temp and calculates the corresponding number of dead time with, for example, the following formula:
-
GateNumber=GateEnd*Prf_temp - If GateNumber>=MaxGateNumber, the Prf_temp fails to meet the requirement and the total score of Prf_temp is set with a worst score (i.e., set to zero, that is Score=0); if GateNumber<MaxGateNumber, scoring the Prf_temp (step 106). The scoring may contain three scoring items: the scoring distance to the previous adjacent dead time, the scoring distance to the next adjacent dead time, and the scoring similarity with the standard PRF value corresponding to current PRF level. The three steps will be further specifically described herein.
- 1. Scoring Distance to the Previous Adjacent Dead Time
- The time distance from a beginning time to the previous dead time may be as follow, where the beginning time is the one at which the signal at the beginning position of sampling gate returns to the probe
-
- The score may be function of the time distance. The larger the time distance, the better the score.
- 2. Scoring Distance to the Next Adjacent Dead Time
- The time distance from an end time to the next dead time may be as follows, wherein the end time is the one at which the signal at the end position of sampling gate returns to the probe:
-
- The score may be function of the time distance. The larger the time distance, the better the score.
- 3. The Scoring Similarity with the Standard PRF Value Corresponding to Current PRF Level
- Comparing the difference between the PRF evaluated currently and the standard PRF value corresponding to the current PRF level, the similarity is the function of difference. The smaller the difference, that is, Prf_temp is closer to use PRF, the better the score.
- 4. Total Score
- After obtaining three scoring items, the final score may be function of one of the three items (including the case where the score is one item itself, which is also called as the function of the item in the disclosure), or any two of the three items (i.e., the product of any two items), or the function of the three items (i.e., the product of the three scoring items); wherein the final score function can be set flexibly, for example, it can be each item itself, or can be other functions.
- Among the aforesaid selectable PRF range, the PRF with best total score is taken as the best PRF of the current PRF level. The size of the sampling gate after adjustment is output at the same time. The scoring of PRF is then finished.
- It will be understood by those skilled in the art that the definition of “the best” or “the worst” may be defined flexibly and is not limited in the disclosure. For example, “the best” can be defined as the higher score and “the worst” as the lower score; alternatively, “the best” can be defined as the lower score and “the worst” as the higher score.
- Additionally, besides the above-mentioned three score items, the score items may be other which can measure the performance of current Prf_temp as understood by those skilled in the art. The scoring methods are not limited in the three score items. They may be any score items which can be used to evaluate the performance of current Prf_temp.
- After obtaining the best PRF, the number of virtual sampling gates and their corresponding SVD can be calculated with the current real SVD as follows:
-
- If the lightest SVD of the virtual sampling gate is smaller than a threshold DepthThre, the current PRF level is not suitable. The next PRF level is should be scored to search for the best PRF corresponding to next PRF level (step 108). The threshold DepthThre is preset and may be flexibly adjusted based on actual image effect.
- In another embodiment, as shown in
FIG. 4 , the reinforcingjudgment step 108 and the evaluating step ofPRF 106 may be interchanged. That is, the reinforcing judgment step may be performed by utilizing the standard PRF value of the current level so as to determine whether to adjust the level. After searching the suitable level, the method may proceed to the evaluating step to search for the best PRF value corresponding to the level, and set the best PRF to be the PRF utilized in transmitting ultrasonic pulses in the system. In this case, the selected PRF is adopted in the aforesaid formulas about the depth of each virtual sampling gate. - In each of the aforesaid embodiments, there are a plurality of preset parameters. These parameters may be set flexibly according to actual situations. They may be directly set by a user through a user interface, or may be set by default by an ultrasonic system.
- By means of the aforesaid process of PRF automatic optimization, the time required for the signal to return to the probe from target position is effectively not fallen into or close to the range of dead time, thus avoiding a low SNR of the Doppler spectrum.
FIG. 5 shows the image before and after optimization, whereinFIG. 5 a is the Doppler spectrum before optimization andFIG. 5 b is the optimized Doppler spectrum. FromFIG. 5 , after HPRF is optimized, the Doppler signals are markedly enhanced, the SNR is improved, and the contour of the velocity information is more clearly shown. -
FIG. 6 illustrates anapparatus 600 for HPRF Doppler ultrasonic imaging according to one embodiment. Theapparatus 600 may include acollecting unit 620, asetting unit 630, and an evaluatingunit 650. In other embodiments, theapparatus 600 optionally further includes a judging unit ofsampling gate size 640 and/or a reinforcingjudgment unit 660. In one embodiment, the collectingunit 620 is configured to implementstep 100 ofFIG. 1 ; thesetting unit 630 is configured to implementstep 102; the judging unit ofsampling gate size 640 is configured to implementstep 104; the evaluatingunit 650 is configured to implementstep 104; and the reinforcingjudgment unit 660 is configured to implementstep 108. - In one embodiment, the
apparatus 600 for HPRF Doppler ultrasonic imaging optionally further includes aninterface 680 and anoutput unit 690. Theinterface 680 may be configured to collect the information selected by the user, i.e., PRF level, SVD, and the sampling gate size. Theoutput unit 690 may be configured to output the best PRF value. - According to another embodiment shown in
FIG. 7 , theapparatus 600 for HPRF Doppler ultrasonic imaging comprises aprocessing unit 613, i.e., DSP or CPU. Theprocessing unit 613 may be a single unit or a plurality of units for implementing the above-mentioned steps. Additionally, the apparatus optionally further includes aninterface 680 and anoutput unit 690. Moreover, theapparatus 600 further comprises at least one non-transitory computer-readable medium 610 which may be in form of non-volatile memory, i.e., EEPROM, flash or disk drive. The computer-readable medium 610 may include acomputer program 611 including program codes, which, when executed by theprocessing unit 613, implements the steps shown inFIG. 1 and/orFIG. 4 . - The
computer program 611 in theapparatus 600 may include acollecting module 611 a for implementingstep 100, asetting module 611 b for implementingstep 102, a judging module ofsampling gate size 611 c for implementingstep 104, an evaluatingmodule 611 d for implementingstep 106, and a reinforcingjudgment module 611 e for implementingstep 108. In other words, when running themodules 611 a-611 e on theprocessing unit 613, the modules respectively correspond to theunits FIG. 6 . - In one embodiment, the
apparatus 600 for HPRF Doppler ultrasonic imaging can be realized in a Doppler ultrasonic imaging system with software, hardware, firmware, or their combination, as shown inFIG. 8 . - Though present disclosure has been described in details by way of aforesaid embodiments, the invention is not limited by these embodiments. It can be understood by those skilled in the art that various modifications, equivalent substitutions, and changes can be made, which should be considered as within the protective scope of the invention defined by the attached claims. For example, the various steps or modules mentioned above may be divided into two or more steps or modules, or multiple steps or modules may be merged into one step or module. Additionally, the terms mentioned in the claims and/or descriptions are used for convenience and should not be construed as limiting. Further, the phrase “one embodiment” may represents different embodiments.
Claims (18)
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010044456.X | 2010-01-21 | ||
CN201010044456.XA CN102133107B (en) | 2010-01-21 | 2010-01-21 | Method and device for improving HPRF (High-Pulse Repetition Frequency) performance in Doppler ultrasonic imaging |
CN201010044456 | 2010-01-21 |
Publications (2)
Publication Number | Publication Date |
---|---|
US20110178406A1 true US20110178406A1 (en) | 2011-07-21 |
US8647277B2 US8647277B2 (en) | 2014-02-11 |
Family
ID=44278047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/007,437 Active 2032-02-10 US8647277B2 (en) | 2010-01-21 | 2011-01-14 | Methods and apparatus for HPRF doppler ultrasonic imaging |
Country Status (2)
Country | Link |
---|---|
US (1) | US8647277B2 (en) |
CN (1) | CN102133107B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120157848A1 (en) * | 2010-12-16 | 2012-06-21 | Samsung Medison Co., Ltd. | Performing receive-focusing based on mid-point algorithm in ultrasound system |
US20160220227A1 (en) * | 2015-01-30 | 2016-08-04 | Samsung Medison Co., Ltd. | Ultrasound system and method of providing guide for improved hprf doppler image |
JP6145202B1 (en) * | 2016-07-15 | 2017-06-07 | 株式会社日立製作所 | Ultrasonic diagnostic equipment |
CN112690826A (en) * | 2019-10-23 | 2021-04-23 | 通用电气精准医疗有限责任公司 | Method and system for automatically providing artifact warning in pulsed wave doppler imaging |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103356242A (en) * | 2012-03-29 | 2013-10-23 | 深圳市蓝韵实业有限公司 | Pulse repetition frequency adjustment method and high pulse Doppler imaging device |
CN102631224B (en) * | 2012-04-24 | 2014-07-02 | 无锡祥生医学影像有限责任公司 | Doppler automatic optimization method and device for touch screen ultrasonic diagnostic apparatus |
CN106236148B (en) * | 2016-07-29 | 2019-06-28 | 飞依诺科技(苏州)有限公司 | A kind of pulse recurrence frequency determines method and device |
CN106859701A (en) * | 2017-02-13 | 2017-06-20 | 深圳安盛生物医疗技术有限公司 | The velocity scale adjusting method and device of a kind of high impulse doppler imaging |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5961462A (en) * | 1998-05-18 | 1999-10-05 | Atl Ultrasound | Ultrasonic doppler imaging at high frame rates of display |
US6447455B2 (en) * | 2000-07-08 | 2002-09-10 | Medison Co., Ltd. | Ultrasound diagnostic apparatus and method for measuring blood flow velocity using doppler effect |
US20070164898A1 (en) * | 2005-12-01 | 2007-07-19 | General Electric Company | Method and apparatus for automatically adjusting spectral doppler gain |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2964472D1 (en) * | 1978-06-16 | 1983-02-17 | Hoffmann La Roche | Ultrasonic imaging apparatus and use thereof for medical ultrasonic diagnostics |
US5325858A (en) * | 1992-01-10 | 1994-07-05 | Kabushiki Kaisha Toshiba | Ultrasonic imaging system capable of varying maximum pulse repetition frequency depending upon pulse transmission frequency |
JP3696763B2 (en) | 1999-11-05 | 2005-09-21 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ultrasound imaging device |
WO2006080011A2 (en) * | 2005-01-25 | 2006-08-03 | Ramot At Tel Aviv University Ltd. | Using pulsed-wave ultrasonography for determining an aliasing-free radial velocity spectrum of matter moving in a region |
US20080288218A1 (en) * | 2005-10-27 | 2008-11-20 | Koninklijke Philips Electronics, N.V. | Using Tissue Acceleration to Create Better Dti Waveforms (Doppler Tissue Imaging) for Crt (Cardiac Resynchronization Therapy) |
JP4820210B2 (en) | 2006-05-11 | 2011-11-24 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ultrasonic diagnostic apparatus and ultrasonic diagnostic image generation method |
US8142361B2 (en) | 2007-08-28 | 2012-03-27 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Method and apparatus for automatic optimization of doppler imaging parameters |
-
2010
- 2010-01-21 CN CN201010044456.XA patent/CN102133107B/en active Active
-
2011
- 2011-01-14 US US13/007,437 patent/US8647277B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5961462A (en) * | 1998-05-18 | 1999-10-05 | Atl Ultrasound | Ultrasonic doppler imaging at high frame rates of display |
US6447455B2 (en) * | 2000-07-08 | 2002-09-10 | Medison Co., Ltd. | Ultrasound diagnostic apparatus and method for measuring blood flow velocity using doppler effect |
US20070164898A1 (en) * | 2005-12-01 | 2007-07-19 | General Electric Company | Method and apparatus for automatically adjusting spectral doppler gain |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120157848A1 (en) * | 2010-12-16 | 2012-06-21 | Samsung Medison Co., Ltd. | Performing receive-focusing based on mid-point algorithm in ultrasound system |
US9091760B2 (en) * | 2010-12-16 | 2015-07-28 | Samsung Medison Co., Ltd. | Performing receive-focusing based on mid-point algorithm in ultrasound system |
US20160220227A1 (en) * | 2015-01-30 | 2016-08-04 | Samsung Medison Co., Ltd. | Ultrasound system and method of providing guide for improved hprf doppler image |
KR20160094221A (en) * | 2015-01-30 | 2016-08-09 | 삼성메디슨 주식회사 | Ultrasound System And Method For Providing Guide To Improve HPRF Doppler Image |
US11033247B2 (en) * | 2015-01-30 | 2021-06-15 | Samsung Medison Co., Ltd. | Ultrasound system and method of providing guide for improved HPRF doppler image |
KR102387708B1 (en) | 2015-01-30 | 2022-04-19 | 삼성메디슨 주식회사 | Ultrasound System And Method For Providing Guide To Improve HPRF Doppler Image |
JP6145202B1 (en) * | 2016-07-15 | 2017-06-07 | 株式会社日立製作所 | Ultrasonic diagnostic equipment |
JP2018007952A (en) * | 2016-07-15 | 2018-01-18 | 株式会社日立製作所 | Ultrasonic diagnostic apparatus |
CN112690826A (en) * | 2019-10-23 | 2021-04-23 | 通用电气精准医疗有限责任公司 | Method and system for automatically providing artifact warning in pulsed wave doppler imaging |
Also Published As
Publication number | Publication date |
---|---|
CN102133107A (en) | 2011-07-27 |
US8647277B2 (en) | 2014-02-11 |
CN102133107B (en) | 2014-10-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8647277B2 (en) | Methods and apparatus for HPRF doppler ultrasonic imaging | |
US8142361B2 (en) | Method and apparatus for automatic optimization of doppler imaging parameters | |
US20190254629A1 (en) | Method and system for displaying ultrasonic elastic measurement | |
US4575799A (en) | Ultrasonic living body tissue characterization system | |
CN104125804B (en) | The method of operating of ultrasound observation apparatus and ultrasound observation apparatus | |
US20160074017A1 (en) | Ultrasound optimization method and ultrasonic medical device therefor | |
EP3153884A1 (en) | Detection apparatus, fish finder, and radar | |
US10324174B2 (en) | Two dimension and three dimension imaging with coded pulses based on speed changes of sound/ultrasound | |
CN104411250A (en) | Ultrasonic observation apparatus, operation method of the same, and computer readable recording medium | |
US20220206111A1 (en) | Non-contact method for detectiing and distinguishing human and animal based on ir-uwb bio-radar signal | |
CN106333707A (en) | Ultrasonic Doppler spectrum calibration method, device and ultrasonic diagnostic system | |
US10548572B2 (en) | Ultrasound processing device | |
EP2103954A1 (en) | Adaptive clutter signal filtering in an ultrasound system | |
CN110507360B (en) | Shear wave imaging method and system | |
US20150216500A1 (en) | Ultrasonic measurement apparatus and ultrasonic measurement method | |
CN109259801A (en) | A kind of shearing wave elastograph imaging method and device | |
CN106859701A (en) | The velocity scale adjusting method and device of a kind of high impulse doppler imaging | |
JP2020536668A (en) | Ultrasonic shear wave imaging with patient-adaptive shear wave generation | |
CN105911537B (en) | A method of reducing active sonar detection blind area | |
US11369341B2 (en) | Method for characterising bone using ultrasonic waves | |
RU2003101179A (en) | METHOD FOR AUTOMATIC SUPPORT OF A MANEUVERING GOAL IN THE ACTIVE LOCATION OF A HYDROACOUSTIC OR RADAR COMPLEX | |
US20220296209A1 (en) | Fetal ultrasound processing unit | |
JP6793077B2 (en) | Ultrasonic diagnostic equipment and waveform processing method | |
JP2001147267A (en) | Target identification device | |
RU2390724C2 (en) | Method for light-range finding |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS CO., LTD, Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:SHUANGSHUANG, LI;LEI, LI;REEL/FRAME:025647/0604 Effective date: 20110113 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
FEPP | Fee payment procedure |
Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1552); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 8 |
|
AS | Assignment |
Owner name: SHENZHEN MINDRAY ANIMAL MEDICAL TECHNOLOGY CO., LTD., CHINA Free format text: LICENSE;ASSIGNOR:SHENZHEN MINDRAY BIOMEDICAL ELECTRONICS CO., LTD.;REEL/FRAME:060440/0796 Effective date: 20220601 |