[go: up one dir, main page]

WO2008010366A1 - Ultrasonographic device - Google Patents

Ultrasonographic device Download PDF

Info

Publication number
WO2008010366A1
WO2008010366A1 PCT/JP2007/062021 JP2007062021W WO2008010366A1 WO 2008010366 A1 WO2008010366 A1 WO 2008010366A1 JP 2007062021 W JP2007062021 W JP 2007062021W WO 2008010366 A1 WO2008010366 A1 WO 2008010366A1
Authority
WO
WIPO (PCT)
Prior art keywords
interpolation
data
ultrasonic
imaging apparatus
depth
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
Application number
PCT/JP2007/062021
Other languages
English (en)
French (fr)
Inventor
Takashi Azuma
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Ltd
Hitachi Medical Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd, Hitachi Medical Corp filed Critical Hitachi Ltd
Priority to US12/374,065 priority Critical patent/US8506483B2/en
Priority to JP2008525809A priority patent/JP4837037B2/ja
Publication of WO2008010366A1 publication Critical patent/WO2008010366A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52044Scan converters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52077Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging with means for elimination of unwanted signals, e.g. noise or interference
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52079Constructional features
    • G01S7/52084Constructional features related to particular user interfaces

Definitions

  • the present invention relates to an apparatus for displaying an ultrasonic tomographic image.
  • a method of performing sector-shaped beam scanning called sector scanning or convex scanning is widely used in order to ensure a narrow, wide, wide, and visual field.
  • a force with a wider field of view is more advantageous for finding lesions.
  • a narrower aperture may be advantageous as follows.
  • the surface area of the probe is limited when the area where the probe can be brought into contact with the subject is limited, such as when the probe is pressed against the gap between the ribs, such as when imaging the heart. It may be advantageous to have a small.
  • the surface area of the probe is reduced, the number of scan lines decreases, and if the scan conversion to a video display device is performed under the condition of a small number of scan lines, image degradation occurs.
  • Patent Document 1 Japanese Patent Laid-Open No. 11-9603
  • An object of the present invention is to provide an ultrasonic diagnostic apparatus to which an optimal interpolation process can be applied in order to contribute to the removal of noise mixed on a scanning line.
  • the dynamic range of the signal is measured for the acquired signal, and whether to use linear interpolation force or sine function interpolation adaptively for the signal, or weighted average of both I do.
  • the depth position for switching between linear interpolation and sine function interpolation may be set in the device in advance, and the two interpolation methods may be switched on the basis of the interpolation method switching depth. Interpolation values obtained by weighting and adding the interpolation values of the two interpolation methods may be used instead of switching the two interpolation methods completely at the boundary of the interpolation method switching depth.
  • FIG. 1 is a functional block diagram showing an example of an ultrasonic diagnostic apparatus according to the present invention.
  • FIG. 3 is an explanatory diagram of linear interpolation.
  • FIG. 5 is an explanatory diagram regarding the beam width.
  • FIG. 6 is an example of experimental data showing the interpolation result.
  • FIG. 7 is an explanatory diagram showing the relationship between the signal dynamic range and the beam width.
  • FIG. 8 is an explanatory diagram showing the relationship between the signal dynamic range and the beam width.
  • FIG. 9 is a diagram for explaining that the influence of noise differs depending on the interpolation method.
  • FIG. 10 is an explanatory diagram showing how to switch between a signal dynamic range and an interpolation method.
  • FIG. 11 An illustration of the different attenuation factors depending on the mode and the corresponding interpolation point switching position.
  • FIG. 12 is a diagram for explaining how to change the weighting function for switching the interpolation method.
  • FIG. 13 is a schematic diagram of a switching position control input unit.
  • FIG. 14 is a schematic diagram of a switching position control input unit.
  • FIG. 15 is a functional block diagram showing an example of an ultrasonic diagnostic apparatus according to the present invention.
  • FIG. 16 is a functional block diagram showing an example of an ultrasonic diagnostic apparatus according to the present invention.
  • FIG. 17 is an explanatory diagram of a noise removal filter.
  • FIG. 18 is a control flow diagram of the noise removal filter.
  • FIG. 19 is a diagram for explaining a weighting function of a noise removal filter.
  • FIG. 1 is a functional block diagram showing an example of an ultrasonic diagnostic apparatus according to the present invention.
  • a transmission electric pulse is transmitted from the transmission beamformer 3 to the ultrasonic probe 1 installed on the surface of the subject via the transmission / reception switching switch 2 under the control of the control unit 4.
  • the transmission beamformer controls the delay time between each channel of the probe 1 to be in an appropriate state so that the ultrasonic beam travels on a desired scanning line.
  • This transmit beamformer The electrical signal from 3 is converted into an ultrasonic signal by the ultrasonic probe 1 and an ultrasonic pulse is transmitted into the subject.
  • a part of the ultrasonic pulse scattered in the subject is received again by the ultrasonic probe 1 as an echo signal and converted into an electric signal.
  • This received signal is sent to the receiving beamformer 20 via the transmission / reception switching switch 2, and the echo signal from the desired depth on the desired scanning line is selectively enhanced as the memory.
  • the data interpolation unit 30 interpolates the data between the scanning lines that actually obtained the data so as to increase the number of scanning lines. As will be described later, the most suitable interpolation method is selected from a plurality of methods.
  • the interpolated data is sent to the scan con- verter 23, and scan convolution is performed.
  • the data after the scan conversion is sent to the display unit 24 and displayed as an ultrasonic tomographic image.
  • the scanning lines are arranged in the horizontal direction, the frame rate is lowered. Therefore, the scanning lines are usually arranged in the horizontal direction at intervals of one-fourth of the half power of the wavelength at the center frequency.
  • the horizontal sampling interval is, for example, When the center frequency is 3MHz, it is about 250 / zm to 125 / zm.
  • the depth direction sufficient sampling is performed with respect to the carrier frequency.
  • FIG. 3 is an explanatory diagram of linear interpolation for obtaining post-interpolation data from the interpolation of the luminance of the two closest points across the interpolated data
  • Fig. 4 is a sine function that uses data of four or more points. It is explanatory drawing of interpolation. Sine function interpolation considers data points as representative points in a space of finite width, regards this finite width as a rectangular function, and determines the interpolation coefficient using the sine function, which is the Fourier transform of the rectangular function.
  • the coefficient 2 ⁇ from the sine function (solid line and dotted line in Fig. 4) having a peak at the adjacent pixel is added to the luminance value at the adjacent pixel.
  • the sum of the product multiplied by the luminance value at the pixel outside it and the coefficient from the sine function (dashed line and dashed line in Fig. 4) with a peak at the pixel position – 3 ⁇ 2 ⁇ is the interpolated value. .
  • Taylor expansion up to the third order of the sine function is Sin (a) / a + (x— a) (cos (a) / a— sin (a) / a) + (x— a) 2 / 2x (— sin (a) / a— 2cos (a) / a 2 + 2cos (a) / a) + (x— a) 3 / 6x (— cos (a) / a + 3sin (a) / a 2 + 4cos (a) / a 3 — 2sin (a) / a 3 — 6cos (a) / a 4 )
  • An ultrasound image is a convolution of a scatterer distribution of a subject and a point response function specific to the imaging conditions of the apparatus determined by beam forming and post-processing. Beam forming is performed for both transmission and reception, but in the current ultrasound imaging, focusing on one point on one scan line is given priority to not reducing the frame rate for transmission. And it is off-focus elsewhere.
  • dynamic focus in which the focus position is switched continuously according to the receiving time, and the on-focus is uniform in the depth direction.
  • the spread of the point response function in the deep part of the living body where cross flow becomes a problem is determined by the received beam width.
  • the beam width BW at the focal position is determined by the diffraction effect of the ultrasonic wave, as shown in Fig. 5, rather than the geometric beam width.
  • the diffraction angle ⁇ can be approximated by the following equation. Where ⁇ is the wavelength and D is the aperture width.
  • the center frequency in the deep part is 2 MHz and the transmission aperture weight is approximated as half the aperture width, the aperture width is 12.5 mm, and the distance from the probe surface to the imaging site is In the case of Z, the diffraction angle is 0.06 rad, so the beam width is 0.06 X Zmm.
  • the width of the scanning line is typically about 0. Ol x Zmm. In other words, since the sampling is about 6 times more powerful than the beam width, linear interpolation is usually considered sufficient.
  • Fig. 6 shows the experimental results.
  • the vertical axis indicates the luminance
  • the horizontal axis indicates the position in the horizontal direction.
  • Fig. 6 (a) relates to the echo signal from the point reflector, and the black circles connected by the solid line are double-density scans and dots.
  • the white circles connected by lines are the result of double-density by linear interpolation after data acquisition by normal scanning. All results overlap well with almost no difference. In this case, there is no image degradation due to linear interpolation. As already mentioned, this is as expected when the beam width is sufficiently wider than the scanning line interval.
  • FIG. 6 (b) relates to a speckle signal, and shows the result of double density scanning after data acquisition with a solid line for double density scanning and a dotted line for normal scanning. It can be seen that the signal with higher spatial frequency components is lost in the result of linear interpolation. This result gives the impression that the tomographic image causes a lateral flow.
  • Fig. 9 shows the results of evaluation by principle experiments regarding the effect of noise on the data after interpolation.
  • the ideal beam before adding noise was approximated by a sine function, and evaluation was performed by virtually adding delta function noise.
  • the solid line shows the ideal beam data
  • the dotted line shows the linear interpolation data
  • the broken line force nc function interpolation data.
  • sine function interpolation is closer to the original data than linear interpolation.
  • Figure 9 (b) shows the result of adding noise at the position of the arrow in the figure. It can be seen that the range affected by noise differs between linear interpolation and sine function interpolation.
  • linear interpolation is used when the effective dynamic range of the signal is wide, and sine function interpolation is used when the effective dynamic range of the signal is narrow. Therefore, an adaptive interpolation method is considered effective. Since the signal size and the electrical noise level specific to the device do not change significantly when the target, target part, and transmission focus position are determined, the switching depth may be determined in advance for each mode. I can do it.
  • a region to be imaged is selected.
  • optimized imaging parameters are selected for each target region such as the liver, kidney, uterus, fetus, and aorta. Switching this parameter set selection is referred to here as switching modes.
  • the attenuation rate associated with the propagation of ultrasonic waves varies greatly depending on the target site. For example, in the case of a fetus, most of the propagation medium is amniotic fluid, so Does not decay.
  • the propagation route is mostly occupied by the subcutaneous fat and the real organ called the liver, so the attenuation rate is much larger than when the amniotic fluid is transmitted.
  • the slope with respect to the depth of the echo signal changes depending on the mode, so that the depth within the dynamic range where beam width> raster interval changes.
  • the attenuation rate varies depending on the type of disease and the degree of progression (eg, cirrhosis) even at the same imaging site, so it is desirable to be able to finely adjust the mode including the lesion and the degree of progression.
  • the ultrasonic imaging apparatus of the present embodiment has a depth storage unit 21, and the depth storage unit 21 ⁇ This is the interpolation method switching depth, that is, the depth at which linear interpolation and sine function interpolation should be switched.
  • the data interpolation unit 30 refers to the interpolation method switching depth stored in the depth storage unit 21.
  • the data interpolation unit 30 performs data interpolation using the linear interpolation method. If it is deeper than the stored switching depth, data interpolation is performed using the sine function interpolation method.
  • the switching position may appear as an artifact.
  • a switching position control input unit 33 may be provided in the operation panel of the diagnostic machine 32 so that the operator can control the position of the switching point according to his / her preference.
  • the switching position may be changed in the depth direction, or the distribution of wl and ws may be changed.
  • Fig. 13 shows the input of one parameter with the switching position as the depth position.
  • Fig. 14 shows the switching of the distribution of wl and ws to the switching point (a in the figure).
  • the weight distribution change slope (b in the figure) and two parameters are input.
  • the position of the switching point of the interpolation method is fixed when the imaging mode is determined.
  • the progression of the medical condition such as cirrhosis
  • the position of the optimal switching point varies greatly depending on the variation between subjects such as the condition of the skin and the thickness of the subcutaneous fat. Therefore, it is preferable to set the optimal switching point position according to the acquired signal.
  • the ultrasound diagnostic device is called time gain control, and the gain can be adjusted for each signal depth, so parameters that determine the imaging mode of the device (target part, transmission focus position, transmission / reception) The frequency of the signal alone may not completely determine the dynamic range of the signal.
  • Fig. 15 shows an example of the configuration of an apparatus corresponding to such a case. Since the apparatus of the present embodiment is the same as the apparatus shown in FIG. 1 except that the interpolation method switching determination unit 22 is provided, description of common parts is omitted.
  • the interpolation method switching determination unit 22 obtains a change profile of the signal intensity in the depth direction as shown in FIG. Find the position below the dynamic range. This point is set as an interpolation method switching position.
  • the solid line in Fig. 10 is the actual data, and the dotted line is the data obtained by calculating the trend in the depth direction.
  • FIG. 16 shows an example of the configuration of an apparatus in which interpolation is performed after SZN is improved by a signal processing technique. Specific processing contents in the noise removal filter 31 in FIG. 16 will be described with reference to FIGS.
  • FIG. 17 is a diagram for explaining a target pixel for noise removal processing and a range in which a weight is calculated.
  • FIG. 18 is a flowchart of noise removal processing.
  • the memory in the noise removal filter 31 includes a one-dimensional image relating to 1, 2,..., N scanning lines, which respectively change in the direction of the time axis t. Two-dimensional data in which the data is arranged in the direction in which the scanning lines are arranged is stored. [0033] Range for calculating the weight by surrounding the target pixel (luminance I) for noise removal processing (weight calculation range)
  • the shape of the weight function is set in the parameter setting step 101 in FIG.
  • the width of the frequency distribution of the luminance difference obtained by obtaining the difference from the luminance I of all the target pixels in the above processing is obtained.
  • the width of the frequency distribution of this luminance difference is used to set the weighting function described later.
  • Pixels in the calculation range are set.
  • the weight is calculated based on a weight function described later. If this is performed for all the pixels within the weight calculation range, when it is determined in the weight calculation end determination step 104, the luminance value to be set for the extraction target pixel is obtained in the luminance value replacement step 105. It is done.
  • the target pixel end determination step 106 When the calculation is performed for all the pixels of the two-dimensional data to be subjected to the extraction process of the structure of the living tissue by shifting the position of the target pixel of the extraction process, it is determined in the target pixel end determination step 106 At this time, the noise removal process ends.
  • FIG. 19 is a diagram for explaining the weight function.
  • FIG. 19 (a) is a diagram showing the frequency distribution of the luminance difference described above in a typical ultrasonic image, and the horizontal axis is the luminance difference (I ⁇ I) described above, and the vertical axis Is the frequency of appearance of the luminance difference (I-I).
  • Figure 19 (b) shows the luminance difference (
  • the weight W is calculated according to the value of the horizontal axis (I I).
  • the integral value of the absolute value w of the function in the range of infinity is finite. Using this weight function, the brightness is calculated as follows.
  • This denoising filter works by continuity of pixel brightness to surrounding pixels. Is different. When the difference between the luminance I and the luminance I of the surrounding pixels is small, the weight is almost constant.
  • the weight is increased for the pixels along the interface between the two structures, so in the direction along the interface between the two structures, a one-dimensional low-pass filter, two structures In the direction perpendicular to the interface, it works as a global filter, so the sharpness of the interface will not be dulled.
  • non-linear filters that have different functions depending on the luminance distribution shape of the pixels in this way, the effective dynamic range of the signal is reduced by eliminating noise while minimizing the effect of smoothing the edges in the image. You can zoom in and change the switching position from linear interpolation to sine function interpolation deeper.
  • an ultrasonic tomographic image having a good spatial resolution and a signal-to-noise ratio can be obtained.

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Remote Sensing (AREA)
  • Medical Informatics (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)

Description

明 細 書
超音波撮像装置
技術分野
[0001] 本発明は、超音波断層像を表示する装置に関する。
背景技術
[0002] 超音波撮像にお!、ては、狭 、口径から、広 、視野を確保するために、セクタスキヤ ンもしくはコンベックススキャンと呼ばれる扇形のビーム走査を行う方法が広く利用さ れている。視野が広い方が病変を見つけるのに有利である力 以下のように口径が 狭い方が有利な場合がある。すなわち、心臓の撮像などのように、肋骨の隙間に探 触子を押し当てる場合など、被検体に探触子を接触させることの出来る範囲が限定 されている場合などには探触子の表面積が小さい方が有利な場合がある。しかし、探 触子の表面積を小さくすると走査線本数が減少してしま 、、少な 、走査線本数の条 件のままビデオ表示装置へのスキャンコンバージョンを行うと画像の劣化が起こる。 画質劣化を防止するために、特許文献 1に開示されるような、受信信号を複素信号 に変換し、複素信号の実数部と虚数部について個別に sine関数を用いた補間処理 を行う方法がある。
特許文献 1:特開平 11― 9603号公報
発明の開示
発明が解決しょうとする課題
[0003] しかし、 sine関数を用いる場合、一つの補間データを求めるのに用いる走査線の数 力 単純な線形補間に比べて増えるため、ある走査線上のデータにノイズが混入した 場合、その影響が遠くまで及んでしまう。ノイズが入った場合には、必ずしも sine関数 補間が最適とは言えない場合があり、 sine関数の適用の良否について判定するという 未解決な問題が依然としてあった。
[0004] 本発明の目的は、走査線上に混入するノイズの除去に寄与するために最適の補間 処理を適用可能とした超音波診断装置を提供することにある。
課題を解決するための手段 [0005] 本発明では、取得した信号に対して、信号のダイナミックレンジを測定し、信号に対 して適応的に線形補間力もしくは sine関数補間を用いるかの選択、もしくは両者の重 み付け平均を行う。線形補間と sine関数補間を切り替える深さ位置を予め装置に設 定しておき、その補間方式切り替え深さを境に 2つの補間方式を切り替えてもよい。 補間方式切り替え深さを境に 2つの補間方式を完全に切り替えるのではなぐ 2つの 補間方式による補間値を重み付け加算したものを補間値としてもよい。
図面の簡単な説明
[0006] [図 1]本発明による超音波診断装置の一例を示す機能ブロック図。
[図 2]スキャンコンバージョンの原理説明図。
[図 3]線形補間の説明図。
[図 4]sinc関数 間の説明図。
[図 5]ビーム幅に関する説明図。
[図 6]補間結果を示す実験データの一例の図。
[図 7]信号のダイナミックレンジとビーム幅の関係を示す説明図。
[図 8]信号のダイナミックレンジとビーム幅の関係を示す説明図。
[図 9]ノイズの影響が補間方式によって異なることを説明する図。
[図 10]信号のダイナミックレンジと補間方式の切り替え方を示す説明図。
[図 11]モードによって減衰率が異なることと、それに対応した補間点切り替え位置に 関する説明図。
[図 12]補間方式切り替えの重み関数の変化の仕方を説明する図。
[図 13]切り替え位置制御入力部の模式図。
[図 14]切り替え位置制御入力部の模式図。
[図 15]本発明による超音波診断装置の一例を示す機能ブロック図。
[図 16]本発明による超音波診断装置の一例を示す機能ブロック図。
[図 17]ノイズ除去フィルタの説明図。
[図 18]ノイズ除去フィルタの制御フロー図。
[図 19]ノイズ除去フィルタの重み関数を説明する図。
符号の説明 [0007] 1…超音波探触子
2…送受切り替えスィッチ
3…送波ビームフォーマ
4…制御系
20…受波ビームフォーマ
21 · ··深さ記憶部
22· ··補間方式切り替え判定部
23· ··スキャンコンバータ
24…表示部
29· ··メモリ
30· ··データ補間部
31…ノイズ除去フィルタ
32…診断機本体
33· ··切り替え位置制御入力部
101· ··パラメータ設定工程
102· ··重み計算範囲設定工程
103…重み値計算工程
104· ··終了判定工程
105…輝度値変換工程
106· ··判定工程
発明を実施するための最良の形態
[0008] 以下、図面を参照して本発明の実施の形態を説明する。
[0009] 図 1は、本発明による超音波診断装置の一例を示す機能ブロック図である。最初に 、超音波診断装置における、画像化のための信号処理の流れを説明する。被検体の 表面に設置された超音波探触子 1に対して、送受切り替えスィッチ 2を介して、制御 部 4の制御のもと送信ビームフォーマ 3から送波電気パルスが送られる。このとき送信 ビームフォーマは、所望の走査線上に超音波ビームが進むように、探触子 1の各チヤ ネル間の遅延時間が適した状態になるように制御して 、る。この送信ビームフォーマ 3からの電気信号は超音波探触子 1にお 、て超音波信号に変換され、被検体内に 超音波パルスが送波される。被検体内において散乱された超音波パルスは、一部が エコー信号として再び超音波探触子 1によって受信され、電気信号に変換される。こ の受信された信号は送受切り替えスィッチ 2を介して受波ビームフォーマ 20に送られ 、所望の走査線上の所望の深さからのエコー信号が選択的に増強された、データと して、メモリ 29に格納される。データ補間部 30は、走査線の数を増やすように、実際 にデータを所得した走査線と走査線の間のデータを補間する。補間の方式は後に説 明するように、複数の方式の中から最適なものが選択される。この補間されたデータ がスキャンコンノ ータ 23に送られ、スキャンコンノ一ジョンが行われる。スキャンコン バージョン後のデータは、表示部 24に送られ、超音波断層像として表示される。
[0010] 以下、本発明の補間と、スキャンコンバージョンに関して説明を行う。まず、図 2を用 いて、セクタスキャンもしくはコンベックススキャンで得られたエコーデータから直交座 標系のビデオ画像データへのスキャンコンバージョン方法にっ 、て説明する。スキヤ ンコンバージョンでは、スキャンコンバージョン後のデータ点(図中の黒丸)の周囲の 複数のスキャンコンバージョン前のデータ点(同じく白丸)での輝度値を用いて、スキ ヤンコンバージョン後のデータ点での輝度値を算出する。このとき、スキャンコンパ一 ジョン後の空間サンプリング間隔がスキャンコンバージョン前の空間サンプリング間隔 に比べ、十分に広い場合には、スキャンコンバージョン方法に対して注意を払う必要 はないが、スキャンコンバージョン後の空間サンプリング間隔がスキャンコンパージョ ン前の空間サンプリング間隔に比べ、十分に広くない場合には、スキャンコンパージ ヨン方法によって出力画像が大きな影響を受ける。スキャンコンバージョンを、(1)補間 によりスキャンコンバージョン前のサンプリングを十分に細力べして、(2)スキャンコンパ 一ジョンする、というように 2段階に分けて取り扱うことにすると、スキャンコンパージョ ン方式の検討は、補間方式の検討に帰着する。
[0011] 超音波画像においては、送受波それぞれでフォーカスを行うが、回折効果によって 横方向へのビーム絞込みには限界がある。また、横方向に細かく走査線を配すると フレームレートが低下するため、通常は中心周波数での波長の半分力 4分の一の 間隔で横方向に走査線を配する。この結果、横方向のサンプリング間隔は、例えば 中心周波数が 3MHzの場合には 250 /z mから 125 /z m程度となる。一方、深さ方向 に関しては、搬送波の周波数に対して十分にサンプリングを行うため、通常 30から 4 OMHzで AZD変翻のサンプリングを行い、音速を 1500mZsとして計算すると、 2 0から 25 /z mのサンプリング間隔となる。このように、深さ方向に比べ、方位方向には サンプリングが極めて粗いため、二次元スキャンコンバージョンにおける補間は、方 位方向への一次元補間の問題として取り扱うことが出来る。
[0012] 一次元の補間方式としては、図 3に示す補間方法と図 4に示す補間方法が一般的 である。図 3は、補間後のデータを挟む最近接な二点の輝度の内挿から補間後のデ ータを求める線形補間の説明図であり、図 4は、 4点以上のデータを用いる sine関数 補間の説明図である。データ点を有限の幅の空間の代表点と考え、この有限幅を矩 形関数とみなし、矩形関数のフーリエ変換である sine関数を用いて補間係数を決定 するのが、 sine関数補間である。図 4に示す破線位置での補間データの輝度を求め るには、隣接する画素での輝度値に、隣接する画素にピークを持つ sine関数(図 4の 実線及び点線)からの係数 2Ζ πを乗じたものと、その外側の画素での輝度値に、そ の画素位置にピークを持つ sine関数(図 4の破線及び一点鎖線)からの係数— 3Ζ2 πを乗じたものの総和が補間値となる。
[0013] データのサンプリング周波数力 データの周波数に比べ十分に細かい場合はいず れの方法を用いても結果は同等である力 十分に細力べない場合は、 sine関数を用 いた補間の方が、補間による誤差が少ないことが知られている。補間位置が不変で ある場合は、 sine関数を使った、あらかじめ求めた係数を使えばよいので、(上記のよ うに、 2Ζ π、 - 3Ζ2 πなど) sine関数を用いることによる、計算上の負荷は発生しない 。しかし、補間位置が一定でない場合は、補間位置ごとに sine関数を用いて補間係 数を求めることは計算負荷が大きくなつてしまう。実際上は、必ずしも sine関数ではな くても、 sine関数を taylor展開して、有限の項で打ち切った近似的な関数を用いても、 実用上は問題ない。特に、実際上は DSP上で超音波診断装置のフレームレートに 間に合うように処理をする必要があるので、有限の項で打ち切る方が、むしろ実用的 には好ましい。この場合、例えば、位置 aを中心として、 sine関数の 3次までの Taylor 展開は、 Sin(a)/ a+(x— a)(cos(a)/ a— sin(a)/ a )+(x— a)2/2x(— sin(a)/ a— 2cos(a)/ a2+2cos(a)/ a )+(x— a)3/ 6x(— cos(a)/ a+3sin(a)/ a2+4cos(a)/ a3— 2sin(a)/ a3— 6cos(a)/ a4)
となるので、 π Ζ2の位置を中心に展開すると、
(2/ π 2—8/3 π 3) χ3+(4/ π 2— 4/ π )χ2+(— 4/ π 2— 2/ π +5/2)χ+(4/ π— π /2+1/3) と高々 3次の多項式で近似できる。隣接点からの補間、二つとなりの点からの補間な ど、補間位置によって、展開の中心値異なるようにしたほうが、展開係数で必要な次 数を下げることが出来る。
[0014] 超音波画像は、ビームフォーミング及びポストプロセス処理で決定される装置の撮 像条件固有の点応答関数と、被写体の散乱体分布のコンボリューシヨンである。ビー ムフォーミングは送波と受波それぞれで行われるが、現状の超音波撮像においては 、送波に関してはフレームレートを低下させないことを優先して、一つの走査線上で は 1箇所にフォーカスするのみで、それ以外の場所ではオフフォーカスとなっている。 一方受波では、ダイナミックフォーカスと呼ばれる、受波の時間に応じて連続的にフ オーカス位置を切り替え、深度方向で一様にオンフォーカスとなるような装置構成を 実現している。
[0015] 横流れが問題となるような生体深部での点応答関数の広がり、すなわち方位方向 の分解能は、受波ビーム幅で決まる。焦点位置でのビーム幅 BWは、幾何学的なビ ーム幅よりは、図 5に示すように、超音波の回折効果で決まる。回折角 Θは次式で近 似出来る。ここで、 λは波長、 Dは口径幅である。
[0016] Θ =sin_1 ( l /D)
以下典型的な例で計算すると、深部での中心周波数 2MHz、送波口径重みを近 似して半分の口径幅として扱うと、口径幅 12. 5mm,探触子表面から撮像部位まで の距離を Zとした場合、回折角は 0. 06radとなるので、ビーム幅は 0. 06 X Zmmであ る。一方、走査線の幅は、典型的には 0. Ol X Zmm程度である。すなわち、ビーム 幅に対して約 6倍細力べサンプリングしていることになるので、通常は線形補間で十分 であると考えられる。
[0017] 実験結果を図 6に示す。縦軸は輝度、横軸は横方向の位置を示す。図 6 (a)は点反 射体からのエコー信号に関するもので、実線で結ばれている黒丸が倍密度走査、点 線で結ばれている白丸が通常走査でデータの取得を行ったあと、線形補間により倍 密度にした結果である。いずれの結果もほとんど差がなぐよく重なっており、このケ ースでは線形補間によって画像の劣化は起きていない。これは既に述べたように、ビ ーム幅が走査線間隔より十分に広い場合として予測されたとおりである。
[0018] 一方、図 6 (b)はスペックル信号に関するもので、実線が倍密度走査、点線が通常 走査でデータの取得を行ったあと、線形補間により倍密度にした結果である。線形補 間された結果の方が、高い空間周波数成分の信号が失われていることがわかる。こ の結果が、断層像が横流れを起こして 、る印象を与えて 、る。
[0019] まず、当初の見積もりと異なり、スペックル信号に対しては線形補間では横流れが 起き、 sine関数補間では横流れが起きな力つた原因に関して考察を行う。これは、信 号のダイナミックレンジが有限であることを考えると説明できる。現状の超音波診断装 置は、 AZD変翻でおよそ 150dBと有限なダイナミックレンジを持っている。そのた め、図 7に示すように、散乱体の反射強度が大きい場合には、ビーム幅と反射信号の 空間的な幅は一致する。しかし、スペックル信号のように反射強度が比較的小さぐ かつ深部力 のエコーのように超音波の生体中の伝播に伴う減衰によって信号強度 力 S小さくなつてしまった場合には、図 8に示すように、有効ダイナミックレンジ内での信 号の幅は狭くなつてしまう。このように、図 8に示すようなケースでは、補間方式は線形 補間よりむしろ sine関数補間を用いる必要がある。
[0020] ノイズが補間後のデータに与える影響に関して、原理実験によって評価した結果を 図 9に示す。ノイズが加わる前の理想ビームを sine関数で近似し、仮想的にデルタ関 数的なノイズを加えて評価を行った。図 9 (a)はノイズが無い場合で、実線が理想ビ ームデータ、点線が線形補間データ、破線力 nc関数補間データを示す。ノイズが無 い場合は、 sine関数補間の方が線形補間より元データに近い。図 9 (b)は、図中の矢 印の位置にノイズをカ卩えた結果である。ノイズの影響の出る範囲が線形補間と sine関 数補間で異なることが解る。つまり、線形補間の場合は、ノイズが入ったときに、その ノイズが補間に影響を与える範囲は、ノイズの入った位置に隣接する補間点のみで あるが、 sine関数補間の場合にはノイズの入った位置に隣接する補間点及びそのもう 一つ隣の点にまで影響がでる。 この補間方式がノイズに与える影響のシミュレーション評価の結果を表 1にまとめる 。評価は 201og (SinC関数補間での誤差 Z線形補間での誤差)を評価しており、 dBが 正の場合は線形補間が適しており、負の場合は sine関数補間が適していることを示し ている。
[表 1]
Figure imgf000010_0001
[0022] この結果は、 sine関数補間を行う場合は隣接 4点のデータを使うため、隣接 2点しか 用いない線形補間に比べノイズがより多くの点に対して影響を与えてしまうことを示し ており、信号の周波数成分に対して十分細力べサンプリング出来ているときは、 sine関 数補間よりむしろ線形補間が優れていることを示している。
[0023] この結果、図 10に示すように、信号の有効なダイナミックレンジが広い場合は線形 補間を用い、信号の有効なダイナミックレンジが狭 、場合には sine関数補間を用いる という、信号に対して適応的な補間方法が有効であると考えられる。信号大きさ及び 装置固有の電気ノイズのレベルは、対象物、対象部位、送波フォーカス位置が決ま つた場合には大きな変化がないので、予めモードごとに、切り替え深さを決めておくこ とが出来る。
[0024] 通常超音波診断装置では、接続する超音波探触子を選んだときに、撮像対象部位 を選択する。例えば低周波コンペックススキャン型探触子の場合、肝臓、腎臓、子宮 、胎児、大動脈など対象部位毎に、それぞれに最適化した撮像用パラメータを選択 するようになつている。このパラメータセットの選択を切り替えることを、ここではモード を切り替えるという。このような対象部位によって、超音波の伝播に伴う減衰率が大き く異なる。例えば、胎児の場合は、伝播媒質のほとんどは、羊水であるため、ほとんど 減衰しない。一方肝臓の場合などは、伝播経路はほとんど、皮下脂肪と、肝臓という 実質臓器で占められているため、羊水を伝播する場合に比べ、減衰率がはるかに大 きい。このため、図 11に示すように、モードによって、エコー信号の深さに対する傾き が変わってくるので、ビーム幅 >ラスター間隔となるダイナミックレンジになる深さが変 化する。同一撮像部位でも病気の種類や進行度合いによって、(例えば肝硬変の場 合など)減衰率が変わってくるので、病変や、進行度合いも含めたモードの細かい調 節が可能であることが望ま 、。
[0025] このため、本実施例の超音波撮像装置は深さ記憶部 21を有し、深さ記憶部 21〖こ 補間方式の切り替え深さ、すなわち線形補間と sine関数補間を切り替えるべき深さの データを格納しておく。データ補間部 30は、深さ記憶部 21に格納された補間方式の 切り替え深さを参照し、データ点の深さが記憶している補間方式の切り替え深さより 浅いときは線形補間方式でデータ補間を行い、記憶している切り替え深さより深いと きは sine関数補間方式によってデータ補間を行う。また、切り替え深さで、急峻に線 形補間から sine関数補間に切り替えると、切り替え位置がアーチファクトとして現れて しまう場合がある。これを防ぐ方法として、以下のように重み値を連続的に変えていく 方法がある。線形補間値を II、これに対する重み値 wl、 sine関数補間値を Is、これに 対する重み値 wsを使って、補間結果の出力 Iを次式のように表し、 wlと wsを図 12に 示すように切り替え点で連続的に切り替わるように構成することが出来る。
[0026] I=wl X Il+ws X Is
なお、この切り替え点の位置は、オペレータが好みによって制御出来るように、診断 機 32の操作パネルの中に切り替え位置制御入力部 33を設けてもよい。その場合、 深さ方向に切り替え位置が変わるようにしてもよいし、前記の wlと wsの配分が変更さ れるようにしても良い。図 13は、切り替え位置を深さの位置として、一つのパラメータ を入力する形にしたものであり、図 14は、 wlと wsの配分の切り替えを、切り替え点(図 中の a)と切り替え時の重み配分変更の傾き(図中の b)と二つのパラメータで入力す る形にしたものである。
[0027] 図 1に示した装置の場合、補間方法の切り替え点の位置は、撮像モードが決まった ら固定であった。し力 実際には、例えば肝臓の場合などでも、肝硬変など病状の進 み具合や、皮下脂肪の厚みなど、被検体間のばらつきによって、最適な切り替え点 の位置は大きく異なる。そのため、最適な切り替え点の位置は、取得した信号に応じ て設定するのが好ましい。特に、超音波診断装置は、 time gain controlといって、信 号の深さ毎にゲインを調節出来るようになっているので、装置の撮像モードを決める パラメータ (対象部位、送波フォーカス位置、送受信の周波数など)だけでは、信号の ダイナミックレンジが完全に決定されるとは限らない。
[0028] このような場合に対応した装置の構成例を図 15に示す。本実施例の装置は、補間 方式切り替え判定部 22があること以外は図 1に示した装置と同じであるので、共通す る部分の説明は省略する。補間方式切り替え判定部 22では、補間の対象となる走査 線に沿った 1次元データが入力されると、図 10に示すような、深さ方向への信号強度 の変化プロファイルを求め、予め設定したダイナミックレンジを下回る位置を求める。 この点を補間方式切り替え位置とする。図 10の実線が実際のデータで、点線が深さ 方向のトレンドを算出したデータである。このトレンドの算出にあたっては、実線のデ ータに対して、ローパスフィルタやメディアンフィルタのようなフィルタを用いてもよいし 、深さ Xでのエコー強度 I(x)力 I(x) = ax + bで表現できると仮定して、係数 a及び bを 最小自乗フィッティングで求めてもよ!、。
[0029] 一つの走査線毎に補間方式切り替え点が大きく異なると、画像に縞々状のアーチ ファクトが出てしまうので、複数走査線の平均データに対して信号強度の変化プロフ アイルを求める方法や、走査線毎に求まった補間方式切り替え位置に対して、走査 線間での平均値を求めて、これを使う方法などがある。
[0030] 次に、信号処理技術によって、 SZNを改善してから補間を行う場合の装置の構成 例を図 16に示す。図 16のノイズ除去フィルタ 31での具体的な処理内容について、 図 17, 18, 19を用いて説明する。
[0031] 図 17は、ノイズ除去処理の対象画素と、重みを計算する範囲を説明する図である。
図 18はノイズ除去処理のフロー図である。
[0032] ノイズ除去フィルタ 31内のメモリには、図 17に示すように、時間軸 tの方向にそれぞ れ変化する、 1, 2, ···, N本の走査線に関する 1次元の画像データを走査線の並ぶ 方向に並べた 2次元データが記憶されて 、る。 [0033] ノイズ除去処理の対象画素 (輝度 I )を囲み、重みを計算する範囲 (重み計算範囲)
0
を設定する。重み計算範囲は、 i= l, 2, · ··, i 、及び、 j = l, 2,…, j からなる対
max max
象画素の周辺画素 (輝度 Iをもつ、 i Xj 個の画素)を含む。重み計算範囲が大き
ij max max
いほどノイズ除去フィルタの効果は大きいが、その分、演算速度は遅くなる。 i と j
max max
、重み関数の形状が、図 18のパラメータ設定工程 101で設定される。パラメータ設定 工程 101では、対象画素の輝度 Iと、この対象画素を囲む重み計算範囲の内の画素
0
の輝度 Iとの差を、上記の処理の全ての対象画素に対して求めて得られる輝度差の 度数分布の幅を求める。この輝度差の度数分布の幅を後述の重み関数を設定する ために用いる。
[0034] 重み計算範囲の設定工程 102で、対象画素の位置と、 i tj によって定まる重み
max max
計算範囲の画素が設定される。重み計算工程 103で、後述の重み関数に基づいて 重みの計算が行なわれる。これが重み計算範囲内の全ての画素に対して行われたと 、重み計算の終了判定工程 104で判定された時、輝度値置換工程 105で、上記の 抽出処理の対象画素に設定する輝度値が求められる。上記の抽出処理の対象画素 の位置をシフトして、生体組織の構造の抽出処理を行なうとする 2次元データの全て の画素に関して計算が行なわれたと、対象画素の終了判定工程 106で判定された 時、ノイズ除去処理が終了する。
[0035] 図 19は、前記の重み関数を説明する図である。図 19 (a)は、典型的な超音波画像 における、先に説明した輝度差の度数分布を示す図であり、横軸は先に説明した輝 度差 (I—I )であり、縦軸は輝度差 (I—I )の出現度数である。図 19 (b)は、輝度差 (
0 ij 0 ij
I I )の絶対値が大きくなるほど重みが単調に小さくなる重み関数 Wの例を示して
0 ij
おり、横軸 (I I )の値に対応して、重み Wが計算される。重み関数 Wの例として偶
0 ij ij
数次の多項式はじめ、その他、ガウシアン関数、 l/ (x2 + a2)のような関数等の多様 な偶関数が有り得る。関数 Wは (I I ) =0の時に極大点をもち、負の無限大から正
0 ij
の無限大の範囲での関数の絶対値 wの積分値は有限である。この重み関数を用い て、輝度の計算は次式のように行われる。
[0036] I +∑{ (I -I )W }/∑W
0 ij 0 ij ij
このノイズ除去フィルタは、画素の輝度の周囲画素への連続性によって、その機能 が異なる。輝度 Iと周囲画素の輝度 Iの差が小さい時は、重みがほぼ一定の値になる
0 ij
ので、 2次元の低域フィルタとして働き、輝度 Iをもつ画素が、 2つの構造物 (組織)の
0
界面に位置する場合は、 2つの構造物の界面に沿った画素に対して重みが大きくな るので、 2つの構造物の界面に沿った方向では、 1次元の低域フィルタ、 2つの構造 物の界面に直交する方向では、全域フィルタとして働くので、界面の鋭さを鈍らせる ことはな 、。このように画素の輝度の分布形状によって機能が異なる非線形フィルタ を用いることにより、画像中のエッジをなまらせる効果を最小限に抑えて、ノイズ除去 を行うことによって、信号の実効的なダイナミックレンジが拡大し、線形補間から sine 関数補間への切り替え位置をより深部に変えることが出来る。
産業上の利用の可能性
本発明によると、空間分解能と信号対雑音比の良好な超音波断層像を得ることが できる。

Claims

請求の範囲
[1] 超音波ビームの送波と超音波エコー信号の受波を行う超音波探触子と、
前記超音波探触子から扇形に走査される超音波ビームが発生されるように前記超 音波探触子に信号を送信する送波ビームフォーマと、
前記超音波探触子から扇状に走査される各走査線上の複数のデータ点からのェコ 一信号を取得する受波ビームフォーマと、
前記複数のデータ点からのエコー信号を各データ点のデータとして格納するメモリ と、
前記メモリに格納された走査線間のデータを補間するデータ補間部とを備え、 前記データ補間部は、データ点の深さに応じた複数の補間方式を用いてデータ補 間を行うことを特徴とする超音波撮像装置。
[2] 請求項 1に記載の超音波撮像装置において、前記複数の補間方式を切り替える深 さを入力する入力部を有することを特徴とする超音波撮像装置。
[3] 請求項 1に記載の超音波撮像装置にお!、て、前記複数の補間方式は線形補間と 二次以上の高次関数補間を含むことを特徴とする超音波撮像装置。
[4] 請求項 1に記載の超音波撮像装置にお!、て、前記複数の補間方式は線形補間と si nc関数補間であることを特徴とする超音波撮像装置。
[5] 請求項 3に記載の超音波撮像装置において、前記線形補間の結果と前記二次以 上の高次関数による補間の結果を重み付け加算したものを補間データとし、データ 点の深さが予め決められた補間方式切り替え深さより浅いときは前記線形補間の結 果に対する重みを大きくし、データ点の深さが前記補間方式切り替え深さより深いと きは前記二次以上の高次関数による補間の結果に対する重みを大きくすることを特 徴とする超音波撮像装置。
[6] 請求項 1に記載の超音波撮像装置において、前記データ補間を行う前にノイズ除 去を行うことを特徴とする超音波撮像装置。
[7] 請求項 1記載の超音波撮像装置において、モード毎に前記複数の補間方式の切り 替え深さを記憶して ヽることを特徴とする超音波撮像装置。
[8] 超音波ビームの送波と超音波エコー信号の受波を行う超音波探触子と、 前記超音波探触子から扇形に走査される超音波ビームが発生されるように前記超 音波探触子に信号を送信する送波ビームフォーマと、
前記超音波探触子から扇状に走査される各走査線上の複数のデータ点からのェコ 一信号を取得する受波ビームフォーマと、
前記複数のデータ点からのエコー信号を各データ点のデータとして格納するメモリ と、
前記走査線に沿った深さ方向への前記エコー信号強度の変化プロファイルを求め 、当該信号強度が予め定めた信号強度を下回る深さ位置を補間方式切り替え位置と する補間方式切り替え判定部と、
前記メモリに格納された走査線間のデータを、前記補間方式切り替え位置より浅 ヽ データ点のデータに関しては線形補間し、前記補間方式切り替え位置より深 、デー タ点のデータに関しては二次以上の高次関数で補間して走査線間のデータを補間 するデータ補間部と
を有することを特徴とする超音波撮像装置。
[9] 請求項 8に記載の超音波撮像装置において、前記二次以上の高次関数は sine関 数であることを特徴とする超音波撮像装置。
[10] 請求項 8に記載の超音波撮像装置において、前記線形補間の結果と前記二次以 上の高次関数による補間の結果を重み付け加算したものを補間データとし、データ 点の深さが前記補間方式切り替え位置より浅いときは前記線形補間の結果に対する 重みを大きくし、データ点の深さが前記補間方式切り替え位置より深いときは前記二 次以上の高次関数による補間の結果に対する重みを大きくすることを特徴とする超 音波撮像装置。
[11] 請求項 8に記載の超音波撮像装置において、前記線形補間の結果と前記二次以 上の高次関数により補間した結果を重みを着けて合成する手段と、重み値を制御す る入力部を有する超音波撮像装置。
PCT/JP2007/062021 2006-07-18 2007-06-14 Ultrasonographic device Ceased WO2008010366A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US12/374,065 US8506483B2 (en) 2006-07-18 2007-06-14 Ultrasonographic device
JP2008525809A JP4837037B2 (ja) 2006-07-18 2007-06-14 超音波撮像装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2006195674 2006-07-18
JP2006-195674 2006-07-18

Publications (1)

Publication Number Publication Date
WO2008010366A1 true WO2008010366A1 (en) 2008-01-24

Family

ID=38956701

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2007/062021 Ceased WO2008010366A1 (en) 2006-07-18 2007-06-14 Ultrasonographic device

Country Status (3)

Country Link
US (1) US8506483B2 (ja)
JP (1) JP4837037B2 (ja)
WO (1) WO2008010366A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103776907A (zh) * 2014-01-03 2014-05-07 浙江大学 基于sinc插值的超声相控阵接收信号精延时方法
WO2016152602A1 (ja) * 2015-03-23 2016-09-29 富士フイルム株式会社 音響波画像生成装置およびその制御方法

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101055577B1 (ko) * 2007-11-23 2011-08-23 삼성메디슨 주식회사 초음파 시스템
EP2605035B1 (en) 2010-11-10 2016-10-26 Samsung Medison Co., Ltd. Enhancing quality of ultrasound image in an ultrasound system via image filtering
KR102605151B1 (ko) * 2016-02-16 2023-11-23 삼성메디슨 주식회사 빔포밍을 수행하는 방법 및 빔포머.
US11980495B2 (en) * 2020-04-28 2024-05-14 GE Precision Healthcare LLC Method and system for providing enhanced color flow doppler and pulsed wave doppler ultrasound images by applying clinically specific flow profiles

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05154153A (ja) * 1991-12-04 1993-06-22 Fuji Electric Co Ltd 超音波診断装置
JPH0795978A (ja) * 1993-09-30 1995-04-11 Shimadzu Corp 超音波診断装置
JPH11276477A (ja) * 1998-01-28 1999-10-12 Hitachi Medical Corp 超音波装置
JP2001286472A (ja) * 2000-04-10 2001-10-16 Toshiba Corp 超音波診断装置
JP2005087266A (ja) * 2003-09-12 2005-04-07 Fuji Photo Film Co Ltd 超音波撮像装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6443238A (en) * 1987-08-12 1989-02-15 Toshiba Corp Ultrasonic blood flow imaging apparatus
US4787252A (en) * 1987-09-30 1988-11-29 Panametrics, Inc. Differential correlation analyzer
JP3806229B2 (ja) * 1997-06-27 2006-08-09 アロカ株式会社 超音波診断装置
US7347548B2 (en) * 2003-05-01 2008-03-25 The Cleveland Clinic Foundation Method and apparatus for measuring a retinal sublayer characteristic

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05154153A (ja) * 1991-12-04 1993-06-22 Fuji Electric Co Ltd 超音波診断装置
JPH0795978A (ja) * 1993-09-30 1995-04-11 Shimadzu Corp 超音波診断装置
JPH11276477A (ja) * 1998-01-28 1999-10-12 Hitachi Medical Corp 超音波装置
JP2001286472A (ja) * 2000-04-10 2001-10-16 Toshiba Corp 超音波診断装置
JP2005087266A (ja) * 2003-09-12 2005-04-07 Fuji Photo Film Co Ltd 超音波撮像装置

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103776907A (zh) * 2014-01-03 2014-05-07 浙江大学 基于sinc插值的超声相控阵接收信号精延时方法
WO2016152602A1 (ja) * 2015-03-23 2016-09-29 富士フイルム株式会社 音響波画像生成装置およびその制御方法
US10925578B2 (en) 2015-03-23 2021-02-23 Fujifilm Corporation Acoustic wave image generating apparatus and control method thereof

Also Published As

Publication number Publication date
JP4837037B2 (ja) 2011-12-14
US20100049042A1 (en) 2010-02-25
US8506483B2 (en) 2013-08-13
JPWO2008010366A1 (ja) 2009-12-17

Similar Documents

Publication Publication Date Title
JP4757307B2 (ja) 超音波画像処理装置
US10729407B2 (en) Ultrasonic diagnostic apparatus, image processing apparatus, and image processing method
KR101205107B1 (ko) 스페클 감소 필터의 구현 방법, 스페클 감소 필터링 장치 및 초음파 촬상 시스템
JP6342212B2 (ja) 超音波診断装置
US6790181B2 (en) Overlapped scanning for multi-directional compounding of ultrasound images
JP4632685B2 (ja) 超音波診断装置及び画像データ処理装置
EP2905633B1 (en) Ultrasonic diagnosis apparatus, image processing apparatus, and image processing method
US20130258805A1 (en) Methods and systems for producing compounded ultrasound images
CN106469461B (zh) 视图方向自适应体积超声成像
JP6309340B2 (ja) 超音波診断装置及び超音波イメージングプログラム
JP2004129773A (ja) 超音波イメージング装置及び超音波信号処理方法
WO2008010366A1 (en) Ultrasonographic device
CN101088468B (zh) 超声波图像取得装置、图像处理装置和超声波图像处理方法
KR101051555B1 (ko) 개선된 3차원 초음파 영상을 형성하는 초음파 영상 장치 및방법
JP2022023982A (ja) 超音波画像を処理するための方法及びシステム
JP2019097794A (ja) 超音波診断装置、医用画像処理装置及びそのプログラム
JP2019097795A (ja) 超音波診断装置、医用画像処理装置及びそのプログラム
JP2008212678A (ja) 周波数領域でデジタルスキャン変換を実施して超音波映像を形成する方法
JP6553140B2 (ja) 超音波診断装置、画像処理装置及び画像処理方法
JP2008220652A (ja) 超音波診断装置、及び超音波画像生成プログラム
JP7280713B2 (ja) 超音波診断装置
JP7171625B2 (ja) 超音波画像を処理するための方法及びシステム
JP2006212054A (ja) 超音波観測装置、画像処理装置、及び、画像処理プログラム
WO2019145188A1 (en) Ultrasound system with improved noise performance by persistence processing
JP2019025247A (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: 07745278

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2008525809

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

NENP Non-entry into the national phase

Ref country code: RU

WWE Wipo information: entry into national phase

Ref document number: 12374065

Country of ref document: US

122 Ep: pct application non-entry in european phase

Ref document number: 07745278

Country of ref document: EP

Kind code of ref document: A1