WO2011089898A1 - 変位推定方法、変位推定装置 - Google Patents
変位推定方法、変位推定装置 Download PDFInfo
- Publication number
- WO2011089898A1 WO2011089898A1 PCT/JP2011/000259 JP2011000259W WO2011089898A1 WO 2011089898 A1 WO2011089898 A1 WO 2011089898A1 JP 2011000259 W JP2011000259 W JP 2011000259W WO 2011089898 A1 WO2011089898 A1 WO 2011089898A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- displacement
- signal
- window
- window size
- displacement estimation
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0609—Display arrangements, e.g. colour displays
- G01N29/0645—Display representation or displayed parameters, e.g. A-, B- or C-Scan
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/22—Details, e.g. general constructional or apparatus details
- G01N29/225—Supports, positioning or alignment in moving situation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4454—Signal recognition, e.g. specific values or portions, signal events, signatures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/50—Processing the detected response signal, e.g. electronic circuits specially adapted therefor using auto-correlation techniques or cross-correlation techniques
-
- 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/8977—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
-
- 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/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52038—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
-
- 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/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/42—Details of probe positioning or probe attachment to the patient
- A61B8/4245—Details of probe positioning or probe attachment to the patient involving determining the position of the probe, e.g. with respect to an external reference frame or to the patient
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/044—Internal reflections (echoes), e.g. on walls or defects
Definitions
- the present invention introduces a method for estimating displacement using ultrasonic signals. This can be used in applications where displacement estimation is required, either in the final result or in an intermediate step for further processing, and used in medical and industrial ultrasound equipment Can do.
- An ultrasonic device is a non-invasive diagnosis of a medium by transmitting and receiving high-frequency mechanical sound waves to and from a medium to be scanned.
- the transducer of such a device transmits ultrasonic waves to the medium to be scanned.
- Ultrasonic waves interact with the lower layer structure of the medium by scattering and reflecting.
- the lower layer structure refers to, for example, a structure inside the medium rather than the surface of the medium.
- the lower layer structure is, for example, a blood vessel in a human body.
- Scattered and reflected ultrasound contains useful information of the underlying structure and is received by the transducer. Thereafter, it is processed by an ultrasonic device and presented to the user as image information, for example.
- FIG. 8 is a diagram showing an ultrasonic RF signal.
- RF Radiofrequency Signal
- B mode luminance mode
- Doppler images etc.
- other types of data such as luminance mode (B mode) images, Doppler images, etc. can be derived from ultrasonic RF signals.
- One of many applications is to analyze the motion of the underlying structure of the scanned media.
- the Doppler effect is used as a simple method for estimating the direction and strength of structure movement.
- Doppler ultrasound is suitable for applications where precision is not strictly required, such as blood flow monitoring in medical ultrasound.
- the estimated displacement quality largely depends on the quality and resolution of the B-mode image.
- the cross-correlation technique increases the amount of calculation and can only estimate the displacement corresponding to a plurality of sampling points.
- a micrometer displacement usually corresponds to only a small part of one sampling interval.
- Autocorrelation depends on the phase information of the quadrature demodulated signal (also known as baseband signal) of the received RF signal.
- This method has an advantage that the displacement corresponding to the sub-sample of the RF signal described in [2] can be estimated.
- the estimated displacement is biased more than the high signal power region.
- the “coarse to fine” approach uses different window regions for different estimation stages, so that the displacement is roughly estimated in the first stage and more precisely estimated in the second stage. Improve accuracy.
- the second estimation stage is executed and combined with the estimation result of the first stage to perform signal warping based on the estimation result of the first stage.
- Patent Documents 1 to 6 are known as conventional examples.
- the current method has problems in terms of accuracy and displacement resolution.
- the method based on cross-correlation cannot estimate a small displacement corresponding to the sub-sampling interval without signal interpolation.
- the cross-correlation-based method has high accuracy but low displacement resolution.
- Signal interpolation can be applied for the purpose of improving the displacement resolution of the cross-correlation method.
- the accuracy of this method depends on the quality of the interpolation algorithm.
- the method based on autocorrelation can estimate small displacements, but it is very easy to generate noise.
- a two-stage estimation method can be used to improve accuracy.
- a comprehensive iterative technique needs to incorporate techniques that lead to the convergence of iterative estimation.
- An object of the present invention includes providing a displacement estimation device (displacement estimation method) capable of accurately estimating a tissue displacement, and thus distinguishing a malignant tumor from a benign tumor or a normal tissue based on the estimated displacement.
- a displacement estimation device capable of accurately estimating a tissue displacement, and thus distinguishing a malignant tumor from a benign tumor or a normal tissue based on the estimated displacement.
- Another object of the present invention includes providing a displacement estimation apparatus that can converge an iterative displacement estimation process with a small number of iterations and with high accuracy.
- a displacement estimation method of the present invention is a displacement estimation method that repeatedly estimates displacement using an ultrasonic signal, and scans at least one ultrasonic signal to obtain a medium. Transmitting the ultrasonic signal reflected from the scanned medium; calculating a window size; calculating a window boundary based on the calculated window size; Using the window by the calculated boundary to estimate a displacement at each depth of the ultrasound signal; warping the ultrasound signal based on the estimated displacement; and warped Using the ultrasonic signal, a step for leading the convergence of the displacement estimation method so that the correlation value of the ultrasonic signal is increased.
- a displacement estimation method comprising and.
- the inventor teaches a method for estimating displacement using ultrasonic waves using an iterative estimation method.
- the estimated displacement does not converge if the estimated window size remains the same or is not properly configured for all iterative estimation rounds.
- the window calculation method calculates the window boundary to overcome the effects of noise and bias caused by the non-uniform distribution of signal power and the part of calculating the window size to guide the convergence of the iterative estimation And two parts.
- signal power is used to determine the boundary of the estimated window at each estimated position.
- the window is usually symmetric around its position.
- the present invention guarantees the accuracy of displacement estimation using a comprehensive iterative technique, including a method for inducing convergence, a method for overcoming elements that limit accuracy, and a method for evaluating the quality of results.
- FIG. 10 is a diagram showing an improvement result of the displacement estimation after the backtrack algorithm is applied to the simulation.
- FIG. 11 is a diagram showing the same result in the phantom experiment.
- FIG. 10 is a diagram showing that the estimated displacement (displacement 1001) in the first round of estimation deviates from the simulated displacement profile (profile 1000).
- the estimated displacement fits very well with the simulated profile.
- FIG. 11 is a diagram showing an estimation result by a phantom experiment in which a global displacement is generated.
- a single line is selected (see profile 1100).
- the final round estimation result (see data 1103) represents the global movement much better than the first round estimation result.
- An object of the present invention includes providing a displacement estimation device (displacement estimation method) capable of accurately estimating a tissue displacement, and thus distinguishing a malignant tumor from a benign tumor or a normal tissue based on the estimated displacement.
- a displacement estimation device capable of accurately estimating a tissue displacement, and thus distinguishing a malignant tumor from a benign tumor or a normal tissue based on the estimated displacement.
- Another object of the present invention includes providing a displacement estimation apparatus that can converge an iterative displacement estimation process with a small number of iterations and with high accuracy.
- FIG. 1 is a diagram showing a displacement estimation method of the present invention.
- FIG. 2 is a diagram showing a window calculation method according to the present invention.
- FIG. 3 is a diagram illustrating an example of window size calculation.
- FIG. 4 is a diagram illustrating an example of a symmetric window and an energy equalization window.
- FIG. 5 is a diagram illustrating an RF signal warping method according to the present invention.
- FIG. 6 is a diagram illustrating RF signal warping and RF sample delay.
- FIG. 7 is a diagram illustrating a method that leads to convergence of the present invention.
- FIG. 8 is a diagram showing an ultrasonic RF signal.
- FIG. 9 is a diagram illustrating displacement estimation using autocorrelation.
- FIG. 9 is a diagram illustrating displacement estimation using autocorrelation.
- FIG. 10 is a diagram showing an example of a displacement estimation result using the present invention for the simulation.
- FIG. 11 is a diagram showing an example of a displacement estimation result using the present invention for the phantom experiment.
- FIG. 12 is a diagram illustrating the displacement of the lower layer structure and the influence on the RF signal.
- FIG. 13 is a block diagram of the apparatus.
- the displacement estimation method of the embodiment uses a displacement estimation method (displacement estimation apparatus in FIG. 13) that repeatedly estimates displacement (see displacement 1207 in FIG. 12) using an ultrasonic signal (ultrasonic signal 1201 s in FIG. 13). 1X).
- the displacement magnitude is identified with relatively high accuracy.
- At least one ultrasonic signal is scanned and transmitted to the medium (the medium 1200M (FIG. 12) and the lower layer structure (displacement measurement target) 1203x (see FIGS. 12 and 13) in the medium 1200M). (Refer to the transmission unit 1X1 in the ultrasonic processing unit 1Xa in FIG. 13).
- ultrasonic signals are transmitted to a plurality of positions in the line direction (direction 81L in FIG. 8), and scanning of the plurality of positions is performed.
- the ultrasonic signal (ultrasonic signal 1201s) reflected from the scanned medium is received (see the receiving unit 1X2).
- the ultrasonic processing unit 1Xa is, for example, a probe (probe).
- window size 1X3M in FIG. 13 for example, window size (winSize) in columns (a) and (b) in FIG. 3) is calculated (size calculation included in information processing unit 1Xb in FIG. 13). Part 1X3).
- the window size in each round is determined.
- the size calculation unit 1X3 may include, for example, the window calculation unit 102 (FIG. 1), or may be at least a part of the window calculation unit 102.
- the window size in the n-th round is determined (n ⁇ 2) based on the correlation value in the (n ⁇ 1) -th round (correlation value 1 ⁇ 7M in FIG. 13).
- the window boundary (boundary information 1X4M: the window start edge (winStart) and the window end (winEnd) shown in the fields (a) and (b) of FIG. 4) Etc.) (see the boundary calculation unit 1X4).
- boundary calculation unit 1X4 may be at least a part of the window calculation unit 102 (FIG. 1), for example.
- start and end points are specified in which the width between the start and end points is the calculated window size.
- the displacement (displacement 1XbM) at each depth of the ultrasonic signal is estimated using the calculated window based on the boundary (see the estimation unit 1X5).
- the estimation unit 1X5 may be, for example, at least a part of the displacement estimation unit 101 (FIG. 1), for example.
- the ultrasonic signal is warped (see warping unit 1X6).
- the warping unit 1X6 may be at least a part of the RF signal warping unit 100 (FIG. 1), for example.
- the convergence of the displacement estimation method is guided so that the correlation value (correlation value 1X7M) of the ultrasonic signal becomes large (see the convergence control unit 1X7).
- the convergence control unit 1X7 may be at least a part of the convergence guide unit 103, for example.
- the magnitude of the displacement estimated in the round with the maximum correlation value calculated from the warped ultrasonic signal is specified as the most accurate magnitude.
- the calculated window size is changed so that the window size used for a plurality of consecutive rounds gradually decreases.
- the window size calculated in each round is smaller than the window size calculated in the previous round, and the calculated window size is changed to a smaller size.
- the window size is then calculated so that the signal energy in all (depth) windows is equal.
- window size in the n-th round is calculated as the correlation value (correlation value in the (n-1) -th round) of the ultrasonic signal increases.
- the relationship between the evaluation value of the degree of convergence in each round, such as the correlation value, and the window size can be determined with reference to a predetermined relational expression based on experimental values in the living body, or a table indicating the correspondence between the two. .
- the window size that the determined correspondence corresponds to the correlation value in the (n ⁇ 1) -th round is specified as the window size in the n-th round.
- the ultrasonic signal 1201s may be transmitted to the lower layer structure 1203x (FIG. 12) that is a measurement target (for example, a malignant tumor with new blood vessels), and the transmitted ultrasonic signal 1201s may be received.
- a measurement target for example, a malignant tumor with new blood vessels
- the identified delay time 1208 is a delay time between the position of the first pulse (see the RF signal 1204 in the column (b) of FIG. 12) and the position of the second pulse (see the RF signal 1205). It is.
- the first pulse is a pulse before the lower layer structure 1203x is displaced 1207 in the received ultrasonic signal 1201s.
- the second pulse is a pulse after the measurement object has displaced 1207.
- the displacement 1207 is determined based on the specified delay time 1208.
- the lower layer structure 1203x is determined as a malignant tumor (cancer), and when it is determined that the movement is not the movement, it is determined as a benign tumor or a normal tissue. May be.
- the processing based on the window size may be processing in which, for example, the correlation value in the n-th round is calculated from data of the window portion of the window size in the received ultrasonic signal 1201s.
- FIG. 1 is a diagram showing a displacement estimation method.
- FIG. 1 The main embodiment of the present invention is shown in FIG. 1
- RFSig (d, l, f) is an RF signal obtained from the ultrasonic unit, d is the depth direction (direction 81D in FIG. 8), l is the line direction (direction 81L), and f is the frame direction (direction 81F). Represents.
- One line (see 801 shown) in one frame of such RF signal is generated by a mechanism in which the ultrasonic transducer 800 transmits a pulse in the scanning direction (direction 802) at the position of the line.
- the pulse interacts with the underlying structure along the path while being reflected and scattered, and the ultrasonic transducer 800 receives the reflected and scattered signal.
- the signal is converted into a corresponding line by the ultrasonic unit.
- the same mechanism generates multiple lines at various positions and multiple frames at various time instances.
- FIG. 12 shows the mechanism by which the displacement of the underlying structure of the scanned medium is reflected as the delay of the RF signal.
- the ultrasonic transducer 1200 transmits a pulse toward the scanned medium (see pulse 1201).
- This pulse propagates through the medium until it reaches the reflection boundary of the lower layer structure 1203 (1203x).
- This pulse causes the pulse to be reflected back towards the transducer (see reflected pulse 1202).
- the transducer receives this pulse and converts it into an RF signal 1204 (see the first row of the table in FIG. 12 (b)).
- the position of the pulse in the entire RF signal 1204 indicates the time required for the transmitted pulse 1201 to travel toward the lower layer structure 1203 and reflect back.
- the RF signal 1205 (the third row of the table in the column (b) in FIG. 12) is created by the same mechanism. Is displaced to a new position (see the lower layer structure 1206 after the displacement).
- Displacement x causes the transmitted pulse (pulse 1201) to propagate and take longer to reflect at the boundary of the underlying structure, resulting in RF signal delay (delay time 1208). reference).
- the direction of the displacement x is the direction of the direction 1203d shown in FIG.
- the displacement (displacement 1207) can be estimated.
- a plurality of reflected pulses are generated.
- the delay of the reflected pulse in the RF signal is different because the amount of displacement of the plurality of lower layer structures is different from each other.
- the present invention presents a method for accurately estimating the displacement of the lower layer structure by repeatedly estimating.
- the present invention comprises the following main blocks.
- the RF signal warping unit 100 (see FIG. 5), the displacement estimation unit 101 (see FIG. 9), the window calculation unit 102 (see FIG. 2), and the convergence guide unit 103 (see FIG. 7) And an output switch 104.
- the RF signal warping unit 100 warps in the target estimation round (n-th round) based on the displacement estimation result dispRound (d, l, round) of the previous estimation round (n-1 round). In this warping, the selected line in the selected frame of the RF signal is warped for the target estimation round.
- the output of this block is the warped RF signal RFSigWarp (d, l, f, round) used in the estimation round.
- dispRound (d, l, 0) is initialized to zero, so that it is obtained from RFSigWarp (d, l, f, 1) and the ultrasonic unit.
- the RF signals RFSig (d, l, f) ⁇ are substantially the same.
- a pre-defined set of frames can be selected as input for this block to estimate displacement.
- two frames are selected and are represented by f1 and f2.
- RF signal warping The purpose of RF signal warping is to modify one RF signal to match another RF signal based on the displacement estimated in the previous round.
- a new displacement is created by estimating the residual displacement from the corrected RF signal and adding it to the previous round displacement.
- the estimated displacement will almost match the modified RF signal and the residual displacement will converge to zero.
- the displacement estimation unit 101 estimates the displacement from the RF signal RFSigWarp (d, l, f, round) by deriving from the time delay of the RF signal.
- a preferred method for performing this task is, but is not limited to, an autocorrelation method where the estimation window at each depth is described by winStart (d, l) and winEnd (d, l). .
- the window calculation unit 102 calculates the estimated window parameters winStart (d, l) and winEnd (d, l) used for the selected line at each depth.
- the window calculation unit 102 acquires the RF signal RFSigWarp (d, l, f, round), and obtains the target round number curRound and the RF signal difference residue (d, l, round) from the convergence guide unit 103. Get as input.
- the convergence guide unit 103 performs processing for guiding convergence according to the present technology.
- the difference between the warped RF signals is calculated to determine the estimated displacement quality.
- the output of this block includes the target estimated round number curRound, the round number minResidueRound (d, l) with the smallest RF signal difference, and the RF signal difference residue (d, l, round).
- the output switch 104 acquires curRound and minResidueRound (d, l) as inputs.
- the most accurate estimation result is selected as the final output.
- the most accurate estimation result is selected as the input of the RF signal warping unit 100 for use in subsequent estimation rounds, and residue (d, l, round) and curRound are input as the displacement estimation unit. 101.
- FIG. 5 is a diagram showing an RF signal warping method.
- the RF signal warping unit 100 of the main embodiment is realized by the method shown in FIG.
- the most accurate displacement estimation result dispRound (d, l, minResidueRound (d, l)) is used to adapt RFSig (d, l, f1) to RFSig (d, l, f2).
- the delay value (see the delay time 1208 in FIG. 12) is calculated by the delay value calculation unit 500 as follows.
- fs is the sampling frequency of the RF signal
- c is the speed of sound in the scanned medium.
- the displacement in the scanned medium is reflected as a delay in the resulting RF signal as a result.
- the delay in the RF sample delay unit 501 is performed using signal interpolation.
- the delay in the RF sample delay unit 501 is performed using a fractional delay filter.
- FIG. 6 is a diagram showing RF signal warping and RF sample delay.
- the RF signal warping process is shown in FIG.
- RF signal warping The purpose of RF signal warping is to nullify this effect by delaying one RF signal in the opposite direction to the estimated displacement.
- the warped RF signals are expected to be compatible with each other.
- the delayValues (d, l) calculated from the delay value calculation unit 500 represents a value by which each sample in one RF signal is delayed to achieve this.
- Two RF signals 600 and 601 are used for this illustration.
- Data 602 indicates delayValues (d, l) calculated from the delay value calculation unit 500 after a certain estimation round.
- a value at a specific depth of delayValues (d, l) represents an amount necessary for delay for the sample at the same depth in the RF signal 600 to be adapted to the delay value calculation unit 601.
- delayValues (d, l) is not always an integer.
- the displacement estimation unit 101 By applying the displacement estimation unit 101 to the RF signal, the displacement used in the first round and the residual displacement used in the subsequent round are obtained, and these are combined to obtain the final displacement.
- FIG. 9 A preferred embodiment is shown in FIG. 9, but is not limited thereto.
- FIG. 9 shows displacement estimation using autocorrelation.
- the IQ demodulator 900 converts the RF signal RFSigWarp (d, l, f, round) into a baseband signal IQSig (d, l, f).
- the autocorrelation calculation unit 901 calculates autocorr (d, l) according to the following equation.
- conj () represents a complex conjugate
- the autocorr (d, l) is converted into a displacement by the displacement calculator 902 according to the following equation.
- arg () is a function for calculating the argument of a complex number.
- DispRound (d, l, 0) is initialized to zero by a round starting from 1 (representing the first round).
- FIG. 2 is a diagram showing a window calculation method.
- the window calculation unit 102 determines different estimation windows for each depth of each line in each estimation round.
- FIG. 3 is a diagram showing a specific example of window size calculation.
- the window size calculation unit 200 determines the window size for each estimation round.
- the window size for each round changes as the function of rounds decreases.
- the window size is directly related to the estimated number of rounds by the decreasing function winSizeRound (round).
- the window sizes for all depths are specified based on the following.
- the “fixed windowing scheme” is a process of selecting the same value of winSize (d, l) that becomes an estimated window size at all depths used for displacement estimation.
- the advantage of the “fixed windowing scheme” is that it requires less computation (relatively less) and is suitable when the RF signal energy is regularly distributed along the depth direction. .
- the window size at each depth is based on the upper limit maxWin (round), the lower limit minWin (round), and the RF signal power sigPow (d, l). Is calculated.
- MaxWin (round) and minWin (round) are round decreasing functions.
- the “instantaneous power windowing scheme” requires more computation than the “fixed windowing scheme”, but it is more especially when the RF signal energy is fairly non-uniform along the depth direction. High accuracy.
- a larger estimation window is specified for regions with low signal power, and for regions with large signal power to distribute displacement in more detail.
- a smaller estimation window is specified.
- the window size at all depths may be calculated such that the signal energy in all windows is kept the same, which is the “constant energy window” Ing scheme ".
- the signal power calculator 201 calculates the average power of the RF signal along the frame direction.
- the signal power calculation unit 201 acquires the warped RF signal RFSigWarp (d, l, f, round) and the target round number curRound as inputs.
- N is the total number of frames selected for use in the calculation.
- FIG. 4 is a diagram showing examples of a symmetric window and an energy equalization window.
- a window boundary at each depth of each line is calculated by the window boundary calculation unit 202 (FIG. 2).
- winSize (d, l) calculated from the window size calculation unit 200 is an input, and winStart (d, l) and winEnd (d, l) are expressed by the following equations:
- the window is calculated to be symmetric around the corresponding depth (referred to as a “symmetric window”).
- ROUND () indicates a process of rounding to the nearest integer.
- the advantage of the symmetric window is that it does not require a large amount of computation.
- the operation is as follows.
- winSize (d, l) calculated from the window size calculation unit 200 is an input to the window boundary calculation unit 202.
- winStart (d, l) and winEnd (d, l) are calculated so that the signal energies on both sides of each depth are equal (referred to as “energy equalization window”) according to the following equation.
- the “Energy Equal Window” requires more computation than the “Symmetric Window”, but the signal power is non-uniformly distributed along the depth direction, which caused the bias in the conventional windowing technology. The influence of can be suppressed.
- FIG. 7 is a diagram showing a method for leading to convergence.
- the convergence guide unit 103 acquires the warped RF signal RFRFSigWarp (d, l, f, round) as an input of the convergence guide unit 103 as shown in FIG.
- the RF signal difference calculation unit 700 calculates the difference between the warped RF signals, and outputs the difference in each estimation round as residue (d, l, round) for each depth of each line.
- the main purpose of the RF signal warping unit 100 is to modify one of the RF signals so as to be compatible with the other, this residual (d, l, round) It is a measure of whether or not it is compatible.
- the RF signal warping unit 100 is used as an input to the RF signal warping unit 100. , DispRound (d, l, round) is acquired.
- warped signals are also used as inputs for displacement estimation in subsequent rounds and are therefore represented as RFSigWarp (d, l, f, round + 1).
- residue (d, l, round) is calculated from RFSigWarp (d, l, f, round + 1).
- f1 and f2 represent two frames selected for use in estimation.
- the same signal difference value is designated for all the depths of one line, and this indicates that the RF signal of the corresponding line is generally compatible.
- f1 and f2 represent two frames selected for use in estimation.
- each depth of each line has a separate signal difference value, which indicates that the RF signal is locally adapted.
- the minimum residual calculation unit 701 confirms the number of rounds with the minimum residual for each depth and line in residue (d, l, round) based on the following equation.
- the outputs of the convergence guide unit 103 are the target round number curRound, the round number minResidueRound (d, l) with the smallest residual, and the signal difference residue (d, l, round).
- the output switch 104 acquires curRound and minResidueRound (d, l) as inputs.
- the output switch 104 selects the most accurate displacement estimation result dispRound (d, l, minResidueRound (d, l)) as the final output dispOut.
- dispRound (d, l, curRound) instead of subjecting the target round dispRound (d, l, curRound) to displacement estimation, dispRound (d, l, minResoundRound () used as the input of the RF signal warping unit 100 for the next estimation round d, l)) is estimated for displacement.
- the pulsation direction of biological tissues such as blood vessels and tumor tissues is not limited to one dimension, and the direction varies depending on the location.
- the displacement in a plurality of directions at a specific position may be obtained, and the two-dimensional or three-dimensional displacement of the tissue may be measured.
- the ultrasonic transducer changes the direction of the probe using a linear probe arranged in one dimension, or the ultrasonic transducer is on a two-dimensional array.
- a method of performing transmission / reception in a plurality of directions by beam forming using a matrix probe arranged in the position is possible.
- the position and orientation of the probe is acquired by tracking the position of the probe with a position sensor such as a magnetic sensor, acceleration sensor, gyroscope, or camera, and the direction during transmission / reception is determined. Can be determined.
- a position sensor such as a magnetic sensor, acceleration sensor, gyroscope, or camera
- Displacement measurement is performed only in one direction, such as the ultrasonic depth direction, and this displacement vector is decomposed into components in the respective axial directions of a predetermined two-dimensional or three-dimensional coordinate system in a multidimensional manner. You may evaluate.
- strain measurement is a time derivative of displacement and serves as an index of tissue hardness.
- the displacement estimation accuracy is improved, and as a result, strain measurement accuracy is improved.
- a differential value is used.
- a time-change waveform of the displacement of a specific part of the tissue can also be used for determination of the tissue properties. For example, it is known that angiogenesis occurs in and around cancer during the growth process of cancer.
- the pulsation of these blood vessels or the pulsation of the surrounding tissue accompanying the pulsation of the blood vessels is regarded as a temporal change in the displacement amount, and it is determined whether the change pattern and the amplitude are peculiar to cancer. Thus, cancer detection may be performed.
- a displacement estimator for repetitively estimating displacement using an ultrasound signal, means for transmitting at least one ultrasound signal from a transducer to a scanned medium; Means for receiving an ultrasonic signal reflected from the scanned medium; means for calculating an estimated window size; means for calculating a boundary of the estimated window based on the estimated window size; Means for estimating displacement along the depth direction of the ultrasonic signal using the estimated window, means for warping the ultrasonic signal based on the estimated displacement, and the convergence tendency of the method
- a displacement estimation device is disclosed that comprises means for inducing convergence by calculating as a difference between warped ultrasound signals.
- the disclosed method includes a method for inducing convergence, a method for overcoming the factors that limit accuracy, and a method for evaluating the quality of the results, and uses a comprehensive iterative method to ensure the accuracy of the displacement estimation.
- the prior art provides a method for increasing the accuracy and resolution of estimation.
- the present invention can be realized not only as a device, a system, an integrated circuit, etc., but also as a method that uses processing means constituting the device as steps, or as a program that causes a computer to execute these steps, It can be realized as a computer-readable recording medium such as a CD-ROM in which the program is recorded, or can be realized as information, data or a signal indicating the program.
- These programs, information, data, and signals may be distributed via a communication network such as the Internet.
- This technology introduces a method for estimating displacement using ultrasonic signals. This can be used in applications where displacement estimation is required, either in the final result or in an intermediate step for further processing, and used in medical and industrial ultrasound equipment Can do.
- a displacement estimation device that can accurately estimate the displacement can be provided.
- a displacement estimation device that can provide information suitable for discriminating between a malignant tumor and a benign tumor or a normal tissue can be provided based on the estimated displacement. .
Abstract
Description
winSize(d,l) = winSizeRound(curRound)
1X1 送信部
1X2 受信部
1X3 サイズ計算部
1X4 境界計算部
1X5 推定部
1X6 ワーピング部
1X7 収束制御部
1201s 超音波信号
1203x 下層構造
1X3M ウインドウサイズ
1X4M 境界の情報
1XbM 変位
1X7M 相関値
100 RF信号ワーピング部
103 収束ガイド部
101 変位推定部
102 ウインドウ計算部
103 収束ガイド部
104 出力スイッチ
600、1204、1205 RF信号
800、1200 超音波トランスデューサ
1200M 媒体
1208 遅延時間
100a、101a、1021a、103a、104a、104b、104c 情報
201n、201m、201a、200m、200a、202a 情報
500m、500a、501n、501a 情報
700m、700a、701a、700b 情報
900m、900a、901m、901a、902a 情報
Claims (10)
- 超音波信号を用いて、変位を反復的に推定する変位推定方法であって、
少なくとも1つの超音波信号を走査して、媒体に送信するステップと、
走査された前記媒体から反射した前記超音波信号を受信するステップと、
ウインドウサイズを計算するステップと、
計算された前記ウインドウサイズに基づいて、ウインドウの境界を計算するステップと、
計算された前記境界による前記ウインドウを用いて、前記超音波信号のそれぞれの深度における変位を推定するステップと、
推定された前記変位に基づいて、前記超音波信号をワーピングするステップと、
ワーピングされた前記超音波信号を用いて、前記超音波信号の相関値が大きくなるように、この変位推定方法の収束を導くステップとを含む
変位推定方法。 - 計算される前記ウインドウサイズは、連続する複数のラウンドに用いられる前記ウインドウサイズが、徐々に減少するように、変化される
請求項1記載の変位推定方法。 - 連続する複数の前記ラウンドの前記ウインドウサイズの上限と下限とのそれぞれを、徐々に減少するように変化させるステップと、
各深度での信号電力に基づいて、それぞれの深度における前記ウインドウサイズとして、その深度における上限と下限との間のウインドウサイズを計算するステップとを含む
請求項1記載の変位推定方法。 - 前記ウインドウサイズは、全てのウインドウにおける信号エネルギーが等しくなるように計算される
請求項1記載の変位推定方法。 - 超音波信号の相関値が大きくなるほど、より小さい前記ウインドウサイズが計算される
請求項1記載の変位推定方法。 - 前記境界は、当該境界による前記ウインドウが、対応する深度を中心に、対称のウインドウになるように、当該深度の両側に延びる
請求項1記載の変位推定方法。 - 前記境界は、対応する深度での、当該境界による前記ウインドウの両側における信号エネルギーが互いに等しくなるように、対応する前記深度の両側に延びる
請求項1記載の変位推定方法。 - 前記超音波信号を、計算された遅延値により、前記超音波信号内の各サンプルを遅延させることにより、ワーピングするステップを含む
請求項1記載の変位推定方法。 - 各ラウンドの後に、ワーピングされたRF信号間の信号差分を計算するステップと、
計算された前記信号差分が最小になるラウンドを決定するステップとを含む
請求項1記載の変位推定方法。 - 超音波信号を用いて、変位を反復的に推定する変位推定装置であって、
少なくとも1つの超音波信号を走査して、媒体に送信する送信部と、
走査された前記媒体から反射した前記超音波信号を受信する受信部と、
ウインドウサイズを計算するサイズ計算部と、
計算された前記ウインドウサイズに基づいて、ウインドウの境界を計算する境界計算部と、
計算された前記境界による前記ウインドウを用いて、前記超音波信号のそれぞれの深度における変位を推定する推定部と、
推定された前記変位に基づいて、前記超音波信号をワーピングするワーピング部と、
ワーピングされた前記超音波信号を用いて、前記超音波信号の相関値が大きくなるように、当該変位推定装置が実行する変位推定方法の収束を導く収束制御部とを含む
変位推定装置。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011533902A JP5692079B2 (ja) | 2010-01-20 | 2011-01-19 | 変位推定方法、変位推定装置 |
EP11734505.8A EP2526872A4 (en) | 2010-01-20 | 2011-01-19 | OFFSET METHOD AND DISPLACEMENT METER |
US13/256,642 US8675450B2 (en) | 2010-01-20 | 2011-01-19 | Displacement estimating method and displacement estimating apparatus |
CN2011800014600A CN102361597A (zh) | 2010-01-20 | 2011-01-19 | 位移估计方法、位移估计装置 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010010336 | 2010-01-20 | ||
JP2010-010336 | 2010-01-20 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2011089898A1 true WO2011089898A1 (ja) | 2011-07-28 |
Family
ID=44306701
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2011/000259 WO2011089898A1 (ja) | 2010-01-20 | 2011-01-19 | 変位推定方法、変位推定装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8675450B2 (ja) |
EP (1) | EP2526872A4 (ja) |
JP (1) | JP5692079B2 (ja) |
CN (1) | CN102361597A (ja) |
WO (1) | WO2011089898A1 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109765540A (zh) * | 2019-02-26 | 2019-05-17 | 南京莱斯电子设备有限公司 | 一种频率步进体制米波雷达目标提取方法 |
US10314964B2 (en) | 2013-07-25 | 2019-06-11 | Fresenius Medical Care Deutschland Gmbh | Sterile tube covering for a medical tubing system |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014187270A (ja) * | 2013-03-25 | 2014-10-02 | Sony Corp | 固体撮像装置およびその製造方法、並びに電子機器 |
JP6180882B2 (ja) * | 2013-10-31 | 2017-08-16 | ソニーセミコンダクタソリューションズ株式会社 | 固体撮像装置、信号処理装置、および電子機器 |
US9878664B2 (en) * | 2015-11-04 | 2018-01-30 | Zoox, Inc. | Method for robotic vehicle communication with an external environment via acoustic beam forming |
US9494940B1 (en) | 2015-11-04 | 2016-11-15 | Zoox, Inc. | Quadrant configuration of robotic vehicles |
US9804599B2 (en) | 2015-11-04 | 2017-10-31 | Zoox, Inc. | Active lighting control for communicating a state of an autonomous vehicle to entities in a surrounding environment |
CN111127349B (zh) * | 2019-12-05 | 2022-08-12 | 河海大学 | 一种数字图像相关的非均匀位移后处理方法及装置 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6270459B1 (en) | 1998-05-26 | 2001-08-07 | The Board Of Regents Of The University Of Texas System | Method for estimating and imaging of transverse displacements, transverse strains and strain ratios |
US6277074B1 (en) | 1998-10-02 | 2001-08-21 | University Of Kansas Medical Center | Method and apparatus for motion estimation within biological tissue |
US20010034485A1 (en) | 2000-04-10 | 2001-10-25 | Tetsuya Kawagishi | Ultrasonic diagnosis apparatus |
JP2003180686A (ja) * | 2001-12-21 | 2003-07-02 | Chikayoshi Sumi | 変位ベクトル計測方法および装置並びに歪テンソル計測装置 |
JP2007029703A (ja) * | 2005-07-27 | 2007-02-08 | Medison Co Ltd | 弾性映像信号の非相関度を減少させて超音波映像を形成する方法 |
JP2007152074A (ja) * | 2005-01-21 | 2007-06-21 | Chikayoshi Sumi | 変位又は歪計測方法及び装置、速度計測方法、弾性率・粘弾性率計測装置、及び、超音波診断装置 |
US20080019609A1 (en) | 2006-07-20 | 2008-01-24 | James Hamilton | Method of tracking speckle displacement between two images |
WO2008038615A1 (fr) | 2006-09-25 | 2008-04-03 | Tohoku University | dispositif ÉCHOGRAPHIQUE |
US20090221916A1 (en) | 2005-12-09 | 2009-09-03 | The Trustees Of Columbia University In The City Of New York | Systems and Methods for Elastography Imaging |
JP2009537273A (ja) * | 2006-05-23 | 2009-10-29 | ケンブリッジ・エンタープライズ・リミテッド | 画像データ処理システム |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6571016B1 (en) * | 1997-05-05 | 2003-05-27 | Microsoft Corporation | Intra compression of pixel blocks using predicted mean |
US6117081A (en) * | 1998-10-01 | 2000-09-12 | Atl Ultrasound, Inc. | Method for correcting blurring of spatially compounded ultrasonic diagnostic images |
US6813315B1 (en) * | 2001-04-24 | 2004-11-02 | Vweb Corporation | Motion estimation using multiple search windows |
US20040034304A1 (en) | 2001-12-21 | 2004-02-19 | Chikayoshi Sumi | Displacement measurement method and apparatus, strain measurement method and apparatus elasticity and visco-elasticity constants measurement apparatus, and the elasticity and visco-elasticity constants measurement apparatus-based treatment apparatus |
US20050096538A1 (en) * | 2003-10-29 | 2005-05-05 | Siemens Medical Solutions Usa, Inc. | Image plane stabilization for medical imaging |
CN1313054C (zh) * | 2004-08-20 | 2007-05-02 | 清华大学 | 一种多尺度的生物组织位移估计方法 |
US20060047462A1 (en) * | 2004-08-31 | 2006-03-02 | Picciotto Carl E | Displacement estimation system and method |
US8211019B2 (en) | 2005-01-21 | 2012-07-03 | Chikayoshi Sumi | Clinical apparatuses |
US20070038101A1 (en) | 2005-07-27 | 2007-02-15 | Medison Co., Ltd. | Method of forming strain images |
US8012093B2 (en) | 2005-07-27 | 2011-09-06 | Medison Co., Ltd. | Ultra sound system for forming strain images with instantaneous frequency |
EP2281231B1 (en) * | 2008-03-18 | 2013-10-30 | Elliptic Laboratories AS | Object and movement detection |
-
2011
- 2011-01-19 JP JP2011533902A patent/JP5692079B2/ja active Active
- 2011-01-19 EP EP11734505.8A patent/EP2526872A4/en not_active Withdrawn
- 2011-01-19 CN CN2011800014600A patent/CN102361597A/zh active Pending
- 2011-01-19 US US13/256,642 patent/US8675450B2/en not_active Expired - Fee Related
- 2011-01-19 WO PCT/JP2011/000259 patent/WO2011089898A1/ja active Application Filing
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6270459B1 (en) | 1998-05-26 | 2001-08-07 | The Board Of Regents Of The University Of Texas System | Method for estimating and imaging of transverse displacements, transverse strains and strain ratios |
US6277074B1 (en) | 1998-10-02 | 2001-08-21 | University Of Kansas Medical Center | Method and apparatus for motion estimation within biological tissue |
US20010034485A1 (en) | 2000-04-10 | 2001-10-25 | Tetsuya Kawagishi | Ultrasonic diagnosis apparatus |
US6506158B2 (en) | 2000-04-10 | 2003-01-14 | Kabushiki Kaisha Toshiba | Ultrasonic diagnosis apparatus |
JP2003180686A (ja) * | 2001-12-21 | 2003-07-02 | Chikayoshi Sumi | 変位ベクトル計測方法および装置並びに歪テンソル計測装置 |
JP2007152074A (ja) * | 2005-01-21 | 2007-06-21 | Chikayoshi Sumi | 変位又は歪計測方法及び装置、速度計測方法、弾性率・粘弾性率計測装置、及び、超音波診断装置 |
JP2007029703A (ja) * | 2005-07-27 | 2007-02-08 | Medison Co Ltd | 弾性映像信号の非相関度を減少させて超音波映像を形成する方法 |
US20090221916A1 (en) | 2005-12-09 | 2009-09-03 | The Trustees Of Columbia University In The City Of New York | Systems and Methods for Elastography Imaging |
JP2009537273A (ja) * | 2006-05-23 | 2009-10-29 | ケンブリッジ・エンタープライズ・リミテッド | 画像データ処理システム |
US20080019609A1 (en) | 2006-07-20 | 2008-01-24 | James Hamilton | Method of tracking speckle displacement between two images |
WO2008038615A1 (fr) | 2006-09-25 | 2008-04-03 | Tohoku University | dispositif ÉCHOGRAPHIQUE |
Non-Patent Citations (4)
Title |
---|
C.PELLOT-BARAKAT ET AL.: "Ultrasound Elastography Based on Multiscale Estimations of Regularized Displacement Fields", IEEE TRANS. MED. IMAGING, vol. 23, no. 2, February 2004 (2004-02-01), pages 153 - 163, XP011106351 * |
C.SUMI ET AL.: "Fine Elasticity Imaging Utilizing the Iterative RF-echo Phase Matching Method", IEEE TRANS. ON UFFC., vol. 46, no. 1, January 1999 (1999-01-01), pages 158 - 166, XP011437871 * |
CHIKAYOSHI SUMI ET AL.: "Improvement of the measurement accuracy of the local tissue displacement by making a local region smaller during the iterative rf-echo phase matching", JOURNAL OF MEDICAL ULTRASONICS, vol. 24, no. 3, 1997, pages 346, XP008168138 * |
CHIKAYOSHI SUMI ET AL.: "Nijigen Cross Specter no Iso o Mochiita Hen'i Vector Keisokuho", THE JAPAN SOCIETY OF ULTRASONICS IN MEDICINE, DAI 64 KAI KENKYU HAPPYOKAI KOEN SHOROKUSHU, vol. 21, pages 279, XP008168236 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10314964B2 (en) | 2013-07-25 | 2019-06-11 | Fresenius Medical Care Deutschland Gmbh | Sterile tube covering for a medical tubing system |
CN109765540A (zh) * | 2019-02-26 | 2019-05-17 | 南京莱斯电子设备有限公司 | 一种频率步进体制米波雷达目标提取方法 |
Also Published As
Publication number | Publication date |
---|---|
US20120002506A1 (en) | 2012-01-05 |
JP5692079B2 (ja) | 2015-04-01 |
US8675450B2 (en) | 2014-03-18 |
JPWO2011089898A1 (ja) | 2013-05-23 |
EP2526872A4 (en) | 2014-04-30 |
CN102361597A (zh) | 2012-02-22 |
EP2526872A1 (en) | 2012-11-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5692079B2 (ja) | 変位推定方法、変位推定装置 | |
KR102332464B1 (ko) | 분석 데이터로부터 전단파 추정 | |
US9622711B2 (en) | System and method for measurement of shear wave speed from multi-directional wave fields | |
KR101586998B1 (ko) | 의료용 초음파 이미징에서의 전단파 정보의 해석 | |
JP5539218B2 (ja) | 超音波診断システムおよび超音波診断装置 | |
JP5121823B2 (ja) | 画像データ処理システム | |
KR101983126B1 (ko) | 음향 방사력 이미징에서의 적응식 모션 추정 | |
EP2881041A1 (en) | Apparatus and method for ultrasonic diagnosis | |
US10324063B2 (en) | Methods and systems for measuring properties with ultrasound | |
US9642600B2 (en) | Shear wave attenuation from k-space analysis system | |
US11761928B2 (en) | Method and system for ultrasonic characterization of a medium | |
CN114176639A (zh) | 用于介质的超声表征的方法和系统 | |
KR20220036863A (ko) | 매질의 초음파 특성화 방법 및 시스템 | |
US20210007712A1 (en) | Ultrasonic diagnostic apparatus and method for controlling ultrasonic diagnostic apparatus | |
US20110245676A1 (en) | Method and apparatus for ultrasound signal acquisition and processing | |
GB2436655A (en) | Deformation determination in pulse echo or ultrasonic imaging | |
US20220082693A1 (en) | Method and system for ultrasonic characterization of a medium | |
JP2023540954A (ja) | 媒体の超音波特性評価用の方法及びシステム | |
US9918698B2 (en) | System and method for gradient-based k-space search for shear wave velocity dispersion estimation | |
Rodriguez-Molares et al. | Reconstruction of specular reflectors by iterative image source localization | |
JP2669204B2 (ja) | 探査装置 | |
Nguyen et al. | Improvements to ultrasonic beamformer design and implementation derived from the task-based analytical framework | |
US10928515B2 (en) | Method for generating ultrasound image and computer readable medium | |
WO2023278445A9 (en) | System and method for non-invasive determination of pressure in a biological compartment | |
Cui et al. | A robust phase zero estimator for ultrasonic elastography using quality-guided seeding strategy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 201180001460.0 Country of ref document: CN |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2011533902 Country of ref document: JP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2011734505 Country of ref document: EP |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 11734505 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13256642 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |