JP6825474B2 - 超音波診断装置、及び超音波信号処理方法。 - Google Patents

超音波診断装置、及び超音波信号処理方法。 Download PDF

Info

Publication number
JP6825474B2
JP6825474B2 JP2017085699A JP2017085699A JP6825474B2 JP 6825474 B2 JP6825474 B2 JP 6825474B2 JP 2017085699 A JP2017085699 A JP 2017085699A JP 2017085699 A JP2017085699 A JP 2017085699A JP 6825474 B2 JP6825474 B2 JP 6825474B2
Authority
JP
Japan
Prior art keywords
acoustic line
line signal
unit
data
ultrasonic
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.)
Active
Application number
JP2017085699A
Other languages
English (en)
Other versions
JP2018183297A (ja
Inventor
津島 峰生
峰生 津島
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 JP2017085699A priority Critical patent/JP6825474B2/ja
Priority to US15/942,634 priority patent/US10702246B2/en
Publication of JP2018183297A publication Critical patent/JP2018183297A/ja
Application granted granted Critical
Publication of JP6825474B2 publication Critical patent/JP6825474B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/54Control of the diagnostic device
    • 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/52047Techniques for image enhancement involving transmitter or receiver for elimination of side lobes or of grating lobes; for increasing resolving power
    • 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/52087Details related to the ultrasound signal acquisition, e.g. scan sequences using synchronization techniques
    • 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/4477Constructional features of the ultrasonic, sonic or infrasonic diagnostic device using several separate ultrasound transducers or 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/5284Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving retrospective matching to a physiological signal

Landscapes

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

Description

本開示は、超音波信号処理方法、及び、それを用いた超音波診断装置に関し、特に、超音波信号処理方法における受信ビームフォーミング処理方法に関する。
超音波診断装置は、超音波プローブより被検体内に超音波を送信し、被検体組織内の音響インピーダンスの差異から生じる超音波の反射波を受信し、得られた電気信号に基づいて被検体の内部組織の形状を示す超音波断層画像を生成して表示する装置である。
従来の超音波診断装置では、受信した反射波に基づく受信ビームフォーミング方法として、整相加算法と呼ばれる方法が使用されている(例えば、非特許文献1)。この方法では、複数の振動子から被検体のある深さでフォーカスを結ぶよう超音波ビームを送信して超音波ビームの中心軸上にある音響線信号を生成する。そのため、中心軸上以外に照射した超音波の利用効率が悪く、送信フォーカス点近傍以外では得られる音響線信号の空間分解能や信号S/N比が低いという課題を有していた。
これに対し、合成開口法(Synthetic Aperture Method)を用いて、送信フォーカス点近傍以外の領域においても空間分解能の高い高画質な画像を得る受信ビームフォーミング方法が提案されている(例えば、非特許文献2)。合成開口法によれば、超音波送信波の伝播経路と、その伝播経路による反射波の振動子への到達時間の両方を加味した遅延制御を行うことで、1回の超音波送信から送信フォーカス点近傍以外を含めた超音波主照射領域全体に対して音響線信号を生成することができる。さらに、合成開口法では、複数回の超音波送信から得た同一観測点に対する複数の音響線信号を重ね合わせて、空間分解能及びS/N比を向上することができる。
伊東正安、望月剛共著「超音波診断装置」コロナ社出版、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
ところが、合成開口法においては、1回の超音波送信で音響線信号が生成される対象領域の面積が大きくなると領域に内在する観測点の数が増加するため、整相加算における演算量、整相加算後の音響線信号のデータを格納するために必要なメモリ容量、音響線信号のデータを転送するために必要なデータ伝送能力が増加する。そのため、整相加算の演算処理を高速に行うため演算処理能力の高いハードウェアと、音響線信号を格納するメモリ容量や転送するためのデータ伝送能力の拡大が必要となり、超音波診断装置のコストが増加するという課題があった。
本開示は、上記課題に鑑みてなされたものであり、空間分解能及びS/N比の低下を抑制しながら整相加算の演算量、生成される音響線信号のデータ量を削減でき、超音波診断装置の受信ビームフォーマにおける必要な内部メモリ容量やデータ伝送能力を縮減できる超音波信号処理方法、及び、それを用いた超音波診断装置を提供することを目的とする。
本開示の一態様に係る超音波診断装置は、複数の振動子を備えた超音波プローブを用いて被検体に超音波ビームを送信し、被検体から得られた反射波に基づいて音響線信号サブフレームデータを生成する超音波診断装置であって、被検体内における超音波ビームの集束点を規定するフォーカス点を設定し、前記複数の振動子から選択される複数の送信振動子の列から前記集束点に集束する超音波ビームを送信させる送信部と、前記超音波プローブが前記被検体から受波した反射波に基づいて、前記超音波プローブの振動子各々に対する受信信号列を生成する受信部と、前記超音波ビームの主照射領域内の位置に対応する複数の領域内観測点のうち前記フォーカス点を通過する対象線群上に存在する複数の線上観測点について、前記主照射領域内から得られた反射波に基づく前記受信信号列を整相加算して音響線信号サブフレームデータを生成する整相加算部と、を備え、前記受信部は、前記送信振動子の列を分割した複数の部分振動子列の各々について、その部分振動子列に含まれる振動子各々が被検体から受波した反射波に基づき受信信号列を生成する複数のパート受信部と、前記整相加算部は、前記複数の部分振動子列に対応する受信信号列を整相加算することにより前記複数の線上観測点について音響線信号を生成して、前記複数の部分振動子列の各々に対応する音響線信号部分サブフレームデータを生成する複数のパート整相加算部と、前記複数の部分振動子列の各々に対応する前記音響線信号部分サブフレームデータから、前記複数の線上観測点に対する音響線信号の列を抽出し、前記対象線群に含まれる対象線間の間隔を詰めて配列することにより音響線信号部分サブフレーム折り畳みデータを生成する複数のパート折り畳み部と、前記複数の部分振動子列の各々に対応する前記音響線信号部分サブフレーム折り畳みデータを前記配列された観測点の位置を基準に加算して音響線信号サブフレーム折り畳みデータを生成する主加算部と、前記音響線信号フレーム折り畳みデータにおける音響線信号を、前記主照射領域における前記線上観測点の位置に再配列して前記音響線信号サブフレームデータを生成する再配列部とを有することを特徴とする。
本開示の一態様に係る超音波信号処理方法、及び、それを用いた超音波診断装置によれば、空間分解能及びS/N比の低下を抑制しながら整相加算の演算量、及び、生成される音響線信号のデータ量を削減することができる。その結果、超音波診断装置の受信ビームフォーマにおける必要な内部メモリ容量やデータ伝送能力を抑制できる。
実施の形態に係る超音波診断装置100の構成を示す機能ブロック図である。 送信ビームフォーマ部103による超音波送信波の伝播経路を示す模式図である。 受信ビームフォーマ部104の構成を示す機能ブロック図である。 受信部40、整相加算部41の構成を示す機能ブロック図である。 (a)(b)は、受信ビームフォーマ部104における、音響線信号の生成方法の概要を示す模式図である。 (a)から(c)は、部分振動子列101a_lの分割数nを2としたときの、パート受信部401、及び整相加算部41のパート受信部401、パート整相加算部411における分割処理の概要を示す模式図、(c)は、音響線信号サブフレームデータds_sf生成の概要を示す模式図である。 実施の形態に係るパート整相加算部411の構成を示す機能ブロック図である。 対象線群Bxを示す模式図である。 受信開口設定部4113により設定された受信開口Rxと送信開口Txとの関係を示す模式図である。 (a)(b)は、送信開口Txから観測点Pijを経由して受信振動子Rkに到達する超音波の伝播経路を示す模式図である。 実施の形態に係るパート折り畳み部412の構成を示す機能ブロック図である。 (a)から(d)は、パート折り畳み部412の動作を示す説明図である。 再配列部414の動作を示す説明図である。 実施の形態に係る合成部42の構成を示す機能ブロック図である。 加算処理部421における合成音響線信号を合成する処理を示す模式図である。 合成音響線信号における最大重畳数と、増幅処理部422における増幅処理の概要を示す模式図である。 実施の形態に係るビームフォーミング処理動作を示すフローチャートである。 図17におけるステップS224の詳細を示すフローチャートである。 観測点Pijについての音響線信号生成動作を説明するための模式図である。 図17におけるステップS240の詳細を示すフローチャートである。 図17におけるステップS270の詳細を示すフローチャートである。 (a)〜(d)は、実施例及び比較例1の受信ビームフォーミングにより得た超音波画像である。 実施例に係る超音波診断装置100Aにおける受信ビームフォーマ部104の模式実装図である。 別の実施例に係る超音波診断装置100Bにおける受信ビームフォーマ部104の模式実装図である。 比較例に係る受信ビームフォーマ部104における観測点Pijについての音響線信号生成動作を説明するための模式図である。
≪発明を実施するための形態に至った経緯≫
発明者は、合成開口法を用いる超音波診断装置において、音響線信号の空間分解能及びS/N比の低下を抑止しながら演算量を削減するために各種の検討を行った。
一般に、集束型の送信ビームフォーミングでは、被検体のある深さで超音波ビームがフォーカスを結ぶよう波面を集束させる。そのため、1度の超音波の送信によって、超音波送信に用いられる複数の振動子から、超音波主照射領域に主として超音波が照射される。送信フォーカス点が1点である場合には、超音波主照射領域は、送信振動子列を底辺とし、底辺の両端のそれぞれから送信フォーカス点を通る2つの直線で囲まれる砂時計形状の領域となり、波面は、送信フォーカス点を中心とした円弧状となる。すなわち、超音波主照射領域は、フォーカス点のある深さではフォーカス点またはその近傍に集束し、それ以外の深さでは、フォーカス深さまでの距離が遠いほど列方向の幅が広がる形状となる。
合成開口法では、1回の送信イベントにおいて、超音波主照射領域の全域に対して観測点を設定できるため、超音波主照射領域全域を、対象領域とするのが好ましい。
しかしながら、対象領域に含まれる観測点の数は対象領域の面積に比例するため、必然的に、整相加算の演算量、整相加算後の音響線信号を格納するために必要なメモリ容量、音響線信号を転送するために必要なデータ伝送能力は、対象領域の面積に比例することとなる。したがって、対象領域の面積の増大は、そのまま、超音波診断装置が必要とするメモリ容量の増大に直結する。時間解像度の低下やユーザビリティの低下を抑制するためには、整相加算の演算を高速に行えるような処理能力の高いプロセッサ、例えば高性能のGPUなどが必要となり、超音波診断装置のコスト増加を招くこととなる。
演算量を削減するには、対象領域に含まれる観測点の数を削減することが考えられる。発明者は、観測点の数を削減するための対象領域として、フォーカス点またはその近傍を通過する複数の対象線からなる対象線群を設定するという着想を得た。このようにすることで、深さ方向には観測点の数も密度を維持したうえで、対象線と垂直方向には観測点の密度を低下させることで距離分解能の縮小しない構成を採るとともに、演算対象となる観測点の数を削減することで、整相加算の演算量を削減することができる。
さらに、整相加算処理の分散処理化によりハードウエアコストを低減することが考えられる。その場合、整相加算処理を複数の演算器に分散化して1演算器あたりの演算量を削減することが必要である。例えば、超音波プローブに列設された複数の振動子の列を複数の部分振動子列に分割し、それぞれの部分振動子列毎に得られる受信信号に基づいて前段の整相加算処理を行い、それぞれの整相加算処理から生成された音響線信号を後段で合成するという方法が考えられる。この構成では、部分振動子列毎に対応した前段の整相加算処理と後段の合成処理とのそれぞれにおいて、例えば、小規模のFPGA(Field Programmable Gate Array)を複数用いることでハードウェアを構成できるので、ハードウェアにおける演算器のコストを大幅に低減することができる。すなわち、1石のFPGA等を用いた集中処理により機能を実装した場合に、演算規模が大きいことより生じるコスト増大を回避できる。すなわち、大規模演算を可能とするFPGA等のハードウエアは非常に高価であり、同じ演算量を複数の集積回路に分割して処理した方が、コスト低減を低減できる。
しかしながら、処理を複数の集積回路に分割しで実現する場合には、集積回路間でデータ転送を行う必要があり、総演算量は賄えてもデータ転送が非常に多く、安価なハードウエアでは実現できない場合がある。具体的には、前段の整相加算により得た音響線信号を一旦格納するために必要なメモリ容量や、前段の処理から後段の処理へ音響線信号を転送するために必要なデータ伝送能力がコスト低減における律速となる。これらを低減するためには、上記した演算対象となる観測点の数を削減する方法だけでは不十分であり、整相加算後の音響線信号の外形的なデータ量を縮減することが必要となる。
そこで、発明者は、特別な演算負荷を伴うことなく音響線信号の外形的なデータ量を縮減できる方法について鋭意検討を行い、以下に示す実施の形態に想到するに至った。
以下、実施の形態に係る超音波画像処理方法及びそれを用いた超音波診断装置について図面を用いて詳細に説明する。
≪実施の形態≫
1.全体構成
以下、実施の形態に係る超音波診断装置100について、図面を参照しながら説明する。
図1は、実施の形態に係る超音波診断装置100の構成を示す機能ブロック図である。図1に示すように、超音波診断システム1000は、被検体に向けて超音波を送信しその反射波の受信する複数の振動子101aを有する超音波プローブ101(以後、「プローブ101」とする)、プローブ101に超音波の送受信を行わせプローブ101からの出力信号に基づき超音波画像を生成する超音波診断装置100、超音波画像を画面上に表示する表示部106を有する。プローブ101、表示部106は、それぞれ、超音波診断装置100に各々接続可能に構成されている。図1は超音波診断装置100に、プローブ101、表示部106が接続された状態を示している。なお、プローブ101と、表示部106とは、超音波診断装置100の内部にあってもよい。
2.超音波診断装置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、受信ビームフォーマ部104、超音波画像生成部105、制御部108は、それぞれ、例えば、FPGA(Field Programmable Gate Array)、ASIC(Application Specific Integrated Circuit)などのハードウェア回路により実現される。あるいは、プロセッサなどのプログラマブルデバイスとソフトウェアにより実現される構成であってもよい。プロセッサとしてはCPU(Central Processing Unit)やGPGPUを用いることができ、GPUを用いる構成はGPGPU(General-Purpose computing on Graphics Processing Unit)と呼ばれる。これらの構成要素は一個の回路部品とすることができるし、複数の回路部品の集合体にすることもできる。また、複数の構成要素を組合せて一個の回路部品とすることができるし、複数の回路部品の集合体にすることもできる。
データ格納部107は、コンピュータ読み取り可能な記録媒体であり、例えば、フレキシブルディスク、ハードディスク、MO、DVD、DVD−RAM、BD、半導体メモリ等を用いることができる。また、データ格納部107は、超音波診断装置100に外部から接続された記憶装置であってもよい。
3.超音波診断装置100の主要部の構成
実施の形態1に係る超音波診断装置100は、プローブ101の各振動子101aから超音波送信を行わせる送信ビームフォーマ部103と、プローブ101での超音波反射波の受信から得た電気信号を演算して超音波画像を生成するための音響線信号を生成する受信ビームフォーマ部104に特徴を有する。そのため、本明細書では、主に、送信ビームフォーマ部103及び受信ビームフォーマ部104について、その構成及び機能を説明する。なお、送信ビームフォーマ部103及び受信ビームフォーマ部104以外の構成については、公知の超音波診断装置に使われるものと同じ構成を適用可能であり、公知の超音波診断装置のビームフォーマ部に本実施の形態に係るビームフォーマ部を置き換えて使用することが可能である。
以下、送信ビームフォーマ部103と、受信ビームフォーマ部104の構成について説明する。
3.1 送信ビームフォーマ部103
送信ビームフォーマ部103は、マルチプレクサ部102を介してプローブ101と接続され、プローブ101から超音波の送信を行うためにプローブ101に存する複数の振動子101aから選択され、複数の振動子101aの全てもしくは一部に当たる送信振動子の列(送信開口Tx)に含まれる複数の振動子の各々に対する高電圧印加のタイミングを制御する。送信ビームフォーマ部103は送信部1031から構成される。
送信部1031は、制御部108からの送信制御信号に基づき、プローブ101に存する複数の振動子101a中、送信開口Txに含まれる各振動子に超音波ビームを送信させるためのパルス状の送信信号を供給する送信処理を行う。具体的には、送信部1031は、例えば、クロック発生回路、パルス発生回路、遅延回路を備えている。クロック発生回路は、超音波ビームの送信タイミングを決定するクロック信号を発生させる回路である。パルス発生回路は、各振動子を駆動するパルス信号を発生させるための回路である。遅延回路は、超音波ビームの送信タイミングを振動子毎に遅延時間を設定し、遅延時間だけ超音波ビームの送信を遅延させて超音波ビームのフォーカシングを行うための回路である。
送信部1031は、超音波送信ごとに送信開口Txを列方向に順次移動させながら超音波送信を繰り返し、プローブ101に存する全ての振動子101aから超音波送信を行う。すなわち、本実施の形態では、超音波送信ごとに、送信開口Txが振動子1つ分ずつ移動していく。送信開口Txに含まれる振動子の位置を示す情報は制御部108を介してデータ格納部107に出力される。例えば、プローブ101に存する振動子101a全数を192としたとき、送信開口Txを構成する振動子列の数として、例えば20〜100を選択してもよく、超音波送信毎に振動子1つ分だけ移動させる構成としてもよい。以後、送信部1031により同一の送信開口Txから行われる超音波送信を「送信イベント」と称呼する。
図2は、送信ビームフォーマ部103による超音波送信波の伝播経路を示す模式図である。ある送信イベントにおいて、超音波送信に寄与するアレイ状に配列された振動子101aの列(送信振動子列)を送信開口Txとして図示している。また、送信開口Txの列長を送信開口長と呼ぶ。
送信ビームフォーマ部103において、送信開口Txの中心に位置する振動子ほど送信タイミングを遅らせるように各振動子の送信タイミングを制御する。これにより、送信開口Tx内の振動子列から送信された超音波送信波は、被検体のある深度(Focal depth)において、波面がある一点、すなわち送信フォーカス点F(Focal point)で、フォーカスがあう(集束する)状態となる。
必ずしも超音波ビームが1点でフォーカスを結ぶとは限らず、例えば、振動子1.5個分から数個分程度にフォーカスした領域に集束するだけの場合もあるが、この場合、超音波主照射領域はフォーカス深さまでは列方向の幅が狭まり、フォーカス深さでフォーカス領域の列方向の幅となり、フォーカス深さより深い領域では再び列方向が広がる形状となる。なお、この場合においては、フォーカス深さにおける、フォーカス領域の中心点を便宜上「フォーカス点」と規定する。
送信フォーカス点Fの深さ(Focal depth)(以下、「フォーカス深さ」とする)は、任意に設定することができる。送信フォーカス点Fで合焦した波面は、再び拡散し、送信開口Txを底とし送信フォーカス点Fを節とする交差する2つの直線で区切られた砂時計型の空間内を超音波送信波が伝播する。すなわち、送信開口Txで放射された超音波は、次第にその空間上での幅(図中の横軸方向)を小さくし、送信フォーカス点Fでその幅を最小化し、それよりも深部(図中では下部)に進行するにしたがって、再び、その幅を大きくしながら拡散し、伝播することとなる。この砂時計型の領域が超音波主照射領域Axである。
ここで、「超音波主照射領域」Axとは、領域内のすべての点において、送信振動子列を構成する各振動子から送波される超音波の位相が揃っている領域を指す。なお、上述したように、超音波主照射領域Axは、1点の送信フォーカス点Fの近傍に集束するように超音波送信波を送信してもよい。
3.2 受信ビームフォーマ部104の構成
受信ビームフォーマ部104は、プローブ101で受信した超音波の反射波に基づき、複数の振動子101aで得られた電気信号から音響線信号を生成する。なお、「音響線信号」とは、ある観測点に対する、整相加算処理がされた後の信号である。整相加算処理については後述する。図3は、受信ビームフォーマ部104の構成を示す機能ブロック図である。図4は、受信部40、整相加算部41の構成を示す機能ブロック図である。図3に示すように、受信ビームフォーマ部104は、受信部40、整相加算部41、合成部42を備える。
以下、受信ビームフォーマ部104を構成する各部の構成の概要について説明する。
3.2.1 受信部40
受信部40は、マルチプレクサ部102を介してプローブ101と接続され、送信イベントに同期してプローブ101での超音波反射波の受信から得た電気信号を増幅した後AD変換した受信信号(RF信号)を生成する回路である。送信イベントの順に時系列に受信信号を生成しデータ格納部107に出力し、データ格納部107に受信信号を保存する。
ここで、受信信号(RF信号)とは、各振動子にて受信された反射超音波から変換された電気信号をA/D変換したデジタル信号であり、各振動子にて受信された超音波の送信方向(被検体の深さ方向)に連なった信号の列を形成している。
送信イベントでは、上述のとおり、送信部1031は、プローブ101に存する複数の振動子101a中、送信開口Txに含まれる複数の振動子の各々に超音波ビームを送信させる。これに対し、受信部40は、送信イベントに同期してプローブ101に存する複数の振動子101aの一部又は全部にあたる振動子の各々が得た反射超音波に基づいて、各振動子に対する受信信号の列を生成する。ここで、反射超音波を受波する振動子を「受波振動子」と称呼する。受波振動子の数は、送信開口Txに含まれる振動子の数よりも多いことが好ましい。また、受波振動子の数はプローブ101に存する振動子101aの全数としてもよい。
受信部40は、図4に示すように、プローブ101に存在する複数の振動子101aから選択されたm個(mは1より大きい自然数)の送信振動子の列101aTxをn分割(nはm以下であり1より大きい自然数)した複数の部分振動子列101a_l(l=1〜n)の各々について、その部分振動子列101a_lに含まれる振動子各々が被検体から受波した反射波に基づき受信信号列を生成する複数(n)のパート受信部401_lに分かれている。各部分振動子列101a_lには、m/n個の振動子が含まれる。例えば、送信振動子の列101aTxの個数mを128、nを2、部分振動子列101a_lに含まれる振動子の個数を64としてもよい。なお、本明細書において、パート受信部401_lは、順番lを区別しない場合、「パート受信部401」と表記する。
送信部1031は、送信イベントに同期して送信開口Txを列方向に順次移動させながら超音波送信を繰り返し、プローブ101に存する複数の振動子101a全体から超音波送信を行う。受信部40は、送信イベントに同期して各受波振動子に対する受信信号の列を生成し、生成された受信信号はデータ格納部107に保存される。
3.2.2 整相加算部41
整相加算部41は、超音波ビームの主照射領域Ax内の位置に対応する複数の領域内観測点Pijのうちフォーカス点Fを通過する対象線群Bx上に存在する複数の線上観測点PBxijについて、超音波主照射領域Ax内から得られた反射波に基づく受信信号列を整相加算して音響線信号サブフレームデータdssfを生成する回路である。なお、「対象線群」の詳細については、後述する。
図4に示すように、整相加算部41は、複数(n)のパート整相加算部411_l(l=1〜n)、複数(n)のパート折り畳み部412_l、主加算部413、再配列部414を備える。
3.2.2.1 パート整相加算部411
(1)概 要
パート整相加算部411は、図4に示すように、複数(n)のパート受信部401_lからそれぞれ入力した複数の部分振動子列101a_lに対応する受信信号列を整相加算することにより複数の線上観測点PBxijについて音響線信号を生成する。複数(n)のパート整相加算部411_lは、それぞれ、複数の部分振動子列101a_lの各々に対応して、複数の領域内観測点Pijのうち対象線群Bx上の全ての観測点PBxijに対する音響線信号部分サブフレームデータds_lを生成する。なお、本明細書において、パート整相加算部411_lは、順番lを区別しない場合、「パート整相加算部411」と表記する。
図5(a)(b)は、受信ビームフォーマ部104における、音響線信号の生成方法の概要を示す模式図である。同図では、簡単のために、振動子列(101a)を2つの部分振動子列101a_1、a_2に分割し、整相加算部41が2つのパート整相加算部411_1、411_2より構成される場合を示した。
図5(a)(b)に示すように、パート整相加算部411_lでは、部分振動子列101a_lに対する整相加算処理として、超音波主照射領域Axに存在する複数の域内観測点Pijのうち、フォーカス点Fを通過する対象線群Bx上に存在する複数の線上観測点PBxijについて整相加算処理が行われる。このとき、パート整相加算部411_lにおける整相加算処理の対象となる受信信号参照範囲は、部分振動子列101a_lに含まれる受信振動子Rk(kは1より大きい自然数)の列(受信開口Rx)が被検体から受波した反射波に基づく受信信号列rfである。
パート整相加算部411における分割処理について、さらに、具体的に説明する。
i)受信部40における、受信信号列rfの生成
図6(a)から(c)は、部分振動子列101a_lの分割数nを2としたときの、パート受信部401、及び整相加算部41のパート整相加算部411における分割処理の概要を示す模式図である。(a)はパート受信部401_1、401_2における受受信号列rfの取得範囲、(b)はパート整相加算部411_1、(c)はパート整相加算部411_2による整相加算処理における、受受信号列rfの参照範囲と音響線信号部分サブフレームデータds_1、ds_2を生成すべき線上観測点PBxijの存在範囲を示したものである。
送信イベントでは、上述のとおり、送信部1031は、プローブ101に存する複数の振動子101a中、全てもしくは一部に当たる振動子の列を送信開口Txとして設定する。図6(a)では、送信開口Txが、部分振動子列101a_1、101a_2に跨って位置している状態を示している。送信部1031は、送信開口Txに含まれる複数の振動子の各々に送信フォーカス点Fに集束するよう超音波ビームを送信させる。超音波ビームは、送信フォーカス点Fで波面が合焦した後拡散し、超音波主照射領域Axは送信開口Txを上辺とする砂時計型となる。そして、超音波主照射領域Ax内の組織で反射した反射波ecは各振動子101aに到達し、パート受信部401_1、401_2は、図6(a)に示すように、それぞれ、部分振動子列101a_1、101a_2に含まれる全ての振動子が受波した反射波ecに基づき受信信号列をrf生成する。
ii)パート整相加算部411_1における、音響線信号部分サブフレームデータデータds_1の生

次に、図6(b)に示すように、パート整相加算部411_1は、パート受信部401_1が出力した受信信号列rfに基づき、超音波主照射領域Axに含まれる対象線群Bx上に位置する線上観測点PBxij全てについて、それぞれに受信開口Rxを設定して整相加算処理を行う。受信開口Rxについては後に詳述するが、例えば、受信開口Rxを構成する振動子列の列中心が線上観測点PBxijに最も近接する振動子と合致するよう選択する構成としてもよい。
例えば、部分振動子列101a_1下方の超音波主照射領域Ax内に位置する線上観測点PBxij(A)_1については受信開口Rx(A)が選択される。この場合には、線上観測点PBxij(A)_1については、パート受信部401_1により生成された受信信号列rfのうち、受信開口Rx(A)内の振動子により受波された受信信号列rfに基づき整相加算が行われ、線上観測点PBxij(A)_1に対する音響線信号部分データds_1が生成される。図6(b)に示すように、受信開口Rx(A)内の振動子により受波されパート受信部401_1により算出された受信信号列rfは、受信開口Rx(A)の全部の振動子に基づくものであるため、線上観測点PBxij(A)_1に対する音響線信号部分データds_1は反射率に応じて最大の信号強度を得ることができる。
他方、部分振動子列101a_2下方の超音波主照射領域Ax内に位置する線上観測点PBxij(B)_1については受信開口Rx(B)が選択される。この場合には、線上観測点PBxij(B)_1については、パート受信部401_1により生成された受信信号列rfのうち、受信開口Rx(B)内の振動子により受波された受信信号列rfにのみ基づき整相加算が行われ、線上観測点PBxij(B)_1に対する音響線信号部分データds_1が生成される。図6(b)に示すように、受信開口Rx(B)内の振動子により受波されパート受信部401_1により算出された受信信号列rfは、受信開口Rx(B)の一部の振動子に基づくものであるため、線上観測点PBxij(B)_1に対する音響線信号部分データds_1は小さな信号強度しか得られない。
このように、パート整相加算部411_1は、パート受信部401_1が出力した受信信号列rfに基づき、対象線群Bx上に位置する線上観測点PBxij全てについて、音響線信号部分サブフレームデータデータds_1を生成する。
iii)パート整相加算部411_2における、音響線信号部分サブフレームデータデータds_2の生

図6(c)に示すように、パート整相加算部411_2は、パート受信部401_2が出力した受信信号列rfに基づき、対象線群Bx上に位置する線上観測点PBxij全てについて、それぞれに受信開口Rxを設定して整相加算処理を行う。
例えば、部分振動子列101a_2下方に位置する線上観測点PBxij(B)_2については受信開口Rx(B)が選択される。この場合には、線上観測点PBxij(B)_2については、パート受信部401_2により生成された受信信号列rfのうち、受信開口Rx(B)内の振動子により受波された受信信号列rfに基づき整相加算が行われ、線上観測点PBxij(B)_2に対する音響線信号部分データds_2が生成される。図6(c)に示すように、受信開口Rx(B)内の振動子により受波されパート受信部401_2により算出された受信信号列rfは、受信開口Rx(A)の大部分の振動子に基づく受信信号列rfであるため、線上観測点PBxij(B)_2に対する音響線信号部分データds_2は十分な信号強度を得ることができる。
他方、部分振動子列101a_1下方に位置する線上観測点PBxij(A)_2については受信開口Rx(A)が選択される。この場合には、線上観測点PBxij(A)_2については、パート受信部401_2により生成された受信信号列rfのうち、受信開口Rx(A)内の振動子により受波された受信信号列rfにのみ基づき整相加算が行われ、線上観測点PBxij(A)_2に対する音響線信号部分データds_2が生成される。図6(c)に示すように、受信開口Rx(A)内の振動子により受波されパート受信部401_2により算出された受信信号列rfは存在しないため、線上観測点PBxij(A)_2に対する音響線信号部分データds_2の信号強度はゼロとなる。
このように、パート整相加算部411_2は、パート受信部401_2が出力した受信信号列rfに基づき、対象線群Bx上に位置する線上観測点PBxij全てについて、音響線信号部分サブフレームデータデータds_2を生成する。
パート整相加算部411_1、411_2での整相加算方法は、受信アポダイゼーションによる重み係数をap、受信信号の強度をrf、複数の振動子101aを示す変数をk(kは1より大きい自然数)としたとき、それぞれ、次式により示される。
(2) パート整相加算部411の各構成
図7は、パート整相加算部411の構成を示す機能ブロック図である。図7に示すように、パート整相加算部411は、対象線群設定部4112、受信開口設定部4113、送信時間算出部4114、受信時間算出部4115、遅延量算出部4116、遅延処理部4117、重み算出部4118、及び加算部4119を備える。
以下、パート整相加算部411を構成する各部の構成について説明する。
i)対象線群設定部4112
対象線群設定部4112は、被検体内において音響線信号サブフレームデータの生成を行う対象線群Bxを設定する。「対象線群」とは、送信イベントに同期して被検体内において音響線信号サブフレームデータの生成が行われるべき信号上の領域であり、対象線群Bx上の観測点Pijについて音響線信号が生成される。対象線群Bxは、音響線信号の生成が行われる観測対象点の集合として、1回の送信イベントに同期して計算の便宜上設定される。
なお、上述のとおり、パート整相加算部411_lは、図4に示すように、複数(n)のパート受信部401_lからそれぞれ入力した複数の部分振動子列101a_lに対応して分割されている。そのため、本実施の形態では、パート整相加算部411_lは、それぞれ、部分振動子列101a_lが受波した信号列から出力される受信信号列rfを入力信号として取得する。しかしながら、パート整相加算部411_lにおける整相加算処理では、超音波主照射領域Ax全体に位置する対象線群Bx上の全ての観測点PBxijに対する音響線信号部分サブフレームデータds_lを生成する構成を採る。
なお、本実施の形とは異なる態様として、パート整相加算部411_lにおける整相加算処理として、超音波主照射領域Ax全体に位置する対象線群Bx上の観測点PBxijのうち部分振動子列101a_lに対応する範囲内に位置する観測点について音響線信号部分サブフレームデータds_lを生成する構成を採ることも可能である。
ここで、「音響線信号サブフレームデータ」とは、1回の送信イベントから生成される対象線群Bx上に存在する全ての観測点Pijに対する音響線信号の集合である。なお、「サブフレーム」とは、1回の送信イベントで得られ、対象線群Bx上に存在する全ての観測点Pijに対応するまとまった信号を形成する単位をさす。取得時間の異なる複数のサブフレームを合成したものがフレームとなる。
対象線群設定部4112は、送信イベントに同期して、送信ビームフォーマ部103から取得する送信開口Txの位置を示す情報に基づき対象線群Bxを設定する。
図8は、対象線群Bxを示す模式図である。図8に示すように、対象線群Bxは、超音波主照射領域Ax内に存在し、複数の対象線BL1〜BL7(以後、対象線の区別を要しない場合は、「対象線BL」と総称する)からなる。各対象線は、フォーカス点Fまたはその近傍を通過する直線である。なお、対象線BL1およびBL7は、それぞれ超音波主照射領域Axの外郭線に該当し、対象線BL4は送信開口Tx中心軸上に存在する。なお、超音波主照射領域Axの外郭線は便宜上、送信開口Txの一端とフォーカス点Fとを通過する直線、および、送信開口Txの他端とフォーカス点Fとを通過する直線、の2つであるものとする。また、各対象線BL1〜BL7において、隣接する対象線がなす角は概ね等しい。すなわち、フォーカス点Fを中心とする円弧上に存在する観測点は等間隔に並ぶ。また、例えば、iを対象線の識別番号を示す変数、jを対象線に沿った位置を示す変数としたとき、Pijとフォーカス点Fとの距離が所定距離未満である場合を除き、同一の対象線BL2上の隣接する観測点PijとPi(j+1)との間の距離djは、隣接する対象線BL2、BL3上のそれぞれの観測点PijとP(i+1)jとの間の距離diよりも小さい。なお、距離djは、少なくとも距離diの2倍であり、好ましくは4倍以上、より好ましくは8倍以上である。このようにすることで、観測点は深さ方向に高密度、振動子の並ぶ方向(フォーカス点Fを中心とした円周方向)には低密度となるように、超音波主照射領域Axのほぼ全域に均等に配置される。ここで、所定距離とは、対象線上の点と隣接する対象線上の点との距離が、対象線に沿った観測点の間隔より小さくなる範囲を示す。例えば、隣接する対象線のなす角をθとした場合、所定距離dpは以下の式を満たす。
di=2・dp・sin(θ/2)
なお、対象線群Bxの形状は、上述の場合に限られず、例えば、各対象線BL1〜BL7において、送信振動子列と接する位置の距離が等間隔である、としてもよい。また、対象線群Bxは7本の対象線からなるが、これは一例であり、対象線の数は、3以上の範囲内であれば、任意に設定してよい。また、本例では、簡単のために奇数の場合を示して説明を行っているが、言うまでもなく対象線の数は奇数に限定されない。
また、各観測点Pijは対象線群Bx上に存在するとしたが、例えば、観測点の一部または全部が、対象線に近接する位置に設定されてもよい。例えば、i、jを直交座標を示す変数としたとき、観測点Pijは、振動子の並ぶ方向(x方向)と深さ方向(z方向)とで示され、各振動子の中心を格子点とする直交座標における格子点上に存在する、とすることができる。このようにすることで、各観測点Pijに対してx座標が同一となる振動子が必ず存在するため、音響線信号の品質を向上させることができるからである。この場合、対象線は必ずしもz方向と平行でないため、観測点を設けようとする位置で、対象線が格子点を通過しないことがある。このような場合、観測点は対象線の直上ではなく、対象線に近接する格子点上に設ける。例えば、対象線上に設けるべき観測点の座標について、所定の桁数において四捨五入などの丸め処理を行うことで、観測点の実際の座標を規定する。具体的には、以下のように観測点を定める。振動子の並ぶ方向(x方向)について、プローブの振動子の数を192、プローブの振動子列の一端の振動子の位置をx=0、他端の振動子の位置をx=191とする。深さ方向について、振動子列の位置をz=0とし、プローブの振動子1つ分だけ深い点をz=1とする。フォーカス点Fの座標を(64,1000)として、座標(31,0)を通る対象線を設定すると、対象線は以下の数式で示される。
z=(1000/33)・(x−31)
この時、例えば、z=1500の深さに観測点を設けようとすると、(80.5,1500)となる。この場合、観測点を(81,1500)としてもよい。このようにすることで、整相加算においてx=81の振動子を基準とすることができ、音響線信号の品質を向上させることができる。なお、観測点を実際に設ける位置は上述した場合に限られず、対象線群Bx上に設けるべき観測点について、その座標の値を丸めて計算上好ましい近接点を実際の観測点とする限りにおいて、任意に行ってよい。
設定された対象線群Bxは送信時間算出部4114、受信時間算出部4115、遅延処理部4117に出力される。
ii)受信開口設定部4113
受信開口設定部4113は、制御部108からの制御信号と、送信ビームフォーマ部103からの送信開口Txの位置を示す情報とに基づき、プローブ101に存する複数の振動子の一部に当たる振動子列(受信振動子列)を受信振動子として選択して受信開口Rxを設定する回路である。ここで、「受信開口Rx」とは、線上観測点PBxijの音響線信号を生成するために整相加算の対象となる受信信号列を受波した振動子の列である。本実施の形態では、受信開口設定部4113は、列中心が観測点Pijに最も空間的に近接する振動子Xkと合致するよう受信開口Rx振動子列を選択する。図9は、受信開口設定部4113により設定された受信開口Rxと送信開口Txとの関係を示す模式図である。図9に示すように、受信開口Rx振動子列の列中心が、観測点Pijに最も空間的に近接する振動子Xkと合致するように受信開口Rx振動子列が選択される。そのため、受信開口Rxの位置は、観測点Pijの位置によって定まり、送信イベントに同期して変動する送信開口Txの位置に基づいては変化しない。すなわち、異なる送信イベントであっても、同一位置にある観測点Pijについての音響線信号を生成する処理においては、同一の受信開口Rx内の受信振動子Rkによって取得された受信信号に基づき整相加算が行われる。
また、超音波主照射領域全体からの反射波を受信するために、受信開口Rxに含まれる振動子の数は、対応する送信イベントにおける送信開口Txに含まれる振動子の数以上に設定することが好ましい。受信開口Rxを構成する振動子列の数は、例えば32、64、96、128、192等としてもよい。
受信開口Rxの設定は、送信イベントに対応して、少なくとも送信イベントと同じ回数だけ行われる。また、受信開口Rxの設定は、送信イベントに同期して漸次行われる構成であってもよく、あるいは、全ての送信イベントが終了した後に、各送信イベントに対応した受信開口Rxの設定が送信イベントの回数分まとめて行われる構成であってもよい。
選択された受信開口Rxの位置を示す情報は制御部108を介してデータ格納部107に出力される。
データ格納部107は、受信開口Rxの位置を示す情報と受信振動子に対応する受信信号とを、送信時間算出部4114、受信時間算出部4115、遅延処理部4117、重み算出部4118に出力する。
iii)送信時間算出部4114
送信時間算出部4114は、送信された超音波が被検体中の観測点Pijに到達する送信時間を算出する回路である。送信イベントに対応して、データ格納部107から取得した、送信開口Txに含まれる振動子の位置を示す情報と、対象線群設定部4112から取得した超音波主照射領域Axを含む対象線群Bxの位置を示す情報とに基づき対象線群Bx上に存在する任意の観測点Pijについて、送信された超音波が被検体中の観測点Pijに到達する送信時間を算出する。
図10(a)(b)は、送信開口Txから放射され対象線群Bx上の任意の位置にある観測点Pijにおいて反射され受信開口Rx内に位置する受信振動子Rkに到達する超音波の伝播経路を説明するための模式図であり、図10(a)は観測点Pijの深さがフォーカス深さ以上の場合、図10(b)は観測点Pijがフォーカス深さより浅い場合を示している。
送信開口Txから放射された送信波は、経路401を通って送信フォーカス点Fにて波面が集束し、再び、拡散する。送信波が集束または拡散する途中で観測点Pijに到達し、観測点Pijで音響インピーダンスに変化があれば反射波を生成し、その反射波がプローブ101における受信開口Rx内の受信振動子Rkに戻っていく。送信フォーカス点Fは送信ビームフォーマ部103の設計値として規定されているので、送信フォーカス点Fと任意の観測点Pijとの間の経路402の長さは幾何学的に算出することができる。
送信時間の算出方法を、以下、さらに詳細に説明する。
まず、観測点Pijの深さがフォーカス深さ以上の場合について、図10(a)を用いて説明する。観測点Pijの深さがフォーカス深さ以上の場合は、送信開口Txから放射された送信波が、経路401を通って送信フォーカス点Fに到達し、送信フォーカス点Fから経路402を通って観測点Pijに到達したものとして算出する。したがって、送信波が経路401を通過する時間と、経路402を通過する時間を合算した値が、送信時間となる。具体的な算出方法としては、例えば、経路401の長さと経路402の長さとを加算した全経路長を、被検体内における超音波の伝搬速度で除算することで求められる。
一方、観測点Pijがフォーカス深さより浅い場合について、図10(b)を用いて説明する。観測点Pijがフォーカス深さより浅い場合は、送信開口Txから放射された送信波が、経路401を通って送信フォーカス点Fに到達する時刻と、経路404を通って観測点Pijに到達した後、観測点Pijから経路402を通って送信フォーカス点Fに到達する時刻とが同一であるものとして算出する。つまり、送信波が経路401を通過する時間から、経路402を通過する時間を差し引いた値が、送信時間となる。具体的な算出方法としては、例えば、経路401の長さから経路402の長さを減算した経路長差を、被検体内における超音波の伝搬速度で除算することで求められる。
なお、観測点Pijがフォーカス深さである場合の送信時間は、観測点Pijがフォーカス深さより深い場合と同じ算出方法、すなわち、送信波が経路401を通過する時間と、経路402を通過する時間を合算する計算方法を用いるとした。しかしながら、観測点Pijがフォーカス深さより浅い場合と同じ算出方法、すなわち、送信波が経路401を通過する時間から、経路402を通過する時間を差し引く計算方法を用いるとしてもよい。経路402の長さが0となるため、いずれで算出しても経路401を通過する時間と一致するためである。
送信時間算出部4114は、1回の送信イベントに対し、対象線群Bx上の全ての観測点Pijについて、送信された超音波が被検体中の観測点Pijに到達する送信時間を算出して遅延量算出部4116に出力する。
iv)受信時間算出部4115
受信時間算出部4115は、観測点Pからの反射波が、受信開口Rxに含まれる受信振動子Rkの各々に到達する受信時間を算出する回路である。送信イベントに対応して、データ格納部107から取得した受信振動子Rkの位置を示す情報と、対象線群設定部4112から取得した対象線群Bxの位置を示す情報とに基づき対象線群Bx上に存在する任意の観測点Pijについて、送信された超音波が被検体中の観測点Pijで反射され受信開口Rxの各受信振動子Rkに到達する受信時間を算出する。
上述のとおり、観測点Pijに到達した送信波は、観測点Pijで音響インピーダンスに変化があれば反射波を生成し、その反射波がプローブ101における受信開口Rx内の各受信振動子Rkに戻っていく。受信開口Rx内の各受信振動子Rkの位置情報はデータ格納部107から取得されるので、任意の観測点Pijから各受信振動子Rkまでの経路403の長さは幾何学的に算出することができる。
受信時間算出部4115は、1回の送信イベントに対し、対象線群Bx上に存在する全ての観測点Pijについて、送信された超音波が観測点Pijで反射して各受信振動子Rkに到達する受信時間を算出して遅延量算出部4116に出力する。
v)遅延量算出部4116
遅延量算出部4116は、送信時間と受信時間とから受信開口Rx内の各受信振動子Riへの総伝播時間を算出し、当該総伝播時間に基づいて、各受信振動子Rkに対する受信信号の列に適用する遅延量を算出する回路である。遅延量算出部4116は、送信時間算出部4114から送信された超音波が観測点Pijに到達する送信時間と、観測点Pijで反射して各受信振動子Rkに到達する受信時間を取得する。そして、送信された超音波が各受信振動子Rkへ到達するまでの総伝播時間を算出し、各受信振動子Rkに対する総伝播時間の差異により、各受信振動子Rkに対する遅延量を算出する。遅延量算出部4116は、対象線群Bx上に存在する全ての観測点Pijについて、各受信振動子Riに対する受信信号の列に適用する遅延量を算出して遅延処理部4117に出力する。
vi)遅延処理部4117
遅延処理部4117は、受信開口Rx内の受信振動子Rkに対する受信信号の列から、各受信振動子Rkに対する遅延量に相当する受信信号を、観測点Pijからの反射超音波に基づく各受信振動子Rkに対応する受信信号として同定する回路である。
遅延処理部4117は、送信イベントに対応して、受信開口設定部4113から受信振動子Rkの位置を示す情報、データ格納部107から受信振動子Rkに対応する受信信号、対象線群設定部4112から取得した対象線群Bxの位置を示す情報、遅延量算出部4116から各受信振動子Rkに対する受信信号の列に適用する遅延量を入力として取得する。そして、各受信振動子Rkに対応する受信信号の列から、各受信振動子Rkに対する遅延量を差引いた時間に対応する受信信号を観測点Pijからの反射波に基づく受信信号として同定し、加算部4119に出力する。
vii)重み算出部4118
重み算出部4118は、受信開口Rxの列方向の中心に位置する振動子に対する重みが最大となるよう各受信振動子Rkに対する重み数列(受信アボダイゼーション)を算出する回路である。
図9に示すように、重み数列は受信開口Rx内の各振動子に対応する受信信号に適用される重み係数の数列である。重み数列は、送信フォーカス点Fを中心として対称な分布をなす。重み数列の分布の形状は、ハミング窓、ハニング窓、矩形窓などを用いることができ、分布の形状は特に限定されない。重み数列は、受信開口Rxの列方向の中心に位置する振動子に対する重みが最大となるように設定され、重みの分布の中心軸は、受信開口中心軸Rxoと一致し、分布は中心軸に対し対称な形状をなす。重み算出部4118は、受信開口設定部4113から出力される受信振動子Rkの位置を示す情報を入力として、各受信振動子Rkに対する重み数列を算出し加算部4119に出力する。
viii)加算部4119
加算部4119は、遅延処理部4117から出力される各受信振動子Rkに対応して同定された受信信号を入力として、それらを加算して、観測点Pijに対する整相加算された音響線信号を生成する回路である。あるいは、さらに、重み算出部4118から出力される各受信振動子Rkに対する重み数列を入力として、各受信振動子Rkに対応して同定された受信信号に、各受信振動子Rkに対する重みを乗じて加算して、観測点Pijに対する音響線信号を生成する構成としてもよい。遅延処理部4117において受信開口Rx内に位置する各受信振動子Rkが検出した受信信号の位相を整えて加算部4119にて加算処理をすることにより、観測点Pijからの反射波に基づいて各受信振動子Rkで受信した受信信号を重ね合わせてその信号S/N比を増加し、観測点Pijからの受信信号を抽出することができる。
1回の送信イベントとそれに伴うパート整相加算部411_lでの処理から、複数の部分振動子列101a_lの各々に対応して、超音波主照射領域Ax全体に位置する複数の領域内観測点Pijのうち対象線群Bx上の全ての観測点PBxijに対する音響線信号部分サブフレームデータds_lを生成することができる。
3.2.2.2 パート折り畳み部412
パート折り畳み部412_lは、パート整相加算部411_lが生成する複数の部分振動子列101a_lの各々に対応する音響線信号部分サブフレームデータds_lを入力として、複数の線上観測点PBxijに対する音響線信号ds(i,j)を、対象線群Bxを構成する対象線BL間の間隔を詰めて配列することにより音響線信号部分サブフレーム折り畳みデータdsc_lsを生成する回路である。
図11は、実施の形態に係るパート折り畳み部412の構成を示す機能ブロック図である。図11に示すように、パート折り畳み部412は、線上観測点音響線信号抽出部4121、配列決定部4122、対象線配列部4123を備える。
i)線上観測点音響線信号抽出部4121
線上観測点音響線信号抽出部4121は、音響線信号部分サブフレームデータds_lを入力として、線上観測点PBxijに対する音響線信号ds(i,j)を抽出する。
図12(a)から(d)は、パート折り畳み部412の動作を示す説明図である。パート整相加算部411_l(例として、l=1or2)では、図12(a)に示すように、部分振動子列101a_1、a_2からの超音波主照射領域Ax全体に位置する対象線群Bx(BL1〜BL7)上に線上観測点PBxijが設定され整相加算処理が行われる。このとき、パート整相加算部411_l又は411_2により生成される音響線信号部分サブフレームデータds_l(l=1or2)の態様は、図12(b)の模式図のようになる。音響線信号部分サブフレームデータds_lは、振動子列方向(x)と深さ(z)を2軸とする直交座量(i,j)(iが振動子列方向(x)に対応、jが深さ(z)に対応)上にマッピングされた、本例では13×29の音響線信号データの集合である。このうち、図12(b)の対象線群Bx近傍に位置する線上観測点PBxijには算出された121個の音響線信号ds(i,j)(図12(b)中の「●」)が内挿されており、線上観測点PBxij以外の観測点には256個のゼロを示す数固定値ds0(図12(b)中の「○」)が内挿されている。各音響線信号ds(i,j)には、音響線強度を示す情報の他、振動子列方向(x)、深さ(z)を示す座標情報(i,j)が含まれている。
線上観測点音響線信号抽出部4121は、例えば、音響線信号部分サブフレームデータds_lを、jの各位置でそれぞれiについてスキャンすることにより、音響線信号部分サブフレームデータds_l内の音響線信号ds(i,j)をサーチし、音響線信号ds(i,j)(図12(b)中の「●」)が検出されたときには、対象線群Bx(BL1〜BL7)中の対象線BLを識別する識別番号sを付与することにより、線上観測点PBxijに対する音響線信号ds(i,j)を抽出する。図12(a)(b)に示す例では、対象線BL1、2、3、4、5、6、7に、それぞれs=−2、−1、0、1、2の識別番号を付与する。そして、音響線信号部分サブフレームデータds_l上の各音響線信号ds(i,j)に対し、識別番号sを付加情報として追加する。したがって、各音響線信号dsc(s,j,i)は、音響線強度を示す情報の他、振動子列方向(x)、深さ(z)を示す座標情報(i,j)、対象線群Bx中の対象線BLの識別番号sを示す情報が含む信号となる。
ii)対象線配列部4123
対象線配列部4123は、抽出された線上観測点PBxijに対する音響線信号ds(i,j)を、対象線群Bx中の対象線BL間の間隔を詰めて配列する。具体的には、例えば、図12(c)に示すように、識別番号sを付加した音響線信号dsc(s,j,i)を、識別番号sとjを2軸とする直交座量(s,j)、sが振動子列方向(x)に対応、jが深さ(z)に対応)に配列することにより対象線BL間の間隔を詰めた音響線信号部分サブフレーム折り畳みデータdsc_lを生成する。あるいは、図12(d)に示すように、識別番号sの絶対値とjとを2軸とする直交座標に配列して音響線信号部分サブフレーム折り畳みデータdsc_lを生成してもよい。直交座量(s,j)上において、線上観測点PBxijに対する121個の音響線信号dsc(s,j,i)が対象線BL間の間隔を詰めて配列されており、音響線信号dsc(s,j,i)以外の部分には24個のゼロを示す数値ds0が内挿されている。
このように、音響線信号部分サブフレーム折り畳みデータdsc_lは、線上観測点PBxijに対する音響線信号ds(i,j)を対象線BL間の間隔を詰めて配列したものであるので、データ量を縮小することができ、必要な内部メモリ容量やデータ伝送能力を縮減できる。具体的には、図12(c)(d)に示す例では、音響線信号部分サブフレーム折り畳みデータdsc_lは、直交座量(s,j)上において、5×29の音響線信号データの集合となり、音響線信号部分サブフレームデータds_lに比べて、データ量が5/13(約38.5%)に縮小されていることがわかる。
iii)配列決定部4122
配列決定部4122は、制御部108からの情報に基づき、超音波診断装置100における超音波計測の条件を示す各種情報に基づき、対象線配列部4123における対象線群上の音響線信号の配列方法を規定する配列情報を決定して対象線配列部4123に出力し、対象線配列部4123は配列情報基づき音響線信号部分サブフレーム折り畳みデータdsc_lを生成する。このとき、配列決定部4122は、超音波計測の条件として、帯域(送受信周波数範囲)、チャネル数(振動子数)、チャネルピッチ(振動子間隔)、形状情報(コンベックス、リニア等)、送信開口(送信振動子列)、送信周波数、フォーカス位置、送信間隔、ステアリング角度、受信開口(受信振動子列)、最大深度、設定音速、モード種別(Bモード、Cモード、Dモード、Eモード、THI)、サンプリング周波数、解像度、PRF、送信間隔、密度、ハードウエアの演算能力、伝送能力、動作モードの種別(画質優先動作モード、フレームレート優先動作モード、バランス動作モード、低電力動作モード)から選択される1以上の情報を用いることができる。これにより、対象線配列部4123は、超音波計測条件に適合するような態様にて象線毎に音響線信号の列を配列することができ、超音波計測条件に応じて空間分解能及びS/N比の低下を抑制しながら整相加算の演算量及び、生成される音響線信号のデータ量を削減できる。
3.2.2.3 主加算部413
主加算部413は、複数の部分振動子列101a_lの各々に対応する音響線信号部分サブフレーム折り畳みデータdsc_lを、識別番号sの絶対値と深さ(z)とを2軸とする直交座量(s,j)上における観測点の位置を基準に加算して音響線信号サブフレーム折り畳みデータdscを生成してデータ格納部107に出力する。
3.2.2.4 再配列部414
再配列部414は、音響線信号サブフレーム折り畳みデータdsc_sfにおける音響線信号dsc(s,j,i)を、超音波主照射領域Axにおける線上観測点PBxijの位置に再配列して音響線信号サブフレームデータds_sfを生成する回路である。上述のとおり、音響線信号部分サブフレーム折り畳みデータdsc_l上の音響線信号dsc(s,j,i)には、音響線強度を示す情報の他、振動子列方向(x)、深さ(z)を示す座標情報(i,j)、対象線BLの識別番号sを示す情報が含まれている。したがって、音響線信号部分サブフレーム折り畳みデータdsc_lを加算して得た音響線信号サブフレーム折り畳みデータdsc上の音響線信号dsc(s,j,i)も、同様に音響線強度を示す情報の他、振動子列方向(x)、深さ(z)を示す座標情報(i,j)、対象線BLの識別番号sを示す情報を含む信号とすることができる。
図13は、再配列部414の動作を示す説明図である。したがって、図13に示すように、再配列部414は、音響線信号サブフレーム折り畳みデータdsc_sfを入力として、音響線信号dsc(s,j,i)を直交座標(i,j)上の線上観測点PBxijの位置に再配置し、線上観測点PBxij以外にはゼロを示す数固定値ds0内挿して音響線信号サブフレームデータds_sfを生成する。
1回の送信イベントとそれに伴う処理から、整相加算部41により、送信イベントに同期して対象線群Bx上に存在する全ての観測点Pijに対する音響線信号サブフレームデータds_sfが生成される。生成された音響線信号サブフレームデータds_sfは、データ格納部107に出力され保存される。
図6(d)は、音響線信号サブフレームデータds_sf生成の概要を示す模式図である。パート折り畳み部412_1、412_2、主加算部413、再配列部414における折り畳み・加算・展開処理(以後、「折り畳み・加算・展開処理」の前後のデータを比較したとき、パート整相加算部411_1で生成された音響線信号部分データds_1と、パート整相加算部411_2で生成された音響線信号部分データds_2とを線上観測点PBxijの位置を基準に算術的に加算して、検出波照射領域Ax内に存在する全ての線上観測点PBxijについて音響線信号サブフレームデータds_sfが生成される。
具体的には、図6(c)に示すように、パート整相加算部411_1で生成された線上観測点PBxij(A)_1及び(B)_1に対する音響線信号部分ムデータds_1と、パート整相加算部411_2で生成された線上観測点PBxij(A)_2及び(B)_2に対する音響線信号部分データds_2とを算術的に加算して線上観測点PBxij(A)及び(B)に対する音響線信号データds_1+2を生成する。これにより、線上観測点PBxij(A)及び(B)の両方において、十分な信号強度を有する音響線信号データds_1+2がそれぞれで得られる。そして、このような折り畳み・加算・展開処理を、全ての線上観測点PBxijについて行うことで、音響線信号サブフレームデータds_sfが生成される。加算処理は、次式により示される。
3.2.3 合成部42
送信イベントに同期して送信開口Txを列方向に順次移動させながら超音波送信を繰り返し、プローブ101に存する全ての振動子101aから超音波送信を行うことにより1フレームの合成された音響線信号である音響線信号フレームデータds_fを生成する。音響線信号フレームデータを構成する観測点ごとの合成された音響線信号を、以後、「合成音響線信号」と称呼する。
合成部42は、送信イベントに同期して生成される音響線信号サブフレームデータds_sfから音響線信号フレームデータds_fを合成する回路である。図14は、合成部42の構成を示す機能ブロック図である。図14に示すように、合成部42は、加算処理部421、増幅処理部422を備える。
以下、合成部42を構成する各部の構成について説明する。
i)加算処理部421
加算処理部421は、音響線信号フレームデータを合成するための一連の音響線信号サブフレームデータds_sfの生成が終了したのち、データ格納部107に保持されている複数の音響線信号サブフレームデータds_sfを読み出す。そして、各音響線信号サブフレームデータds_sfに含まれる音響線信号が取得された観測点Pijの位置を指標として複数の音響線信号サブフレームデータds_sfを加算することにより、各観測点に対する合成音響線信号を生成して音響線信号フレームデータds_fを合成する。そのため、複数の音響線信号サブフレームデータds_sfに含まれる同一位置の観測点に対する音響線信号は加算されて合成音響線信号が生成される。
図15は、加算処理部421における合成音響線信号を合成する処理を示す模式図である。上述のとおり、送信イベントに同期して送信振動子列(送信開口Tx)に用いる振動子を振動子列方向に振動子1つ分だけ異ならせて超音波送信が順次行われる。そのため、異なる送信イベントに基づく対象線群Bxも送信イベントごとに同一方向に振動子1つ分だけ位置が異なる。複数の音響線信号サブフレームデータds_sfを、各音響線信号サブフレームデータds_sfに含まれる音響線信号が取得された観測点Pijの位置を指標として加算することにより、全ての対象線群Bxを網羅した音響線信号フレームデータが合成される。
また、位置の異なる複数の対象線群Bxにまたがって存在する観測点Pijについては、各音響線信号サブフレームデータds_sfにおける音響線信号の値が加算されるので、合成音響線信号は、跨りの程度に応じて大きな値を示す。以後、観測点Pijが異なる対象線群Bxに含まれる回数を「重畳数」、振動子列方向における重畳数の最大値を「最大重畳数」と称する。
図16(a)は、合成音響線信号における最大重畳数と、増幅処理部422における増幅処理の概要を示す模式図である。本実施の形態では、対象線群Bxは砂時計形状の領域内に存在する。そのため、図16(a)に示すように、重畳数及び最大重畳数は被検体の深さ方向において変化するので、合成音響線信号の値も同様に深さ方向において変化する。但し、本実施例の形態では対象線群の数が例えば11の場合、最大重畳数は結果的に高々11に制限される。
なお、各音響線信号サブフレームデータds_sfに含まれる音響線信号が取得された観測点Pijの位置を指標として加算する際に、観測点Pijの位置を指標として重みづけしながら加算してもよい。
合成された音響線信号フレームデータds_fは増幅処理部422に出力される。
ii)増幅処理部422
上述のとおり、合成音響線信号の値は被検体の深さ方向において変化する。これを補うために、増幅処理部422は、音響線信号フレームデータds_fに含まれる合成音響線信号の合成において、加算が行われた回数に応じて決定した増幅率を各合成音響線信号に乗じる増幅処理を行う。
図16(b)は、増幅処理部422における増幅処理の概要を示す模式図である。図16(b)に示すように、最大重畳数は被検体の深さ方向において変化するので、この変化を補うように、最大重畳数に応じて決定された被検体深さ方向において変化する増幅率が合成音響線信号に乗じられる。これにより、深さ方向における重畳数の変化に伴う合成音響線信号の変動要因は解消され、増幅処理後の合成音響線信号の値は深さ方向において均一化が図られる。
また、重畳数に応じて決定された振動子列方向において変化する増幅率を合成音響線信号に乗じる処理を行ってもよい。振動子列方向において重畳数が変化する場合に、その変動要因を解消し、振動子列方向において増幅処理後の合成音響線信号の値の均一化が図られる。
なお、生成した各観測点に対する合成音響線信号に増幅処理を施した信号を音響線信号フレームデータとしてもよい。
4.動 作
以上の構成からなる超音波診断装置100の動作について説明する。
図17は、受信ビームフォーマ部104のビームフォーミング処理動作を示すフローチャートである。
先ず、ステップS101において、送信部1031は、プローブ101に存する複数の振動子101a中送信開口Txに含まれる各振動子に超音波ビームを送信させるための送信信号を供給する送信処理(送信イベント)を行う。
次に、ステップS102において、受信部40は、プローブ101での超音波反射波の受信から得た電気信号に基づき受信信号を生成しデータ格納部107に出力し、データ格納部107に受信信号を保存する。プローブ101に存する全ての振動子101aから超音波送信が完了したか否かを判定する(ステップS103)。そして、完了していない場合(ステップS103において、No)にはステップS101に戻り、送信開口Txを列方向に振動子1つ分だけ移動させながら送信イベントを行い、完了している場合(ステップS103、Yes)にはステップS201に進む。
次に、ステップS210において、対象線群設定部4112は、送信イベントに同期して、送信開口Txの位置を示す情報に基づき対象線群Bxを設定する。1回目のループでは初回の送信イベントにおける送信開口Txから求められる対象線群Bxが設定される。
次に、部分振動子列101a_lの識別番号lを1に設定する(ステップS215)。
次に、観測点同期型ビームフォーミング処理(ステップS220(S221〜S228))に進む。ステップS220では、まず、観測点Pijの位置を示す座標ijを対象線群Bx上の最小値に初期化し(ステップS221、S222)、受信開口設定部4113は、列中心が観測点Pijに最も空間的に近接する振動子Xkと合致するよう受信開口Rx振動子列を選択する(ステップS223)。これにより、次のビームフォーミング処理(ステップS220)において、パート整相加算部411_lは、部分振動子列101a_lのうち受信開口Rx振動子列と重複する範囲に含まれる振動子が被検体から受波した反射波に基づく受信信号列rfを入力信号として整相加算処理を行う。
次に、ビームフォーミング処理(ステップS220)により、観測点Pijについて音響線信号を生成する(ステップS224)。
ここで、ステップS224における、観測点Pijについて音響線信号を生成する動作について説明する。図18は、図17におけるステップS224の詳細を示すフローチャートであり、受信ビームフォーマ部104における観測点Pijについての音響線信号生成動作を示す。図19は、受信ビームフォーマ部104における観測点Pijについての音響線信号生成動作を説明するための模式図である。
まず、ステップS2241において、送信時間算出部4114は、対象線群Bx上に存在する任意の観測点Pijについて、送信された超音波が被検体中の観測点Pijに到達する送信時間を算出する。送信時間は、(1)観測点Pijの深さがフォーカス深さ以上の場合には、幾何学的に定まる受信開口Rx内の受信振動子Rkから送信フォーカス点Fを経由して観測点Pijに至る経路(401+402)の長さを超音波の音速csで除することにより、(2)観測点Pijがフォーカス深さより浅い場合には、幾何学的に定まる、受信開口Rx内の受信振動子Rkから送信フォーカス点Fに至る経路と観測点Pijからフォーカス点に至る経路との差分(401−402)の長さを、超音波の音速csで除することにより算出できる。
次に受信開口Rxから求められる受信開口Rx内の受信振動子Rkの位置を示す座標kを受信開口Rx内の最小値に初期化し(ステップS2242)、受信時間算出部4115は、送信された超音波が被検体中の観測点Pijで反射され受信開口Rxの受信振動子Rkに到達する受信時間を算出する(ステップS2243)。受信時間は、幾何学的に定まる観測点Pijから受信振動子Rkまでの経路403の長さを超音波の音速csで除することにより算出できる。さらに、遅延量算出部4116は、送信時間と受信時間の合計から、送信開口Txから送信された超音波が観測点Pijで反射して受信振動子Rkに到達するまでの総伝播時間を算出し(ステップS2244)、受信開口Rx内の各受信振動子Rkに対する総伝播時間の差異により、各受信振動子Rkに対する遅延量を算出する(ステップS2245)。
受信開口Rx内に存在する全ての受信振動子Rkについて遅延量の算出を完了したか否かを判定し(ステップS2246)、完了していない場合には座標kをインクリメント(ステップS2247)して、更に受信振動子Rkについて遅延量の算出を行い(ステップS2243)、完了している場合にはステップS2248に進む。この段階では、受信開口Rx内に存在する全ての受信振動子Rkについて観測点Pijからの反射波到達の遅延量が算出されている。
ステップS2248において、遅延処理部4117は、受信開口Rx内の受信振動子Rkに対応する受信信号の列から、各受信振動子Rkに対する遅延量を差引いた時間に対応する受信信号を観測点Pijからの反射波に基づく受信信号として同定する。
次に、重み算出部4118は、受信開口Rxの列方向の中心に位置する振動子に対する重みが最大となるよう各受信振動子Rkに対する重み数列を算出する(ステップS2249)。加算部4119は、各受信振動子Rkに対応して同定された受信信号に、各受信振動子Rkに対する重みを乗じて加算して、観測点Pijに対する音響線信号を生成し(ステップS2250)、生成された観測点Pijについて音響線信号はデータ格納部107に出力され保存される(ステップS2251)。
次に、図17に戻り、座標i、jをインクリメントしてステップS224を繰り返すことにより、対象線群Bx上の座標ijに位置する全ての観測点Pij(図19中の「・」)について音響線信号が生成される。対象線群Bx上に存在する全ての観測点Pijについて音響線信号の生成を完了したか否かを判定し(ステップS225、S227)、完了していない場合には座標ijをインクリメント(ステップS226、S228)して、観測点Pijについて音響線信号を生成し(ステップS224)、完了した場合にはステップS228に進む。この段階では、1回の送信イベントに伴う対象線群Bx上に存在する全ての観測点Pijについての音響線信号部分サブフレームデータds_lが生成され、データ格納部107に出力され保存されている。
次に、パート折り畳み部412_lは、音響線信号部分サブフレーム折り畳みデータの生成処理を行う(ステップS250)。図20は、図17におけるステップS240の詳細を示すフローチャートである。
対象線群Bx上の線上観測点PBxijの位置を示すj、jを最小値に初期化し(ステップS2411、S2412)、線上観測点音響線信号抽出部4121は、音響線信号部分サブフレームデータds_l上の座標(j、j)に位置する信号ds(i,j)をデータ格納部107から読み込み(ステップS2413)、ds(i,j)がゼロを示す固定値ds0(0)であるか否かを判定する(ステップS2414)。固定値ds0(0)である場合(ステップS2414において、No)には、ステップS2418に進む。ds(i,j)が、固定値ds0(0)でない場合(ステップS2414において、Yes)には、ds(i,j)が対象線群Bx上の線上観測点PBxijにおける音響線信号ds(i,j)であると判定して、音響線信号ds(i,j)を、音響線信号ds(i,j)に対象線群Bx中の対象線BLを識別する識別番号sを付加した音響線信号dsc(s,j,i)に置き換え(ステップS2415)、データ格納部107に出力する(ステップS2416)とともに、識別番号sをインクリメントする(ステップS2417)。
次に、位置を示すj、jのうち、全てのiについて処理を完了したか否か(ステップS2418)、超音波主照射領域Ax内の全てのjについて処理を完了したか否か(ステップS2420)について判定し、完了していない場合はi、jをインクリメント(ステップS2419、S2421)して、ds(i,j)が固定値ds0(0)であるか否かを判定し(ステップS2414)、全てのi、jについて処理を完了している場合には、図17におけるステップ250に進む。
この段階では、部分振動子列101a_lに関するイタレーションについて、超音波主照射領域Axに含まれる対象線群Bx上の線上観測点PBxijにおける全ての音響線信号ds(i,j)が、音響線信号dsc(s,j,i)に置き換えて抽出されている。これにより、部分振動子列101a_lに関する、音響線信号部分サブフレーム折り畳みデータdsc_lが生成される。図12(c)(d)において、説明したように、音響線信号部分サブフレーム折り畳みデータdsc_lは、線上観測点PBxijに対する音響線信号ds(i,j)を対象線BL間の間隔を詰めて配列したものであるので、データ量を縮小することができ、必要な内部メモリ容量やデータ伝送能力を縮減できる。
次に、全ての部分振動子列101a_lについて、音響線信号部分サブフレームデータds_lの生成が完了したか否かを判定し(ステップS250)、完了していない場合にはlをインクリメント(ステップS251)してステップS221に戻り、完了している場合にはステップS260に進む。
この段階では、1回の送信イベントに伴う全ての部分振動子列101a_lに対して音響線信号部分サブフレーム折り畳みデータdsc_lが生成され、データ格納部107に出力され保存されている。
次に、主加算部413は、複数の部分振動子列101a_lの各々に対応する音響線信号部分サブフレーム折り畳みデータdsc_lを、データ格納部107から読み出し、識別番号sの絶対値と深さ(z)とを2軸とする直交座量(s,j)上における観測点の位置を基準に、それらを加算して音響線信号サブフレーム折り畳みデータdscを生成する(ステップS260)。
続けて、再配列部414は、音響線信号サブフレーム折り畳みデータdsc_sfにおける音響線信号dsc(s,j,i)を、超音波主照射領域Axにおける線上観測点PBxijの位置に再配列して音響線信号サブフレームデータds_sfを生成する(ステップS270)図21は、図17におけるステップS270の詳細を示すフローチャートである。
対象線群Bx上の線上観測点PBxijの位置を示すj、jを最小値に初期化し(ステップS271、S272)、再配列部414は、音響線信号サブフレーム折り畳みデータdsc_sf上の座標(s,j)に位置する音響線信号dsc(s,j,i)をデータ格納部107から読み込み、対象線BLを識別する識別番号sを付加した音響線信号dsc(s,j,i)を音響線信号ds(i,j)に置き換え(ステップS273)、データ格納部107に出力する(ステップS274)。このとき、再配列部414は、音響線信号dsc(s,j,i)を直交座標(i,j)上の線上観測点PBxijの位置に再配列するとともに、線上観測点PBxij以外にはゼロを示す数固定値ds0内挿する。
次に、識別番号s、位置を示すjのうち、全てのsについて処理を完了したか否か(ステップS275)、超音波主照射領域Ax内の全てのjについて処理を完了したか否か(ステップS277)について判定し、完了していない場合はs、jをインクリメント(ステップS276、S277)して、音響線信号dsc(s,j,i)を音響線信号ds(i,j)に置き換え(ステップS273)、全てのs、jについて処理を完了している場合には、図17におけるステップ290に進む。
この段階では、再配列部414は、部分振動子列101a_lに関するイタレーションについて、音響線信号サブフレーム折り畳みデータdsc_sfを入力として、音響線信号サブフレームデータds_sfを生成している。
次に、図17に戻り、全ての送信イベントについて、音響線信号サブフレームデータds_sfの生成が終了したか否かを判定し(ステップS290)、終了していない場合には、ステップS210に戻り、観測点Pijの位置を示す座標ijを、次の送信イベントでの送信開口Txから求められる対象線群Bx上の最小値に初期化し(ステップS221、S222)、受信開口Rxを設定する(ステップS223)し、終了している場合にはステップS301に進む。
次に、ステップS301において、加算処理部421は、データ格納部107に保持されている複数の音響線信号サブフレームデータds_sfを読み出し、観測点Pijの位置を指標として複数の音響線信号サブフレームデータds_sfを加算して各観測点Pijに対する合成音響線信号を生成して音響線信号フレームデータds_fを合成する。次に、増幅処理部422は、音響線信号フレームデータds_fに含まれる各合成音響線信号の加算回数に応じて決定された増幅率を各合成音響線信号に乗じ(ステップS302)、増幅された音響線信号フレームデータds_fを、超音波画像生成部105及びデータ格納部107に出力し(ステップS303)処理を終了する。
5.画質評価について
実施の形態に係る受信ビームフォーマ部104を実施例に相当する例と比較例を用いて、生成されるBモード画像の画質評価を行った。
(1)供試サンプル
図22(a)〜(d)は、実施例及び比較例1の受信ビームフォーミングにより得た超音波画像である。サンプルの具体的な仕様には、以下のとおりである。
i)図22(a)
部分振動子列101a_lの分割数nを2としたときの、パート整相加算部411_1により生成された音響線信号部分サブフレームデータds_lに基づくBモード画像
ii)図22(b)
パート整相加算部411_2により生成された音響線信号部分サブフレームデータds_2に基づくBモード画像。
iii)図22(c)
音響線信号部分サブフレームデータds_lと音響線信号部分サブフレームデータds_2とを算術加算して得た音響線信号サブフレームデータds_sfに基づくBモード画像である。
(2)評価結果
図22(a)では、左半分に画像が描画され、(b)では右半分に画像が描画され、(c)では、良好な品質の画像が描画されている。
6.受信ビームフォーマ部104による効果
本実施の形態に係る超音波診断装置100の実施例について、主に、受信ビームフォーマ部104について、その構成及び機能を説明する。なお、受信ビームフォーマ部104以外の構成については、公知の超音波診断装置に使われるものと同じ構成を適用可能である。
(実施例1)
図23は、本実施の形態の実施例に係る超音波診断装置100Aにおける受信ビームフォーマ部104の模式実装図である。部分振動子列101a_lの分割数nを2としたときの、受信ビームフォーマ部104を構成する複数の集積回路において、受信部40、整相加算部41、合成部42における各機能がどのように配置されて実現されているかを示す。図23に示すように、超音波診断装置100Aでは、部分振動子列101a_1、101a_2は、128個の振動子からなる送信振動子101aの列をn分割(n=2)したものであり、それぞれ64個の振動子(64チャンネル)に接続された2個の集積回路501_1、501_2と、集積回路501_3とを備える。複数のパート受信部401_1、401_2、複数のパート整相加算部411_1、411_2、複数のパート折り畳み部412_1、412_2は、それぞれ各1づつ2個の集積回路501_1、501_2に含まれ、主加算部413、再配列部414は集積回路501_3に含まれる。そして、2個の集積回路501_1、501_2のそれぞれから、集積回路501_3に音響線信号部分サブフレーム折り畳みデータdsc_lが伝送される構成を採る。さらに、超音波診断装置100Aでは、合成部42も、集積回路501_3に含まれる構成としてもよい。
上述のとおり、音響線信号部分サブフレーム折り畳みデータdsc_lは、線上観測点PBxijに対する音響線信号ds(i,j)を対象線BL間の間隔を詰めて配列したものであるので、音響線信号部分サブフレームデータds_lに比べて、データ量を約20%から80%の範囲に縮減することができる(上述の例では、データ量が5/13(約38.5%)に縮減した)。集積回路501_1、501_2から集積回路501_3へのデータ転送に際し、パート折り畳み部412_1、412_2により音響線信号部分サブフレームデータds_lを音響線信号部分サブフレーム折り畳みデータdsc_lに変換してからデータ転送を行うことにより、転送に係るデータ量を縮減することができる。その結果、集積回路501_1、501_2、501_3の必要な内部メモリ容量や、集積回路501_1及び501_2〜集積回路501_3間の必要なデータ伝送能力の抑制することができる。
また、受信振動子101aの列をn分割して受信振動子101aに関する整相加算処理を複数の集積回路501_1、501_2、501_3に分散して行う。すなわち、整相加算処理を複数の演算器に分散化して1演算器あたりの演算量を削減することができる。これにより、部分振動子列毎に対応した前段の整相加算処理、折り畳み処理を行う集積回路501_1、501_2と、後段の主加算処理、再配列処理、合成処理とを行う集積回路501_3とのそれぞれにおいて、例えば、小規模のFPGA(Field Programmable Gate Array)を複数用いることでハードウェアを構成できるので、ハードウェアにおける演算器のコストを大幅に低減することができる。
合成部42は、送信イベントに同期して生成される音響線信号サブフレームデータds_sfを再配列部414から取得して音響線信号フレームデータds_fを合成し、データ格納部107に出力する回路である。超音波診断装置100Aでは、集積回路501_3がデータ格納部107を内部メモリーとして備える構成を採る。データ格納部107(内部メモリ)には、前送信イベントに同期して合成された音響線信号フレームデータds_fが保存されている。合成部42は、送信イベントごとに再配列部414から新たな音響線信号サブフレームデータds_sfを取得して新たに音響線信号フレームデータds_fを合成し、データ格納部107保存された音響線信号フレームデータds_fを新たに作成されたデータに更新する。同時に、合成部42は、前送信イベント時に作成された音響線信号フレームデータds_fと現信イベント時に作成された音響線信号フレームデータds_fとの差分に相当する音響線信号dsを超音波画像生成部105に出力する。送信イベント間の差分に相当する音響線信号dsとは1個の振動子101aの相当する幅の深さ方向に伸長した音響線信号dsの列であり、音響線信号フレームデータds_fのデータ量を送信イベント数で除したデータ量になる。本例では、送信イベント数を振動子101aの総数を128としたとき、音響線信号フレームデータds_fのデータ量に対し1/128のデータ量となる。
係る構成により、送信イベント毎に集積回路501_3から転送されるデータは音響線信号フレームデータds_fのうち送信イベント間の差分に相当する音響線信号dsのみとすることができ、集積回路501_3より後段へのデータ伝送に必要なデータ伝送能力を縮減することができる。
(実施例2)
図24は、本実施の形態の別の実施例に係る超音波診断装置100Bにおける受信ビームフォーマ部104の模式実装図である。ここでも、部分振動子列101a_lの分割数nは2としている。図24に示すように、超音波診断装置100Bでは、部分振動子列101a_1、101a_2は、受信振動子101aの列をn分割(n=2)したものであり、n個の集積回路511、511を備え、複数のパート受信部401_1、401_2、複数のパート整相加算部411_1、411_2、複数のパート折り畳み部412_1、412_2は、それぞれ各1づつn個の集積回路511、512に含まれ、主加算部413及び再配列部414は、n個の集積回路511、512の何れかに含まれ、主加算部413及び再配列部414を含む集積回路に対し、それ以外の集積回路のそれぞれから、音響線信号部分サブフレーム折り畳みデータdsc_lが伝送される構成を採る。さらに、超音波診断装置100Bでは、合成部42も、主加算部413及び再配列部414を含む集積回路に含まれる構成としてもよい。
超音波診断装置100Bでは、上記超音波診断装置100Aで示した効果に加え、さらに、以下の効果を得ることができる。超音波診断装置100Aでは、超音波プローブに列設された複数の振動子の列を複数の部分振動子列に分割し、それぞれの部分振動子列毎に得られる受信信号に基づいて前段で整相加算処理、折り畳み処理を行い、それぞれの整相加算処理から生成された音響線信号を後段で、主加算処理、再配列処理、合成処理を行う構成を採る。このうち、集積回路501_3にて行われる後段の主加算処理、再配列処理及び合成処理は、集積回路501_1、501_2にて行われる前段の整相加算処理、折り畳み処理に比べて演算量が少なく、集積回路501_1、501_2の何れかにて行う構成としてもよい。超音波診断装置100Bでは、後段の主加算処理、再配列処理及び合成処理を、前段の整相加算処理、折り畳み処理を行う集積回路511、512の何れかにて行わせる構成とした点に特徴がある。係る構成により、超音波診断装置100Bでは、集積回路501_3を削減することができ、ハードウェアにおける演算器のコストをさらに低減することができる。
なお、図24では、簡単のために、振動子列(101a)を2つの部分振動子列101a_1、a_2に分割し、整相加算部41が2つのパート整相加算部411_1、411_2より構成される場合を示したが、分割数nは2に限られないことは言うまでもない。
7.変形例
以上、実施の形態に係る超音波信号処理装置を説明したが、本発明は、その本質的な特徴的構成要素を除き、以上の実施の形態に何ら限定を受けるものではない。例えば、各実施の形態に対して当業者が思いつく各種変形を施して得られる形態や、本発明の趣旨を逸脱しない範囲で各実施の形態における構成要素及び機能を任意に組み合わせることで実現される形態も本発明に含まれる。以下では、そのような形態の一例として、変形例に係る超音波信号処理装置について説明する。
(1)実施の形態1に係る超音波診断装置100では、受信開口設定部4113は、列中心が線上観測点PBxijに最も空間的に近接する振動子と合致するよう受信開口Rxを選択する構成とした。しかしながら、受信開口Rxの構成は適宜変更することができる。
変形例1では、列中心が送信開口Tx振動子列の列中心と合致する受信開口Rx振動子列を選択する送信同期型受信開口設定部(以後、「Tx受信開口設定部」)を備えた点で実施の形態1と相違する。Tx受信開口設定部以外の構成については、実施の形態1に示した各要素と同じであり、同じ部分については説明を省略する。
図25は、Tx受信開口設定部により設定された受信開口Rxと送信開口Txとの関係を示す模式図である。変形例1では、受信開口Rx振動子列の列中心が送信開口Tx振動子列の列中心と合致するように受信開口Rx振動子列が選択される。受信開口Rxの中心軸の位置は、送信開口Txの中心軸の位置と同一であり、受信開口Rxは、送信フォーカス点Fを中心として対称な開口である。したがって、送信イベントにごとに列方向に移動する送信開口Txの位置変化に同期して、受信開口Rxの位置も移動する。
また、受信開口Rxの中心軸及び送信開口Txの中心軸上に位置する振動子に対する重みが最大となるよう受信開口Rxの各受信振動子Rkに対する重み数列(受信アボダイゼーション)は算出される。重み数列は、振動子Xkを中心として対称な分布をなす。重み数列の分布の形状は、ハミング窓、ハニング窓、矩形窓などを用いることができ、分布の形状は特に限定されない。
変形例1に係る超音波診断装置では、実施の形態1において示した効果のうち観測点同期型の受信開口に関する部分を除いた効果に変えて、以下の効果を奏する。すなわち、変形例1では、Tx受信開口設定部は送信イベントに対応して列中心が送信開口Txに含まれる振動子列の列中心と合致する振動子列を受信振動子として選択して受信開口Rxを設定する。そのため、受信開口Rxの中心軸の位置は、送信開口Txの中心軸の位置と同一であり、送信イベントにごとに列方向に移動する送信開口Txの位置変化に同期して、受信開口Rxの位置も変化(移動)する。よって、送信イベントに同期してそれぞれ異なる受信開口にて整相加算を行うことができ、複数の送信イベントにわたって受信時刻は異なるものの、結果としてより一層広い受信開口を用いた受信処理4の効果が得られ、広い観測領域で空間分解能を均一にすることができる。
(2)実施の形態では、対象線群Bxが、7本の対象線を等角度に配置した場合を例示したが、対象線群Bxの対象線の数は3以上であればよい。また、対象線の位置関係は、隣接する2本の対象線のなす角が一定である場合に限られず、例えば、同じ深さにおいて、対象線上の観測点の間隔が等間隔である、としてもよい。ここで、直線Lは振動子の並ぶ方向と平行な直線であり、直線Lと各対象線の交点の間隔は一定である。
なお、対象線群Bxは上述の例に限られず、送信フォーカス点Fを通過する複数の対象線からなり、隣接する2つの対象線上にあり、フォーカス点Fから等しい距離にある、または、同じ深さにある、2つの観測点の間の距離が、同一の対象線上にある2つの隣接する観測点の間の距離より大きい限り、任意の形状でよい。ただし、上述したように、対象線の向きは互いに類似しないほうが好ましいため、角度または振動子の並ぶ方向に均等であることが好ましい。また、2つの対象線がなす角の最大値は大きいことが好ましく、2つの対象線が超音波主照射領域Axの外郭線に一致していることが好ましい。さらに、1つの対象線が送信開口中心軸上にあることが好ましい。これは、最も信号強度の高い音響線信号を取得することができるからである。
(3)本実施の形態に係る超音波診断装置100は、図1で示した構成の超音波診断装置に限定されない。例えば、マルチプレクサ部102がなく、送信ビームフォーマ部103と受信ビームフォーマ部104とが直接、プローブ101の各振動子101aに接続されていてもよい。また、プローブ101に送信ビームフォーマ部103や受信ビームフォーマ部104、またその一部などが内蔵される構成であってもよい。これは、本実施の形態に係る超音波診断装置100に限られず、後に説明する他の実施の形態や変形例に係る超音波診断装置でも同様である。
(4)本開示を上記実施の形態に基づいて説明してきたが、本発明は、上記の実施の形態に限定されず、以下のような場合も本発明に含まれる。
例えば、本発明は、マイクロプロセッサとメモリを備えたコンピュータシステムであって、上記メモリは、上記コンピュータプログラムを記憶しており、上記マイクロプロセッサは、上記コンピュータプログラムにしたがって動作するとしてもよい。例えば、本発明の超音波診断装置の診断方法のコンピュータプログラムを有しており、このプログラムに従って動作する(又は接続された各部位に動作を指示する)コンピュータシステムであってもよい。
また、上記超音波診断装置の全部、もしくは一部、またビームフォーミング部の全部又は一部を、マイクロプロセッサ、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等のプロセッサがプログラムを実行することにより実現してもよい。上記超音波診断装置の診断方法や、ビームフォーミング方法を実施させるプログラムが記録された非一時的なコンピュータ読み取り可能な記録媒体であってもよい。プログラムや信号を記録媒体に記録して移送することにより、プログラムを独立した他のコンピュータシステムにより実施するとしてもよい、また、上記プログラムは、インターネット等の伝送媒体を介して流通させることができるのは言うまでもない。
また、上記実施形態に係る超音波診断装置の各構成要素は、CPU(Central Processing Unit)やGPU(Graphics Processing Unit)やプロセッサなどのプログラマブルデバイスとソフトウェアにより実現される構成であってもよい。後者の構成は、いわゆるGPGPU(General-Purpose computing on Graphics Processing Unit)である。これらの構成要素は一個の回路部品とすることができるし、複数の回路部品の集合体にすることもできる。また、複数の構成要素を組合せて一個の回路部品とすることができるし、複数の回路部品の集合体にすることもできる。
上記実施形態に係る超音波診断装置では、記憶装置であるデータ格納部を超音波診断装置内に含む構成としたが、記憶装置はこれに限定されず、半導体メモリ、ハードディスクドライブ、光ディスクドライブ、磁気記憶装置、等が、超音波診断装置に外部から接続される構成であってもよい。
また、ブロック図における機能ブロックの分割は一例であり、複数の機能ブロックを一つの機能ブロックとして実現したり、一つの機能ブロックを複数に分割したり、一部の機能を他の機能ブロックに移してもよい。また、類似する機能を有する複数の機能ブロックの機能を単一のハードウェア又はソフトウェアが並列又は時分割に処理してもよい。
また、上記のステップが実行される順序は、本発明を具体的に説明するために例示するためのものであり、上記以外の順序であってもよい。また、上記ステップの一部が、他のステップと同時(並列)に実行されてもよい。
また、超音波診断装置には、プローブ及び表示部が外部から接続される構成としたが、これらは、超音波診断装置内に一体的に具備されている構成としてもよい。
また、プローブは、送受信部の一部の機能をプローブに含んでいてもよい。例えば、送受信部から出力された送信電気信号を生成するための制御信号に基づき、プローブ内で送信電気信号を生成し、この送信電気信号を超音波に変換する。併せて、受信した反射超音波を受信電気信号に変換し、プローブ内で受信電気信号に基づき受信信号を生成する構成を採ることができる。
また、各実施の形態に係る超音波診断装置、及びその変形例の機能のうち少なくとも一部を組み合わせてもよい。更に上記で用いた数字は、全て本発明を具体的に説明するために例示するものであり、本発明は例示された数字に制限されない。
さらに、本実施の形態に対して当業者が思いつく範囲内の変更を施した各種変形例も本発明に含まれる。
≪まとめ≫
実施の形態に係る超音波信号処理装置は、複数の振動子を備えた超音波プローブを用いて被検体に超音波ビームを送信し、被検体から得られた反射波に基づいて音響線信号サブフレームデータを生成する超音波診断装置であって、被検体内における超音波ビームの集束点を規定するフォーカス点を設定し、前記複数の振動子から選択される複数の送信振動子の列から前記集束点に集束する超音波ビームを送信させる送信部と、前記超音波プローブが前記被検体から受波した反射波に基づいて、前記超音波プローブの振動子各々に対する受信信号列を生成する受信部と、前記超音波ビームの主照射領域内の位置に対応する複数の領域内観測点のうち前記フォーカス点を通過する対象線群上に存在する複数の線上観測点について、前記主照射領域内から得られた反射波に基づく前記受信信号列を整相加算して音響線信号サブフレームデータを生成する整相加算部と、を備え、前記受信部は、前記送信振動子の列を分割した複数の部分振動子列の各々について、その部分振動子列に含まれる振動子各々が被検体から受波した反射波に基づき受信信号列を生成する複数のパート受信部と、前記整相加算部は、前記複数の部分振動子列に対応する受信信号列を整相加算することにより前記複数の線上観測点について音響線信号を生成して、前記複数の部分振動子列の各々に対応する音響線信号部分サブフレームデータを生成する複数のパート整相加算部と、前記複数の部分振動子列の各々に対応する前記音響線信号部分サブフレームデータから、前記複数の線上観測点に対する音響線信号の列を抽出し、前記対象線群に含まれる対象線間の間隔を詰めて配列することにより音響線信号部分サブフレーム折り畳みデータを生成する複数のパート折り畳み部と、前記複数の部分振動子列の各々に対応する前記音響線信号部分サブフレーム折り畳みデータを前記配列された観測点の位置を基準に加算して音響線信号サブフレーム折り畳みデータを生成する主加算部と、前記音響線信号フレーム折り畳みデータにおける音響線信号を、前記主照射領域における前記線上観測点の位置に再配列して前記音響線信号サブフレームデータを生成する再配列部とを有することを特徴とする。
係る構成により、空間分解能及びS/N比の低下を抑制しながら整相加算の演算量、及び、生成される音響線信号のデータ量を削減することができ、超音波診断装置の受信ビームフォーマにおける必要な内部メモリ容量やデータ伝送能力を縮減できる。
また、別の態様では、上記何れかの態様において、前記対象線群に含まれる各対象線は直線であり、1つの対象線上であって前記フォーカス点からの距離が所定距離以上である1つの観測点について、前記1つの対象線上にある最近接の観測点との距離は、前記1つの対象線に隣接する対象線上に存在する最近接の観測点との距離より小さい構成としてもよい。
係る構成により、空間分解能及びS/N比の低下を抑制しながら整相加算の演算量及び、生成される音響線信号のデータ量を削減できる。
また、別の態様では、上記何れかの態様において、さらに、前記音響線信号サブフレームデータを複数合成して音響線信号フレームデータを生成する合成部を備え、前記送信部は、前記超音波ビームを送信する送信イベントを前記フォーカス点を変更しながら複数回行い、前記受信部は、各送信イベントに同期して前記受信信号列を生成し、前記整相加算部は、各送信イベントに同期して前記音響線信号サブフレームデータを複数生成し、前記合成部は、各送信イベントに同期して生成された複数の前記音響線信号サブフレームデータを、観測点の位置を基準に合成して音響線信号フレームデータを生成する構成としてもよい。
係る構成により、合成開口法を用いて空間分解能及びS/N比の低下を抑制できる。
また、別の態様では、上記何れかの態様において、前記部分振動子列は、前記送信振動子の列をn分割したものであり、n個の第1の集積回路と第2の集積回路とを備え、前記複数のパート受信部、前記複数のパート整相加算部、前記複数のパート折り畳み部は、それぞれ各1づつn個の第1の集積回路に含まれ、前記主加算部及び前記再配列部は、前記第2の集積回路に含まれ、前記n個の集積回路のそれぞれから、前記第2の集積回路に前記音響線信号部分サブフレーム折り畳みデータが伝送される構成としてもよい。
係る構成により、集積回路501_1、501_2から集積回路501_3へのデータ転送に際し、転送に係るデータ量を縮減することができ、集積回路501_1、501_2、501_3の必要な内部メモリ容量や、集積回路501_1及び501_2〜集積回路501_3間の必要なデータ伝送能力の抑制することができる。
また、受信振動子101aの列をn分割して受信振動子101aに関する整相加算処理を複数の集積回路501_1、501_2、501_3に分散して行い、整相加算処理を複数の演算器に分散化して1演算器あたりの演算量を削減することができる。これにより、部分振動子列毎に対応した前段の整相加算処理、折り畳み処理を行う集積回路501_1、501_2と、後段の主加算処理、再配列処理、合成処理とを行う集積回路501_3とのそれぞれにおいて、例えば、小規模のFPGA(Field Programmable Gate Array)を複数用いることでハードウェアを構成できるので、ハードウェアにおける演算器のコストを大幅に低減することができる。
また、別の態様では、上記何れかの態様において、前記部分振動子列は、前記送信振動子の列をn分割したものであり、n個の集積回路を備え、前記複数のパート受信部、前記複数のパート整相加算部、前記複数のパート折り畳み部は、それぞれ各1づつn個の集積回路に含まれ、前記主加算部及び前記再配列部は、前記n個の集積回路の何れかに含まれ、前記主加算部及び前記再配列部が含まれる集積回路に対し、それ以外の集積回路のそれぞれから、前記音響線信号部分サブフレーム折り畳みデータが伝送される構成としてもよい。
係る構成により、超音波診断装置100Bでは、集積回路を削減することができ、ハードウェアにおける演算器のコストをさらに低減することができる。
また、別の態様では、上記何れかの態様において、さらに、超音波計測条件を示す各種に基づき、前記音響線信号部分サブフレーム折り畳みデータにおける前記対象線上の音響線信号の配列方法を規定する配列情報を決定する配列決定部を備え、前記パート折り畳み部は、前記配列情報に基づき対象線毎に音響線信号の列を配列して音響線信号部分サブフレーム折り畳みデータを生成する構成としてもよい。
係る構成により、超音波計測条件に適合するような態様にて象線毎に音響線信号の列を配列することができ、超音波計測条件に応じて空間分解能及びS/N比の低下を抑制しながら整相加算の演算量及び、生成される音響線信号のデータ量を削減できる。
また、別の態様では、上記何れかの態様において、前記合成部は、各送信イベントに同期して合成された音響線信号フレームデータ間の差分を算出し、各送信イベントに同期して前記差分に相当する音響線信号を出力する構成を採る。
送信イベント間の差分に相当する音響線信号dsとは1個の振動子101aの相当する幅の深さ方向に伸長した音響線信号dsの列であり、音響線信号フレームデータds_fのデータ量を送信イベント数で除したデータ量になる。
係る構成により、送信イベント毎に集積回路501_3から転送されるデータは音響線信号フレームデータds_fのうち送信イベント間の差分に相当する音響線信号dsのみとすることができ、集積回路501_3より後段へのデータ伝送に必要なデータ伝送能力を縮減することができる。
また、本実施の形態に係る超音波信号処理方法は、複数の振動子を備えた超音波プローブを用いて被検体に超音波ビームを送信し、被検体から得られた反射波に基づいて音響線信号サブフレームデータを生成する超音波信号処理方法であって、被検体内における超音波ビームの集束点を規定するフォーカス点を設定し、前記複数の振動子から選択される複数の送信振動子の列から前記集束点に集束する超音波ビームを送信させる送信ステップと、前記超音波プローブが前記被検体から受波した反射波に基づいて、前記超音波プローブの振動子各々に対する受信信号列を生成する受信ステップと、前記超音波ビームの主照射領域内の位置に対応する複数の領域内観測点のうち前記フォーカス点を通過する対象線群上に存在する複数の線上観測点について、前記主照射領域内から得られた反射波に基づく前記受信信号列を整相加算して音響線信号サブフレームデータを生成する整相加算ステップと、を有し、前記受信ステップは、前記送信振動子の列を分割した複数の部分振動子列の各々について、その部分振動子列に含まれる振動子各々が被検体から受波した反射波に基づき受信信号列を生成する複数のパート受信サブステップと、前記整相加算ステップは、前記複数の部分振動子列に対応する受信信号列を整相加算することにより前記複数の線上観測点について音響線信号を生成して、前記複数の部分振動子列の各々に対応する音響線信号部分サブフレームデータを生成する複数のパート整相加算サブステップと、前記複数の部分振動子列の各々に対応する前記音響線信号部分サブフレームデータから、前記複数の線上観測点に対する音響線信号の列を抽出し、前記対象線群に含まれる対象線間の間隔を詰めて配列することにより音響線信号部分サブフレーム折り畳みデータを生成する複数のパート折り畳みサブステップと、前記複数の部分振動子列の各々に対応する前記音響線信号部分サブフレーム折り畳みデータを前記配列された観測点の位置を基準に加算して音響線信号サブフレーム折り畳みデータを生成する主加算サブステップと、前記音響線信号フレーム折り畳みデータにおける音響線信号を、前記主照射領域における前記線上観測点の位置に再配列して前記音響線信号サブフレームデータを生成する再配列サブステップとを有することを特徴とする。
また、別の態様では、上記何れかの態様において、前記対象線群に含まれる各対象線は直線であり、1つの対象線上であって前記フォーカス点からの距離が所定距離以上である1つの観測点について、前記1つの対象線上にある最近接の観測点との距離は、前記1つの対象線に隣接する対象線上に存在する最近接の観測点との距離より小さい構成としてもよい。
係る構成により、超音波診断装置の受信ビームフォーマにおける必要な内部メモリ容量やデータ伝送能力を縮減できる超音波信号処理方法、及び、それを用いた超音波診断装置を実現に供することができる。
係る構成により、空間分解能及びS/N比の低下を抑制しながら整相加算の演算量及び、生成される音響線信号のデータ量を削減できる。
また、別の態様では、上記何れかの態様において、さらに、前記音響線信号サブフレームデータを複数合成して音響線信号フレームデータを生成する合成ステップを有し、
前記送信ステップでは、前記超音波ビームを送信する送信イベントを前記フォーカス点を変更しながら複数回行い、前記受信ステップでは、各送信イベントに同期して前記受信信号列を生成し、前記整相加算ステップでは、各送信イベントに同期して前記音響線信号サブフレームデータを複数生成し、前記合成ステップでは、各送信イベントに同期して生成された複数の前記音響線信号サブフレームデータを、観測点の位置を基準に合成して音響線信号フレームデータを生成する
構成としてもよい。
係る構成により、合成開口法を用いて空間分解能及びS/N比の低下を抑制できる。
≪補足≫
以上で説明した実施の形態は、いずれも本発明の好ましい一具体例を示すものである。実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置及び接続形態、工程、工程の順序などは一例であり、本発明を限定する主旨ではない。また、実施の形態における構成要素のうち、本発明の最上位概念を示す独立請求項に記載されていない工程については、より好ましい形態を構成する任意の構成要素として説明される。
また、上記の工程が実行される順序は、本発明を具体的に説明するために例示するためのものであり、上記以外の順序であってもよい。また、上記工程の一部が、他の工程と同時(並列)に実行されてもよい。
また、発明の理解の容易のため、上記各実施の形態で挙げた各図の構成要素の縮尺は実際のものと異なる場合がある。また本発明は上記各実施の形態の記載によって限定されるものではなく、本発明の要旨を逸脱しない範囲において適宜変更可能である。
本開示にかかる超音波信号処理方法、超音波診断装置は、従来の超音波診断装置の性能向上、特に、演算装置のコスト削減や演算負荷軽減によるフレームレート向上、内部メモリ容量やデータ伝送能力の抑制に有用である。また本開示は超音波への適用のみならず、複数のアレイ素子を用いたセンサ等の用途にも応用できる。
1000 超音波診断システム
100 超音波診断装置
101 プローブ
101a 振動子
101a_l(l=1〜m) 部分振動子列
102 マルチプレクサ部
103 送信ビームフォーマ部
1031 送信部
104 受信ビームフォーマ部
40 受信部
401_l(l=1〜m) パート受信部
41 整相加算部
411_l(l=1〜m) パート整相加算部
4112 対象線群設定部
4113 受信開口設定部
4114 送信時間算出部
4115 受信時間算出部
4116 遅延量算出部
4117 遅延処理部
4118 重み算出部
4119 加算部
412_l(l=1〜m) パート折り畳み部
4121 線上観測点音響線信号抽出部
4122 配列決定部
4123 対象線配列部
413 主加算部
414 再配列部
42 合成部
421 加算処理部
422 増幅処理部
105 超音波画像生成部
106 表示部
107 データ格納部
108 制御部
150 超音波信号処理装置

Claims (10)

  1. 複数の振動子を備えた超音波プローブを用いて被検体に超音波ビームを送信し、被検体から得られた反射波に基づいて音響線信号サブフレームデータを生成する超音波診断装置であって、
    被検体内における超音波ビームの集束点を規定するフォーカス点を設定し、前記複数の振動子から選択される複数の送信振動子の列から前記集束点に集束する超音波ビームを送信させる送信部と、
    前記超音波プローブが前記被検体から受波した反射波に基づいて、前記超音波プローブの振動子各々に対する受信信号列を生成する受信部と、
    前記超音波ビームの主照射領域内の位置に対応する複数の領域内観測点のうち前記フォーカス点を通過する対象線群上に存在する複数の線上観測点について、前記主照射領域内から得られた反射波に基づく前記受信信号列を整相加算して音響線信号サブフレームデータを生成する整相加算部と、
    を備え、
    前記受信部は、前記送信振動子の列を分割した複数の部分振動子列の各々について、その部分振動子列に含まれる振動子各々が被検体から受波した反射波に基づき受信信号列を生成する複数のパート受信部と、
    前記整相加算部は、前記複数の部分振動子列に対応する受信信号列を整相加算することにより前記複数の線上観測点について音響線信号を生成して、前記複数の部分振動子列の各々に対応する音響線信号部分サブフレームデータを生成する複数のパート整相加算部と、
    前記複数の部分振動子列の各々に対応する前記音響線信号部分サブフレームデータから、前記複数の線上観測点に対する音響線信号の列を抽出し、前記対象線群に含まれる対象線間の間隔を詰めて配列することにより音響線信号部分サブフレーム折り畳みデータを生成する複数のパート折り畳み部と、
    前記複数の部分振動子列の各々に対応する前記音響線信号部分サブフレーム折り畳みデータを前記配列された観測点の位置を基準に加算して音響線信号サブフレーム折り畳みデータを生成する主加算部と、
    前記音響線信号フレーム折り畳みデータにおける音響線信号を、前記主照射領域における前記線上観測点の位置に再配列して前記音響線信号サブフレームデータを生成する再配列部とを有する
    超音波診断装置。
  2. 前記対象線群に含まれる各対象線は直線であり、1つの対象線上であって前記フォーカス点からの距離が所定距離以上である1つの観測点について、前記1つの対象線上にある最近接の観測点との距離は、前記1つの対象線に隣接する対象線上に存在する最近接の観測点との距離より小さい
    請求項1に記載の超音波診断装置。
  3. さらに、前記音響線信号サブフレームデータを複数合成して音響線信号フレームデータを生成する合成部を備え、
    前記送信部は、前記超音波ビームを送信する送信イベントを前記フォーカス点を変更しながら複数回行い、
    前記受信部は、各送信イベントに同期して前記受信信号列を生成し、
    前記整相加算部は、各送信イベントに同期して前記音響線信号サブフレームデータを複数生成し、
    前記合成部は、各送信イベントに同期して生成された複数の前記音響線信号サブフレームデータを、観測点の位置を基準に合成して音響線信号フレームデータを生成する
    請求項1又は2に記載の超音波診断装置。
  4. 前記部分振動子列は、前記送信振動子の列をn分割(nは1より大きい自然数)したものであり、
    n個の第1の集積回路と第2の集積回路とを備え、
    前記複数のパート受信部、前記複数のパート整相加算部、前記複数のパート折り畳み部は、それぞれ各1づつn個の第1の集積回路に含まれ、
    前記主加算部及び前記再配列部は、前記第2の集積回路に含まれ、
    前記n個の集積回路のそれぞれから、前記第2の集積回路に前記音響線信号部分サブフレーム折り畳みデータが伝送される
    請求項1から3の何れか1項に記載の超音波診断装置。
  5. 前記部分振動子列は、前記送信振動子の列をn分割(nは1より大きい自然数)したものであり、
    n個の集積回路を備え、
    前記複数のパート受信部、前記複数のパート整相加算部、前記複数のパート折り畳み部は、それぞれ各1づつn個の集積回路に含まれ、
    前記主加算部及び前記再配列部は、前記n個の集積回路の何れかに含まれ、
    前記主加算部及び前記再配列部が含まれる集積回路に対し、それ以外の集積回路のそれぞれから、前記音響線信号部分サブフレーム折り畳みデータが伝送される
    請求項1から3の何れか1項に記載の超音波診断装置。
  6. さらに、超音波計測条件を示す各種に基づき、前記音響線信号部分サブフレーム折り畳みデータにおける前記対象線上の音響線信号の配列方法を規定する配列情報を決定する配列決定部を備え、
    前記パート折り畳み部は、前記配列情報に基づき、対象線毎に音響線信号の列を配列して音響線信号部分サブフレーム折り畳みデータを生成する
    請求項1から5の何れか1項に記載の超音波診断装置。
  7. 前記合成部は、各送信イベントに同期して合成された音響線信号フレームデータ間の差分を算出し、各送信イベントに同期して前記差分に相当する音響線信号を出力する
    請求項4に記載の超音波診断装置。
  8. 複数の振動子を備えた超音波プローブを用いて被検体に超音波ビームを送信し、被検体から得られた反射波に基づいて音響線信号サブフレームデータを生成する超音波信号処理方法であって、
    被検体内における超音波ビームの集束点を規定するフォーカス点を設定し、前記複数の振動子から選択される複数の送信振動子の列から前記集束点に集束する超音波ビームを送信させる送信ステップと、
    前記超音波プローブが前記被検体から受波した反射波に基づいて、前記超音波プローブの振動子各々に対する受信信号列を生成する受信ステップと、
    前記超音波ビームの主照射領域内の位置に対応する複数の領域内観測点のうち前記フォーカス点を通過する対象線群上に存在する複数の線上観測点について、前記主照射領域内から得られた反射波に基づく前記受信信号列を整相加算して音響線信号サブフレームデータを生成する整相加算ステップと、
    を有し、
    前記受信ステップは、前記送信振動子の列を分割した複数の部分振動子列の各々について、その部分振動子列に含まれる振動子各々が被検体から受波した反射波に基づき受信信号列を生成する複数のパート受信サブステップと、
    前記整相加算ステップは、前記複数の部分振動子列に対応する受信信号列を整相加算することにより前記複数の線上観測点について音響線信号を生成して、前記複数の部分振動子列の各々に対応する音響線信号部分サブフレームデータを生成する複数のパート整相加算サブステップと、
    前記複数の部分振動子列の各々に対応する前記音響線信号部分サブフレームデータから、前記複数の線上観測点に対する音響線信号の列を抽出し、前記対象線群に含まれる対象線間の間隔を詰めて配列することにより音響線信号部分サブフレーム折り畳みデータを生成する複数のパート折り畳みサブステップと、
    前記複数の部分振動子列の各々に対応する前記音響線信号部分サブフレーム折り畳みデータを前記配列された観測点の位置を基準に加算して音響線信号サブフレーム折り畳みデータを生成する主加算サブステップと、
    前記音響線信号フレーム折り畳みデータにおける音響線信号を、前記主照射領域における前記線上観測点の位置に再配列して前記音響線信号サブフレームデータを生成する再配列サブステップとを有する
    超音波信号処理方法。
  9. 前記対象線群に含まれる各対象線は直線であり、1つの対象線上であって前記フォーカス点からの距離が所定距離以上である1つの観測点について、前記1つの対象線上にある最近接の観測点との距離は、前記1つの対象線に隣接する対象線上に存在する最近接の観測点との距離より小さい
    請求項8に記載の超音波信号処理方法。
  10. さらに、前記音響線信号サブフレームデータを複数合成して音響線信号フレームデータを生成する合成ステップを有し、
    前記送信ステップでは、前記超音波ビームを送信する送信イベントを前記フォーカス点を変更しながら複数回行い、
    前記受信ステップでは、各送信イベントに同期して前記受信信号列を生成し、
    前記整相加算ステップでは、各送信イベントに同期して前記音響線信号サブフレームデータを複数生成し、
    前記合成ステップでは、各送信イベントに同期して生成された複数の前記音響線信号サブフレームデータを、観測点の位置を基準に合成して音響線信号フレームデータを生成する
    請求項8又は9に記載の超音波信号処理方法。
JP2017085699A 2017-04-24 2017-04-24 超音波診断装置、及び超音波信号処理方法。 Active JP6825474B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2017085699A JP6825474B2 (ja) 2017-04-24 2017-04-24 超音波診断装置、及び超音波信号処理方法。
US15/942,634 US10702246B2 (en) 2017-04-24 2018-04-02 Ultrasound diagnostic apparatus and an ultrasound signal processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017085699A JP6825474B2 (ja) 2017-04-24 2017-04-24 超音波診断装置、及び超音波信号処理方法。

Publications (2)

Publication Number Publication Date
JP2018183297A JP2018183297A (ja) 2018-11-22
JP6825474B2 true JP6825474B2 (ja) 2021-02-03

Family

ID=63852526

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017085699A Active JP6825474B2 (ja) 2017-04-24 2017-04-24 超音波診断装置、及び超音波信号処理方法。

Country Status (2)

Country Link
US (1) US10702246B2 (ja)
JP (1) JP6825474B2 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6747108B2 (ja) * 2016-07-05 2020-08-26 コニカミノルタ株式会社 超音波信号処理装置、超音波信号処理方法、及び、超音波診断装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3910860B2 (ja) * 2002-02-05 2007-04-25 株式会社日立メディコ 超音波撮像装置
US9060669B1 (en) * 2007-12-20 2015-06-23 Zonare Medical Systems, Inc. System and method for providing variable ultrasound array processing in a post-storage mode
CN103635829B (zh) * 2011-06-30 2016-04-27 皇家飞利浦有限公司 具有两个波束成形器阶段的二维超声诊断成像系统
JP6380194B2 (ja) * 2015-03-27 2018-08-29 コニカミノルタ株式会社 超音波信号処理装置、及び、超音波診断装置

Also Published As

Publication number Publication date
JP2018183297A (ja) 2018-11-22
US10702246B2 (en) 2020-07-07
US20180303461A1 (en) 2018-10-25

Similar Documents

Publication Publication Date Title
US11360199B2 (en) Ultrasound signal processing device, ultrasound diagnostic device
CN108209971B (zh) 超声波信号处理装置和方法以及超声波诊断装置
JP6406019B2 (ja) 超音波信号処理装置、及び超音波診断装置
JP2010057730A (ja) 光音響トモグラフィの受信データ処理装置
JP6579231B2 (ja) 超音波画像生成方法、及び超音波診断装置
CN107569254B (zh) 超声波信号处理装置、超声波信号处理方法以及超声波诊断装置
US20180055486A1 (en) Ultrasound signal processing device, ultrasound diagnostic device, and ultrasound signal processing method
JP6825474B2 (ja) 超音波診断装置、及び超音波信号処理方法。
JP7147399B2 (ja) 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法
US11744555B2 (en) Ultrasound signal processing device, ultrasound diagnostic device, and ultrasound signal processing method
US11484292B2 (en) Ultrasound signal processing device that uses synthetic aperture method and delay and sum method
JP6848793B2 (ja) 超音波信号処理方法、及び超音波信号処理装置。
JP2018082835A (ja) 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法
US11129595B2 (en) Ultrasonic diagnostic apparatus, interpolation processing unit, and interpolation processing method
JP6562122B2 (ja) 超音波診断装置、及び、超音波画像生成方法
JP6863817B2 (ja) 超音波撮像装置
JP7020052B2 (ja) 超音波信号処理装置、超音波診断装置、超音波信号処理方法、および、超音波画像表示方法
JP5627756B2 (ja) 光音響装置の制御方法
JP6933102B2 (ja) 超音波信号処理装置、及び超音波信号処理方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200318

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: 20201215

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20201216

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20201228

R150 Certificate of patent or registration of utility model

Ref document number: 6825474

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150