WO2005065547A1 - 超音波診断装置 - Google Patents

超音波診断装置 Download PDF

Info

Publication number
WO2005065547A1
WO2005065547A1 PCT/JP2005/000074 JP2005000074W WO2005065547A1 WO 2005065547 A1 WO2005065547 A1 WO 2005065547A1 JP 2005000074 W JP2005000074 W JP 2005000074W WO 2005065547 A1 WO2005065547 A1 WO 2005065547A1
Authority
WO
WIPO (PCT)
Prior art keywords
filter coefficient
data
filter
reception
diagnostic apparatus
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/JP2005/000074
Other languages
English (en)
French (fr)
Inventor
Akiko Uchikawa
Yoshihiko Ito
Yasuhiro Nakamura
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.)
Panasonic Holdings Corp
Original Assignee
Matsushita Electric Industrial Co Ltd
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 Matsushita Electric Industrial Co Ltd filed Critical Matsushita Electric Industrial Co Ltd
Priority to US10/585,388 priority Critical patent/US8246544B2/en
Priority to EP05703342A priority patent/EP1709908A4/en
Priority to JP2005516882A priority patent/JP4290699B2/ja
Publication of WO2005065547A1 publication Critical patent/WO2005065547A1/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/52025Details of receivers for pulse systems
    • G01S7/52026Extracting wanted echo signals
    • 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/52085Details related to the ultrasound signal acquisition, e.g. scan sequences
    • 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/52085Details related to the ultrasound signal acquisition, e.g. scan sequences
    • G01S7/52095Details related to the ultrasound signal acquisition, e.g. scan sequences using multiline receive beamforming
    • 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
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8979Combined Doppler and pulse-echo imaging systems
    • 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

Definitions

  • the present invention relates to an ultrasonic diagnostic apparatus having a parallel reception function of generating a plurality of reception beams having the same transmission beam power.
  • an ultrasonic diagnostic apparatus is based on two-dimensional information obtained by scanning an ultrasonic beam with an ultrasonic probe 30 that transmits and receives ultrasonic waves to and from a subject. It displays sound wave images.
  • the signal received by the ultrasonic probe 30 is supplied to a spatial filter processing circuit 32 via a signal processing unit 31.
  • the spatial filter processing circuit 32 is arranged at a stage before the image display conversion unit 33, and is configured by a spatial filter for reducing noise on a received signal.
  • the control circuit 34 controls the filter coefficient of the spatial filter processing circuit 32 so as to change from near the transducer plane of the acoustic data to the far point direction. In other words, in the conventional ultrasonic diagnostic apparatus, when performing spatial filter processing on acoustic data before coordinate transformation, the filter coefficient is changed according to the distance on the acoustic line data! Reference 1).
  • Patent Document 1 Japanese Unexamined Patent Publication No. 2001-340338 (FIG. 1)
  • the signal difference between a plurality of reception beams obtained with the same transmission beam power is the signal difference between reception beams obtained from different transmission beams.
  • the filter coefficient for reducing the image quality difference between adjacent beams is fixed. That is, between the received beam signals having the same transmitted beam power, the filter is too effective, and the image becomes uniform due to the smoothing effect, and the image cannot be displayed in detail.
  • the present invention has been made to solve the above conventional problem, and an object of the present invention is to optimize a filtering process for a signal between a plurality of reception beams obtained from the same transmission beam. Another object of the present invention is to provide an ultrasonic diagnostic apparatus capable of suppressing the occurrence of stripes in the direction in which acoustic lines are arranged and displaying an ultrasonic image of good image quality capable of expressing details.
  • an ultrasonic diagnostic apparatus stores reception beam data obtained by digitally converting a reception beam formed from an ultrasonic reception signal.
  • First control means memory control unit, first memory control unit for controlling reading and writing of data to and from the first storage means, and the positional relationship between the transmission beam and the reception beam.
  • the filter coefficient calculation unit that calculates the filter coefficient based on the included information, and the image quality difference between adjacent beams based on the filter coefficient for a plurality of reception beam data including the parallel reception beam obtained from the same transmission beam
  • a first spatial filter operation unit that performs a filtering process for reducing the noise, and scans the image data output from the first spatial filter operation unit on a display monitor. It is configured to display the image after converting it to.
  • the ultrasonic diagnostic apparatus further includes A two-dimensional Doppler signal processing unit that performs two-dimensional Doppler processing on received beam data, a second storage unit (two-dimensional Doppler memory) that stores two-dimensional Doppler data output from the two-dimensional Doppler signal processing unit, And a second control unit (second memory control unit) for controlling reading and writing of data to and from the second storage unit, and obtaining from the same transmission beam according to the filter coefficient supplied from the filter coefficient calculation unit. And a second spatial filter operation unit that performs filtering processing on a plurality of received two-dimensional Doppler data including parallel received beams to reduce image quality differences between adjacent beams. Can be done.
  • the filter coefficient calculation unit may be configured to be able to control the filter coefficient according to the reception depth.
  • the filter coefficient calculation section may be configured to be able to control the filter coefficient according to the reception beam angle.
  • the filter coefficient calculation unit may be configured to be able to control the filter coefficient according to the focal position of the transmission beam.
  • the filter coefficient can be optimized such that the filter coefficient has a weak correlation near the focal position of the transmission beam, and the filter coefficient has a strong correlation when the focal position of the transmission beam is far away. Accordingly, it is possible to display a favorable ultrasonic image having a uniform image quality regardless of the focal position of the transmission beam.
  • the present invention by optimizing a filtering process on a signal between a plurality of reception beams that have also obtained the same transmission beam power, it is possible to suppress the occurrence of fringes in the arrangement direction of acoustic lines, It is possible to provide an ultrasonic diagnostic apparatus capable of displaying an ultrasonic image having a high image quality that can be expressed to a high degree.
  • FIG. 1 is a functional block diagram showing a configuration example of an ultrasonic diagnostic apparatus according to a first embodiment of the present invention.
  • FIG. 2A is a schematic diagram showing image data at an arbitrary depth when a lateral filter is used and one receiving beam is generated from one transmitting beam in the first embodiment.
  • FIG. 2B is a schematic diagram showing filter coefficients corresponding to the image data in FIG. 2A.
  • FIG. 3A is a schematic diagram showing image data at an arbitrary depth when a diamond-shaped filter is used in the case where one transmission beam is generated from one transmission beam in the first embodiment. is there.
  • FIG. 3B is a schematic diagram showing filter coefficients corresponding to the image data in FIG. 3A.
  • FIG. 4A is a schematic diagram showing image data at an arbitrary depth in a case where two reception beams are generated from one transmission beam and a lateral filter is used in the first embodiment. is there.
  • FIG. 4B is a schematic diagram showing a filter coefficient when data OutL of an arbitrary sampling point on reception beam RX1L is generated for the image data of FIG. 4A.
  • FIG. 4C is a schematic diagram showing filter coefficients when generating data OutR of an arbitrary sampling point on reception beam RX1R with respect to the image data of FIG. 4A.
  • FIG. 5A shows two reception beams from one transmission beam in the first embodiment.
  • FIG. 9 is a schematic diagram showing image data of an arbitrary depth when a diamond-shaped filter is used in the case of generation.
  • FIG. 5B is a schematic diagram showing a filter coefficient when data OutL of an arbitrary sampling point on reception beam RX1L is generated for the image data of FIG. 5A.
  • FIG. 5C is a schematic diagram showing a filter coefficient when data OutR of an arbitrary sampling point on reception beam RX1R is generated for the image data of FIG. 5A.
  • FIG. 6 is a functional block diagram showing one configuration example of an ultrasonic diagnostic apparatus according to a second embodiment of the present invention.
  • FIG. 7 is a functional block diagram showing one configuration example of an ultrasonic diagnostic apparatus according to a third embodiment of the present invention.
  • FIG. 8 is a schematic diagram for explaining a filter coefficient control function in a case where a difference occurs in acoustic line intervals depending on depth in the third embodiment.
  • FIG. 9 is a functional block diagram showing a configuration example of an ultrasonic diagnostic apparatus according to a fourth embodiment of the present invention.
  • FIG. 10A is a schematic diagram for explaining a filter coefficient control function when transmitting and receiving an acoustic ray at an angle in the fourth embodiment.
  • FIG. 10B is a schematic diagram for explaining a filter coefficient control function in the case where there is a difference in deflection angle between acoustic lines in the fourth embodiment.
  • FIG. 11 is a functional block diagram showing one configuration example of an ultrasonic diagnostic apparatus according to a fifth embodiment of the present invention.
  • FIG. 12 is a schematic diagram for describing a filter coefficient control function based on a focal position of a transmission beam in a fifth embodiment.
  • FIG. 13 is a functional block diagram showing a configuration example of a conventional ultrasonic diagnostic apparatus. Explanation of symbols
  • Second memory control unit (second control means)
  • FIG. 1 is a functional block diagram showing one configuration example of the ultrasonic diagnostic apparatus according to the first embodiment of the present invention.
  • This ultrasonic diagnostic apparatus includes an ultrasonic reception data processing unit 1 that forms a reception beam of an ultrasonic reception signal and converts the reception beam into reception beam data, a spatial filter 9, and information on the reception beam.
  • An acoustic beam control unit 6 that supplies the image data to the ultrasonic reception data processing unit 1, a scan conversion unit 7 that converts output data from the spatial filter 9 into image data, and a monitor 8 (display unit) that displays the image data. Consists of
  • the spatial filter 9 includes a memory 3 (first storage unit) that stores the reception beam data supplied from the ultrasonic reception data processing unit 1, and a memory control that controls reading and writing of data from and to the memory 3.
  • Unit 2 first control means and read by memory control unit 2.
  • Filter operation unit 4 (first spatial filter operation unit) that performs a filtering process on the received data, and a filter coefficient operation unit 5 that calculates the filter coefficient and supplies it to the spatial filter operation unit 4 .
  • the filter coefficient calculator 5 is configured to calculate a filter coefficient based on information including a positional relationship between a reception beam and a transmission beam.
  • the positional relationship with the transmission beam is, for example, information on whether or not the reception beam is a force generated from the same transmission beam.
  • the spatial filter operation unit 4 performs a filtering process on the plurality of reception beam data based on the filter coefficient supplied from the filter coefficient operation unit 5 in order to reduce the image quality difference between adjacent beams.
  • the output data is supplied to the scan converter 7. Note that the information on the reception beam may be directly supplied to the memory control unit 2 and the filter coefficient calculation unit 5.
  • the filter coefficient is a value that cancels out the difference in correlation between image data without making the value bilaterally symmetric.
  • a filter coefficient to be applied is reduced, and a filter applied to image data on reception beams generated from different reception beams is used. Increase the coefficient. Thereby, the image quality difference that may occur between adjacent beams is reduced.
  • the following description relates to such a filtering process.
  • FIG. 2A is a schematic diagram showing image data at an arbitrary depth when a lateral filter is used.
  • FIG. 2B is a schematic diagram showing filter coefficients corresponding to the image data of FIG. 2A.
  • FIGS. 2A and 2B show how receive beams RX1, RX2, RX3, RX4, and RX5 are generated from transmit beams TX1, ⁇ 2, ⁇ 3, ⁇ 4, and ⁇ 5, respectively.
  • Dl, D2, D3, D4, and D5 indicate image data at arbitrary depths on the reception beams RX1, RX2, RX3, RX4, and RX5, respectively.
  • a filter coefficient for generating data at an arbitrary sampling point on the reception beam RX3 is C1.
  • the coefficient of the filter corresponding to the receive beam RX2 is C21
  • the coefficient of the filter corresponding to the receive beam RX4 is C22
  • the coefficient of the filter corresponding to the receive beam RX1 is C31
  • the coefficient of the filter corresponding to the receive beam RX5. Is C32.
  • the output data Out of the filtering result is represented by the following equations (1) and (2).
  • the filter coefficient in the filtering process for reducing the image quality difference between adjacent beams is such that the value of the filter coefficient C1 at the same position as an arbitrary sampling point (Out) is obtained. It is set to the maximum, and the other filter coefficients are set to values symmetrical to C1. That is, the filter coefficients C21, C22, C31, and C32 can be expressed by the following equation (3).
  • the output data Out of the filtering result is represented by the following equations (4), (5), and (6).
  • FIGS. 3A and 3B show setting of filter coefficients when a diamond-shaped filter is used instead of the lateral filter.
  • 3A and 3B are schematic diagrams showing image data at an arbitrary depth and a corresponding filter coefficient when a diamond-shaped filter is used, respectively.
  • the filter coefficient C1 at the same position as an arbitrary sampling point (Out) is set to the maximum, and the other filter coefficients are set. Is set to a value that is symmetrical left and right and up and down with respect to the filter coefficient C1.
  • FIG. 4A is a schematic diagram showing image data at an arbitrary depth.
  • Figures 4B and 4C show the filter coefficients for generating data OutL at any sampling point on the receive beam RX1L and the arbitrary sampling on the receive beam RX1R, respectively, when using a lateral filter.
  • FIG. 9 is a schematic diagram illustrating a filter coefficient when generating point data OutR.
  • the transmit beam TX0 generates-receive beams RX0L, RX0R
  • the transmit beam TX1 generates-receive beams RX1L, RX1R
  • the transmit beam From TX2-one receive beam RX2L, RX2R is generated.
  • Dl, D2, D3, D4, D5, and D6 indicate image data at arbitrary depths on the receive beams RX0L, RX0R, RX1L, RX1R, RX2L, and RX2R, respectively.
  • the image data Dl and D2 are image data obtained from the same transmission beam TXO, they often have closer correlation values than D2 and D3, which have strong correlation. This relationship is expressed by the following equation (10).
  • D3 and D4 and D5 and D6 also have strong correlations.
  • the value of the filter coefficient is not bilaterally symmetric, and as described above, the difference between the correlation between image data obtained from the same transmission beam and the correlation between image data obtained from different transmission beams. Asymmetrically set to negate. For example, when the filter coefficient is a positive value, the filter coefficient is set so as to satisfy the following relational expression (13).
  • the filtering process for mitigating the image quality difference between the adjacent beams is optimized, and it is possible to display an ultrasonic image with good image quality in which the generation of fringes in the arrangement direction of acoustic lines is suppressed. .
  • the coefficients of the filters that generate data OutR at arbitrary sampling points on the receive beam RX1R are Cl, C21, C22, C31, and C32, and the filter coefficients are If the value is a positive value, the output data OutR of the filtering result is represented by the following equations (14), (15), and (16).
  • FIG. 5A is a schematic diagram showing image data at an arbitrary depth.
  • 5B and 5C show the filter coefficients for generating data at any sampling point on the receive beam RX1L and the arbitrary sampling on the receive beam RX1R, respectively, when a diamond-shaped filter is used.
  • FIG. 9 is a schematic diagram illustrating filter coefficients when generating point data.
  • the filter coefficients are set to be asymmetrical and asymmetrical about an arbitrary sampling point (OutL or OutR), as in the case of the lateral filter, so that the same transmission beam power can be obtained.
  • the difference between the correlation between the data and the correlation between the image data obtained with different transmitted beam powers is set to cancel.
  • the output data OutL of the filtering result uses the filter coefficient shown in FIG.5B
  • the output data OutR uses the filter coefficient shown in FIG. Expressions (17), (18), (19) and (20) are given.
  • a signal between a plurality of reception beams including a parallel reception beam obtained from the same transmission beam depends on the position of the reception beam with respect to the transmission beam.
  • the filter coefficient can be optimally controlled. As a result, It is possible to suppress the occurrence of stripes in the column direction and to display an ultrasonic image with good image quality that can be expressed in detail. Further, beam distortion can be corrected by optimizing the filter coefficient.
  • FIG. 6 is a functional block diagram showing one configuration example of the ultrasonic diagnostic apparatus according to the second embodiment of the present invention.
  • the present embodiment relates to an ultrasonic diagnostic apparatus having a two-dimensional Doppler (color Doppler) function.
  • This ultrasonic diagnostic apparatus is obtained by adding elements for a two-dimensional Doppler function to the configuration of the first embodiment.
  • the memory control unit 2 in FIG. 1 is referred to as a first memory control unit 2
  • the spatial filter operation unit 4 in FIG. 1 is referred to as a first spatial filter operation unit 4.
  • a brightness signal processing unit 10 and a two-dimensional Doppler signal processing unit 11 are provided as elements interposed between the ultrasonic reception data processing unit 1 and the spatial filter 15.
  • the spatial filter 15 includes a first memory control unit 2, a memory 3, a first spatial filter operation unit 4, and a filter coefficient operation unit 5, similar to the spatial filter 9 according to the first embodiment.
  • a two-dimensional Doppler memory 13 (second storage unit), a second memory control unit 12 (second control unit), and a second spatial filter operation unit 14 are provided.
  • the data supplied via the brightness signal processing unit 10 is processed by the first memory control unit 2, the memory 3, and the first spatial filter operation unit 4 in the same processing as in the configuration of the first embodiment. And supplied to the scan conversion unit 7.
  • the two-dimensional Doppler signal processing unit 11 performs two-dimensional Doppler processing on the reception beam data from the ultrasonic reception data processing unit 1.
  • the two-dimensional Doppler memory 13 stores the two-dimensional Doppler data output from the two-dimensional Doppler signal processing unit 11.
  • the second memory control unit 12 controls reading and writing of data from and to the two-dimensional Doppler memory 13.
  • the second spatial filter operation unit 14 performs a filtering process on the received two-dimensional Doppler data read by the second memory control unit 12 according to the filter coefficient from the filter coefficient operation unit 5, and performs a scan conversion unit. Supply output data to 7.
  • the filtering process is performed on a plurality of received two-dimensional Doppler data to reduce the image quality difference between adjacent beams.
  • a signal between a plurality of reception beams including a parallel reception beam obtained from the same transmission beam can be used.
  • the filter coefficient can be optimally controlled according to the position of the reception beam with respect to the transmission beam. As a result, it is possible to suppress the occurrence of remarkable stripes in the arrangement direction of acoustic lines in the two-dimensional Doppler, and to display a supersonic image of good image quality that can be expressed in detail.
  • beam distortion can be corrected by optimizing the filter coefficients.
  • FIG. 7 is a functional block diagram showing one configuration example of the ultrasonic diagnostic apparatus according to the third embodiment of the present invention.
  • the present embodiment relates to an ultrasonic diagnostic apparatus having a filter coefficient control function according to a reception depth.
  • This ultrasonic diagnostic apparatus is provided with a depth information generation unit 17 that supplies reception depth information to a filter coefficient calculation unit 16 in addition to the configuration of the first embodiment.
  • the filter coefficient it is possible to set the filter coefficient to be different depending on the depth of the image data. For example, as shown in FIG. 8, there is a difference in the distance between the acoustic lines 18 depending on the depth. Therefore, when the distance between the image data D on the adjacent acoustic lines 18 differs depending on the depth, the degree of the correlation different depending on the depth is determined. It is possible to mitigate. For this purpose, a shallower part is set to a stronger correlation, a filter coefficient, and a deeper part is set to a weaker correlation and a filter coefficient. As a result, it is possible to display an ultrasonic image with less lateral flow and good image quality.
  • FIG. 9 is a functional block diagram showing one configuration example of the ultrasonic diagnostic apparatus according to the fourth embodiment of the present invention.
  • the present embodiment relates to an ultrasonic diagnostic apparatus having a filter coefficient control function according to a reception beam angle.
  • This ultrasonic diagnostic apparatus is provided with a beam angle information generation unit 20 that supplies received beam angle information to a filter coefficient calculation unit 19, in addition to the configuration of the first embodiment.
  • the filter coefficient it is possible to set the filter coefficient to be different depending on the reception beam angle. For example, as shown in FIG. 10A, when transmitting and receiving an acoustic line 21 at an angle, By optimizing the filter coefficient such as weak correlation between the same depths and using a filter coefficient, it becomes possible to display an ultrasonic image of good image quality with less lateral flow. Further, as shown in FIG. 10B, when there is a difference in the deflection angle due to the acoustic line 22, the beam distortion can be corrected by optimizing the filter coefficient.
  • FIG. 11 is a functional block diagram showing one configuration example of the ultrasonic diagnostic apparatus according to the fifth embodiment of the present invention.
  • the present embodiment relates to an ultrasonic diagnostic apparatus having a filter coefficient control function according to the focal position of a transmission beam.
  • This ultrasonic diagnostic apparatus includes a transmission beam focal position information generation unit 24 that supplies focal position information of a transmission beam to a filter coefficient calculation unit 23 in addition to the configuration of the first embodiment.
  • the filter coefficient it is possible to control the filter coefficient according to the focal position of the transmission beam. For example, if the focus of the transmit beam is narrowed to the focus position F shown in Fig. 12, a filter coefficient with a weak correlation is used near the focus position F of the transmit beam, and the focus position F of the transmit beam is also far away. Then, the strength of the correlation is set to the filter coefficient. By optimizing the filter coefficients in this way, it is possible to display a uniform and good-quality ultrasound image regardless of the focal position of the transmission beam.
  • the ultrasonic diagnostic apparatus of the present invention suppresses the generation of remarkable fringes in the direction of arrangement of acoustic lines in a two-dimensional Doppler for a signal between a plurality of reception beams having the same transmission beam power. This is useful for displaying an ultrasonic image with good image quality that can be expressed up to the maximum.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Processing (AREA)

Abstract

 超音波受信信号から形成した受信ビームをディジタル変換して得られる受信ビームデータを格納する第1の記憶手段(3)と、第1の記憶手段に対するデータの読み出しおよび書き込みを制御する第1の制御手段(2)と、受信ビームに関する送信ビームとの位置関係を含む情報に基づきフィルタ係数を演算するフィルタ係数演算部(5)と、同一の送信ビームから得られる並列受信ビームを含む複数本の受信ビームデータに対して、フィルタ係数に基づき、隣り合うビーム間の画質差を緩和するためのフィルタリング処理を施す第1の空間フィルタ演算部(4)とを備え、第1の空間フィルタ演算部から出力された画像データを表示用モニタ(8)の走査に変換して画像を表示する。並列受信機能を有する超音波診断装置において、音響線の配列方向における縞の発生を抑制し、細部まで表現できる画質の良好な超音波画像を表示することができる。

Description

明 細 書
超音波診断装置
技術分野
[0001] 本発明は、同一の送信ビーム力 複数の受信ビームを生成する並列受信機能を有 する超音波診断装置に関する。
背景技術
[0002] 超音波診断装置は、図 13に示すように、被検体に対して超音波の送受信を行う超 音波探触子 30により超音波ビームを走査して得られる 2次元情報にもとづき、超音波 画像を表示するものである。超音波探触子 30により受信された信号は、信号処理部 31を介して空間フィルタ処理回路 32に供給される。空間フィルタ処理回路 32は、画 像表示変換部 33の前段に配置され、受信信号に対してノイズを低減するための空 間フィルタで構成されている。制御回路 34は、空間フィルタ処理回路 32のフィルタ係 数を、音響データの振動子面近傍から、遠点方向にかけて変化させるように制御す る。つまり、従来の超音波診断装置では、座標変換前の音響データに対し空間フィ ルタ処理を行う際に、音響線データ上の距離に応じてフィルタ係数を変化させて!/、る (例えば、特許文献 1を参照)。
特許文献 1 :特開 2001— 340338号公報 (第 1図)
発明の開示
発明が解決しょうとする課題
[0003] 上記従来の超音波診断装置においては、並列受信機能における、同一の送信ビ ーム力 得られる複数の受信ビーム間の信号差は、異なる送信ビームから得られる 受信ビーム間の信号差と比較すると小さい。このため、隣り合うビーム間の画質差を 緩和するためのフィルタ係数が固定されていることにより、以下のような問題を生じて いた。すなわち、同一の送信ビーム力も得られた受信ビーム信号間では、フィルタが 効きすぎて、スムージング効果により画像が均一化して、画像の細部が表示できなく なる。一方、異なる送信ビーム力 得られる受信ビーム信号間では、フィルタが効き 難く、エッジ強調効果により、受信ビーム間の境界が表示されてしまう。 [0004] 特に、二次元ドプラでは、同一音響線上で送受信を 10回程度繰り返すことによって 時間変化を検出するため、異なる送信ビーム間での時間経過が大きい。そのため、 受信ビームが同一送信ビームから得られたか、異なる送信ビームカゝら得られたかによ る変化が著しぐ画像の均一化や受信ビーム間の境界の表示が顕著になり、音響線 の配列方向に縞が発生し易い。また、二次元ドプラにおいて、時間変化を色付表示 する場合は、変化がゼロで色表示をしな 、時と変化があって色付けされる境界が歴 然とするため、データの有無による音響線方向の縞が発生し易 、。
[0005] 本発明は、上記従来の問題を解決するためになされたもので、その目的は、同一 の送信ビームから得られた複数の受信ビーム間での信号に対するフィルタリング処理 を最適化することで、音響線の配列方向における縞の発生を抑制し、細部まで表現 できる画質の良好な超音波画像を表示することが可能な超音波診断装置を提供す ることにめる。
課題を解決するための手段
[0006] 前記の目的を達成するため、本発明に係る超音波診断装置は、超音波受信信号 から形成した受信ビームをディジタル変換して得られる受信ビームデータを格納する 第 1の記憶手段 (メモリ)と、第 1の記憶手段に対するデータの読み出しおよび書き込 みを制御する第 1の制御手段 (メモリ制御部、第 1のメモリ制御部)と、受信ビームに関 する送信ビームとの位置関係を含む情報に基づきフィルタ係数を演算するフィルタ係 数演算部と、同一の送信ビームから得られる並列受信ビームを含む複数本の受信ビ ームデータに対して、フィルタ係数に基づき、隣り合うビーム間の画質差を緩和する ためのフィルタリング処理を施す第 1の空間フィルタ演算部とを備え、第 1の空間フィ ルタ演算部から出力された画像データを表示用モニタの走査に変換して画像を表示 する構成を有する。
[0007] この構成により、同一の送信ビームから得られる並列受信ビームを含む複数の受信 ビーム間の信号に対して、送信ビームに対する受信ビームの位置に応じてフィルタ 係数を最適に制御することができ、それにより、音響線の配列方向における縞の発生 を抑制し、細部まで表現できる画質の良好な超音波画像を表示することができる。
[0008] また、上記構成の超音波診断装置にお!、て更に、超音波受信データ処理部からの 受信ビームデータに対して二次元ドプラ処理を施す二次元ドプラ信号処理部と、二 次元ドプラ信号処理部から出力される二次元ドプラデータを格納する第 2の記憶手 段(二次元ドプラメモリ)と、第 2の記憶手段に対するデータの読み出しおよび書き込 みを制御する第 2の制御手段 (第 2のメモリ制御部)と、フィルタ係数演算部から供給 されるフィルタ係数に従って、同一の送信ビームから得られる並列受信ビームを含む 複数本の受信二次元ドプラデータに対して、隣り合うビーム間の画質差を緩和するた めのフィルタリング処理を施す第 2の空間フィルタ演算部とを備えた構成とすることが できる。
[0009] この構成により、二次元ドプラ (カラードプラ)機能を有する超音波診断装置におい て、同一の送信ビーム力 得られる並列受信ビームを含む複数の受信ビーム間の信 号に対して、送信ビームに対する受信ビームの位置に応じてフィルタ係数を最適に 制御することができ、それにより、二次元ドプラにおいて顕著な音響線の配列方向に おける縞の発生を抑制し、細部まで表現できる画質の良好な超音波画像を表示する ことができる。
[0010] また、上記構成の超音波診断装置にお!、て、フィルタ係数演算部は、フィルタ係数 を受信深度に応じて制御可能である構成とすることができる。
[0011] この構成により、深度によって音響線間隔に差の生じる場合、浅い部分ではより相 関の強いフィルタ係数、深い部分では相関の弱いフィルタ係数とすることができ、そ れにより、横流れの少ない画質の良好な超音波画像を表示することが可能になる。
[0012] また、上記構成の超音波診断装置にお!、て、フィルタ係数演算部は、フィルタ係数 を受信ビーム角度に応じて制御可能である構成とすることができる。
[0013] この構成により、音響線に角度をつけて送受信する場合、同一深度間での相関の 弱いフィルタ係数とするなどフィルタ係数を最適化することができ、それにより、横流 れの少ない画質の良好な超音波画像を表示することが可能になる。また、音響線に よって偏向角度に差がある場合、フィルタ係数を最適化することで、ビームの歪みを 補正することが可能になる。
[0014] また、上記構成の超音波診断装置にお!、て、フィルタ係数演算部は、フィルタ係数 を送信ビームの焦点位置に応じて制御可能である構成とすることができる。 [0015] この構成により、送信ビームの焦点位置近くでは相関の弱いフィルタ係数とし、送信 ビームの焦点位置力も離れると相関の強 、フィルタ係数とするなど、フィルタ係数を 最適化することができ、それにより、送信ビームの焦点位置によらず均質な画質の良 好な超音波画像を表示することが可能になる。
発明の効果
[0016] 本発明によれば、同一の送信ビーム力も得られた複数の受信ビーム間での信号に 対するフィルタリング処理を最適化することで、音響線の配列方向における縞の発生 を抑制し、細部まで表現できる画質の良好な超音波画像を表示することができる超 音波診断装置を提供することが可能になる、という格別な効果を奏する。
図面の簡単な説明
[0017] [図 1]図 1は、本発明の第 1の実施の形態に係る超音波診断装置の一構成例を示す 機能ブロック図である。
[図 2A]図 2Aは、第 1の実施の形態における、一送信ビームから一受信ビームを生成 する場合で、ラテラルフィルタを用いた場合の任意の深さの画像データを示す模式 図である。
[図 2B]図 2Bは、図 2Aの画像データに対応するフィルタ係数を示す模式図である。
[図 3A]図 3Aは、第 1の実施の形態における、一送信ビームから一受信ビームを生成 する場合で、ダイアモンド形状のフィルタを用いた場合の任意の深さの画像データを 示す模式図である。
[図 3B]図 3Bは、図 3Aの画像データに対応するフィルタ係数を示す模式図である。
[図 4A]図 4Aは、第 1の実施の形態における、一送信ビームから二本の受信ビームを 生成する場合で、ラテラルフィルタを用いた場合の任意の深さの画像データを示す 模式図である。
[図 4B]図 4Bは、図 4Aの画像データに対して、受信ビーム RX1L上の任意のサンプリ ング点のデータ OutLを生成する場合のフィルタ係数を示す模式図である。
[図 4C]図 4Cは、図 4Aの画像データに対して、受信ビーム RX1R上の任意のサンプ リング点のデータ OutRを生成する場合のフィルタ係数を示す模式図である。
[図 5A]図 5Aは、第 1の実施の形態における、一送信ビームから二本の受信ビームを 生成する場合で、ダイアモンド形状のフィルタを用いた場合の任意の深さの画像デ ータを示す模式図である。
[図 5B]図 5Bは、図 5Aの画像データに対して、受信ビーム RX1L上の任意のサンプリ ング点のデータ OutLを生成する場合のフィルタ係数を示す模式図である。
[図 5C]図 5Cは、図 5Aの画像データに対して、受信ビーム RX1R上の任意のサンプ リング点のデータ OutRを生成する場合のフィルタ係数を示す模式図である。
[図 6]図 6は、本発明の第 2の実施の形態に係る超音波診断装置の一構成例を示す 機能ブロック図である。
[図 7]図 7は、本発明の第 3の実施の形態に係る超音波診断装置の一構成例を示す 機能ブロック図である。
[図 8]図 8は、第 3の実施の形態における、深度によって音響線間隔に差の生じる場 合のフィルタ係数制御機能を説明するための模式図である。
[図 9]図 9は、本発明の第 4の実施の形態に係る超音波診断装置の一構成例を示す 機能ブロック図である。
[図 10A]図 10Aは、第 4の実施の形態における、音響線に角度をつけて送受信する 場合のフィルタ係数制御機能を説明するための模式図である。
[図 10B]図 10Bは、第 4の実施の形態における、音響線によって偏向角度に差がある 場合のフィルタ係数制御機能を説明するための模式図である。
[図 11]図 11は、本発明の第 5の実施の形態に係る超音波診断装置の一構成例を示 す機能ブロック図である。
[図 12]図 12は、第 5の実施の形態における、送信ビームの焦点位置によるフィルタ係 数制御機能を説明するための模式図である。
[図 13]図 13は、従来の超音波診断装置の一構成例を示す機能ブロック図である。 符号の説明
1 超音波受信データ処理部
2 メモリ(第 1の記憶手段)
3 メモリ制御部、第 1のメモリ制御部 (第 1の制御手段)
4 空間フィルタ演算部、第 1の空間フィルタ演算部 5、 16、 19、 23 フィルタ係数演算部
6 音響ビーム制御部
7 走査変換部
8 モニタ (表示手段)
9、 15、 16 空間フィルタ
10 ブライトネス信号処理部
11 二次元ドプラ信号処理部
12 二次元ドプラメモリ(第 2の記憶手段)
13 第 2のメモリ制御部 (第 2の制御手段)
14 第 2の空間フィルタ演算部
17 深度情報生成部
18、 21、 22 音響線
20 ビーム角度情報生成部
24 送信ビーム焦点位置情報生成部
発明を実施するための最良の形態
[0019] 以下、本発明の好適な実施の形態について、図面を参照しながら詳細に説明する
[0020] (第 1の実施の形態)
図 1は、本発明の第 1の実施の形態に係る超音波診断装置の一構成例を示す機能 ブロック図である。
[0021] この超音波診断装置は、超音波受信信号の受信ビームを形成し、受信ビームを受 信ビームデータにディジタル変換する超音波受信データ処理部 1と、空間フィルタ 9と 、受信ビームに関する情報を超音波受信データ処理部 1に供給する音響ビーム制御 部 6と、空間フィルタ 9からの出力データを画像データに変換する走査変換部 7と、画 像データを表示するモニタ 8 (表示手段)とから構成される。
[0022] 空間フィルタ 9は、超音波受信データ処理部 1から供給される受信ビームデータを 格納するメモリ 3 (第 1の記憶手段)と、メモリ 3に対するデータの読み出しおよび書き 込みを制御するメモリ制御部 2 (第 1の制御手段)と、メモリ制御部 2により読み出され たデータに対してフィルタリング処理を施す空間フィルタ演算部 4 (第 1の空間フィル タ演算部)と、フィルタ係数を演算して空間フィルタ演算部 4に供給するフィルタ係数 演算部 5とから構成される。
[0023] フィルタ係数演算部 5は、受信ビームに関する送信ビームとの位置関係を含む情報 に基づきフィルタ係数を演算するように構成される。送信ビームとの位置関係とは、例 えば、受信ビームが同一の送信ビームから生成されたものである力否かに関する情 報である。空間フィルタ演算部 4は、フィルタ係数演算部 5から供給されるフィルタ係 数に基づき、複数本の受信ビームデータに対し、隣り合うビーム間の画質差を緩和 するためのフィルタリング処理を施して、その出力データを走査変換部 7に供給する 。なお、受信ビームに関する情報は、メモリ制御部 2およびフィルタ係数演算部 5に直 接供給されてもよい。
[0024] 次に、以上のように構成された超音波診断装置の動作にっ 、て、図 2A、 2B、図 3 A、 3B、図 4A— 4C、および図 5A— 5Cを参照して説明する。特に本発明の主要部 である空間フィルタ 9に重点を置いて説明する。
[0025] 本実施の形態においては、同一の送信ビーム力 複数の受信ビームを生成する場 合のフィルタ係数は、値を左右対称にはせずに、画像データ間の相関の差異を打ち 消す値に設定する。つまり、同一送信ビームから生成される複数の受信ビーム上の 画像データは互いに相関が強いため、掛けられるフィルタ係数は小さくし、異なる受 信ビームから生成される受信ビーム上の画像データに掛けられるフィルタ係数は大き くする。それにより、隣り合うビーム間で発生する可能性のある画質差を緩和する。以 下の説明は、このようなフィルタリング処理に関するものである。
[0026] まず、一送信ビーム力 一受信ビームを生成する場合の、隣り合うビーム間の画質 差を緩和するためのフィルタリング処理について、ラテラルフィルタを例に挙げ、図 2 Aおよび図 2Bを参照して説明する。この処理は、従来技術と同様である。図 2Aは、 ラテラルフィルタを用いた場合における、任意の深さの画像データを示す模式図であ る。図 2Bは、図 2Aの画像データに対応するフィルタ係数を示す模式図である。
[0027] 図 2Aおよび図 2Bには、送信ビーム TX1、 ΤΧ2、 ΤΧ3、 ΤΧ4および ΤΧ5から、受 信ビーム RX1、 RX2、 RX3、 RX4および RX5がそれぞれ生成される様子が示される 。図 2Aにおいて、 Dl、 D2、 D3、 D4および D5はそれぞれ受信ビーム RX1、 RX2、 RX3、 RX4および RX5上の任意の深さの画像データを示す。ここでは、ラテラルフィ ルタが、 5つの受信ビームから 1つのビームを生成する場合を例として説明する。図 2 Bに示すように、受信ビーム RX3上の任意のサンプリング点のデータを生成するフィ ルタの係数を C1とする。同様に、受信ビーム RX2に対応するフィルタの係数を C21 、受信ビーム RX4に対応するフィルタの係数を C22、受信ビーム RX1に対応するフ ィルタの係数を C31、および受信ビーム RX5に対応するフィルタの係数を C32とする 。この場合、フィルタリング結果の出力データ Outは、以下の式(1)、式(2)で表され る。
[0028] Out=(DlXC31 + D2XC21 + D3XCl + D4XC22
+ D5XC32) …式(1)
C31 + C21 + C1 + C22 + C32=1 …式(2)
[0029] 隣り合うビーム間の画質差を緩和するためのフィルタリング処理におけるフィルタ係 数は、一般に図 2Bに示すように、任意のサンプリング点(Out)と同じ位置でのフィル タ係数 C1の値が最大に設定され、他のフィルタ係数は C1に対して左右対称な値に 設定される。つまり、フィルタ係数 C21、 C22、 C31、 C32は、以下の関係を有する式 (3)で表せる。
[0030] C21 = C22 = C2、 C31 = C32 = C3 …式(3)
[0031] ここで、フィルタ係数を正の値とした場合、フィルタリング結果の出力データ Outは、 以下の式 (4)、式(5)、式 (6)で表される。
[0032] Out=(DlXC3 + D2XC2 + D3XCl + D4XC2 + D5XC3)
…式 (4)
C3 + C2 + C1 + C2 + C3 = 1 …式(5)
C1>C2>C3 …式(6)
[0033] ラテラルフィルタに代えて、ダイアモンド形状のフィルタを用いた場合におけるフィル タ係数の設定について、図 3Aおよび図 3Bに示す。図 3Aおよび図 3Bは、それぞれ 、ダイアモンド形状のフィルタを用いた場合における任意の深さの画像データ、およ び対応するフィルタ係数を示す模式図である。 [0034] ダイアモンド形状のフィルタにおいても、フィルタ係数は、前記ラテラルフィルタと同 様に、任意のサンプリング点(Out)と同じ位置でのフィルタ係数 C1の値が最大に設 定され、他のフィルタ係数はフィルタ係数 C1に対して左右対称および上下対称な値 に設定される。ここで、受信ビーム RX3上の任意のサンプリング点のデータを生成す るものとし、フィルタ係数を正の値とした場合、フィルタリング結果の出力データ Out は、以下の式(7)、式 (8)、式(9)で表される。
[0035] Out=(DlXC3 + D2XC2 + D3XCl + D4XC2 + D5XC3)
+ (D2aXC4 + D3aXC5 + D4aXC4)
+ (D21XC4 + D31XC5 + D41XC4)
+ D3bXC6 + D32XC6 …式(7)
C1 + 2XC2 + 2XC3+4XC4 + 2XC5 + 2XC6 = 1
…式 (8)
C1>C2>C3力つ C1>C5>C6力つ C5>C4力つ C2>C4
…式 (9)
[0036] 次に、同一の送信ビーム力も複数の受信ビームを生成する場合の、隣り合うビーム 間の画質差を緩和するためのフィルタリング処理について、ラテラルフィルタにより二 本の受信ビームを生成する場合を例に挙げ、図 4A— 4Cを参照して説明する。この 処理は、本実施の形態に特有の条件で行われる。
[0037] 図 4Aは、任意の深さの画像データを示す模式図である。図 4Bおよび図 4Cは、そ れぞれ、ラテラルフィルタを用いた場合における、受信ビーム RX1L上の任意のサン プリング点のデータ OutLを生成する場合のフィルタ係数、および受信ビーム RX1R 上の任意のサンプリング点のデータ OutRを生成する場合のフィルタ係数を示す模 式図である。
[0038] 図 4A— 4Cに示されるように、送信ビーム TX0から-本の受信ビーム RX0L、 RX0 Rが生成され、送信ビーム TX1から-本の受信ビーム RX1L、 RX1Rが生成され、ま た送信ビーム TX2から-本の受信ビーム RX2L、RX2Rが生成される。
[0039] 図 4Aにおいて、 Dl、 D2、 D3、 D4、 D5および D6はそれぞれ、受信ビーム RX0L 、 RX0R、 RX1L、 RX1R、 RX2Lおよび RX2R上の任意の深さの画像データを示す 。画像データ Dlと D2は、同一の送信ビーム TXOから得られた画像データであるた め相関が強ぐ D2と D3よりも近い値となることが多い。この関係は以下の式(10)で 表される。また、 D3と D4、 D5と D6も同様に相関が強い。
[0040] I D1-D2 I <|D3-D2 | …式(10)
[0041] 5つの受信ビームから 1つのビームを生成するものとし、図 4Bに示されるように、受 信ビーム RX1L上の任意のサンプリング点のデータ OutLを生成するフィルタの係数 を Cl、 C21、 C22、 C31および C32とした場合、フィルタリング結果の出力データ Ou tLは、以下の式(11)、式(12)で表される。
[0042] OutL=(DlXC31 + D2XC21+D3XCl + D4XC22
+ D5XC32) …式(11)
C31 + C21 + C1 + C22 + C32=1 …式(12)
[0043] ここで、フィルタ係数の値は左右対称ではなぐ上述のように、同一送信ビームから 得られる画像データ間の相関と、異なる送信ビームから得られる画像データ間の相 関の間の差異を打ち消すように非対称に設定される。例えば、フィルタ係数を正の値 とした場合、フィルタ係数は、以下の関係式(13)を満たすように設定される。
[0044] C22<C21 …式(13)
[0045] この結果、隣り合うビーム間の画質差を緩和するためのフィルタリング処理が最適 化され、音響線の配列方向における縞の発生を抑制した画質の良好な超音波画像 を表示することができる。
[0046] 同様に、図 4Cに示されるように、受信ビーム RX1R上の任意のサンプリング点のデ ータ OutRを生成するフィルタの係数を Cl、 C21、 C22、 C31および C32とし、フィル タ係数を正の値とした場合、フィルタリング結果の出力データ OutRは、以下の式(14 )、式(15)、式(16)で表される。
[0047] OutR=(D2XC32 + D3XC22 + D4XCl + D5XC21
+ D6XC31) …式(14)
C32 + C22 + C1 + C21 + C31 = 1 …式(15)
C22<C21 …式(16)
[0048] ラテラルフィルタに代えて、ダイアモンド形状のフィルタを用いた場合におけるフィル タ係数の設定について、図 5A、図 5Bおよび図 5Cに示す。図 5Aは、任意の深さの 画像データを示す模式図である。図 5Bおよび図 5Cは、それぞれ、ダイアモンド形状 のフィルタを用いた場合における、受信ビーム RX1L上の任意のサンプリング点のデ ータを生成する場合のフィルタ係数、および受信ビーム RX1R上の任意のサンプリン グ点のデータを生成する場合のフィルタ係数を示す模式図である。
[0049] ダイアモンド形状のフィルタにおいても、フィルタ係数は、前記ラテラルフィルタと同 様に、任意のサンプリング点(OutLまたは OutR)を中心に、左右非対称および上下 非対称にして、同一送信ビーム力も得られる画像データ間の相関と、異なる送信ビー ム力 得られる画像データ間の相関の間の差異を打ち消すように設定される。例えば 、フィルタ係数を正の値とした場合、フィルタリング結果の出力データ OutLは、図 5B に示されるフィルタ係数を用いて、出力データ OutRは、図 5Cに示されるフィルタ係 数を用いて、以下の式(17)、式(18)、式(19)および式(20)で表される。
[0050] OutL=(DlXC31 + D2XC21+D3XCl + D4XC22
+ D5XC32)
+ (D2a XC41 + D3aXC5 + D4a X C42)
+ (D21XC41 + D31XC5 + D41XC42)
+ D3bXC6 + D32XC6 …式(17)
OutR=(D2XC32 + D3XC22 + D4XCl + D5XC21
+ D6XC31)
+ (D3aXC42 + D4aXC5 + D5aXC41)
+ (D31XC42 + D41XC5 + D51XC41)
+ D4bXC6 + D42XC6 …式(18)
C1 + C21 + C22 + C31 + C32 + 2XC41 + 2XC42 + 2XC5
+ 2XC6 = 1 …式(19)
C21>C22力つ C41>C42 …式(20)
[0051] 以上のように、本実施の形態によれば、同一の送信ビームから得られる並列受信ビ ームを含む複数の受信ビーム間の信号に対して、送信ビームに対する受信ビームの 位置に応じてフィルタ係数を最適に制御することができる。それにより、音響線の配 列方向における縞の発生を抑制し、細部まで表現できる画質の良好な超音波画像を 表示することが可能となる。またフィルタ係数の最適化により、ビームの歪みも補正す ることが可能となる。
[0052] (第 2の実施の形態)
図 6は、本発明の第 2の実施の形態に係る超音波診断装置の一構成例を示す機能 ブロック図である。本実施の形態は、二次元ドプラ (カラードプラ)機能を有する超音 波診断装置に関する。この超音波診断装置は、第 1の実施の形態の構成に対して、 二次元ドプラ機能ための要素をカ卩えたものである。なお、以下の説明では、図 1のメ モリ制御部 2を第 1のメモリ制御部 2と称し、図 1の空間フィルタ演算部 4を第 1の空間 フィルタ演算部 4と称する。
[0053] 本実施の形態により追加される要素に関して、以下に説明する。まず、超音波受信 データ処理部 1と空間フィルタ 15との間に介在する要素として、ブライトネス信号処理 部 10および二次元ドプラ信号処理部 11が設けられる。空間フィルタ 15は、第 1の実 施の形態の空間フィルタ 9と同様の、第 1のメモリ制御部 2、メモリ 3、第 1の空間フィル タ演算部 4、およびフィルタ係数演算部 5からなる構成に加えて、二次元ドプラメモリ 1 3 (第 2の記憶手段)と、第 2のメモリ制御部 12 (第 2の制御手段)と、第 2の空間フィル タ演算部 14とが設けられる。
[0054] ブライトネス信号処理部 10を介して供給されるデータは、第 1のメモリ制御部 2、メモ リ 3、第 1の空間フィルタ演算部 4により第 1の実施の形態の構成と同様の処理を施さ れて、走査変換部 7に供給される。
[0055] 二次元ドプラ信号処理部 11は、超音波受信データ処理部 1からの受信ビームデー タに対して二次元ドプラ処理を施す。二次元ドプラメモリ 13は、二次元ドプラ信号処 理部 11から出力される二次元ドプラデータを格納する。第 2のメモリ制御部 12は、二 次元ドプラメモリ 13に対するデータの読み出しおよび書き込みを制御する。第 2の空 間フィルタ演算部 14は、フィルタ係数演算部 5からのフィルタ係数に従って、第 2のメ モリ制御部 12により読み出される受信二次元ドプラデータに対してフィルタリング処 理を施し、走査変換部 7に出力データを供給する。フィルタリング処理は、複数本の 受信二次元ドプラデータに対して隣り合うビーム間の画質差を緩和するように施され る。
[0056] このように構成することにより、二次元ドプラ (カラードプラ)機能を有する超音波診 断装置において、同一の送信ビームから得られる並列受信ビームを含む複数の受信 ビーム間の信号に対して、送信ビームに対する受信ビームの位置に応じてフィルタ 係数を最適に制御することができる。それにより、二次元ドプラにおいて顕著な音響 線の配列方向における縞の発生を抑制し、細部まで表現できる画質の良好な超音 波画像を表示すること可能になる。また、フィルタ係数の最適化により、ビームの歪み を補正することも可能になる。
[0057] (第 3の実施の形態)
図 7は、本発明の第 3の実施の形態に係る超音波診断装置の一構成例を示す機能 ブロック図である。本実施の形態は、受信深度に応じたフィルタ係数制御機能を有す る超音波診断装置に関する。この超音波診断装置は、第 1の実施の形態の構成に加 えて、フィルタ係数演算部 16に受信深度情報を供給する深度情報生成部 17を設け たものである。
[0058] この構成によれば、画像データの深度に応じてフィルタ係数を異ならせる設定が可 能である。例えば、図 8に示すように深度によって音響線 18の間隔に差が生じ、従つ て、隣り合う音響線 18上の画像データ Dの間隔が深度によって異なる場合に、深度 によって異なる相関の程度を緩和することが可能である。そのためには、浅い部分で はより相関の強 、フィルタ係数、深 、部分では相関の弱 、フィルタ係数に設定する。 その結果、横流れの少な 、画質の良好な超音波画像を表示することが可能になる。
[0059] (第 4の実施の形態)
図 9は、本発明の第 4の実施の形態に係る超音波診断装置の一構成例を示す機能 ブロック図である。本実施の形態は、受信ビーム角度に応じたフィルタ係数制御機能 を有する超音波診断装置に関する。この超音波診断装置は、第 1の実施の形態の構 成に加えて、フィルタ係数演算部 19に受信ビーム角度情報を供給するビーム角度情 報生成部 20を設けたものである。
[0060] この構成によれば、受信ビーム角度に応じてフィルタ係数を異ならせる設定が可能 である。例えば、図 10Aに示すように、音響線 21に角度をつけて送受信する場合、 同一深度間での相関の弱 、フィルタ係数とするなど、フィルタ係数を最適化すること で、横流れの少ない画質の良好な超音波画像を表示することが可能になる。また、 図 10Bに示すように、音響線 22によって偏向角度に差がある場合、フィルタ係数を 最適化することで、ビームの歪みを補正することが可能になる。
[0061] (第 5の実施の形態)
図 11は、本発明の第 5の実施の形態に係る超音波診断装置の一構成例を示す機 能ブロック図である。本実施の形態は、送信ビームの焦点位置に応じたフィルタ係数 制御機能を有する超音波診断装置に関する。この超音波診断装置は、第 1の実施の 形態の構成に加えて、フィルタ係数演算部 23に送信ビームの焦点位置情報を供給 する送信ビーム焦点位置情報生成部 24を設けたものである。
[0062] この構成によれば、送信ビームの焦点位置に応じてフィルタ係数を制御することが 可能である。例えば、図 12に示す焦点位置 Fに、送信ビームの焦点が絞られている 場合に、送信ビームの焦点位置 Fの近くでは相関の弱いフィルタ係数とし、送信ビー ムの焦点位置 F力も離れた位置では相関の強 、フィルタ係数に設定する。そのように フィルタ係数を最適化することで、送信ビームの焦点位置によらず均質で画質の良 好な超音波画像を表示することが可能になる。
産業上の利用可能性
[0063] 本発明の超音波診断装置は、同一の送信ビーム力 得られる複数の受信ビーム間 の信号に対して、二次元ドプラにおいて顕著な音響線の配列方向における縞の発生 を抑制し、細部まで表現できる画質の良好な超音波画像を表示するために有用であ る。

Claims

請求の範囲
[1] 超音波受信信号力 形成した受信ビームをディジタル変換して得られる受信ビーム データを格納する第 1の記憶手段と、
前記第 1の記憶手段に対するデータの読み出しおよび書き込みを制御する第 1の 制御手段と、
前記受信ビームに関する送信ビームとの位置関係を含む情報に基づきフィルタ係 数を演算するフィルタ係数演算部と、
同一の送信ビーム力 得られる並列受信ビームを含む複数本の前記受信ビームデ ータに対して、前記フィルタ係数に基づき、隣り合うビーム間の画質差を緩和するた めのフィルタリング処理を施す第 1の空間フィルタ演算部とを備え、
前記第 1の空間フィルタ演算部力 出力された画像データを表示用モニタの走査 に変換して画像を表示することを特徴とする超音波診断装置。
[2] 前記超音波受信データ処理部力 の受信ビームデータに対して二次元ドプラ処理 を施す二次元ドプラ信号処理部と、
前記二次元ドプラ信号処理部から出力される二次元ドプラデータを格納する第 2の 記憶手段と、
前記第 2の記憶手段に対するデータの読み出しおよび書き込みを制御する第 2の 制御手段と、
前記フィルタ係数演算部力 供給されるフィルタ係数に従って、同一の送信ビーム 力 得られる並列受信ビームを含む複数本の前記受信二次元ドプラデータに対して 、隣り合うビーム間の画質差を緩和するためのフィルタリング処理を施す第 2の空間フ ィルタ演算部とを更に備えた請求項 1記載の超音波診断装置。
[3] 前記フィルタ係数演算部は、前記フィルタ係数を受信深度に応じて制御可能である 請求項 1または 2記載の超音波診断装置。
[4] 前記フィルタ係数演算部は、前記フィルタ係数を受信ビーム角度に応じて制御可 能である請求項 1または 2記載の超音波診断装置。
[5] 前記フィルタ係数演算部は、前記フィルタ係数を送信ビームの焦点位置に応じて 制御可能である請求項 1または 2記載の超音波診断装置。
PCT/JP2005/000074 2004-01-08 2005-01-06 超音波診断装置 Ceased WO2005065547A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US10/585,388 US8246544B2 (en) 2004-01-08 2005-01-06 Ultrasonic diagnosis apparatus
EP05703342A EP1709908A4 (en) 2004-01-08 2005-01-06 ULTRASONIC DIAGNOSTIC APPARATUS
JP2005516882A JP4290699B2 (ja) 2004-01-08 2005-01-06 超音波診断装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2004003310 2004-01-08
JP2004-003310 2004-01-08

Publications (1)

Publication Number Publication Date
WO2005065547A1 true WO2005065547A1 (ja) 2005-07-21

Family

ID=34747072

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2005/000074 Ceased WO2005065547A1 (ja) 2004-01-08 2005-01-06 超音波診断装置

Country Status (5)

Country Link
US (1) US8246544B2 (ja)
EP (1) EP1709908A4 (ja)
JP (1) JP4290699B2 (ja)
CN (1) CN100539950C (ja)
WO (1) WO2005065547A1 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101919706A (zh) * 2009-06-12 2010-12-22 深圳迈瑞生物医疗电子股份有限公司 抽取滤波的方法及装置
JP2016135187A (ja) * 2015-01-23 2016-07-28 コニカミノルタ株式会社 超音波画像診断装置、超音波画像生成方法及びプログラム
JP6041957B1 (ja) * 2015-10-01 2016-12-14 株式会社日立製作所 超音波診断装置
JP2020069304A (ja) * 2018-11-02 2020-05-07 コニカミノルタ株式会社 超音波診断装置、超音波診断装置の制御方法、及び、超音波診断装置の制御プログラム
US12350100B2 (en) 2020-11-16 2025-07-08 Ablic Inc. Ultrasonic device and control method of ultrasonic device

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6685645B1 (en) 2001-10-20 2004-02-03 Zonare Medical Systems, Inc. Broad-beam imaging
CN101645167B (zh) * 2008-08-06 2012-08-15 深圳迈瑞生物医疗电子股份有限公司 帧相关系数控制及帧相关处理方法、装置及超声成像系统
JP5683213B2 (ja) 2009-11-17 2015-03-11 キヤノン株式会社 画像形成装置及び画像形成方法
KR101888649B1 (ko) * 2011-11-17 2018-08-16 삼성전자주식회사 빔포밍 방법, 이를 수행하는 장치 및 의료영상시스템
KR102185415B1 (ko) * 2013-01-11 2020-12-02 삼성전자주식회사 빔 포밍 모듈, 상기 빔 포밍 모듈을 이용하는 초음파 이미징 장치, 상기 빔 포밍 모듈을 이용한 빔 포밍 방법 및 상기 빔 포밍 모듈을 이용한 초음파 이미징 장치의 제어 방법

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0614929A (ja) * 1992-04-17 1994-01-25 Hewlett Packard Co <Hp> 音響イメージング装置
JPH078492A (ja) * 1993-06-28 1995-01-13 Toshiba Corp 超音波診断装置
JPH10118063A (ja) * 1996-10-24 1998-05-12 Aloka Co Ltd 超音波診断装置
JPH10328185A (ja) * 1997-06-04 1998-12-15 Aloka Co Ltd 超音波診断装置および診断用超音波送受信方法
JP2001340338A (ja) * 2000-06-05 2001-12-11 Olympus Optical Co Ltd 超音波診断装置

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5779640A (en) * 1996-04-26 1998-07-14 Acuson Corporation Ultrasound multi-beam distortion correction system and method
US6416475B1 (en) * 1996-06-28 2002-07-09 Sonosite, Inc. Ultrasonic signal processor for a hand held ultrasonic diagnostic instrument
US6283919B1 (en) 1996-11-26 2001-09-04 Atl Ultrasound Ultrasonic diagnostic imaging with blended tissue harmonic signals
US6074347A (en) 1996-12-04 2000-06-13 Acuson Corporation Method and apparatus for controlling acoustic signal bandwidth in an ultrasonic diagnostic imaging system
US6193663B1 (en) * 1997-12-18 2001-02-27 Acuson Corporation Diagnostic ultrasound imaging method and system with improved frame rate
US6511426B1 (en) * 1998-06-02 2003-01-28 Acuson Corporation Medical diagnostic ultrasound system and method for versatile processing
US6142942A (en) 1999-03-22 2000-11-07 Agilent Technologies, Inc. Ultrasound imaging system and method employing an adaptive filter

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0614929A (ja) * 1992-04-17 1994-01-25 Hewlett Packard Co <Hp> 音響イメージング装置
JPH078492A (ja) * 1993-06-28 1995-01-13 Toshiba Corp 超音波診断装置
JPH10118063A (ja) * 1996-10-24 1998-05-12 Aloka Co Ltd 超音波診断装置
JPH10328185A (ja) * 1997-06-04 1998-12-15 Aloka Co Ltd 超音波診断装置および診断用超音波送受信方法
JP2001340338A (ja) * 2000-06-05 2001-12-11 Olympus Optical Co Ltd 超音波診断装置

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101919706A (zh) * 2009-06-12 2010-12-22 深圳迈瑞生物医疗电子股份有限公司 抽取滤波的方法及装置
JP2016135187A (ja) * 2015-01-23 2016-07-28 コニカミノルタ株式会社 超音波画像診断装置、超音波画像生成方法及びプログラム
JP6041957B1 (ja) * 2015-10-01 2016-12-14 株式会社日立製作所 超音波診断装置
WO2017056898A1 (ja) * 2015-10-01 2017-04-06 株式会社日立製作所 超音波診断装置及び受信データ処理方法
JP2020069304A (ja) * 2018-11-02 2020-05-07 コニカミノルタ株式会社 超音波診断装置、超音波診断装置の制御方法、及び、超音波診断装置の制御プログラム
JP7192404B2 (ja) 2018-11-02 2022-12-20 コニカミノルタ株式会社 超音波診断装置、超音波診断装置の制御方法、及び、超音波診断装置の制御プログラム
US12350100B2 (en) 2020-11-16 2025-07-08 Ablic Inc. Ultrasonic device and control method of ultrasonic device

Also Published As

Publication number Publication date
JP4290699B2 (ja) 2009-07-08
CN100539950C (zh) 2009-09-16
JPWO2005065547A1 (ja) 2007-12-20
US8246544B2 (en) 2012-08-21
US20090198138A1 (en) 2009-08-06
EP1709908A4 (en) 2009-07-29
EP1709908A1 (en) 2006-10-11
CN1909836A (zh) 2007-02-07

Similar Documents

Publication Publication Date Title
EP2107388A1 (en) Ultrasound diagnosis apparatus
US20120157850A1 (en) Ultrasound diagnosis apparatus
WO2005065547A1 (ja) 超音波診断装置
JP2008080104A (ja) 超音波映像を形成する超音波システム及び方法
US20080087089A1 (en) Apparatus and method for forming an ultrasound image
CN102695456B (zh) 波束形成方法、超声波诊断装置、集成电路
JP2009125593A (ja) 超音波システム
CN114027873A (zh) 一种超声成像方法、装置及计算机可读存储介质
JP2008018111A (ja) ビームフォーミング装置及び方法
CN110267599B (zh) 一种超声成像方法及装置、计算机可读存储介质
JP7581682B2 (ja) 超音波診断装置、超音波信号処理方法、及びプログラム
CN115177290A (zh) 复合超声成像方法及超声成像设备
JP2011050742A (ja) ステアリング角度を考慮して超音波合成映像を提供する超音波システムおよび方法
WO2012067391A2 (ko) 컬러 플로우 영상 촬상 방법과 그를 위한 초음파 장치
CN115137392B (zh) 一种抑制杂波的超声波束合成方法、系统及存储介质
JP3987388B2 (ja) 超音波診断装置
JP6041957B1 (ja) 超音波診断装置
US10925578B2 (en) Acoustic wave image generating apparatus and control method thereof
US12343212B2 (en) Ultrasonic diagnostic apparatus and image processing apparatus
JP4864532B2 (ja) 超音波診断装置、画像データ表示装置及び3次元画像データ生成方法
JP4664209B2 (ja) 超音波診断装置およびその撮像を実行する超音波イメージングプログラム
JP2020069304A (ja) 超音波診断装置、超音波診断装置の制御方法、及び、超音波診断装置の制御プログラム
KR20010077539A (ko) 초음파 영상 시스템용 인터레이싱 다중 빔 집속 장치 및집속 방법
JP5634817B2 (ja) 超音波診断装置
JP2010081965A (ja) 超音波診断装置

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): BW GH GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2005516882

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 10585388

Country of ref document: US

Ref document number: 200580002088.X

Country of ref document: CN

NENP Non-entry into the national phase

Ref country code: DE

WWW Wipo information: withdrawn in national office

Ref document number: DE

WWE Wipo information: entry into national phase

Ref document number: 2005703342

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 2005703342

Country of ref document: EP