JP7467026B2 - 医用情報処理装置、医用情報処理プログラム、医用情報処理システム - Google Patents

医用情報処理装置、医用情報処理プログラム、医用情報処理システム Download PDF

Info

Publication number
JP7467026B2
JP7467026B2 JP2018191548A JP2018191548A JP7467026B2 JP 7467026 B2 JP7467026 B2 JP 7467026B2 JP 2018191548 A JP2018191548 A JP 2018191548A JP 2018191548 A JP2018191548 A JP 2018191548A JP 7467026 B2 JP7467026 B2 JP 7467026B2
Authority
JP
Japan
Prior art keywords
blood flow
myocardial
coronary artery
display
flow rate
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
JP2018191548A
Other languages
English (en)
Other versions
JP2020058591A (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.)
Canon Medical Systems Corp
Original Assignee
Canon Medical Systems 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 Canon Medical Systems Corp filed Critical Canon Medical Systems Corp
Priority to JP2018191548A priority Critical patent/JP7467026B2/ja
Priority to US16/598,294 priority patent/US20200113450A1/en
Publication of JP2020058591A publication Critical patent/JP2020058591A/ja
Priority to JP2023060713A priority patent/JP2023082164A/ja
Application granted granted Critical
Publication of JP7467026B2 publication Critical patent/JP7467026B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/0035—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A61B5/0044—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the heart
    • 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/02028—Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7278—Artificial waveform generation or derivation, e.g. synthesising signals from measured signals
    • 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/743—Displaying an image simultaneously with additional graphical information, e.g. symbols, charts, function plots
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46—Arrangements for interfacing with the operator or the patient
    • A61B6/461—Displaying means of special interest
    • A61B6/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48—Diagnostic techniques
    • A61B6/486—Diagnostic techniques involving generating temporal series of image data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/503—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/507—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • 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
    • G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/60—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46—Arrangements for interfacing with the operator or the patient
    • A61B6/461—Displaying means of special interest
    • A61B6/466—Displaying means of special interest adapted to display 3D data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46—Arrangements for interfacing with the operator or the patient
    • A61B6/467—Arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B6/469—Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/56—Details of data transmission or power supply, e.g. use of slip rings
    • A61B6/563—Details of data transmission or power supply, e.g. use of slip rings involving image data transmission via a network
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20092—Interactive image processing based on input by user
    • G06T2207/20104—Interactive definition of region of interest [ROI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30172—Centreline of tubular or elongated structure

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)
  • Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Cardiology (AREA)
  • Physiology (AREA)
  • Optics & Photonics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Vascular Medicine (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Dentistry (AREA)
  • Hematology (AREA)
  • Artificial Intelligence (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Human Computer Interaction (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Quality & Reliability (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Description

本発明の実施形態は、医用情報処理装置、医用情報処理プログラム、医用情報処理システムに関する。
従来、心筋の虚血の状態を評価する場合に、冠動脈の血流量及び心筋の血流量に加えて、心筋における毛細血管の血管抵抗（微小循環抵抗）を評価することが重要であると知られている。ここで、微小循環抵抗を評価する方法としては、指先などにレーザー光を当てて、指などの毛細血管における血流の評価を行う方法が知られている。
特開２０１５－０３９４４８号公報
本発明が解決しようとする課題は、毛細血管における抵抗を評価することができる指標を提供することである。
実施形態に係る医用情報処理装置は、取得部と、算出部とを備える。取得部は、冠動脈における血管血流量と、前記冠動脈によって血液が供給される心筋領域における心筋血流量とを取得する。算出部は、前記血管血流量と前記心筋血流量とを組み合わせて、前記心筋領域に対して前記血液を供給する毛細血管における毛細血管抵抗量を示す指標を算出する。
図１は、第１の実施形態に係る医用情報処理装置の構成の一例を示す図である。 図２は、第１の実施形態に係る医用情報処理装置による毛細血管抵抗指標の算出の一例を説明するための図である。 図３は、第１の実施形態に係る医用情報処理装置による毛細血管抵抗指標の算出の一例を説明するための図である。 図４は、第１の実施形態に係る医用情報処理装置による毛細血管抵抗指標の算出の一例を説明するための図である。 図５Ａは、第１の実施形態に係る医用情報処理装置による毛細血管抵抗指標の算出の一例を説明するための図である。 図５Ｂは、第１の実施形態に係る医用情報処理装置による毛細血管抵抗指標の算出の一例を説明するための図である。 図６は、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。 図７Ａは、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。 図７Ｂは、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。 図８Ａは、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。 図８Ｂは、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。 図９は、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。 図１０は、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。 図１１は、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。 図１２は、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。 図１３は、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。 図１４は、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。 図１５は、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。 図１６は、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。 図１７は、第１の実施形態に係る治療方針を判定するための判定基準の一例を説明するための図である。 図１８は、第１の実施形態に係る医用情報処理装置による処理手順を示すフローチャートである。
以下に添付図面を参照して、本願に係る医用情報処理装置、医用情報処理プログラム、医用情報処理システムの実施形態を詳細に説明する。なお、本願に係る医用情報処理装置、医用情報処理プログラム、医用情報処理システムは、以下に示す実施形態によって限定されるものではない。
（第１の実施形態）
図１は、第１の実施形態に係る医用情報処理装置３００の構成の一例を示す図である。図１に示すように、第１の実施形態に係る医用情報処理装置３００は、ネットワーク４００を介して、医用画像診断装置１００と、サーバ装置２００に接続される。なお、図１に示す例は、あくまでも一例であり、ネットワーク４００にその他種々の装置（例えば、端末装置等）が接続される場合であってもよい。
医用画像診断装置１００は、Ｘ線診断装置、Ｘ線ＣＴ（Computed Tomography）装置、ＭＲＩ（Magnetic Resonance Imaging）装置、超音波診断装置、ＳＰＥＣＴ（Single Photon Emission Computed Tomography）装置、ＰＥＴ（Positron Emission computed Tomography）装置、ＳＰＥＣＴ装置とＸ線ＣＴ装置とが一体化されたＳＰＥＣＴ－ＣＴ装置、ＰＥＴ装置とＸ線ＣＴ装置とが一体化されたＰＥＴ－ＣＴ装置、又はこれらの装置群等である。また、第１の実施形態に係る医用画像診断装置１００は、３次元の医用画像データ（ボリュームデータ）を生成可能である。
ここで、医用画像診断装置１００は、冠動脈の血流量又は心筋の血流量を定量することが可能な医用画像データを収集する。例えば、医用画像診断装置１００であるＸ線ＣＴ装置は、造影剤を投与した被検体の心臓を略中心にＸ線管及びＸ線検出器を旋回移動させ、被検体を透過したＸ線を検出して投影データを収集する。そして、Ｘ線ＣＴ装置は、収集した投影データに基づいて、時系列の３次元のＣＴ画像データ（ボリュームデータ）を生成する。一例を挙げると、Ｘ線ＣＴ装置は、流体解析によって冠動脈の血流量を算出するための冠動脈造影ＣＴ画像データや、パーフュージョン解析によって心筋の血流量を算出するための心筋造影ＣＴ画像データを収集する。
そして、医用画像診断装置１００は、医用情報処理装置３００からの要求に応じて、収集した医用画像データを医用情報処理装置３００に送信する。また、医用画像診断装置１００は、収集した医用画像データに対して各種解析を行った結果を、医用情報処理装置３００に送信することもできる。
サーバ装置２００は、医用画像診断装置によって収集された医用画像データ（例えば、Ｘ線ＣＴ装置によって収集されたＣＴ画像データ及びＣＴ画像等）や、各種検査情報（例えば、プレッシャーワイヤによって計測された血管内圧力の情報等）等を記憶したり、医用画像データに対して各種画像処理を行ったりする装置である。ここで、サーバ装置２００は、ネットワーク４００を介して医用画像診断装置１００から取得した医用画像データや、各種検査情報等を装置内又は装置外に設けられた記憶回路に記憶させる。そして、サーバ装置２００は、医用情報処理装置３００からの要求に応じて、記憶回路に記憶させた医用画像データや、各種検査情報等を医用情報処理装置３００に送信する。
医用情報処理装置３００は、ネットワーク４００を介して医用画像診断装置１００やサーバ装置２００から医用画像データを取得し、取得した画像データを処理する。また、医用情報処理装置３００は、ネットワーク４００を介してサーバ装置２００から各種検査情報を取得し、取得した検査情報を用いて各種処理を実行する。例えば、医用情報処理装置３００は、ワークステーション等のコンピュータ機器によって実現される。
例えば、図１に示すように、医用情報処理装置３００は、通信インターフェース３１０と、記憶回路３２０と、入力インターフェース３３０と、ディスプレイ３４０と、処理回路３５０とを有する。
通信インターフェース３１０は、処理回路３５０に接続され、ネットワーク４００を介して接続された医用画像診断装置１００又はサーバ装置２００との間で行われる各種データの伝送及び通信を制御する。例えば、通信インターフェース３１０は、ネットワークカードやネットワークアダプタ、ＮＩＣ（Network Interface Controller）等によって実現される。一例を挙げると、通信インターフェース３１０は、医用画像診断装置１００又はサーバ装置２００から医用画像データや、検査情報を受信し、受信した医用画像データや検査情報を処理回路３５０に出力する。
記憶回路３２０は、処理回路３５０に接続され、各種データを記憶する。また、記憶回路３２０は、処理回路３５０の処理に用いられる種々の情報や、処理回路３５０による処理結果、処理回路３５０が読み出して実行することで各種機能を実現するための種々のプログラム等を記憶する。例えば、記憶回路３２０は、ＲＡＭ（Random Access Memory）、フラッシュメモリ等の半導体メモリ素子や、ハードディスク、光ディスク等によって実現される。本実施形態では、記憶回路３２０は、医用画像診断装置１００又はサーバ装置２００から受信した医用画像データや検査情報等を記憶する。
入力インターフェース３３０は、処理回路３５０に接続されており、操作者から受け付けた入力操作を電気信号に変換して処理回路３５０に出力する。なお、本明細書において入力インターフェース３３０は、マウス、キーボードなどの物理的な操作部品を備えるものだけに限られない。例えば、装置とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を制御回路へ出力する電気信号の処理回路も入力インターフェースの例に含まれる。
ディスプレイ３４０は、処理回路３５０に接続され、処理回路３５０から出力される各種情報及び各種画像を表示する。例えば、ディスプレイ３４０は、液晶モニタやＣＲＴ（Cathode Ray Tube）モニタ、タッチパネル等によって実現される。例えば、ディスプレイ３４０は、操作者の指示を受け付けるためのＧＵＩ（Graphical User Interface）や、種々の表示画像、処理回路３５０による種々の処理結果を表示する。
処理回路３５０は、入力インターフェース３３０を介して操作者から受け付けた入力操作に応じて、医用情報処理装置３００が有する各構成要素を制御する。例えば、処理回路３５０は、プロセッサによって実現される。本実施形態では、処理回路３５０は、通信インターフェース３１０から出力される医用画像データや、検査情報等を記憶回路３２０に記憶させる。また、処理回路３５０は、記憶回路３２０から医用画像データや、検査情報を読み出して各種処理を実行し、処理結果をディスプレイ３４０に表示させる。
このような構成のもと、本実施形態に係る医用情報処理装置３００は、毛細血管における抵抗を評価することができる指標を提供することを可能にする。具体的には、医用情報処理装置３００は、冠動脈における血管血流量と心筋における心筋血流量とを用いて、冠動脈から流入した血液を心筋に供給する毛細血管における抵抗（微小循環抵抗）を評価するための毛細血管抵抗指標を算出する。
本実施形態に係る処理回路３５０は、図１に示すように、制御機能３５１と、取得機能３５２と、画像生成機能３５３と、算出機能３５４と、判定機能３５５とを実行する。ここで、制御機能３５１は、表示制御部の一例である。また、取得機能３５２は、取得部の一例である。また、算出機能３５４は、算出部の一例である。
制御機能３５１は、医用情報処理装置３００の全体制御を実行する。具体的には、制御機能３５１は、入力インターフェース３３０を介して入力された各種要求に応じた処理を実行するように制御する。例えば、制御機能３５１は、通信インターフェース３１０を介した医用画像データ等の送受信、記憶回路３２０への情報の格納、ディスプレイ３４０への情報（例えば、表示画像や、解析結果）の表示などを制御する。
取得機能３５２は、冠動脈における血管血流量と、冠動脈によって血液が供給される心筋領域における心筋血流量とを取得する。具体的には、取得機能３５２は、医用画像診断装置１００や、サーバ装置２００から取得された医用画像データや、検査情報に基づいて、冠動脈における血管血流量と、冠動脈によって血液が供給される心筋領域における心筋血流量とを取得する。
例えば、取得機能３５２は、医用画像診断装置１００であるＸ線ＣＴ装置、或いは、サーバ装置２００から取得した冠動脈造影ＣＴ画像データに対して流体解析を実行することにより、冠動脈における血管血流量を取得する。また、取得機能３５２は、流体解析による血管血流量の取得以外にも、プレッシャーワイヤによって計測されたＦＦＲ（Fractional Flow Reserve）や、瞬時ＦＦＲ、冠動脈造影ＭＲ画像データに対する流体解析等を用いることにより、冠動脈における血管血流量を取得することができる。以下、冠動脈造影ＣＴ画像データに対して流体解析を実行することにより、冠動脈における血管血流量を取得する場合を一例に挙げて説明する。
なお、ＦＦＲ及び瞬時ＦＦＲは、血流の阻害の程度を推測する指標であり、血管の岐始部側の血流量と末梢側の血流量との比で定義されるものである。ここで、実際のＦＦＲの計測においては、血管内の血流量と圧力との関係を比例関係にすることで、血流量を圧力に置き換えて計測される。例えば、ＦＦＲの計測においては、アデノシンを投与して最大充血状態（ストレス状態）とすることによって血管内の血流量と圧力との関係を比例関係にし、ＦＦＲを定義する血流量を圧力に置き換える。また、例えば、瞬時ＦＦＲの計測においては、アデノシンを投与せず、安静状態で血管内の血流量と圧力との関係が比例関係となる時相を用いることで、瞬時ＦＦＲを定義する血流量を圧力に置き換える。そこで、取得機能３５２は、プレッシャーワイヤによって計測された圧力の値を血流量の値に変換することで、冠動脈における血管血流量を取得する。
冠動脈造影ＣＴ画像データに基づく流体解析を実行して血管血流量を取得する場合、取得機能３５２は、記憶回路３２０から経時的に収集された複数時相の冠動脈造影ＣＴ画像データを読み出し、読み出した複数時相の冠動脈造影ＣＴ画像データに対して画像処理を行うことで、時系列の血管形状データを抽出する。
ここで、取得機能３５２は、冠動脈造影ＣＴ画像データに含まれる血管領域に指標値を算出する対象領域を設定する。具体的には、取得機能３５２は、操作者による入力インターフェース３３０を介した指示又は画像処理によって、血管領域に対象領域を設定する。なお、対象領域の設定の詳細については、後に詳述する。そして、取得機能３５２は、設定した対象領域の血管形状データとして、例えば、血管の芯線（芯線の座標情報）、芯線に垂直な断面での血管及び内腔の断面積、芯線に垂直な断面での円柱方向の、芯線から内壁までの距離及び芯線から外壁までの距離などを冠動脈造影ＣＴ画像データから抽出する。
さらに、取得機能３５２は、流体解析の解析条件を設定する。具体的には、取得機能３５２は、解析条件として、血液の物性値、反復計算の条件、解析の初期値などを設定する。例えば、取得機能３５２は、血液の物性値として、血液の粘性、密度などを設定する。また、取得機能３５２は、反復計算の条件として、反復計算における最大反復回数、緩和係数、残差の許容値などを設定する。また、取得機能３５２は、解析の初期値として、血流量、圧力、流体抵抗、圧力境界の初期値などを設定する。なお、取得機能３５２によって用いられる各種値は、システムに予め組み込んでおいてもよいし、操作者が対話的に定義してもよい。
そして、取得機能３５２は、冠動脈造影ＣＴ画像データを用いた流体解析により血管の血流に関する指標値を算出する。具体的には、取得機能３５２は、血管形状データと解析条件とを用いた流体解析を実行し、血管の対象領域における血流に関する指標値を算出する。例えば、取得機能３５２は、血管の内腔や外壁の輪郭、血管の断面積及び芯線などの血管形状データと、血液の物性値、反復計算の条件及び解析の初期値などの設定条件に基づいて、血管の所定の位置ごとに、圧力、血流量、血液の流速、ベクトル及びせん断応力などの指標値を算出する。これにより、取得機能３５２は、冠動脈の位置ごとの血流量を取得することができる。
また、例えば、取得機能３５２は、医用画像診断装置１００であるＸ線ＣＴ装置、或いは、サーバ装置２００から取得した心筋造影ＣＴ画像データに対してパーフュージョン解析を実行することにより、心筋における心筋血流量を取得する。なお、取得機能３５２は、心筋造影ＣＴ画像データに基づくパーフュージョン解析による心筋血流量の取得以外にも、ＳＰＥＣＴ装置を用いた心筋シンチグラフィや、ＰＥＴ装置を用いた心筋血流量の定量、心筋造影ＭＲ画像データに対するパーフュージョン解析等を用いることにより、心筋における心筋血流量を取得することができる。以下、心筋造影ＣＴ画像データに対してパーフュージョン解析を実行することにより、心筋における心筋血流量を取得する場合を一例に挙げて説明する。
心筋造影ＣＴ画像データに基づくパーフュージョン解析を実行して心筋血流量を取得する場合、取得機能３５２は、記憶回路３２０から複数心拍分の期間において経時的に収集された複数時相の心筋造影ＣＴ画像データを読み出し、読み出した複数時相の心筋造影ＣＴ画像データに対して画像処理を行うことで心筋血流量を算出する。例えば、取得機能３５２は、複数心拍分の期間において経時的に収集された複数時相の心筋造影ＣＴ画像データに基づいて、左心室内腔（或いは、大動脈）及び心筋に含まれる画素ごとのＴＤＣ（Time Density Curve）を算出する。そして、取得機能３５２は、算出した左心室内腔（或いは、大動脈）におけるＴＤＣ（Time Density Curve）に対する心筋の画素ごとのＴＤＣの対応に基づいて、心筋の画素ごとの心筋血流量を算出する。これにより、取得機能３５２は、心筋の位置ごとの血流量を取得することができる。なお、取得機能３５２は、心筋における所定の領域に含まれる画素群の心筋血流量の平均値を所定の領域における心筋血流量として算出することもできる。
上述したように、取得機能３５２は、ＣＴ画像データを用いて血管血流量及び心筋血流量を算出することで、冠動脈の位置ごとの血管血流量及び心筋の位置ごとの心筋血流量をそれぞれ取得することができる。ここで、取得機能３５２は、操作者によって指定された位置の血管血流量及び心筋血流量を取得する。具体的には、取得機能３５２は、毛細血管抵抗指標を算出する対象領域に対応する冠動脈位置の血管血流量及び心筋領域の心筋血流量を取得する。例えば、取得機能３５２は、狭窄が認められる冠動脈の末端領域や、虚血が認められる心筋領域、或いは、診療の対象となる心筋領域などを、毛細血管抵抗指標を算出する対象領域として、対応する位置の血管血流量及び心筋血流量を取得する。
一例を挙げると、取得機能３５２は、画像生成機能３５３によって生成された表示画像に含まれる冠動脈のうち、操作者によって選択された冠動脈における血管血流量を取得する。そして、取得機能３５２は、選択された冠動脈によって血液が供給される心筋領域を特定し、特定した心筋領域の心筋血流量を取得する。ここで、操作者は、冠動脈において血管血流量を取得する位置を指定することもできる。すなわち、操作者は、血管血流量を取得する冠動脈を選択するとともに、選択した血管のうち血管血流量を算出する位置をさらに指定することできる。
なお、血管血流量の取得する位置が指定されていない場合（冠動脈が選択されただけの場合）、取得機能３５２は、選択された冠動脈の岐始部の血管血流量を算出するようにしてもよい。また、選択された冠動脈の支配領域は、例えば、ボロノイ法によって特定される。すなわち、取得機能３５２は、選択された冠動脈の形状に基づいた領域拡張を行うことにより、選択された冠動脈の支配領域を特定する。
また、例えば、操作者によって心筋領域が指定される場合、取得機能３５２は、まず、操作者によって指定された心筋領域に血液を供給する冠動脈を特定する。一例を挙げると、取得機能３５２は、冠動脈と心筋との形態的な位置関係に基づいて、指定された心筋領域に血液を供給する冠動脈を特定する。ここで、取得機能３５２は、指定された心筋領域に血液を供給する冠動脈の範囲をさらに特定することができる。そして、取得機能３５２は、特定した冠動脈の形状に基づいた領域拡張を行うことにより、選択された冠動脈の支配領域をさらに特定し、特定した支配領域における心筋血流量を取得する。
なお、指定された心筋領域に血液を供給する冠動脈が特定された場合、取得機能３５２は、例えば、特定された冠動脈の岐始部の血管血流量を算出する。また、指定された心筋領域に血液を供給する冠動脈の範囲がさらに特定された場合、取得機能３５２は、例えば、特定された冠動脈の範囲における岐始部側の端部の血管血流量を算出する。
画像生成機能３５３は、記憶回路３２０によって記憶された医用画像データを読み出し、読み出した医用画像データから表示画像を生成する。例えば、画像生成機能３５３は、ＣＴ画像データを読み出し、読み出したＣＴ画像データに対して種々の画像処理を施すことにより、心臓全体や、冠動脈、部分的な心筋領域を示す表示画像を生成する。一例を挙げると、画像生成機能３５３は、ＣＴ画像データに対して画像処理を施すことにより、ボリュームレンダリング画像や、ＣＰＲ（Curved Multi Planer Reconstruction）画像、ＭＰＲ（Multi Planer Reconstruction）画像、ＳＰＲ（Stretched Multi Planer Reconstruction）画像、ポーラーマップ（Polar Map）などを生成する。
また、画像生成機能３５３は、取得機能３５２によって取得された血管血流量や心筋血流量、ＦＦＲ等の検査結果、算出機能３５４によって算出された毛細血管抵抗指標などを用いた表示情報を生成する。例えば、画像生成機能３５３は、心筋血流量を示すパーフュージョン画像や、毛細血管抵抗指標を示すグラフなどを生成する。
算出機能３５４は、血管血流量と心筋血流量とを組み合わせて、心筋領域に対して血液を供給する毛細血管における毛細血管抵抗量を示す毛細血管抵抗指標を算出する。例えば、算出機能３５４は、血管血流量と心筋血流量との比率に基づく毛細血管抵抗指標を算出する。図２は、第１の実施形態に係る医用情報処理装置３００による毛細血管抵抗指標の算出の一例を説明するための図である。
ここで、図２においては、冠動脈の位置Ｐ１における血管血流量と、心筋の領域Ｒ１における心筋血流量が取得された場合について示す。なお、上述したように、位置Ｐ１及び領域Ｒ１は、操作者によって冠動脈が指定されることにより決定される場合であってもよく、或いは、操作者によって心筋領域が指定されることによって決定される場合であってもよい。
算出機能３５４は、冠動脈の位置Ｐ１における血管血流量と、心筋の領域Ｒ１における心筋血流量を用いて、冠動脈の位置Ｐ１及び心筋領域Ｒ１に対応する毛細血管抵抗指標を算出する。例えば、算出機能３５４は、以下の式により毛細血管抵抗指標を算出する。
毛細血管抵抗指標＝心筋血流量（Ｑmyo）／血管血流量（Ｑvessel）
すなわち、算出機能３５４は、冠動脈の血管血流量と、当該冠動脈に支配された心筋領域の心筋血流量との比率を算出することで、冠動脈から心筋領域への血液の供給能（冠動脈と心筋との間にある毛細血管における血液の流れの阻害度合い）を算出することができる。ここで、例えば、毛細血管抵抗が高い場合には冠動脈から心筋領域への血液の供給が減少するため、毛細血管抵抗指標が「１」よりも小さくなる。一方、毛細血管抵抗が低い場合には冠動脈から心筋領域への血液の供給が減少しないため、毛細血管抵抗指標が「１」に近づくこととなる。
例えば、算出機能３５４は、取得機能３５２によって取得された位置Ｐ１における「血管血流量（Ｑvessel）」と心筋領域Ｒ１における「心筋血流量（Ｑmyo）」とから、上記した式に基づいて、位置Ｐ１及び心筋領域Ｒ１に対応する毛細血管抵抗指標を算出する。
上述した実施形態では、単一の冠動脈及び領域を対象に毛細血管抵抗指標を算出する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、医用情報処理装置３００は、複数の冠動脈によって支配される心筋領域を対象に毛細血管抵抗指標を算出することもできる。かかる場合には、例えば、取得機能３５２は、複数本の冠動脈における複数の血管血流量と、複数本の冠動脈によって血液がそれぞれ供給される複数の心筋領域における複数の心筋血流量をまとめた単一の心筋血流量とを取得する。そして、算出機能３５４は、複数の血管血流量と単一の心筋血流量とを組み合わせて毛細血管抵抗指標を算出する。
図３は、第１の実施形態に係る医用情報処理装置３００による毛細血管抵抗指標の算出の一例を説明するための図である。ここで、図３においては、２本の冠動脈によって支配される心筋領域の毛細血管抵抗指標を算出する場合を示す。例えば、図３の例では、毛細血管抵抗指標を算出する対象領域が指定されたことにより、取得機能３５２が、指定された対象領域に対して血液を供給する２本の冠動脈の範囲を特定する。そして、取得機能３５２は、特定した２本の冠動脈の範囲それぞれにおいて、岐始部側の端部の位置Ｐ２の血管血流量「Ｑvessel1」及び位置Ｐ３における血管血流量「Ｑvessel2」を取得する。
さらに、取得機能３５２は、特定した２本の冠動脈の範囲それぞれの支配領域をボロノイ法によってそれぞれ特定し、特定した支配領域をまとめた単一の領域Ｒ２を抽出する。そして、取得機能３５２は、抽出した領域Ｒ２における心筋血流量「Ｑmyo」を取得する。なお、上述した位置Ｐ２、位置Ｐ３及び領域Ｒ２の決定はあくまでも一例であり、例えば、操作者によって位置Ｐ２及び位置Ｐ３が指定され、取得機能３５２が、指定された位置Ｐ２及び位置Ｐ３に基づいて、領域Ｒ２を特定する場合でもよい。
上述したように、取得機能３５２が、位置Ｐ２の血管血流量「Ｑvessel1」及び位置Ｐ３における血管血流量「Ｑvessel2」と、領域Ｒ２における心筋血流量「Ｑmyo」とを取得すると、算出機能３５４は、例えば、「心筋血流量（Ｑmyo）／（血管血流量（Ｑvessel1）＋血管血流量（Ｑvessel2））により、領域Ｒ２における毛細血管抵抗指標を算出する。
上述した例では、心筋領域に血液を供給する範囲を特定する際に、岐始部側の端部のみを特定する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、岐始部側の端部及び末梢側の端部を特定する場合であってもよい。かかる場合には、取得機能３５２は、心筋領域に対して血液を供給する冠動脈の範囲の上流側端部における血管血流量及び下流側端部における血管血流量と、指定された心筋領域における心筋血流量とを取得する。そして、算出機能３５４は、上流側端部における血管血流量と下流側端部における血管血流量との差分と、指定された心筋領域における心筋血流量とを組み合わせて指標を算出する。
図４は、第１の実施形態に係る医用情報処理装置３００による毛細血管抵抗指標の算出の一例を説明するための図である。ここで、図４においては、１本の冠動脈における所定の範囲によって支配される心筋領域の毛細血管抵抗指標を算出する場合を示す。例えば、図４の例では、毛細血管抵抗指標を算出する対象領域が、１本の冠動脈における所定の範囲内に指定されたことにより、取得機能３５２が、指定された対象領域に対して血液を供給する冠動脈の範囲（冠動脈の岐始部側の端部及び末梢側の端部）を特定する。そして、取得機能３５２は、特定した冠動脈の範囲それぞれにおいて、岐始部側の端部の位置Ｐ４の血管血流量「Ｑvessel3」及び位置Ｐ５における血管血流量「Ｑvessel4」を取得する。さらに、取得機能３５２は、位置Ｐ４から位置Ｐ５までの冠動脈における血管血流量として、「Ｑvessel3－Ｑvessel4」を算出する。
そして、取得機能３５２は、特定した冠動脈の範囲に支配された心筋領域Ｒ３をボロノイ法によって特定する。そして、取得機能３５２は、特定した領域Ｒ３における心筋血流量「Ｑmyo」を取得する。なお、上述した位置Ｐ４、位置Ｐ５及び領域Ｒ３の決定はあくまでも一例であり、例えば、操作者によって位置Ｐ４及び位置Ｐ５が指定され、取得機能３５２が、指定された位置Ｐ４及び位置Ｐ５に基づいて、領域Ｒ３を特定する場合でもよい。
上述したように、取得機能３５２が、位置Ｐ４から位置Ｐ５までの冠動脈における血管血流量「Ｑvessel3－Ｑvessel4」と、領域Ｒ３における心筋血流量「Ｑmyo」とを取得すると、算出機能３５４は、例えば、「領域Ｒ３における心筋血流量（Ｑmyo）／（血管血流量（Ｑvessel3）－血管血流量（Ｑvessel4））により、位置Ｐ４から位置Ｐ５までの冠動脈及び領域Ｒ３に対応する毛細血管抵抗指標を算出する。
上述した例では、単一の心筋領域を対象に毛細血管抵抗指標を算出する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、心筋領域を部分領域に分割し、分割した各部分領域における毛細血管抵抗指標を算出する場合であってもよい。かかる場合には、取得機能３５２は、冠動脈及び冠動脈によって血液が供給される心筋領域を複数の範囲に分割し、分割した範囲ごとに、範囲の上流側端部における血管血流量及び下流側端部における血管血流量と、当該範囲に対応する心筋領域における心筋血流量とをそれぞれ取得する。そして、算出機能３５４は、上流側端部における血管血流量と下流側端部における血管血流量との差分と、範囲に対応する心筋領域における心筋血流量とを組み合わせて、分割した範囲ごとに毛細血管抵抗指標をそれぞれ算出する。
図５Ａは、第１の実施形態に係る医用情報処理装置３００による毛細血管抵抗指標の算出の一例を説明するための図である。ここで、図５Ａにおいては、冠動脈によって支配される心筋領域を複数の部分領域に分割し、部分領域ごとの毛細血管抵抗指標を算出する場合を示す。例えば、図５Ａの例では、毛細血管抵抗指標を算出する対象領域が指定されたことにより、取得機能３５２が、指定された対象領域に対して血液を供給する冠動脈の範囲（冠動脈の岐始部側の端部の位置Ｐ６）を特定する。そして、取得機能３５２は、特定した冠動脈の範囲によって支配された心筋領域Ｒ４をボロノイ法によって特定する。
その後、取得機能３５２は、特定した範囲の冠動脈の芯線を取得して、芯線に直交する方向で冠動脈及び心筋領域Ｒ４を複数の範囲に分割する。一例を挙げると、取得機能３５２は、冠動脈の芯線に直交し、位置Ｐ６１を通過する線分によって、領域Ｒ４から部分領域Ｒ４１１及び部分領域Ｒ４１２を分割する。すなわち、取得機能３５２は、位置Ｐ６１を通過する線分によって、位置Ｐ６から位置Ｐ６１までの冠動脈の領域と、部分領域Ｒ４１１及び部分領域Ｒ４１２とを含む範囲を、冠動脈及び領域Ｒ４から分割する。
また、取得機能３５２は、冠動脈の芯線に直交し、位置Ｐ６２を通過する線分によって、領域Ｒ４から部分領域Ｒ４２１及び部分領域Ｒ４２２を分割し、かつ、冠動脈から位置Ｐ６１から位置Ｐ６２までの領域を分割する。同様に、取得機能３５２は、位置Ｐ６３～Ｐ６１０を通過する各線分によって、部分領域Ｒ４３１～Ｒ４９２及び部分領域Ｒ４２を分割する。
そして、取得機能３５２は、分割した範囲それぞれについて、血管血流量と心筋血流量とを取得する。一例を挙げると、取得機能３５２は、位置Ｐ６から位置Ｐ６１までの冠動脈の領域と、部分領域Ｒ４１１及び部分領域Ｒ４１２とを含む範囲について、血管血流量と心筋血流量とを取得する。例えば、取得機能３５２は、位置Ｐ６の血管血流量「Ｑvessel6」及び位置Ｐ６１における血管血流量「Ｑvessel61」を取得する。さらに、取得機能３５２は、位置Ｐ６から位置Ｐ６１までの冠動脈における血管血流量として、「Ｑvessel6－Ｑvessel61」を算出する。また、取得機能３５２は、部分領域Ｒ４１１における心筋血流量「Ｑmyo411」と、部分領域Ｒ４２１における心筋血流量「Ｑmyo421」とを取得する。
同様に、取得機能３５２は、位置Ｐ６１から位置Ｐ６２までの冠動脈の領域と、部分領域Ｒ４２１及び部分領域Ｒ４２２とを含む範囲～位置Ｐ６１０から冠動脈末端までの冠動脈の領域と、部分領域Ｒ４２とを含む範囲までの各範囲について、血管血流量と心筋血流量とをそれぞれ取得する。
算出機能３５４は、取得機能３５２によって取得された範囲ごとの血管血流量及び部分領域における心筋血流量を用いて、部分領域ごとの毛細血管抵抗指標を算出する。例えば、算出機能３５４は、部分領域Ｒ４１１における毛細血管抵抗指標として、「Ｑmyo411／（Ｑvessel6－Ｑvessel61）」を算出する。また、算出機能３５４は、部分領域Ｒ４２１における毛細血管抵抗指標として、「Ｑmyo421／（Ｑvessel6－Ｑvessel61）」を算出する。なお、算出機能３５４は、位置Ｐ６から位置Ｐ６１までの冠動脈に対応する毛細血管抵抗指標として、例えば、部分領域Ｒ４１１における毛細血管抵抗指標と部分領域Ｒ４２１における毛細血管抵抗指標との平均を算出する場合でもよい。
同様に、算出機能３５４は、位置Ｐ６１から位置Ｐ６２までの冠動脈の領域と、部分領域Ｒ４２１及び部分領域Ｒ４２２とを含む範囲～位置Ｐ６１０から冠動脈末端までの冠動脈の領域と、部分領域Ｒ４２とを含む範囲までの各範囲について、部分領域ごとの毛細血管抵抗指標をそれぞれ算出する。なお、位置Ｐ６１０から冠動脈末端までの冠動脈の領域と、部分領域Ｒ４２とを含む範囲については、位置Ｐ６１０から冠動脈末端までの冠動脈の領域における心筋血流量と、部分領域Ｒ４２における心筋血流量とから、１つの毛細血管抵抗指標が算出される。
ここで、図５Ａに示す分割の例はあくまでも一例であり、実施形態はこれに限定されるものではない。すなわち、領域Ｒ４を分割する際の部分領域の数は、図５Ａに図示される数に限定されるものではなく、例えば、図５Ａに示す数以上の部分領域に分割される場合でもよく、或いは、図５Ａに示す数以下の部分領域に分割される場合でもよい。
上述した例では、１本の冠動脈における支配領域を分割する場合について説明した。以下では、冠動脈が分枝する領域の支配領域を分割する場合について説明する。図５Ｂは、第１の実施形態に係る医用情報処理装置３００による毛細血管抵抗指標の算出の一例を説明するための図である。ここで、図５Ｂにおいては、冠動脈が分枝する領域について拡大したものを示す。すなわち、図５Ｂに示す場合にも、実際には、図５Ａと同様に、冠動脈の範囲と、その範囲の冠動脈によって支配される心筋領域が特定されている。
例えば、図５Ｂの例では、取得機能３５２は、まず、位置Ｐ７１から位置Ｐ７２までの冠動脈の領域と、部分領域Ｒ５１１及び部分領域Ｒ５１２とを含む範囲、及び、位置Ｐ７２から位置Ｐ７３までの冠動脈の領域と、部分領域Ｒ５２１及び部分領域Ｒ５２２とを含む範囲などを分割する。そして、取得機能３５２は、分枝が含まれる部分領域Ｒ５２２について、さらに複数の範囲に分割する。例えば、取得機能３５２は、部分領域Ｒ５２２から、位置Ｐ７２１から位置Ｐ７２２までの冠動脈の領域と、部分領域Ｒ５２２１及び部分領域Ｒ５２２２とを含む範囲を分割する。同様に、取得機能３５２は、分枝が含まれる部分領域Ｒ５２２から複数の範囲を分割する。
上述したように分割を行うと、取得機能３５２は、各範囲について、血管血流量と心筋血流量とをそれぞれ取得する。すなわち、取得機能３５２は、分割した範囲ごとに、冠動脈の血管血流量と、部分領域ごとの心筋血流量を取得する。例えば、取得機能３５２は、位置Ｐ７１から位置Ｐ７２までの冠動脈の領域と、部分領域Ｒ５１１及び部分領域Ｒ５１２とを含む範囲については、図５Ａでの説明と同様に、血管血流量と心筋血流量とを取得する。
また、取得機能３５２は、位置Ｐ７２から位置Ｐ７３までの冠動脈の領域と、部分領域Ｒ５２１及び分枝を含む部分領域Ｒ５２２とを含む範囲については、例えば、以下のように、各領域の血流量を取得する。まず、分枝を含まない部分領域Ｒ５２１における毛細血管抵抗指標を算出するための血流量として、取得機能３５２は、図５Ａでの説明と同様に、位置Ｐ７２の血管血流量「Ｑvessel72」及び位置Ｐ７３における血管血流量「Ｑvessel73」を取得する。さらに、取得機能３５２は、位置Ｐ７２から位置Ｐ７３までの冠動脈における血管血流量として、「Ｑvessel72－Ｑvessel73」を算出する。また、取得機能３５２は、部分領域Ｒ５２１における心筋血流量「Ｑmyo521」を取得する。
一方、分枝を含む部分領域Ｒ５２２における毛細血管抵抗指標を算出するための血流量として、取得機能３５２は、位置Ｐ７２の血管血流量「Ｑvessel72」と、分枝位置Ｐ７２１における血管血流量「Ｑvessel721」と、位置Ｐ７３における血管血流量「Ｑvessel73」と、位置Ｐ７２２における血管血流量「Ｑvessel722」と、位置Ｐ７２３における血管血流量「Ｑvessel723」と、位置Ｐ７２４における血管血流量「Ｑvessel724」とを取得する。そして、取得機能３５２は、位置Ｐ７２から位置Ｐ７２１までの冠動脈における血管血流量「Ｑvessel72－Ｑvessel721」、位置Ｐ７２１から位置Ｐ７３までの冠動脈における血管血流量「Ｑvessel721－Ｑvessel73」、位置Ｐ７２１から位置Ｐ７２２までの冠動脈における血管血流量「Ｑvessel721－Ｑvessel722」、位置Ｐ７２２から位置Ｐ７２３までの冠動脈における血管血流量「Ｑvessel722－Ｑvessel723」、位置Ｐ７２３から位置Ｐ７２４までの冠動脈における血管血流量「Ｑvessel723－Ｑvessel724」を算出する。
さらに、取得機能３５２は、部分領域Ｒ５２２１における心筋血流量「Ｑmyo5221」及び部分領域Ｒ５２２１における心筋血流量「Ｑmyo5221」を含む部分領域Ｒ５２２から分割された各部分領域における心筋血流量をそれぞれ取得する。
上述したように、各領域における血流量が取得されると、算出機能３５４は、各血流量を用いて、分割された範囲ごとに、毛細血管抵抗指標を算出する。例えば、算出機能３５４は、位置Ｐ７１から位置Ｐ７２までの冠動脈の領域と、部分領域Ｒ５１１及び部分領域Ｒ５１２とを含む範囲については、図５Ａでの説明と同様に、毛細血管抵抗指標を算出する。
一方、分枝を含む部分領域Ｒ５２２を含む範囲については、算出機能３５４は、例えば、以下のように、毛細血管抵抗指標を算出する。まず、分枝を含まない部分領域Ｒ５２１における毛細血管抵抗指標については、算出機能３５４は、図５Ａでの説明と同様に、「Ｑmyo521／（Ｑvessel72－Ｑvessel73）」を算出する。
そして、分枝を含む部分領域Ｒ５２２に含まれる各部分領域における毛細血管抵抗指標については、算出機能３５４は、例えば、部分領域Ｒ５２２１における毛細血管抵抗指標として、「Ｑmyo5221／（（Ｑvessel72－Ｑvessel721）＋（Ｑvessel721－Ｑvessel722））」を算出する。すなわち、算出機能３５４は、位置Ｐ７２から位置Ｐ７２１までの冠動脈からの血液の供給と、位置Ｐ７２１から位置Ｐ７２２までの冠動脈からの血液の供給とを考慮した毛細血管抵抗指標を算出する。
同様に、算出機能３５４は、例えば、部分領域Ｒ５２２２における毛細血管抵抗指標として、位置Ｐ７２１から位置Ｐ７３までの冠動脈からの血液の供給と、位置Ｐ７２１から位置Ｐ７２２までの冠動脈からの血液の供給とを考慮した毛細血管抵抗指標「Ｑmyo5222／（（Ｑvessel721－Ｑvessel73）＋（Ｑvessel721－Ｑvessel722））」を算出する。なお、それ以降の範囲については、上記と同様に算出される。すなわち、算出機能３５４は、例えば、位置Ｐ７２２から位置Ｐ７２３までの範囲については、両脇の部分領域の心筋血流量に対する、位置Ｐ７２２から位置Ｐ７２３までの冠動脈における血管血流量「Ｑvessel722－Ｑvessel723」をそれぞれ算出することで、両脇の部分領域ごとの毛細血管抵抗指標を算出する。
なお、図５Ｂに示す分割の例はあくまでも一例であり、実施形態はこれに限定されるものではない。すなわち、分割する際の部分領域の数は、図５Ｂに図示される数に限定されるものではなく、例えば、図５Ｂに示す数以上の部分領域に分割される場合でもよく、或いは、図５Ｂに示す数以下の部分領域に分割される場合でもよい。
上述したように、本実施形態に係る取得機能３５２及び算出機能３５４は、冠動脈における血管血流量と心筋における心筋血流量から毛細血管抵抗指標を算出する。ここで、上述した毛細血管抵抗指標の算出は、心臓全体に対して適用することができる。すなわち、取得機能３５２及び算出機能３５４は、心臓全体の心筋と、分枝される全ての血管を含む右冠動脈（Right coronary artery：ＲＣＡ)、左冠動脈前下行枝（Left anterior descending coronary artery：ＬＡＤ)及び左冠動脈回旋枝（Left circumflex coronary artery：ＬＣＸ)とを対象として、毛細血管抵抗指標を算出することができる。
なお、毛細血管抵抗指標を算出する対象領域の設定及び、対象領域を分割するか否か、分割する場合の分割数などは、操作者によって任意に設定することができる。例えば、毛細血管抵抗指標を算出するごとに、任意に設定することができる。また、毛細血管抵抗指標を算出する対象領域ごとに条件が予め設定される場合であってもよい。かかる場合、例えば、心尖部における心筋領域を対象領域として設定された場合、ＲＣＡ、ＬＡＤ及びＬＣＸの３本の血管における血管血流量を用いて毛細血管抵抗指標を算出するように、予め設定されている場合でもよい。
また、上述した血管血流量及び心筋血流量は、ストレス状態及び安静状態のいずれの状態で取得される場合であってもよい。すなわち、本実施形態に係る取得機能３５２及び算出機能３５４は、ストレス状態での毛細血管抵抗指標及び安静状態での毛細血管抵抗指標を算出することができる。
図１に戻って、判定機能３５５は、算出機能３５４によって算出された毛細血管抵抗指標に基づいて、対象領域の状態を判定する。例えば、判定機能３５５は、算出された毛細血管抵抗指標値と閾値とを比較することによって、心筋に対する血液の供給において微小循環が十分に保たれているか否かを判定する。一例を挙げると、判定機能３５５は、毛細血管抵抗指標値が所定の閾値よりも低い場合、毛細血管において血液の供給が阻害されていると判定する。
また、例えば、判定機能３５５は、ストレス状態及び安静状態で算出した毛細血管抵抗指標を比較して、対象領域の状態を判定することもできる。一例を挙げると、判定機能３５５は、虚血が生じている領域におけるストレス状態及び安静状態で算出した毛細血管抵抗指標を比較して両者に違いがない場合には、冠動脈に障害が生じていると判定する。
また、例えば、判定機能３５５は、毛細血管抵抗指標と、他の検査情報とを用いて、治療方針を判定することもできる。なお、治療方針の判定については、後に詳述する。
上述したように、本実施形態に係る医用情報処理装置３００は、心臓全体について、毛細血管抵抗指標を算出することができる。さらに、医用情報処理装置３００は、算出した毛細血管抵抗指標を種々の形態で表示させることができる。すなわち、制御機能３５１が、算出機能３５４によって算出された毛細血管抵抗指標に関する情報をディスプレイ３４０にて表示させることができる。以下、図６～図１６を用いて、本実施形態に係る毛細血管抵抗指標の表示形態の例を説明する。図６～図１６は、第１の実施形態に係る毛細血管抵抗指標の表示形態の例を示す図である。
例えば、制御機能３５１は、心筋領域及び冠動脈の少なくとも一方を含む医用画像上に算出機能３５４によって算出された毛細血管抵抗指標に関する情報を示した表示画像を表示させるように制御する。一例を挙げると、制御機能３５１は、図６に示すように、心臓全体を示すボリュームレンダリング画像上の領域Ｒ１を、算出された毛細血管抵抗指標の値に応じた色で示したカラー画像を表示させる。
かかる場合には、画像生成機能３５３が、例えば、取得されたＣＴ画像データから心臓全体のボリュームレンダリング画像を生成する。制御機能３５１は、生成されたボリュームレンダリング画像の領域Ｒ１を、算出された毛細血管抵抗指標の値に応じた色で示したカラー画像を表示させる。なお、毛細血管抵抗指標に対する色の割り当ては任意に行うことができる。なお、図６に示した例はあくまでも一例であり、表示画像としてその他種々の画像を用いることができる。例えば、サーフェスレンダリング画像が用いられる場合でもよい。
また、制御機能３５１は、分割した範囲ごとにそれぞれ算出された毛細血管抵抗指標を色で示し、心筋領域又は冠動脈を示す３次元画像上に割り当てたカラー画像を表示させるように制御する。例えば、対象領域が部分領域に分割され、部分領域ごとの毛細血管抵抗指標が算出されと、制御機能３５１は、図７Ａに示すように、心筋を示すボリュームレンダリング画像上の領域Ｒ１について、分割された部分領域ごとに算出された毛細血管抵抗指標の値に応じた色で各部分領域を示したカラー画像を表示させる。
また、例えば、対象領域が部分領域に分割され、部分領域ごとの毛細血管抵抗指標が算出されと、制御機能３５１は、図７Ｂに示すように、冠動脈を示すボリュームレンダリング画像上の冠動脈の各範囲に対応する領域について、毛細血管抵抗指標の値に応じた色で各領域を示したカラー画像を表示させる。かかる場合には、画像生成機能３５３が、例えば、取得されたＣＴ画像データから冠動脈のボリュームレンダリング画像を生成する。制御機能３５１は、生成されたボリュームレンダリング画像における冠動脈の各領域を、毛細血管抵抗指標の値に応じた色で示したカラー画像を表示させる。
ここで、複数の部分領域を含む範囲を対象とする場合（例えば、図５Ａにおける部分領域Ｒ４１１及び部分領域Ｒ４１２を含む範囲を対象とする場合）、当該範囲における冠動脈の領域（例えば、図５Ａにおける位置Ｐ１から位置Ｐ６１までの領域）に対して割り当てる毛細血管抵抗指標の値は、各部分領域における毛細血管抵抗指標の平均であってもよく、或いは、部分領域ごとの毛細血管抵抗指標のうち最も低い毛細血管抵抗指標であってもよい。
また、細かい分割によって冠動脈に沿った複数の毛細血管抵抗指標が細かく算出された場合、制御機能３５１は、冠動脈の位置ごとの毛細血管抵抗指標を示すグラフを表示させることができる。例えば、制御機能３５１は、冠動脈を長軸方向に沿って２次元平面上に示した表示画像と、当該冠動脈の各範囲で算出された指標の変化を示すグラフとを、表示画像における位置とグラフにおける位置とを対応付けて表示させるように制御する。
一例を挙げると、制御機能３５１は、図８Ａに示すように、縦軸に毛細血管抵抗指標を示し、横軸に冠動脈の位置を示したグラフを表示させる。さらに、制御機能３５１は、横軸の冠動脈の位置を合わせたＳＰＲ画像をグラフに対応付けて表示させる。かかる場合には、画像生成機能３５３が、例えば、毛細血管抵抗指標が算出された冠動脈のＳＰＲ画像をＣＴ画像データから生成する。さらに、画像生成機能３５３は、生成したＳＰＲ画像の冠動脈における位置と、グラフの横軸における位置とが対応づく縮尺で毛細血管抵抗指標のグラフを生成する。制御機能３５１は、画像生成機能３５３によって生成されたグラフ及びＳＰＲ画像を、図８Ａに示すように対応付けて表示させる。
また、さらに、制御機能３５１は、表示画像に心筋パーフュージョンの結果を表示させるように制御することもできる。例えば、制御機能３５１は、図８Ｂに示すように、ＳＰＲ画像の心筋領域を、取得機能３５２が算出した心筋血流量の結果に対応する色で示したカラー画像を表示させる。
また、医用情報処理装置３００は、心筋全体を表示させるために、ポーラーマップを用いることができる。例えば、制御機能３５１は、心筋を極座標表示させた表示画像上で、冠動脈及び心筋領域を識別可能に表示させ、分割した範囲ごとにそれぞれ算出された指標を色で示したカラー画像を表示させるように制御する。一例を挙げると、制御機能３５１は、図９に示すように、心筋全体を極座標表示させたポーラーマップに対して冠動脈と毛細血管抵抗指標を算出した対象領域とを重畳させた表示画像を表示させる。
かかる場合には、画像生成機能３５３が、ポーラーマップ上の位置に合わせた形状の冠動脈の画像と対象領域の画像とを生成する。そして、制御機能３５１が、生成された冠動脈の画像と対象領域の画像とを、ポーラーマップ上の対応する位置に重畳させて表示させる。ここで、制御機能３５１は、図９に示すように、対象領域における各部分領域を、算出機能３５４によって算出された毛細血管抵抗指標の値に応じた色で示すことで、カラー画像を表示させる。
上述したように、制御機能３５１は、毛細血管抵抗指標を種々の表示形態で表示させることできる。ここで、制御機能３５１は、上述した種々の形態の表示画像を適宜組み合わせて表示させることができる。例えば、制御機能３５１は、図６～図９に示す各表示画像全てを並べて表示させたり、各表示画像から選択された複数の表示画像を並べて表示させたりすることができる。
ここで、制御機能３５１は、複数の表示画像を並べて表示させる場合、画像間での位置関係を示すマーカを表示させることができる。例えば、制御機能３５１は、カラー画像や、グラフ及びＳＰＲ画像、ポーラーマップなどを並べて表示させる際に、それら画像に対して略同一位置を指示するマーカを配置して表示させる。
一例を挙げると、制御機能３５１は、図１０に示すように、縦軸に毛細血管抵抗指標を示し、横軸に冠動脈の位置を示したグラフ及び横軸の冠動脈の位置を合わせたＳＰＲ画像と、冠動脈のボリュームレンダリング画像を用いたカラー画像とを表示させる際に、冠動脈上の略同一位置を示すマーカＭ１及びマーカＭ２を表示させる。
すなわち、制御機能３５１は、画像生成機能３５３によるＳＰＲ画像及びボリュームレンダリング画像の生成時の座標情報に基づいて、グラフ（ＳＰＲ画像）の横軸の位置と、ボリュームレンダリング画像内の冠動脈の位置との位置関係の情報を取得する。そして、制御機能３５１は、取得した位置関係の情報に基づいて、冠動脈上の略同一位置を示すマーカＭ１をグラフ（或いは、ＳＰＲ画像）の横軸に配置し、マーカＭ２をボリュームレンダリング画像における冠動脈に配置する。
ここで、マーカの表示は、複数の表示画像が表示されると同時に開始される場合でもよく、或いは、入力インターフェース３３０を介した操作者によるマーカの表示開始指示に応じて開始される場合でもよい。また、マーカは、入力インターフェース３３０を介した操作者による移動指示に応じて移動される。すなわち、制御機能３５１は、入力インターフェース３３０を介したマーカの移動操作に応じた位置にマーカを移動させて表示させる。このとき、制御機能３５１は、各表示画像に配置された複数のマーカを連動して移動させる。例えば、操作者がマーカＭ１を移動させる操作を実行した場合に、制御機能３５１は、マーカＭ２も連動して移動させる。
さらに、制御機能３５１は、冠動脈において指定された位置の短軸断面画像を表示させることもできる。例えば、制御機能３５１は、冠動脈に対してマーカを配置した際に、配置した位置における冠動脈の短軸断面画像をさらに表示させるように制御する。一例を挙げると、図１０に示すように、冠動脈におけるマーカＭ１が配置されると、まず、画像生成機能３５３が、マーカＭ１が配置された位置のクロスカット画像（芯線に直交する断面画像）を生成する。制御機能３５１は、ＳＰＲ画像（或いは、グラフ）の横軸においてマーカＭ１に対応する位置に、生成された短軸断面画像を表示させる。
上述した例では、毛細血管抵抗指標を種々の表示形態で表示させる場合について説明した。しかしながら、医用情報処理装置３００は、毛細血管抵抗指標以外の指標値をさらに表示させることができる。
例えば、制御機能３５１は、冠動脈を示す３次元画像と、心筋領域を示す３次元画像上に範囲ごとの指標を色で示したカラー画像とを重畳表示させ、冠動脈における位置を指示するマーカと、当該マーカの位置における冠血流予備量比とをさらに表示させるように制御する。一例を挙げると、制御機能３５１は、図１１に示すように、冠動脈を示すボリュームレンダリング画像と毛細血管抵抗指標を算出した領域Ｒ６のボリュームレンダリング画像とを重畳した重畳画像にマーカＭ３を表示させる。そして、制御機能３５１は、マーカＭ３におけるＦＦＲの値「０．８」をさらに表示させる。
すなわち、制御機能３５１は、取得機能３５２による流体解析の結果からマーカＭ３の位置におけるＦＦＲの値を取得して、重畳画像のマーカＭ３の位置に表示させる。ここで、マーカＭ３の位置は、入力インターフェース３３０を介した操作者による移動操作によって移動される。制御機能３５１は、移動操作によってマーカＭ３が移動されるごとに、移動後のマーカＭ３の位置のＦＦＲの値を取得して、重畳画像のマーカＭ３の位置に表示させる。
また、例えば、制御機能３５１は、冠動脈を示す３次元画像と、冠動脈における位置を指示するマーカと、当該マーカの位置における毛細血管抵抗指標及び冠血流予備量比とを表示させるように制御する。一例を挙げると、制御機能３５１は、図１２に示すように、冠動脈を示すボリュームレンダリング画像にマーカＭ３とマーカＭ４とを表示させる。そして、制御機能３５１は、マーカＭ３におけるＦＦＲの値「０．８」とマーカＭ４における毛細血管抵抗指標に値「０．９」をさらに表示させる。
ここで、マーカＭ３及びマーカＭ４の位置は、入力インターフェース３３０を介した操作者による移動操作によって移動される。制御機能３５１は、移動操作によってマーカＭ３が移動されるごとに、移動後のマーカＭ３の位置のＦＦＲの値を取得して、ボリュームレンダリング画像のマーカＭ３の位置に表示させる。また、制御機能３５１は、移動操作によってマーカＭ４が移動されるごとに、移動後のマーカＭ４の位置の毛細血管抵抗指標の値、ボリュームレンダリング画像のマーカＭ４の位置に表示させる。なお、マーカＭ３及びマーカＭ４は、どちらか一方の移動操作に連動して他方も移動するように制御させることが可能である。すなわち、マーカＭ３及びマーカＭ４は、それぞれ別々に移動させることができるとともに、連動して移動させることもできる。
また、例えば、制御機能３５１は、心筋を極座標表示させた表示画像上に心筋パーフュージョンの結果を表示させるとともに、当該表示画像に含まれる心筋領域において毛細血管抵抗指標が閾値より低い領域を識別可能に表示させるように制御する。一例を挙げると、制御機能３５１は、図１３に示すように、心筋パーフュージョンによって得られた心筋血流量の値に応じた色をポーラーマップに割り当て、さらに、毛細血管抵抗指標が閾値より低い領域を識別可能に表示させる。ここで、制御機能３５１は、例えば、閾値より低い領域を枠で囲んだり、別の色でマスクしたりすることによって、識別可能に表示させる。
また、例えば、制御機能３５１は、心筋を示す３次元画像上に心筋パーフュージョンの結果を表示させるとともに、当該３次元画像に含まれる心筋領域において指標が閾値より低い領域を識別可能に表示させるように制御する。一例を挙げると、制御機能３５１は、図１４に示すように、心筋パーフュージョンによって得られた心筋血流量の値に応じた色をポーラーマップに割り当て、さらに、毛細血管抵抗指標が閾値より低い領域を識別可能に表示させる。ここで、制御機能３５１は、例えば、閾値より低い領域を枠で囲んだり、別の色でマスクしたりすることによって、識別可能に表示させる。
また、例えば、制御機能３５１は、冠動脈を示す３次元画像を表示させ、３次元画像に対する冠動脈における位置を指定する指定操作を受け付けた場合に、受け付けた位置における冠血流予備量比、受け付けた位置の冠動脈によって血液が供給される心筋領域における心筋パーフュージョンの結果、受け付けた位置における毛細血管抵抗指標、及び、受け付けた位置の画素値に基づく値を並べて表示させるように制御する。すなわち、制御機能３５１は、冠動脈上で位置が指定された場合に、種々の指標を同時に表示させることができる。
一例を挙げると、制御機能３５１は、図１５に示すように、冠動脈のボリュームレンダリング画像内の冠動脈に対して、冠動脈における位置を指定するためのマーカＭ５を配置して表示する。そして、制御機能３５１は、マーカＭ５の位置における「ＦＦＲ：ａａａ」、「ＭＢＦ：ｂｂｂ」、「毛細血管抵抗指標：ｃｃｃ」及び「ＴＡＧ：ｄｄｄ」を並べて表示させる。
なお、ＭＢＦ（Myocardial Blood Flow）は、心筋血流量であり、マーカＭ５の位置に隣接する心筋における心筋血流量が表示される。また、ＴＡＧ（Transluminal Attenuation Gradient）は、冠動脈における画素値（ＨＵ値）に基づく値であり、冠動脈の芯線に沿った上流側から下流側への距離を横軸にとり、各距離におけるＨＵ値を縦軸にとったグラフにおける傾き（ＨＵ値の空間的な変化率）である。
例えば、造影剤が注入されて冠動脈上の２点間の距離を造影剤が移動する際に血流が遅いほど時間がかかるため、病変によって血流が低下している血管の方が造影効果の差が大きくなり、血流が遅いほどＴＡＧが大きくなる。なお、ＴＡＧは、冠動脈造影ＣＴ画像データを用いて、取得機能３５２によって算出される。
また、例えば、制御機能３５１は、心筋を示す３次元画像を表示させ、３次元画像に対する心筋領域を指定する指定操作を受け付けた場合に、受け付けた心筋領域に対して血液を供給する冠動脈における冠血流予備量比、受け付けた心筋領域における心筋パーフュージョンの結果、及び、受け付けた心筋領域における毛細血管抵抗指標を並べて表示させるように制御する。
一例を挙げると、制御機能３５１は、図１６に示すように、心臓のボリュームレンダリング画像内の心筋に対して毛細血管抵抗指標を算出する領域Ｒ１が指定された場合、領域Ｒ１に血液を供給する冠動脈の位置Ｐ１における「ＦＦＲ：ｅｅｅ」、領域Ｒ１における「ＭＢＦ：ｆｆｆ」及び領域Ｒ１における「毛細血管抵抗指標：ｇｇｇ」を並べて表示させる。
また、例えば、制御機能３５１は、並べて表示させた各表示項目と、表示項目ごとに設定された閾値との比較結果に基づいて、受け付けた位置及び受け付けた心筋領域の少なくとも一方に対する推奨の治療方針を表示させる。例えば、制御機能３５１は、図１５や、図１６に示すように、複数の表示項目（ＦＦＲや、ＭＢＦなど）を並べて表示させた際に、推奨の治療方針を表示させることができる。
かかる場合には、まず、判定機能３５５が、各表示項目をそれぞれ閾値と比較することで、推奨の治療方針を判定する。例えば、判定機能３５５は、図１７に示す判定基準に基づいて、推奨の治療方針を判定する。図１７は、第１の実施形態に係る治療方針を判定するための判定基準の一例を説明するための図である。ここで、図１７においては、ＦＦＲと、ＭＢＦと、毛細血管抵抗指標とを用いた治療方針の判定に用いられる判定基準を示す。また、図１７における上向き矢印（↑）は閾値よりも高いことを示し、下向き矢印（↓）は閾値よりも低いことを示す。
なお、図１７に示すような判定基準は、予め設定されて記憶回路３２０に記憶される。すなわち、判定機能３５５は、判定に際して、記憶回路３２０に記憶された判定基準を参照して判定を行う。また、各表示項目の判定に用いられる各閾値は、操作者などによって予め設定される。
例えば、判定機能３５５は、図１７に示すように、「ＦＦＲ：↓、ＭＢＦ：↓、毛細血管抵抗指標：↑」である場合、「推奨治療」として「カテーテル治療」を判定する。また、判定機能３５５は、図１７に示すように、「ＦＦＲ：↑、ＭＢＦ：↓、毛細血管抵抗指標：↓」である場合、「推奨治療」として「薬物治療」を判定する。また、判定機能３５５は、図１７に示すように、「ＦＦＲ：↓、ＭＢＦ：↑、毛細血管抵抗指標：↑」である場合、「推奨治療」として「治療不要」を判定する。また、判定機能３５５は、図１７に示すように、「ＦＦＲ：↓、ＭＢＦ：↓、毛細血管抵抗指標：↓」である場合、「推奨治療」として「カテーテル治療」及び「薬物治療」を判定する。
制御機能３５１は、判定機能３５５によって判定された結果（推奨治療）を、表示画像における各表示項目にさらに対応付けて表示させる。上述した推奨治療の判定は、冠動脈の位置ごと、或いは、心筋領域ごとに実行することができる。すなわち、操作者が、入力インターフェース３３０を用いてマーカＭ５を移動させたり、毛細血管抵抗指標を算出する対象領域を指定したりするごとに、判定機能３５５が、移動後の位置ごとに推奨の治療方針を判定し、制御機能３５１が、判定結果を表示させる。
なお、図１７に示す判定基準はあくまでも一例であり、実施形態はこれに限定されるものではない。例えば、閾値と比較する項目としてＴＡＧや、その他の指標が用いられる場合でもよい。
次に、第１の実施形態に係る医用情報処理装置３００による処理の手順について説明する。図１８は、第１の実施形態に係る医用情報処理装置３００による処理手順を示すフローチャートである。ここで、図１８におけるステップＳ１０１、ステップＳ１１１は、例えば、処理回路３５０が制御機能３５１に対応するプログラムを記憶回路３２０から呼び出して実行することにより実現される。また、ステップＳ１０３～ステップＳ１０８は、例えば、処理回路３５０が取得機能３５２に対応するプログラムを記憶回路３２０から呼び出して実行することにより実現される。また、ステップＳ１０９は、例えば、処理回路３５０が算出機能３５４に対応するプログラムを記憶回路３２０から呼び出して実行することにより実現される。また、ステップＳ１０２、ステップＳ１１０は、例えば、処理回路３５０が制御機能３５１及び画像生成機能３５３に対応するプログラムを記憶回路３２０から呼び出して実行することにより実現される。
本実施形態に係る医用情報処理装置３００では、まず、処理回路３５０が、ＣＴ画像データを取得する（ステップＳ１０１）。そして、処理回路３５０が、ＣＴ画像データから心臓の画像を生成して表示する（ステップＳ１０２）。そして、処理回路３５０は、冠動脈に対する位置の指定を受け付けたか否かを判定する（ステップＳ１０３）。ここで、冠動脈に対する位置の指定を受け付けた場合（ステップＳ１０３肯定）、処理回路３５０は、指定された冠動脈が血液を供給する心筋領域を特定し（ステップＳ１０４）、指定された冠動脈における血管血流量及び特定した心筋領域における心筋血流量を算出する（ステップＳ１０５）。
一方、ステップＳ１０３において、冠動脈に対する位置の指定を受けていない場合（ステップＳ１０３否定）、処理回路３５０は、心筋に対する領域の指定を受け付けたか否かを判定する（ステップＳ１０６）。ここで、心筋に対する領域の指定を受け付けた場合（ステップＳ１０６肯定）、処理回路３５０は、指定された心筋領域に対して血液を供給する冠動脈を特定し（ステップＳ１０７）、指定された心筋領域における心筋血流量及び特定した冠動脈における血管血流量を算出する（ステップＳ１０８）。なお、ステップＳ１０６において、心筋に対する領域の指定を受け付けていない場合（ステップＳ１０６否定）、処理回路３５０は、ステップＳ１０３に戻って、判定を継続する。
ステップＳ１０５又はステップＳ１０８において、血管血流量及び心筋血流量を算出すると、処理回路３５０は、血管血流量及び心筋血流量を用いて毛細血管抵抗指標を算出する（ステップＳ１０９）。そして、処理回路３５０は、算出した毛細血管抵抗指標に関する情報を生成して、表示する（ステップＳ１１０）。その後、処理回路３５０は、毛細血管抵抗指標の算出が終了されたか否かを判定して（ステップＳ１１１）、終了された場合に（ステップＳ１１１肯定）、処理を終了する。一方、終了されていない場合（ステップＳ１１１否定）、処理回路３５０は、ステップＳ１０３に戻って、判定を継続する。
上述したように、第１の実施形態によれば、取得機能３５２は、冠動脈における血管血流量と、冠動脈によって血液が供給される心筋領域における心筋血流量とを取得する。算出機能３５４は、血管血流量と心筋血流量とを組み合わせて、心筋領域に対して血液を供給する毛細血管における毛細血管抵抗量を示す指標を算出する。従って、第１の実施形態に係る医用情報処理装置３００は、冠動脈から心筋領域への血液の供給能（冠動脈と心筋との間にある毛細血管における血液の流れの阻害度合い）を算出することができ、毛細血管における抵抗を評価することができる指標を提供することを可能にする。
例えば、心筋虚血を診断するための指標として、ＣＦＲ（Coronary Flow Reserve）が知られている。本願にかかる毛細血管抵抗指標は、臨床的にこのＣＦＲと似た意味合いをもつものである。しかしながら、ＣＦＲは、ストレス状態の心筋血流量と安静状態の心筋血流量との比であるため、ＣＴ画像データからこれらを算出するためには、ストレス状態と安静状態で２回の撮影を行うこととなり、被ばく量が多くなる。また、ＣＦＲを算出するためには、血管をストレス状態とするために血管拡張剤を用いるため、被検体に対して負担がかかる。
これに対して、本願にかかる毛細血管抵抗指標は、複数心拍分の複数時相のＣＴ画像データを取得することで、算出させることができ、被ばく量を抑制することができる。また、本願にかかる毛細血管抵抗指標は、安静状態で算出することができるため、被検体に対して薬物投与の負担をかけずに済む。すなわち、本願に係る毛細血管抵抗指標は、血管拡張剤を投与することができない被検体でも算出することができる。
また、本願にかかる毛細血管抵抗指標は、種々の手法によって算出することができる冠動脈の血管血流量と心筋血流量とを用いるため、様々な被検体において算出させることができる。
また、第１の実施形態によれば、取得機能３５２は、複数本の冠動脈における複数の血管血流量と、複数本の冠動脈によって血液がそれぞれ供給される複数の心筋領域における複数の心筋血流量をまとめた単一の心筋血流量とを取得する。算出機能３５４は、複数の血管血流量と単一の心筋血流量とを組み合わせて毛細血管抵抗指標を算出する。従って、第１の実施形態に係る医用情報処理装置３００は、冠動脈から心筋への血液の供給をより正確に反映させた毛細血管抵抗指標を提供することを可能にする。
また、第１の実施形態によれば、取得機能３５２は、心筋領域に対して血液を供給する冠動脈の範囲の上流側端部における血管血流量及び下流側端部における血管血流量と、指定された心筋領域における心筋血流量とを取得する。算出機能３５４は、上流側端部における血管血流量と下流側端部における血管血流量との差分と、指定された心筋領域における心筋血流量とを組み合わせて毛細血管抵抗指標を算出する。従って、第１の実施形態に係る医用情報処理装置３００は、より精度の高い毛細血管抵抗指標を提供することを可能にする。
また、第１の実施形態によれば、取得機能３５２は、冠動脈及び冠動脈によって血液が供給される心筋領域を複数の範囲に分割し、分割した範囲ごとに、範囲の上流側端部における血管血流量及び下流側端部における血管血流量と、当該範囲に対応する心筋領域における心筋血流量とをそれぞれ取得する。算出機能３５４は、上流側端部における血管血流量と下流側端部における血管血流量との差分と、範囲に対応する心筋領域における心筋血流量とを組み合わせて、分割した範囲ごとに毛細血管抵抗指標をそれぞれ算出する。従って、第１の実施形態に係る医用情報処理装置３００は、より細かい領域について毛細血管抵抗指標を提供することを可能にする。
また、第１の実施形態によれば、制御機能３５１は、心筋領域及び冠動脈の少なくとも一方を含む医用画像上に算出機能３５４によって算出された毛細血管抵抗指標に関する情報を示した表示画像を表示させるように制御する。従って、第１の実施形態に係る医用情報処理装置３００は、毛細血管抵抗指標について、よりわかりやすい情報を提供することを可能にする。
また、第１の実施形態によれば、制御機能３５１は、分割した範囲ごとにそれぞれ算出された毛細血管抵抗指標を色で示し、心筋領域又は冠動脈を示す３次元画像上に割り当てたカラー画像を表示させるように制御する。従って、第１の実施形態に係る医用情報処理装置３００は、各領域における毛細血管抵抗指標について、よりわかりやすい情報を提供することを可能にする。
また、第１の実施形態によれば、制御機能３５１は、冠動脈を長軸方向に沿って２次元平面上に示した表示画像と、当該冠動脈の各範囲で算出された毛細血管抵抗指標の変化を示すグラフとを、表示画像における位置とグラフにおける位置とを対応付けて表示させるように制御する。従って、第１の実施形態に係る医用情報処理装置３００は、心筋の性状と毛細血管抵抗指標との関係をよりわかりやすく表示させることを可能にする。
また、第１の実施形態によれば、制御機能３５１は、表示画像に心筋パーフュージョンの結果を表示させるように制御する。従って、第１の実施形態に係る医用情報処理装置３００は、心筋血流量と毛細血管抵抗指標との関係について、よりわかりやすい情報を提供することを可能にする。
また、第１の実施形態によれば、制御機能３５１は、心筋を極座標表示させた表示画像上で、冠動脈及び心筋領域を識別可能に表示させ、分割した範囲ごとにそれぞれ算出された指標を色で示したカラー画像を表示させるように制御する。従って、第１の実施形態に係る医用情報処理装置３００は、心筋全体を把握させながら、対象領域の毛細血管抵抗指標について、よりわかりやすい情報を提供することを可能にする。
また、第１の実施形態によれば、制御機能３５１は、分割した範囲ごとにそれぞれ算出された毛細血管抵抗指標を色で示し、心筋領域又は冠動脈を示す３次元画像上に割り当てた第１のカラー画像、冠動脈を長軸方向に沿って２次元平面上に示した表示画像及び当該冠動脈の各範囲で算出された毛細血管抵抗指標の変化を示すグラフ、及び、心筋を極座標表示させた表示画像上で、冠動脈及び心筋領域を識別可能に表示させ、分割した範囲ごとにそれぞれ算出された毛細血管抵抗指標を色で示した第２のカラー画像のうち少なくとも２つを表示させるように制御する。さらに、制御機能３５１は、表示させた第１のカラー画像、表示画像及びグラフ、及び、第２のカラー画像のうち少なくとも２つに対して、略同一位置を指示するマーカを配置して表示させるように制御する。従って、第１の実施形態に係る医用情報処理装置３００は、複数の表示画像を比較して診断することができるとともに、複数の画像間において位置関係を明瞭にすることを可能にする。また、心筋と冠動脈の形態と毛細血管抵抗指標とを比較して観察することができるため、虚血と心筋梗塞とを見極めて治療方針を決定することを可能にする。
また、第１の実施形態によれば、制御機能３５１は、冠動脈に対してマーカを配置した際に、配置した位置における冠動脈の短軸断面画像をさらに表示させるように制御する。従って、第１の実施形態に係る医用情報処理装置３００は、より詳細な形態情報を提供することを可能にする。
また、第１の実施形態によれば、制御機能３５１は、冠動脈を示す３次元画像と、心筋領域を示す３次元画像上に範囲ごとの毛細血管抵抗指標を色で示したカラー画像とを重畳表示させ、冠動脈における位置を指示するマーカと、当該マーカの位置における冠血流予備量比とをさらに表示させるように制御する。従って、第１の実施形態に係る医用情報処理装置３００は、毛細血管抵抗指標とＦＦＲの値とを比較して観察することを可能にする。
また、第１の実施形態によれば、制御機能３５１は、冠動脈を示す３次元画像と、冠動脈における位置を指示するマーカと、当該マーカの位置における毛細血管抵抗指標及び冠血流予備量比とを表示させるように制御する。従って、第１の実施形態に係る医用情報処理装置３００は、毛細血管抵抗指標とＦＦＲの値とを比較して観察することを可能にする。
また、第１の実施形態によれば、制御機能３５１は、心筋を極座標表示させた表示画像上に心筋パーフュージョンの結果を表示させるとともに、当該表示画像に含まれる心筋領域において毛細血管抵抗指標が閾値より低い領域を識別可能に表示させるように制御する。また、制御機能３５１は、心筋を示す３次元画像上に心筋パーフュージョンの結果を表示させるとともに、当該３次元画像に含まれる心筋領域において毛細血管抵抗指標が閾値より低い領域を識別可能に表示させるように制御する。従って、第１の実施形態に係る医用情報処理装置３００は、心筋血流量の結果と、毛細血管抵抗指標とを比較観察できるため、再灌流によって復活する心筋領域を特定して治療方針を決定することを可能にする。
また、第１の実施形態によれば、制御機能３５１は、冠動脈を示す３次元画像を表示させ、３次元画像に対する冠動脈における位置を指定する指定操作を受け付けた場合に、受け付けた位置における冠血流予備量比、受け付けた位置の冠動脈によって血液が供給される心筋領域における心筋パーフュージョンの結果、受け付けた位置における毛細血管抵抗指標、及び、受け付けた位置の画素値に基づく値を並べて表示させるように制御する。また、制御機能３５１は、心筋を示す３次元画像を表示させ、３次元画像に対する心筋領域を指定する指定操作を受け付けた場合に、受け付けた心筋領域に対して血液を供給する冠動脈における冠血流予備量比、受け付けた心筋領域における心筋パーフュージョンの結果、及び、受け付けた心筋領域における指標を並べて表示させるように制御する。従って、第１の実施形態に係る医用情報処理装置３００は、特定の指標だけでは判断することが難しい虚血の原因を特定し、治療方針を決定することを可能にする。
また、第１の実施形態によれば、制御機能３５１は、並べて表示させた各表示項目と、表示項目ごとに設定された閾値との比較結果に基づいて、受け付けた位置及び受け付けた心筋領域の少なくとも一方に対する推奨の治療方針を表示させる。従って、第１の実施形態に係る医用情報処理装置３００は、治療方針の決定を支援することを可能にする。
また、第１の実施形態によれば、算出機能３５４は、血管血流量と心筋血流量との比率に基づく毛細血管抵抗指標を算出する。従って、第１の実施形態に係る医用情報処理装置３００は、毛細血管の抵抗を評価するための指標を容易に提供することを可能にする。
（第２の実施形態）
さて、これまで第１の実施形態について説明したが、上述した第１の実施形態以外にも、種々の異なる形態にて実施されてよいものである。
上述した実施形態では、血管血流量と心筋血流量との比率に基づいて毛細血管抵抗指標を算出する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、対数比をとる、分母や分子に係数を加算・乗算するなど、２つの値を使用して算出するものであればよい。
一例を挙げると、血液の流速や、血液の粘性、心臓全体の駆出率・駆出量などを補正係数として掛けあわせてもよい。かかる場合には、算出機能３５４は、取得機能３５２による流体解析の結果から流速や粘性の情報を取得する。また、算出機能３５４は、１心拍以上における複数時相で収集された経時的なＣＴ画像データに基づいて、心臓全体の駆出率・駆出量を算出する。
ここで、算出機能３５４は、血液の流速を用いる場合、例えば、流速が速い場合に毛細血管抵抗指標が低くなるように、血液の流速を補正係数として用いる。また、算出機能３５４は、血液の粘性を用いる場合、例えば、粘性が高い場合に毛細血管抵抗指標が低くなるように、血液の粘性を補正係数として用いる。ここで、補正係数として血液の粘性を用いた場合、例えば、生活習慣病の改善などによって血液の粘性が低くなった場合の心筋血流量の変化も計算することができる。また、算出機能３５４は、心臓全体の駆出率・駆出量を用いる場合、例えば、心臓全体の駆出率・駆出量の変化に伴って毛細血管抵抗指標が規則的に変化するように、心臓全体の駆出率・駆出量を補正係数として用いる。
また、上述した実施形態では、毛細血管抵抗指標をそのまま用いた情報を表示する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、毛細血管抵抗指標に対して他の係数などを掛けた値を用いる場合でもよい。一例を挙げると、算出機能３５４は、算出した毛細血管抵抗指標に対して「細胞間の血液の伝わりやすさ」を示す指標値を掛けた値を用いて、上述した種々の表示情報を表示してもよい。これにより、医用情報処理装置３００は、心筋梗塞領域とその進行具合を示す情報を提供することを可能にする。
また、上述した実施形態では、ＳＰＲ画像やポーラーマップに対して心筋パーフュージョンの結果を用いる場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、代わりに心筋の機能指標を示す他の値を用いる場合でもよい。例えば、心臓を超音波撮像することによって心筋の局所的な動きを抽出し、動きの大きさを心筋パーフュージョンの結果の代わりに用いる場合でもよい。
また、本願に係る医用情報処理装置３００は、算出した血管血流量及び心筋血流量を用いて機械学習を行うこともできる。例えば、医用情報処理装置３００は、血管血流量及び心筋血流量と、対象領域におけるＣＦＲの値との関係を機械学習することにより、識別器を構築する。構築された識別器は、新たな血管血流量及び心筋血流量の入力に対して、ＣＦＲの値を出力する。これにより、血管血流量及び心筋血流量からＣＦＲを推測することを可能にする。すなわち、被ばく量の増大や、血管拡張剤の投与を回避しつつ、ＣＦＲの値を取得することを可能にする。なお、上述した機械学習では、ＣＦＲの代わりに、毛細血管抵抗試験などによって計測された値が用いられる場合でもよい。
また、上述した実施形態では、医用情報処理装置３００単体で各種処理を実行する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、処理回路３５０は、ネットワークを介して接続された外部装置のプロセッサを利用して、機能を実現することとしてもよい。例えば、処理回路３５０は、記憶回路３２０から各機能に対応するプログラムを読み出して実行するとともに、医用情報処理装置３００とネットワークを介して接続されたサーバ群（クラウド）を計算資源として利用することにより、図１に示す各機能を実現する。また、例えば、記憶回路３２０は、医用情報処理装置３００とネットワークを介して接続されたサーバ群（クラウド）により実現されることとしてもよい。
上述した実施形態では、単一の処理回路（処理回路３５０）によって各処理機能が実現される場合の例を説明したが、実施形態はこれに限られない。例えば、処理回路３５０は、複数の独立したプロセッサを組み合わせて構成され、各プロセッサが各プログラムを実行することにより各処理機能を実現するものとしても構わない。また、処理回路３５０が有する各処理機能は、単一又は複数の処理回路に適宜に分散又は統合されて実現されてもよい。
また、上述した各実施形態の説明で用いた「プロセッサ」という文言は、例えば、ＣＰＵ（Central Processing Unit）、ＧＰＵ(Graphics Processing Unit)、或いは、特定用途向け集積回路（Application Specific Integrated Circuit：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Simple Programmable Logic Device：ＳＰＬＤ）、複合プログラマブル論理デバイス（Complex Programmable Logic Device：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Field Programmable Gate Array：ＦＰＧＡ））等の回路を意味する。ここで、記憶回路にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むように構成しても構わない。この場合には、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。また、本実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて一つのプロセッサとして構成され、その機能を実現するようにしてもよい。
ここで、プロセッサによって実行されるプログラムは、ＲＯＭ（Read Only Memory）や記憶回路等に予め組み込まれて提供される。なお、このプログラムは、これらの装置にインストール可能な形式又は実行可能な形式のファイルでＣＤ（Compact Disk）－ＲＯＭ、ＦＤ（Flexible Disk）、ＣＤ－Ｒ（Recordable）、ＤＶＤ（Digital Versatile Disk）等のコンピュータで読み取り可能な記憶媒体に記録されて提供されてもよい。また、このプログラムは、インターネット等のネットワークに接続されたコンピュータ上に格納され、ネットワーク経由でダウンロードされることにより提供又は配布されてもよい。例えば、このプログラムは、後述する各機能部を含むモジュールで構成される。実際のハードウェアとしては、ＣＰＵが、ＲＯＭ等の記憶媒体からプログラムを読み出して実行することにより、各モジュールが主記憶装置上にロードされて、主記憶装置上に生成される。
以上説明した少なくとも一つの実施形態によれば、毛細血管における抵抗を評価することができる指標を提供することができる。
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
３００ 医用情報処理装置
３５１ 制御機能
３５２ 取得機能
３５３ 画像生成機能
３５４ 算出機能
３５５ 判定機能

Claims (21)

  1. 最大充血状態又は安静状態での冠動脈における血管血流量を流体解析によって取得し、前記血管血流量の取得時と同一の状態での前記冠動脈によって血液が供給される心筋領域における心筋血流量をパーフュージョン解析によって取得する取得部と、
    前記血管血流量と前記心筋血流量との比較に基づいて、前記心筋領域に対して前記血液を供給する毛細血管における毛細血管抵抗量を示す指標を算出する算出部と、
    を備える、医用情報処理装置。
  2. 前記取得部は、複数本の冠動脈における複数の血管血流量と、前記複数本の冠動脈によって血液がそれぞれ供給される複数の心筋領域における複数の心筋血流量をまとめた単一の心筋血流量とを取得し、
    前記算出部は、前記複数の血管血流量と前記単一の心筋血流量とを組み合わせて前記指標を算出する、請求項１に記載の医用情報処理装置。
  3. 前記取得部は、心筋領域に対して血液を供給する冠動脈の範囲の上流側端部における血管血流量及び下流側端部における血管血流量と、前記心筋領域における心筋血流量とを取得し、
    前記算出部は、前記上流側端部における血管血流量と前記下流側端部における血管血流量との差分と、前記心筋領域における心筋血流量とを組み合わせて前記指標を算出する、請求項１に記載の医用情報処理装置。
  4. 前記取得部は、前記冠動脈及び前記冠動脈によって血液が供給される心筋領域を複数の範囲に分割し、分割した範囲ごとに、範囲の上流側端部における血管血流量及び下流側端部における血管血流量と、当該範囲に対応する心筋領域における心筋血流量とをそれぞれ取得し、
    前記算出部は、前記上流側端部における血管血流量と前記下流側端部における血管血流量との差分と、前記範囲に対応する心筋領域における心筋血流量とを組み合わせて、前記分割した範囲ごとに前記指標をそれぞれ算出する、請求項１に記載の医用情報処理装置。
  5. 前記心筋領域及び前記冠動脈の少なくとも一方を含む医用画像上に前記算出部によって算出された前記指標に関する情報を示した表示画像を表示させるように制御する表示制御部をさらに備える、請求項１～４のいずれか１つに記載の医用情報処理装置。
  6. 前記分割した範囲ごとにそれぞれ算出された前記指標を色で示し、前記心筋領域又は前記冠動脈を示す３次元画像上に割り当てたカラー画像を表示させるように制御する表示制御部をさらに備える、請求項４に記載の医用情報処理装置。
  7. 前記冠動脈を長軸方向に沿って２次元平面上に示した表示画像と、当該冠動脈の各範囲で算出された前記指標の変化を示すグラフとを、前記表示画像における位置と前記グラフにおける位置とを対応付けて表示させるように制御する表示制御部をさらに備える、請求項４に記載の医用情報処理装置。
  8. 前記表示制御部は、前記表示画像に心筋パーフュージョンの結果を表示させるように制御する、請求項７に記載の医用情報処理装置。
  9. 心筋を極座標表示させた表示画像上で、前記冠動脈及び前記心筋領域を識別可能に表示させ、前記分割した範囲ごとにそれぞれ算出された前記指標を色で示したカラー画像を表示させるように制御する表示制御部をさらに備える、請求項４に記載の医用情報処理装置。
  10. 前記分割した範囲ごとにそれぞれ算出された前記指標を色で示し、前記心筋領域又は前記冠動脈を示す３次元画像上に割り当てた第１のカラー画像、前記冠動脈を長軸方向に沿って２次元平面上に示した表示画像及び当該冠動脈の各範囲で算出された前記指標の変化を示すグラフ、及び、心筋を極座標表示させた表示画像上で、前記冠動脈及び前記心筋領域を識別可能に表示させ、前記分割した範囲ごとにそれぞれ算出された前記指標を色で示した第２のカラー画像のうち少なくとも２つを表示させるように制御する表示制御部をさらに備え、
    前記表示制御部は、表示させた前記第１のカラー画像、前記表示画像及び前記グラフ、及び、前記第２のカラー画像のうち少なくとも２つに対して、略同一位置を指示するマーカを配置して表示させるように制御する、請求項４に記載の医用情報処理装置。
  11. 前記表示制御部は、前記冠動脈に対して前記マーカを配置した際に、配置した位置における冠動脈の短軸断面画像をさらに表示させるように制御する、請求項１０に記載の医用情報処理装置。
  12. 前記冠動脈を示す３次元画像と、前記心筋領域を示す３次元画像上に前記範囲ごとの指標を色で示したカラー画像とを重畳表示させ、前記冠動脈における位置を指示するマーカと、当該マーカの位置における冠血流予備量比とをさらに表示させるように制御する表示制御部をさらに備える、請求項４に記載の医用情報処理装置。
  13. 前記冠動脈を示す３次元画像と、前記冠動脈における位置を指示するマーカと、当該マーカの位置における前記指標及び冠血流予備量比とを表示させるように制御する表示制御部をさらに備える、請求項４に記載の医用情報処理装置。
  14. 心筋を極座標表示させた表示画像上に心筋パーフュージョンの結果を表示させるとともに、当該表示画像に含まれる前記心筋領域において前記指標が閾値より低い領域を識別可能に表示させるように制御する表示制御部をさらに備える、請求項４に記載の医用情報処理装置。
  15. 心筋を示す３次元画像上に心筋パーフュージョンの結果を表示させるとともに、当該３次元画像に含まれる前記心筋領域において前記指標が閾値より低い領域を識別可能に表示させるように制御する表示制御部をさらに備える、請求項４に記載の医用情報処理装置。
  16. 前記冠動脈を示す３次元画像を表示させ、前記３次元画像に対する前記冠動脈における位置を指定する指定操作を受け付けた場合に、受け付けた位置における冠血流予備量比、前記受け付けた位置の冠動脈によって血液が供給される心筋領域における心筋パーフュージョンの結果、前記受け付けた位置における前記指標、及び、前記受け付けた位置の画素値に基づく値を並べて表示させるように制御する表示制御部をさらに備える、請求項４に記載の医用情報処理装置。
  17. 心筋を示す３次元画像を表示させ、前記３次元画像に対する心筋領域を指定する指定操作を受け付けた場合に、受け付けた心筋領域に対して血液を供給する冠動脈における冠血流予備量比、前記受け付けた心筋領域における心筋パーフュージョンの結果、及び、前記受け付けた心筋領域における前記指標を並べて表示させるように制御する表示制御部をさらに備える、請求項４に記載の医用情報処理装置。
  18. 前記表示制御部は、並べて表示させた各表示項目と、表示項目ごとに設定された閾値との比較結果に基づいて、受け付けた位置及び受け付けた心筋領域の少なくとも一方に対する推奨の治療方針を表示させる、請求項１６又は１７に記載の医用情報処理装置。
  19. 前記算出部は、前記血管血流量と前記心筋血流量との比率に基づく前記指標を算出する、請求項１～１８のいずれか１つに記載の医用情報処理装置。
  20. 最大充血状態又は安静状態での冠動脈における血管血流量を流体解析によって取得し、前記血管血流量の取得時と同一の状態での前記冠動脈によって血液が供給される心筋領域における心筋血流量をパーフュージョン解析によって取得し、
    前記血管血流量と前記心筋血流量との比較に基づいて、前記心筋領域に対して前記血液を供給する毛細血管における毛細血管抵抗量を示す指標を算出する、
    各処理をコンピュータに実行させる、医用情報処理プログラム。
  21. 最大充血状態又は安静状態での冠動脈における血管血流量を流体解析によって取得し、前記血管血流量の取得時と同一の状態での前記冠動脈によって血液が供給される心筋領域における心筋血流量をパーフュージョン解析によって取得する取得部と、
    前記血管血流量と前記心筋血流量との比較に基づいて、前記心筋領域に対して前記血液を供給する毛細血管における毛細血管抵抗量を示す指標を算出する算出部と、
    を備える、医用情報処理システム。
JP2018191548A 2018-10-10 2018-10-10 医用情報処理装置、医用情報処理プログラム、医用情報処理システム Active JP7467026B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2018191548A JP7467026B2 (ja) 2018-10-10 2018-10-10 医用情報処理装置、医用情報処理プログラム、医用情報処理システム
US16/598,294 US20200113450A1 (en) 2018-10-10 2019-10-10 Medical information processing system and medical information processing method
JP2023060713A JP2023082164A (ja) 2018-10-10 2023-04-04 医用情報処理装置、医用情報処理プログラム、医用情報処理システム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018191548A JP7467026B2 (ja) 2018-10-10 2018-10-10 医用情報処理装置、医用情報処理プログラム、医用情報処理システム

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2023060713A Division JP2023082164A (ja) 2018-10-10 2023-04-04 医用情報処理装置、医用情報処理プログラム、医用情報処理システム

Publications (2)

Publication Number Publication Date
JP2020058591A JP2020058591A (ja) 2020-04-16
JP7467026B2 true JP7467026B2 (ja) 2024-04-15

Family

ID=70161863

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2018191548A Active JP7467026B2 (ja) 2018-10-10 2018-10-10 医用情報処理装置、医用情報処理プログラム、医用情報処理システム
JP2023060713A Pending JP2023082164A (ja) 2018-10-10 2023-04-04 医用情報処理装置、医用情報処理プログラム、医用情報処理システム

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2023060713A Pending JP2023082164A (ja) 2018-10-10 2023-04-04 医用情報処理装置、医用情報処理プログラム、医用情報処理システム

Country Status (2)

Country Link
US (1) US20200113450A1 (ja)
JP (2) JP7467026B2 (ja)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019134757A1 (en) * 2018-01-08 2019-07-11 Brainlab Ag Optimizing an atlas
JP7467026B2 (ja) 2018-10-10 2024-04-15 キヤノンメディカルシステムズ株式会社 医用情報処理装置、医用情報処理プログラム、医用情報処理システム
EP3671649A1 (en) * 2018-12-19 2020-06-24 Siemens Healthcare GmbH Method and computer system for generating a combined tissue-vessel representation
CN113837985B (zh) * 2020-06-24 2023-11-07 上海博动医疗科技股份有限公司 用于血管造影图像处理的训练方法和装置、自动处理方法和装置
JP7548755B2 (ja) 2020-08-27 2024-09-10 キヤノンメディカルシステムズ株式会社 医用画像処理装置、システム及び方法
CN115869003A (zh) * 2022-12-30 2023-03-31 杭州脉流科技有限公司 基于ct图像的冠状动脉微循环阻力指数计算方法和装置
KR102611037B1 (ko) * 2023-07-04 2023-12-08 주식회사 휴런 뇌와 관련된 정보로부터의 혈관 함수를 추출하기 위한 장치 및 방법

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014100540A (ja) 2012-10-22 2014-06-05 Fujita Gakuen 医用画像診断装置、医用画像処理装置及び医用画像処理プログラム
JP2016511682A (ja) 2013-03-04 2016-04-21 ハートフロー， インコーポレイテッド 血流特性のモデリングにおける感度解析方法及びシステム
US20170095292A1 (en) 2015-10-02 2017-04-06 Heartflow, Inc. System and method for diagnosis and assessment of cardiovascular disease by comparing arterial supply capacity to end-organ demand
JP2018503419A (ja) 2014-12-08 2018-02-08 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． インタラクティブな心臓テストデータならびに関連するデバイス、システム、及び方法
US20200113450A1 (en) 2018-10-10 2020-04-16 Canon Medical Systems Corporation Medical information processing system and medical information processing method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014100540A (ja) 2012-10-22 2014-06-05 Fujita Gakuen 医用画像診断装置、医用画像処理装置及び医用画像処理プログラム
JP2016511682A (ja) 2013-03-04 2016-04-21 ハートフロー， インコーポレイテッド 血流特性のモデリングにおける感度解析方法及びシステム
JP2018503419A (ja) 2014-12-08 2018-02-08 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． インタラクティブな心臓テストデータならびに関連するデバイス、システム、及び方法
US20170095292A1 (en) 2015-10-02 2017-04-06 Heartflow, Inc. System and method for diagnosis and assessment of cardiovascular disease by comparing arterial supply capacity to end-organ demand
JP2018534970A (ja) 2015-10-02 2018-11-29 ハートフロー， インコーポレイテッド 動脈の供給能力と終末器官の必要量とを比較することにより心血管疾患の診断及び評価をするためのシステム及び方法
US20200113450A1 (en) 2018-10-10 2020-04-16 Canon Medical Systems Corporation Medical information processing system and medical information processing method

Also Published As

Publication number Publication date
US20200113450A1 (en) 2020-04-16
JP2023082164A (ja) 2023-06-13
JP2020058591A (ja) 2020-04-16

Similar Documents

Publication Publication Date Title
JP7467026B2 (ja) 医用情報処理装置、医用情報処理プログラム、医用情報処理システム
US11896416B2 (en) Method for calculating coronary artery fractional flow reserve on basis of myocardial blood flow and CT images
US20200375565A1 (en) Medical information processing apparatus, x-ray ct apparatus, and medical information processing method
CN105792738B (zh) 用于改善功能狭窄分析的局部ffr估计和可视化
US11636635B2 (en) Medical information processing apparatus and medical information processing method
KR101939778B1 (ko) 필요 혈류량 결정 방법 및 장치, 혈류 영상 생성 방법 및 장치, 심근 관류 영상 처리 방법 및 장치
EP2932469B1 (en) Method of determining the blood flow through coronary arteries
US10052032B2 (en) Stenosis therapy planning
JP6929925B2 (ja) 動脈網における流量および圧力勾配を患者特定コンピュータ断層撮影アルゴリズムに基づくコントラスト分布から判断するための方法
US10278662B2 (en) Image processing apparatus and medical image diagnostic apparatus
US20200214659A1 (en) Medical-information processing apparatus and x-ray ct apparatus
JP6981807B2 (ja) 医用情報処理装置、ｘ線ｃｔ装置、医用情報処理プログラム、医用情報処理方法及び医用情報処理システム
JP7032111B2 (ja) 医用画像処理装置、ｘ線ｃｔ装置及び医用画像処理プログラム
US10157490B2 (en) Angulation planning for a three dimensional angiography
CN107773243B (zh) 利用不同记录模态的组合来确定临床特征参量
JP2019000628A (ja) 画像処理装置、医用画像診断装置及び画像処理プログラム
US9905001B2 (en) Image processing apparatus and image processing method
JP2017140365A (ja) 画像処理装置及び医用画像診断装置
JP2014046109A (ja) 医用画像処理装置
KR20170098643A (ko) 관상동맥 분지관의 체적을 이용한 혈류 역학 모델링 방법
JP7548755B2 (ja) 医用画像処理装置、システム及び方法
US20240215937A1 (en) Determination of hemodynamic indices
JP2022037581A (ja) 医用画像処理装置、システム及び方法
KR20180023359A (ko) 심근 체적의 분할을 이용한 혈류 역학 시뮬레이션 방법

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210830

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220825

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220906

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20230104

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230404

C60 Trial request (containing other claim documents, opposition documents)

Free format text: JAPANESE INTERMEDIATE CODE: C60

Effective date: 20230404

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20230412

C21 Notice of transfer of a case for reconsideration by examiners before appeal proceedings

Free format text: JAPANESE INTERMEDIATE CODE: C21

Effective date: 20230418

A912 Re-examination (zenchi) completed and case transferred to appeal board

Free format text: JAPANESE INTERMEDIATE CODE: A912

Effective date: 20230502

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240111

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240403

R150 Certificate of patent or registration of utility model

Ref document number: 7467026

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150