JPWO2018159708A1 - 血流解析装置および方法並びにプログラム - Google Patents

血流解析装置および方法並びにプログラム Download PDF

Info

Publication number
JPWO2018159708A1
JPWO2018159708A1 JP2019503082A JP2019503082A JPWO2018159708A1 JP WO2018159708 A1 JPWO2018159708 A1 JP WO2018159708A1 JP 2019503082 A JP2019503082 A JP 2019503082A JP 2019503082 A JP2019503082 A JP 2019503082A JP WO2018159708 A1 JPWO2018159708 A1 JP WO2018159708A1
Authority
JP
Japan
Prior art keywords
blood vessel
blood flow
blood
backflow
flow velocity
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2019503082A
Other languages
English (en)
Other versions
JP6690053B2 (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 JPWO2018159708A1 publication Critical patent/JPWO2018159708A1/ja
Application granted granted Critical
Publication of JP6690053B2 publication Critical patent/JP6690053B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, 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/026—Measuring blood flow
    • A61B5/0263—Measuring blood flow using NMR
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computerised tomographs
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50—Clinical applications
    • A61B6/504—Clinical applications involving diagnosis of blood vessels, e.g. by angiography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50—Clinical applications
    • A61B6/507—Clinical applications involving determination of haemodynamic parameters, e.g. perfusion CT
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06—Measuring blood flow
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0891—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461—Displaying means of special interest
    • A61B8/466—Displaying means of special interest adapted to display 3D data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/483—Diagnostic techniques involving the acquisition of a 3D volume of data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/488—Diagnostic techniques involving Doppler signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223—Devices 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
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, 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/026—Measuring blood flow
    • A61B5/0285—Measuring or recording phase velocity of blood waves
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays
    • A61B5/7425—Displaying combinations of multiple images regardless of image source, e.g. displaying a reference anatomical image with a live image

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Hematology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Physiology (AREA)
  • Optics & Photonics (AREA)
  • Cardiology (AREA)
  • Vascular Medicine (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Pulmonology (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

想定外の箇所に逆流があった場合でも、その逆流箇所をとらえることができ、逆流箇所の見落としリスクを減らすことができる血流解析装置および方法並びにプログラムを提供する。血管を含む被写体を撮影した血管画像から血管領域を抽出する血管領域抽出部（１１）と、血管領域内における血流速度および血流方向を表す流速ベクトルを取得する流速ベクトル取得部（１２）と、血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する逆流特定部（１３）とを備える。

Description

本発明は、医用画像から抽出された血管領域を用いて、その血管領域内を流れる血流の解析を行う血流解析装置および方法並びにプログラムに関するものである。
近年、心臓などを撮影した医用画像を用いて、心臓などの血管内の血流を解析することが行われている。
このような医用画像を用いた血流解析方法としては、たとえば３次元シネ位相コントラスト磁気共鳴法によって撮影されたＭＲＩ（Magnetic Resonance Imaging）画像を用いて、ボクセル毎、ピクセル毎または領域毎に流速ベクトルを算出し、これを表示する方法が提案されている。
また、特開２０１４−１８８３２３号公報においては、医用画像から血管モデルを生成し、その血管モデルに基づいて、ＣＦＤ（Computational Fluid Dynamics）を行って流速ベクトルを求めることが提案されている。
ここで、血流解析においては、血流が逆流している箇所を特定することが、診断上重要である。上述したように流速ベクトルを表示した場合、ユーザが、その向きと大きさを観察することによって逆流箇所を特定することができる。
また、医用画像から抽出された血管領域に関心領域を設定し、その関心領域に流れる血流を表す流線を表示させて観察することにより、逆流を確認することもできる。
また、超音波診断の分野でも逆流を確認することが行われている。たとえば特開２００９−１５３９１９号公報および特開平８−１５０１４２号公報においては、超音波診断装置を用いてカラードップラ計測を行って、逆流を可視化することが提案されている。
しかしながら、上述した３次元シネ位相コントラスト磁気共鳴法およびＣＦＤを用いて逆流箇所を確認する場合、ユーザが関心のある箇所付近の流速ベクトルの向きを確認したり、関心領域を設定した際の流線の形成具合を確認したりすることによって、逆流箇所を確認することができるが、この場合、ユーザが目的を持って観察しようとした場所以外で逆流があった場合、逆流を見落とす可能性がある。
また、超音波診断によって逆流を確認する場合、ユーザが関心のある領域にプローブを当ててカラードップラ計測が行われるが、超音波画像は２次元的であり、２次元画像上での流れの正負を表示させることはできるが、血管形状に沿った流れを抽出することは困難であり、逆流として画像化することは難しい。
本発明は、上記事情に鑑み、想定外の箇所に逆流があった場合でも、その逆流箇所をとらえることができ、逆流箇所の見落としリスクを減らすことができる血流解析装置および方法並びにプログラムを提供することを目的とするものである。
本発明の血流解析装置は、血管を含む被写体を撮影した血管画像から血管領域を抽出する血管領域抽出部と、血管領域内における血流速度および血流方向を表す流速ベクトルを取得する流速ベクトル取得部と、血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する逆流特定部とを備える。
また、上記本発明の血流解析装置において、逆流特定部は、積分値が単調増加または単調減少でない場合に、逆流が有ると特定することができる。
また、上記本発明の血流解析装置においては、逆流の箇所を血管領域の座標軸上にマッピングしたマップ画像を表示部に表示させる表示制御部を備えることができる。
また、上記本発明の血流解析装置において、表示制御部は、流速ベクトルを予め設定された期間だけ積分した積分値を用いてマップ画像を生成することができる。
また、上記本発明の血流解析装置において、上記予め設定された期間は、１心拍の期間とすることができる。
また、上記本発明の血流解析装置において、表示制御部は、流速ベクトルを時間軸方向に積分する毎にマップ画像を時系列に生成して動画として表示部に表示させることができる。
また、上記本発明の血流解析装置において、表示制御部は、マップ画像を２次元または３次元の血管画像上に重畳表示させることができる。
また、上記本発明の血流解析装置において、表示制御部は、積分値に基づいて、逆流の量の情報を表示部にさせることができる。
また、上記本発明の血流解析装置において、表示制御部は、マップ画像以外の血流を示す指標を表示部に表示させることができる。
また、上記本発明の血流解析装置においては、血流を示す指標として、流速ベクトル、流線または流跡線を表示させることができる。
また、上記本発明の血流解析装置において、逆流特定部は、血管領域のボクセル毎または複数のボクセル毎に積分値を算出することができる。
また、上記本発明の血流解析装置において、逆流特定部は、血管領域の中心軸方向に直交する断面毎に積分値を算出することができる。
また、上記本発明の血流解析装置において、流速ベクトル取得部は、時系列に撮影された３次元の血管画像または任意の時点において撮影された３次元の血管画像に基づいて、流速ベクトルを取得することができる。
また、上記本発明の血流解析装置において、流速ベクトル取得部は、血管領域の数値流体力学(Computational Fluid Dynamics)モデル、３次元シネ位相コントラスト磁気共鳴法または３次元の超音波画像を用いて流速ベクトルを取得することができる。
本発明の血流解析方法は、血管を含む被写体を撮影した血管画像から血管領域を抽出し、血管領域内における血流速度および血流方向を表す流速ベクトルを取得し、血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する。
本発明の血流解析プログラムは、コンピュータを、血管を含む被写体を撮影した血管画像から血管領域を抽出する血管領域抽出部と、血管領域内における血流速度および血流方向を表す流速ベクトルを取得する血流ベクトル取得部と、血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する逆流特定部として機能させる。
本発明の他の血流解析装置は、コンピュータに実行させるための命令を記憶するメモリと、記憶された命令を実行するよう構成されたプロセッサを備え、プロセッサが、血管を含む被写体を撮影した血管画像から血管領域を抽出する処理と、血管領域内における血流速度および血流方向を表す流速ベクトルを取得する処理と、血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する処理とを実行するように構成されている。
本発明の血流解析装置および方法並びにプログラムによれば、血管を含む被写体を撮影した血管画像から血管領域を抽出し、血管領域内における血流速度および血流方向を表す流速ベクトルを取得する。そして、血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する。このように流速ベクトルの中心軸方向成分の積分値に基づいて逆流の有無を特定することによって、想定外の箇所に逆流があった場合でも、その逆流箇所をとらえることができ、逆流箇所の見落としリスクを減らすことができる。
本発明の血流解析装置の一実施形態を用いた血流解析システムの概略構成を示すブロック図 血管領域内における流速ベクトルを模式的に示す図 逆流箇所の特定方法を説明するための図 積分値が単調増加または単調減少しない例を示す図 断面毎に積分値を算出した場合の逆流箇所の表示例を示す図 本発明の血流解析装置の一実施形態を用いた血流解析システムの作用を説明するためのフローチャート 複数のボクセル毎に積分値を算出した場合の逆流箇所の表示例を示す図
以下、図面を参照して本発明の血流解析装置の一実施形態を用いた血流解析システムについて説明する。図１は、本実施形態の血流解析システムの概略構成を示すブロック図である。
本実施形態の血流解析システムは、具体的には、図１に示すように、血流解析装置１と、医用画像保管サーバ２と、表示部３と、入力部４とを備えている。
血流解析装置１は、コンピュータに本実施形態の血流解析プログラムをインストールしたものである。
血流解析装置１は、中央処理装置（ＣＰＵ（central processing unit））、半導体メモリ、およびハードディスクやＳＳＤ(Solid State Drive)等のストレージデバイスを備えている。ストレージデバイスには、本実施形態の血流解析プログラムがインストールされており、この血流解析プログラムが中央処理装置によって実行されることによって、図１に示す画像取得部１０、血管領域抽出部１１、流速ベクトル取得部１２、逆流特定部１３および表示制御部１４が動作する。
血流解析プログラムは、ＤＶＤ（Digital Versatile Disc）およびＣＤ−ＲＯＭ（Compact Disc Read Only Memory）などの記録媒体に記録されて配布され、その記録媒体からコンピュータにインストールされる。または、血流解析プログラムは、ネットワークに接続されたサーバコンピュータの記憶装置もしくはネットワークストレージに対して、外部からアクセス可能な状態で記憶される。そして、外部からの要求に応じてコンピュータにダウンロードされ、インストールされる。
画像取得部１０は、予め撮影された患者の医用画像を取得するものである。３次元医用画像６は、本実施形態においては、たとえばＣＴ（Computed Tomography）装置、ＭＲＩ
（Magnetic Resonance Imaging）装置および超音波撮影装置などによって撮影された３次元医用画像６を取得する。本実施形態では、患者の心臓（本発明の被写体に相当する）の３次元医用画像６（本発明の血管画像に相当する）を取得する場合について説明するが、これに限らず、肺、肝臓および頭部などその他の臓器でもよい。
３次元医用画像６は、医用画像保管サーバ２に患者の識別情報とともに予め保管されており、画像取得部１０は、入力部４などを用いてユーザによって入力された患者の識別情報に基づいて、その識別情報を有する３次元医用画像６を医用画像保管サーバ２から読み出して一時記憶するものである。
血管領域抽出部１１は、３次元医用画像６から血管領域を抽出するものである。具体的には、本実施形態の血管領域抽出部１１は、心臓の３次元医用画像６から冠動脈の領域を血管領域として抽出する。血管領域抽出部１１は、たとえば心臓の３次元医用画像６に対して多重解像度変換を行い、各解像度の画像に対してヘッセ行列の固有値解析を行い、各解像度の画像における解析結果を統合することによって、心臓領域中の様々なサイズの線構造（血管）の集合体として、冠動脈の領域を抽出する（たとえばY Sato, et al.、「Three-dimensional multi-scale line filter for segmentation and visualization of curvilinear structures in medical images.」、Medical Image Analysis、1998年6月、Vol.2、No.2、p.p.143-168など参照）。また、さらに最小全域木アルゴリズム等を用いて、抽出された各線構造の中心点を連結することにより、冠動脈を表す木構造のデータを生成し、抽出された冠動脈の中心点を結ぶ芯線上の各点（木構造データの各ノード）において、芯線に直交する断面を求め、各断面において、グラフカット法などの公知のセグメンテーション手法を用いて冠動脈の輪郭を認識し、その輪郭を表す情報を木構造データの各ノードに関連づけることによって、冠動脈の領域を抽出するようにしてもよい。
なお、冠動脈の領域の抽出方法としては上記の方法に限らず、領域拡張法などのその他の公知な手法を用いるようにしてもよい。
流速ベクトル取得部１２は、血管領域内における血流速度および血流方向を表す流速ベクトルを取得するものである。流速ベクトルの取得方法としては、種々の方法を採用することができる。具体的には、たとえば血管領域抽出部１１によって抽出された血管領域を用いてＣＦＤ（Computational Fluid Dynamics）による血流解析を行うことによって流速ベクトルを取得することができる。なお、この際、任意の時点において撮影された３次元医用画像６から抽出された血管領域を用いて血流解析を行うようにしてもよいし、時系列に撮影された３次元医用画像６のそれぞれから抽出された血管領域を用いて血流解析を行うようにしてもよい。
また、画像取得部１０において３次元シネ位相コントラスト磁気共鳴法（3D cine PC MRI）によって撮影されたＭＲＩ画像を取得し、そのＭＲＩ画像に基づいて取得された血管領域内の速度情報を用いて流速ベクトルを取得するようにしてもよい。
また、画像取得部１０においてドップラー計測によって時系列に撮影された３次元の超音波画像を取得し、その超音波画像に基づいて取得された血管領域内の速度情報を用いて流速ベクトルを取得するようにしてもよい。図２は、血管領域Ｂ内における流速ベクトルＶを模式的に示した図である。
逆流特定部１３は、流速ベクトル取得部１２によって取得された血管領域内の流速ベクトルに基づいて、血管領域内の血流の逆流の有無を特定するものである。
逆流特定部１３は、具体的には、図２に示すように、血管領域Ｂの延伸方向に延びる中心軸Ｃを設定し、流速ベクトルＶの中心軸Ｃの方向成分を時間軸方向に積分した積分値を算出し、その積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する。
流速ベクトルＶの中心軸Ｃの方向成分については、たとえば図２および図３に示すように、中心軸Ｃ上の各点について、中心軸Ｃに直交する断面Ｓを設定し、その断面Ｓ内の流速ベクトルＶの中心軸方向成分Ｖｃを算出する。そして、中心軸方向成分Ｖｃを時間軸方向に積分した積分値を算出し、その積分値が、単調増加するかまたは単調減少するかを判定する。すなわち、たとえば血管領域内において逆流が生じていない場合には、上述した積分値は単調増加するか、または単調減少する。これに対し、血管領域内において逆流が生じている場合には、流速ベクトルＶの中心軸方向成分Ｖｃが逆方向に変化するので、中心軸方向成分Ｖｃの積分値は、図４ＩまたはＩＩに示すように、単調増加または単調減少とはならない。
したがって、逆流特定部１３は、中心軸方向成分Ｖｃの積分値が単調増加するかまたは単調減少する場合には、逆流は無いと判定し、中心軸方向成分Ｖｃの積分値が単調増加しないまたは単調減少しない場合には、逆流が有ると判定する。なお、中心軸方向成分Ｖｃの積分値としては、断面Ｓ内の各流速ベクトルＶの中心軸方向成分Ｖｃの積分値の平均値を用いるようにしてもよいし、中央値を用いるようにしてもよい。
なお、血管領域Ｂの中心軸Ｃに直交する断面Ｓの中心軸方向の設定間隔については、ユーザによって任意の間隔に設定可能である。また、積分値を算出する際の積分期間についても、ユーザによって任意の間隔に設定可能である。本実施形態のように、心臓の血管内を流れる血流の逆流の有無を特定する場合には、積分期間を１心拍の期間とすることが好ましい。
図１に戻り、表示制御部１４は、逆流特定部１３によって特定された逆流箇所の情報を表示部３に表示させるものである。具体的には、表示制御部１４は、逆流箇所を血管領域の座標軸上にマッピングしたマップ画像を生成し、表示部３に表示させる。図５は、血管領域Ｂの座標軸上に逆流箇所Ｒをマッピングしたマッピング画像の一例を示す図である。本実施形態においては、上述したように血管領域の中心軸に直交する断面を設定し、断面毎に積分値を算出するので、図５に示すように、逆流箇所Ｒは断面単位で表示される。
また、断面毎に算出される積分値は逆流量を表す。したがって、積分値に応じて表示を変更することによって、逆流量を表す画像とすることができる。この際、全積分期間の積分値に基づいて逆流量を表す画像を生成して表示することによって、静止画として表示するようにしてもよいし、流速ベクトルを時間軸方向に積分する毎に、すなわち積分期間内における各積分過程の積分値毎に、逆流量を表す画像を順次生成して表示することによって、動画として表示するようにしてもよい。
このように逆流量を表す画像を表示させることによって、直観的に逆流量をとらえることができ、診断や治療の際に役立てることができる。
図５に示すマップ画像は、積分値が大きいほど濃度を大きくしたマップ画像の一例である。マップ画像としては、カラー表示することが望ましく、その場合には、積分値に応じて色相を変化させたり、彩度を変化させたりすることが望ましい。
また、マップ画像を３次元医用画像６上に重畳表示させるようにしてもよい。なお、マップ画像を重ね合わせる画像としては、３次元医用画像６に限らず、２次元の断面画像上にマップ画像を重畳表示するようにしてもよい。
図１に戻り、医用画像保管サーバ２は、各種データを保存して管理するコンピュータであり、大容量記憶装置およびデータベース管理用プログラムを備えている。医用画像保管サーバ２は、予め撮影された３次元医用画像６などをネットワーク経由でＣＴ装置またはＭＲＩ装置などの撮影装置から取得し、大容量記憶装置に保存して管理する。
入力部４は、ユーザによる種々の設定入力を受け付けるものであり、キーボードやマウスなどの入力デバイスを備えたものである。入力部４は、たとえば患者の識別情報の設定入力などを受け付けるものである。
次に、本実施形態の血流解析システムの作用について、図６に示すフローチャートを参照しながら説明する。
まず、ユーザによる患者の識別情報などの設定入力に応じて、その患者の医用画像が画像取得部１０によって取得される（Ｓ１０）。本実施形態においては、上述したように心臓の３次元医用画像６が取得される。
そして、画像取得部１０によって取得された心臓の３次元医用画像６は、血管領域抽出部１１に入力され、血管領域抽出部１１は、入力された３次元医用画像６から血管領域を抽出する（Ｓ１２）。
そして、流速ベクトル取得部１２において、たとえば血管領域を用いてＣＦＤにより血流解析が行われ、血管領域内の流速ベクトルが取得される（Ｓ１４）。そして、逆流特定部１３において、血管領域内の流速ベクトルに基づいて、血流の逆流の有無が特定される。具体的には、逆流特定部１３は、上述したように血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、その積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する（Ｓ１６）。
そして、逆流特定部１３によって特定された逆流箇所の情報が表示制御部１４に出力され、表示制御部１４は、入力された情報に基づいてマップ画像を生成して表示部３に表示させる（Ｓ１８）。
上記実施形態の血流解析システムによれば、血管を含む被写体を撮影した血管画像から血管領域を抽出し、血管領域内における血流速度および血流方向を表す流速ベクトルを取得する。そして、血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する。このように流速ベクトルの中心軸方向成分の積分値に基づいて逆流の有無を特定することによって、想定外の箇所に逆流があった場合でも、その逆流箇所をとらえることができ、逆流箇所の見落としリスクを減らすことができる。
また、上記実施形態の血流解析システムにおいては、血管領域の中心軸方向に直交する断面毎に流速ベクトルの積分値を算出するようにしたが、これに限らず、血管領域のボクセル毎または複数のボクセル毎に積分値を算出するようにしてもよい。図７は、複数のボクセル毎に積分値を算出し、逆流箇所Ｒ１，Ｒ２を特定してマップ画像を生成した例を示す図である。複数のボクセル毎の積分値としては、たとえば複数のボクセルのそれぞれについて算出された積分値の平均値を用いるようにしてもよいし、中央値などを用いるようにしてもよい。このようにボクセル毎または複数のボクセル毎に積分値を算出することによって、逆流箇所をより詳細に特定することができる。
また、上記実施形態の血流解析システムにおいては、表示制御部１４が、マップ画像以外の血流を示す指標を表示部３に表示させるようにしてもよい。具体的には、血流を示す指標として、流速ベクトルを表す矢印をマップ画像に重畳表示したり、血流を示す流線または流跡線などをマップ画像に重畳表示したりしてもよい。
また、上記実施形態の血流解析システムにおいては、３次元医用画像６を用いて逆流箇所を特定するようにしたが、必ずしも３次元的な画像でなくてもよく、２次元の断面画像を用いて、上記実施形態と同様にして流速ベクトルを取得し、その血管領域の中心軸方向成分の積分値の変化に基づいて逆流箇所を特定するようにしてもよい。
１ 血流解析装置
２ 医用画像保管サーバ
３ 表示部
４ 入力部
６ ３次元医用画像
１０ 画像取得部
１１ 血管領域抽出部
１２ 流速ベクトル取得部
１３ 逆流特定部
１４ 表示制御部
Ｂ 血管領域
Ｃ 中心軸
Ｒ 逆流箇所
Ｒ１，Ｒ２ 逆流箇所
Ｓ 断面
Ｖ 流速ベクトル
Ｖｃ 流速ベクトルの中心軸方向成分

Claims (16)

  1. 血管を含む被写体を撮影した血管画像から血管領域を抽出する血管領域抽出部と、
    前記血管領域内における血流速度および血流方向を表す流速ベクトルを取得する流速ベクトル取得部と、
    前記血管領域の延伸方向に延びる中心軸を設定し、前記流速ベクトルの前記中心軸の方向成分を時間軸方向に積分した積分値を算出し、前記積分値の変化に基づいて、前記血管領域内の血流の逆流の有無を特定する逆流特定部とを備えた血流解析装置。
  2. 前記逆流特定部が、前記積分値が単調増加または単調減少でない場合に、前記逆流が有ると特定する請求項１記載の血流解析装置。
  3. 前記逆流の箇所を前記血管領域の座標軸上にマッピングしたマップ画像を表示部に表示させる表示制御部を備えた請求項１または２記載の血流解析装置。
  4. 前記表示制御部が、前記流速ベクトルを予め設定された期間だけ積分した前記積分値を用いて前記マップ画像を生成する請求項３記載の血流解析装置。
  5. 前記予め設定された期間が、１心拍の期間である請求項４記載の血流解析装置。
  6. 前記表示制御部が、前記流速ベクトルを時間軸方向に積分する毎に前記マップ画像を時系列に生成して動画として前記表示部に表示させる請求項３記載の血流解析装置。
  7. 前記表示制御部が、前記マップ画像を２次元または３次元の血管画像に重畳表示させる請求項３から６いずれか１項記載の血流解析装置。
  8. 前記表示制御部が、前記積分値に基づいて、前記逆流の量の情報を前記表示部にさせる請求項３から７いずれか１項記載の血流解析装置。
  9. 前記表示制御部が、前記マップ画像以外の血流を示す指標を前記表示部に表示させる請求項３から８いずれか１項記載の血流解析装置。
  10. 前記血流を示す指標が、前記流速ベクトル、流線または流跡線である請求項９記載の血流解析装置。
  11. 前記逆流特定部が、前記血管領域のボクセル毎または複数のボクセル毎に前記積分値を算出する請求項１から１０いずれか１項記載の血流解析装置。
  12. 前記逆流特定部が、前記血管領域の中心軸方向に直交する断面毎に前記積分値を算出する請求項１から１０いずれか１項記載の血流解析装置。
  13. 前記流速ベクトル取得部が、時系列に撮影された３次元の血管画像または任意の時点において撮影された３次元の血管画像に基づいて、前記流速ベクトルを取得する請求項１から１２いずれか１項記載の血流解析装置。
  14. 前記流速ベクトル取得部が、前記血管領域の数値流体力学(Computational Fluid Dynamics)モデル、３次元シネ位相コントラスト磁気共鳴法または３次元の超音波画像を用いて前記流速ベクトルを取得する請求項１から１２いずれか１項記載の血流解析装置。
  15. 血管を含む被写体を撮影した血管画像から血管領域を抽出し、
    前記血管領域内における血流速度および血流方向を表す流速ベクトルを取得し、
    前記血管領域の延伸方向に延びる中心軸を設定し、前記流速ベクトルの前記中心軸の方向成分を時間軸方向に積分した積分値を算出し、前記積分値の変化に基づいて、前記血管領域内の血流の逆流の有無を特定する血流解析方法。
  16. コンピュータを、
    血管を含む被写体を撮影した血管画像から血管領域を抽出する血管領域抽出部と、
    前記血管領域内における血流速度および血流方向を表す流速ベクトルを取得する血流ベクトル取得部と、
    前記血管領域の延伸方向に延びる中心軸を設定し、前記流速ベクトルの前記中心軸の方向成分を時間軸方向に積分した積分値を算出し、前記積分値の変化に基づいて、前記血管領域内の血流の逆流の有無を特定する逆流特定部として機能させる血流解析プログラム。
JP2019503082A 2017-02-28 2018-02-28 血流解析装置および方法並びにプログラム Active JP6690053B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2017037060 2017-02-28
JP2017037060 2017-02-28
PCT/JP2018/007610 WO2018159708A1 (ja) 2017-02-28 2018-02-28 血流解析装置および方法並びにプログラム

Publications (2)

Publication Number Publication Date
JPWO2018159708A1 true JPWO2018159708A1 (ja) 2019-07-25
JP6690053B2 JP6690053B2 (ja) 2020-04-28

Family

ID=63371317

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019503082A Active JP6690053B2 (ja) 2017-02-28 2018-02-28 血流解析装置および方法並びにプログラム

Country Status (3)

Country Link
US (1) US11266322B2 (ja)
JP (1) JP6690053B2 (ja)
WO (1) WO2018159708A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110309849B (zh) * 2019-05-10 2021-08-06 腾讯医疗健康(深圳)有限公司 血管图像处理方法、装置、设备及存储介质
JP2021186008A (ja) 2020-05-26 2021-12-13 キヤノンメディカルシステムズ株式会社 体液解析装置、体液解析装置の制御方法、およびプログラム

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004201730A (ja) * 2002-12-24 2004-07-22 Hitachi Ltd 複数方向の投影映像を用いた３次元形状の生成方法
JP2005533593A (ja) * 2002-07-25 2005-11-10 コーニンクレッカ　フィリップス　エレクトロニクス　エヌ　ヴィ 最適化された視野マップ
JP2010075636A (ja) * 2008-09-29 2010-04-08 Toshiba Corp 超音波診断装置及び超音波画像表示方法
US20180140420A1 (en) * 2015-04-29 2018-05-24 Kephalios S.A.S. An annuloplasty system and a method for monitoring the effectiveness of an annuloplasty treatment

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3517002B2 (ja) 1994-11-30 2004-04-05 株式会社東芝 超音波診断装置
JP5196994B2 (ja) 2007-12-27 2013-05-15 株式会社東芝 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム
JP6162452B2 (ja) 2013-03-28 2017-07-12 東芝メディカルシステムズ株式会社 血流解析装置及び血流解析プログラム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005533593A (ja) * 2002-07-25 2005-11-10 コーニンクレッカ　フィリップス　エレクトロニクス　エヌ　ヴィ 最適化された視野マップ
JP2004201730A (ja) * 2002-12-24 2004-07-22 Hitachi Ltd 複数方向の投影映像を用いた３次元形状の生成方法
JP2010075636A (ja) * 2008-09-29 2010-04-08 Toshiba Corp 超音波診断装置及び超音波画像表示方法
US20180140420A1 (en) * 2015-04-29 2018-05-24 Kephalios S.A.S. An annuloplasty system and a method for monitoring the effectiveness of an annuloplasty treatment

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
WANG,WUGANG ET AL.: "Quantification of Mitral Regurgitation by General Imaging Three-Dimensional Quantification:Feasibili", J AM SOC ECHOCARDIOGR., vol. 第27巻、第3号, JPN6018014370, March 2014 (2014-03-01), pages 268 - 276, ISSN: 0004235580 *
宮本憲行ほか: "僧帽弁逸脱による僧帽弁逆流の定量におけるproximal isovelocity surfece area法の臨床的有用性に関する検", 超音波医学（日本超音波医学会第６３回研究発表会講演論文集), vol. 第20巻, SupplementII, JPN6018014368, 15 October 1993 (1993-10-15), pages 203 - 204, ISSN: 0004235579 *

Also Published As

Publication number Publication date
WO2018159708A1 (ja) 2018-09-07
US20190374114A1 (en) 2019-12-12
US11266322B2 (en) 2022-03-08
JP6690053B2 (ja) 2020-04-28

Similar Documents

Publication Publication Date Title
JP5820117B2 (ja) データセットを選択式及び対話式に処理する方法
JP2014161734A (ja) 医療画像を整合する方法及びその装置
US10888290B2 (en) Medical-image processing apparatus and medical-image diagnostic apparatus
JP6653673B2 (ja) 血流解析装置および方法並びにプログラム
CN108471994B (zh) 移动ffr模拟
EP1542589B1 (en) Display of image data information
JP2015512292A (ja) マルチプルイメージデータループを獲得して分析する方法及びシステム
JP6690053B2 (ja) 血流解析装置および方法並びにプログラム
US20200170519A1 (en) Fluid analysis apparatus, method for operating fluid analysis apparatus, and fluid analysis program
JP6898047B2 (ja) 時変データの定量的評価
JP7059391B2 (ja) 流体解析装置、方法およびプログラム
JP7368500B2 (ja) 流体解析装置、方法およびプログラム
JP6811872B2 (ja) 流体解析装置および流体解析装置の作動方法並びに流体解析プログラム
JP7195458B2 (ja) 流体解析装置、方法およびプログラム
JP6807981B2 (ja) 画像位置合わせ装置および方法並びにプログラム
Al-Ani et al. 3D visualization based on surface estimation techniques
JP2022171345A (ja) 医用画像処理装置、医用画像処理方法及びプログラム
Ogiela et al. Cognitive methods for semantic image analysis in medical imaging applications
Grabner et al. Aorta cross-section calculation and 3D visualization from CT or MRT data using VRML

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190408

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200408

R150 Certificate of patent or registration of utility model

Ref document number: 6690053

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250