JP7095177B2 - 超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ - Google Patents

超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ Download PDF

Info

Publication number
JP7095177B2
JP7095177B2 JP2021514934A JP2021514934A JP7095177B2 JP 7095177 B2 JP7095177 B2 JP 7095177B2 JP 2021514934 A JP2021514934 A JP 2021514934A JP 2021514934 A JP2021514934 A JP 2021514934A JP 7095177 B2 JP7095177 B2 JP 7095177B2
Authority
JP
Japan
Prior art keywords
blood vessel
period
diameter
frame
data
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
JP2021514934A
Other languages
English (en)
Other versions
JPWO2020213562A1 (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.)
Fujifilm Corp
Original Assignee
Fujifilm Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujifilm Corp filed Critical Fujifilm Corp
Publication of JPWO2020213562A1 publication Critical patent/JPWO2020213562A1/ja
Application granted granted Critical
Publication of JP7095177B2 publication Critical patent/JP7095177B2/ja
Active 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/54Control of the diagnostic device
    • A61B8/543Control of the diagnostic device involving acquisition triggered by a physiological signal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0891Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/488Diagnostic techniques involving Doppler signals
    • 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
    • A61B8/5223Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • 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
    • A61B8/5238Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5246Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
    • 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

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Public Health (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Hematology (AREA)
  • Vascular Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

本発明は、Bモードデータとドプラデータを取得する超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサに関する。
従来から、被検体の内部の画像を得るものとして、超音波診断装置が知られている。超音波診断装置は、一般的に、複数の素子が配列された振動子アレイが備えられた超音波プローブを備えている。この超音波プローブを被検体の体表に接触させた状態において、振動子アレイから被検体内に向けて超音波ビームが送信され、被検体からの超音波エコーを振動子アレイにおいて受信して素子データが取得される。さらに、超音波診断装置は、得られた素子データを電気的に処理して、被検体の当該部位に対する超音波画像を生成する。
例えば、特許文献1には、ドプラデータに基づいて、被検体の血管が写るBモード(Brightness mode:輝度モード)画像における血管領域を検出する超音波診断装置が開示されている。特許文献1の超音波診断装置は、検出された血管領域に基づいてBモード画像における血管をセグメント化し、セグメント化された血管の直径をBモード画像に基づいて算出する。
特表2017-524455号公報
ここで、一般的に、血管の弾性を表す弾性指標等を算出するために、血管の最大径および最小径が測定されることがある。特許文献1に記載されているような従来の超音波診断装置を用いて、血管の最小径および最大径を測定する場合には、連続的に取得された複数フレームのそれぞれのBモード画像から算出された血管の直径を比較して、その最小値および最大値を求める必要があるため、超音波診断装置における負担が大きく、血管の最小径および最大径を算出するために多大な時間を要するという問題があった。また、Bモード画像のみに基づいて血管の最小径および最大径を測定しようとすると、ノイズの影響等に起因して輝度値に変動が生じた場合に、誤った最小値および最大値を測定してしまうおそれがあった。
本発明は、このような従来の問題点を解決するためになされたものであり、血管の直径が最大となるBモードデータおよび血管の直径が最小となるBモードデータの少なくとも一方を容易に且つ正確に特定することができる超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサを提供することを目的とする。
上記目的を達成するために、本発明に係る超音波診断装置は、被検体の血管が含まれる領域に対して定められた期間にわたってBモードデータとドプラデータとが連続的に取得される超音波診断装置であって、ドプラデータに基づいて各心拍期間における時相サーチ期間を特定する時相サーチ期間特定部と、時相サーチ期間特定部により特定された時相サーチ期間内における複数フレームのBモードデータを解析することにより、各心拍期間において血管の直径が最大となるフレームのBモードデータおよび血管の直径が最小となるフレームのBモードデータの少なくとも一方を特定するフレーム特定部とを備えることを目的とする。
フレーム特定部により特定された、血管の直径が最大となるフレームのBモードデータおよび血管の直径が最小となるフレームのBモードデータの少なくとも一方に基づいて血管の最大径および最小径の少なくとも一方を算出する血管径算出部をさらに備えることが好ましい。
また、血管径算出部により算出された血管の最大径および最小径の少なくとも一方を用いて血管の断面積を算出する断面積算出部と、各心拍期間におけるドプラデータを取得するドプラ処理部と、ドプラ処理部により取得されたドプラデータに基づいて血流速度を算出する血流速度算出部と、断面積算出部により算出された断面積と血流速度算出部により算出された血流速度に基づいて血流量を計測する血流量計測部とをさらに備えることが好ましい。
時相サーチ期間特定部は、各心拍期間におけるドプラデータが最小値を有する時点を含む第1期間を時相サーチ期間として特定し、フレーム特定部は、第1期間内における複数フレームのBモードデータに基づいて血管の直径が最小となるフレームのBモードデータを特定し、血管径算出部は、フレーム特定部により特定された、血管の直径が最小となるフレームのBモードデータに基づいて、血管の最小径を算出し、断面積算出部は、血管径算出部により算出された血管の最小径を用いて血管の断面積を算出することができる。
もしくは、時相サーチ期間特定部は、各心拍期間におけるドプラデータが最大値を有する時点を含む第2期間を時相サーチ期間として特定し、フレーム特定部は、第2期間内における複数フレームのBモードデータに基づいて血管の直径が最大となるフレームのBモードデータを特定し、血管径算出部は、フレーム特定部により特定された、血管の直径が最大となるフレームのBモードデータに基づいて、血管の最大径を算出し、断面積算出部は、血管径算出部により算出された血管の最大径を用いて血管の断面積を算出することもできる。
もしくは、時相サーチ期間特定部は、各心拍期間におけるドプラデータが最小値を有する時点を含む第1期間および各心拍期間におけるドプラデータが最大値を有する時点を含む第2期間を有する時相サーチ期間を特定し、フレーム特定部は、第1期間内における複数フレームのBモードデータに基づいて血管の直径が最小となるフレームのBモードデータを特定し、且つ、第2期間内における複数フレームのBモードデータに基づいて血管の直径が最大となるフレームのBモードデータを特定し、血管径算出部は、フレーム特定部により特定された、血管の直径が最小となるフレームのBモードデータに基づいて、血管の最小径を算出し、且つ、フレーム特定部により特定された、血管の直径が最大となるフレームのBモードデータに基づいて、血管の最大径を算出し、断面積算出部は、血管径算出部により算出された血管の最小径および血管の最大径とから算出された各心拍期間における血管の平均径を用いて血管の断面積を算出することもできる。
第1期間は、各心拍期間におけるドプラデータが最小値を有する時点から最大値を有する時点までの期間であることが好ましい。
また、第2期間は、各心拍期間におけるドプラデータが最大値を有する期間から各心拍期間よりも短い時間に設定された定められた時間が経過するまでの期間であることが好ましい。
さらに、第1期間および第2期間は、それぞれ、各心拍期間の10%以上20%以下の時間幅を有することが好ましい。
また、血管の最大径と血管の最小径との差に基づいて血管の弾性指標を算出する弾性指標算出部をさらに備えることができる。
また、超音波プローブと、超音波プローブを介して被検体内に超音波ビームを送信し且つ被検体内から超音波エコーを受信して受信データを生成する送受信回路と、送受信回路により生成された受信データに基づいてBモードデータを生成するBモード処理部とをさらに備えることが好ましい。
本発明の超音波診断装置の制御方法は、被検体の血管が含まれる領域に対して定められた期間にわたってBモードデータとドプラデータとが連続的に取得される超音波診断装置の制御方法であって、ドプラデータに基づいて各心拍期間における時相サーチ期間を特定し、特定された時相サーチ期間内における複数フレームのBモードデータを解析することにより、各心拍期間において血管の直径が最大となるフレームのBモードデータおよび血管の直径が最小となるフレームのBモードデータの少なくとも一方を特定することを特徴とする。
本発明の超音波診断装置用プロセッサは、被検体の血管が含まれる領域に対して定められた期間にわたってBモードデータとドプラデータとを連続的に取得し、ドプラデータに基づいて各心拍期間における時相サーチ期間を特定し、特定された時相サーチ期間内における複数フレームのBモードデータを解析することにより、各心拍期間において血管の直径が最大となるフレームのBモードデータおよび血管の直径が最小となるフレームのBモードデータの少なくとも一方を特定することを特徴とする。
本発明によれば、ドプラデータに基づいて各心拍期間における時相サーチ期間を特定する時相サーチ期間特定部と、時相サーチ期間特定部により特定された時相サーチ期間内における複数フレームのBモードデータを解析することにより、各心拍期間において血管の直径が最大となるフレームのBモードデータおよび血管の直径が最小となるフレームのBモードデータの少なくとも一方を特定するフレーム特定部とを備えるため、血管の直径が最大となるBモードデータと血管の直径が最小となるBモードデータの少なくとも一方を容易に且つ正確に特定することができる。
本発明の実施の形態1に係る超音波診断装置の構成を示すブロック図である。 本発明の実施の形態1における受信回路の内部構成を示すブロック図である。 本発明の実施の形態1におけるBモード処理部の内部構成を示すブロック図である。 本発明の実施の形態1におけるドプラ処理部の内部構成を示すブロック図である。 本発明の実施の形態1におけるドプラ波形画像を模式的に示す図である。 本発明の実施の形態1においてドプラゲートが設定されたBモード画像の例を模式的に示す図である。 本発明の実施の形態1における理想的なドプラ波形画像と血管径の理想的な時間変化の関係を模式的に示す図である。 本発明の実施の形態1における実際のドプラ波形画像の例を示す図である。 本発明の実施の形態1において血管径を算出する方法を模式的に示す図である。 本発明の実施の形態1に係る超音波診断装置の動作を示すフローチャートである。 本発明の実施の形態1において表示装置に表示されたBモード画像とドプラ波形画像を模式的に示す図である。 本発明の実施の形態1において表示装置に表示されたBモード画像、ドプラ波形画像、血流量の計測値を模式的に示す図である。 本発明の実施の形態2に係る超音波診断装置を示すブロック図である。 本発明の実施の形態3に係る超音波診断装置を示すブロック図である。
以下、この発明の実施の形態を添付図面に基づいて説明する。
以下に記載する構成要件の説明は、本発明の代表的な実施態様に基づいてなされるが、本発明はそのような実施態様に限定されるものではない。
なお、本明細書において、「~」を用いて表される数値範囲は、「~」の前後に記載される数値を下限値および上限値として含む範囲を意味する。
また、本明細書において、「垂直」および「平行」とは、本発明が属する技術分野において許容される誤差の範囲を含むものとする。例えば、「垂直」および「平行」とは、厳密な垂直あるいは平行に対して±10°未満の範囲内であることなどを意味し、厳密な垂直あるいは平行に対しての誤差は、5°以下であることが好ましく、3°以下であることがより好ましい。
本明細書において、「同一」、「同じ」は、技術分野で一般的に許容される誤差範囲を含むものとする。また、本明細書において、「全部」、「いずれも」または「全面」などというとき、100%である場合のほか、技術分野で一般的に許容される誤差範囲を含み、例えば99%以上、95%以上、または90%以上である場合を含むものとする。
実施の形態1
図1に、本発明の実施の形態1に係る超音波診断装置1の構成を示す。図1に示すように、超音波診断装置1は、振動子アレイ2を備えており、送信回路3および受信回路4がそれぞれ接続されている。ここで、送信回路3と受信回路4により、送受信回路5が構成されている。受信回路4には、Bモード(Brightness mode:輝度モード)処理部6およびドプラ処理部7が接続され、これらのBモード処理部6およびドプラ処理部7に表示制御部8を介して表示装置9が接続されている。
また、Bモード処理部6にゲート設定部10が接続され、ゲート設定部10に、ドプラ処理部7および表示制御部8が接続されている。また、Bモード処理部6に、画像メモリ11が接続され、画像メモリ11に、フレーム特定部12が接続されている。また、フレーム特定部12に、血管径算出部13、断面積算出部14、血流量計測部15が順次接続されている。また、ドプラ処理部7に、時相サーチ期間特定部16と血流速度算出部17が接続されている。時相サーチ期間特定部16に、フレーム特定部12が接続され、血流速度算出部17に、血流量計測部15が接続されている。
また、送受信回路5、Bモード処理部6、ドプラ処理部7、表示制御部8、ゲート設定部10、フレーム特定部12、血管径算出部13、断面積算出部14、血流量計測部15、時相サーチ期間特定部16、血流速度算出部17に、装置制御部18が接続されており、装置制御部18に、入力装置19および格納部20が接続されている。装置制御部18と格納部20は、互いに双方向の情報の受け渡しが可能に接続されている。
また、振動子アレイ2は、超音波プローブ21に含まれており、Bモード処理部6、ドプラ処理部7、表示制御部8、ゲート設定部10、フレーム特定部12、血管径算出部13、断面積算出部14、血流量計測部15、時相サーチ期間特定部16、血流速度算出部17および装置制御部18により、超音波診断装置1用のプロセッサ22が構成されている。
図1に示す超音波プローブ21の振動子アレイ2は、1次元または2次元に配列された複数の振動子を有している。これらの振動子は、それぞれ送信回路3から供給される駆動信号に従って超音波を送信すると共に、被検体からの超音波エコーを受信して、超音波エコーに基づく信号を出力する。各振動子は、例えば、PZT(Lead Zirconate Titanate:チタン酸ジルコン酸鉛)に代表される圧電セラミック、PVDF(Poly Vinylidene Di Fluoride:ポリフッ化ビニリデン)に代表される高分子圧電素子およびPMN-PT(Lead Magnesium Niobate-Lead Titanate:マグネシウムニオブ酸鉛-チタン酸鉛固溶体)に代表される圧電単結晶等からなる圧電体の両端に電極を形成することにより構成される。
送信回路3は、例えば、複数のパルス発生器を含んでおり、装置制御部18からの制御信号に応じて選択された送信遅延パターンに基づいて、振動子アレイ2の複数の振動子から送信される超音波が超音波ビームを形成するようにそれぞれの駆動信号を、遅延量を調節して複数の振動子に供給する。このように、振動子アレイ2の振動子の電極にパルス状または連続波状の電圧が印加されると、圧電体が伸縮し、それぞれの振動子からパルス状または連続波状の超音波が発生して、それらの超音波の合成波から、超音波ビームが形成される。
送信された超音波ビームは、例えば、被検体の部位等の対象において反射され、超音波プローブ21の振動子アレイ2に向かって伝搬する。このように振動子アレイ2に向かって伝搬する超音波は、振動子アレイ2を構成するそれぞれの振動子により受信される。この際に、振動子アレイ2を構成するそれぞれの振動子は、伝搬する超音波エコーを受信することにより伸縮して電気信号を発生させ、これらの電気信号を受信回路4に出力する。
受信回路4は、装置制御部18からの制御信号に従い、振動子アレイ2から出力される信号の処理を行って、いわゆるRF(Radio Frequency:高周波)データである、受信データを生成する。図2に示すように、受信回路4は、増幅部23、AD(Analog Digital:アナログデジタル)変換部24およびビームフォーマ25が直列に接続された構成を有している。
増幅部23は、振動子アレイ2を構成するそれぞれの振動子から入力された信号を増幅し、増幅した信号をAD変換部24に送信する。AD変換部24は、増幅部23から送信された信号をデジタルデータに変換し、これらのデータをビームフォーマ25に送信する。ビームフォーマ25は、装置制御部18からの制御信号に応じて選択された受信遅延パターンに基づいて設定される音速または音速の分布に従い、AD変換部24により変換された各データに対してそれぞれの遅延を与えて加算することにより、いわゆる受信フォーカス処理を行う。この受信フォーカス処理により、AD変換部24により変換された各データが整相加算され且つ超音波エコーの焦点が絞り込まれた受信データが取得される。
Bモード処理部6は、図3に示されるように、信号処理部26、DSC(Digital Scan Converter:デジタルスキャンコンバータ)27および画像処理部28が順次直列に接続された構成を有している。
信号処理部26は、受信回路4により生成された受信データに対し、超音波の反射位置の深度に応じて距離による減衰の補正を施した後、包絡線検波処理を施すことにより、被検体内の組織に関する断層画像情報であるBモード画像信号を生成する。
DSC27は、信号処理部26で生成されたBモード画像信号を通常のテレビジョン信号の走査方式に従う画像信号に変換(ラスター変換)する。
画像処理部28は、DSC27から入力されるBモード画像信号に階調処理等の各種の必要な画像処理を施した後、Bモード画像信号を表示制御部8に出力する。
ドプラ処理部7は、後述するゲート設定部10により血管領域に設定されるドプラゲート内のドプラデータを取得し、ドプラ波形画像を生成するものであり、図4に示されるように、直交検波部29とハイパスフィルタ30と高速フーリエ変換部(Fast Fourier Transformer)31とドプラ波形画像生成部32が順次直列に接続されると共に直交検波部29の出力端にデータメモリ33が接続された構成を有している。
直交検波部29は、受信回路4で生成された受信データに参照周波数のキャリア信号を混合することで、受信データを直交検波して複素データに変換する。
ハイパスフィルタ30は、いわゆるウォールフィルタ(Wall Filter)として機能するもので、直交検波部29で生成された複素データから被検体の体内組織の運動に由来する周波数成分を除去する。
高速フーリエ変換部31は、複数のサンプル点の複素データをフーリエ変換することにより周波数解析して血流速度を求め、スペクトル信号を生成する。
ドプラ波形画像生成部32は、高速フーリエ変換部31で生成されたスペクトル信号を時間軸上に揃えつつ各周波数成分の大きさを輝度で表すことによりドプラ波形画像信号を生成する。ここで、図5に、ドプラ波形画像信号に基づく理想的なドプラ波形画像UDの例を示す。ドプラ波形画像UDは、横軸に時間軸を示し、縦軸にドプラシフト周波数すなわち血流速度を示し、波形の輝度が各周波数成分におけるパワーを表すものであり、ドプラ波形画像UDにおける血流速度の値は、心拍期間HCに従って周期的に変化する。心拍期間HCについては、血流速度が最小値V1を有する時点T1から次に最小値V1を有する時点T4までの期間を心拍期間HCと定義することとする。
また、データメモリ33は、直交検波部29で受信データから変換された複素データを保存する。
装置制御部18は、格納部20等に予め記憶されているプログラムおよび入力装置19を介したユーザによる入力操作に基づいて、超音波診断装置1の各部の制御を行う。
表示制御部8は、装置制御部18の制御の下、Bモード処理部6により生成されたBモード画像信号、ドプラ処理部7により生成されたドプラ波形画像信号等に所定の処理を施して、Bモード画像、ドプラ波形画像等を表示装置9に表示する。
超音波診断装置1の表示装置9は、表示制御部8により生成された画像を表示するものであり、例えば、LCD(Liquid Crystal Display:液晶ディスプレイ)、有機ELディスプレイ(Organic Electroluminescence Display)等のディスプレイ装置を含む。
超音波診断装置1の入力装置19は、ユーザが入力操作を行うためのものであり、キーボード、マウス、トラックボール、タッチパッドおよびタッチパネル等を備えて構成することができる。
ゲート設定部10は、入力装置19を介したユーザの入力操作に基づいて、図6に示すように、Bモード画像UB上の血管領域BR内にドプラゲートDGを設定する。例えば、ユーザが、表示装置9に表示されたBモード画像UBを見ながら、入力装置19を介してBモード画像UB上の血管領域BR内における適切な位置を指定すると、ゲート設定部10は、ユーザにより指定された位置にドプラゲートDGを配置する。このようにして設定されたドプラゲートDGは、Bモード画像UBに重畳して表示装置9に表示される。
ドプラ処理部7は、ゲート設定部10により設定されたドプラゲートDG内のドプラデータを取得する。
ところで、一般的に、超音波診断装置を用いて被検体内の血管の最小径および最大径を算出する場合には、各心拍期間HCにわたって連続的に生成された複数フレームのBモード画像のすべてに対して、算出された血管の直径を比較する必要があった。そのため、血管の最小径および最大径を算出する際に、超音波診断装置における負担が大きく、多大な時間を要するという問題があった。
そこで、本発明者らは、ドプラ波形画像UDにおける血流速度の時間変化と血管の直径の時間変化に着目して、図7に示すような関係を見出した。血管の直径は、血流速度と同様に、心拍期間HCに従って最小値D1と最大値D2の間を周期的に変化するが、血流速度が最小値V1を有する時点T1を含む期間、すなわち、血流速度が最小値V1を有する時点T1から最大値V2を有する時点T2までの第1期間FPにおいて血管の直径が最小値D1を有する。また、血流速度が最大値V2を有する時点T2を含む期間、すなわち、血流速度が最大値V2を有する時点T2から、各心拍期間HCよりも短い時間、例えば0.25秒等に設定された定められた時間が経過する時点T3までの第2期間SPにおいて血管の直径が最大値D2を有する。
そのため、第1期間FPおよび第2期間SPを特定することにより、超音波診断装置1の負担を軽減して、短時間に且つ容易に、血管の直径が最小値D1を有するBモード画像および血管の直径が最大値D2を有するBモード画像を特定することができる。さらに、血管の直径が最小値D1を有するBモード画像と血管の直径が最大値D2を有するBモード画像を特定する際に、第1期間FP内および第2期間SP内のBモード画像における血管の直径を互いに比較すればよく、ノイズの影響等に起因して輝度値に変動が生じたBモード画像を誤って特定するおそれを低減することができる。
時相サーチ期間特定部16は、ドプラ処理部7により生成されたドプラデータに基づいて、各心拍期間HCにおける第1期間FPと第2期間SPの少なくとも一方を時相サーチ期間として特定する。ここで、時相サーチ期間とは、各心拍期間HCにおいて血管の直径が最大となるBモード画像信号と血管の直径が最小となるBモード画像信号の少なくとも一方を特定するための期間である。
なお、時相サーチ期間特定部16は、図8に示すように、各心拍期間HCにおいてドプラ波形WDの包絡線Eを設定し、設定された包絡線Eに基づいて時相サーチ期間を特定するものとする。
画像メモリ11は、Bモード処理部6により連続的に生成された複数フレームのBモード画像信号を保存するメモリである。画像メモリ11としては、フラッシュメモリ、HDD(Hard Disc Drive:ハードディスクドライブ)、SSD(Solid State Drive:ソリッドステートドライブ)、FD(Flexible Disc:フレキシブルディスク)、MOディスク(Magneto-Optical disc:光磁気ディスク)、MT(Magnetic Tape:磁気テープ)、RAM(Random Access Memory:ランダムアクセスメモリ)、CD(Compact Disc:コンパクトディスク)、DVD(Digital Versatile Disc:デジタルバーサタイルディスク)、SDカード(Secure Digital card:セキュアデジタルカード)、USBメモリ(Universal Serial Bus memory:ユニバーサルシリアルバスメモリ)等の記録メディア、またはサーバ等を用いることができる。
フレーム特定部12は、画像メモリ11に保存された複数フレームのBモード画像信号のうち、時相サーチ期間特定部16により特定された時相サーチ期間内における複数フレームのBモード画像信号を解析することにより、各心拍期間HCにおいて血管の直径が最大となるフレームのBモード画像信号および血管の直径が最小となるフレームのBモード画像信号の少なくとも一方を特定する。この際に、フレーム特定部12は、例えば、時相サーチ期間内の各フレームのBモード画像信号において血管の直径を算出し、算出された血管の直径に基づいて、血管の直径が最大となるフレームのBモード画像信号および血管の直径が最小となるフレームのBモード画像信号の少なくとも一方を特定することができる。
フレーム特定部12は、各フレームのBモード画像信号において血管の直径を算出する際に、例えば、図9に示すように、入力装置19を介したユーザの入力操作等によりBモード画像UBに設定された縦方向の直線SL上におけるBモード画像UBの輝度が、一定の値よりも高い2点の位置を血管前壁W1の位置および血管後壁W2の位置として特定し、特定された2点間のBモード画像UB上の距離Lを算出することができる。なお、例えば、直線SLとして、Bモード画像UBの中心を通る垂直線を使用してもよい。
フレーム特定部12は、時相サーチ期間の第1期間FP内において、このようにして算出された血管の直径が最小となるフレームのBモード画像信号を特定する。また、フレーム特定部12は、時相サーチ期間の第2期間SPにおいて、このようにして算出された血管の直径が最大となるフレームのBモード画像信号を特定する。
血管径算出部13は、フレーム特定部12により特定された、血管の直径が最小となるフレームのBモード画像信号および血管の直径が最大となるフレームのBモード画像信号の少なくとも一方に基づいて、血管の最大径および最小径の少なくとも一方を算出する。例えば、血管径算出部13は、血管の直径が最大となるフレームのBモード画像信号において、フレーム特定部12により血管の直径として算出されたBモード画像UB上の距離Lを、血管の実際の直径に換算する等により、血管の最大径を算出することができる。また、例えば、血管径算出部13は、血管の直径が最小となるフレームのBモード画像信号において、フレーム特定部12により血管の直径として算出されたBモード画像UB上の距離Lに基づいて、血管の最小径を算出することができる。
断面積算出部14は、血管が円形の断面を有するものとして、血管径算出部13により算出された血管の直径に基づいて、血管の断面積を算出する。例えば、血管径算出部13により、血管の最小径が算出された場合に、断面積算出部14は、血管の最小径を用いて血管の断面積を算出する。また、例えば、血管径算出部13により、血管の最大径が算出された場合に、断面積算出部14は、血管の最大径を用いて血管の断面積を算出する。また、例えば、血管径算出部13により、血管の最小径および最大径の双方が算出された場合には、断面積算出部14は、血管の最小径と最大径に基づいて血管の平均径を算出し、算出された血管の平均径を用いて血管の断面積を算出することができる。
血流速度算出部17は、ドプラ処理部7により取得されたドプラデータに基づいて、いわゆるパルスドプラ法により、血流速度を算出する。なお、血流速度算出部17は、各心拍期間HCにおける平均血流速度を算出することもできる。
血流量計測部15は、断面積算出部14により算出された血管の断面積と、血流速度算出部17により算出された血流速度とに基づいて、血管内を流れる血液の単位時間当たりの体積を表す血流量を計測する。
格納部20は、超音波診断装置1の動作プログラム等を格納するもので、フラッシュメモリ、HDD、SSD、FD、MOディスク、MT、RAM、CD、DVD、SDカード、USBメモリ等の記録メディア、またはサーバ等を用いることができる。
なお、Bモード処理部6、ドプラ処理部7、表示制御部8、ゲート設定部10、フレーム特定部12、血管径算出部13、断面積算出部14、血流量計測部15、時相サーチ期間特定部16、血流速度算出部17および装置制御部18を有するプロセッサ22は、CPU(Central Processing Unit:中央処理装置)、および、CPUに各種の処理を行わせるための制御プログラムから構成されるが、FPGA(Field Programmable Gate Array:フィードプログラマブルゲートアレイ)、DSP(Digital Signal Processor:デジタルシグナルプロセッサ)、ASIC(Application Specific Integrated Circuit:アプリケーションスペシフィックインテグレイテッドサーキット)、GPU(Graphics Processing Unit:グラフィックスプロセッシングユニット)、その他のIC(Integrated Circuit:集積回路)を用いて構成されてもよく、もしくはそれらを組み合わせて構成されてもよい。
また、プロセッサ22のBモード処理部6、ドプラ処理部7、表示制御部8、ゲート設定部10、フレーム特定部12、血管径算出部13、断面積算出部14、血流量計測部15、時相サーチ期間特定部16、血流速度算出部17および装置制御部18は、部分的にあるいは全体的に1つのCPU等に統合させて構成することもできる。
次に、図10に示すフローチャートを用いて、実施の形態1における超音波診断装置1の動作を詳細に説明する。
まず、ステップS1において、図6に示すように、入力装置19を介したユーザの入力操作に基づいて、ゲート設定部10により、Bモード画像UB上にドプラゲートDGが設定される。
次に、ステップS2において、Bモード処理部6とドプラ処理部7により、心拍期間HCよりも長い定められた期間にわたって、Bモード画像信号とドプラ波形画像信号が生成される。この際に、Bモード処理部6は、少なくとも血管領域BRが撮像されているBモード画像信号を順次、連続的に生成し、生成されたBモード画像信号に基づくBモード画像UBを表示装置9に表示させる。また、ドプラ処理部7は、ステップS1で設定されたドプラゲートDG内のドプラデータを取得し、取得されたドプラデータに基づいて、ドプラ波形画像信号を順次、連続的に生成し、生成されたドプラ波形画像信号に基づくドプラ波形画像UDを表示装置9に表示させる。これにより、例えば図11に示すように、Bモード画像UBとドプラ波形画像UDが表示装置9に表示される。また、ステップS2で生成された複数フレームのBモード画像信号は、画像メモリ11に保存される。
続くステップS3において、時相サーチ期間特定部16は、ステップS2で生成されたドプラ波形画像UDに基づいて時相サーチ期間を特定する。この際に、時相サーチ期間特定部16は、図7に示すように、時相サーチ期間として、ドプラ波形画像UDにおける血流速度が最小値V1を有する時点T1から最大値V2を有する時点T2までの第1期間FPと、血流速度が最大値V2を有する時点T2から各心拍期間HCよりも短い時間、例えば0.25秒等に設定された定められた時間が経過する時点T3までの第2期間SPの少なくとも一方を特定する。図7に示すように、第1期間FPにおいて、血管の直径が最小値D1を有し、第2期間SPにおいて、血管の直径が最大値D2を有する。
ステップS4において、フレーム特定部12は、ステップS2で生成された複数フレームのBモード画像信号のうち、ステップS3で特定された時相サーチ期間における複数フレームのBモード画像信号を解析することにより、各心拍期間HCにおいて血管の直径が最大となるフレームのBモード画像信号および血管の直径が最小となるフレームのBモード画像信号の少なくとも一方を特定する。
例えば、ステップS3で時相サーチ期間として第1期間FPが特定された場合に、フレーム特定部12は、第1期間FP内の複数フレームのBモード画像信号において血管壁を検出し、検出された血管壁に基づいて、血管の直径を算出し、算出された血管の直径が最小値D1を有するフレームのBモード画像信号を特定する。
ここで、フレーム特定部12は、各フレームのBモード画像信号において血管の直径を算出する際に、例えば、図9に示すように、入力装置19を介したユーザの入力操作等によりBモード画像UBに設定された、縦方向の直線SL上におけるBモード画像UBの輝度が一定の値よりも高い2点の位置を、血管前壁W1の位置および血管後壁W2の位置として特定し、特定された2点間のBモード画像UB上の距離Lを算出することができる。
また、例えば、ステップS3で時相サーチ期間として第2期間SPが特定された場合に、フレーム特定部12は、第2期間SP内の複数フレームのBモード画像信号において血管壁を検出し、検出された血管壁に基づいて、血管の直径を算出し、算出された血管の直径が最大値D2を有するフレームのBモード画像信号を特定する。
なお、ステップS3で時相サーチ期間として第1期間FPと第2期間SPの双方が特定された場合には、第1期間FP内において血管の直径が最小値D1を有するフレームのBモード画像信号を特定し、且つ、第2期間SP内において血管の直径が最大値D2を有するフレームのBモード画像信号を特定する。
このように、フレーム特定部12は、ステップS3で特定された時相サーチ期間内において血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号の少なくとも一方を特定するため、超音波診断装置1における負担を軽減し、血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号の少なくとも一方を短時間に且つ容易に特定することができる。また、フレーム特定部12は、ノイズの影響等に起因して輝度値に変動が生じたBモード画像信号を誤って特定するおそれを低減することができる。
ステップS5において、血管径算出部13は、ステップS4で特定された、血管の直径が最小となるフレームのBモード画像信号および血管の直径が最大となるフレームのBモード画像信号の少なくとも一方に基づいて、血管の最大径および最小径の少なくとも一方を算出する。例えば、ステップS3で時相サーチ期間として第1期間FPが特定され且つステップS4で血管の直径が最小となるフレームのBモード画像信号が特定された場合に、血管径算出部13は、血管の直径が最小となるフレームのBモード画像信号において、ステップS4で血管の直径として算出されたBモード画像UB上の距離Lを血管の実際の直径に換算する等により、血管の最小径を算出することができる。
また、例えば、ステップS3で時相サーチ期間として第2期間SPが特定され且つステップS4で血管の直径が最大となるフレームのBモード画像信号が特定された場合に、血管径算出部13は、血管の直径が最大となるフレームのBモード画像信号において、ステップS4で血管の直径として算出されたBモード画像UB上の距離Lを血管の実際の直径に換算する等により、血管の最大径を算出することができる。
また、例えば、ステップS3で時相サーチ期間として第1期間FPと第2期間SPの双方が特定され且つステップS4で血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号の双方が特定された場合に、血管径算出部13は、血管の最小径と最大径の双方を算出することができる。
ステップS6において、断面積算出部14は、血管が円形の断面を有するものとして、ステップS5で算出された血管の最小径と最大径の少なくとも一方を用いて、血管の断面積を算出する。例えば、ステップS5で血管の最小径が算出された場合に、断面積算出部14は、血管の最小径を用いて血管の断面積を算出する。また、例えば、ステップS5で血管の最大径が算出された場合に、断面積算出部14は、血管の最大径を用いて血管の断面積を算出する。また、例えば、ステップS5で血管の最小径と最大径の双方が算出された場合には、断面積算出部14は、血管の最小径と最大径に基づいて血管の平均径を算出し、算出された血管の平均径を用いて血管の断面積を算出することができる。
ステップS7において、血流速度算出部17は、ステップS2で生成されたドプラデータに基づいて、血流速度を算出する。
ステップS8において、血流量計測部15は、血管の最小径、最大径または平均径に基づいてステップS6で算出された血管の断面積とステップS7で算出された血流速度とに基づいて、血管内を流れる血液の単位時間当たりの体積を表す血流量を計測し、図12に示すように、血流量の計測値MVを表示装置9に表示する。図9に示す例では、Bモード画像UBおよびドプラ波形画像UDと共に、血流量の計測値MVが表示装置9に表示されている。
このようにして、ステップS8で血流量が計測されると、超音波診断装置1の動作が終了する。
ここで、一連のフレームのBモード画像の中から、例えば、画像が鮮明であるという理由でユーザにより選択されたフレームのBモード画像に基づいて血流量の計測が行われると、血流量を計測する毎に、血管が最小径を有する条件または最大径を有する条件が異なってしまい、同一の被検体に対して、過去に計測された血流量の計測値と新たに計測された血流量の値を比較する際に、正確な比較が困難となってしまうという問題があった。
本発明の実施の形態1に係る超音波診断装置1では、血管の最小径、最大径または平均径に基づいてステップS6で算出された血管の断面積とステップS7で算出された血流速度とに基づいて、ステップS8で血流量が計測されるため、一定の基準に従って、血流量が算出される。そのため、例えば、過去に計測された血流量の計測値と新たに計測された血流量の値を比較する際に、一定の条件に従う正確な比較を行うことが可能である。
また、以上から、本発明の実施の形態1に係る超音波診断装置1によれば、時相サーチ期間特定部16により特定された時相サーチ期間内において血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号の少なくとも一方が特定されるため、定められた期間内にドプラ波形画像信号と共に生成されたすべてのBモード画像信号を用いる場合よりも小さい計算量で血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号の少なくとも一方を特定することができる。そのため、超音波診断装置1における負担を軽減し、血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号の少なくとも一方を短時間に且つ容易に特定することができる。
さらに、超音波診断装置1によれば、血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるBモード画像信号の少なくとも一方を特定する際に、第1期間FP内と第2期間SP内の少なくとも一方のBモード画像信号における血管の直径を互いに比較すればよく、ノイズの影響等に起因して輝度値に変動が生じたBモード画像信号を誤って特定するおそれを低減することができる。そのため、血管径算出部13により、血管の最小径および最大径の少なくとも一方が正確に算出され、血流量計測部15により算出される血流量についても、正確な値が計測されることが可能である。
なお、時相サーチ期間特定部16は、ドプラ波形画像UDにおける血流速度が最小値V1を有する時点T1から最大値V2を有する時点T2までの第1期間FPに対して、第2期間SPを、血流速度が最大値V2を有する時点T2から各心拍期間HCよりも短い時間、例えば0.25秒等に設定された定められた時間が経過する時点T3までの期間として特定するが、第2期間SPの特定方法は、これに限定されない。
例えば、時相サーチ期間特定部16は、第1期間FPと第2期間SPとの和が、心拍期間HCの半分の時間幅を有するように、血流速度が最大値V2を有する時点T2から定められた時間が経過する時点T3までの期間を第2期間SPとして特定することができる。
また、例えば、時相サーチ期間特定部16は、第2期間SPが心拍期間HCの10%以上20%以下の時間幅を有するように、血流速度が最大値V2を有する時点T2から定められた時間が経過する時点T3までの期間を第2期間SPとして特定することもできる。この際に、第1期間FPについても、心拍期間HCの10%以上20%以下の時間幅を有することが好ましい。
また、例えば、時相サーチ期間特定部16は、第2期間SPが第1期間FPと同一の時間幅を有するように、血流速度が最大値V2を有する時点T2から定められた時間が経過する時点T3までの期間を第2期間SPとして特定することもできる。
また、Bモード処理部6とドプラ処理部7により、定められた期間にわたってBモード画像信号とドプラ波形画像信号が連続的に生成されているが、Bモード画像信号とドプラ波形画像信号が生成される代わりに、定められた期間にわたって生成されたBモード画像信号と、同一の定められた期間にわたって生成されたドプラ波形画像信号が予め保存された図示しないメモリが超音波診断装置1に備えられていてもよい。この場合には、時相サーチ期間特定部16により、図示しないメモリに保存されたドプラ波形画像信号に基づいて時相サーチ期間が特定され、フレーム特定部12により、図示しないメモリに保存された複数フレームのBモード画像信号のうち時相サーチ期間内の複数フレームのBモード画像信号に基づいて、血管径が最小となるフレームのBモード画像信号と血管径が最大となるフレームのBモード画像信号の少なくとも一方が特定される。
そのため、Bモード画像信号とドプラ波形画像信号が生成される代わりに、定められた期間にわたって生成されたBモード画像信号と、同一の定められた期間にわたって生成されたドプラ波形画像信号が予め保存された、図示しないメモリが超音波診断装置1に備えられている場合でも、超音波診断装置1は、Bモード画像信号とドプラ波形画像信号が生成される場合と同様に、血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号の少なくとも一方を特定する際の負担を軽減し、血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号の少なくとも一方を短時間に且つ容易に特定することができる。さらに、超音波診断装置1は、ノイズの影響等に起因して輝度値に変動が生じたBモード画像信号を誤って特定するおそれを低減することができる。
また、図10のフローチャートにおけるステップS2で、Bモード処理部6とドプラ処理部7により、定められた期間においてBモード画像信号とドプラ波形画像信号が生成された後、ステップS3において時相サーチ期間が特定されるが、ステップS2の定められた期間内で順次生成されるドプラ波形画像信号に基づいて、リアルタイムに時相サーチ期間を特定する処理が行われてもよい。この場合に、時相サーチ期間特定部16は、例えば、心拍期間HCに満たない期間においてドプラ処理部7がドプラ波形画像信号を生成する毎に、生成されたドプラ波形画像信号に基づいて、時相サーチ期間を特定する処理を行う。この場合には、Bモード処理部6とドプラ処理部7により、第1期間FP以降または第2期間SP以降においてBモード画像信号とドプラ波形画像信号が取得される必要がなくなるため、より短時間に且つ容易に、血管の直径が最小となるBモード画像信号および血管の直径が最大となるBモード画像信号の少なくとも一方が特定される。
また、フレーム特定部12は、Bモード処理部6により生成されたBモード画像信号に基づいて、血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号を特定する代わりに、受信回路4により生成されたRFデータに基づいて、血管の直径が最小となるフレームのRFデータと血管の直径が最大となるフレームのRFデータを特定することもできる。この場合でも、フレーム特定部12により、Bモード画像信号に基づいて、血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号が特定される場合と同様に、超音波診断装置1における負担を軽減して、血管の直径が最小となるフレームのRFデータと血管の直径が最大となるフレームのRFデータが容易に且つ正確に特定されることができる。
実施の形態2
実施の形態1では、フレーム特定部12により、血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号の双方が特定され、血管径算出部13により、各心拍期間HCにおける血管の最小径と最大径の双方が算出されることができるが、例えば、このようにして算出された血管の最小径と最大径を用いて、血管の弾性を表す弾性指標を算出することができる。
図13に示すように、本発明の実施の形態2に係る超音波診断装置1Aは、図1に示す実施の形態1の超音波診断装置1において、装置制御部18の代わりに装置制御部18Aを備え、弾性指標算出部41が追加されたものである。
超音波診断装置1Aにおいて、血管径算出部13に弾性指標算出部41が接続され、弾性指標算出部41に表示制御部8および装置制御部18Aが接続されている。また、Bモード処理部6、ドプラ処理部7、表示制御部8、ゲート設定部10、フレーム特定部12、血管径算出部13、断面積算出部14、血流量計測部15、時相サーチ期間特定部16、血流速度算出部17、装置制御部18Aおよび弾性指標算出部41により、超音波診断装置1A用のプロセッサ22Aが構成されている。
弾性指標算出部41は、血管径算出部13により算出された血管の最大径と最小径との差に基づいて、血管の弾性指標を算出する。血管の弾性指標とは、血管の弾性を表す指標である。血管の最大径と最小径との差が大きいほど、血管の直径の変化が大きく、血管の弾性率が低いと判断でき、血管の最大径と最小径との差が小さいほど、血管の直径の変化が小さく、血管の弾性率が高いと判断できる。そのため、弾性指標算出部41は、例えば、血管の最大径と最小径との差を、血管の弾性指標として算出することができる。また、弾性指標算出部41は、例えば、血管の最大径と最小径との差を血管の最小径で除すことにより規格化したものを弾性指標として算出することもできる。
また、図示しない血圧計を用いて、血管の直径が最小となる時点における被検体の血圧P1と、血管の直径が最大となる時点における被検体の血圧P2とを計測することにより、弾性指標算出部41は、血圧P1、P2、心拍期間HCにおける血管の直径の最小値D1、心拍期間HCにおける血管の直径の最大値D2を用いて、特許第5384919号公報に記載されるスティフネスパラメータB={Log(P2/P1)}/{(D2/D1)-1}を弾性指標として算出することもできる。
以上から、本発明の実施の形態2に係る超音波診断装置1Aによれば、実施の形態1に係る超音波診断装置1と同様に、血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるBモード画像信号を特定する際に、超音波診断装置1Aにおける負担を軽減し、血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるBモード画像信号を短時間に且つ容易に特定し、さらに、ノイズの影響等に起因して輝度値に変動が生じたBモード画像信号を誤って特定するおそれを低減することができるため、弾性指標についても、短時間に且つ容易に且つ正確に算出されることができる。
実施の形態3
実施の形態1の超音波診断装置1は、表示装置9、入力装置19、超音波プローブ21がプロセッサ22に直接的に接続される構成を有しているが、例えば、表示装置9、入力装置19、超音波プローブ21、プロセッサ22がネットワークを介して間接的に接続されることもできる。
図14に示すように、実施の形態3における超音波診断装置1Bは、表示装置9、入力装置19、超音波プローブ21がネットワークNWを介して超音波診断装置本体51に接続されたものである。超音波診断装置本体51は、図1に示す実施の形態1の超音波診断装置1において、表示装置9、入力装置19、超音波プローブ21を除いたものであり、送受信回路5、格納部20およびプロセッサ22により構成されている。
超音波診断装置1Bがこのような構成を有している場合でも、実施の形態1の超音波診断装置1と同様に、時相サーチ期間特定部16により特定された時相サーチ期間内において血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号の少なくとも一方が特定されるため、定められた期間内にドプラ波形画像信号と共に生成されたすべてのBモード画像信号を用いる場合よりも小さい計算量で血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号の少なくとも一方を特定することができる。そのため、超音波診断装置1における負担を軽減し、血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるフレームのBモード画像信号の少なくとも一方を短時間に且つ容易に特定することができる。
さらに、超音波診断装置1Bによれば、実施の形態1の超音波診断装置1と同様に、血管の直径が最小となるフレームのBモード画像信号と血管の直径が最大となるBモード画像信号の少なくとも一方を特定する際に、第1期間FP内と第2期間SP内の少なくとも一方のBモード画像信号における血管の直径を互いに比較すればよく、ノイズの影響等に起因して輝度値に変動が生じたBモード画像信号を誤って特定するおそれを低減することができる。
また、表示装置9、入力装置19、超音波プローブ21がネットワークNWを介して超音波診断装置本体51と接続されているため、超音波診断装置本体51を、いわゆる遠隔サーバとして使用することができる。これにより、例えば、ユーザは、表示装置9、入力装置19、超音波プローブ21をユーザの手元に用意することにより、被検体の診断を行うことができるため、超音波診断の際の利便性を向上させることができる。
また、例えば、いわゆるタブレットと呼ばれる携帯型の薄型コンピュータが表示装置9および入力装置19として使用される場合には、ユーザは、より容易に被検体の超音波診断を行うことができ、超音波診断の利便性をさらに向上させることができる。
なお、表示装置9、入力装置19、超音波プローブ21がネットワークNWを介して超音波診断装置本体51に接続されているが、この際に、表示装置9、入力装置19、超音波プローブ21は、ネットワークNWに有線接続されていてもよく、無線接続されていてもよい。
また、実施の形態3の態様は、実施の形態1に適用されることが説明されているが、実施の形態2についても、同様に適用されることができる。
1,1A,1B 超音波診断装置、2 振動子アレイ、3 送信回路、4 受信回路、5 送受信回路、6 Bモード処理部、7 ドプラ処理部、8 表示制御部、9 表示装置、10 ゲート設定部、11 画像メモリ、12 フレーム特定部、13 血管径算出部、14 断面積算出部、15 血流量計測部、16 時相サーチ期間特定部、17 血流速度算出部、18,18A 装置制御部、19 入力装置、20 格納部、21 超音波プローブ、22,22A プロセッサ、23 増幅部、24 AD変換部、25 ビームフォーマ、26 信号処理部、27 DSC、28 画像処理部、29 直交検波部、30 ハイパスフィルタ、31 高速フーリエ変換部、32 ドプラ波形画像生成部、33 データメモリ、41 弾性指標算出部、51 超音波診断装置本体、52 送受信回路、BR 血管領域、D1,V1 最小値、D2,V2 最大値、DG ドプラゲート、E 包絡線、FP 第1期間、HC 心拍期間、L 距離、MV 計測値、SL 直線、SP 第2期間、T1,T2,T3,T4 時点、UB Bモード画像、UD ドプラ波形画像、W1 血管前壁、W2 血管後壁、WD ドプラ波形。

Claims (13)

  1. 被検体の血管が含まれる領域に対して定められた期間にわたってBモードデータとドプラデータとが連続的に取得される超音波診断装置であって、
    前記ドプラデータに基づいて各心拍期間における時相サーチ期間を特定する時相サーチ期間特定部と、
    前記時相サーチ期間特定部により特定された前記時相サーチ期間内における複数フレームの前記Bモードデータを解析することにより、各心拍期間において前記血管の直径が最大となるフレームのBモードデータおよび前記血管の直径が最小となるフレームのBモードデータの少なくとも一方を特定するフレーム特定部と
    を備える超音波診断装置。
  2. 前記フレーム特定部により特定された、前記血管の直径が最大となるフレームのBモードデータおよび前記血管の直径が最小となるフレームのBモードデータの少なくとも一方に基づいて前記血管の最大径および最小径の少なくとも一方を算出する血管径算出部をさらに備える請求項1に記載の超音波診断装置。
  3. 前記血管径算出部により算出された前記血管の最大径および最小径の少なくとも一方を用いて前記血管の断面積を算出する断面積算出部と、
    各心拍期間における前記ドプラデータを取得するドプラ処理部と、
    前記ドプラ処理部により取得された前記ドプラデータに基づいて血流速度を算出する血流速度算出部と、
    前記断面積算出部により算出された前記断面積と前記血流速度算出部により算出された前記血流速度に基づいて血流量を計測する血流量計測部と
    をさらに備える請求項2に記載の超音波診断装置。
  4. 前記時相サーチ期間特定部は、各心拍期間におけるドプラデータが最小値を有する時点を含む第1期間を前記時相サーチ期間として特定し、
    前記フレーム特定部は、前記第1期間内における複数フレームの前記Bモードデータに基づいて前記血管の直径が最小となるフレームのBモードデータを特定し、
    前記血管径算出部は、前記フレーム特定部により特定された、前記血管の直径が最小となるフレームのBモードデータに基づいて、前記血管の最小径を算出し、
    前記断面積算出部は、前記血管径算出部により算出された前記血管の最小径を用いて前記血管の断面積を算出する請求項3に記載の超音波診断装置。
  5. 前記時相サーチ期間特定部は、各心拍期間におけるドプラデータが最大値を有する時点を含む第2期間を前記時相サーチ期間として特定し、
    前記フレーム特定部は、前記第2期間内における複数フレームの前記Bモードデータに基づいて前記血管の直径が最大となるフレームのBモードデータを特定し、
    前記血管径算出部は、前記フレーム特定部により特定された、前記血管の直径が最大となるフレームのBモードデータに基づいて、前記血管の最大径を算出し、
    前記断面積算出部は、前記血管径算出部により算出された前記血管の最大径を用いて前記血管の断面積を算出する請求項3に記載の超音波診断装置。
  6. 前記時相サーチ期間特定部は、各心拍期間におけるドプラデータが最小値を有する時点を含む第1期間および各心拍期間におけるドプラデータが最大値を有する時点を含む第2期間を有する前記時相サーチ期間を特定し、
    前記フレーム特定部は、前記第1期間内における複数フレームの前記Bモードデータに基づいて前記血管の直径が最小となるフレームのBモードデータを特定し、且つ、前記第2期間内における複数フレームの前記Bモードデータに基づいて前記血管の直径が最大となるフレームのBモードデータを特定し、
    前記血管径算出部は、前記フレーム特定部により特定された、前記血管の直径が最小となるフレームのBモードデータに基づいて、前記血管の最小径を算出し、且つ、前記フレーム特定部により特定された、前記血管の直径が最大となるフレームのBモードデータに基づいて、前記血管の最大径を算出し、
    前記断面積算出部は、前記血管径算出部により算出された前記血管の最小径および前記血管の最大径とから算出された各心拍期間における前記血管の平均径を用いて前記血管の断面積を算出する請求項3に記載の超音波診断装置。
  7. 前記第1期間は、各心拍期間におけるドプラデータが最小値を有する時点から最大値を有する時点までの期間である請求項4または6に記載の超音波診断装置。
  8. 前記第2期間は、各心拍期間におけるドプラデータが最大値を有する期間から各心拍期間よりも短い時間に設定された定められた時間が経過するまでの期間である請求項5または6に記載の超音波診断装置。
  9. 前記第1期間および前記第2期間は、それぞれ、各心拍期間の10%以上20%以下の時間幅を有する請求項4~6のいずれか一項に記載の超音波診断装置。
  10. 前記血管の最大径と前記血管の最小径との差に基づいて前記血管の弾性指標を算出する弾性指標算出部をさらに備える請求項6に記載の超音波診断装置。
  11. 超音波プローブと、
    前記超音波プローブを介して前記被検体内に超音波ビームを送信し且つ前記被検体内から超音波エコーを受信して受信データを生成する送受信回路と、
    前記送受信回路により生成された前記受信データに基づいて前記Bモードデータを生成するBモード処理部と
    をさらに備える請求項1~10のいずれか一項に記載の超音波診断装置。
  12. 被検体の血管が含まれる領域に対して定められた期間にわたってBモードデータとドプラデータとが連続的に取得される超音波診断装置の制御方法であって、
    前記ドプラデータに基づいて各心拍期間における時相サーチ期間を特定し、
    特定された前記時相サーチ期間内における複数フレームの前記Bモードデータを解析することにより、各心拍期間において前記血管の直径が最大となるフレームのBモードデータおよび前記血管の直径が最小となるフレームのBモードデータの少なくとも一方を特定する
    超音波診断装置の制御方法。
  13. 被検体の血管が含まれる領域に対して定められた期間にわたってBモードデータとドプラデータとを連続的に取得し、前記ドプラデータに基づいて各心拍期間における時相サーチ期間を特定し、特定された前記時相サーチ期間内における複数フレームの前記Bモードデータを解析することにより、各心拍期間において前記血管の直径が最大となるフレームのBモードデータおよび前記血管の直径が最小となるフレームのBモードデータの少なくとも一方を特定する超音波診断装置用プロセッサ。
JP2021514934A 2019-04-17 2020-04-13 超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ Active JP7095177B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2019078241 2019-04-17
JP2019078241 2019-04-17
PCT/JP2020/016285 WO2020213562A1 (ja) 2019-04-17 2020-04-13 超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ

Publications (2)

Publication Number Publication Date
JPWO2020213562A1 JPWO2020213562A1 (ja) 2020-10-22
JP7095177B2 true JP7095177B2 (ja) 2022-07-04

Family

ID=72837458

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021514934A Active JP7095177B2 (ja) 2019-04-17 2020-04-13 超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ

Country Status (3)

Country Link
US (1) US20220022849A1 (ja)
JP (1) JP7095177B2 (ja)
WO (1) WO2020213562A1 (ja)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001218768A (ja) 2000-02-10 2001-08-14 Aloka Co Ltd 超音波診断装置
WO2006011504A1 (ja) 2004-07-28 2006-02-02 Matsushita Electric Industrial Co., Ltd. 超音波診断装置および超音波診断装置の制御方法
JP2006102489A (ja) 2004-09-07 2006-04-20 Toshiba Corp 超音波ドプラ診断装置及び診断パラメータ計測方法
JP2016027835A (ja) 2012-12-11 2016-02-25 日立アロカメディカル株式会社 超音波撮像装置及び方法
JP2017524455A (ja) 2015-04-03 2017-08-31 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 血管を識別する超音波システム及び方法
US20190099153A1 (en) 2016-03-22 2019-04-04 Imperial Innovations Limited Fluid flow analysis

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001218768A (ja) 2000-02-10 2001-08-14 Aloka Co Ltd 超音波診断装置
WO2006011504A1 (ja) 2004-07-28 2006-02-02 Matsushita Electric Industrial Co., Ltd. 超音波診断装置および超音波診断装置の制御方法
JP2006102489A (ja) 2004-09-07 2006-04-20 Toshiba Corp 超音波ドプラ診断装置及び診断パラメータ計測方法
JP2016027835A (ja) 2012-12-11 2016-02-25 日立アロカメディカル株式会社 超音波撮像装置及び方法
JP2017524455A (ja) 2015-04-03 2017-08-31 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 血管を識別する超音波システム及び方法
US20190099153A1 (en) 2016-03-22 2019-04-04 Imperial Innovations Limited Fluid flow analysis

Also Published As

Publication number Publication date
JPWO2020213562A1 (ja) 2020-10-22
WO2020213562A1 (ja) 2020-10-22
US20220022849A1 (en) 2022-01-27

Similar Documents

Publication Publication Date Title
JP5230106B2 (ja) 超音波診断装置、imt計測方法及びimt計測プログラム
JP7476376B2 (ja) 超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ
JP6484389B2 (ja) 超音波診断装置および超音波診断装置の制御方法
US11039777B2 (en) Ultrasonic diagnostic apparatus and control method
JPWO2017038162A1 (ja) 超音波診断装置および超音波診断装置の制御方法
JP5588924B2 (ja) 超音波診断装置
JP5663552B2 (ja) 超音波検査装置、超音波検査装置の信号処理方法およびプログラム
US20210267569A1 (en) Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus
WO2021014926A1 (ja) 超音波診断装置および超音波診断装置の制御方法
WO2021161683A1 (ja) 超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ
JPWO2019058753A1 (ja) 超音波診断装置および超音波診断装置の制御方法
JP5844175B2 (ja) 超音波診断装置および超音波画像生成方法
JP6856816B2 (ja) 超音波診断装置および超音波診断装置の制御方法
JP7095177B2 (ja) 超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ
JP2013244164A (ja) 超音波画像生成方法および超音波画像診断装置
JP5869411B2 (ja) 超音波診断装置および超音波画像生成方法
US20200359989A1 (en) Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus
JP5851345B2 (ja) 超音波診断装置およびデータ処理方法
JP2010124852A (ja) 超音波診断装置
WO2017122411A1 (ja) 超音波診断装置および音速定量化方法
JP2013244162A (ja) 超音波診断装置
JP2006149701A (ja) 断層画像計測装置および超音波診断装置
WO2017126209A1 (ja) 超音波診断装置および音速定量化方法
WO2017122414A1 (ja) 超音波診断装置および音速定量化方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20211011

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220622

R150 Certificate of patent or registration of utility model

Ref document number: 7095177

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150