[go: up one dir, main page]

JP5934071B2 - 管状構造の最短経路探索装置および方法並びにプログラム - Google Patents

管状構造の最短経路探索装置および方法並びにプログラム Download PDF

Info

Publication number
JP5934071B2
JP5934071B2 JP2012213408A JP2012213408A JP5934071B2 JP 5934071 B2 JP5934071 B2 JP 5934071B2 JP 2012213408 A JP2012213408 A JP 2012213408A JP 2012213408 A JP2012213408 A JP 2012213408A JP 5934071 B2 JP5934071 B2 JP 5934071B2
Authority
JP
Japan
Prior art keywords
path
route
divided
search
tubular structure
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
JP2012213408A
Other languages
English (en)
Other versions
JP2014064824A (ja
JP2014064824A5 (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
Priority to JP2012213408A priority Critical patent/JP5934071B2/ja
Priority to US14/033,200 priority patent/US9198603B2/en
Publication of JP2014064824A publication Critical patent/JP2014064824A/ja
Publication of JP2014064824A5 publication Critical patent/JP2014064824A5/ja
Application granted granted Critical
Publication of JP5934071B2 publication Critical patent/JP5934071B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • A61B5/1075Measuring physical dimensions, e.g. size of the entire body or parts thereof for measuring dimensions by non-invasive methods, e.g. for determining thickness of tissue layer
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0073Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by tomography, i.e. reconstruction of 3D images from 2D projections
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/02007Evaluating blood vessel condition, e.g. elasticity, compliance
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
    • A61B5/055Detecting, 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
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30048Heart; Cardiac
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00Indexing scheme for image generation or computer graphics
    • G06T2210/41Medical
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03Recognition of patterns in medical or anatomical images
    • G06V2201/031Recognition of patterns in medical or anatomical images of internal organs

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Vascular Medicine (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Pulmonology (AREA)
  • Optics & Photonics (AREA)

Description

本発明は、被検体の画像から血管等の管状構造を抽出し、抽出した管状構造の最短経路を探索する装置および方法並びにプログラムに関するものである。
近年、医療機器(例えば多検出器型CT等)の進歩により質の高い3次元画像が画像診断に用いられるようになってきている。しかしながら、3次元画像は多数の2次元画像から構成され情報量が多いため、医師が所望の観察部位を見つけ診断することに時間を要する場合がある。そこで、注目する臓器を抽出しMIP、VR、CPR等の表示を行うことにより、臓器全体や病変の視認性を高め診断の効率化を図ることが行われている。
一方、心臓CT画像に対する解析、特に冠動脈解析あるいは脳血管の解析を行う際に、血管経路に沿って展開したCPR画像、狭窄部および狭窄率といった診断に有効な情報を得るために、管状構造を含む3次元画像から、冠動脈等の血管の中心経路をグラフ構造として抽出して表示することが行われている。例えば、特許文献1には、血管等の管状構造の始点から目標とする部位との間において、目標とする部位に通じる複数の管状構造の中心を通る中心線を抽出し、複数の中心線から最短となる最短経路を探索し、最短経路を対応する管状構造の表示画面に重ねて表示する手法が提案されている。
しかしながら、特許文献1に記載された手法においては、プラーク等により血管が細くなっている部分、あるいは脂肪等の管状構造とは異なる組織が付与してしまった部分等において、中心線が分断されてしまう場合がある。このように中心線が分断されてしまうと、最短経路を探索できなくなる。このため、中心線が分断されている場合に、ユーザによる途切れた部分の指定を受け付けることにより、途切れた部分を接続して、最短経路を探索できるようにする手法が提案されている(特許文献2参照)。
特開平11−120327号公報 特開2004−283373号公報
しかしながら、特許文献2に記載された手法においては、分断されている箇所をユーザが指定する必要がある。このため、管状構造の最短経路の確認の作業が煩雑なものとなり、作業を行うユーザの負担が大きい。
本発明は、上記事情に鑑みなされたものであり、ユーザに負担をかけることなく、3次元画像から管状構造の最短経路を探索できるようにすることを目的とする。
本発明による管状構造の最短経路探索装置は、管状構造を含む3次元画像から管状構造の経路を抽出する経路抽出手段と、
経路において、分断されている分断経路を検出する分断経路検出手段と、
分断経路を接続する検索用経路を生成する検索用経路生成手段と、
経路および検索用経路に基づいて、経路の始点および終点の間の最短経路を探索する探索手段とを備えたことを特徴とするものである。
なお、本発明による最短経路探索装置においては、分断経路検出手段を、経路の端点と端点の周囲にある他の経路との距離が特定距離以内である場合に、端点と他の経路との間を分断経路として検出する手段としてもよい。
また、本発明による最短経路探索装置においては、分断経路検出手段を、経路の端点における管状構造の径と、分断された経路の周囲にある複数の他の経路における管状構造の径とにより定まる距離を算出し、距離が最も短い他の経路と端点との間を分断経路として検出する手段としてもよい。
また、本発明による最短経路探索装置においては、経路の始点および終点を設定する始点終点設定手段をさらに備えるものとしてもよい。
また、本発明による最短経路探索装置においては、経路抽出手段を、抽出された経路に対する結合および/または分断の指示に基づいて、経路を修正する手段としてもよい。
また、本発明による最短経路探索装置においては、分断経路検出手段を、分断の指示がなされて修正された経路においては、分断経路としての検出を禁止する手段としてもよい。
また、本発明による最短経路探索装置においては、探索手段を、検索用経路に対して優先度を低くして最短経路を探索する手段としてもよい。
また、本発明による最短経路探索装置においては、経路抽出手段を、分断経路において経路の再度の抽出を行う手段とし、
検索用経路生成手段を、経路の再度の抽出結果に基づいて、検索用経路を生成する手段としてもよい。
本発明による最短経路探索方法は、管状構造を含む3次元画像から管状構造の経路を抽出し、
経路において、分断されている分断経路を検出し、
分断経路を接続する検索用経路を生成し、
経路および検索用経路に基づいて、経路の始点および終点の間の最短経路を探索することを特徴とするものである。
なお、本発明による管状構造の最短経路探索方法をコンピュータに実行させるためのプログラムとして提供してもよい。
本発明によれば、管状構造を含む3次元画像から管状構造の経路が抽出され、経路において、分断されている分断経路が検出され、分断経路を接続する検索用経路が生成され、経路および検索用経路に基づいて、経路の始点および終点の間の最短経路が探索される。このため、ユーザは何ら手を加えなくても、分断経路が接続されて、最短経路が探索されることとなる。したがって、3次元画像から管状構造の最短経路を探索する際のユーザの負担を軽減できる。
また、経路の端点と端点の周囲にある他の経路との距離が特定距離以内である場合に、端点と他の経路との間を分断経路として検出することにより、本来連続した経路でない端点が、経路として探索されてしまうことを防止できるため、最短経路探索の精度を向上できる。
また、経路の端点における管状構造の径と、分断された経路の周囲にある複数の他の経路における管状構造の径とにより定まる距離を算出し、この距離が最も短い他の経路と端点との間を分断経路として検出することによっても、本来連続した経路でない端点が、経路として探索されてしまうことを防止できるため、最短経路探索の精度を向上できる。
また、経路の始点および終点を検出することにより、ユーザは経路の始点および終点を指定する必要が無くなるため、ユーザの負担をさらに軽減することができる。
また、抽出された経路に対する結合および/または分断の指示を受け付け、この指示に基づいて、経路を修正することにより、ユーザの意図を反映させて最短経路を探索することができる。
この場合、分断の指示がなされて修正された経路においては、分断経路としての検出を禁止することにより、ユーザが意図的に分断した経路が検索用経路により接続されてしまうことを防止できる。
また、検索用経路に対して優先度を低くして最短経路を探索することにより、検索用経路が優先的に最短経路として探索されることが無くなるため、誤った経路が探索されてしまう可能性を低減できる。
また、分断経路の端点間において経路の再度の抽出を行うことにより、本来存在すべき経路が探索される場合がある。このため、経路の再度の抽出結果に基づいて、検索用経路を生成することにより、本来存在すべき経路の可能性が高い検索用経路を用いて最短経路を探索できるため、最短経路探索の精度を向上できる。
本発明の実施形態による管状構造の最短経路探索装置を適用した3次元画像表示装置の構成を示す概略ブロック図 抽出された冠動脈の経路を示す図 分断経路の検出を説明するための図(その1) 分断経路の検出を説明するための図(その2) 検索用経路の生成を説明するための図 始点および終点の設定を説明するための図 最短経路がマッピングされたVR画像を示す図 本実施形態において行われる処理を示すフローチャート
以下、図面を参照して本発明の実施形態について説明する。図1は本発明の実施形態による管状構造の最短経路探索装置を適用した3次元画像表示装置の構成を示す概略ブロック図である。なお、図1に示す3次元画像表示装置1の構成は、補助記憶装置に読み込まれた3次元画像表示プログラムをコンピュータ上で実行することにより実現される。このプログラムは、CD−ROM等の記録媒体に記録され、もしくはインターネット等のネットワークを介して配布され、コンピュータにインストールされる。なお、以降の説明においては、3次元画像を人体の胸部を撮影することにより取得されたものとし、管状構造を冠動脈とする。
本実施形態による3次元画像表示装置1は、画像取得部10、記憶部12、構造物抽出部14、経路抽出部16、分断経路検出部18、検索用経路生成部20、始点終点設定部22、経路探索部24、表示制御部26、入力部28および表示部30を備える。
画像取得部10は、マルチスライスCT装置またはMRI装置等のモダリティ2において、被写体の胸部を撮影して得られた3次元画像V0を取得する、通信インターフェースの機能を有する。なお、本実施形態においては、モダリティ2はマルチスライスCT装置であるものとする。3次元画像0はLAN経由でモダリティ2から送信される。
ここで、3次元画像V0は、診断対象となる胸部を断層面に垂直な方向に沿って順に得られる2次元の断層画像を積層することによって取得されるものであり、本実施形態においては、モダリティ2において撮影された複数の断層画像を重ね合わせることにより生成される。なお、CT装置を用いて取得した3次元画像は、3次元空間上での格子点を構成するボクセル(すなわち画素位置)毎にX線の吸収量を蓄えたデータとなり、各画素位置に対して1つの信号値(CT装置で撮影した場合は、X線の吸収量を示す値)が与えられたデータとなる。
なお、3次元画像V0には、DICOM(Digital Imaging and Communications in Medicine)規格で規定された付帯情報が付加される。付帯情報は、例えば、3次元画像を識別するための画像ID、被写体を識別するための患者ID、検査を識別するための検査ID、画像情報毎に割り振られるユニークなID(UID)、その画像情報が生成された検査日、検査時刻、その画像情報を取得するための検査で使用されたモダリティの種類、患者氏名、年齢、性別等の患者情報、検査部位(撮影部位、本実施形態においては胸部)、撮影条件(造影剤の使用有無や、放射線量等)、1回の検査で複数の画像を取得したときのシリーズ番号あるいは採取番号等の情報が含まれうる。
記憶部12は、ハードディスク等の大容量の記憶装置であり、3次元画像V0が記憶される。
構造物抽出部14は、3次元画像V0から心臓を抽出する。心臓を抽出する手法としては、例えば特開平2008−259682号公報に記載されたように、部位認識を行い、その結果に基づいて心臓を抽出する手法を用いることができる。特開平2008−259682号公報に記載された手法は、入力された3次元画像V0を構成する断層画像を正規化し、正規化された断層画像から多数の特徴量を算出し、正規化した断層画像毎に算出された特徴量を、AdaBoost手法によって得られた判別器に入力して、部位らしさを表す部位毎のスコアを算出し、算出された部位スコアを入力として、動的計画法を用いて、人体の体部の並び順が保たれるように断層画像に表された部位(すなわち心臓)を決定する手法である。また、テンプレートマッチングによる方法(例えば特開2002−253539号公報参照)および、各部位(すなわち心臓)の固有画像を用いた方法(例えば特開2003−10166号公報参照)等を用いることもできる。
経路抽出部16は、心臓における冠動脈の経路を抽出する。具体的には、構造物抽出部14が抽出した心臓内の局所領域ごとに、3×3のヘシアン(Hessian)行列の固有値を算出することにより線状構造の探索を行う。線状構造が含まれる領域では、ヘシアン行列の3つの固有値のうち1つは0に近い値となり、他の2つは相対的に大きな値となる。また、値が0に近い固有値に対応する固有ベクトルは、線状構造の主軸方向を示すものとなる。経路抽出部16は、この関係を利用して、局所領域ごとに、ヘシアン行列の固有値に基づいて線状構造らしさを判定し、線状構造が識別された局所領域については、その中心点を候補点として検出する。
そして、探索により検出された候補点を、所定のアルゴリズムに基づいて連結する。これにより、候補点および候補点同士を連結する血管枝(エッジ)からなる木構造が構築され、この木構造を冠動脈の経路として抽出する。なお、経路上における複数の候補点の座標情報や、血管枝の方向を示すベクトル情報は、候補点や血管枝の識別子とともに記憶部12に記憶される。続いて、検出された候補点ごとに、周辺のボクセルの値(CT値)に基づき、冠動脈の経路に垂直な断面において、血管の輪郭(血管の外壁)を識別する。形状の識別は、Graph-Cutsに代表される公知のセグメンテーション手法を用いて行う。以上の処理により、図2に示すような管状構造である冠動脈の経路K0および冠動脈が抽出され、その抽出された経路K0および冠動脈の特定に必要な情報が生成され、記憶部12に記憶される。
経路抽出部16は、上述したように冠動脈の中心線を経路として抽出するが、プラーク等により血管が細くなっている部分、あるいは脂肪等の血管とは異なる組織が付与してしまった部分等において、本来あるべき経路が分断されてしまう場合がある。分断経路検出部18は、このような本来あるべき経路が分断されてしまっている分断経路を検出する。このため、分断経路検出部18は、まず経路抽出部16が抽出した経路の端点を検出する。この端点の検出は、経路K0が途切れている点を検出することにより行う。
そして、検出した端点を中心とした所定距離Th0以内に、抽出された経路K0が存在するか否かを判定し、存在する場合にその端点と抽出された経路との間は分断されているものとして、そこを分断経路として検出する。図3は分断経路の検出を説明するための図である。分断経路検出部18は、図3(a)に示すように端点Ptを中心とした所定距離Th0の円形領域A0を設定し、円形領域A0内に他の経路が存在する場合には、端点Ptとその経路との間を分断経路として検出する。一方、図3(b)に示すように、円形領域A0内に他の経路が存在しない場合、端点Ptとその経路との間は分断経路としては検出しない。また、図3(c)に示すように、円形領域A0内に他の複数の経路K1,K2が存在する場合には、端点Ptと端点Ptに最も近い方の経路K1との間を分断経路として検出する。ここで、所定距離Th0はユーザが任意の値に設定することが可能である。
なお、このような分断経路の検出に代えて、経路の端点とその周囲にある他の経路との距離を算出し、この距離が最も短くなる他の経路と端点との間を分断経路として検出してもよい。この際、端点が位置する経路を抽出した冠動脈の径、および距離を算出する対象となる他の経路を抽出した冠動脈の径を用いて、端点とその周囲の他の経路との距離を判断してもよい。すなわち、図4に示すように、経路K3の端点Ptと経路K4との距離を算出するに際し、経路K3を抽出した血管の半径R3および他の経路K4を抽出した血管の半径R4を、端点Ptと経路K4との距離L0から減算した距離L1を算出し、その距離L1が最小となる他の経路と端点Ptとの間を分断経路として検出してもよい。なお、上記距離L0は、端点Ptと他の経路K4との最短距離である。
検索用経路生成部20は、分断経路検出部18が検出した分断経路を接続する検索用経路を生成する。すなわち、図5に示すように、経路K5の端点Ptと端点Ptに近い最も近い経路K6との間に、破線で示す検索用経路Kk0を生成する。この際、検索用経路Kk0は、端点Ptから経路K6に向かう経路K5の延長線を通るように生成すればよい。また、経路抽出部16において、分断経路すなわち端点Ptと経路K6との間の領域に対して、再度の経路抽出を行うようにしてもよい。このように再度の抽出処理を行う際には、端点Ptから経路K6に向かう直線の進め方を小さくする、さらには、端点Ptから経路K6の方向にのみ探索を行う等して、抽出の精度を向上させることが好ましい。このように抽出精度を高めることにより、経路が抽出される可能性があり、このように経路が抽出された場合には、その経路を検索用経路として用いればよい。
始点終点設定部22は、経路抽出部16が抽出した冠動脈の経路K0に始点および終点を設定する。具体的には、図6に示すように、表示制御部26により、心臓に経路K0をマッピングしたボリュームレンダリング画像(VR画像)を表示部30に表示し、ユーザによる入力部28を操作しての始点Psおよび終点Peの指定を受け付けることにより、始点Psおよび終点Peを設定する。なお、VR画像の表示については後述する。ここで、本実施形態においては、始点Psは冠動脈と大動脈との接続部分であり、終点Peは冠動脈の末梢端とする。なお、始点Psおよび終点Peを自動で検出するようにしてもよい。これにより、ユーザは経路の始点Psおよび終点Peを指定する必要が無くなるため、ユーザの負担を軽減することができる。
経路探索部24は、経路抽出部16が抽出した経路および検索用経路に基づいて、始点Psおよび終点Peの間の最短経路を探索する。すなわち、経路探索部24は、経路抽出部16が抽出した経路および検索用経路に沿って、始点Psから終点Peまでの長さ(距離)を検出するトラッキング処理を行う。そして、このトラッキング処理により始点Psから終点Peまでの最短経路を検出する。具体的には、始点Psから終点Peまでの長さが検出された経路のうち、最短の経路(中心線の長さが一番小さくなる経路の探索ベクトルデータ)を検出し、この経路を始点Psから終点Peまでの最短経路として検出する。
なお、このように最短経路を探索するに際し、検索用経路についてはトラッキングの際に、優先度を低くすることが好ましい。具体的には、抽出した経路よりも長さを長く見せるように重みづけをすることが好ましい。
表示制御部26は、構造物抽出部14が3次元画像V0から抽出した心臓および経路抽出部16が抽出した冠動脈を、ボリュームレンダリング法を用いてボリュームレンダリング(VR)表示する。すなわち、抽出した心臓および冠動脈に対して投影面から仮想的な光線を照射し、3次元画像V0の各信号値に対応した色(R,G,B)および不透明度(オパシティ)に基づいて、物体内部からの仮想的な反射光による3次元画像を作成し、この3次元画像から投影面に物体内部の3次元構造を透視する投影画像を生成して、これをVR画像として表示する。さらに、表示制御部26は、経路探索部24が探索した冠動脈上に最短経路をマッピング表示する。図7は最短経路がマッピングされたVR画像を示す図である。図7に示すように、冠動脈上に最短経路Ksがマッピングされている。
入力部28は、キーボードおよびマウス等の公知の入力装置からなる。
表示部30は、液晶、CRT等の公知の表示装置からなる。
次いで、本実施形態において行われる処理について説明する。図8は本実施形態において行われる処理を示すフローチャートである。なお、3次元画像V0は画像取得部10により取得されて記憶部12に記憶され、さらに構造物抽出部14により、3次元画像V0から心臓が抽出されているものとする。まず、経路抽出部16が、冠動脈の経路および冠動脈を抽出し(ステップST1)、さらに分断経路検出部18が、本来あるべき経路が分断されてしまっている分断経路を検出する(ステップST2)。
そして、検索用経路生成部20が、分断経路を接続する検索用経路を生成し(ステップST3)、始点終点設定部22が、経路抽出部16が抽出した冠動脈に始点Psおよび終点Peを設定する(ステップST4)。そして、経路探索部24が、経路抽出部16が抽出した経路および検索用経路に基づいて、始点Psおよび終点Peの間の最短経路を探索する(ステップST5)。さらに、表示制御部26が、心臓および冠動脈のVR画像に最短経路をマッピングし、最短経路がマッピングされたVR画像を表示部30に表示し(ステップST6)、処理を終了する。
このように、本実施形態においては、3次元画像V0から冠動脈の経路を抽出し、経路において分断経路を検出し、分断経路に検索用経路を生成し、抽出した経路および検索用経路に基づいて、始点Psおよび終点Peの間の最短経路を探索するようにしたものである。このため、ユーザは何ら手を加えなくても、分断経路が接続されて、最短経路が探索されることとなる。したがって、3次元画像V0から冠動脈等の管状構造の最短経路を探索する際のユーザの負担を軽減できる。
また、経路の端点と端点の周囲にある他の経路との距離が所定距離Th0以内である場合に、端点と他の経路との間を分断経路として検出しているため、本来連続した経路でない端点が経路として探索されてしまうことを防止でき、これにより、最短経路探索の精度を向上できる。
また、検索用経路に対して優先度を低くして最短経路を探索しているため、検索用経路が優先的に最短経路として探索されることが無くなり、これにより、誤った経路が探索されてしまう可能性を低減できる。
また、分断経路において経路の再度の抽出を行うことにより、本来存在すべき経路が探索される場合がある。このため、経路の再度の抽出結果に基づいて、検索用経路を生成することにより、本来存在すべき経路の可能性が高い検索用経路を用いて最短経路を探索できるため、最短経路探索の精度を向上できる。
なお、上記実施形態においては、最短経路を探索した後に、VR画像とともに最短経路を表示しているが、経路抽出部16が抽出した経路を、VR画像にマッピングして表示部30に表示し、ユーザによる入力部28を用いての経路の結合および/または分断の指示を受け付け、経路抽出部16がこの指示に応じて抽出した経路を修正するようにしてもよい。これにより、ユーザの意図を反映させて最短経路を探索することができる。また、この場合において、ユーザは経路を分断する指示を行う場合があるが、ユーザによる分断の指示がなされた経路においては、分断経路としての検出を禁止することが好ましい。これにより、ユーザが意図的に分断した経路が検索用経路により接続されてしまうことを防止できる。
また、上記実施形態においては、管状構造として、心臓の冠動脈を用いているが、これに限定されるものではなく、脳血管、肺の気管支、肝臓の血管(動脈、静脈、門脈)等、人体を構成する管状構造を有する構造物の最短経路の探索に、本発明を適用することができる。
1 3次元画像表示装置
10 画像取得部
12 記憶部
14 構造物抽出部
16 経路抽出部
18 分断経路検出部
20 検索用経路生成部
22 始点終点設定部
24 経路探索部
26 表示制御部
28 入力部
30 表示部

Claims (15)

  1. 管状構造を含む3次元画像から該管状構造の経路を抽出する経路抽出手段と、
    前記経路において、分断されている分断経路を検出する分断経路検出手段と、
    該分断経路を接続する検索用経路を生成する検索用経路生成手段と、
    前記経路および前記検索用経路に基づいて、前記経路の始点および終点の間の最短経路を探索する探索手段とを備え
    前記探索手段は、前記検索用経路に対して優先度を低くして前記最短経路を探索する手段であることを特徴とする管状構造の最短経路探索装置。
  2. 前記経路抽出手段は、前記抽出された経路に対する結合および/または分断の指示に基づいて、前記経路を修正する手段であることを特徴とする請求項記載の最短経路探索装置。
  3. 前記分断経路検出手段は、分断の指示がなされて修正された経路においては、前記分断経路としての検出を禁止する手段であることを特徴とする請求項記載の最短経路探索装置。
  4. 管状構造を含む3次元画像から該管状構造の経路を抽出する経路抽出手段と、
    前記経路において、分断されている分断経路を検出する分断経路検出手段と、
    該分断経路を接続する検索用経路を生成する検索用経路生成手段と、
    前記経路および前記検索用経路に基づいて、前記経路の始点および終点の間の最短経路を探索する探索手段とを備え、
    前記経路抽出手段は、前記抽出された経路に対する結合および/または分断の指示に基づいて、前記経路を修正する手段であり、
    前記分断経路検出手段は、分断の指示がなされて修正された経路においては、前記分断経路としての検出を禁止する手段であることを特徴とする管状構造の最短経路探索装置。
  5. 前記分断経路検出手段は、前記経路の端点と該端点の周囲にある他の経路との距離が特定距離以内である場合に、前記端点と前記他の経路との間を分断経路として検出する手段であることを特徴とする請求項1から4のいずれか1項記載の最短経路探索装置。
  6. 前記分断経路検出手段は、前記経路の端点における前記管状構造の径と、該端点の周囲にある複数の他の経路における前記管状構造の径とにより定まる距離を算出し、該距離が最も短い他の経路と前記端点との間を前記分断経路として検出する手段であることを特徴とする請求項1から4のいずれか1項記載の最短経路探索装置。
  7. 管状構造を含む3次元画像から該管状構造の経路を抽出する経路抽出手段と、
    前記経路において、分断されている分断経路を検出する分断経路検出手段と、
    該分断経路を接続する検索用経路を生成する検索用経路生成手段と、
    前記経路および前記検索用経路に基づいて、前記経路の始点および終点の間の最短経路を探索する探索手段とを備え、
    前記分断経路検出手段は、前記経路の端点における前記管状構造の径と、該端点の周囲にある複数の他の経路における前記管状構造の径とにより定まる距離を算出し、該距離が最も短い他の経路と前記端点との間を前記分断経路として検出する手段であることを特徴とする最短経路探索装置。
  8. 前記経路の始点および終点を設定する始点終点設定手段をさらに備えたことを特徴とする請求項1からのいずれか1項記載の最短経路探索装置。
  9. 前記経路抽出手段は、前記分断経路において前記経路の再度の抽出を行う手段であり、
    前記検索用経路生成手段は、前記経路の再度の抽出結果に基づいて、前記検索用経路を生成する手段であることを特徴とする請求項1からのいずれか1項記載の最短経路探索装置。
  10. 管状構造を含む3次元画像から該管状構造の経路を抽出し、
    前記経路において、分断されている分断経路を検出し、
    該分断経路を接続する検索用経路を生成し、
    前記経路および前記検索用経路に基づいて、前記経路の始点および終点の間の最短経路を探索するに際し、
    前記検索用経路に対して優先度を低くして前記最短経路を探索することを特徴とする管状構造の最短経路探索方法。
  11. 管状構造を含む3次元画像から該管状構造の経路を抽出し、
    前記経路において、分断されている分断経路を検出し、
    該分断経路を接続する検索用経路を生成し、
    前記経路および前記検索用経路に基づいて、前記経路の始点および終点の間の最短経路を探索するに際し、
    前記抽出された経路に対する結合および/または分断の指示に基づいて、前記経路を修正し、
    分断の指示がなされて修正された経路においては、前記分断経路としての検出を禁止することを特徴とする管状構造の最短経路探索方法。
  12. 管状構造を含む3次元画像から該管状構造の経路を抽出し、
    前記経路において、分断されている分断経路を検出し、
    該分断経路を接続する検索用経路を生成し、
    前記経路および前記検索用経路に基づいて、前記経路の始点および終点の間の最短経路を探索するに際し、
    前記経路の端点における前記管状構造の径と、該端点の周囲にある複数の他の経路における前記管状構造の径とにより定まる距離を算出し、該距離が最も短い他の経路と前記端点との間を前記分断経路として検出することを特徴とする管状構造の最短経路探索方法。
  13. 管状構造を含む3次元画像から該管状構造の経路を抽出する手順と、
    前記経路において、分断されている分断経路を検出する手順と、
    該分断経路を接続する検索用経路を生成する手順と、
    前記経路および前記検索用経路に基づいて、前記経路の始点および終点の間の最短経路を探索する手順とを有し、
    前記探索する手順は、前記検索用経路に対して優先度を低くして前記最短経路を探索する手順である、管状構造の最短経路探索方法をコンピュータに実行させることを特徴とするプログラム。
  14. 管状構造を含む3次元画像から該管状構造の経路を抽出する手順と、
    前記経路において、分断されている分断経路を検出する手順と、
    該分断経路を接続する検索用経路を生成する手順と、
    前記経路および前記検索用経路に基づいて、前記経路の始点および終点の間の最短経路を探索する手順とを有し、
    前記経路を抽出する手順は、前記抽出された経路に対する結合および/または分断の指示に基づいて、前記経路を修正する手順であり、
    前記分断経路を検出する手順は、分断の指示がなされて修正された経路においては、前記分断経路としての検出を禁止する手順である、管状構造の最短経路探索方法をコンピュータに実行させることを特徴とするプログラム。
  15. 管状構造を含む3次元画像から該管状構造の経路を抽出する手順と、
    前記経路において、分断されている分断経路を検出する手順と、
    該分断経路を接続する検索用経路を生成する手順と、
    前記経路および前記検索用経路に基づいて、前記経路の始点および終点の間の最短経路を探索する手順とを有し、
    前記分断経路を検出する手順は、前記経路の端点における前記管状構造の径と、該端点の周囲にある複数の他の経路における前記管状構造の径とにより定まる距離を算出し、該距離が最も短い他の経路と前記端点との間を前記分断経路として検出する手順である、管状構造の最短経路探索方法をコンピュータに実行させることを特徴とするプログラム。
JP2012213408A 2012-09-27 2012-09-27 管状構造の最短経路探索装置および方法並びにプログラム Active JP5934071B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2012213408A JP5934071B2 (ja) 2012-09-27 2012-09-27 管状構造の最短経路探索装置および方法並びにプログラム
US14/033,200 US9198603B2 (en) 2012-09-27 2013-09-20 Device, method and program for searching for the shortest path in a tubular structure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2012213408A JP5934071B2 (ja) 2012-09-27 2012-09-27 管状構造の最短経路探索装置および方法並びにプログラム

Publications (3)

Publication Number Publication Date
JP2014064824A JP2014064824A (ja) 2014-04-17
JP2014064824A5 JP2014064824A5 (ja) 2015-01-29
JP5934071B2 true JP5934071B2 (ja) 2016-06-15

Family

ID=50339534

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012213408A Active JP5934071B2 (ja) 2012-09-27 2012-09-27 管状構造の最短経路探索装置および方法並びにプログラム

Country Status (2)

Country Link
US (1) US9198603B2 (ja)
JP (1) JP5934071B2 (ja)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140253544A1 (en) * 2012-01-27 2014-09-11 Kabushiki Kaisha Toshiba Medical image processing apparatus
US10360672B2 (en) * 2013-03-15 2019-07-23 University Of Iowa Research Foundation Automated separation of binary overlapping trees
US9895803B1 (en) * 2015-06-19 2018-02-20 X Development Llc Calculating trajectory corridor for robot end effector
EP3565469B1 (en) * 2017-01-05 2022-08-24 Koninklijke Philips N.V. Visualizing a course of a vasculature structure with an occlusion
JP7144408B2 (ja) * 2017-06-02 2022-09-29 テルモ株式会社 経路選択支援システム、経路選択支援プログラムを記録した記録媒体、および経路選択支援方法
JP7127212B2 (ja) * 2019-04-25 2022-08-29 富士フイルム株式会社 画像の向き設定装置、方法およびプログラム
KR102356341B1 (ko) * 2020-02-10 2022-01-28 주식회사 메디픽셀 혈관 내 의료 도구 삽입 장치를 위한 가이드 데이터 제공 방법 및 장치
US20210401500A1 (en) * 2020-06-29 2021-12-30 Biosense Webster (Israel) Ltd. Efficient automatic finding of minimal ear-nose-throat (ent) path for probe
CN115953495B (zh) * 2023-03-14 2023-08-01 北京唯迈医疗设备有限公司 基于二维造影图像的智能路径规划装置、系统和存储介质

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0765154A (ja) * 1993-08-31 1995-03-10 Toshiba Corp 血管像の定量解析装置及びその定量解析方法
JP4053117B2 (ja) 1997-10-17 2008-02-27 東芝医用システムエンジニアリング株式会社 画像処理装置
JP4545971B2 (ja) 2001-03-05 2010-09-15 日本電信電話株式会社 医用画像識別システム,医用画像識別処理方法,医用画像識別用プログラムおよびその記録媒体
JP3766033B2 (ja) 2001-04-25 2006-04-12 富士写真フイルム株式会社 画像処理方法および装置並びにプログラム
DE60237719D1 (de) * 2001-10-16 2010-10-28 Koninkl Philips Electronics Nv Verfahren zur auswahl von ästen für das ausrichten von sonden
JP4421203B2 (ja) * 2003-03-20 2010-02-24 株式会社東芝 管腔状構造体の解析処理装置
WO2007065221A1 (en) * 2005-12-07 2007-06-14 Commonwealth Scientific And Industrial Research Organisation Linear feature detection method and apparatus
DE102006001884A1 (de) * 2006-01-13 2007-07-19 Siemens Ag Verfahren zur bildlichen Unterstützung einer Führung eines in einen Hohlraum eines Untersuchungsobjekts mindestens teilweise eingeführten medizinischen Instruments
WO2010058398A2 (en) * 2007-03-08 2010-05-27 Sync-Rx, Ltd. Image processing and tool actuation for medical procedures
JP2008259682A (ja) 2007-04-12 2008-10-30 Fujifilm Corp 部位認識結果修正装置、方法、およびプログラム
US8126232B2 (en) * 2008-03-05 2012-02-28 Siemens Aktiengesellschaft System and method for 3D vessel segmentation with minimal cuts
JP5242235B2 (ja) * 2008-05-08 2013-07-24 株式会社東芝 Mri画像診断装置及びmr画像形成方法
US9679389B2 (en) * 2009-05-19 2017-06-13 Algotec Systems Ltd. Method and system for blood vessel segmentation and classification
US8934686B2 (en) * 2009-11-26 2015-01-13 Algotec Systems Ltd. User interface for selecting paths in an image
JP5263995B2 (ja) * 2011-03-18 2013-08-14 富士フイルム株式会社 ネットワーク構築装置および方法ならびにプログラム
DE102011076188B4 (de) * 2011-05-20 2013-03-21 Siemens Aktiengesellschaft Vorrichtung und Verfahren zur Planung einer endovaskulären Prozedur mit einem medizinischen Instrument
CN103957832B (zh) * 2011-10-26 2016-09-28 皇家飞利浦有限公司 血管树图像的内窥镜配准

Also Published As

Publication number Publication date
US9198603B2 (en) 2015-12-01
JP2014064824A (ja) 2014-04-17
US20140088416A1 (en) 2014-03-27

Similar Documents

Publication Publication Date Title
JP5934071B2 (ja) 管状構造の最短経路探索装置および方法並びにプログラム
JP7309986B2 (ja) 医用画像処理方法、医用画像処理装置、医用画像処理システム及び医用画像処理プログラム
US8630467B2 (en) Diagnosis assisting system using three dimensional image data, computer readable recording medium having a related diagnosis assisting program recorded thereon, and related diagnosis assisting method
JP5584006B2 (ja) 投影画像生成装置、投影画像生成プログラムおよび投影画像生成方法
JP4937397B2 (ja) 医用画像診断支援装置および方法、並びにプログラム
JP6080249B2 (ja) 3次元画像表示装置および方法並びにプログラム
JP6080248B2 (ja) 3次元画像表示装置および方法並びにプログラム
JP5572470B2 (ja) 診断支援装置および方法並びにプログラム
US10052032B2 (en) Stenosis therapy planning
CN102007515B (zh) 用于对肺动脉进行分割的方法和系统
JP5812841B2 (ja) 画像診断支援装置、画像診断支援方法、画像診断支援プログラム
JP5700964B2 (ja) 医用画像処理装置、方法およびプログラム
JP2009028515A (ja) 医用画像処理装置及び医用画像診断装置
JP2010528750A (ja) 管状構造の検査
JP5415245B2 (ja) 医用画像表示装置および方法並びにプログラム
JP6301277B2 (ja) 診断補助画像生成装置および診断補助画像生成方法、並びに、診断補助画像生成プログラム
JP5380231B2 (ja) 医用画像表示装置および方法ならびにプログラム
JP2007135858A (ja) 画像処理装置
JP6734111B2 (ja) 所見情報作成装置及びシステム
JP6934948B2 (ja) 流体解析装置および流体解析装置の作動方法並びに流体解析プログラム
WO2014155917A1 (ja) 手術支援装置、方法およびプログラム
JP5159195B2 (ja) 医用画像処理装置
JP2022047359A (ja) 医用画像処理装置、医用画像処理システム、及び、医用画像処理方法

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20141205

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20141205

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150819

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150908

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20151106

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160506

R150 Certificate of patent or registration of utility model

Ref document number: 5934071

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

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