JP4441664B2 - 血管形状測定装置および血流量測定装置 - Google Patents

血管形状測定装置および血流量測定装置 Download PDF

Info

Publication number
JP4441664B2
JP4441664B2 JP2004219866A JP2004219866A JP4441664B2 JP 4441664 B2 JP4441664 B2 JP 4441664B2 JP 2004219866 A JP2004219866 A JP 2004219866A JP 2004219866 A JP2004219866 A JP 2004219866A JP 4441664 B2 JP4441664 B2 JP 4441664B2
Authority
JP
Japan
Prior art keywords
blood vessel
section
array
cross
vessel wall
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.)
Expired - Fee Related
Application number
JP2004219866A
Other languages
English (en)
Other versions
JP2006034667A (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.)
UNEX CORPORATION
University of Tokushima
Original Assignee
UNEX CORPORATION
University of Tokushima
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 UNEX CORPORATION, University of Tokushima filed Critical UNEX CORPORATION
Priority to JP2004219866A priority Critical patent/JP4441664B2/ja
Priority to PCT/JP2005/013813 priority patent/WO2006011544A1/ja
Priority to US10/565,675 priority patent/US7429244B2/en
Priority to EP05767223A priority patent/EP1772102A4/en
Publication of JP2006034667A publication Critical patent/JP2006034667A/ja
Application granted granted Critical
Publication of JP4441664B2 publication Critical patent/JP4441664B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4887Locating particular structures in or on the body
    • A61B5/489Blood vessels
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/02007Evaluating blood vessel condition, e.g. elasticity, compliance
    • 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/13Tomography
    • A61B8/14Echo-tomography
    • 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/523Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for generating planar views from image data in a user selectable plane not corresponding to the acquisition plane

Landscapes

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

Description

本発明は、生体の表示下に位置する血管の形状を超音波を用いて測定する装置、その血管内の血流速度を測定する血流速度測定装置、およびそれを用いてその血管内の血流量を測定する血流量測定装置に関するものである。
生体の状態を把握するために、非侵襲でその生体の血管の径や断面積、血流速度、血流量を正確に測定することが要請されている。たとえば、血管の内皮機能検査に際しては、内径4mmφの動脈において、1%の精度で血管径や断面積、血流速度、血流量を測定することが求められる。
これに対し、特許文献1には、非侵襲で生体表皮下の血管の長手方向の2位置に、一方向に配列された複数の第1超音波素子を有する第1アレイと、該第1超音波素子の配列方向と平行な方向に配列された複数の第2超音波素子を有する第2アレイとを配置し、それら第1アレイおよび第2アレイの直下の2位置において血管の断面形状を検出する技術が提案されている。これによれば、超音波を用いているため、非侵襲にて血管の断面形状が検出される。
特開平11−76233号公報
しかしながら、上記従来の血管の断面形状を検出する装置は、第1アレイの複数の第1超音波素子の配列方向を含む測定断面内および第2アレイの複数の第1超音波素子の配列方向を含む測定断面内における血管の断面形状を検出するものであることから、必ずしも血管の長手方向とその測定断面とは直交するものではないので、十分な精度が得られないという問題があった。生体の皮膚下の血管上に上記第1アレイおよび第2アレイが固定された超音波プローブが装着されるとき、動脈の位置は目視し難いことから第1アレイおよび第2アレイがその動脈に直交するように装着したとしても必ずしも正確ではなく、しかも、動脈は皮膚に対して平行でない場合が殆どであるので、動脈血管の長手方向と上記第1アレイや第2アレイの測定断面とは直交しないので、その測定断面内で測定される血管の径や断面積等の断面形状、それに基づく血流速度や血流量を正確に測定することができなかった。
本発明は以上の事情を背景として為されたもので、その目的とするところは、高精度で血管の形状、その血管内の血流速度や血流量を測定することができる血管形状測定装置および血流量測定装置を提供することにある。
上記目的を達成するための請求項1に係る発明の血管形状測定装置の要旨とするところは、一方向に配列された複数の第1超音波素子を有する第1アレイと、その第1超音波素子の配列方向と平行な方向に配列された複数の第2超音波素子を有する第2アレイとを備え、生体の表皮下に位置する血管に交差してその表皮上に配置されたその第1アレイおよび第2アレイから検出されるエコー信号に基づいてその血管の形状を測定する血管形状測定装置であって、(a) 前記第1アレイの各第1超音波素子から放射される放射信号とその各第1超音波素子により検出される前記血管壁からの反射信号との時間差に基づいてその血管壁までの距離をそれぞれ算出し、その距離に基づいて前記第1アレイの測定断面内における血管壁の位置を決定する第1血管壁位置演算手段と、(b) 前記第2アレイの各第2超音波素子から放射される放射信号と各第2超音波素子により検出される前記血管壁からの反射信号との時間差に基づいてその血管壁までの距離をそれぞれ算出し、その距離に基づいて前記第2アレイの測定断面内における血管壁の位置を決定する第2血管壁位置演算手段と、(c) 前記第1血管壁位置演算手段により決定された各第1超音波素子に対応する血管壁部位の位置と前記第2血管壁位置演算手段により決定された各第2超音波素子に対応する血管壁部位の位置とに基づいて、前記血管の直交断面における形状を算出する血管形状算出手段とを、含み、その血管形状算出手段は、(d) 前記第1血管壁位置演算手段により決定された各第1超音波素子に対応する血管壁部位の位置と前記第2血管壁位置演算手段により決定された各第2超音波素子に対応する血管壁部位の位置とに基づいて、前記第1アレイの測定断面内と前記第2アレイの測定断面内とにおける血管壁の断面の中心点と長軸径および/または短軸径とをそれぞれ算出する測定断面内形状算出手段と、(e) その測定断面内形状算出手段により算出された第1アレイの測定断面と前記第2アレイの測定断面とにおける血管壁の断面の中心点に基づいてその血管の中心軸を算出する中心軸算出手段と、(f) その中心軸算出手段により算出された血管の中心軸に基づいて、その血管の直交断面と前記測定断面との交差角度を算出する交差角度算出手段と、(g) その交差角度算出手段により算出された交差角度に基づいて、前記血管の直交断面の長軸径および/または短軸径となるように、前記形状算出手段により算出された長軸径および/または短軸径を補正する補正手段とを、含むことにある。
また、請求項にかかる発明では、上記請求項に係る発明において、前記血管形状算出手段は、前記補正手段によって補正された長軸径および短軸径に基づいて前記血管の直交断面における断面積を算出する直交断面積算出手段を、さらに含むものである。
前記目的を達成するための請求項に係る発明の血流量測定装置の要旨とするところは、(a) 一方向に配列された複数の第1超音波素子を有する第1アレイと、その第1超音波素子の配列方向と平行な方向に配列された複数の第2超音波素子を有する第2アレイと、ドップラ用超音波素子とを備え、その第1アレイおよび第2アレイが前記生体の表皮下に位置する血管と交差し且つそのドップラ用超音波素子からの超音波を放射方向がその血管に対して鋭角を成すように装着される超音波プローブと、(b) 前記ドップラ用超音波素子から前記血管に向かって超音波を放射したときに得られる、その血管内の血流速度に基づくドップラ効果により変化させられるドップラ反射波に基づいてその血流速度を算出する血流速度算出手段と、(c) 前記第1アレイの各第1超音波素子から検出される反射信号に基づいて、前記血管のその第1アレイの直下に位置する部分の血管壁のうちその各第1超音波素子に対応する部位の位置をそれぞれ算出する第1血管壁位置演算手段と、(d) 前記第2アレイの各第2超音波素子から検出されるエコー信号に基づいて、前記血管のその第2アレイの直下に位置する部分の血管壁のうちその各第2超音波素子に対応する部位の位置をそれぞれ算出する第2血管壁位置演算手段と、(e) 前記第1血管壁位置演算手段により算出された各第1超音波素子に対応する血管壁部位の位置と前記第2血管壁位置演算手段により算出された各第2超音波素子に対応する血管壁部位の位置とに基づいて、前記第1アレイの測定断面内と前記第2アレイの測定断面内とにおける血管壁の断面の中心点を算出し、その第1アレイの測定断面と前記第2アレイの測定断面とにおける血管壁の断面の中心点に基づいてその血管の中心軸を算出する中心軸算出手段と、(f) その中心軸算出手段により算出された前記血管の中心軸と前記ドップラ用超音波素子からその血管に向かって超音波を放射する方向との実際の相対角度を算出する相対角度算出手段と、(g) その相対角度算出手段により算出された実際の相対角度に基づいて前記血流速度算出手段により算出された血流速度を補正する血流速度補正手段とを、含む血流速度測定装置を備え、(h) 前記第1血管壁位置演算手段により算出された各第1超音波素子に対応する血管壁部位の位置と前記第2血管壁位置演算手段により算出された各第2超音波素子に対応する血管壁部位の位置とに基づいて、前記第1アレイの測定断面内と前記第2アレイの測定断面内とにおける血管壁の断面の中心点と長軸径および/または短軸径とをそれぞれ算出する測定断面内形状算出手段と、(i) 前記中心軸算出手段により算出された血管の中心軸に基づいて、その血管の直交断面と前記測定断面との交差角度を算出する交差角度算出手段と、(j) その交差角度算出手段により算出された交差角度に基づいて、前記血管の直交断面の長軸径および/または短軸径となるように、前記形状算出手段により算出された長軸径および/または短軸径を補正する補正手段と、(k) その補正手段によって補正された長軸径および短軸径に基づいて前記血管の直交断面における断面積を算出する直交断面積算出手段と、(l) その直交断面積算出手段により算出された前記血管の直交断面における断面積と前記血流速度補正手段により補正された血流速度とに基づいて、前記血管内の血流量を算出する血流量算出手段とを、含むことにある。
前記請求項1に係る発明によれば、生体の表皮下に位置する血管に交差してその表皮上に配置されたその第1アレイおよび第2アレイから検出されるエコー信号に基づいて、その第1アレイおよび第2アレイの直下に位置する部分の血管壁のうち各第1超音波素子および第2超音波素子に対応する部位の位置がそれぞれ算出され、それら血管壁の第1超音波素子および第2超音波素子に対応する部位の位置に基づいて、血管の直交断面における形状が算出されるので、生体皮膚下の血管の方向が上記第1アレイおよび第2アレイと直交していなくても、或いは皮膚に平行でなくても、その血管の血管径、血管断面形状、或いは断面積等の断面形状が正確に得られる。また、第1血管壁位置演算手段および第2血管壁位置演算手段によって、各第1超音波素子および各第2超音波素子から放射される放射信号とその各第1超音波素子および各第2超音波素子により検出される前記血管壁からの反射信号との時間差に基づいてその血管壁までの距離がそれぞれ算出され、その距離に基づいて前記第1アレイの測定断面内における血管壁の位置が決定されるので、第1アレイおよび第2アレイの測定断面内において各血管壁の位置が正確に得られる。また、前記血管形状算出手段において、第1アレイおよび第2アレイの測定断面内の血管壁の部位の位置に基づいて、その第1アレイの測定断面内と前記第2アレイの測定断面内とにおける血管壁の断面の中心点と長軸径および/または短軸径とがそれぞれ算出され、その第1アレイの測定断面内と前記第2アレイの測定断面内とにおける血管壁の断面の中心点に基づいて血管の中心軸が算出され、その血管の中心軸に基づいて前記血管の直交断面と前記測定断面との交差角度が算出され、その交差角度に基づいて前記血管の直交断面の長軸径および/または短軸径となるように、長軸径および/または短軸径が補正されるので、生体皮膚下の血管の方向が上記第1アレイおよび第2アレイと直交していなくても、或いは皮膚に平行でなくても、その血管の断面形状が正確に得られる。
また、請求項に係る発明によれば、前記血管形状算出手段は、前記補正手段によって補正された長軸径および短軸径に基づいて前記血管の直交断面における断面積を算出する直交断面積算出手段を含むものであるので、生体皮膚下の血管の方向が上記第1アレイおよび第2アレイと直交していなくても、或いは皮膚に平行でなくても、その血管の断面積が正確に得られる。
また、請求項に係る発明の血流量測定装置によれば、(a) 一方向に配列された複数の第1超音波素子を有する第1アレイと、その第1超音波素子の配列方向と平行な方向に配列された複数の第2超音波素子を有する第2アレイと、ドップラ用超音波素子とを備え、その第1アレイおよび第2アレイが前記生体の表皮下に位置する血管と交差し且つそのドップラ用超音波素子からの超音波を放射方向がその血管に対して鋭角を成すように装着される超音波プローブと、(b) 前記ドップラ用超音波素子から前記血管に向かって超音波を放射したときに得られる、その血管内の血流速度に基づくドップラ効果により変化させられるドップラ反射波に基づいてその血流速度を算出する血流速度算出手段と、(c) 前記第1アレイの各第1超音波素子から検出される反射信号に基づいて、前記血管のその第1アレイの直下に位置する部分の血管壁のうちその各第1超音波素子に対応する部位の位置をそれぞれ算出する第1血管壁位置演算手段と、(d) 前記第2アレイの各第2超音波素子から検出されるエコー信号に基づいて、前記血管のその第2アレイの直下に位置する部分の血管壁のうちその各第2超音波素子に対応する部位の位置をそれぞれ算出する第2血管壁位置演算手段と、(e) 前記第1血管壁位置演算手段により算出された各第1超音波素子に対応する血管壁部位の位置と前記第2血管壁位置演算手段により算出された各第2超音波素子に対応する血管壁部位の位置とに基づいて、前記第1アレイの測定断面内と前記第2アレイの測定断面内とにおける血管壁の断面の中心点を算出し、その第1アレイの測定断面と前記第2アレイの測定断面とにおける血管壁の断面の中心点に基づいてその血管の中心軸を算出する中心軸算出手段と、(f) その中心軸算出手段により算出された前記血管の中心軸と前記ドップラ用超音波素子からその血管に向かって超音波を放射する方向との実際の相対角度を算出する相対角度算出手段と、(g) その相対角度算出手段により算出された実際の相対角度に基づいて前記血流速度算出手段により算出された血流速度を補正する血流速度補正手段とを、含む血流速度測定装置を備え、(h) 前記第1血管壁位置演算手段により算出された各第1超音波素子に対応する血管壁部位の位置と前記第2血管壁位置演算手段により算出された各第2超音波素子に対応する血管壁部位の位置とに基づいて、前記第1アレイの測定断面内と前記第2アレイの測定断面内とにおける血管壁の断面の中心点と長軸径および/または短軸径とをそれぞれ算出する測定断面内形状算出手段と、(i) 前記中心軸算出手段により算出された血管の中心軸に基づいて、その血管の直交断面と前記測定断面との交差角度を算出する交差角度算出手段と、(j) その交差角度算出手段により算出された交差角度に基づいて、前記血管の直交断面の長軸径および/または短軸径となるように、前記形状算出手段により算出された長軸径および/または短軸径を補正する補正手段と、(k) その補正手段によって補正された長軸径および短軸径に基づいて前記血管の直交断面における断面積を算出する直交断面積算出手段と、(l) その直交断面積算出手段により算出された前記血管の直交断面における断面積と前記血流速度補正手段により補正された血流速度とに基づいて、前記血管内の血流量を算出する血流量算出手段とが、設けられているので、生体皮膚下の血管の方向が上記第1アレイおよび第2アレイと直交していなくても、或いは皮膚に平行でなくても、その血管内の血流量が正確に得られる。
ここで、好適には、前記血管は、生体の皮膚下に位置する動脈たとえば上腕動脈、撓骨動脈、足背動脈、頸動脈、浅側頭動脈等である。FMD(flow-medeated dilation: 内皮依存性血管拡張反応)検査の場合には、たとえば上腕動脈、撓骨動脈、足背動脈が対象となり、頭部への血流把握の場合には、たとえば頸動脈、浅側頭動脈が対象となる。
また、好適には、前記第1アレイおよび第2アレイを構成する複数の超音波素子は、振動子と受信子とが共通の振動子から構成されたものであってもよいし、一対の発信子と受信子とからそれぞれ構成されていてもよい。
また、好適には、ドップラ用超音波素子は、送信子と受信子とが共通の1個の振動子であってもよいし、一対の発信子と受信子とから構成されていてもよい。また、このドップラ用超音波素子は、前記血管の方向に対して交差する方向に配列された複数個の振動子アレイから構成されてもよい。
以下、本発明の好適な実施の形態について図面を参照しつつ詳細に説明する。図1は、本発明の一実施例の超音波測定装置10に備えられた超音波プローブ12の生体14に対する装着状態を示している。この超音波プローブ12は、図2或いは図3に示すように、装着バンド16によって生体の一部たとえば上腕部において上腕動脈18の真上に位置するように装着される。この上腕動脈18は、上腕二頭筋の下端部下側から表皮20に向かって接近する形状を有している。
上記超音波プローブ12は、その表皮20と接触する接触面22を有し、アクリル樹脂等の合成樹脂、セラミックス、金属などの超音波透過性材料から構成された本体24と、その本体24内に所定の間隔Lを隔てて互いに平行となるように埋設され且つ上記上腕動脈18と交差する長手方向の第1アレイ26および第2アレイ28と、その第2アレイ28から下流側に位置するように配設され且つ上記上腕動脈18と交差する長手方向の第3アレイ30とを備えている。第1アレイ26および第2アレイ28は、上記一方向すなわち上記長手方向に配列された複数個の超音波素子26n および28n (nは整数)をそれぞれ備えている。上記第3アレイ30も、上記の第1アレイ26および第2アレイ28と平行な一方向すなわち上記長手方向に配列された複数個の超音波素子30n を備えている。上記超音波素子26n および28n は、測定対象となる上腕動脈18の径よりも十分に小さい間隔たとえば0.2乃至0.5mm程度となるように配置されている。超音波素子30n もそれと同等であるかそれよりも大きい間隔で配置されている。
上記第1アレイ26の超音波素子26n および第2アレイ28の超音波素子28n は、圧電素子などから構成されることによって超音波の発信子および受信子として機能するものであり、真下に向かって超音波を発信するように接触面22に対して垂直に設けられている。上記第3アレイ30の超音波素子30n も、圧電素子などから構成されることによって超音波の発信子および受信子として機能するものであり、上流部位の上腕動脈18に向かって超音波を発信するように接触面22に対して超音波放射面が斜めにたとえば45度程度傾斜して設けられている。上記第1アレイ26の超音波素子26n および第2アレイ28の超音波素子28n は、X−Z平面として設定される接触面22内に略位置するように配置され、第1アレイ26の超音波素子26n のうちの端部に位置する素子261 が上記X−Y−Z三次元直交座標の原点として設定されている。
第1アレイ26の複数個の超音波素子26n および第2アレイ28の複数個の超音波素子28n は、それぞれ、上記接触面22に対して垂直方向に超音波を放射し、その超音波の伝播過程に存在する界面からの反射波をエコー信号として受信するので、図2に示すように、第1アレイ26および第2アレイ28の直下には、互いに平行であって上記上腕動脈18がそれぞれ貫通させられる測定断面Aおよび測定断面Bが形成されるようになっている。それらの測定断面Aおよび測定断面Bは前記X−Y−Z三次元直交座標のX−Y平面に平行な面となり、それら測定断面Aおよび測定断面Bの面内の位置は座標値により表される。
図3に示すように、上記超音波測定装置10は、第1アレイ26、第2アレイ28、第3アレイ30を駆動制御するための超音波駆動制御回路32と、アナログ信号およびデジタル信号の一方から他方へ変換するための信号変換器34と、電子制御装置36と、数字および画像を表示するための表示器38とを備えている。上記電子制御装置36はCPU40、ROM42、RAM44、図示しない記憶装置やインターフェース等から成る所謂マイクロコンピュータから構成されており、CPU40はRAM44の一時記憶機能を利用しつつ予めROM42に記憶されたプログラムにしたがって入力信号を処理して、血管形状測定装置、血流速度測定装置、および血流量を算出し、演算結果を表示器38に表示させ、或いは他の機器へ出力する。このため、超音波測定装置10は、血管形状測定装置、血流速度測定装置、および血流量測定装置として機能している。
図4は、上記電子制御装置36の演算制御機能の要部を説明する機能ブロック線図である。図4において、第1管壁位置演算手段48および第2管壁位置演算手段50は、第1アレイ26の直下の測定断面Aおよび第2アレイ28の直下の測定断面Bにおける上腕動脈18の断面形状を、各超音波素子26n および各超音波素子28n がそれぞれ受信したエコー信号に基づいて算出する。図5は、第1アレイ26の各超音波素子26n 毎に示す送信信号Sn と上腕動脈18の管壁からのエコー信号En の波形を示すタイムチャートであり、その送信信号Sn とエコー信号En との時間差が第1アレイ26からの深さ寸法に対応している。このため、上記第1管壁位置演算手段48は、上記送信信号Sn とエコー信号En との間の時間差と、生体内の音速とに基づいて各超音波素子26n に対応する管壁の深さ位置をそれぞれ算出する。第2管壁位置演算手段50も同様にして、測定断面Bにおける各超音波素子28n に対応する管壁の深さ位置をそれぞれ逐次算出する。
次いで、測定断面形状算出手段52は、上記測定断面AおよびBにおいて、上記各超音波素子26n および28n に対応する管壁の深さ位置から、図6に示すようにX−Y座標内の点で特定し、さらに、それらの各点から曲線補完を用いてそれら各点を結ぶ閉曲線Kを算出して上腕動脈18の内腔の形状とし、その閉曲線Kの長径軸寸法2b、短軸径寸法2a、X軸に対する長径軸の傾斜角度θをそれぞれ逐次算出する。図6は測定断面Aにおける座標を示している。
中心軸算出手段54は、測定断面形状算出手段52により測定断面AおよびBにおいて求められた閉曲線Kから、その各測定断面AおよびBにおける閉曲線Kの中心点CA (X0A,Y0A,ZA )およびCB (X0B,Y0B,ZB )をそれぞれ算出する。そして、それら各測定断面AおよびBにおける閉曲線Kの中心点CA (X0A,Y0A,ZA )およびCB (X0B,Y0B,ZB )を結ぶ直線を、上腕動脈18の中心軸CLとして逐次決定する。図7は接触面22に平行な面(水平面)内の中心軸CLを示し、図8は接触面22に垂直な面(垂直面)内の中心軸CLを示している。たとえば、上記中心点CA は、閉曲線Kに近似する楕円の式(1) を利用した次式(1) から、最小自乗法を用いて未知数であるX0 およびY0 を算出する。
[(Xi −X0 )cos θ+( Yi −Y0 )sin θ]2/a2
+[ ( Xi −X0 )sin θ+( Yi −Y0 )cos θ]2/b2 =1 ・・・(1)
交差角算出手段56は、図7に示す接触面22に平行な面内において、測定断面A或いはBについて、中心軸算出手段54により算出された上腕動脈18の中心軸CLに基づいて中心点CA を通ってその中心軸CLと直交する上腕動脈18の直交断面A’或いは中心点CB を通ってその中心軸CLと直交する上腕動脈18の直交断面B’を決定し、その直交断面A’或いはB’と測定断面A或いはBとの交差角度β(度)を逐次算出する。同様に、図8に示す接触面22に垂直な面内において、測定断面A或いはBについて、中心軸算出手段54により算出された上腕動脈18の中心軸CLに基づいて中心点CA を通ってその中心軸CLと直交する上腕動脈18の直交断面A’或いは中心点CB を通ってその中心軸CLと直交する上腕動脈18の直交断面B’を決定し、その直交断面A’或いはB’と測定断面A或いはBとの交差角度γ(度)を逐次算出する。
補正手段58は、上記交差角度算出手段56により算出された交差角度βおよびγに基づいて、測定断面形状算出手段52により算出された形状を直交断面A’およびB’における形状となるように補正する。すなわち、上腕動脈18の直交断面A’およびB’における長軸径2b’および短軸径2a’となるように、補正式(2) 、(3) により測定断面形状算出手段52により算出された長軸径2bおよび短軸径2aをそれぞれ逐次補正する。
a’=a/cos β ・・・(2)
b’=b/cos γ ・・・(3)
但し、cos β=( X0 A−X0 B)/L
cos γ=( X0 A−X0 B)/L
直交断面積算出手段60は、直交断面A’およびB’における上腕動脈18の内腔面積SA ’およびSB ’を、上記補正手段58によって補正された長軸径2b’および短軸径2a’に基づいて逐次算出する。たとえば、直交断面A’およびB’における値に補正された長軸径2b’および短軸径2a’から特定される楕円の式(4) を用いて積分することにより、直交断面A’およびB’における面積SA ’およびSB ’が算出される。本実施例では、上記測定断面形状算出手段52、中心軸算出手段54、交差角算出手段56、補正手段58、および直交断面積算出手段60が、血管形状算出手段62に対応している。
(( Xi ’−X0 )/a’)2+(( Yi ’−Y0 )/b’)2=1 ・・・(4)
アクティブ素子選択手段64は、第3アレイ30の複数の超音波素子30n のうち上腕動脈18に最も近い素子或いは超音波放射方向線が上腕動脈18の中心軸CLと最も近いアクティブ素子を、中心軸算出手段54により求められた上腕動脈18の中心軸CLの位置に基づいて選択する。
血流速度算出手段66は、予め記憶された式(5) から、上記アクティブ素子から放射された超音波が血流によるドップラ効果によって位相変化或いは周波数変化させられたドップラ反射波の周波数fdに基づいて1拍毎に発生する最大瞬時血流速度Umax を算出する。式(5) において、fcは放射される超音波の波数、cは生体中の音速である。式(5) において、θ2 は図9に示す接触面22に垂直な面内のアクティブ素子からの超音波放射方向線USLと中心軸CLとの角度、θ1 は図10に示す接触面22に垂直な面内のアクティブ素子からの超音波放射方向線USLと中心軸CLとの角度、θ3 は図10に示す接触面22に垂直な面内の受波素子の受波方向線RSLと中心軸CLとの角度である。アクティブ素子が受波素子を兼ねる場合にはθ3 =0となる。それらの角度θ1 、θ2 、θ3 は、既知のアクティブ素子の幾何的位置と前記中心軸算出手段54により算出された中心軸CLとに基づいて予め算出されるとともに、上記式(5) は、それらの角度θ1 、θ2 、θ3 による補正が加味されているので、本実施例の血流速度算出手段66は、式(5) を用いることにより超音波放射方向線USLと中心軸CLとの間の相対角度θ1 、θ2 、θ3 に起因する最大瞬時血流速度Umax のずれを補正した最大瞬時血流速度Umax を一挙に算出しているので、超音波放射方向線USLと中心軸CLとの間の相対角度θ1 、θ2 、θ3 を算出する相対角度算出手段68、およびその相対角度θ1 、θ2 、θ3 に起因する最大瞬時血流速度Umax のずれを補正する血流補正手段70を兼ねているが、血流速度算出手段66は最大瞬時血流速度Umax を基本的に算出し、相対角度算出手段68は既知のアクティブ素子の幾何的位置と前記中心軸算出手段54により算出された中心軸CLとに基づいて超音波放射方向線USLと中心軸CLとの間の相対角度θ1 、θ2 、θ3 を算出し、血流補正手段70はその相対角度θ1 、θ2 、θ3 に起因する最大瞬時血流速度Umax のずれを補正するように役割を分担するように構成されてもよい。
fd=−(fc/c)(cos θ1cosθ2 +cos θ3cosθ2 ) ・・・(5)
血流量算出手段72は、前記直交断面積算出手段60により求められた直交断面B’における上腕動脈18のSB ’と、上記血流速度算出手段66において求められた補正後の最大瞬時血流速度Umax とに基づいて、直交断面B’における血流量QB (=SB ’×Umax /2)を逐次算出する。そして、出力手段74は、上記のようにして求められた、補正後の長軸径2b’および短軸径2a’、補正後の血流速度U、血流量QB をそれぞれ数字、グラフにて表示器38に画像表示或いは印字表示させるとともに、図示しない記憶装置に記憶させる。
図11および図12は、前記電子制御装置36の制御作動の要部を説明するフローチャートであって、図11は測定開始操作判定後に実行される測定制御ルーチンであり、図12はその測定制御ルーチン内の直交断面内の血管形状算出ルーチンを示している。
図11において、測定開始操作が判定されると、ステップ(以下、ステップを省略する)S1において、実際の測定実行時間Tが予め設定された最大測定時間Tmよりも小さいか否かが判断される。当初はこのS1の判断が肯定されるので、S2において走査素子番号nがクリアされた後、S3においてその走査素子番号nに「1」が加算される。次に、S4において、走査素子番号nがその最大番号nmax よりも低いか否かが判断される。この走査素子番号nは、第1アレイ或いは第2アレイの超音波素子数に対応した値に設定される。当初はこのS4の判断が肯定されるので、S5においてn番目の超音波素子26n および28n から超音波は発信され、S6においてn番目の超音波素子26n および28n によりエコー信号が受信される。次いで、前記第1血管壁位置演算手段48および第2血管壁位置演算手段50に対応するS7において、上腕動脈18のうちの上記n番目の超音波素子26n および28n に対応するに部位の血管壁の深さ位置が算出される。そして、上記S4の判断が否定されるまで上記S3以下が繰り返し実行される。
1番目の超音波素子261 および281 からnmax 番目の超音波素子26nmaxおよび28nmaxまでの送信受信が終了すると上記S4の判断が否定されるので、前記アクティブ素子選択手段64に対応するS8において、測定断Bにおける最大血管径位置或いは管壁の最大深さ位置が上記S7において算出された各超音波素子28n に対応するに部位の血管壁の深さ位置に基づいて決定されるとともに、第3アレイ30のうちの上記最大血管径位置或いは管壁の最大深さ位置に対応して位置するドップラ用超音波素子30n が血流検出に用いられるアクティブ素子として決定される。次いで、前記血流算出手段66に対応するS9において、よく知られた前記の式(5) から、上記アクティブ素子から放射された超音波が血流によるドップラ効果によって位相変化或いは周波数変化させられたドップラ反射波の周波数fdに基づいて最大瞬時血流速度Umax が算出される。
次のS10では、図12に示す直交断面の血管形状算出ルーチンが実行される。図12において、前記測定断面内形状算出手段52に対応するSA1では、測定断面AおよびBにおいて各超音波素子26n および28n に対応する管壁の深さ位置が図6に示すようにX−Y座標内の点で特定され、それらの各点から曲線補完を用いてそれら各点を結ぶ閉曲線Kが算出されて上腕動脈18の内腔の形状とされ、測定断面Aの場合の図6に示すように、その閉曲線Kの長径軸寸法2b、短軸径寸法2aが測定断面AおよびBにおいてそれぞれ算出される。次に、前記中心軸算出手段54に対応するSA2では、SA1により求められた各測定断面AおよびB内の閉曲線Kから、その閉曲線Kの中心点CA (X0A,Y0A,ZA )およびCB (X0B,Y0B,ZB )がそれぞれ算出され、それら各測定断面AおよびBにおける閉曲線Kの中心点CA (X0A,Y0A,ZA )およびCB (X0B,Y0B,ZB )を結ぶ直線が、上腕動脈18の中心軸CLとして決定される。続いて、前記交差角算出手段56に対応するSA3では、図7に示す接触面22に平行な面内において、測定断面Bについて、SA2により算出された上腕動脈18の中心軸CLに基づいて中心点CB を通ってその中心軸CLと直交する上腕動脈18の直交断面B’が決定され、その直交断面B’と測定断面Bとの交差角度β(度)が算出される。また、図8に示す接触面22に垂直な面内において、測定断面Bについて、中心軸算出手段54により算出された上腕動脈18の中心軸CLに基づいて中心点CB を通ってその中心軸CLと直交する上腕動脈18の直交断面B’が決定され、その直交断面B’と測定断面Bとの交差角度γ(度)が算出される。次に、前記補正手段58に対応するSA4では、上記SA3により算出された交差角度βおよびγに基づいて、SA1により算出された形状が直交断面B’における形状となるように補正する。すなわち、上腕動脈18の直交断面B’における長軸径2b’および短軸径2a’となるように、補正式(2) 、(3) により測定断面形状算出手段52により算出された長軸径2bおよび短軸径2aがそれぞれ補正される。そして、前記直交断面積算出手段60に対応するSA5では、直交断面B’における上腕動脈18の内腔面積SB ’が上記SA4によって補正された長軸径2b’および短軸径2a’に基づいて算出する。たとえば、直交断面B’における値に補正された長軸径2b’および短軸径2a’から特定される楕円の式(4) を用いて積分することにより、直交断面B’における面積SB ’が算出される。本実施例では、上記SA1乃至SA5が、血管形状算出手段62に対応している。
図11に戻って、前記相対角度算出手段68に対応するS11および前記血流速度補正手段70に対応するS12では、既知のアクティブ素子の幾何的位置と前記中心軸算出手段54により算出された中心軸CLとに基づいて超音波放射方向線USLと中心軸CLとの間の相対角度θ1 、θ2 、θ3 が算出されるとともに、その相対角度θ1 、θ2 、θ3 に起因する最大瞬時血流速度Umax のずれが補正される。本実施例では、S9において角度θ1 、θ2 、θ3 による補正が加味された式(5) が用いられることによって最大瞬時血流速度Umax が算出されているので、上記S11およびS12は、S9において実行されていることになる。
以上のようにして直交断面B’における上腕動脈18の内腔の直交断面積SB ’と最大瞬時血流速度Umax とが算出されると、前記血流量算出手段72に対応するS13において、上腕動脈18の血流量QB (=SB ’×Umax /2)が算出される。そして、前記出力手段74に対応するS14では、上記のようにして求められた、補正後の長軸径2b’および短軸径2a’、補正後の最大瞬時血流速度Umax 、血流量QB がそれぞれ数字、グラフにて表示器38に画像表示出力或いは印字表示出力させられるとともに、図示しない記憶装置に記憶させられる。このようにして、直交断面A’およびB’における上腕動脈18の内腔の長軸径2b’および短軸径2a’、その内腔の断面形状や、直交断面B’における直交断面積SB 、最大瞬時血流速度Umax 、および血流量QB が逐次求められるので、血流により血管壁に加えられるシェアストレスを計算できる。これは内皮機機能の評価に用いられ得る。
上述のように、本実施例によれば、血管形状算出手段62(S10)により、生体の一部14である上腕部の表皮20下に位置する上腕動脈18に交差してその表皮20上に配置された第1アレイ26および第2アレイ28から検出されるエコー信号に基づいて、その第1アレイ26および第2アレイ28の直下に位置する部分の血管壁のうち各第1超音波素子26n および第2超音波素子28n に対応する部位の位置がそれぞれ算出され、それら血管壁の第1超音波素子26n および第2超音波素子28n に対応する部位の位置に基づいて、上腕動脈18の直交断面における形状が算出されるので、生体皮膚20下の上腕動脈18の方向が上記第1アレイ26および第2アレイ28と直交していなくても、或いは皮膚20に平行でなくても、その上腕動脈18の断面形状が正確に得られる。
また、本実施例によれば、第1血管壁位置演算手段48(S7)および第2血管壁位置演算手段50(S7)によって、各第1超音波素子26n および各第2超音波素子28n から放射される放射信号とその各第1超音波素子26n および各第2超音波素子28n により検出される血管壁からのエコー(反射)信号との時間差に基づいてその血管壁までの距離がそれぞれ算出され、血管形状算出手段62(S10)により、その距離に基づいて前記第1アレイ26の測定断面A内および第2アレイ28の測定断面B内における血管壁の位置が決定されるので、第1アレイ26および第2アレイ28の測定断面A内および測定断面B内において各血管壁の位置が正確に得られる。
また、本実施例によれば、血管形状算出手段62(S10)において、第1アレイ26および第2アレイ28の測定断面内の血管壁の部位の位置に基づいて、その第1アレイ26の測定断面A内および前記第2アレイ28の測定断面B内における血管壁の断面の中心点CA およびCB と長軸径2bおよび短軸径2aとがそれぞれ算出され、その第1アレイ26の測定断面A内と前記第2アレイ28の測定断面B内とにおける血管壁の断面の中心点CA とCB とに基づいて上腕動脈18の中心軸CLが算出され、その上腕動脈18の中心軸CLに基づいてその血管の直交断面B’と測定断面Bとの交差角度βおよびγが算出され、その交差角度βおよびγに基づいて上記血管の直交断面B’内の長軸径2b’および短軸径2a’となるように、長軸径2bおよび短軸径2aが補正されるので、生体皮膚20下の上腕動脈18の方向が上記第1アレイ26および第2アレイ28と直交していなくても、或いは皮膚20に平行でなくても、その血管18の断面形状が正確に得られる。
また、本実施例によれば、血管形状算出手段62(S10)は、補正手段58(SA4)によって補正された長軸径2b’および短軸径2a’に基づいて血管18の直交断面B’における断面積SB ’を算出する直交断面積算出手段60(SA5)を含むものであるので、生体皮膚20下の血管18の方向が上記第1アレイ26および第2アレイ28と直交していなくても、或いは皮膚20に平行でなくても、その血管18の断面積SB ’が正確に得られる。
また、本実施例によれば、(a) 一方向に配列された複数の第1超音波素子26n を有する第1アレイ26と、その第1超音波素子26n の配列方向と平行な方向に配列された複数の第2超音波素子28n を有する第2アレイ28と、ドップラ用超音波素子30n を有する第3アレイ30とを備え、その第1アレイ26および第2アレイ28が生体の表皮20下に位置する上腕動脈18と交差し且つそのドップラ用超音波素子30n からの超音波を放射方向USLが上腕動脈18に対して鋭角を成すように装着される超音波プローブ12と、(b) ドップラ用超音波素子30n から上腕動脈18に向かって超音波を放射したときに得られる、その上腕動脈18内の血流速度に基づくドップラ効果により変化させられるドップラ反射波に基づいてその血流速度Uを算出する血流速度算出手段66(S9)と、(c) 第1アレイ26の各第1超音波素子26n から検出されるエコー(反射)信号に基づいて、上腕動脈18の第1アレイ26の直下に位置する部分の血管壁のうち各第1超音波素子26n に対応する部位の位置をそれぞれ算出する第1血管壁位置演算手段48(S7)と、(d) 第2アレイ28の各第2超音波素子28n から検出されるエコー信号に基づいて、上腕動脈18第2アレイ28の直下に位置する部分の血管壁のうち各第2超音波素子28n に対応する部位の位置をそれぞれ算出する第2血管壁位置演算手段50(S7)と、(e) 第1血管壁位置演算手段48により算出された各第1超音波素子26n に対応する血管壁部位の位置と第2血管壁位置演算手段50により算出された各第2超音波素子28n に対応する血管壁部位の位置とに基づいて、第1アレイ26の測定断面A内と第2アレイ28の測定断面B内とにおける血管壁の断面の中心点CA とCB とを算出し、第1アレイ26の測定断面Aと第2アレイ28の測定断面Bとにおける血管壁の断面の中心点CA とCB に基づいてその血管18の中心軸CLを算出する中心軸算出手段54(SA2)と、(f) その中心軸算出手段54により算出された血管18の中心軸CLとドップラ用超音波素子30n からその血管18に向かって超音波を放射する方向USLとの実際の相対角度θ1 、θ2 、θ3 を算出する相対角度算出手段68(S11)と、(g) その相対角度算出手段68により算出された実際の相対角度θ1 、θ2 、θ3 に基づいて血流速度算出手段66により算出された血流速度を補正する血流速度補正手段70(S12)とが、設けられているので、生体皮膚20下の血管18の方向が上記第1アレイ26および第2アレイ28と直交していなくても、或いは皮膚20に平行でなくても、その血管18内の最大瞬時血流速度Umax が正確に得られる。
また、本実施例によれば、(a) 上記の血流速度測定のための構成と、(b) 第1血管壁位置演算手段48により算出された各第1超音波素子26n に対応する血管壁部位の位置と第2血管壁位置演算手段50により算出された各第2超音波素子28n に対応する血管壁部位の位置とに基づいて、第1アレイ26の測定断面A内と第2アレイ28の測定断面B内とにおける血管壁の断面の中心点CA およびCB と長軸径2bおよび短軸径2とをそれぞれ算出する測定断面内形状算出手段52(SA1)と、(c) 中心軸算出手段54により算出された血管18の中心軸CLに基づいて、その血管18の直交断面B’と測定断面Bとの交差角度β、γを算出する交差角度算出手段56(SA3)と、(d) その交差角度算出手段56により算出された交差角度β、γに基づいて、血管18の直交断面B’内の長軸径2/b’および短軸径2/a’となるように、測定断面形状算出手段52により算出された長軸径2/bおよび短軸径2/aを補正する補正手段58(SA4)と、(e) その補正手段58によって補正された長軸径2/bおよび短軸径2/aに基づいて血管18の直交断面B’における断面積S’を算出する直交断面積算出手段60(SA5)と、(f) その直交断面積算出手段60により算出された血管18の直交断面B’における断面積S’と血流速度補正手段70により補正された最大瞬時血流速度Umax とに基づいて、血管18内の血流量QB を算出する血流量算出手段72(S13)とが、設けられているので、生体皮膚20下の血管18の方向が上記第1アレイ26および第2アレイ28と直交していなくても、或いは皮膚20に平行でなくても、その血管18内の血流量QB が正確に得られる。
以上、本発明の実施例を図面に基づいて詳細に説明したが、本発明はその他の態様においても適用され得る。
たとえば、前述の実施例では、血流量算出手段72により直交断面B’における上腕動脈18の血流量QB (=SB ’×Umax /2)が算出されていたが、直交断面A’における上腕動脈18の血流量QA (=SA ’×Umax /2)が、直交断面積算出手段60により求められた直交断面A’における上腕動脈18の直交断面積SA ’に基づいて算出されてもよい。また、直交断面A’と直交断面B’との平均値に対する血流量が求められてもよい。
また、直交断面A’およびB’における上腕動脈18の内腔の長軸径2b’および短軸径2a’、その内腔の断面形状や、直交断面B’における直交断面積SB 等の変化の時間差に基づいて脈波伝播速度が算出されてもよい。
また、前述の実施例において、血流量算出手段72により直交断面B’における上腕動脈18の平均の血流量QB (=SB ’×Umax /2)が算出されていたが、血流速度算出手段66によって算出される血流速度Uが平均速度である場合には、上記血流量QB は式(QB =SB ’×U)に基づいて算出される。また、血流量算出手段72により、最大血流量QBmax(=SB ’×Umax )が求められてもよい。
また、前述の実施例において、第3アレイ30は複数のドップラ用超音波素子30n から構成されていたが、単一のドップラ用超音波素子から構成されていたもよい。
また、前述の実施例では、たとえば第1アレイ26では、各超音波素子26n が、図5に示すように、個々に超音波の発信および受信を行っていたが、位相が異なる駆動信号を用いて数個の超音波素子から超音波を送信して細い超音波ビームを送信し、受信する場合も数個の超音波素子を用いて受信する所謂ビームフォーミング技術を採用することもできる。これによれば、所定の距離で収束するように超音波ビームを形成することができるので、検出精度が高められる。
なお、上述したのは、あくまでも一実施形態であり、本発明は当業者の知識に基づいて種々の変更、改良を加えた態様で実施することができる。
本発明の一実施例の超音波装置に備えられた超音波プローブが上腕に装着された状態を示す斜視図である。 図1の実施例の超音波装置の超音波プローブに設けられた第1アレイの測定断面および第2アレイの測定断面と上腕動脈との関係を示す斜視図である。 図1の実施例の超音波装置の要部構成を概略説明する図である。 図3の電子制御装置の制御機能の要部を説明する機能ブロック線図である。 図1の実施例において、第1アレイの測定断面における各超音波素子の送信波および受信波を説明する図である。 図5の第1血管壁位置演算手段により図5の送信波と受信波との時間差に基づいて算出された第1アレイの測定断面における各測定点の座標と、それらを接続した閉曲線である上腕動脈の内腔の形状を説明する図である。 図5の中心線算出手段により求められた中心線CLに直交する直交断面B’と第2アレイの測定断面Bとの、接触面に平行な平面内の交差角度βを示す図である。 図5の中心線算出手段により求められた中心線CLに直交する直交断面B’と第2アレイの測定断面Bとの、接触面に垂直な面内の交差角度γを示す図である。 図5の中心線算出手段により求められた中心線CLとドップラ用超音波素子の放射方向線USLとの、接触面に垂直な面内の相対角度θ2 を示す図である。 図5の中心線算出手段により求められた中心線CLとドップラ用超音波素子の放射方向線USLとの、接触面に平行な平面内の相対角度θ1 と、中心線CLとドップラ用超音波素子の受波方向線RSLとの、接触面に平行な平面内の相対角度θ3 とを示す図である。 図3の電子制御装置の制御作動の要部を説明するフローチャートである。 図11のS10の直交断面内における血管形状を算出するためのルーチンを示す図である。
符号の説明
10:超音波測定装置(血管形状測定装置、血流速度測定装置、血流量測定装置)
12:超音波プローブ
14:生体の一部(上腕部)
18:上腕動脈(血管)
26:第1アレイ
26n :第1超音波素子
28:第1アレイ
28n :第2超音波素子
48:第1血管壁位置演算手段
50:第2血管壁位置演算手段
52:血管断面内形状算出手段
54:中心軸算出手段
56:交差角算出手段
58:形状補正手段
60:直交断面積算出手段
62:血管形状算出手段
66:血流速度算出手段
68:相対角度算出手段
70:血流速度補正手段
72:血流量算出手段
























Claims (3)

  1. 一方向に配列された複数の第1超音波素子を有する第1アレイと、該第1超音波素子の配列方向と平行な方向に配列された複数の第2超音波素子を有する第2アレイとを備え、生体の表皮下に位置する血管に交差して該表皮上に配置された該第1アレイおよび第2アレイから検出されるエコー信号に基づいて該血管の形状を測定する血管形状測定装置であって、
    前記第1アレイの各第1超音波素子から放射される放射信号と該各第1超音波素子により検出される前記血管壁からの反射信号との時間差に基づいて該血管壁までの距離をそれぞれ算出し、該距離に基づいて前記第1アレイの測定断面内における血管壁の位置を決定する第1血管壁位置演算手段と、
    前記第2アレイの各第2超音波素子から放射される放射信号と各第2超音波素子により検出される前記血管壁からの反射信号との時間差に基づいて該血管壁までの距離をそれぞれ算出し、該距離に基づいて前記第2アレイの測定断面内における血管壁の位置を決定する第2血管壁位置演算手段と、
    前記第1血管壁位置演算手段により決定された各第1超音波素子に対応する血管壁部位の位置と前記第2血管壁位置演算手段により決定された各第2超音波素子に対応する血管壁部位の位置とに基づいて、前記血管の直交断面における形状を算出する血管形状算出手段と
    を、含み、
    該血管形状算出手段は、
    前記第1血管壁位置演算手段により決定された各第1超音波素子に対応する血管壁部位の位置と前記第2血管壁位置演算手段により決定された各第2超音波素子に対応する血管壁部位の位置とに基づいて、前記第1アレイの測定断面内と前記第2アレイの測定断面内とにおける血管壁の断面の中心点と長軸径および/または短軸径とをそれぞれ算出する測定断面内形状算出手段と、
    該測定断面内形状算出手段により算出された第1アレイの測定断面と前記第2アレイの測定断面とにおける血管壁の断面の中心点に基づいて該血管の中心軸を算出する中心軸算出手段と、
    該中心軸算出手段により算出された血管の中心軸に基づいて、該血管の直交断面と前記測定断面との交差角度を算出する交差角度算出手段と、
    該交差角度算出手段により算出された交差角度に基づいて、前記血管の直交断面の長軸径および/または短軸径となるように、前記形状算出手段により算出された長軸径および/または短軸径を補正する補正手段と
    を、含むことを特徴とする血管形状測定装置。
  2. 前記血管形状算出手段は、前記補正手段によって補正された長軸径および短軸径に基づいて前記血管の直交断面における断面積を算出する直交断面積算出手段を、さらに含むものである請求項の血管形状測定装置。
  3. 一方向に配列された複数の第1超音波素子を有する第1アレイと、該第1超音波素子の配列方向と平行な方向に配列された複数の第2超音波素子を有する第2アレイと、ドップラ用超音波素子とを備え、該第1アレイおよび第2アレイが前記生体の表皮下に位置する血管と交差し且つ該ドップラ用超音波素子からの超音波を放射方向が該血管に対して鋭角を成すように装着される超音波プローブと、
    前記ドップラ用超音波素子から前記血管に向かって超音波を放射したときに得られる、該血管内の血流速度に基づくドップラ効果により変化させられるドップラ反射波に基づいて該血流速度を算出する血流速度算出手段と、
    前記第1アレイの各第1超音波素子から検出される反射信号に基づいて、前記血管の該第1アレイの直下に位置する部分の血管壁のうち該各第1超音波素子に対応する部位の位置をそれぞれ算出する第1血管壁位置演算手段と、
    前記第2アレイの各第2超音波素子から検出されるエコー信号に基づいて、前記血管の該第2アレイの直下に位置する部分の血管壁のうち該各第2超音波素子に対応する部位の位置をそれぞれ算出する第2血管壁位置演算手段と、
    前記第1血管壁位置演算手段により算出された各第1超音波素子に対応する血管壁部位の位置と前記第2血管壁位置演算手段により算出された各第2超音波素子に対応する血管壁部位の位置とに基づいて、前記第1アレイの測定断面内と前記第2アレイの測定断面内とにおける血管壁の断面の中心点を算出し、該第1アレイの測定断面と前記第2アレイの測定断面とにおける血管壁の断面の中心点に基づいて該血管の中心軸を算出する中心軸算出手段と、
    該中心軸算出手段により算出された前記血管の中心軸と前記ドップラ用超音波素子から該血管に向かって超音波を放射する方向との実際の相対角度を算出する相対角度算出手段と、
    該相対角度算出手段により算出された実際の相対角度に基づいて前記血流速度算出手段により算出された血流速度を補正する血流速度補正手段と
    を、含む血流速度測定装置を備えた血流量測定装置であって、
    前記第1血管壁位置演算手段により算出された各第1超音波素子に対応する血管壁部位の位置と前記第2血管壁位置演算手段により算出された各第2超音波素子に対応する血管壁部位の位置とに基づいて、前記第1アレイの測定断面内と前記第2アレイの測定断面内とにおける血管壁の断面の中心点と長軸径および/または短軸径とをそれぞれ算出する測定断面内形状算出手段と、
    前記中心軸算出手段により算出された血管の中心軸に基づいて、該血管の直交断面と前記測定断面との交差角度を算出する交差角度算出手段と、
    該交差角度算出手段により算出された交差角度に基づいて、前記血管の直交断面の長軸径および/または短軸径となるように、前記形状算出手段により算出された長軸径および/または短軸径を補正する補正手段と、
    該補正手段によって補正された長軸径および短軸径に基づいて前記血管の直交断面における断面積を算出する直交断面積算出手段と、
    該直交断面積算出手段により算出された前記血管の直交断面における断面積と前記血流速度補正手段により補正された血流速度とに基づいて、前記血管内の血流量を算出する血流量算出手段と
    を、含むことを特徴とする血流量測定装置。
JP2004219866A 2004-07-28 2004-07-28 血管形状測定装置および血流量測定装置 Expired - Fee Related JP4441664B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2004219866A JP4441664B2 (ja) 2004-07-28 2004-07-28 血管形状測定装置および血流量測定装置
PCT/JP2005/013813 WO2006011544A1 (ja) 2004-07-28 2005-07-28 血管形状測定装置、血流速度測定装置、および血流量測定装置
US10/565,675 US7429244B2 (en) 2004-07-28 2005-07-28 Blood-vessel-shape measuring apparatus, blood-flow-velocity measuring apparatus, and blood-flow-amount measuring apparatus
EP05767223A EP1772102A4 (en) 2004-07-28 2005-07-28 INSTRUMENT FOR MEASURING THE FORM OF A BLOOD VESSEL, INSTRUMENT FOR MEASURING THE SPEED OF BLOOD FLOW AND APPARATUS FOR MEASURING BLOOD FLOW

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004219866A JP4441664B2 (ja) 2004-07-28 2004-07-28 血管形状測定装置および血流量測定装置

Publications (2)

Publication Number Publication Date
JP2006034667A JP2006034667A (ja) 2006-02-09
JP4441664B2 true JP4441664B2 (ja) 2010-03-31

Family

ID=35786292

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004219866A Expired - Fee Related JP4441664B2 (ja) 2004-07-28 2004-07-28 血管形状測定装置および血流量測定装置

Country Status (4)

Country Link
US (1) US7429244B2 (ja)
EP (1) EP1772102A4 (ja)
JP (1) JP4441664B2 (ja)
WO (1) WO2006011544A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012085789A (ja) * 2010-10-19 2012-05-10 Seiko Epson Corp 血管径測定装置
JP2015107174A (ja) * 2013-12-03 2015-06-11 株式会社ユネクス 生体血管状態測定装置

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102004043676B4 (de) * 2004-09-09 2014-01-09 Siemens Aktiengesellschaft Verfahren zur Visualisierung von Plaqueablagerungen aus 3D-Bilddatensätzen von Gefäßstrukturen
US7940974B2 (en) * 2006-11-21 2011-05-10 General Electric Company Method and system for adjusting 3D CT vessel segmentation
US8160395B2 (en) 2006-11-22 2012-04-17 General Electric Company Method and apparatus for synchronizing corresponding landmarks among a plurality of images
US8126238B2 (en) * 2006-11-22 2012-02-28 General Electric Company Method and system for automatically identifying and displaying vessel plaque views
US8244015B2 (en) * 2006-11-22 2012-08-14 General Electric Company Methods and apparatus for detecting aneurysm in vasculatures
US7957574B2 (en) * 2006-11-22 2011-06-07 General Electric Company Methods and apparatus for generating a risk metric for soft plaque in vessels
FR2908976B1 (fr) * 2006-11-24 2009-02-20 Gen Electric Procede de mesures de dimensions d'un vaisseau.
JP5499938B2 (ja) 2010-06-25 2014-05-21 セイコーエプソン株式会社 超音波センサー、測定装置、プローブ、および測定システム
JP5499939B2 (ja) * 2010-06-25 2014-05-21 セイコーエプソン株式会社 測定装置、生体検査装置、流速測定方法、および圧力測定方法
JP2015166024A (ja) * 2015-07-02 2015-09-24 セイコーエプソン株式会社 超音波センサーおよび超音波測定装置
US11918413B2 (en) 2015-11-03 2024-03-05 Koninklijke Philips N.V. Ultrasound probe, system and method for measuring arterial parameters using non-imaging ultrasound
US20180092622A1 (en) * 2016-09-30 2018-04-05 Robert Bosch Gmbh Phased Array for Detecting Artery Location to Measure Blood Velocity
US20180092621A1 (en) * 2016-09-30 2018-04-05 Robert Bosch Gmbh Single piezoelectric transmitter and receiver to detect blood velocities
WO2020175515A1 (ja) * 2019-02-27 2020-09-03 株式会社村田製作所 血管位置検出装置、および、血流量計測装置
JP2023028891A (ja) * 2021-08-20 2023-03-03 日本光電工業株式会社 生体情報処理方法、生体情報処理装置及びプログラム
CN115568876B (zh) * 2022-11-24 2023-05-09 苏州圣泽医疗科技有限公司 血流速度测量值的修正方法及多普勒血流检测装置

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1991011146A1 (en) * 1990-01-25 1991-08-08 Commonwealth Scientific And Industrial Research Organisation Ultrasonic method and apparatus for determination of vessel location and size
JPH0556971A (ja) 1991-08-30 1993-03-09 Toshiba Corp 超音波診断装置
US5339816A (en) * 1991-10-23 1994-08-23 Aloka Co., Ltd. Ultrasonic doppler blood flow monitoring system
JP3381968B2 (ja) 1993-07-09 2003-03-04 株式会社東芝 超音波診断装置
US6261233B1 (en) * 1996-01-05 2001-07-17 Sunlight Medical Ltd. Method and device for a blood velocity determination
FR2759892A1 (fr) 1996-12-31 1998-08-28 Philips Electronics Nv Systeme d'echographie ultrasonore pour l'examen des arteres
JP3857788B2 (ja) 1997-09-01 2006-12-13 テルモ株式会社 循環器情報計測システム
AUPP227898A0 (en) 1998-03-11 1998-04-09 Commonwealth Scientific And Industrial Research Organisation Improvements in ultrasound techniques
JP3378941B2 (ja) 2000-06-29 2003-02-17 徳島県 位置自動追尾型超音波血流計測装置
US6780155B2 (en) * 2001-12-18 2004-08-24 Koninklijke Philips Electronics Method and system for ultrasound blood flow imaging and volume flow calculations
JP2004129797A (ja) * 2002-10-09 2004-04-30 Matsushita Electric Ind Co Ltd 超音波診断装置

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012085789A (ja) * 2010-10-19 2012-05-10 Seiko Epson Corp 血管径測定装置
US9642594B2 (en) 2010-10-19 2017-05-09 Seiko Epson Corporation Blood vessel diameter measurement device
JP2015107174A (ja) * 2013-12-03 2015-06-11 株式会社ユネクス 生体血管状態測定装置

Also Published As

Publication number Publication date
US7429244B2 (en) 2008-09-30
WO2006011544A1 (ja) 2006-02-02
EP1772102A1 (en) 2007-04-11
EP1772102A4 (en) 2010-11-10
US20060241427A1 (en) 2006-10-26
JP2006034667A (ja) 2006-02-09

Similar Documents

Publication Publication Date Title
JP4441664B2 (ja) 血管形状測定装置および血流量測定装置
US6673020B2 (en) Ultrasonic diagnostic apparatus
JP5499939B2 (ja) 測定装置、生体検査装置、流速測定方法、および圧力測定方法
EP1273267B1 (en) Ultrasonic diagnostic apparatus
WO2019205167A1 (zh) 一种超声瞬时弹性测量设备及方法
EP1079242B1 (en) Ultrasound diagnostic apparatus
JP7456151B2 (ja) 超音波診断装置、超音波診断装置の制御方法、及び、超音波診断装置の制御プログラム
EP1934629B1 (en) A method for generating hardness and/or strain information of a tissue
JP2001218768A (ja) 超音波診断装置
JP6861624B2 (ja) 超音波送受信装置および超音波送受信方法
JP5199157B2 (ja) 超音波診断装置、圧較差測定方法及び血管弾性測定方法
JP4768100B2 (ja) 超音波診断装置
JP2018102589A (ja) 脈波伝播速度測定装置、血圧測定装置、および脈波伝播速度測定方法
JP2005074146A (ja) 超音波測定方法、超音波発生機構
CN109561880B (zh) 用于确定心输出量的系统
JP4620423B2 (ja) 平滑筋弛緩状態評価装置
JP2012192133A (ja) 超音波診断装置および超音波画像生成方法
JP2002224063A (ja) 脈波伝播速度計測装置及び超音波診断装置
JP2970884B2 (ja) 血管弾性率測定用プローブ装置
JP4230904B2 (ja) 超音波診断装置
JP7224016B2 (ja) 血管径センサおよび血管径測定装置
JP2021003214A (ja) 超音波診断装置、超音波診断装置の制御方法、及び、超音波診断装置の制御プログラム
US20150374330A1 (en) Blood vessel search device, ultrasonic measurement apparatus, and blood vessel search method
WO2014129203A1 (ja) 超音波診断装置および超音波診断装置の制御方法
JP2017192479A (ja) 画像生成装置及び画像生成方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20061031

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090811

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20091009

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20091009

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20091009

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

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20091215

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20091215

R150 Certificate of patent or registration of utility model

Ref document number: 4441664

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130122

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20160122

Year of fee payment: 6

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees