[go: up one dir, main page]

JP2018110784A - 超音波診断装置、および、その制御方法 - Google Patents

超音波診断装置、および、その制御方法 Download PDF

Info

Publication number
JP2018110784A
JP2018110784A JP2017004351A JP2017004351A JP2018110784A JP 2018110784 A JP2018110784 A JP 2018110784A JP 2017004351 A JP2017004351 A JP 2017004351A JP 2017004351 A JP2017004351 A JP 2017004351A JP 2018110784 A JP2018110784 A JP 2018110784A
Authority
JP
Japan
Prior art keywords
ultrasonic
acoustic line
unit
line signal
transmission
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2017004351A
Other languages
English (en)
Other versions
JP6724797B2 (ja
Inventor
津島 峰生
Mineo Tsushima
峰生 津島
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.)
Konica Minolta Inc
Original Assignee
Konica Minolta Inc
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 Konica Minolta Inc filed Critical Konica Minolta Inc
Priority to JP2017004351A priority Critical patent/JP6724797B2/ja
Priority to EP17205915.6A priority patent/EP3349034A3/en
Priority to US15/861,936 priority patent/US20180199918A1/en
Priority to CN201810018290.0A priority patent/CN108294778B/zh
Publication of JP2018110784A publication Critical patent/JP2018110784A/ja
Application granted granted Critical
Publication of JP6724797B2 publication Critical patent/JP6724797B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4427Device being portable or laptop-like
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4483Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B06GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS IN GENERAL
    • B06BMETHODS OR APPARATUS FOR GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS OF INFRASONIC, SONIC, OR ULTRASONIC FREQUENCY, e.g. FOR PERFORMING MECHANICAL WORK IN GENERAL
    • B06B1/00Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency
    • B06B1/02Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy
    • B06B1/06Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy operating with piezoelectric effect or with electrostriction
    • B06B1/0607Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy operating with piezoelectric effect or with electrostriction using multiple elements
    • B06B1/0622Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy operating with piezoelectric effect or with electrostriction using multiple elements on one surface
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4444Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4444Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • A61B8/4472Wireless probes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/22Details, e.g. general constructional or apparatus details
    • G01N29/24Probes
    • G01N29/2437Piezoelectric probes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/22Details, e.g. general constructional or apparatus details
    • G01N29/24Probes
    • G01N29/2481Wireless probes, e.g. with transponders or radio links
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/34Generating the ultrasonic, sonic or infrasonic waves, e.g. electronic circuits specially adapted therefor
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/36Detecting the response signal, e.g. electronic circuits specially adapted therefor
    • 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/8909Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • 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/8997Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using synthetic aperture techniques
    • 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/52046Techniques for image enhancement involving transmitter or receiver
    • G01S7/52049Techniques for image enhancement involving transmitter or receiver using correction of medium-induced phase aberration
    • 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/52082Constructional features involving a modular construction, e.g. a computer with short range imaging equipment
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06GANALOGUE COMPUTERS
    • G06G7/00Devices in which the computing operation is performed by varying electric or magnetic quantities
    • G06G7/12Arrangements for performing computing operations, e.g. operational amplifiers
    • G06G7/14Arrangements for performing computing operations, e.g. operational amplifiers for addition or subtraction 
    • HELECTRICITY
    • H10SEMICONDUCTOR DEVICES; ELECTRIC SOLID-STATE DEVICES NOT OTHERWISE PROVIDED FOR
    • H10NELECTRIC SOLID-STATE DEVICES NOT OTHERWISE PROVIDED FOR
    • H10N39/00Integrated devices, or assemblies of multiple devices, comprising at least one piezoelectric, electrostrictive or magnetostrictive element covered by groups H10N30/00 – H10N35/00
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4405Device being mounted on a trolley
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B06GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS IN GENERAL
    • B06BMETHODS OR APPARATUS FOR GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS OF INFRASONIC, SONIC, OR ULTRASONIC FREQUENCY, e.g. FOR PERFORMING MECHANICAL WORK IN GENERAL
    • B06B1/00Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency
    • B06B1/02Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy
    • B06B1/0207Driving circuits
    • B06B1/0215Driving circuits for generating pulses, e.g. bursts of oscillations, envelopes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/01Indexing codes associated with the measuring variable
    • G01N2291/012Phase angle
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/10Number of transducers
    • G01N2291/106Number of transducers one or more transducer arrays

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Acoustics & Sound (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Chemical & Material Sciences (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Mathematical Physics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Gynecology & Obstetrics (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

【課題】集束型の送信ビームフォーミングを用いた合成開口法を行う超音波診断装置であって、必要な演算能力に応じて構成及び処理内容を変更することのできる超音波診断装置を提供する。【解決手段】少なくとも対象領域の一部を含む第1対象領域内の複数の観測点について、受信信号列を整相加算し第1音響線信号を生成する第1整相加算部と、第1音響線信号に基づいてパラメータを算出するパラメータ算出部と、対象領域の全部または一部である第2対象領域内の複数の観測点について、パラメータに基づき、受信信号列を整相加算しサブフレーム音響線信号を生成する第2整相加算部と、サブフレーム音響線信号を合成してフレーム音響線信号を生成する合成部と、第1音響線信号とフレーム音響線信号のいずれに基づいて超音波画像を生成するかを決定する制御部と、超音波画像を生成する超音波画像生成部とを備える超音波診断装置。【選択図】図24

Description

本開示は、超音波診断装置に関し、特に、超音波診断装置における送受信ビームフォーミング処理方法、制御方法、および装置システム構成に関する。
超音波診断装置は、超音波プローブ(以後、「プローブ」とする)により被検体内部に超音波を送信し、被検体組織の音響インピーダンスの差異により生じる超音波反射波(エコー)を受信する。さらに、この受信から得た電気信号に基づいて、被検体の内部組織の構造を示す超音波断層画像を生成して表示するものである。超音波診断装置は、被検体への侵襲が少なく、リアルタイムに体内組織の状態を断層画像などで観察できるため、生体の形態診断に広く用いられている。
従来、受信した反射超音波に基づく信号の受信ビームフォーミング方法として、一般的に整相加算法と呼ばれる方法が使用されている(例えば、非特許文献1)。この方法では、一般に、被検体への超音波送信が行われる際、被検体のある深さで超音波ビームがフォーカスを結ぶよう送信ビームフォーミングがなされる。また、送信超音波ビームの中心軸上またはその近傍に観測点を設定する。そのため、超音波主照射領域の面積に比べて観測点数が少なく超音波の利用効率が悪い。また、観測点が送信フォーカス点近傍から離れた位置にある場合には、得られる音響線信号の空間分解能及び信号S/N比が低くなる課題も有している。なお、超音波主照射領域とは、超音波ビームが伝播する領域を指す。
これに対して、合成開口法(Synthetic Aperture Method)により、送信フォーカス点近傍以外の領域においても空間分解能の高い、高画質な画像を得る受信ビームフォーミング方法が考案されている(例えば、非特許文献2)。この方法によれば、超音波送信波の伝播経路と、その伝播経路による反射波の振動子への到達時間の両方を加味した遅延制御を行うことで、送信フォーカス点近傍以外に位置する超音波主照射領域からの反射超音波も反映した受信ビームフォーミングを行うことができる。その結果、送信超音波ビームの中心軸のみならず超音波主照射領域全体に対して音響線信号(受信ビームフォーミングにより生成される、観測点からの反射超音波に基づく信号)を生成することができる。また、合成開口法では、複数の送信イベントから得た同一観測点に対する複数の受信信号をもとに仮想的に送信フォーカスを合わせることで、非特許文献1記載の受信ビームフォーミング方法と比較して、空間分解能及びS/N比の高い超音波画像を得ることが可能となる。
特表2008−536578号公報
伊東正安、望月剛共著「超音波診断装置」コロナ社出版、2002年8月26日(P42−P45) "Virtual ultrasound sources in high resolution ultrasound imaging", S.I.Nikolov and J.A.Jensen, in Proc, SPIE − Progress in biomedical optics and imaging, vol. 3, 2002, P. 395−405 "Synthetic Aperture Sequential Beamforming", Jacob Kortbek, et.al., IEEE Ultrasonics Symposium, 2−5 Nov. 2008 pp 966−969
一方で、合成開口法においては、超音波利用効率と解像度向上の観点から、1回の超音波送信イベントで音響線信号を生成する領域(以下、「対象領域」と呼ぶ)の面積が大きいことが好ましく、超音波主照射領域全域を対象領域とすることがより好ましい。しかしながら、対象領域の面積が大きくなると、その内部に存在する観測点(受信ビームフォーミングの演算対象となる場所)の数が対象領域の面積に比例して増加するため、送信と受信の遅延を考慮した整相加算の演算量が増加する。そのため、超音波主照射領域の面積が大きくなると、整相加算の演算処理を高速に行うため演算処理能力の高いハードウエアが必要となる。しかしながら、可搬型の超音波診断装置では、例えば、装置の大きさ、排熱、バッテリーによる駆動可能時間などにより、演算能力の向上について制約が生じることが少なくない。一方で、用途に応じて演算能力の異なる超音波診断装置を用いる場合、超音波診断装置の台数が増加し非効率である。
本発明は、上記課題に鑑みてなされたものであり、集束型の送信ビームフォーミングを用いた合成開口法を行う超音波診断装置であって、必要な演算能力に応じて構成及び処理内容を変更することのできる超音波診断装置を提供することを目的とする。
本発明の一態様に係る超音波診断装置は、複数の振動子を備えた超音波プローブを用いて被検体に集束型の超音波ビームを送信する送信イベントを複数回繰り返すとともに、各送信イベントに同期して被検体から反射超音波を受波して受信信号列を生成し、前記受信信号列に基づいて超音波画像を生成する超音波診断装置であって、超音波ビームの集束する位置を規定するフォーカス点を送信イベントごとに変更しながら、送信イベントごとに、前記超音波プローブに列設された複数の振動子から送信振動子列を選択して当該送信振動子列から被検体内の対象領域に超音波ビームを送信させる送信部と、各送信イベントに同期して、前記超音波プローブが前記対象領域から受波した反射超音波に基づいて、振動子各々に対する受信信号列を生成する受信部と、前記送信イベントごとに、少なくとも前記対象領域の一部を含む第1対象領域内の複数の観測点について、前記受信信号列を整相加算し第1音響線信号を生成する第1整相加算部と、前記第1音響線信号に基づいて、サブフレーム音響線信号を生成するためのパラメータを算出するパラメータ算出部と、前記対象領域の全部または一部である第2対象領域内の複数の観測点について、前記パラメータに基づき、前記受信信号列を整相加算しサブフレーム音響線信号を生成する第2整相加算部と、前記サブフレーム音響線信号を合成してフレーム音響線信号を生成する合成部と、前記第1音響線信号と前記フレーム音響線信号のいずれに基づいて超音波画像を生成するかを決定する制御部と、前記制御部の決定に基づき、前記第1音響線信号と前記フレーム音響線信号のいずれか一方から超音波画像を生成する超音波画像生成部とを含む超音波信号処理回路を備えることを特徴とする。
本発明の一態様に係る超音波診断装置によれば、必要な演算能力に応じて、第2整相加算部を用いる動作と、第2整相加算部を用いない動作とを切り替えることができる。
実施の形態1に係る超音波診断システム1000の構成を示す機能ブロック図である。 実施の形態1に係る送信ビームフォーマ部103による超音波送信波の伝播経路を示す模式図である。 実施の形態1に係る受信ビームフォーマ部104の構成を示す機能ブロック図である。 実施の形態1に係る第1整相加算部1041の構成を示す機能ブロック図である。 実施の形態1に係る対象領域Bxと観測線BL、代表点Qkを示す模式図である。 (a)、(b)はいずれも、実施の形態1に係る、送信開口Txから代表点Qkを経由して受信振動子Rmに到達する超音波の伝播経路を示す模式図である。 受信振動子Rmごとの遅延時間の関係を示す模式図である。 実施の形態1に係る第2整相加算部1042の構成を示す機能ブロック図である。 実施の形態1に係る、観測点Pijの位置と重み算出部1247が算出する重み数列との関係を示す模式図である。 (a)、(b)はいずれも、実施の形態1に係る、送信開口Txから観測点Pijを経由して受信振動子Rijに到達する超音波の伝播経路を示す模式図である。 (a)、(b)はいずれも、実施の形態に係る、遅延量と超音波速度との関係を示す模式図である。 (a)は実施の形態1に係る第1合成部1044の、(b)は実施の形態1に係る第2合成部1045の、それぞれの構成を示す機能ブロック図である。 実施の形態1に係る、音響線信号展開部1341におけるサブフレーム音響線信号生成処理を示す模式図である。 実施の形態1に係る、重みづけ合成部1342における重みづけ合成処理と重み数列を示す模式図である。 実施の形態1に係る、音響線信号における最大重畳数と、増幅処理の概要を示す模式図である。 実施の形態1に係る、第2整相加算部1042を使用する場合における、受信ビームフォーマ部104のビームフォーミング処理を示すフローチャートである。 実施の形態1に係る、線領域音響線信号の生成動作を示すフローチャートである。 実施の形態1に係る、代表点Qkについての音響線信号生成動作を説明するための模式図である。 実施の形態1に係る、サブフレーム音響線信号の生成動作を示すフローチャートである。 実施の形態1に係る、観測点Pijについての音響線信号生成動作を説明するための模式図である。 実施の形態1に係る、第2整相加算部1042を使用しない場合における、受信ビームフォーマ部104のビームフォーミング処理を示すフローチャートである。 変形例に係る、ゲイン補正値を変更した場合の、観測点の深さと音響線信号のS/N比の関係を示す模式図である。 変形例に係る、対象領域Cxの幅を変更する制御を示す模式図である。 実施の形態2に係る超音波診断システム2000、可搬本体部200、カート203のそれぞれを示す模式図である。 実施の形態2に係る受信ビームフォーマ部304の構成を示す機能ブロック図である。
≪実施の形態1≫
<全体構成>
以下、実施の形態1に係る超音波診断装置100について、図面を参照しながら説明する。
図1は、実施の形態1に係る超音波診断システム1000の機能ブロック図である。図1に示すように、超音波診断システム1000は、被検体に向けて超音波を送信しその反射波を受信する複数の振動子101aを有するプローブ101、プローブ101に超音波の送受信を行わせプローブ101からの出力信号に基づき超音波画像を生成する超音波診断装置100、超音波画像を画面上に表示する表示部106を有する。プローブ101、表示部106は、それぞれ、超音波診断装置100に各々接続可能に構成されている。図1は超音波診断装置100に、プローブ101、表示部106が接続された状態を示している。なお、プローブ101と、表示部106とは、超音波診断装置100の内部にあってもよい。
<超音波診断装置100の構成>
超音波診断装置100は、プローブ101の複数ある振動子101aのうち送信又は受信の際に用いる振動子のそれぞれについて入出力を確保するマルチプレクサ部102、超音波の送信を行うためにプローブ101の各振動子101aに対する高電圧印加のタイミングを制御する送信ビームフォーマ部103と、プローブ101で受信した超音波の反射波に基づき、複数の振動子101aで得られた電気信号を増幅し、A/D変換し、受信ビームフォーミングして音響線信号を生成する受信ビームフォーマ部104を有する。また、受信ビームフォーマ部104からの出力信号に基づいて超音波画像(Bモード画像)を生成する超音波画像生成部105、受信ビームフォーマ部104が出力する音響線信号および超音波画像生成部105が出力する超音波画像を保存するデータ格納部107と、各構成要素を制御する制御部108を備える。
このうち、マルチプレクサ部102、送信ビームフォーマ部103、受信ビームフォーマ部104、超音波画像生成部105は、超音波信号処理装置150を構成する。
超音波診断装置100を構成する、マルチプレクサ部102、送信ビームフォーマ部103、超音波画像生成部105、制御部108は、それぞれ、例えば、FPGA(Field Programmable Gate Array)、ASIC(Application Specific Ingegrated Circuit)などのハードウエア回路により実現される。あるいは、CPU(Central Processing Unit)に代表されるプロセッサなどのプログラマブルデバイスとソフトウェアにより実現される構成であってもよい。これらの構成要素は一個の回路部品とすることができるし、複数の回路部品の集合体にすることもできる。また、複数の構成要素を組合せて一個の回路部品とすることができるし、複数の回路部品の集合体にすることもできる。なお、受信ビームフォーマ部104については後述する。
データ格納部107は、コンピュータ読み取り可能な記録媒体であり、例えば、フレキシブルディスク、ハードディスク、光磁気ディスク、光学ディスク、半導体メモリ等を用いることができる。また、データ格納部107は、超音波診断装置100に外部から接続された記憶装置であってもよい。
なお、本実施の形態1に係る超音波診断装置100は、図1で示した構成の超音波診断装置に限定されない。例えば、マルチプレクサ部102がなく、送信ビームフォーマ部103と受信ビームフォーマ部104とが直接、プローブ101の各振動子101aに接続されていてもよい。また、プローブ101に送信ビームフォーマ部103や受信ビームフォーマ部104、またその一部などが内蔵される構成であってもよい。これは、本実施の形態に係る超音波診断装置100に限られず、後に説明する他の実施の形態や変形例に係る超音波診断装置でも同様である。
<超音波診断装置100の主要部の構成>
実施の形態1に係る超音波診断装置100は、プローブ101の各振動子101aから超音波ビームの送信を行わせる送信ビームフォーマ部103と、プローブ101において超音波反射波の受信から得た電気信号をA/D変換し演算して超音波画像を生成するための音響線信号を生成する受信ビームフォーマ部104に特徴を有する。そのため、本明細書では、主に、送信ビームフォーマ部103及び受信ビームフォーマ部104について、その構成及び機能を説明する。なお、送信ビームフォーマ部103及び受信ビームフォーマ部104以外の構成については、公知の超音波診断装置に使われるものと同じ構成を適用可能であり、公知の超音波診断装置のビームフォーマ部に本実施の形態に係るビームフォーマ部を置き換えて使用することが可能である。
以下、送信ビームフォーマ部103と、受信ビームフォーマ部104の構成について説明する。
1.送信ビームフォーマ部103
送信ビームフォーマ部103は、マルチプレクサ部102を介してプローブ101と接続され、プローブ101から超音波の送信を行うためにプローブ101に存する複数の振動子101aの全てもしくは一部に当たる送信振動子列からなる送信開口Txに含まれる複数の振動子の各々に対する高電圧印加のタイミングを制御する。送信ビームフォーマ部103は送信部1031から構成される。
送信部1031は、制御部108からの送信制御信号に基づき、プローブ101に存する複数の振動子101a中、送信開口Txに含まれる各振動子に超音波ビームを送信させるためのパルス状の送信信号を供給する送信処理を行う。具体的には、送信部1031は、例えば、クロック発生回路、パルス発生回路、遅延回路を備えている。クロック発生回路は、超音波ビームの送信タイミングを決定するクロック信号を発生させる回路である。パルス発生回路は、各振動子を駆動するパルス信号を発生させるための回路である。遅延回路は、超音波ビームの送信タイミングを振動子毎に遅延時間を設定し、遅延時間だけ超音波ビームの送信を遅延させて超音波ビームのフォーカシングを行うための回路である。
送信部1031は、超音波送信ごとに送信開口Txを列方向に移動ピッチMpだけ移動させながら超音波送信を繰り返し、プローブ101に存する全ての振動子101aから超音波送信を行う。本実施の形態では、移動ピッチMpを振動子1つ分とし、超音波送信ごとに、送信開口Txが振動子1つ分ずつ移動していく。なお、移動ピッチMpは振動子1つ分に限られず、例えば、振動子0.5個分としてもよい。送信開口Txに含まれる振動子の位置を示す情報は制御部108を介してデータ格納部107に出力される。例えば、プローブ101に存する振動子101a全数を192としたとき、送信開口Txを構成する振動子列の数として、例えば20〜128を選択してもよく、超音波送信毎に振動子1つ分だけ移動させる構成としてもよい。以後、送信部1031により同一の送信開口Txから行われる超音波送信を「送信イベント」と称呼する。
図2は、送信ビームフォーマ部103による超音波送信波の伝播経路を示す模式図である。ある送信イベントにおいて、超音波送信に寄与するアレイ状に配列された振動子101aの列(送信振動子列)を送信開口Txとして図示している。また、送信開口Txの列長を送信開口長と呼ぶ。
送信ビームフォーマ部103において、送信開口Txの中心に位置する振動子ほど送信タイミングを遅らせるように各振動子の送信タイミングを制御する。これにより、送信開口Tx内の振動子列から送信された超音波送信波は、被検体のある深度(Focal depth)において、波面がある一点、すなわち送信フォーカス点F(Focal point)で、フォーカスがあう(集束する)状態となる。送信フォーカス点Fの深さ(Focal depth)(以下、「フォーカス深さ」とする)は、任意に設定することができる。ここで、フォーカス深さは、超音波送信波が振動子の並ぶ方向(図2におけるx方向)に最も集束する深さ、すなわち、超音波ビームのx方向における幅が最も狭まるy方向の深さである。送信フォーカス点Fは、フォーカス深さにおける超音波ビームのx方向における中心位置である。但し、1フレームに係る複数の送信イベント中ではフォーカス深さは一定とする。すなわち、1フレームに係る複数の送信イベントにおいて、送信開口Txと送信フォーカス点Fとの相対的な関係は変化しない。送信フォーカス点Fで合焦した波面は、再び拡散し、送信開口Txを底とし送信フォーカス点Fを節とする交差する2つの直線で区切られた砂時計型の空間内を超音波送信波が伝播する。すなわち、送信開口Txで放射された超音波は、次第にその空間上での幅(図中の横軸方向)を小さくし、送信フォーカス点Fでその幅を最小化し、それよりも深部(図中では上部)に進行するにしたがって、再び、その幅を大きくしながら拡散し、伝播することとなる。この砂時計型の領域が超音波主照射領域Axである。なお、上述したように、超音波主照射領域Axは、1点の送信フォーカス点Fの近傍に集束するように超音波送信波を送信してもよい。
2.受信ビームフォーマ部104の構成
受信ビームフォーマ部104は、プローブ101で受信した超音波の反射波に基づき、複数の振動子101aで得られた電気信号から音響線信号を生成する。なお、「音響線信号」とは、ある観測点に対する、整相加算処理がされた後の信号である。整相加算処理については後述する。図3は、受信ビームフォーマ部104の構成を示す機能ブロック図である。図3に示すように、受信ビームフォーマ部104は、受信部1040、第1整相加算部1041、第2整相加算部1042、パラメータ算出部1043、第1合成部1044、第2合成部1045、出力部1046を備える。
このうち、第2整相加算部1042、第2合成部1045は、高機能演算回路1047を構成する。
高機能演算回路1047は、例えば、プロセッサなどのプログラマブルデバイスとソフトウェアにより実現される。プロセッサとしては、CPUやGPU(Graphics Processing Unit)を用いることができ、GPUを用いる構成はGPGPU(General−Purpose computing on Graphics Processing Unit)と呼ばれる。また、受信ビームフォーマ部104のうち高機能演算回路1047を除く各要素、すなわち、受信部1040、第1整相加算部1041、パラメータ算出部1043、第1合成部1044、出力部1046は、例えば、FPGA(Field Programmable Gate Array)、ASIC(Application Specific Ingegrated Circuit)などのハードウエア回路、あるいは、プロセッサなどのプログラマブルデバイスとソフトウェアにより実現される。なお、第2整相加算部1042、第2合成部1045は、それぞれ、第1整相加算部1041、第1合成部1044より高い演算能力を備えている。
以下、受信ビームフォーマ部104を構成する各部の構成について説明する。
(1)受信部1040
受信部1040は、マルチプレクサ部102を介してプローブ101と接続され、送信イベントに同期してプローブ101での超音波反射波の受信から得た電気信号を増幅した後AD変換した受信信号(RF信号)を生成する回路である。送信イベントの順に時系列に受信信号を生成しデータ格納部107に出力し、データ格納部107に受信信号を一時的に保存する。
ここで、受信信号(RF信号)とは、各振動子にて受信された反射超音波から変換された電気信号をA/D変換したデジタル信号であり、各振動子にて受信された超音波の送信方向(被検体の深さ方向)に連なった信号の列を形成している。
送信イベントでは、上述のとおり、送信部1031は、プローブ101に存する複数の振動子101a中、送信開口Txに含まれる複数の振動子の各々に超音波ビームを送信させる。これに対し、受信部1040は、送信イベントに同期してプローブ101に存する複数の振動子101aの一部又は全部にあたる振動子(以下、「受波振動子」という)の各々が得た反射超音波に基づいて、各振動子に対する受信信号の列を生成する。受波振動子の数は、送信開口Txに含まれる振動子の数よりも多いことが好ましい。また、受波振動子の数はプローブ101に存する振動子101aの全数としてもよい。
送信部1031は、送信イベントに同期して送信開口Txを列方向に移動ピッチMpだけ移動させながら超音波送信を繰り返し、プローブ101に存する複数の振動子101a全体から超音波送信を行う。受信部1040は、送信イベントに同期して各受波振動子に対する受信信号の列を生成し、生成された受信信号はデータ格納部107に保存される。
(2)第1整相加算部1041
第1整相加算部1041は、送信イベントに同期して、被検体においてサブフレーム音響線信号の生成を行う対象領域Bxを設定する。次に、対象領域Bx内に、送信フォーカス点Fを通過する観測線BLを設定する。本実施の形態では、観測線BLは送信フォーカス点Fと送信開口Txの中心線とを通過する直線であり、振動子列に直交している。なお、観測線BLは送信フォーカス点Fと送信開口Tx内の任意の点を通過していればよく、上述の場合に限られない。次に、観測線BL上に存在する複数の代表点Qkの各々について、代表点Qkから各受信振動子Rmが受信した受信信号列を整相加算する。そして、各代表点Qkにおける音響線信号の列を算出することにより線領域音響線信号を生成する回路である。図4は、第1整相加算部1041の構成を示す機能ブロック図である。図4に示すように、第1整相加算部1041は、対象領域設定部1141、送信時間算出部1142、受信時間算出部1143、遅延量算出部1144、遅延処理部1145、重み算出部1146、及び加算部1147を備える。
以下、第1整相加算部1041を構成する各部の構成について説明する。
i)対象領域設定部1141
対象領域設定部1141は、被検体内においてサブフレーム音響線信号の生成を行う対象領域Bxを設定する。「対象領域」とは、送信イベントに同期して被検体内においてサブフレーム音響線信号の生成が行われるべき信号上の領域である。すなわち、対象領域Bxは、音響線信号の生成が行われる観測対象点の集合として、1回の送信イベントに同期して計算の便宜上設定される。ここで、「サブフレーム音響線信号」とは、1回の送信イベントから生成される対象領域Bx内に存在する全ての観測点に対する音響線信号の集合である。「サブフレーム」とは、1回の送信イベントで得られ、対象領域Bx内に存在する全ての観測点に対応するまとまった信号を形成する単位を指し、取得時間の異なる複数のサブフレームを合成したものがフレームとなる。
対象領域設定部1141は、送信イベントに同期して、送信ビームフォーマ部103から取得する送信開口Txの位置を示す情報に基づき対象領域Bxを設定する。
図5は、対象領域Bxを示す模式図である。図5に示すように、対象領域Bxは、超音波主照射領域Ax内に存在する任意の領域であり、本実施の形態では、超音波主照射領域Axの全域である。
さらに、対象領域設定部1141は、対象領域Bxの内部に、線領域音響線信号の生成を行う対象線BLを設定する。本実施の形態では、対象線BLは、フォーカス点Fまたはその近傍を通過する直線である。なお、上述したように、対象線BLはフォーカス点Fまたはその近傍と、送信開口Tx上の任意の点とを通過する直線上の領域であればよい。そして、この対象線BL上に設定される代表点Qkに対し、線領域音響線信号が生成される。なお、受信開口は、送信開口Txをそのまま用いる。
設定された対象領域Bxは制御部103に出力され、対象線BL、および、送信ビームフォーマ部103から取得した送信開口Txは、送信時間算出部1142、受信時間算出部1143、遅延処理部1145、重み算出部1146に出力される。
ii)送信時間算出部1142
送信時間算出部1142は、送信された超音波が被検体中の観測点Pに到達する送信時間を算出する回路である。送信イベントに対応して、対象領域設定部1141から取得した、送信開口Txに含まれる振動子の位置を示す情報と、対象線BLの位置を示す情報とに基づき対象線BL上に存在する任意の代表点Qkについて、送信された超音波が被検体中の代表点Qkに到達する送信時間を算出する。
図6は、送信開口Txから放射され対象線BL上の任意の位置にある代表点Qkにおいて反射され送信開口Tx内に位置する受信振動子Rmに到達する超音波の伝播経路を説明するための模式図である。なお、図6(a)は代表点Qkの深さが送信フォーカス深さ以上の場合、図6(b)は代表点Qkが送信フォーカス深さより浅い場合を示している。
送信開口Txから放射された送信波は、経路401を通って送信フォーカス点Fにて波面が集束し、再び、拡散する。送信波が集束または拡散する途中で代表点Qkに到達し、代表点Qkで音響インピーダンスに変化があれば反射波を生成し、その反射波がプローブ101における送信開口Tx内の受信振動子Rmに戻っていく。送信フォーカス点Fは送信ビームフォーマ部103の設計値として規定されているので、送信フォーカス点Fと任意の代表点Qkとの間の経路402の長さは幾何学的に算出することができる。
送信時間の算出方法を、以下、さらに詳細に説明する。
まず、代表点Qkの深さが送信フォーカス深さ以上の場合について、図6(a)を用いて説明する。この場合、送信開口Txから放射された送信波が、経路401を通って送信フォーカス点Fに到達し、送信フォーカス点Fから経路402を通って代表点Qkに到達したものとして算出する。したがって、送信波が経路401を通過する時間と、経路402を通過する時間を合算した値が、送信時間となる。具体的には、経路401の長さと経路402の長さとを加算した全経路長を、被検体内における超音波の伝搬速度で除算することで求められる。
一方、代表点Qkが送信フォーカス深さより浅い場合について、図6(b)を用いて説明する。この場合は、送信開口Txから放射された送信波が、経路401を通って送信フォーカス点Fに到達する時刻と、経路404を通って代表点Qkに到達した後、代表点Qkから経路405を通って送信フォーカス点Fに到達する時刻とが同一であるものとして算出する。つまり、送信波が経路401を通過する時間から、経路405を通過する時間を差し引いた値が、送信時間となる。具体的には、経路401の長さから経路405の長さを減算した経路長差を、被検体内における超音波の伝搬速度で除算することで求められる。
なお、代表点Qkが送信フォーカス深さである場合の送信時間は、代表点Qkが送信フォーカス深さより浅い場合と同じ算出方法を用いてもよい。経路402の長さと経路405の長さがいずれも0となるため、いずれの算出方法でも、経路401を通過する時間と一致するためである。
送信時間算出部1142は、1回の送信イベントに対し、対象線BL上の全ての代表点Qkについて、送信された超音波が被検体中の観測点Qkに到達する送信時間を算出して遅延量算出部1144に出力する。また、送信時間算出部1142は、1回の送信イベントに対し、対象線BL上の全ての代表点Qkについて、経路402または経路405の長さを、受信時間算出部1143に出力する。
iii)受信時間算出部1143
受信時間算出部1143は、代表点Qからの反射波が、送信開口Txに含まれる受信振動子Rmの各々に到達する受信時間を算出する回路である。送信イベントに対応して、対象領域設定部1141から取得した受信振動子Rkの位置を示す情報と、対象線BLの位置を示す情報とに基づき対象線BL上に存在する任意の代表点Qkについて、送信された超音波が被検体中の代表点Qkで反射され送信開口Txの各受信振動子Rmに到達する受信時間を算出する。
上述のとおり、代表点Qkに到達した送信波は、代表点Qkで音響インピーダンスに変化があれば反射波を生成し、その反射波がプローブ101における送信開口Tx内の各受信振動子Rmに戻っていく。このとき、受信時間算出部1143は、送信超音波ビームと同様、送信フォーカス点Fを基準として代表点Qkから受信振動子Rmまでの経路を算出する。
受信時間の算出方法について、まず、図6を用いて概念を説明する。但し、算出に当たっては、後述するように、単純化を行うことができる。
まず、代表点Qkの深さが送信フォーカス深さ以上の場合について、図6(a)を用いて説明する。この場合は、代表点Qkが反射した反射波が、経路402を通って送信フォーカス点Fに到達し、送信フォーカス点Fから経路403を通って受信振動子Rmに到達したものとして算出する。したがって、経路402を通過する時間と、経路403を通過する時間を合算した値が、受信時間となる。
一方、代表点Qkが送信フォーカス深さより浅い場合について、図6(b)を用いて説明する。この場合は、送信フォーカス点Fが反射した反射波が、経路405を通って代表点Qkに到達した後経路406を通って受信振動子Rmに到達する時刻と、経路403を通って直接受信振動子Rmに到達する時刻とが同一であるものとして算出する。つまり、代表点Qkが反射した反射波が受信振動子Rmに到達する時間は、送信フォーカス点Fが反射した反射波が経路403を通って受信振動子Rmに到達する時間より、経路405の通過に必要な時間だけ短い、とする。したがって、経路403を通過する時間から経路405を通過する時間を差し引いた値が、受信時間となる。
ここで、代表点Qkごとの経路402または経路405の長さは、送信時間算出部1142が送信時間の一部として算出する代表点Qkごとの経路402または経路405の長さと同一となる。そのため、本実施の形態では、送信時間算出部1142が算出した代表点Qkごとの経路402または経路405の長さを、受信時間の算出に用いる。また、経路403の長さは送信フォーカス点Fと受信振動子Rmとの位置関係のみに依存する。言い換えると、同一の代表点Qに対する2つの受信振動子Rm1、Rm2における受信時間の差は、代表点の位置には全く依存しない。すなわち、同一の代表点Qに対する2つの受信振動子Rm1、Rm2における受信時間の差は、代表点Qk1に対しても、代表点Qk2に対しても、代表点Qk3に対しても、一定となる。
以下、図7を用いてより詳細に説明する。経路403の長さは、受信振動子Rmと送信フォーカス点Fとの位置関係によって定まる。受信振動子Rmの受信時間と送信開口Txの中心に位置する受信振動子Rcの受信時間との差は、送信フォーカス点Fを中心として受信振動子Rcに接する円弧410と、受信振動子Rmとの距離412を超音波が伝播するのに要する時間である。
したがって、受信時間算出部1143は、受信振動子Rcの経路403の長さに対応する経路401の長さと、送信時間算出部1142が算出した代表点Qkごとの経路402または経路405の長さとを用いて、受信振動子Rcに対する代表点Qkごとの受信時間を算出する。さらに、受信振動子Rmごとの距離412を超音波の伝播速度で除することで、受信振動子Rmごとの受信振動子Rcに対する受信時間差を算出する。そして、受信振動子Rcに対する代表点Qkごとの受信時間と、受信振動子Rmごとの受信振動子Rcに対する受信時間差を、遅延量算出部1144に出力する。
iv)遅延量算出部1144
遅延量算出部1144は、送信時間と受信時間とから送信開口Tx内の各受信振動子Rmへの総伝播時間を算出し、当該総伝播時間に基づいて、各受信振動子Rmに対する受信信号の列に適用する遅延量を算出する回路である。遅延量算出部1144は、送信された超音波が代表点Qkに到達する送信時間と、代表点Qkで反射して受信振動子Rcに到達する受信時間と、受信振動子Rcと各受信振動子Rmとの受信時間差とを取得する。そして、送信された超音波が各受信振動子Rmへ到達するまでの総伝播時間を算出し、各受信振動子Rmに対する総伝播時間の差異により、各受信振動子Rmに対する遅延量を算出する。代表点Qkごとの受信振動子Rcに対する総伝播時間は、代表点Qkに対する送信時間と、受信振動子Rcに対する受信時間との合計として得ることができる。また、各受信振動子Rmに対する総伝播時間は、同一の代表点Qkに対する受信振動子Rcの総伝播時間に、受信振動子Rcと各受信振動子Rmとの受信時間差を加算することで得ることができる。遅延量算出部1144は、対象線BL上に存在する全ての代表点Qkについて、受信振動子Rcに対する受信信号の列に適用する遅延量を算出し、受信振動子Rcと各受信振動子Rmとの受信時間差と併せて遅延処理部1145に出力する。
v)遅延処理部1145
遅延処理部1145は、送信開口Tx内の受信振動子Rmに対する受信信号の列から、各受信振動子Rmに対する遅延量に相当する受信信号を、代表点Qkからの反射超音波に基づく各受信振動子Rmに対応する受信信号として同定する回路である。
遅延処理部1145は、送信イベントに対応して、対象領域設定部1141から受信振動子Rmの位置および対象線BLの位置を示す情報、データ格納部107から受信振動子Rmに対応する受信信号、遅延量算出部1144から各受信振動子Rmに対する受信信号の列に適用する遅延量を入力として取得する。そして、各受信振動子Rmに対応する受信信号の列から、各受信振動子Rmに対する遅延量を差引いた時間に対応する受信信号を代表点Qkからの反射波に基づく受信信号として同定し、加算部1147に出力する。
より具体的には、遅延処理部1145は、各受信振動子Rmに対する受信信号の列に対し、受信振動子Rcと各受信振動子Rmとの受信時間差を打ち消すように遅延処理を行う。遅延処理後の受信信号の列から同一の時刻に対応する受信信号を取り出すことで、同一の代表点Qkからの反射超音波に基づく受信信号の集合を取り出すことができる。
vi)重み算出部1146
重み算出部1146は、各受信振動子Rmに対する重み数列(受信アポダイゼーション)を算出する回路である。
重み数列は送信開口Tx内の各振動子に対応する受信信号に適用される重み係数の数列である。重み数列は、送信フォーカス点Fを中心として対称な分布をなす。重み数列は、送信開口Txの列方向の中心に位置する振動子に対する重みが最大となるように設定され、重みの分布の中心軸は、送信開口中心軸と一致する。重み数列の形状としては、例えば、ハミング窓、ハニング窓、矩形窓である。
重み算出部1146は、対象領域設定部1141から出力される送信開口Txの位置を示す情報を入力として各受信振動子Rmに対する重み数列を算出し、算出した代表点Qkごとの重み数列を加算部1147に出力する。
vii)加算部1147
加算部1147は、遅延処理部1145から出力される各受信振動子Rmに対応して同定された受信信号と、重み算出部1146から出力される重みづけプロファイルを入力として、各受信振動子Rmに対応して同定された受信信号に、各受信振動子Rmに対する重みを乗じて加算して、代表点Qkに対する音響線信号を生成する。遅延処理部1145において送信開口Tx内に位置する各受信振動子Rmが検出した受信信号の位相を整えて加算部1147にて加算処理をすることにより、代表点Qkからの反射波に基づいて各受信振動子Rmで受信した受信信号を重ね合わせてその信号S/N比を増加し、代表点Qkからの受信信号を抽出することができる。
以上の処理をまとめると以下のようになる。受信振動子Rmに対する受信信号列をRf(m,t)とする。ここで、mは受信振動子を指す識別子であり、tは代表点Qkからの反射超音波を受信振動子Rcが受信した時刻である。また、受信振動子Rmに対する重み係数をA(m)とする。さらに、受信振動子Rmと受信振動子Rcとの受信時間差をd(m)とする。このとき、代表点Qkに対する音響線信号Das(k)は、以下の式で与えられることとなる。
Figure 2018110784
以上の処理により、1回の送信イベントに係る対象線BL上の全ての代表点Qkについて音響線信号を生成することができる。生成された線領域の音響線信号は、データ格納部107に出力され保存される。
(3)第2整相加算部1042
第2整相加算部1042は、送信イベントに同期して、被検体内においてサブフレーム音響線信号の生成を行う対象領域Cxを設定する。次に、対象領域Cx上に存する複数の観測点Pij各々について、観測点から各受信振動子Rkが受信した受信信号列を整相加算する。そして、各観測点における音響線信号の列を算出することによりサブフレーム音響線信号を生成する回路である。図8は、第2整相加算部1042の構成を示す機能ブロック図である。図8に示すように、第2整相加算部1042は、対象領域設定部1241、受信開口設定部1242、送信時間算出部1243、受信時間算出部1244、遅延量算出部1245、遅延処理部1246、重み算出部1247、及び加算部1248を備える。
以下、第2整相加算部1042を構成する各部の構成について説明する。なお、第1整相加算部1041を構成する同名称の構成とほぼ同じ機能を有する構成については、詳細な説明を省略する。
i)対象領域設定部1241
対象領域設定部1241は、被検体内においてサブフレーム音響線信号の生成を行う対象領域Cxを設定する。
対象領域設定部1241は、送信イベントに同期して、第1整相加算部1041の対象領域設定部1241が設定した対象領域Bxに基づき、対象領域Cxを設定する。本実施の形態において、対象領域Cxは、対象領域Bxの全域である。
設定された対象領域Cxは送信時間算出部1243、受信時間算出部1244、遅延処理部1246に出力される。
ii)受信開口設定部1242
受信開口設定部1242は、制御部108からの制御信号と、送信ビームフォーマ部103からの送信開口Txの位置を示す情報とに基づき、プローブ101に存する複数の振動子の一部に当たり、列中心が観測点に最も空間的に近接する振動子と合致する振動子列(受信振動子列)を受信振動子として選択して受信開口Rxを設定する回路である。
図9は、受信開口設定部1242により設定された受信開口Rxと送信開口Txとの関係を示す模式図である。図9に示すように、受信開口Rxは、その中心が観測点Pijに最も空間的に近接する振動子Xkと合致するように設定される。なお、送信開口Rxは、観測点Pijごとにその位置が決定され、送信イベントに同期して変動する送信開口Txの位置に基づいては変化しない。すなわち、異なる送信イベントであっても、同一位置にある観測点Pijについての音響線信号を生成する処理においては、同一の受信開口Rx内の受信振動子Rkによって取得された受信信号に基づき整相加算が行われる。
また、超音波主照射領域Ax全体からの反射波を受信するために、受信開口Rxに含まれる振動子の数は、対応する送信イベントにおける送信開口Txに含まれる振動子の数以上に設定することが好ましい。受信開口Rxを構成する振動子列の数は、例えば32、64、96、128、192等としてもよい。
受信開口Rxの設定は、少なくとも列方向における観測点Pijの最大数と同じ回数だけ行われる。また、受信開口Rxの設定は、送信イベントに同期して漸次行われる構成であってもよく、あるいは、全ての送信イベントが終了した後に、各送信イベントに対応した受信開口Rxの設定が送信イベントの回数分まとめて行われる構成であってもよい。
選択された受信開口Rxの位置を示す情報は制御部108を介してデータ格納部107に出力される。
データ格納部107は、受信開口Rxの位置を示す情報と受信振動子に対応する受信信号とを、送信時間算出部1243、受信時間算出部1244、遅延処理部1246、重み算出部1247に出力する。
iii)送信時間算出部1243
送信時間算出部1243は、送信された超音波が被検体中の観測点Pijに到達する送信時間を算出する回路である。送信イベントに対応して、データ格納部107から取得した送信開口Txに含まれる振動子の位置を示す情報と、対象領域設定部1241から取得した対象領域Cxの位置を示す情報とに基づき、対象領域Cx内に存在する任意の観測点Pijについて、送信された超音波が被検体中の観測点Cijに到達する送信時間を算出する。送信時間の算出方法は、送信時間算出部1142と同じであるので、ここでは説明を省略する。
送信時間算出部1243は、1回の送信イベントに対し、対象領域Cx内の全ての観測点Pijについて、送信された超音波が被検体中の観測点Pijに到達する送信時間を算出して遅延量算出部1245に出力する。
iv)受信時間算出部1244
受信時間算出部1244は、観測点Pijからの反射波が、受信開口Rxに含まれる受信振動子Rkの各々に到達する受信時間を算出する回路である。送信イベントに対応して、データ格納部107から取得した受信振動子Rkの位置を示す情報と、対象領域設定部1141から取得した対象領域Cxの位置を示す情報とに基づき対象領域Cx内に存在する任意の観測点Pijについて、送信された超音波が被検体中の観測点Pijで反射され受信開口Rxの各受信振動子Rkに到達する受信時間を算出する。
受信時間算出部1143とは異なり、受信時間算出部1244では、反射波は、観測点Pijから受信振動子Rkへの最短経路403を経由するとして受信時間を算出する。
受信時間算出部1244は、1回の送信イベントに対し、対象領域Cx内に存在する全ての観測点Pijについて、送信された超音波が観測点Pijで反射して各受信振動子Rkに到達する受信時間を算出して遅延量算出部1245に出力する。
v)遅延量算出部1245
遅延量算出部1245は、送信時間と受信時間とから受信開口Rx内の各受信振動子Rkへの総伝播時間を算出し、当該総伝播時間に基づいて、各受信振動子Rkに対する受信信号の列に適用する遅延量を算出する回路である。遅延量算出部1245は、送信された超音波が観測点Pijに到達する送信時間と、観測点Pijで反射して各受信振動子Rkに到達する受信時間を取得する。そして、送信された超音波が各受信振動子Rkへ到達するまでの総伝播時間を算出し、各受信振動子Rkに対する総伝播時間の差異により、各受信振動子Rkに対する遅延量を算出する。遅延量算出部1245は、対象領域Cx内に存在する全ての観測点Pijについて、各受信振動子Rkに対する受信信号の列に適用する遅延量を算出して遅延処理部1246に出力する。
vi)遅延処理部1246
遅延処理部1246は、受信開口Rx内の受信振動子Rkに対する受信信号の列から、各受信振動子Rkに対する遅延量に相当する受信信号を、観測点Pijからの反射超音波に基づく各受信振動子Rkに対応する受信信号として同定する回路である。
遅延処理部1246は、送信イベントに対応して、受信開口設定部1242から受信振動子Rkの位置を示す情報、データ格納部107から受信振動子Rkに対応する受信信号、対象領域設定部1241から取得した対象領域Cxの位置を示す情報、遅延量算出部1245から各受信振動子Rkに対する受信信号の列に適用する遅延量を入力として取得する。そして、各受信振動子Rkに対応する受信信号の列から、各受信振動子Rkに対する遅延量を差引いた時間に対応する受信信号を観測点Pijからの反射波に基づく受信信号として同定し、加算部1248に出力する。
vii)重み算出部1247
重み算出部1247は、受信開口Rxの列方向の中心に位置する振動子に対する重みが最大となるよう各受信振動子Rkに対する重み数列を算出する回路である。
図9に示すように、重み数列は受信開口Rx内の各振動子に対応する受信信号に適用される重み係数の数列である。重み数列は、送信フォーカス点Fを中心として対称な分布をなす。重み数列の分布の形状は、ハミング窓、ハニング窓、矩形窓などを用いることができ、分布の形状は特に限定されない。重み数列は、受信開口Rxの列方向の中心に位置する振動子に対する重みが最大となるように設定され、重みの分布の中心軸は、受信開口中心軸Rxoと一致する。重み算出部1247は、受信開口設定部1242から出力される受信振動子Rkの位置を示す情報を入力として、各受信振動子Rkに対する重み数列を算出し加算部1248に出力する。
viii)加算部1248
加算部1248は、遅延処理部1246から出力される各受信振動子Rkに対応して同定された受信信号を入力として、それらを加算して、観測点Pijに対する整相加算された音響線信号を生成する回路である。あるいは、さらに、重み算出部1247から出力される各受信振動子Rkに対する重み数列を入力として、各受信振動子Rkに対応して同定された受信信号に、各受信振動子Rkに対する重みを乗じて加算して、観測点Pijに対する音響線信号を生成する構成としてもよい。
加算部1248により、送信イベントに同期して対象領域Cx内に存在する全ての観測点Pijに対するサブフレームの音響線信号が生成される。生成されたサブフレームの音響線信号は、データ格納部107に出力され保存される。
(4)パラメータ算出部1043
パラメータ算出部1043は、第1整相加算部1041が生成した線領域音響線信号、および/または、後述する第1合成部1044が生成する第1フレーム音響線信号をもとに、第2整相加算部1042が受信ビームフォーミングを行うためのパラメータを算出する回路である。ここで、パラメータとは、例えば、被検体内の超音波の伝播速度(以下、「超音波速度」という)、観測点の深さに対する重みづけ(以下、「ゲイン補正値」という)、対象領域Cxの大きさを指定する情報などである。
本実施の形態では、被検体内の超音波速度をパラメータとして算出する。被検体内の超音波速度は、音響線信号をBモード画像に変換する際の座標変換(時差を深さの差に変換)や、上述したように、整相加算における遅延量の算出に用いる。遅延量は、図11(a)の模式図に示すように、観測点Pijからの反射波に基づく音響線信号を生成する際、振動子Rcにおける観測点Pijからの反射波に基づく受信信号と、振動子Rmにおける観測点Pijからの反射波に基づく受信信号とのタイミングが同じになるように設定する。このとき、振動子Rcと観測点Pijとの幾何学的距離dc、振動子Rmと観測点Pijとの幾何学的距離dm、超音波速度vuを用いて、遅延量を算出する。具体的には、Δd=dm−dcを用いて、振動子Rmにおける観測点Pijからの反射波に基づく受信信号とのタイミングが、振動子Rcにおける観測点Pijからの反射波に基づく受信信号のタイミングよりΔd/vuだけ遅い、という前提を用いて整相加算を行う。
このとき、超音波速度vuは、被検体の特性から推定した値を用いる。しかしながら、超音波速度は、被検体内の組織の硬さ等に依存し、推定値が誤差を含むことがある。ここで、超音波速度の推定値vuの誤差が大きい場合、すなわち、観測点Pijからの反射波が振動子Rcに到達する時刻と振動子Rmに到達する時刻との真の時差Δtと、計算上の遅延量差Δd/vuとの差が大きい場合、整相加算の受信フォーカスが合わない(ずれた)状態となる。図11(b)は、計算上の超音波速度vuと、遅延量と真の時差とのずれΔd/vu−Δtとの関係を示した模式図である。しかしながら、超音波速度を事前に正確に推定することは困難である。
そこで、パラメータ算出部1043は、同一の対象線BLに対して、用いる計算上の超音波速度vuのみが異なる複数の線領域音響線信号を第1整相加算部1041に生成させ、超音波速度vuの推定を行う。具体的には、パラメータ算出部1043は、例えば、同一の対象線BLに対して生成した、超音波速度v1、v2、v3、v4に対応する線領域音響線信号を比較し、最も分散が大きい線領域音響線信号に対応する超音波速度を、超音波速度として第2整相加算部1042に出力する。これは、計算上の超音波速度と実際の超音波速度とのずれが小さいほど、整相加算の受信フォーカスが合った状態になる、すなわち、音響線信号のS/N比が大きくなるためである。つまり、計算上の超音波速度と実際の超音波速度とのずれが小さいほど、1つの超音波反射源からの反射超音波に基づく複数の受信信号において、ピークのタイミングの時差が小さくなるため、ピークが鋭く、すなわち、ピークの信号値が大きく、かつ、半値幅の時間が短くなり、距離解像度(時間解像度)の高い音響線信号が得られるからである。
パラメータ算出部1043は、算出したパラメータを第2整相加算部1042に出力する。
(5)第1合成部1044
第1合成部1044は、送信イベントに同期して生成される線領域音響線信号からサブフレーム音響線信号を生成し、生成したサブフレーム音響線信号から第1フレーム音響線信号を合成する回路である。図12(a)は、第1合成部1044の構成を示す機能ブロック図である。図12(a)に示すように、第1合成部1044は、音響線信号展開部1341、重みづけ合成部1342を備える。
以下、第1合成部1044を構成する各部の構成について説明する。
i)音響線信号展開部1341
音響線信号展開部1341は、第1フレーム音響線信号を合成するための一連の線領域音響線信号の生成が終了したのち、データ格納部107に保持されている複数の線領域音響線信号を読み出す。そして、線領域音響線信号のそれぞれについて、対象領域Bx内の観測点Sijと代表点Qkとの位置関係に基づいて、線領域音響線信号からサブフレーム音響線信号を生成する。
図13は、音響線信号展開部1341におけるサブフレーム音響線信号を生成する処理を示す模式図である。まず、代表点Qkの深さと、観測点Sijの深さがいずれも送信フォーカス深さより深い場合を想定する。上述の通り、代表点Qkの送信時間は、送信開口Txから送信フォーカス点Fまでの距離と、送信フォーカス点Fから代表点Qkまでの距離との合算値に依存する。すなわち、送信フォーカス点Fから代表点Qkと観測点Sijとが等距離に存在した場合、代表点Qkの送信時間と、観測点Sijの送信時間とは等しい。同様に、代表点Qkの受信時間は、観測点Qkから送信フォーカス点Fまでの距離と、送信フォーカス点Fから受信振動子Rmまでの距離との合算値に依存する。すなわち、送信フォーカス点Fから代表点Qkと観測点Sijとが等距離に存在した場合、代表点Qkの受信時間と、観測点Sijの受信時間とは等しい。したがって、代表点Qkに対する音響線信号は、送信フォーカス点Fまでの距離が等しい複数の観測点Sijに対する音響線信号を含んでいる。言い換えれば、送信フォーカス点Fまでの距離が等しい複数の観測点Sijに対する音響線信号の総和を、代表点Qkに対する音響線信号として取得している。この関係は、代表点Qkの深さと、観測点Sijの深さがいずれも送信フォーカス深さより浅い場合であっても成立する。
したがって、音響線信号展開部1341は、代表点Qkの音響線信号を、(1)送信フォーカス深さに対し、代表点Qkの深さと観測点Sijの深さは、いずれも深い、または、いずれも浅い、(2)代表点Qkと送信フォーカス点Fとの距離と、観測点Sijと送信フォーカス点Fとの距離とは等しい、の2条件を満たす観測点Sijの音響線信号の値として適用する。具体的には、対象領域Bx内に、送信フォーカス点Fを中心として代表点Qkを通過する円弧を設定し、円弧上に存在する全ての観測点Sijに対し、当該観測点Sijに対応する音響線信号の値として、当該円弧上に存在する代表点Qkの音響線信号の値を適用する。このとき、対象領域Bx内で連続しない円弧、すなわち、送信フォーカス点Fからの距離は等しいが、フォーカス深さより浅い円弧と、フォーカス深さより深い円弧、例えば、円弧521と円弧511は、別の円弧として扱う。例えば、円弧514上の全ての観測点Sijに対し、当該当該観測点Sijに対応する音響線信号の値として、円弧514に存在する代表点Qkの音響線信号の値を適用する。このような処理により、音響線信号展開部1341は、線領域音響線信号からサブフレーム音響線信号を生成する。
音響線信号展開部1341は、生成したサブフレーム音響線信号を、重みづけ合成部1342に出力する。
ii)重みづけ合成部1342
図14は、重みづけ合成部1342における合成音響線信号を合成する処理を示す模式図である。上述のとおり、送信イベントに同期して送信振動子列(送信開口Tx)に用いる振動子を振動子列方向に移動ピッチMpだけ異ならせて超音波送受信が順次行われる。そのため、異なる送信イベントに基づく対象領域Bxも送信イベントごとに同一方向に移動ピッチMpだけ位置が異なる。複数のサブフレーム音響線信号を、各サブフレーム音響線信号に含まれる音響線信号と対応する観測点Sijの位置を指標として加算することにより、全ての対象領域Bxを網羅した第1フレーム音響線信号が合成される。
このとき、重みづけ合成部1342は、観測点Sijの位置を指標としてサブフレーム音響線信号に対する重みづけを行う。重み数列は、観測点Sijに対応する各サブフレーム音響線信号に適用される重み係数の数列である。重み数列は、サブフレーム音響線信号に対応する送信イベントにおける送信フォーカス点Fの位置によって規定される。重み数列は、観測点Sijを中心として対称な分布をなす。重み数列は、観測点Sijと同じX座標(振動子の並び方向における位置)に送信フォーカス点Fが設定される送信イベントにおけるサブフレーム音響線信号に対する重みが最大となるように設定される。重みの分布の中心軸は、観測点Pijを通り振動子列に直交する直線Pijoと一致する。重み数列の形状としては、例えば、ハミング窓、ハニング窓、矩形窓である。
重みづけ合成部11402は、各観測点Pijについて、観測点Pijに対応する各サブフレーム音響線信号を重みづけ加算することで、第1フレーム音響線信号を合成する。
また、位置の異なる複数の対象領域Bxにまたがって存在する観測点Pijについては、各サブフレーム音響線信号における音響線信号の値が加算されるので、合成音響線信号は、跨りの程度に応じて大きな値を示す。以後、観測点Pijが異なる対象領域Bxに含まれる回数を「重畳数」、振動子列方向における重畳数の最大値を「最大重畳数」と称する。
また、本実施の形態では、対象領域Bxは砂時計形状の領域内に存在する。そのため、図15(a)に示すように、重畳数及び最大重畳数は被検体の深さ方向において変化するので、合成音響線信号の値も同様に深さ方向において変化する。これを補うため、重みづけ合成部11402は、フレーム音響線信号に含まれる合成音響線信号の合成において、加算が行われた回数に応じて決定した増幅率を各合成音響線信号に乗じる増幅処理を行う。
図15(b)は、重みづけ合成部11402における増幅処理の概要を示す模式図である。図15(b)に示すように、最大重畳数は被検体の深さ方向において変化するので、この変化を補うように、最大重畳数に応じて決定された被検体深さ方向において変化する増幅率が合成音響線信号に乗じられる。これにより、深さ方向における重畳数の変化に伴う合成音響線信号の変動要因は解消され、増幅処理後の合成音響線信号の値は深さ方向において均一化が図られる。
また、重畳数に応じて決定された振動子列方向において変化する増幅率を合成音響線信号に乗じる処理を行ってもよい。振動子列方向において重畳数が変化する場合に、その変動要因を解消し、振動子列方向において増幅処理後の合成音響線信号の値の均一化が図られる。
なお、生成した各観測点に対する合成音響線信号に増幅処理を施した信号をフレーム音響線信号としてもよい。
(6)第2合成部
第2合成部1045は、送信イベントに同期して生成されるサブフレーム音響線信号からフレーム音響線信号を合成する回路である。図12(b)は、第2合成部1045の構成を示す機能ブロック図である。図12(b)に示すように、合成部1045は、加算処理部1343、増幅処理部1344を備える。
以下、第2合成部1045を構成する各部の構成について説明する。
i)加算処理部1343
加算処理部1343は、フレーム音響線信号を合成するための一連のサブフレーム音響線信号の生成が終了したのち、データ格納部107に保持されている複数のサブフレーム音響線信号を読み出す。そして、各サブフレーム音響線信号に含まれる音響線信号が取得された観測点Pijの位置を指標として複数のサブフレーム音響線信号を加算することにより、各観測点に対する合成音響線信号を生成してフレーム音響線信号を合成する。そのため、複数のサブフレーム音響線信号に含まれる同一位置の観測点に対する音響線信号は加算されて合成音響線信号が生成される。
加算処理部1343は、重みづけ合成部1342と同様、複数のサブフレーム音響線信号を、各サブフレーム音響線信号に含まれる音響線信号が取得された観測点Pijの位置を指標として加算することにより、全ての対象領域Cxを網羅した第2フレーム音響線信号を合成する。
ii)増幅処理部1344
上述のとおり、合成音響線信号の値は被検体の深さ方向において変化する。これを補うために、増幅処理部1344は、フレーム音響線信号に含まれる合成音響線信号の合成において、加算が行われた回数に応じて決定した増幅率を各合成音響線信号に乗じる増幅処理を行う。増幅処理部1344では、重みづけ合成部1342と同様の増幅処理を行う。図15(b)に示すように、最大重畳数は被検体の深さ方向において変化するので、この変化を補うように、最大重畳数に応じて決定された被検体深さ方向において変化する増幅率が合成音響線信号に乗じられる。これにより、深さ方向における重畳数の変化に伴う合成音響線信号の変動要因は解消され、増幅処理後の合成音響線信号の値は深さ方向において均一化が図られる。
(7)出力部
出力部1046は、制御部108の制御に基づき、第1合成部1044で生成されるフレーム音響線信号と、第2合成部1045で生成されるフレーム音響線信号とのいずれかを超音波画像生成部105に出力するセレクタ回路である。制御部108は、高機能演算回路1047を使用する場合には、出力部1046に第2合成部1045で生成されるフレーム音響線信号を出力させる。一方、制御部108は、高機能演算回路1047を使用しない場合には、出力部1046に第1合成部1044で生成されるフレーム音響線信号を出力させる。
<動作1:第2整相加算部を用いる場合>
以上の構成からなる超音波診断装置100の動作について説明する。
図16は、第2整相加算部を用いる場合における、受信ビームフォーマ部104のビームフォーミング処理動作を示すフローチャートである。
先ず、ステップS101において、送信部1031は、プローブ101に存する複数の振動子101a中送信開口Txに含まれる各振動子に超音波ビームを送信させるための送信信号を供給する送信処理(送信イベント)を行う。
次に、ステップS102において、受信部1040は、プローブ101での超音波反射波の受信から得た電気信号に基づき受信信号を生成しデータ格納部107に出力し、データ格納部107に受信信号を保存する。プローブ101に存する全ての振動子101aから超音波送信が完了したか否かを判定する(ステップS103)。そして、完了していない場合にはステップS101に戻り、送信開口Txを列方向に移動ピッチMpだけ移動させながら送信イベントを行い、完了している場合にはステップS201に進む。
次に、ステップS201において、対象領域設定部1141は、送信イベントに同期して、送信開口Txの位置を示す情報に基づき対象領域Bxを設定する。1回目のループでは初回の送信イベントにおける送信開口Txから求められる対象領域Bxが設定される。
次に、ステップS202において、対象領域設定部1141は、設定した対象領域Bx内に、対象線BLを設定する。対象線BLは対象領域Bx内部に存在し、かつ、送信フォーカス点Fを通過する直線状の領域である。
次に、ステップS210において、代表点Qkについて音響線信号を生成する。
ここで、ステップS210における、代表点Qkについて音響線信号を生成する動作について説明する。なお、以下では、代表点Qkについて1つだけ音響線信号を生成する動作を説明するが、上述したように、超音波速度のパラメータのみを変更して同一の代表点Qkについて複数の音響線信号を生成する。つまり、超音波速度の値のみを変更して、ステップS210を複数回実行する。
図17は、受信ビームフォーマ部104における代表点Qkについての音響線信号生成動作を示すフローチャートである。図18は、受信ビームフォーマ部104における代表点Qkについての音響線信号生成動作を説明するための模式図である。
まず、ステップS2111において、送信時間算出部1142は、送信された超音波が送信フォーカス点Fに到達する第1時間を算出する。第1時間は、幾何学的に定まる送信開口Txから送信フォーカス点Fに至る経路(401)の長さを超音波の音速csで除することにより、算出できる。
次に、ステップS2112において、送信時間算出部1142は、送信フォーカス点Fから代表点Qkに到達する第2時間を算出する。第2時間は、幾何学的に定まる送信フォーカス点Fから代表点Qkに至る経路(402)の長さを超音波の音速csで除することにより、算出できる。なお、代表点Qkの深さが送信フォーカス深さより浅い場合には、算出した値を絶対値とする負の値を、第2時間とする。すなわち、送信フォーカス点から等距離の2つの代表点Qx、Qyに対し、深さが送信フォーカス深さより深い代表点Qxの第2時間を1.5μsと仮定すると、深さが送信フォーカス深さより浅い代表点Qyの第2時間を−1.5μsとする。
送信時間算出部1142は、第1時間と第2時間の和を代表点Qkに対する送信時間として遅延量算出部1144に出力し、第2時間を受信時間算出部1143に出力する。
次に、送信開口Txから求められる受信振動子Rmの位置を示す座標mを送信開口Tx内の最小値に初期化し(ステップS2114)、送信された超音波が被検体中の代表点Qkで反射され送信開口Txの受信振動子Rmに到達する受信時間を算出する。ここで、代表点Qkで反射された超音波が送信フォーカス点Fまでに到達する時間はすでに第2時間としてステップS2112で算出されている。そこで、受信時間算出部1143は、反射された超音波が送信フォーカス点Fから送信開口Txの受信振動子Rmに到達する第3時間を算出する(ステップS2114)。第3時間は、幾何学的に定まる送信フォーカス点Fから受信振動子Rmまでの経路403の長さを超音波の音速csで除することにより算出できる。そして、受信時間算出部1143は、第2時間と第3時間との和を、受信時間として遅延量算出部1144に出力する。さらに、遅延量算出部1144は、送信時間と受信時間の合計から、送信開口Txから送信された超音波が代表点Qkで反射して受信振動子Rmに到達するまでの総伝播時間を算出し(ステップS2115)、送信開口Tx内の各受信振動子Rmに対する総伝播時間の差異により、各受信振動子Rmに対する遅延量を算出する(ステップS2116)。
送信開口Tx内に存在する全ての受信振動子Rmについて遅延量の算出を完了したか否かを判定し(ステップS2117)、完了していない場合には座標mをインクリメント(ステップS2118)して、更に受信振動子Rmについて遅延量の算出を行い(ステップS2114)、完了している場合にはステップS2121に進む。この段階では、送信開口Tx内に存在する全ての受信振動子Rmについて代表点Qkからの反射波到達の遅延量が算出されている。
ステップS2121において、遅延処理部1145は、送信開口Tx内の受信振動子Rmに対応する受信信号の列に対し、各受信振動子Rmに対する遅延量に基づいた遅延処理を行い、代表点Qkからの反射波に基づく受信信号の時刻(タイミング)を同期させる。上述したように、総伝播時間のうち、第1時間は送信フォーカス点Fと送信開口Txとの位置関係によって一意に決まり、第2時間は送信フォーカス点Fと代表点Qkとの位置関係によって一意に決まり、第3時間は送信フォーカス点Fと受信振動子Rmとの位置関係によって一意に決まる。ここで、1つの送信イベントにおいて送信フォーカス点Fの位置と送信開口Txの位置とはいずれも一定であるから、第1時間は全ての代表点Qk、全ての受信振動子Rmに対して一定である。また、第2時間は受信振動子Rmの位置に依存しないから、1つの受信振動子Rmにおいて代表点Qkと代表点Q(k+1)との総伝播時間の差は、座標mに依存しない。すなわち、同一の受信振動子に対応する受信信号の列において、代表点Qkからの反射波に基づく受信信号と代表点Q(k+1)からの反射波に基づく受信信号との時差は代表点Qkと代表点Q(k+1)との間の距離のみに依存し、受信振動子Rmに対応する信号と受信振動子R(m+1)に対応する信号との間で差はない。一方、第3時間は代表点Qkの位置には依存しないから、受信振動子Rmに対応する受信信号の列における代表点Qkからの反射波に基づく受信信号と、受信振動子R(m+1)に対応する受信信号の列における代表点Qkからの反射波に基づく受信信号との時差は受信振動子Rm、受信振動子R(m+1)、送信フォーカス点Fの3つの位置関係のみに依存し、代表点Qkからの反射波に基づく受信信号と代表点Q(k+1)からの反射波に基づく受信信号との間で差はない。したがって、各受信振動子Rmに対応する受信信号の列のそれぞれに対して第3時間の差を打ち消す遅延処理を行うと、受信信号の列の間で、代表点Qkからの反射波に基づく受信信号の時刻が揃うだけでなく、代表点Q(k+1)からの反射波に基づく受信信号の時刻も、代表点Q(k−1)からの反射波に基づく受信信号の時刻も、それぞれ揃うことになる。したがって、代表点Qkごとに総伝播時間に基づいて受信信号を同定する必要はなく、各受信信号の列に対して第3時間に基づいた遅延処理を行うことで、第1時間と第2時間に基づいて代表点Qkごとの受信信号を振動子列方向の受信信号列として同定することができる。
次に、重み算出部1146は、各受信振動子Rmに対する重み数列を算出する(ステップS2123)。重み算出部1146は、例えば、ハミング窓、ハニング窓、矩形窓を適用する。加算部1147は、各受信振動子Rmに対応して同定された受信信号列に、各受信振動子Rmに対する重みを乗じて加算して、代表点Qkに対する音響線信号を生成し(ステップS2123)、生成された代表点Qkに対応する音響線信号はデータ格納部107に出力され保存される(ステップS2124)。
次に、図16に戻り、ステップS211において、パラメータ算出部1043は、線領域音響線信号に基づいてパラメータを算出する。ここでは、同一の対象線領域BLに対して生成された複数の線領域音響線信号について、それぞれの、音響線信号の値の分散(音響線信号の値の二乗の相加平均)を算出する。そして、最も分散の大きい音響線信号に対応する超音波速度を、被検体内における超音波速度として算出する。算出されたパラメータである超音波速度は、第2整相加算部1042に出力される。
次に、ステップS220において、パラメータを用いて観測点Pijについて音響線信号を生成する。図19は受信ビームフォーマ部104における観測点Pijについての音響線信号生成動作を示すフローチャートである。図20は、受信ビームフォーマ部104における観測点Pijについての音響線信号生成動作を説明するための模式図である。
まず、ステップS2210において、対象領域設定部1241は、観測点Pijの設定範囲である対象領域Cxを設定する。ここでは、対象領域Bx全体を、対象領域Cxとして設定する。
次に、観測点同期型ビームフォーミング処理(ステップS2211〜2264)に進む。まず、観測点Pijの位置を示す座標ijを対象領域Cxの最小値に初期化し(ステップS2211、S2212)、受信開口設定部1242は、列中心が観測点Pijに最も空間的に近接する振動子Xkと合致するよう受信開口Rx振動子列を選択する(ステップS2231)。
次に、観測点Pijについて音響線信号を生成する。
まず、ステップS2241において、送信時間算出部1243は、対象領域Cx上に存在する任意の観測点Pijについて、送信された超音波が被検体中の観測点Pijに到達する送信時間を算出する。このとき、送信時間算出部1243は、経路401と経路402とで定まる幾何学的な経路長と、パラメータ算出部1043が算出した超音波の速度とを用いて送信時間を算出する。
次に受信開口Rxから求められる受信開口Rx内の受信振動子Rkの位置を示す座標kを受信開口Rx内の最小値に初期化し(ステップS2242)、送信された超音波が被検体中の観測点Pijで反射され受信開口Rxの受信振動子Rkに到達する受信時間を算出する(ステップS2243)。受信時間は、幾何学的に定まる観測点Pijから受信振動子Rkまでの経路403の長さを、パラメータ算出部1043が算出した超音波の速度で除することにより算出できる。さらに、送信時間と受信時間の合計から、送信開口Txから送信された超音波が観測点Pijで反射して受信振動子Rkに到達するまでの総伝播時間を算出し(ステップS2244)、受信開口Rx内の各受信振動子Rkに対する総伝播時間の差異により、各受信振動子Rkに対する遅延量を算出する(ステップS2245)。
受信開口Rx内に存在する全ての受信振動子Rkについて遅延量の算出を完了したか否かを判定し(ステップS2246)、完了していない場合には座標kをインクリメント(ステップS2247)して、更に受信振動子Rkについて遅延量の算出を行い(ステップS2243)、完了している場合にはステップS2248に進む。この段階では、受信開口Rx内に存在する全ての受信振動子Rkについて観測点Pijからの反射波到達の遅延量が算出されている。
ステップS2248において、遅延処理部1047は、受信開口Rx内の受信振動子Rkに対応する受信信号の列から、各受信振動子Rkに対する遅延量を差引いた時間に対応する受信信号を観測点Pijからの反射波に基づく受信信号として同定する。
次に、重み算出部1247は、受信開口Rxの列方向の中心に位置する振動子に対する重みが最大となるよう各受信振動子Rkに対する重み数列を算出する(ステップS2249)。加算部1248は、各受信振動子Rkに対応して同定された受信信号に、各受信振動子Rkに対する重みを乗じて加算して、観測点Pijに対する音響線信号を生成し(ステップS2250)、生成された観測点Pijについて音響線信号はデータ格納部107に出力され保存される(ステップS2251)。
次に、座標ijをインクリメントしてステップS2231〜S2251を繰り返すことにより、対象領域Bx内の座標ijに位置する全ての観測点Pij(図20中の「・」)について音響線信号が生成される。対象領域Cx内に存在する全ての観測点Pijについて音響線信号の生成を完了したか否かを判定し(ステップS2261、S2263)、完了していない場合には座標ijをインクリメント(ステップS2262、S2264)して、観測点Pijについて音響線信号を生成し、完了した場合にはステップS230に進む。この段階では、1回の送信イベントに伴う対象領域Cx内に存在する全ての観測点Pijについてのサブフレームの音響線信号が生成され、データ格納部107に出力され保存されている。
次に、全ての送信イベントについて、サブフレームの音響線信号の生成が終了したか否かを判定し(ステップS230)、終了していない場合には、ステップS210に戻り、終了している場合にはステップS301に進む。
次に、ステップS301において、加算処理部1343は、データ格納部107に保持されている複数のサブフレーム音響線信号を読み出し、観測点Pijの位置を指標として複数のサブフレーム音響線信号を加算して各観測点Pijに対する合成音響線信号を生成してフレーム音響線信号を合成する。次に、増幅処理部1344は、フレーム音響線信号に含まれる各合成音響線信号の加算回数に応じて決定された増幅率を各合成音響線信号に乗じ(ステップS302)、増幅されたフレーム音響線信号を超音波画像生成部105及びデータ格納部107に出力し(ステップS303)、処理を終了する。
<動作2:第2整相加算部を用いない場合>
一方、第2整相加算部1042を用いない場合の超音波診断装置100の動作について説明する。
図21は、第2整相加算部を用いない場合における、受信ビームフォーマ部104のビームフォーミング処理動作を示すフローチャートである。なお、図16で示した、第2整相加算部を用いる場合と同様の動作については同じステップ番号を付し、詳細な説明を省略する。
先ず、ステップS101において、送信部1031は、プローブ101に存する複数の振動子101a中送信開口Txに含まれる各振動子に超音波ビームを送信させるための送信信号を供給する送信処理(送信イベント)を行う。
次に、ステップS102において、受信部1040は、プローブ101での超音波反射波に基づき受信信号を生成しデータ格納部107に出力し、データ格納部107に受信信号を保存する。プローブ101に存する全ての振動子101aから超音波送信が完了したか否かを判定し(ステップS103)、完了していない場合にはステップS101に戻り、送信開口Txを列方向に移動ピッチMpだけ移動させながら送信イベントを行い、完了している場合にはステップS201に進む。
次に、ステップS201において、対象領域設定部1141は、送信イベントに同期して、送信開口Txの位置を示す情報に基づき対象領域Bxを設定する。
次に、ステップS202において、対象領域設定部1141は、設定した対象領域Bx内に、対象線BLを設定する。
次に、ステップS210において、代表点Qkについて音響線信号を生成する。
次に、ステップS212において、音響線信号展開部1341は、線領域音響線信号に基づいて、対象領域Bx内の観測点Sijについてサブフレーム音響線信号を生成する。上述したように、音響線信号展開部1341は、観測点Sijについて、送信フォーカス点Fとの距離が同一である代表点Qkのうち、観測点Sijが送信フォーカス深さより深い場合は送信フォーカス深さより深い代表点Qk、観測点Sijが送信フォーカス深さより浅い場合は送信フォーカス深さより浅い代表点Qkを特定し、特定した代表点Qkに対する音響線信号を当該観測点Sijに対する音響線信号とする。この段階では、1回の送信イベントに伴う対象領域Bx内に存在する全ての観測点Sijについてのサブフレームの音響線信号が生成されている。
次に、全ての送信イベントについて、サブフレームの音響線信号の生成が終了したか否かを判定し(ステップS230)、終了していない場合には、ステップS201に戻り、送信開口Txに基づいて対象領域Bxを設定し(ステップS201)、終了している場合にはステップS311に進む。
次に、ステップS311において、重みづけ合成部1342は、観測点Pijに対するサブフレーム音響線信号に対応する送信イベントにおける送信フォーカス点Fの位置に基づいて、サブフレーム音響線信号に対する重みづけ数列を設定する。
次に、ステップS312において、重みづけ合成部1342は、重み数列を用いて、観測点Pijの位置を指標として複数のサブフレーム音響線信号を重みづけ加算して各観測点Pijに対する合成音響線信号を生成してフレーム音響線信号を合成する。次に、重みづけ合成部1342は、フレーム音響線信号を、超音波画像生成部105及びデータ格納部107に出力し(ステップS303)処理を終了する。
<まとめ>
以上、説明したように本実施の形態に係る超音波診断装置100によれば、合成開口法により、異なる送信イベントにより生成された同一位置にある観測点についての音響線信号を重ね合わせて合成する。これにより、複数の送信イベントに対して送信フォーカス点F以外の深度にある観測点においても、仮想的に送信フォーカスを行った効果が得られ空間分解能と信号S/N比を向上することができる。
また、超音波診断装置100では、線領域音響線信号を生成する際に、全ての観測点Pijについて整相加算を行うのではなく、送信フォーカス点Fを中心とする円弧状の領域のそれぞれに代表点Qkを1つずつ設け、代表点Qkに対してのみ整相加算を行う。これにより、整相加算の対象となる代表点Qkの数は対象領域Bxの面積ではなく対象領域Bxの深さ方向の長さに依存するため、整相加算の演算量を大きく削減することができる。さらに、送信時間のみならず受信時間も送信フォーカス点Fを基準とすることで、受信振動子Rmごとの受信時間算出処理を代表点Qkごとに繰り返す必要がなくなる。したがって、総伝播時間の算出処理も整相加算処理も単純化することができ、この点においても整相加算の演算量を大きく削減することができる。一方で、同一の観測点について異なるサブフレーム音響線信号を合成することによる空間解像度、信号S/N比の向上効果を得ることができるため、演算量の低下度合いに対して空間解像度、信号S/N比の低下を抑えることができる。
さらに、超音波診断装置100では、第2整相加算部1042を用いる場合には、線領域音響線信号に基づいて超音波速度などのパラメータを算出し、これをサブフレーム音響線信号の生成に反映する。これにより、整相加算により生成された音響線信号に基づいて整相加算のパラメータを変更する場合、演算量の大きいサブフレーム音響線信号の生成を繰り返す必要がなくなる。さらに、評価に用いる線領域音響線信号はサブフレーム音響線信号と独立に生成されるため、第2整相加算部と第2合成部の演算能力をフレーム音響線信号の生成のみに振り分けることができる。したがって、第2整相加算部と第2合成部の演算能力を最大限に発揮させながら、パラメータの調整を行うことができる。
また、超音波診断装置100では、第2整相加算部1042を用いない場合には、線領域音響線信号に基づいてフレーム音響線信号を生成する。したがって、第1整相加算部と第1合成部の演算能力に見合った低演算量でフレーム音響線信号を生成することができる。そのため、第1整相加算部の演算能力が低い場合でも、空間解像度、信号S/N比の低下を抑えた上で超音波画像を生成することができる。
≪変形例1≫
実施の形態1に係る超音波診断装置100では、パラメータ算出部1043が、超音波の速度を算出した。しかしながら、パラメータ算出部が算出するパラメータは超音波速度に限らず、音響線信号から事後的に最適化できるパラメータであれば、任意のパラメータであってよい。
パラメータの他の具体例としては、例えば、ゲイン補正値が挙げられる。受信部は、A/D変換前に、プローブ101の各振動子101aが超音波反射波の受信で得た電気信号を増幅する「ゲインコントロール」の処理を行う。A/D変換で混入する量子化ノイズのレベルは一定であるため、量子化ノイズの影響を最小限とするためには、増幅後の電気信号の振幅が最大値を超えない範囲で大きく、かつ、揃っていることが好ましい。しかしながら、電気信号が微弱な場合、特にホワイトノイズ等が混入している場合、過度に増幅すると、ホワイトノイズ等も同時に増幅され、整相加算後の音響線信号にノイズが混入し、品質が却って低下する課題がある。そのため、ゲインコントロールでは、S/N比が一定となるよう、増幅することが好ましい。しかしながら、被検体内の超音波減衰率は組織の硬さ、組織の境界の存在等に影響されるため、必ずしも適正に行われるとは限らない。そこで、線領域音響線信号から観測点の深さとS/N比の関係を算出し、深さに対するS/N比の変化を最小にするために、観測点の深さに応じて増幅を行う、としてもよい。その増幅率がゲイン補正値である。具体的には、超音波速度同様、観測点の深さに対する重みづけ係数であるゲイン補正値分布を複数備え、ゲイン補正値分布のみが異なる、同一の対象線領域BLに対する複数の線領域音響線信号を生成し、図22に示すように、深さに対するS/N比の変化が最も小さいゲイン補正値分布g2を第2整相加算部に出力することができる。これにより、S/N比のムラのない高品質な音響線信号を生成することができる。
また、線領域音響線信号ではなく、線領域音響線信号に基づき第1合成部が生成した第1フレーム音響線信号に基づいて、パラメータを算出してもよい。パラメータの具体例としては、例えば、対象領域Cxの幅があげられる。具体的には、第1フレーム音響線信号の時間変化に基づいて被検体内の組織および/またはプローブの移動を検知し、図23に示すように、被検体内の組織および/またはプローブの移動量が大きくなるほど、対象領域Cxの幅Wが小さくなるようにパラメータを算出する。被検体内の組織および/またはプローブの移動は、例えば、連続する2フレームの第1フレーム音響線信号で相関処理を行い、相関値が閾値を下回るか否かで検出することができる。これにより、プローブの移動時に重畳数を小さくし、モーションアーチファクトの発生を抑制することができる。また、例えば、被検体内の組織および/またはプローブの移動の大きさに対して、対象領域Cxの幅を多段階で狭めてもよい。具体的には、相関値が第1の閾値を下回るときは、対象領域Cxの幅Wを対象領域Bxの幅Txの1/2とし、相関値が第2の閾値を下回るときは、対象領域Cxの幅Wを対象領域Bxの幅Txの1/3とする、としてもよい。
なお、パラメータは上述の例に限定されず、線領域音響線信号または第1フレーム音響線信号から算出可能な任意のパラメータであってよい。また、超音波速度やゲイン補正値を、第1フレーム音響線信号から算出するとしてもよい。
≪実施の形態2≫
実施の形態1および変形例では、第2整相加算部を使用するか否かで超音波診断装置の構成が変化しない場合について説明した。しかしながら、第2整相加算部と第2合成部とを構成する高機能演算回路を脱着可能に構成し、超音波診断装置は、高機能演算回路が装着されているときにのみ第2整相加算部を使用する、としてもよい。
図23(a)は、実施の形態2に係る超音波診断システム2000を示す外観図である。図23(a)に示すように、超音波診断システム2000は、可搬本体部(HCU;Hand Carry Unit)200、プローブ201、カート203、増設部204、入力部205を有する。表示部202は、可搬本体部200に内蔵されている。
可搬本体部200は、プローブ201に超音波の送受信を行わせプローブ101からの出力信号に基づき超音波画像を生成して表示部202に表示する超音波診断装置であり、超音波診断装置100のうち高機能演算回路1047、すなわち、第2整相加算部1042と第2合成部1045とを除く全構成を有している。可搬本体部200は、例えば、PCI Express×16対応のスイッチファブリックと、スイッチファブリックに接続されたパーソナルコンピュータと、スイッチファブリックに接続されFPGAで構成された受信部1040とで実現される。なお、スイッチファブリックの通信仕様はPCI Express×16に限られず、受信部1040とパーソナルコンピュータ、パーソナルコンピュータと後述する増設部204との通信に十分な速度(帯域)を有していれば、任意のバスを用いてよい。
カート203は可搬本体部200を搭載した状態で移動可能な台であり、天板221、増設部204、入力部205を備える。
増設部204は、高機能演算回路1047を備える。具体的には、GPUと、GPUを動作させるための電源ユニットなどから構成される。なお、具体的構成はGPUに限られず、高機能な演算能力を有するプロセッサや、プロセッサや演算回路を含むモジュールであってもよい。
入力部205は、例えば、可搬本体部200と接続可能なキーボードである。
可搬本体部200は、カート203から取り外されているときは、図20のフローチャートで示したように、第2整相加算部を用いない動作を行う。一方、可搬本体部200は、カート203に取り付けられているときは、図16のフローチャートで示したように、第2整相加算部を用いる動作を行う。
<脱着構成>
可搬本体部200の裏面は脱着部を備えている。脱着部は、例えば、図23(b)に示すように、X方向に延びる2つの係合溝211からなる。係合溝211は、L字型の溝であり、可搬本体部200の背面側には、係止爪212が設けられている。係止爪212はばねで係止溝内に付勢されており、可搬本体部200の背面側が傾斜面となっている。そのため、係合溝211に後述する係合凸部222の挿入を許容するが、係合溝211の奥まで挿入された係合凸部222の引き出しを抑止する構成となっている。可搬本体部200の背面には、係止爪212を係止溝から引き上げるための解除釦213と、コネクタ214が設けられている。コネクタ214には、スイッチファブリックから延びる信号線が格納されている。なお、コネクタ214は、いわゆる活線挿抜(ホットプラグ)に対応していることが好ましい。さらに、コネクタ214は、入力部205からの入力を受け付ける信号線を含んでいてもよい。
一方、カート203の天板221は脱着部を備えている。脱着部は、例えば、図23(c)に示すように、X方向に延びる2つの係合凸部222とコネクタ223が設けられている。係合凸部222は、係合溝211と係合可能なL字型の突起である。なお、図23(c)に示すように、係合凸部222はカート203の前面側に近づくほど、Y方向に延びる部分が小さくなったり、Y方向に延びる部分がない状態でZ方向の高さが低くなったりしてもよい。このような形状にすることで、係合凸部222の一部のみが係合溝211に挿入された状態で可搬本体部200が天板221に対してZ方向に相対移動可能となるため、可搬本体部200の着脱が容易となる。コネクタ223は、天板221の係合凸部222と可搬本体部200の係合溝211が係合して係止爪212によって可搬本体部200が天板221に対して相対移動できない状態で、コネクタ214と係合するように設けられている。コネクタ223には、増設部204の高機能演算回路1047から延びる信号線が格納されている。さらに、コネクタ223は、入力部205からの信号線を含んでいてもよい。
可搬本体部200をカート203に取り付ける際には、係合凸部222のカート203の前面側の端部と係合溝211の位置を合わせ、X方向カート後方側に可搬本体部200をスライドさせる。これにより、係合溝211がガイドとなって、係合凸部222と係合溝211が係合する。係合凸部222が係合溝211に完全に格納されると、コネクタ214とコネクタ223が係合し、可搬本体部200と増設部204が電気的に接続される。また、係止爪212が係合溝211に突出し、カート203に対して可搬本体部200が相対移動することを抑止する。これにより、可搬本体部200がカート203に固定される。
一方、可搬本体部200をカート203から取り外す際には、解除釦213を押しながら可搬本体部200をカート203の前面側にスライドさせ、さらに、可搬本体部200をカート203の前面側に引き出す。解除釦213を押すことで係合凸部222と係合溝211がX方向に相対移動可能となり、可搬本体部200をカート203の前面側にスライドさせることでコネクタ214とコネクタ223との接続が解除される。さらに、可搬本体部200をカート203の前面側にスライドさせることで、係合凸部222と係合溝211の係合が完全に解かれる。
<機能構成>
図24に、可搬本体部200の受信ビームフォーマ部304の構成を示す。図24に示すように、受信ビームフォーマ部304は、受信ビームフォーマ部104から高機能演算回路1047、すなわち、第2整相加算部1042と第2合成部1045とを除く構成である。
受信ビームフォーマ部304は、それ単体で<動作2:第2整相加算部を用いない場合>の動作を行う。さらに、可搬本体部200がカート203に固定されているときは、増設部204の高機能演算回路1047と連携する。すなわち、受信ビームフォーマ部304と高機能演算回路1047との組み合わせで、受信ビームフォーマ部104を実現する。
可搬本体部200の制御部は、可搬本体部200がカート203に固定されているか否かを判定し、固定されていない場合は、第2整相加算部を用いないと決定し、固定されている場合は、第2整相加算部を用いると決定する。可搬本体部200がカート203に固定されているか否かを判定する方法としては、例えば、スイッチファブリックと高機能演算回路1047との接続が有効か否かを検出する方法がある。なお、判定方法はこの場合に限られず、例えば、制御部が高機能演算回路1047への通信を試みて、成功するか否かで判定してもよい。または、コネクタ214に、コネクタ223と接続されているか否かを検出する検出信号線を設け、検出信号線の状態で判定してもよい。または、コネクタ214、係合溝211、あるいは、可搬本体部200の裏面または背面に、可搬本体部200がカート203に固定されているか否かを検出するためのセンサを設けてもよい。
<まとめ>
実施の形態2に係る超音波診断装置によれば、可搬本体部のみで超音波画像を生成する場合には、線領域音響線信号を生成して線領域音響線信号から超音波画像を生成する、という演算量の小さい処理を行う。したがって、可搬本体部が高い演算能力を有している必要がない。また、演算量が小さいことで省電力化が可能となるため、可搬本体部の排熱設計や回路の単純化等による可搬本体部の小型化や、内蔵蓄電池による稼働可能時間の延長が容易となる。一方、可搬本体部がカートに搭載されている場合は、線領域音響線信号に基づいて最適化されたパラメータを用いて、高機能演算回路が合成開口法により超音波画像を生成する。したがって、高機能演算回路の演算能力を超音波画像の生成のみに使用した上でパラメータの最適化が可能となり、高品質の超音波画像を生成した上で、対象領域の広域化やフレームレートの向上が可能となる。そのため、可搬本体部のみでも超音波診断装置として使用可能でありながら、カートと一体化した場合には高品質な超音波診断装置として使用することができる。
≪実施の形態に係るその他の変形例≫
(1)各実施の形態および変形例では、第1整相加算部は、対象領域Bxの内部であり、かつ、フォーカス点Fを通過する対象線BLについて線領域音響線信号を生成し、第1合成部は、線領域音響線信号に基づいて第1フレーム音響線信号を生成する場合について説明した。しかしながら、第1整相加算部、第1合成部の受信ビームフォーミングは上述の例に限られない。例えば、第1整相加算部は、フォーカス点Fまたはその近傍を通過し、振動子列101aに直交する1または複数の直線状の領域BLxについて、従来型の受信ビームフォーミングにより線領域音響線信号を生成するとしてもよい。ここで、従来型の受信ビームフォーミングとは、例えば、観測点Qと送信開口Txとの最短距離(すなわち、観測点Qの深さ)を超音波速度で除した値を送信時間、観測点Qと受信振動子Rとの距離を超音波速度で除した値を受信時間、とする、既知のビームフォーミングである。このような構成においても、実施の形態と同様の効果を得ることができる。
(2)各実施の形態および変形例では、第2整相加算部は、観測点同期型の受信ビームフォーミングを行うとした。しかしながら、第2整相加算部は、送信開口同期型の受信ビームフォーミングを行うとしてもよい。この場合、受信開口を、送信開口の位置に基づいて設定する。具体的には、送信開口の中心軸と受信開口の送信軸が一致するように、受信開口を設ける。したがって、受信開口は観測点の位置には依存せず、送信開口の位置に依存する。そのため、送信イベントに同期してそれぞれ異なる受信開口にて整相加算を行うことができ、複数の送信イベントにわたって受信時刻は異なるものの、結果としてより一層広い受信開口を用いた受信処理の効果が得られ、広い観測領域で空間分解能を均一にすることができるという効果を奏する。
(3)実施の形態2では、可搬本体部200とカート203のそれぞれの脱着部の具体的形状について説明した。しかしながら、脱着部の構成は上述の例に限られない。着脱部の形状は、(a)可搬本体部がユーザの意に反してカートから離脱しない、(b)可搬本体部がカートに装着されているときに、可搬本体部と増設部とが電気的に接続される、の2条件を満たしていれば、任意の形状でよい。なお、ユーザが可搬本体部を取り外そうとしているときにカートからの離脱が容易であることが好ましい。
また、実施の形態2では、可搬本体部200とカート203のそれぞれにコネクタを設ける、としたが、例えば、可搬本体部と増設部とは、無線通信または非接触通信を行うとしてもよい。
(4)各実施の形態および変形例では、パラメータ算出部がパラメータを算出した送信イベントにおいて、第2整相加算部が当該パラメータを用いて整相加算する場合について説明した。しかしながら、パラメータの反映は、同一の送信イベントである必要はない。例えば、第2整相加算部は、パラメータ算出部からパラメータを受信すると、次以降の送信イベントにおいてパラメータを反映するとしてもよい。この構成により、第2整相加算部は、第1整相加算部による線領域音響線信号の生成およびパラメータ算出部によるパラメータの算出を待つ必要がなく、第2整相加算部による整相加算と、第1整相加算部による線領域音響線信号の生成およびパラメータ算出部によるパラメータの算出とを並列処理で行うことができる。したがって、第1整相加算部およびパラメータ算出部による演算能力が第2整相加算部および第2合成部によるフレーム音響線信号の生成時間に影響を与えず、第2整相加算部を使用する場合におけるフレームレートを向上させることができる。
(5)なお、本発明を上記実施の形態に基づいて説明してきたが、本発明は、上記の実施の形態に限定されず、以下のような場合も本発明に含まれる。
例えば、上記超音波診断装置の全部、もしくは一部を、マイクロプロセッサ、ROM、RAM等の記録媒体、ハードディスクユニットなどから構成されるコンピュータシステムで構成した場合も本発明に含まれる。上記RAM又はハードディスクユニットには、上記各装置と同様の動作を達成するコンピュータプログラムが記憶されている。上記マイクロプロセッサが、上記コンピュータプログラムにしたがって動作することにより、各装置はその機能を達成する。
また、上記の各装置を構成する構成要素の一部又は全部は、1つのシステムLSI(Large Scale Integration(大規模集積回路))から構成されているとしてもよい。システムLSIは、複数の構成部を1個のチップ上に集積して製造された超多機能LSIであり、具体的には、マイクロプロセッサ、ROM、RAMなどを含んで構成されるコンピュータシステムである。これらは個別に1チップ化されてもよいし、一部又は全てを含むように1チップ化されてもよい。なお、LSIは、集積度の違いにより、IC、システムLSI、スーパーLSI、ウルトラLSIと呼称されることもある。上記RAMには、上記各装置と同様の動作を達成するコンピュータプログラムが記憶されている。上記マイクロプロセッサが、上記コンピュータプログラムにしたがって動作することにより、システムLSIは、その機能を達成する。例えば、本発明のビームフォーミング方法がLSIのプログラムとして格納されており、このLSIがコンピュータ内に挿入され、所定のプログラム(ビームフォーミング方法)を実施する場合も本発明に含まれる。
なお、集積回路化の手法はLSIに限るものではなく、専用回路または汎用プロセッサで実現してもよい。LSI製造後に、プログラムすることが可能なFPGA(Field Programmable Gate Array)や、LSI内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサー(Reconfigurable Processor)を利用してもよい。
さらには、半導体技術の進歩または派生する別技術によりLSIに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行ってもよい。
また、各実施の形態に係る、超音波診断装置の機能の一部又は全てを、CPU等のプロセッサがプログラムを実行することにより実現してもよい。上記超音波診断装置の診断方法や、ビームフォーミング方法を実施させるプログラムが記録された非一時的なコンピュータ読み取り可能な記録媒体であってもよい。プログラムや信号を記録媒体に記録して移送することにより、プログラムを独立した他のコンピュータシステムにより実施するとしてもよい、また、上記プログラムは、インターネット等の伝送媒体を介して流通させることができるのは言うまでもない。
上記実施形態に係る超音波診断装置では、記憶装置であるデータ格納部を超音波診断装置内に含む構成としたが、記憶装置はこれに限定されず、半導体メモリ、ハードディスクドライブ、光ディスクドライブ、磁気記憶装置、等が、超音波診断装置に外部から接続される構成であってもよい。
また、ブロック図における機能ブロックの分割は一例であり、複数の機能ブロックを一つの機能ブロックとして実現したり、一つの機能ブロックを複数に分割したり、一部の機能を他の機能ブロックに移してもよい。また、類似する機能を有する複数の機能ブロックの機能を単一のハードウエア又はソフトウェアが並列又は時分割に処理してもよい。
また、上記のステップが実行される順序は、本発明を具体的に説明するために例示するためのものであり、上記以外の順序であってもよい。また、上記ステップの一部が、他のステップと同時(並列)に実行されてもよい。
また、超音波診断装置には、プローブ及び表示部が外部から接続される構成としたが、これらは、超音波診断装置内に一体的に具備されている構成としてもよい。
また、上記実施の形態においては、プローブは、複数の圧電素子が一次元方向に配列されたプローブ構成を示した。しかしながら、プローブの構成は、これに限定されるものではなく、例えば、複数の圧電変換素子を二次元方向に配列した二次元配列振動子や、一次元方向に配列された複数の振動子を機械的に揺動させて三次元の断層画像を取得する揺動型プローブを用いてもよく、測定に応じて適宜使い分けることができる。例えば、2次元に配列されたプローブを用いた場合、圧電変換素子に電圧を与えるタイミングや電圧の値を個々に変化させることによって、送信する超音波ビームの照射位置や方向を制御することができる。
また、プローブは、送受信部の一部の機能をプローブに含んでいてもよい。例えば、送受信部から出力された送信電気信号を生成するための制御信号に基づき、プローブ内で送信電気信号を生成し、この送信電気信号を超音波に変換する。併せて、受信した反射超音波を受信電気信号に変換し、プローブ内で受信電気信号に基づき受信信号を生成する構成を採ることができる。
また、各実施の形態に係る超音波診断装置、及びその変形例の機能のうち少なくとも一部を組み合わせてもよい。更に上記で用いた数字は、全て本発明を具体的に説明するために例示するものであり、本発明は例示された数字に制限されない。
さらに、本実施の形態に対して当業者が思いつく範囲内の変更を施した各種変形例も本発明に含まれる。
≪まとめ≫
(1)実施の形態に係る超音波診断装置は、複数の振動子を備えた超音波プローブを用いて被検体に集束型の超音波ビームを送信する送信イベントを複数回繰り返すとともに、各送信イベントに同期して被検体から反射超音波を受波して受信信号列を生成し、前記受信信号列に基づいて超音波画像を生成する超音波診断装置であって、超音波ビームの集束位置を規定するフォーカス点を送信イベントごとに変更しながら、送信イベントごとに、前記超音波プローブに列設された複数の振動子から送信振動子列を選択して当該送信振動子列から被検体内の対象領域に超音波ビームを送信させる送信部と、各送信イベントに同期して、前記超音波プローブが前記対象領域から受波した反射超音波に基づいて、振動子各々に対する受信信号列を生成する受信部と、前記送信イベントごとに、少なくとも前記対象領域の一部を含む第1対象領域内の複数の観測点について、前記受信信号列を整相加算し第1音響線信号を生成する第1整相加算部と、前記第1音響線信号に基づいて、サブフレーム音響線信号を生成するためのパラメータを算出するパラメータ算出部と、前記対象領域の全部または一部である第2対象領域内の複数の観測点について、前記パラメータに基づき、前記受信信号列を整相加算しサブフレーム音響線信号を生成する第2整相加算部と、前記サブフレーム音響線信号を合成してフレーム音響線信号を生成する合成部と、前記第1音響線信号と前記フレーム音響線信号のいずれに基づいて超音波画像を生成するかを決定する制御部と、前記制御部の決定に基づき、前記第1音響線信号と前記フレーム音響線信号のいずれか一方から超音波画像を生成する超音波画像生成部とを含む超音波信号処理回路を備えることを特徴とする。
上記構成により、必要な演算能力に応じて、第2整相加算部を用いる動作と、第2整相加算部を用いない動作とを切り替えることができる。

(2)また、上記(1)の超音波診断装置は、前記超音波信号処理回路は、少なくとも前記第1整相加算部と前記パラメータ算出部と前記制御部とを含む第1回路と、少なくとも前記第2整相加算部と前記合成部とを含む第2回路とを含み、前記第2回路は前記第1回路に対して脱着可能に構成され、前記制御部は、前記第2回路が前記第1回路に接続されているときは、前記フレーム音響線信号に基づいて超音波画像を生成すると決定し、前記第2回路が前記第1回路から切り離されているときは、前記第1音響線信号に基づいて超音波画像を生成すると決定する、としてもよい。
上記構成により、演算能力より可搬性を重視する場合には演算量の小さい受信ビームフォーミングにより超音波画像を生成することができる一方、可搬性より演算能力を重視する場合には高品質の超音波画像を生成することができる。したがって、使用状況に応じて最適な制御を行うことができる。
(3)また、上記(1)〜(2)の超音波診断装置は、前記パラメータ算出部は、前記第1音響線信号に基づいて前記被検体内の動きを検出し、前記制御部は、被検体内の動きが大きいほど、前記超音波プローブの振動子の並ぶ方向における第2対象領域の幅が小さくなるように、前記第2整相加算部に指示する、としてもよい。
(4)また、上記(3)の超音波診断装置は、前記パラメータ算出部は、前記第1音響線信号から生成したフレーム音響線信号と、他のフレームにおける第1音響線信号から生成したフレーム音響線信号とから、前記被検体内の動きを検出する、としてもよい。
これら上記構成により、第2整相加算部の演算に影響を与えずにプローブや被検体の動きを検出し、速やかにモーションアーチファクトの発生の抑制を行うことができる。
(5)また、上記(1)〜(4)の超音波診断装置は、前記第1整相加算部は、整相加算に用いる被検体内の超音波速度について複数の推定値を保持しており、前記送信イベントごとに、同一の第1対象領域に対して、前記複数の推定値のそれぞれを用いて、複数の第1音響線信号を生成し、前記パラメータ算出部は、前記複数の第1音響線信号に基づいて、被検体内の超音波速度を推定し、前記制御部は、前記パラメータ算出部が推定した被検体内の超音波速度を用いて整相加算を行うよう、第2整相加算部に指示する、としてもよい。
(6)また、上記(5)の超音波診断装置は、前記パラメータ算出部は、第1対象領域が同一で超音波速度の推定値が異なる、複数の第1音響線信号のそれぞれに対して、信号の値の分散を算出し、算出した分散が最も大きい第1音響線信号に対応する推定値を、被検体内の超音波速度であると推定する、としてもよい。
これら上記構成により、第2整相加算部の演算に影響を与えることなく、超音波速度を最適化し、超音波画像の解像度を向上させることが可能となる。
(7)また、上記(1)〜(6)の超音波診断装置は、前記パラメータ算出部は、前記第1音響線信号における、観測点の深さと信号のS/N比との関係を算出して、観測点の深さに応じた増幅率を決定し、前記制御部は、前記パラメータ算出部が決定した、観測点の深さに応じた増幅率を用いて、第2音響線信号に重みづけを行うよう、第2整相加算部に指示する、としてもよい。
(8)また、上記(7)の超音波診断装置は、前記第1整相加算部は、観測点の深さに対する信号増幅率について複数のプロファイルを保持しており、前記送信イベントごとに、同一の第1対象領域に対して、前記複数のプロファイルのそれぞれを用いて、複数の第1音響線信号を生成し、前記パラメータ算出部は、前記第1音響線信号における、観測点の深さに対する信号のS/N比の分散が最も小さいプロファイルを、観測点の深さに応じた増幅率として決定する、としてもよい。
これら上記構成により、第2整相加算部の演算に影響を与えることなくゲインコントロールを行い、超音波画像のS/N比を均一化させることが可能となる。
(9)また、上記(1)〜(8)の超音波診断装置は、前記第2整相加算部は、送信された超音波が各観測点に到達する送信時間を、前記観測点の深さが前記被検体内で超音波が集束するフォーカス深さ以上である場合には、送信された超音波が前記送信振動子列から前記フォーカス点に到達するまでの第1時間と、前記基準点から前記観測点に到達する第2時間との合計を送信時間として算出し、前記観測点の深さが前記被検体内で超音波が集束するフォーカス深さ未満である場合には、前記第1時間から前記第2時間を減算した結果を送信時間として算出する、としてもよい。
上記構成により、第2整相加算部を用いる場合に高解像度の超音波画像を生成することが可能となる。
(10)また、上記(1)〜(9)の超音波診断装置は、前記画像生成部は、前記第1音響線信号からフレーム音響線信号を生成する第1音響線信号合成部を備え、前記第1対象領域は、前記フォーカス点を通過し、全域が前記対象領域内に含まれる直線状の領域であり、前記第1整相加算部は、前記送信振動子列に含まれる振動子のそれぞれに対応する受信信号列に対して、観測点と前記フォーカス点との距離と、前記フォーカス点と振動子との距離とをそれぞれ被検体内の超音波速度で除算した値に基づく遅延処理を含む整相加算を行って前記第1音響線信号を生成し、前記画像生成部は、前記第1音響線信号から超音波画像を生成する場合、前記第1音響線信号合成部が生成したフレーム音響線信号から超音波画像を生成する、としてもよい。
(11)また、上記(10)の超音波診断装置は、前記第1整相加算部は、送信された超音波が各観測点に到達する送信時間を、前記観測点の深さが前記被検体内で超音波が集束するフォーカス深さ以上である場合には、送信された超音波が前記送信振動子列から前記フォーカス点に到達するまでの第1時間と、前記基準点から前記観測点に到達する第2時間との合計を送信時間として算出し、前記観測点の深さが前記被検体内で超音波が集束するフォーカス深さ未満である場合には、前記第1時間から前記第2時間を減算した結果を送信時間として算出する、としてもよい。
これら上記構成により、第1整相加算部を用いない場合に、演算量を大きく削減しながらも品質劣化を抑えた超音波画像を生成することができる。また、第2整相加算部を用いる場合、特に超音波速度をパラメータとして算出する際に、誤差を低減させることができる。
(12)また、上記(10)〜(11)の超音波診断装置は、前記第1整相加算部は、各観測点からの反射波が各振動子に到達する受信時間を、前記観測点から前記観測点に最近接する振動子に到達するまでの時間を、前記観測点に最近接する振動子に対応する受信時間として算出し、超音波が前記フォーカス点から振動子に到達するまでの時間と前記フォーカス点から前記観測点に最近接する振動子に到達するまでの時間との差を、前記観測点に最近接する振動子に対応する受信時間に加算することで、当該振動子に対応する受信時間として算出する、としてもよい。
上記構成により、観測点ごとに受信信号を同定する必要がなくなり、第1整相加算部の演算量を大きく削減することができる。
(13)また、上記(10)〜(12)の超音波診断装置は、前記第1音響線信号合成部は、前記送信イベントごとに、前記対象領域内の各観測点について、前記フォーカス点との距離が当該観測点と同一であり前記直線上に存在する観測点の前記第1音響線信号を割り当てて当該観測点の音響線信号とすることで、サブフレーム音響線信号を生成し、生成した複数のサブフレーム音響線信号を合成してフレーム音響線信号を生成する、としてもよい。
上記構成により、対象領域内の観測点の数に対して第1整相加算部の演算量を削減することができる。
(14)また、上記(1)〜(9)の超音波診断装置は、前記第1対象領域は、前記フォーカス点またはその近傍を通過し、前記超音波プローブの振動子の並ぶ方向に直交する1以上の直線からなり、前記第1整相加算部は、観測点の深さを被検体内の超音波速度で除算した値を送信時間、観測点から振動子までの距離を被検体内の超音波速度で除算した値を受信時間として整相加算を行い、前記第1音響線信号を生成し、前記画像生成部は、前記第1音響線信号から超音波画像を生成する場合、前記第1整相加算部が生成した複数の第1音響線信号から超音波画像を生成する、としてもよい。
上記構成により、第1整相加算部を用いない場合に、演算量を削減することができる。
(15)また、実施の形態の他の態様に係る超音波診断装置は、複数の振動子を備えた超音波プローブを用いて被検体に集束型の超音波ビームを送信する送信イベントを複数回繰り返すとともに、各送信イベントに同期して被検体から反射超音波を受波して受信信号列を生成し、前記受信信号列に基づいて超音波画像を生成する超音波診断装置であって、超音波ビームの集束位置を規定するフォーカス点を送信イベントごとに変更しながら、送信イベントごとに、前記超音波プローブに列設された複数の振動子から送信振動子列を選択して当該送信振動子列から被検体内の対象領域に超音波ビームを送信させる送信部と、各送信イベントに同期して、前記超音波プローブが前記対象領域から受波した反射超音波に基づいて、振動子各々に対する受信信号列を生成する受信部と、前記送信イベントごとに、少なくとも前記対象領域の一部を含む第1対象領域内の複数の観測点について、前記受信信号列を整相加算し第1音響線信号を生成する第1整相加算部と、前記第1音響線信号に基づいて、サブフレーム音響線信号を生成するためのパラメータを算出するパラメータ算出部と、前記第1音響線信号から超音波画像を生成する超音波画像生成部と、制御部とを含む超音波信号処理回路を備え、前記対象領域の全部または一部である第2対象領域内の複数の観測点について、前記パラメータに基づき、前記受信信号列を整相加算しサブフレーム音響線信号を生成する第2整相加算部と、前記サブフレーム音響線信号を合成してフレーム音響線信号を生成する合成部とを含む演算回路が前記超音波信号処理回路に対して着脱可能に構成され、前記制御部は、前記演算回路が前記超音波信号処理回路に接続されているとき、前記超音波画像生成部に、前記第1音響線信号に替えて前記フレーム音響線信号から超音波画像を生成させることを特徴とする。
上記構成により、第1整相加算部が生成する第1音響線信号に基づき、演算量の小さい受信ビームフォーミングにより超音波画像を生成することができる。一方、高い演算能力が必要な場合は、第2整相加算部と合成部とを装着することで、高品質の超音波画像を生成することができる。さらに、第2整相加算部に対し第1音響線信号に基づいてパラメータ制御を行うことで、演算量の大きい第2整相加算部による受信ビームフォーミングを不必要に繰り返すことなく、第1整相加算部による低演算量の処理によりパラメータの最適化を図ることができる。
本開示にかかる超音波診断装置は、ポータブルな超音波診断装置としても高演算能力の超音波診断装置としても使用可能であり、かつ、使用状況に応じて最適な動作を行うことができる。
100 超音波診断装置
101、201 プローブ
101a 振動子
102 マルチプレクサ部
103 送信ビームフォーマ部
1031 送信部
104、304 受信ビームフォーマ部
1040 受信部
1041 第1整相加算部
1042 第2整相加算部
1043 パラメータ算出部
1044 第1合成部
1045 第2合成部
1046 出力部
1047 高機能演算回路
1141、1241 対象領域設定部
1242 受信開口設定部
1142、1243 送信時間算出部
1143、1244 受信時間算出部
1144、1245 遅延量算出部
1145、1246 遅延処理部
1146、1247 重み算出部
1147、1248 加算部
1341 音響線信号展開部
1342 重み付け合成部
1343 加算処理部
1344 増幅処理部
105 超音波画像生成部
106、202 表示部
107 データ格納部
108 制御部
150 超音波信号処理装置
200 可搬本体部
203 カート
204 増設部
205 入力部
211 係合溝
212 係止爪
213 解除釦
214、223 コネクタ
221 天板
222 係合凸部

Claims (15)

  1. 複数の振動子を備えた超音波プローブを用いて被検体に集束型の超音波ビームを送信する送信イベントを複数回繰り返すとともに、各送信イベントに同期して被検体から反射超音波を受波して受信信号列を生成し、前記受信信号列に基づいて超音波画像を生成する超音波診断装置であって、
    超音波ビームの集束位置を規定するフォーカス点を送信イベントごとに変更しながら、送信イベントごとに、前記超音波プローブに列設された複数の振動子から送信振動子列を選択して当該送信振動子列から被検体内の対象領域に超音波ビームを送信させる送信部と、
    各送信イベントに同期して、前記超音波プローブが前記対象領域から受波した反射超音波に基づいて、振動子各々に対する受信信号列を生成する受信部と、
    前記送信イベントごとに、少なくとも前記対象領域の一部を含む第1対象領域内の複数の観測点について、前記受信信号列を整相加算し第1音響線信号を生成する第1整相加算部と、
    前記第1音響線信号に基づいて、サブフレーム音響線信号を生成するためのパラメータを算出するパラメータ算出部と、
    前記対象領域の全部または一部である第2対象領域内の複数の観測点について、前記パラメータに基づき、前記受信信号列を整相加算しサブフレーム音響線信号を生成する第2整相加算部と、
    前記サブフレーム音響線信号を合成してフレーム音響線信号を生成する合成部と、
    前記第1音響線信号と前記フレーム音響線信号のいずれに基づいて超音波画像を生成するかを決定する制御部と、
    前記制御部の決定に基づき、前記第1音響線信号と前記フレーム音響線信号のいずれか一方から超音波画像を生成する超音波画像生成部と
    を含む超音波信号処理回路を備えることを特徴とする超音波診断装置。
  2. 前記超音波信号処理回路は、少なくとも前記第1整相加算部と前記パラメータ算出部と前記制御部とを含む第1回路と、少なくとも前記第2整相加算部と前記合成部とを含む第2回路とを含み、
    前記第2回路は前記第1回路に対して脱着可能に構成され、
    前記制御部は、前記第2回路が前記第1回路に接続されているときは、前記フレーム音響線信号に基づいて超音波画像を生成すると決定し、前記第2回路が前記第1回路から切り離されているときは、前記第1音響線信号に基づいて超音波画像を生成すると決定する
    請求項1に記載の超音波診断装置。
  3. 前記パラメータ算出部は、前記第1音響線信号に基づいて前記被検体内の動きを検出し、
    前記制御部は、被検体内の動きが大きいほど、前記超音波プローブの振動子の並ぶ方向における第2対象領域の幅が小さくなるように、前記第2整相加算部に指示する
    ことを特徴とする請求項1または2に記載の超音波診断装置。
  4. 前記パラメータ算出部は、前記第1音響線信号から生成したフレーム音響線信号と、他のフレームにおける第1音響線信号から生成したフレーム音響線信号とから、前記被検体内の動きを検出する
    ことを特徴とする請求項3に記載の超音波診断装置。
  5. 前記第1整相加算部は、整相加算に用いる被検体内の超音波速度について複数の推定値を保持しており、前記送信イベントごとに、同一の第1対象領域に対して、前記複数の推定値のそれぞれを用いて、複数の第1音響線信号を生成し、
    前記パラメータ算出部は、前記複数の第1音響線信号に基づいて、被検体内の超音波速度を推定し、
    前記制御部は、前記パラメータ算出部が推定した被検体内の超音波速度を用いて整相加算を行うよう、第2整相加算部に指示する
    ことを特徴とする請求項1から4のいずれか1項に記載の超音波診断装置。
  6. 前記パラメータ算出部は、第1対象領域が同一で超音波速度の推定値が異なる、複数の第1音響線信号のそれぞれに対して、信号の値の分散を算出し、
    算出した分散が最も大きい第1音響線信号に対応する推定値を、被検体内の超音波速度であると推定する
    ことを特徴とする請求項5に記載の超音波診断装置。
  7. 前記パラメータ算出部は、前記第1音響線信号における、観測点の深さと信号のS/N比との関係を算出して、観測点の深さに応じた増幅率を決定し、
    前記制御部は、前記パラメータ算出部が決定した、観測点の深さに応じた増幅率を用いて、第2音響線信号に重みづけを行うよう、第2整相加算部に指示する
    ことを特徴とする請求項1から6のいずれか1項に記載の超音波診断装置。
  8. 前記第1整相加算部は、観測点の深さに対する信号増幅率について複数のプロファイルを保持しており、前記送信イベントごとに、同一の第1対象領域に対して、前記複数のプロファイルのそれぞれを用いて、複数の第1音響線信号を生成し、
    前記パラメータ算出部は、前記第1音響線信号における、観測点の深さに対する信号のS/N比の分散が最も小さいプロファイルを、観測点の深さに応じた増幅率として決定する
    ことを特徴とする請求項7に記載の超音波診断装置。
  9. 前記第2整相加算部は、送信された超音波が各観測点に到達する送信時間を、前記観測点の深さが前記被検体内で超音波が集束するフォーカス深さ以上である場合には、送信された超音波が前記送信振動子列から前記フォーカス点に到達するまでの第1時間と、前記基準点から前記観測点に到達する第2時間との合計を送信時間として算出し、前記観測点の深さが前記被検体内で超音波が集束するフォーカス深さ未満である場合には、前記第1時間から前記第2時間を減算した結果を送信時間として算出する
    ことを特徴とする請求項1から8のいずれか1項に記載の超音波診断装置。
  10. 前記画像生成部は、前記第1音響線信号からフレーム音響線信号を生成する第1音響線信号合成部を備え、
    前記第1対象領域は、前記フォーカス点を通過し、全域が前記対象領域内に含まれる直線状の領域であり、
    前記第1整相加算部は、前記送信振動子列に含まれる振動子のそれぞれに対応する受信信号列に対して、観測点と前記フォーカス点との距離と、前記フォーカス点と振動子との距離とをそれぞれ被検体内の超音波速度で除算した値に基づく遅延処理を含む整相加算を行って前記第1音響線信号を生成し、
    前記画像生成部は、前記第1音響線信号から超音波画像を生成する場合、前記第1音響線信号合成部が生成したフレーム音響線信号から超音波画像を生成する
    ことを特徴とする請求項1から9のいずれか1項に記載の超音波診断装置。
  11. 前記第1整相加算部は、送信された超音波が各観測点に到達する送信時間を、前記観測点の深さが前記被検体内で超音波が集束するフォーカス深さ以上である場合には、送信された超音波が前記送信振動子列から前記フォーカス点に到達するまでの第1時間と、前記基準点から前記観測点に到達する第2時間との合計を送信時間として算出し、前記観測点の深さが前記被検体内で超音波が集束するフォーカス深さ未満である場合には、前記第1時間から前記第2時間を減算した結果を送信時間として算出する
    ことを特徴とする請求項10に記載の超音波診断装置。
  12. 前記第1整相加算部は、各観測点からの反射波が各振動子に到達する受信時間を、前記観測点から前記観測点に最近接する振動子に到達するまでの時間を、前記観測点に最近接する振動子に対応する受信時間として算出し、超音波が前記フォーカス点から振動子に到達するまでの時間と前記フォーカス点から前記観測点に最近接する振動子に到達するまでの時間との差を、前記観測点に最近接する振動子に対応する受信時間に加算することで、当該振動子に対応する受信時間として算出する
    ことを特徴とする請求項10または11に記載の超音波診断装置。
  13. 前記第1音響線信号合成部は、前記送信イベントごとに、前記対象領域内の各観測点について、前記フォーカス点との距離が当該観測点と同一であり前記直線上に存在する観測点の前記第1音響線信号を割り当てて当該観測点の音響線信号とすることで、サブフレーム音響線信号を生成し、生成した複数のサブフレーム音響線信号を合成してフレーム音響線信号を生成する
    ことを特徴とする請求項10から12のいずれか1項に記載の超音波診断装置。
  14. 前記第1対象領域は、前記フォーカス点またはその近傍を通過し、前記超音波プローブの振動子の並ぶ方向に直交する1以上の直線からなり、
    前記第1整相加算部は、観測点の深さを被検体内の超音波速度で除算した値を送信時間、観測点から振動子までの距離を被検体内の超音波速度で除算した値を受信時間として整相加算を行い、前記第1音響線信号を生成し、
    前記画像生成部は、前記第1音響線信号から超音波画像を生成する場合、前記第1整相加算部が生成した複数の第1音響線信号から超音波画像を生成する
    ことを特徴とする請求項1から9のいずれか1項に記載の超音波診断装置。
  15. 複数の振動子を備えた超音波プローブを用いて被検体に集束型の超音波ビームを送信する送信イベントを複数回繰り返すとともに、各送信イベントに同期して被検体から反射超音波を受波して受信信号列を生成し、前記受信信号列に基づいて超音波画像を生成する超音波診断装置であって、
    超音波ビームの集束位置を規定するフォーカス点を送信イベントごとに変更しながら、送信イベントごとに、前記超音波プローブに列設された複数の振動子から送信振動子列を選択して当該送信振動子列から被検体内の対象領域に超音波ビームを送信させる送信部と、
    各送信イベントに同期して、前記超音波プローブが前記対象領域から受波した反射超音波に基づいて、振動子各々に対する受信信号列を生成する受信部と、
    前記送信イベントごとに、少なくとも前記対象領域の一部を含む第1対象領域内の複数の観測点について、前記受信信号列を整相加算し第1音響線信号を生成する第1整相加算部と、
    前記第1音響線信号に基づいて、サブフレーム音響線信号を生成するためのパラメータを算出するパラメータ算出部と、
    前記第1音響線信号から超音波画像を生成する超音波画像生成部と、
    制御部と
    を含む超音波信号処理回路を備え、
    前記対象領域の全部または一部である第2対象領域内の複数の観測点について、前記パラメータに基づき、前記受信信号列を整相加算しサブフレーム音響線信号を生成する第2整相加算部と、
    前記サブフレーム音響線信号を合成してフレーム音響線信号を生成する合成部と
    を含む演算回路が前記超音波信号処理回路に対して着脱可能に構成され、
    前記制御部は、前記演算回路が前記超音波信号処理回路に接続されているとき、前記超音波画像生成部に、前記第1音響線信号に替えて前記フレーム音響線信号から超音波画像を生成させる
    ことを特徴とする超音波診断装置。
JP2017004351A 2017-01-13 2017-01-13 超音波診断装置、および、その制御方法 Expired - Fee Related JP6724797B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2017004351A JP6724797B2 (ja) 2017-01-13 2017-01-13 超音波診断装置、および、その制御方法
EP17205915.6A EP3349034A3 (en) 2017-01-13 2017-12-07 Ultrasound diagnostic apparatus and control method therefor
US15/861,936 US20180199918A1 (en) 2017-01-13 2018-01-04 Ultrasonic diagnostic apparatus and control method therefor
CN201810018290.0A CN108294778B (zh) 2017-01-13 2018-01-09 超声波诊断装置及其控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017004351A JP6724797B2 (ja) 2017-01-13 2017-01-13 超音波診断装置、および、その制御方法

Publications (2)

Publication Number Publication Date
JP2018110784A true JP2018110784A (ja) 2018-07-19
JP6724797B2 JP6724797B2 (ja) 2020-07-15

Family

ID=60627554

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017004351A Expired - Fee Related JP6724797B2 (ja) 2017-01-13 2017-01-13 超音波診断装置、および、その制御方法

Country Status (4)

Country Link
US (1) US20180199918A1 (ja)
EP (1) EP3349034A3 (ja)
JP (1) JP6724797B2 (ja)
CN (1) CN108294778B (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020039604A (ja) * 2018-09-11 2020-03-19 コニカミノルタ株式会社 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法
US11484292B2 (en) 2018-08-10 2022-11-01 Konica Minolta, Inc. Ultrasound signal processing device that uses synthetic aperture method and delay and sum method

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10705210B2 (en) * 2017-05-31 2020-07-07 B-K Medical Aps Three-dimensional (3-D) imaging with a row-column addressed (RCA) transducer array using synthetic aperture sequential beamforming (SASB)
CN110833432B (zh) * 2018-08-15 2023-04-07 深南电路股份有限公司 超声波模拟前端装置及超声波成像设备
US20200214675A1 (en) * 2018-12-10 2020-07-09 Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of National Defence 3d ultrasound system and method
FR3099580B1 (fr) 2019-08-02 2024-11-01 Centre Nat Rech Scient Procédé et système de caractérisation ultrasonore non invasive d’un milieu hétérogène
JP7488710B2 (ja) * 2020-07-21 2024-05-22 富士フイルムヘルスケア株式会社 超音波撮像装置
WO2022049857A1 (ja) * 2020-09-03 2022-03-10 コニカミノルタ株式会社 超音波式検査装置、支持体の検査方法、及び、支持体の検査プログラム

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11221214A (ja) * 1997-11-10 1999-08-17 Medison Co Ltd 超音波映像化システムの超音波信号集束方法及び装置
JP2013121494A (ja) * 2011-11-07 2013-06-20 Toshiba Corp 超音波診断装置
JP2014030715A (ja) * 2012-07-13 2014-02-20 Fujifilm Corp 超音波検査装置
JP2015192712A (ja) * 2014-03-31 2015-11-05 富士フイルム株式会社 音響波処理装置、音響波処理装置の信号処理方法およびプログラム
JP2016127876A (ja) * 2015-01-09 2016-07-14 コニカミノルタ株式会社 超音波信号処理装置、及び超音波診断装置
JP2016154749A (ja) * 2015-02-25 2016-09-01 コニカミノルタ株式会社 超音波診断装置、及び超音波診断装置の制御方法
JP6059782B1 (ja) * 2015-10-01 2017-01-11 株式会社日立製作所 超音波診断装置及び遅延データ生成方法

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3827527B2 (ja) * 2000-12-28 2006-09-27 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波撮影装置
CN103142251B (zh) 2005-04-14 2015-04-29 维拉声学公司 利用面向像素处理的超声成像系统
US9332958B2 (en) * 2008-08-25 2016-05-10 Hitachi Medical Corporation Ultrasonic diagnostic apparatus and method of displaying ultrasonic image
US8968205B2 (en) * 2011-02-10 2015-03-03 Siemens Medical Solutions Usa, Inc. Sub-aperture control in high intensity focused ultrasound
US10201327B2 (en) * 2011-07-01 2019-02-12 Cliff A. Gronseth Organic specimen feature identification in ultrasound data
CN102908163B (zh) * 2011-08-02 2016-01-20 深圳迈瑞生物医疗电子股份有限公司 一种帧相关方法及其装置、以及超声成像系统
JP5972569B2 (ja) * 2011-12-26 2016-08-17 東芝メディカルシステムズ株式会社 超音波診断装置、超音波画像処置装置、医用画像診断装置及び超音波画像処理プログラム
JP6459451B2 (ja) * 2013-12-12 2019-01-30 コニカミノルタ株式会社 超音波診断装置、超音波画像処理方法、及びコンピュータ読み取り可能な非一時的な記録媒体
EP3117774B1 (en) * 2014-03-14 2019-01-30 Alpinion Medical Systems Co., Ltd. Software-based ultrasound imaging system
JP6165089B2 (ja) * 2014-03-25 2017-07-19 富士フイルム株式会社 音響波処理装置、音響波処理装置の信号処理方法およびプログラム
JP6380194B2 (ja) * 2015-03-27 2018-08-29 コニカミノルタ株式会社 超音波信号処理装置、及び、超音波診断装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11221214A (ja) * 1997-11-10 1999-08-17 Medison Co Ltd 超音波映像化システムの超音波信号集束方法及び装置
JP2013121494A (ja) * 2011-11-07 2013-06-20 Toshiba Corp 超音波診断装置
JP2014030715A (ja) * 2012-07-13 2014-02-20 Fujifilm Corp 超音波検査装置
JP2015192712A (ja) * 2014-03-31 2015-11-05 富士フイルム株式会社 音響波処理装置、音響波処理装置の信号処理方法およびプログラム
JP2016127876A (ja) * 2015-01-09 2016-07-14 コニカミノルタ株式会社 超音波信号処理装置、及び超音波診断装置
JP2016154749A (ja) * 2015-02-25 2016-09-01 コニカミノルタ株式会社 超音波診断装置、及び超音波診断装置の制御方法
JP6059782B1 (ja) * 2015-10-01 2017-01-11 株式会社日立製作所 超音波診断装置及び遅延データ生成方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11484292B2 (en) 2018-08-10 2022-11-01 Konica Minolta, Inc. Ultrasound signal processing device that uses synthetic aperture method and delay and sum method
JP2020039604A (ja) * 2018-09-11 2020-03-19 コニカミノルタ株式会社 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法
JP7147399B2 (ja) 2018-09-11 2022-10-05 コニカミノルタ株式会社 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法

Also Published As

Publication number Publication date
CN108294778B (zh) 2021-06-08
EP3349034A2 (en) 2018-07-18
JP6724797B2 (ja) 2020-07-15
CN108294778A (zh) 2018-07-20
EP3349034A3 (en) 2018-10-10
US20180199918A1 (en) 2018-07-19

Similar Documents

Publication Publication Date Title
JP6724797B2 (ja) 超音波診断装置、および、その制御方法
JP6733530B2 (ja) 超音波信号処理装置、超音波信号処理方法、及び、超音波診断装置
CN105997137B (zh) 超声波信号处理装置、以及超声波诊断装置
JP6406019B2 (ja) 超音波信号処理装置、及び超音波診断装置
JP6579231B2 (ja) 超音波画像生成方法、及び超音波診断装置
CN102970935B (zh) 超声波诊断装置以及超声波诊断装置控制方法
JP6747108B2 (ja) 超音波信号処理装置、超音波信号処理方法、及び、超音波診断装置
US20180055486A1 (en) Ultrasound signal processing device, ultrasound diagnostic device, and ultrasound signal processing method
JP6708101B2 (ja) 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法
US10743842B2 (en) Ultrasound signal processor, ultrasound signal processing method, and ultrasound diagnostic device
JP7052385B2 (ja) 超音波信号処理装置、超音波信号処理方法、および、超音波診断装置
US12042330B2 (en) Ultrasound signal processing device, ultrasound diagnostic device, and ultrasound signal processing method
JP7124505B2 (ja) 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法
JP6562122B2 (ja) 超音波診断装置、及び、超音波画像生成方法
JP7020052B2 (ja) 超音波信号処理装置、超音波診断装置、超音波信号処理方法、および、超音波画像表示方法
JP2023104734A (ja) 超音波診断装置、及び画像処理装置
CN118806329A (zh) 用于非衍射声学波束的发射的系统和方法
US20190150888A1 (en) Ultrasound signal processing device and ultrasound signal processing method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190327

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200228

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200310

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200421

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200526

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200608

R150 Certificate of patent or registration of utility model

Ref document number: 6724797

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees