WO2017119301A1 - Dispositif de capture d'image échographique - Google Patents
Dispositif de capture d'image échographique Download PDFInfo
- Publication number
- WO2017119301A1 WO2017119301A1 PCT/JP2016/088169 JP2016088169W WO2017119301A1 WO 2017119301 A1 WO2017119301 A1 WO 2017119301A1 JP 2016088169 W JP2016088169 W JP 2016088169W WO 2017119301 A1 WO2017119301 A1 WO 2017119301A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- doppler
- aliasing
- imaging apparatus
- shift amount
- processing unit
- 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.)
- Ceased
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
Definitions
- the present invention relates to an ultrasonic imaging apparatus, and more particularly to an ultrasonic imaging apparatus having a Doppler calculation function based on a vector Doppler method.
- an ultrasonic Doppler method is widely known as a method for measuring the velocity of a tissue including blood flow in a target portion.
- the ultrasonic Doppler method only the velocity component (Doppler velocity) in the transmission direction of the ultrasonic beam can be directly measured.
- a method has been proposed in which a Doppler velocity is obtained from two or more different directions for one point (measurement point) of a target tissue, and a velocity vector is calculated using the Doppler velocities having different angles. (Non-Patent Document 1).
- This method is called a vector Doppler method.
- the velocity vector of the tissue to be measured that is, the tissue velocity vector
- the velocity of the tissue is calculated, and the velocity of the tissue can be measured.
- the ultrasonic Doppler method calculates the Doppler velocity using the fact that a frequency shift is caused by moving or moving away from the ultrasonic blood flow.
- an RF (Radio Frequency) signal is orthogonally detected, and a frequency shift is detected and calculated as a Doppler phase shift amount (hereinafter simply referred to as a Doppler shift amount).
- a Doppler shift amount the Doppler phase shift amount
- aliasing can be easily corrected because it is visualized by an ultrasonic Doppler image. Further, since the Doppler shift amount is obtained as an average value of the values of a large number of elements constituting the ultrasonic probe, aliasing inherent in the signal of each element does not become a problem.
- the Doppler shift amount is calculated using the received signal obtained by measuring echoes from a predetermined measurement point at a plurality of angles, aliasing exists even in one of the received Doppler shift amounts. In this case, the calculated vector accuracy is significantly reduced.
- the solution of this problem is indispensable when performing various clinical applications using the vector Doppler method.
- This invention makes it a subject to solve the said problem in the case of applying a vector Doppler method. Specifically, it is an object to provide a technique for automatically correcting aliasing for Doppler shift amounts measured from a plurality of directions.
- one of the representative ultrasonic imaging apparatuses of the present invention is configured to individually measure a plurality of openings set in an ultrasonic probe and the plurality of openings and a measurement target region. Based on the geometric relationship with the point, aliasing of the Doppler shift amount calculated for each opening is detected. That is, the present invention is based on detecting aliasing using Doppler shift information measured from a plurality of directions with respect to a velocity vector at one measurement point, and the spatial continuity of the tissue itself with respect to the Doppler shift amount in a specific direction.
- the prior art for example, the invention described in Patent Document 1) that detects aliasing by using the method is completely different from the technical idea.
- an ultrasonic imaging apparatus of the present invention transmits an ultrasonic wave to an inspection object and receives an echo signal (received signal) reflected from the inspection object, and an ultrasonic probe. And a signal processing unit for processing a reception signal received at each of the plurality of set openings.
- the signal processing unit calculates tissue Doppler velocity information that is information regarding the velocity of the tissue to be examined for each opening from the received signal, and detects aliasing included in the calculated tissue Doppler velocity information.
- the tissue velocity can be calculated with high accuracy.
- FIG. 1 is an overall block diagram illustrating an embodiment of an ultrasonic imaging apparatus. It is a figure which shows the outline
- FIGS. 4A and 4B are diagrams illustrating linear fitting. FIG.
- FIG. 4A illustrates a case where fitting is performed including data that causes aliasing
- FIG. 5B illustrates a case where fitting is performed using data in the vicinity of a reference value. It is a figure explaining the threshold value for aliasing determination. It is a figure which shows the flow of the data in a signal processing part. It is a flowchart figure of the process of the signal processing part by the modification 1 of 1st embodiment. It is a flowchart figure of the aliasing correction
- FIGS. 4A and 4B are diagrams showing examples of display images, where FIG. 5A shows a case where the effect of the aliasing process is insufficient, and FIG. 6B shows a case where the effect of the process by the aliasing processing unit is obtained.
- the ultrasonic imaging apparatus (ultrasonic transmitting / receiving apparatus) of the present embodiment includes an ultrasonic probe (2) that transmits an ultrasonic wave to the inspection object (3) and receives an echo reflected from the inspection object as a reception signal; And a signal processing unit (15) for processing a reception signal received by the ultrasonic probe.
- the ultrasonic probe is formed with a plurality of openings each composed of one or more elements, and the signal processing unit processes a received signal for each opening.
- the signal processing unit (15) calculates the frequency change amount or phase change amount (referred to as Doppler shift amount) of the transmitted ultrasonic wave from the echo signal, and uses the Doppler shift amount, the sound velocity, and the wavelength to transmit the transmitted ultrasonic wave of the tissue velocity.
- Doppler shift amount the frequency change amount or phase change amount
- a Doppler velocity calculation unit (153) that calculates a Doppler velocity that is a transmission direction component and a velocity vector of a measurement target tissue, and an aliasing process that detects aliasing based on the Doppler shift amount calculated by the Doppler velocity calculation unit Part (155).
- the Doppler velocity is calculated using the sound velocity and the wavelength as the frequency change or phase change amount (Doppler shift amount) of the transmitted ultrasonic wave. Since this is a known technique, details are omitted.
- the aliasing processing unit (155) updates the data set composed of the measured Doppler shift amount calculated from the received signals at the individual openings to a new data set composed of the Doppler shift amount not including aliasing, and performs the Doppler velocity calculation unit (153). ) Calculates the Doppler velocity and the tissue velocity vector using the new data set updated by the aliasing processing unit.
- the processing by the aliasing processor and the calculation of the tissue velocity by the Doppler velocity calculator can take various forms.
- the ultrasonic imaging apparatus of the present embodiment includes an input unit for a person who performs an inspection using the apparatus (hereinafter referred to as an inspector) to input conditions and instructions, a display unit for displaying processing results to the inspector, and the like May be provided.
- FIG. 1 shows an apparatus configuration example of an ultrasonic imaging apparatus to which the present invention is applied.
- the ultrasonic imaging apparatus 100 of the present embodiment includes an apparatus main body 1 and an ultrasonic probe 2.
- the apparatus main body 1 generates an ultrasonic image while controlling the ultrasonic probe 2, and includes an input unit 10, a control unit 11, an ultrasonic signal generator 12, an ultrasonic reception circuit 13, a display unit 14, and A signal processing unit 15 is provided.
- the ultrasonic probe 2 is in contact with the living body (subject) 3 and irradiates the irradiation area 30 in the living body with ultrasonic waves according to the signal generated by the ultrasonic signal generator 12. An echo signal that is a reflected wave of the signal is received.
- the ultrasonic probe 2 has a large number (for example, about 30 to 10,000) of driving elements arranged in a one-dimensional direction or a two-dimensional direction.
- the driving elements hereinafter also simply referred to as elements
- Ultrasonic waves can be transmitted and received, and some or all of them are used depending on the measurement.
- the ultrasonic probe 2 generates a continuous wave or a pulse wave according to the scanning method.
- a two-dimensional imaging method for imaging a two-dimensional cross section or a three-dimensional imaging method for imaging a three-dimensional region can be appropriately selected according to the scanning method of the ultrasonic probe 2.
- the ultrasonic signal generator 12 includes an oscillator that generates a signal having a predetermined frequency, and sends a drive signal to the ultrasonic probe 2.
- a convergent beam, a plane wave for high-speed imaging, a diffused wave, or the like can be transmitted.
- the ultrasonic reception circuit 13 includes a circuit that performs signal processing such as amplification and delay processing on the reflected echo signal received by the ultrasonic probe 2 in addition to the reception circuit. Furthermore, the ultrasonic reception circuit 13 may include a reception data memory 131 that stores a time-series reception signal in each element.
- the reception signal stored in the reception data memory 131 is a time-series RF signal, and is hereinafter also referred to as reception data.
- the memory When storing the received data for each element, the same number of received data memories as the number of elements are required, so the memory becomes enormous.
- the memory may be reduced by bundling the reception data of a plurality of elements by the conventional delay addition method to obtain the reception data of one opening.
- the number of elements constituting one opening and the number of openings may be set in advance, or may be set by the user via the input unit 10 described later.
- the signal processing unit 15 has a function of generating an ultrasound image from an echo signal from the ultrasound probe 2 and includes an envelope detection unit and a unit that performs log compression. Details of the signal processing unit 15 will be described later.
- the apparatus main body 1 includes a scan converter and an A / D (Analog-to-digital) converter.
- the A / D converter is provided before the signal processing unit 15.
- the sampling frequency is usually between 20 MHz and 50 MHz.
- the scan converter may be included in the ultrasonic receiving circuit 13 or may be provided in the subsequent stage of the signal processing unit 15.
- the ultrasonic receiving circuit 13 includes a scan converter, there is an advantage that the amount of data handled by the signal processing unit 15 is reduced.
- the scan converter is not included in the ultrasonic receiving circuit 13, a large amount of data can be handled by the signal processing unit 15, and a highly accurate measuring device can be realized.
- the control unit 11 controls the ultrasonic signal generator 12, the ultrasonic reception circuit 13, the display unit 14, and the signal processing unit 15 based on the operating conditions of the ultrasonic imaging apparatus 100 set by the input unit 10. It can be constructed in a CPU (Central Processing Unit) of a computer system. Note that some or all of the calculations performed by the signal processing unit 15 may be realized on the same CPU.
- CPU Central Processing Unit
- the input unit 10 includes a keyboard and a pointing device for a doctor or engineer (hereinafter collectively referred to as an examiner) who operates the ultrasonic imaging apparatus 100 to set the operating conditions of the ultrasonic imaging apparatus 100 to the control unit 11. Also, when information from an external device such as an electrocardiogram is used for the examination, it also functions as an external signal input unit.
- an examiner a doctor or engineer who operates the ultrasonic imaging apparatus 100 to set the operating conditions of the ultrasonic imaging apparatus 100 to the control unit 11. Also, when information from an external device such as an electrocardiogram is used for the examination, it also functions as an external signal input unit.
- the display unit 14 outputs information obtained by the signal processing unit 15, for example, an ultrasonic image.
- the display unit 14 may constitute a UI (User Interface) for the inspector to set operation conditions and the like together with the input unit 10 and display a necessary GUI (Graphical User Interface).
- UI User Interface
- GUI Graphic User Interface
- the signal processing unit 15 includes a tomographic image forming unit 151, a Doppler speed calculation unit 153, an aliasing processing unit 155, a display image forming unit 156, and a memory 157 as main elements.
- a part or all of the arithmetic functions of each unit included in the signal processing unit 15 may be realized by the CPU configuring the control unit 11 described above, ASIC (Application Specific Integrated Circuit) or FPGA (Field Programmable Gate Array). It is also possible to realize with hardware such as).
- the tomographic image forming unit 151 uses, for example, a B-mode image, that is, a two-dimensional tissue shape image using a planar imaging method of an ultrasonic irradiation target, from a reception signal (echo signal) output from the ultrasonic reception circuit 13; Alternatively, a three-dimensional tissue shape image using a three-dimensional imaging method is formed.
- a B-mode image that is, a two-dimensional tissue shape image using a planar imaging method of an ultrasonic irradiation target, from a reception signal (echo signal) output from the ultrasonic reception circuit 13
- a reception signal echo signal
- the tissue shape image created by the tomographic image forming unit 151 is an image for presenting the imaging target shape to the examiner, and is created using an echo signal acquired by imaging different from Doppler imaging.
- the tomographic image forming unit 151 is not an essential component, but it is more diagnostic by presenting a tissue shape image together with tissue Doppler velocity information that is information relating to the velocity of the tissue. Suitable information can be provided.
- the Doppler velocity calculation unit 153 performs planar imaging on various tissue Doppler velocity information such as color Doppler, tissue Doppler, pulse wave Doppler, and continuous wave Doppler from the echo signal output from the ultrasonic receiving circuit 13. Two-dimensional tissue Doppler velocity information using the method or three-dimensional tissue Doppler velocity information using the stereoscopic imaging method is extracted. Furthermore, in the present embodiment, the Doppler velocity calculation unit 153 performs reception processing for each element stored in the reception data memory 131 or reception data obtained by bundling it for each opening to perform processing based on the vector Doppler method. The tissue Doppler velocity information is calculated for each of a plurality of elements or openings.
- the aliasing processing unit 155 detects aliasing when calculating the Doppler velocity based on the tissue Doppler velocity information calculated for each opening. As described above, when the Doppler velocity calculation unit 153 calculates the Doppler shift amount, which is one of the tissue Doppler velocity information, for each opening, aliasing may exist in the Doppler shift amount due to calculation restrictions. The aliasing processing unit 155 specifies which aperture Doppler shift amount is aliased based on the relationship between the geometric position of the aperture portion and each Doppler shift amount, and performs a subsequent calculation using each Doppler shift amount. For example, in calculating the tissue velocity vector, a process for eliminating the influence of aliasing is performed.
- the aliasing processing unit 155 may reject the Doppler shift amount in which aliasing has occurred, or may perform aliasing correction processing. Further, the Doppler velocity may be obtained again based on the rejected or corrected Doppler shift amount, or the tissue velocity vector information may be calculated. The specific contents of the aliasing detection method and processing will be described in detail later.
- the display image forming unit 156 forms a display image displayed on the display unit 14, and is calculated by the tomographic image formed by the tomographic image forming unit 151, the Doppler waveform obtained by Doppler measurement, and the aliasing processing unit 155.
- a display image is formed in accordance with a predetermined format or an instruction input from the input unit 10 with the various amounts.
- the memory 157 stores an echo signal, information necessary for calculation in the signal processing unit 15 and a processing result of the signal processing unit 15.
- the irradiation region 30 shown in FIG. 1 is a part including the carotid artery will be described as a specific example.
- the irradiation region 30 may be a blood vessel or other heart desired by the examiner. .
- FIG. 2 An outline of the processing of the ultrasonic imaging apparatus 100 of the present embodiment is shown in FIG. As shown in FIG. 2, a process of acquiring and receiving a reception signal for each element by transmitting and receiving ultrasonic waves in each of the plurality of openings (S0), and each of the plurality of openings from the reception signal for each element , Processing for generating reception data composed of RF signals (S1), processing for calculating Doppler shift amount of each opening using reception data of each opening (S2), and using Doppler shift amount of each opening Then, processing for detecting aliasing (S3) is performed.
- the processing S0 is mainly performed by the ultrasonic signal generator 12 and the ultrasonic receiving circuit 13, and the processing S1 to S3 is performed by the signal processing unit 15.
- the ultrasonic signal generator 12 sends a drive signal having a predetermined frequency to the ultrasonic probe 2.
- a convergent beam may be transmitted as in a general-purpose medical ultrasonic diagnostic apparatus, or a plane wave or a diffused wave for high-speed imaging may be transmitted. May be.
- the apodization may be set to a Hanning window or a Hamming window using a delay addition method.
- the Doppler frequency is defined by the angle between the flow direction (tissue velocity (eg, blood flow velocity)) and the transmitted beam at one point (measurement point) in the measurement target tissue (eg, blood vessel), and the flow direction at the measurement point. Since it is determined by the angle with the receiving element direction, the transmission method is not limited as long as two angles can be defined.
- FIG. 3 shows an example in which the transmitted beam is a convergent beam.
- the angle ⁇ between the tissue velocity V (thick arrow) that is the flow direction at one measurement point P in the measurement target tissue and the transmitted beam 301 is measured.
- the relationship between the tissue velocity V which is the flow direction at the point P and the angle (referred to as a prospective angle) ⁇ i in the receiving element direction is shown.
- the prospective angle ⁇ i varies depending on the position of the elements 21 in the arrangement direction.
- each element 21 of the ultrasonic probe 2 After each of the plurality of openings of the ultrasonic probe 2 transmits an ultrasonic wave, each element 21 of the ultrasonic probe 2 receives the ultrasonic wave and detects an ultrasonic signal.
- the ultrasonic signal detected by each element 21 is stored in the reception data memory 131 for each element as an RF signal that is a time-series reception signal, and the ultrasonic reception circuit 13 acquires the reception signal for each element 21.
- morphological imaging such as B mode may be performed separately from imaging of ultrasonic Doppler (S0).
- the morphological imaging is the same as that of a conventional ultrasonic imaging apparatus, and will be briefly described.
- the ultrasonic frequency of the B-mode image is in the range of 1 MHz to 20 MHz that can be imaged.
- the frame rate is set within a range in which the movement of the heart that fluctuates depending on the heartbeat can be captured. Specifically, the frequency is set to 3 Hz or higher so that the time phase of the heart can be observed.
- the tomographic image forming unit 151 forms, for example, a B-mode image from the echo signal output from the ultrasonic receiving circuit 13.
- the ultrasound biological image may be either a two-dimensional image using a planar imaging method or a three-dimensional image using a stereoscopic imaging method, and acquires data in time series.
- a three-dimensional image can be obtained by scanning in the short axis direction using, for example, a two-dimensional array in the short axis direction or a mechanical probe.
- These element groups used for one measurement may all be used as one opening, or as shown in FIG. 4, a plurality of elements (six elements in FIG. 4) are bundled to form one opening. Also good.
- the number N of openings may be at least 3 or more and may be 3 or more and the number of elements.
- the plurality of openings do not necessarily have to be continuous, but are preferably continuous.
- Such information on the number and position of the openings that bundle the elements of the ultrasonic probe 2 is set in advance in the ultrasonic receiving circuit 13 and the signal processing unit 15 (Doppler velocity calculation unit 153). However, the inspector can also set the number and position of the openings as one parameter of the operating condition.
- the Doppler speed calculation unit 153 stores the bundled openings in the reception data memory 131 based on the information about the openings set in this way, and measures the received signal that is the reception data of the bundled reception elements. By performing delay addition according to points, reception data for each opening is created. At this time, a sample gate is set for the reception signal, and reception data composed of an RF signal only in a specific region (referred to as a measurement region) is created.
- the Doppler velocity calculation unit 153 calculates the Doppler shift amount by performing quadrature detection on the reception data including the RF signal obtained in the processing S1 and detecting the phase difference between the transmissions of the ultrasonic signals transmitted.
- a technique for calculating the Doppler shift amount from the reception data subjected to the quadrature detection is known, and detailed description thereof is omitted here, but the outline is as follows.
- the I (t) signal and Q (t) signal obtained by quadrature detection of received data consisting of RF signals are subjected to LPF processing to remove high frequency components, and the phase is extracted by normalizing the amplitude components.
- the phase difference (Doppler shift amount) between transmissions is calculated by the autocorrelation method or the like.
- the phase difference ⁇ is expressed by Expression (1).
- Expression (2) there is a relationship expressed by Expression (2) between the Doppler shift frequency fd and the carrier center frequency f0.
- the Doppler shift amount ( ⁇ ) is calculated using all received signals, and the tissue velocity V is obtained from the equations (1) and (2).
- the vector Doppler method individual elements (opening portions) are obtained.
- the phase difference is obtained every time.
- equation (2) can be expressed by equation (3) when the angle ⁇ of the direction of blood flow at the measurement point P shown in FIG. 3 and the expected angle ⁇ i of each element are used.
- fd i is a Doppler shift frequency calculated from the received signal of the i-th element. Equation (3) is a description for each element, but for the opening part where the elements are bundled, the Doppler shift frequency at each opening part can be calculated by using the angle of the center part of the opening part as ⁇ i.
- the Doppler shift frequency in Expression (3) is actually calculated as a phase difference (Doppler shift amount) taking into account Expression (1). Data including the Doppler shift amount of each opening is obtained for a plurality of measurement points in the measurement region.
- the Doppler speed calculation unit 153 obtains a data set (actual measurement data set 200) that is a collection of data that is actually measured at each opening and includes the calculated Doppler shift amount that is the calculated Doppler shift amount. Then, the Doppler speed calculation unit 153 sends the actual measurement data set 200 to the aliasing processing unit 155.
- the aliasing processing unit 155 receives the measured data set 200 including the measured Doppler shift amount of each opening calculated in the processing S2 (FIG. 6: S30), and performs processing using the fitting and the threshold value using this, and the aliasing is performed. Detect where it is happening and remove the data.
- Fig. 5 shows an image of the measured Doppler shift amount plotted at each opening.
- the horizontal axis of the graph is the expected angle ⁇ of the element, and the vertical axis is the Doppler shift amount.
- the number of openings is set to 8, and the measured Doppler shift amount for each opening is calculated.
- the measured Doppler shift amount changes almost linearly (501), but when aliasing occurs,
- the calculated actual Doppler shift amount exceeds the phase detection range, the polarity of the phase is reversed and aliasing occurs (509).
- the phase detection range differs depending on the calculation method used when calculating the measured Doppler shift amount.
- FIG. 5 shows the case of ⁇ ⁇ / 2, but it may be ⁇ ⁇ .
- the measured Doppler shift amount (fd) is calculated by using all the received signals. Therefore, even if there is a phase wrapping for each element or each of the plurality of openings, When the measured Doppler shift amount is obtained, individual phase wrap-around becomes obvious. For this reason, in this process S3, the process which excludes the data which has caused the aliasing is performed.
- Process S3 includes a step (S31) of fitting an actual measured Doppler shift amount of each opening with a predetermined opening among the plurality of openings as a reference point, a step (S32) of setting a threshold, an actual Doppler shift amount and a fitting function, It is determined whether the difference is within the threshold and updating the fitting function (S33), and calculating the tissue velocity vector using the finally obtained fitting function (S34).
- the aliasing processing unit 155 performs S31 to S33
- the Doppler speed calculation unit 153 performs Step S34.
- the aliasing processing unit 155 receives the actual measurement data set 200 including the actual measured Doppler shift amount for each opening obtained by the process S2 from the Doppler velocity calculation unit 153.
- ⁇ Step S31 The aliasing processing unit 155 sets a reference and performs a fitting process (for example, a linear fitting process) based on the measured data set 200 including the measured Doppler shift amount for each opening obtained in the process S2. That is, the aliasing processing unit 155 fits the measured Doppler shift amount of each opening (determines a fitting function of the measured Doppler shift amount). If the fitting function of the measured Doppler shift amount for each opening is ⁇ i, the fitting function can be described as in Expression (4).
- This equation corresponds to the straight line 501 in FIG. 5A, where ⁇ is the slope of the straight line and ⁇ CFM is the intercept value.
- the phase difference is a relative value, and the intercept value (position on the vertical axis in FIG. 5) may change.
- a predetermined opening is determined as a reference point (reference opening), and ⁇ CFM is a reference point.
- the position of the opening as the reference point is not particularly limited, but it is customary to set the color Doppler so that it does not fold when the inspector sets the color Doppler. It is premised that no wrapping has occurred.
- ⁇ can be described as follows if information on all openings is used on average.
- ⁇ may be calculated by weighting according to the angle of ⁇ as in equation (7).
- FIG. 7A shows a case where all data including aliasing is used, and the fitting is remarkably displaced.
- FIG. 7B when the vicinity of the reference point is used, good fitting is possible.
- the aliasing processing unit 155 calculates a threshold for determining whether or not the measured Doppler shift amount causes aliasing, and sets it in the ultrasonic imaging apparatus 100.
- the threshold value for example, using a variance ⁇ of a difference between the fitting function and the measured Doppler shift amount, any value about 0.5 to 4 times the variance value may be used.
- the threshold may be preset on the ultrasonic imaging apparatus 100 side in advance, or may be configured to be set by the inspector via the input unit 10.
- the aliasing processing unit 155 has, at each measurement point, a fitting function ( ⁇ i) (equation (4)) (fitting function ( ⁇ i ) Is calculated (S331), and it is determined whether or not the difference is equal to or smaller than the threshold set in step S32 (S332). When the threshold value is exceeded, the aliasing processing unit 155 determines that there is a high possibility that aliasing has occurred, rejects the data, and does not use it for subsequent calculations (S333).
- the aliasing processing unit 155 uses data that is equal to or less than the threshold value (S334), and sets the measured Doppler shift amount of the opening portion that is less than or equal to the threshold value as a new data set 210 as shown in FIG. It passes to the calculating part 153 (S335).
- Figure 8 shows the relationship between the measured Doppler shift amount and the threshold value.
- the threshold values 510 and 520 are set above and below the fitting function 500.
- the aliasing processing unit 155 determines whether or not the measured Doppler shift amount (indicated by ⁇ ) falls within a range between upper and lower dotted lines that are threshold values.
- the Doppler velocity calculation unit 153 uses the inclination of the determined fitting function and the measured Doppler shift amount of the reference point to obtain the gradient between the tissue velocity V at the measurement point and the transmission beam (transmitting beam) according to the following equation (8).
- (Angle) ⁇ is calculated. That is, if the fitting function is determined and the slope ⁇ is obtained, ⁇ can be calculated.
- the tissue velocity V can be calculated from the above formulas (1) and (2) using ⁇ CFM .
- the Doppler velocity calculation unit 153 calculates the tissue velocity vector using the fitting function.
- step S33 may be completed by a single process, but may be repeated twice or more as shown by the dotted arrow in FIG. 6 to update the fitting function and a new data set. .
- the aliasing processing unit 155 returns to step S31 and repeats the processes from step S31 to S33 M times (S36). Accordingly, it is possible to improve the accuracy of ⁇ calculated in S34.
- the repetition of the step ends when, for example, all the data that remains without being rejected falls below a predetermined threshold, or ends when the predetermined number of times is reached.
- the threshold value set in step S32 may be constant, or the threshold value may be changed for each repetition. When changing, for example, the threshold value is gradually changed to a smaller value.
- the aliasing processing unit 155 and the Doppler velocity calculation unit 153 perform the process S3 illustrated in FIG. 6 for all measurement points included in the measurement region. Thereby, a velocity vector (a velocity value and a direction) is calculated for a tissue including blood flow included in the measurement region.
- the calculated velocity vector is stored in the memory 157 and then displayed on the display unit 14 as a display image in which an arrow indicating the vector is superimposed on a morphological image obtained separately or as a numerical value by the display image creation unit 156. be able to.
- the display mode will be described later.
- the Doppler calculation unit 153 calculates the measured Doppler shift amount for each opening using the reception data for each opening, obtains the measured data set 200 including the measured Doppler shift amount for each opening, and performs aliasing.
- the data is sent to the processing unit 155.
- the aliasing processing unit 155 detects aliasing present in the measured Doppler shift amount for each opening included in the measured data set 200, removes (rejects) the Doppler shift amount in which aliasing has been detected from the measured data set 200, and leaves the remaining openings.
- a new data set 210 consisting of the amount of Doppler shift is obtained and sent to the Doppler calculation unit 153.
- the Doppler calculation unit 153 calculates the velocity vector 250 of the measurement target region using the Doppler shift amount of the remaining opening that has not been removed by the aliasing processing unit 155 included in the new data set 210.
- the display image forming unit 156 forms a display image of a new data set 210 or speed vector 250 displayed on the display unit 14.
- the aliased data is removed from the Doppler shift amount for each opening. As a result, the accuracy of velocity vector calculation can be increased.
- step S34 This modification is different from the first embodiment in the processing (step S34) of the Doppler velocity calculation unit 153 after removing data with aliasing.
- the processing flow of this modification is shown in FIG. In FIG. 10, the same steps as those in FIG.
- the details of step S33 in FIG. 10 are the same as S331 to S334 in S33 shown in FIG. Alternatively, it may be the same as S330 (excluding S335) in FIG.
- the angle ⁇ is calculated using the slope ⁇ of the updated fitting function (Equation (8)).
- the Doppler speed is calculated in the same manner as in the conventional color Doppler method (S35). For this reason, first, the Doppler shift amount for each opening is averaged (S351). At that time, a weighted average may be performed as necessary, thereby further improving the accuracy of the color Doppler.
- the weighting for example, a weight inversely proportional to the size of the expected angle of the opening can be considered.
- the Doppler speed is calculated by Equations (1) and (2) (S352).
- the tissue velocity vector can also be calculated using the Doppler velocity calculated as described above. Since this is described in a known technique (Non-Patent Document 1), it can be referred to.
- the calculation of the Doppler speed is the same as in the conventional method, but the data after the average addition / weighting addition used for the calculation of the Doppler speed does not include the data with aliasing.
- the accuracy of color Doppler can be improved compared with the conventional method which does not remove.
- a step of automatically adjusting the speed range is added after the fitting function is obtained.
- This step can be inserted, for example, between steps S33 and S34 in FIG. 6 or after S34, or between steps S33 and S35 in FIG. 10 or after S35.
- the Doppler shift amount assumed in all openings that is, the range of the Doppler shift amount can be estimated.
- the speed range is automatically adjusted so that the Doppler shift amount of all the openings is within the threshold value.
- the phase shift range is inversely proportional to the speed range. Therefore, for example, when the region determined by the fitting function and the threshold (the region sandwiched between the two dotted lines 510 and 520 in FIG. 8) is narrower than ⁇ ⁇ / 2 in the total width of ⁇ , the entire range of ⁇ ⁇ / 2 is obtained. Narrow the speed range to spread out. As a result, the blood flow velocity distributed within the narrowed velocity range can be depicted with high resolution.
- the speed range is set so that the Doppler shift amount of all the openings is within the threshold. To spread. As a result, a relatively slow blood flow to a fast blood flow can be depicted in a wide measurement region without phase wrapping.
- the calculation of the speed range adjustment amount is executed by adding a speed range correction unit as a function of the signal processing unit 15, for example, the aliasing processing unit 153 or the Doppler calculation unit 153.
- the speed range is automatically adjusted based on the determined fitting function, so that the speed range according to the blood flow velocity distribution in the measurement region can be adjusted without adjusting the speed range manually as in the past. Blood flow velocity can be depicted. Note that changing the speed range (changing the speed) is not always recommended because there is a trade-off that information in the low speed range is lost. Therefore, instead of automatically adjusting, a means or GUI for manual adjustment by the inspector may be provided.
- the aliasing processing unit corrects the data (measured Doppler shift amount) determined to be aliasing by the threshold value instead of simply rejecting it.
- Other configurations are the same as those of the first embodiment, and hereinafter, the ultrasonic imaging apparatus of this embodiment will be described focusing on differences from the ultrasonic imaging apparatus 100 of the first embodiment.
- the ultrasonic imaging apparatus of the present embodiment also performs the processes S0 to S3 having the contents shown in FIG. 2 as in the first embodiment, but the contents of the aliasing detection process S3 are changed. That is, in the first embodiment, in step S33 (FIG. 6) of process S3, the measured Doppler shift amount outside the range between the upper and lower thresholds of the fitting function is detected and removed as aliasing, but the second In the embodiment, the aliasing correction is performed on the removed one.
- the aliasing processing unit 155 provides a reference for the measured data set 200 including the measured Doppler shift amount of each opening for each measurement point, performs fitting (S31), and calculates and sets a threshold value. (S32).
- the aliasing processing unit 155 calculates a difference between the measured Doppler shift amount of each opening and the fitting function (S331), and determines whether the difference is equal to or less than a threshold value (S332). Data equal to or less than the threshold value is used for the calculation of the velocity vector as it is as in the first embodiment (S334, S335).
- the fitting function is updated by performing iterative calculation, it is added to the data set for iterative calculation (S334, S335).
- step S331 determines that aliasing has occurred and performs aliasing.
- a correction process is performed (S336). An outline of the aliasing correction will be described with reference to FIG.
- the aliasing processing unit 155 further determines whether or not the corrected Doppler shift data group 303 is equal to or smaller than the threshold (S337). If the corrected Doppler shift data group 303 is within the threshold region, the data is used (S338). If it is outside the threshold region, it is rejected (S339). As a result, a new data set 210 is created (S335).
- the threshold value used in step S337 may be the same as or different from the threshold value used in step S332, and is set in step S32.
- the aliasing processing unit 155 passes the new data set 210 to the Doppler velocity calculation unit 153, and the Doppler velocity calculation unit 153 calculates a tissue velocity vector (S34). Alternatively, the aliasing processing unit 155 repeats the processing from step S31 to S330 M times (S36), improves the accuracy of the data set, and then passes it to the Doppler speed calculation unit 153, as in the first embodiment. . Further, the Doppler velocity calculation unit 153 may calculate the velocity vector by calculating the angle ⁇ in the velocity direction by the calculation of the vector Doppler method, or the conventional color Doppler as in the first modification of the first embodiment. The velocity vector may be calculated by obtaining an average value of the Doppler shift amount of each opening based on the method.
- the present embodiment it is possible to acquire a more accurate data set by performing a series of processing one or more iterations.
- the threshold value may be decreased stepwise.
- the aliasing processing unit corrects the measured Doppler shift amount determined to be aliasing out of the data set including the measured Doppler shift amount, and creates a new one. Get the dataset.
- the aliasing processing unit adds the corrected Doppler shift amount to the new data set and sets the difference outside the threshold value. You may add the processing to reject.
- the Doppler speed calculation can be performed without reducing the number of measured data.
- the ultrasonic imaging apparatus of the present embodiment is characterized by including means for enabling an inspector to select processing contents and conditions by the aliasing processing unit. That is, the ultrasonic imaging apparatus of the present embodiment includes an input unit that receives selection of a processing mode by the aliasing processing unit. In addition, a display unit that displays the processing result of the aliasing processing unit or the calculation result of the Doppler velocity calculation unit reflecting the processing result is provided.
- the display image forming unit of the signal processing unit creates a UI that accepts selection by the examiner and displays it on the display unit.
- Other configurations are the same as those of the first and second embodiments, and the present embodiment will be described by appropriately using the elements in the drawings used in these embodiments.
- FIG. 13 shows a functional block diagram of the signal processing unit 15 of the present embodiment.
- the Doppler velocity calculation unit 153 the aliasing processing unit 155, and the display image creation unit 156 among the components constituting the signal processing unit 15 are illustrated, and other elements are omitted, but FIG. Each part shown and other elements may be included.
- the Doppler speed calculation unit 153 includes a vector Doppler calculation unit 1533 and a color Doppler calculation unit 1535.
- the vector Doppler calculation unit 1533 calculates the velocity vector using the measured Doppler shift amount for each opening from which the influence of aliasing is removed by the aliasing processing unit 155.
- the color Doppler calculation unit 1535 calculates the Doppler speed by using the Doppler shift amount obtained by averaging or weighted addition of the measured Doppler shift amount for each aperture from which the influence of aliasing has been removed.
- the Doppler speed calculation unit 153 may further include a speed range setting unit 1531. A speed range used for speed calculation and display is set in the speed range setting unit 1531. The speed range is set via the UI unit 140. Alternatively, when the signal processing unit 15 includes the speed range automatic adjustment function described in the second modification of the first embodiment, the function is set by this function.
- the aliasing processing unit 155 includes an opening setting unit 1553, a first processing unit 1551, and a second processing unit 1552.
- the opening setting unit 1553 holds information on a preset or received through the UI unit 140.
- the first processing unit 1551 creates a new data set by removing the Doppler shift amount causing aliasing from the measured Doppler shift amount of each opening.
- the second processing unit 1552 creates a new data set by performing aliasing correction processing of the Doppler shift amount in which aliasing has occurred among the measured Doppler shift amounts of the respective openings.
- the UI unit 140 is a user interface device including the input unit 10 and the display unit 14 and operates interactively with the inspector.
- An example of a display screen in the UI unit 140 is shown in FIG.
- the display screen 400 includes a result display block (area) 410 that displays the calculation result of the Doppler velocity calculation unit 153, a menu selection block 420 that displays a menu selected by the examiner, and a menu selection.
- An operation block 430 for displaying objects such as buttons and knobs for setting or adjusting parameters in accordance with the imaging method (imaging technique) and operation mode (processing mode) selected via the block 420 is included.
- a menu screen 450 as shown in FIG. 14B is displayed.
- a GUI for receiving selection of an imaging method, selection of a processing mode in the vector Doppler method, setting of parameters (processing conditions) in the vector Doppler method, and the like is displayed.
- selection of the imaging method for example, selection of color Doppler method or vector Doppler method is accepted.
- the selection of the imaging method may be alternative or combined.
- the Doppler velocity calculation unit 153 includes a new data set (each of which is updated by the aliasing processing unit 155 by the vector Doppler unit 1533).
- a calculation using equations (7) and (8) is performed using a measured data set of the Doppler shift amount of the opening to calculate a velocity vector.
- the color Doppler unit 1535 performs a normal color Doppler method calculation using the new data set updated by the aliasing processing unit 155.
- the color Doppler unit 1535 When the combined use with the color Doppler method is selected, the color Doppler unit 1535 performs the color Doppler method together with the vector Doppler method performed by the vector Doppler unit 1533.
- the blood flow vector directly obtained by the color Doppler method is a component in the direction parallel to the transmission beam.
- the color Doppler unit 1535 is a component in the orthogonal direction by an estimation method using a law of conservation of mass. In this case, a blood flow vector in the beam plane is calculated.
- a technique for calculating a blood flow vector in the color Doppler method is known, and a description thereof is omitted here.
- a GUI for setting conditions such as the number of openings and a threshold value may be activated.
- the value may be displayed.
- the buttons and knobs displayed on the operation block 430 are switched to a function of adjusting, for example, a threshold value and a speed range (width and upper limit value or lower limit value).
- the opening setting unit 1553 and the speed range setting unit 1531 accept the operation of the operation block 430 and set the number of openings and the speed range, respectively.
- the first processing unit 1551 and the second processing unit 1552 perform the respective processes (the processes in FIGS. 6 and 10) using the set threshold value.
- a UI for selecting a processing mode is displayed on another menu screen or a pull-down menu screen as illustrated.
- a processing mode for example, according to the first embodiment, a mode for removing data having a large difference from the fitting result (referred to as a normal mode), or a mode for correcting data having a large difference from the fitting result according to the second embodiment One of (correction mode) can be selected.
- These processing modes are determined by a processing program executed by the aliasing unit 155. Either processing mode may be preset in the apparatus, or the inspector can select any one.
- the first processing unit 1551 performs aliasing detection / removal and velocity vector calculation according to the flow of FIG.
- the second processing unit 1552 performs aliasing detection / correction and velocity vector calculation according to the flow of FIG.
- the speed vector calculation method may be determined depending on whether the color Doppler method or the vector Doppler method is selected when the imaging method is selected.
- the ultrasonic imaging apparatus of the present invention can provide various diagnostic information reflecting the blood flow velocity vector.
- the method of providing diagnostic information is not particularly limited, but a typical providing method is a method of displaying as a display image on the display unit 14 of the ultrasonic imaging apparatus.
- the display image is generated by the display image generation unit 156.
- the diagnostic information image generated by the display image generation unit 156 is displayed on, for example, the result display block 410 of the screen 400 shown in FIG.
- a numerical value block 415 may be provided in the result display block 410 to display the blood flow velocity (average value) and the like.
- the display method is not limited to these, and various combinations and omission of elements that are not essential are possible.
- FIG. 15 shows a state where a color Doppler is displayed and a vector is superimposed on the flow of the blood vessel.
- FIG. 15A the vector direction and the blood vessel running direction are deviated, and it is highly possible that the vector calculation accuracy in aliasing is reduced. In such a case, it is expected that correction is performed as shown in FIG. 15B by using the aliasing correction mode.
- the corrected information may be presented to the inspector.
- a mode for automatically increasing the speed range may be prepared on the apparatus side, or the inspector may be prompted to change the speed range.
- the examiner looks at this display and determines that the aliasing processing has not been performed properly if the blood flow velocity vector or blood flow velocity is unnatural or significantly different from the predicted result.
- the processing mode may be changed, or parameters used for the aliasing processing, such as the number of openings, the threshold value, and the speed range may be reset.
- the ultrasonic Doppler measurement suitable for the measurement target such as changing the processing mode when it is determined that the processing needs to be re-executed by looking at the result obtained in the preset normal mode processing. Can be implemented.
- the ultrasonic imaging apparatus of the present invention calculates the Doppler velocity of the tissue using the received signal received at each opening of the ultrasonic probe.
- an aliasing processing unit for detecting aliasing at the time of calculating the tissue velocity vector at the time of calculating the tissue velocity vector (including the case of calculating based on the Doppler velocity).
- the tissue velocity can be calculated with high accuracy.
- the present invention since the present invention does not assume spatial continuity, the present invention can also be applied to the speed of one measurement point.
- the ultrasonic imaging apparatus of the present invention is not limited to the above-described embodiments, and elements can be added or deleted as appropriate.
- the present invention is not limited to the above-described embodiments, and includes various modifications.
- the above-described embodiments have been described in detail for easy understanding of the present invention, and are not necessarily limited to those having all the configurations described.
- a part of the configuration of one embodiment can be replaced with the configuration of another embodiment, and the configuration of another embodiment can be added to the configuration of one embodiment.
- Each of the above-described configurations, functions, and the like may be realized by software by interpreting and executing a program that realizes each function by the processor.
- Information such as a program for realizing each function can be stored in a memory, a hard disk, a recording device such as an SSD (Solid State Drive), or a recording medium such as an IC (Integrated Circuit) card, an SD card, or a DVD.
- Various information does not depend on the data structure. Further, although the description has been given of recording each information in the memory by the expression “store”, it may be expressed as “register” or “set”.
- control lines and information lines indicate what is considered necessary for the explanation, and not all the control lines and information lines on the product are necessarily shown. Actually, it may be considered that almost all the components are connected to each other.
- the calculation accuracy of the tissue velocity vector can be improved. This can contribute to more reliable diagnosis.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Hematology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
L'objectif de la présente invention est, en mesure Doppler vectorielle, d'identifier automatiquement un repliement par rapport à des décalages Doppler qui sont mesurés depuis une pluralité de directions. L'invention concerne un dispositif de capture d'image échographique, comprenant une unité de traitement de signal qui traite des signaux d'écho qui sont reçus avec une pluralité de parties d'ouverture qui sont formées sur une sonde échographique, l'unité de traitement de signal comprenant en outre une unité de calcul de vitesse Doppler qui calcule la vitesse Doppler d'un tissu à partir des signaux reçus qui sont reçus avec chacune des parties d'ouverture, et une unité de traitement de repliement qui détecte un repliement (enroulement de phase) qui survient avec le calcul de vitesse Doppler. L'unité de traitement de repliement détecte le repliement par ajustement du déphasage Doppler que l'unité de calcul de vitesse Doppler a calculé pour chacune des ouvertures, et exécute un traitement tel que l'élimination ou la correction du repliement.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017560100A JP6515205B2 (ja) | 2016-01-04 | 2016-12-21 | 超音波撮像装置 |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016-000045 | 2016-01-04 | ||
| JP2016000045 | 2016-01-04 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2017119301A1 true WO2017119301A1 (fr) | 2017-07-13 |
Family
ID=59274327
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/JP2016/088169 Ceased WO2017119301A1 (fr) | 2016-01-04 | 2016-12-21 | Dispositif de capture d'image échographique |
Country Status (2)
| Country | Link |
|---|---|
| JP (1) | JP6515205B2 (fr) |
| WO (1) | WO2017119301A1 (fr) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20200029210A (ko) * | 2018-09-10 | 2020-03-18 | 국방과학연구소 | 초음파 카메라의 영상 보정 방법 및 장치 |
| CN113116380A (zh) * | 2019-12-31 | 2021-07-16 | 深圳迈瑞生物医疗电子股份有限公司 | 血流成像的混叠指数、混叠区域的显示方法及显示装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013116372A (ja) * | 2007-03-29 | 2013-06-13 | Hitachi Aloka Medical Ltd | 超音波画像処理方法および超音波画像処理装置 |
| JP2014534852A (ja) * | 2011-10-19 | 2014-12-25 | ヴェラゾニックス,インコーポレーテッド | 平面波送信を使用するベクトルドップラーイメージングのための推定及び表示 |
| JP2015500062A (ja) * | 2011-12-01 | 2015-01-05 | マウイ イマギング,インコーポレーテッド | ピングベース及び多数開口ドップラー超音波を用いた運動の検出 |
-
2016
- 2016-12-21 JP JP2017560100A patent/JP6515205B2/ja active Active
- 2016-12-21 WO PCT/JP2016/088169 patent/WO2017119301A1/fr not_active Ceased
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013116372A (ja) * | 2007-03-29 | 2013-06-13 | Hitachi Aloka Medical Ltd | 超音波画像処理方法および超音波画像処理装置 |
| JP2014534852A (ja) * | 2011-10-19 | 2014-12-25 | ヴェラゾニックス,インコーポレーテッド | 平面波送信を使用するベクトルドップラーイメージングのための推定及び表示 |
| JP2015500062A (ja) * | 2011-12-01 | 2015-01-05 | マウイ イマギング,インコーポレーテッド | ピングベース及び多数開口ドップラー超音波を用いた運動の検出 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20200029210A (ko) * | 2018-09-10 | 2020-03-18 | 국방과학연구소 | 초음파 카메라의 영상 보정 방법 및 장치 |
| KR102093897B1 (ko) * | 2018-09-10 | 2020-03-26 | 국방과학연구소 | 초음파 카메라의 영상 보정 방법 및 장치 |
| CN113116380A (zh) * | 2019-12-31 | 2021-07-16 | 深圳迈瑞生物医疗电子股份有限公司 | 血流成像的混叠指数、混叠区域的显示方法及显示装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP6515205B2 (ja) | 2019-05-15 |
| JPWO2017119301A1 (ja) | 2018-06-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5735718B2 (ja) | 超音波診断装置、及び弾性評価方法 | |
| JP3746115B2 (ja) | 超音波診断装置 | |
| JP4787358B2 (ja) | 超音波診断装置 | |
| JP6734079B2 (ja) | 医用診断装置、および医用解析プログラム | |
| US10959704B2 (en) | Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing method | |
| JP6058295B2 (ja) | 超音波診断装置、医用画像処理装置、医用画像処理方法、および医用画像処理プログラム | |
| JP6987496B2 (ja) | 解析装置 | |
| EP2939598A1 (fr) | Dispositif de diagnostic ultrasonore et procédé d'évaluation d'élasticité | |
| US10101450B2 (en) | Medical image processing apparatus, a medical image processing method and a medical diagnosis apparatus | |
| WO2013157553A1 (fr) | Dispositif et procédé de capture d'image ultrasonore | |
| JP2008136860A (ja) | 超音波診断装置およびその画像処理プログラム | |
| US11766245B2 (en) | Ultrasonic diagnostic apparatus and control method | |
| US11039777B2 (en) | Ultrasonic diagnostic apparatus and control method | |
| JP6253360B2 (ja) | 被検体情報取得装置、被検体情報取得方法、及びプログラム | |
| JP6698511B2 (ja) | 超音波信号処理装置、および、それを用いた超音波送受信装置 | |
| JP6515205B2 (ja) | 超音波撮像装置 | |
| JP2020168232A (ja) | 超音波診断装置、信号処理装置、および、プログラム | |
| JP2018186923A (ja) | 超音波撮像装置及び血流速度算出方法 | |
| JP2021049129A (ja) | 超音波診断装置、及び画像処理装置 | |
| US12433565B2 (en) | Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus | |
| JP6677975B2 (ja) | 医用画像診断装置、及び医用画像処理装置 | |
| JP6444519B2 (ja) | 超音波診断装置、及び超音波撮像方法 | |
| JP2021500168A (ja) | 血管狭窄の超音波測定 | |
| JP2004357876A (ja) | 超音波診断装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 16883800 Country of ref document: EP Kind code of ref document: A1 |
|
| ENP | Entry into the national phase |
Ref document number: 2017560100 Country of ref document: JP Kind code of ref document: A |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 16883800 Country of ref document: EP Kind code of ref document: A1 |