JP2023521773A - 血管映像を自動処理する方法及び装置 - Google Patents

血管映像を自動処理する方法及び装置 Download PDF

Info

Publication number
JP2023521773A
JP2023521773A JP2022561515A JP2022561515A JP2023521773A JP 2023521773 A JP2023521773 A JP 2023521773A JP 2022561515 A JP2022561515 A JP 2022561515A JP 2022561515 A JP2022561515 A JP 2022561515A JP 2023521773 A JP2023521773 A JP 2023521773A
Authority
JP
Japan
Prior art keywords
blood vessel
branch
error
image processing
target
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
JP2022561515A
Other languages
English (en)
Other versions
JP7427293B2 (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.)
Medipixel Inc
Original Assignee
Medipixel Inc
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 Medipixel Inc filed Critical Medipixel Inc
Publication of JP2023521773A publication Critical patent/JP2023521773A/ja
Application granted granted Critical
Publication of JP7427293B2 publication Critical patent/JP7427293B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices 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
    • 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/46Arrangements for interfacing with the operator or the patient
    • A61B6/467Arrangements for interfacing with the operator or the patient characterised by special input means
    • 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/46Arrangements for interfacing with the operator or the patient
    • A61B6/467Arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B6/469Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
    • 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/50Apparatus 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/504Apparatus 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
    • 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/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/77Retouching; Inpainting; Scratch removal
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/187Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/44Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
    • G06V10/457Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components by analysing connectivity, e.g. edge linking, connected component analysis or slices
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/98Detection or correction of errors, e.g. by rescanning the pattern or by human intervention; Evaluation of the quality of the acquired patterns
    • G06V10/993Evaluation of the quality of the acquired pattern
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/40ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • 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/48Diagnostic techniques
    • A61B6/481Diagnostic techniques involving the use of contrast agents
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00Indexing scheme for image data processing or generation, in general
    • G06T2200/24Indexing scheme for image data processing or generation, in general involving graphical user interfaces [GUIs]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20092Interactive image processing based on input by user
    • G06T2207/20096Interactive definition of curve of interest
    • 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
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30168Image quality inspection
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Public Health (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Quality & Reliability (AREA)
  • Human Computer Interaction (AREA)
  • Dentistry (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Physiology (AREA)
  • Multimedia (AREA)
  • Vascular Medicine (AREA)
  • Geometry (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Analysis (AREA)

Abstract

プロセッサによって実行される血管映像から血管映像を処理する方法であって、前記血管映像から対象血管を抽出するステップと、前記対象血管の抽出結果で前記対象血管に関する血管構造データ、前記対象血管の曲率情報、前記対象血管の直径情報、及び前記対象血管の明度情報のうち少なくとも1つに基づいてエラー部位を識別するステップと、前記対象血管でエラー部位が識別される場合に応答して、前記識別されたエラー部位を修正するステップとを含む。

Description

以下、血管映像を処理する方法に関する技術が提供される。
血管造影検査(angiography)映像は、主な血管を観察して血管内の問題部位を診断して必要な施術及び措置を行うために広範囲に活用されている。従来には、便宜性を高めて診断結果を定量化するために、血管映像から主な血管を自動識別する映像処理エンジンが利用されている。例えば、Piメディカル社のCaas QCAエンジンは、血管造影映像から主な血管を探して表示することができる。一方、エンジンから抽出された主な血管は、実際に抽出しなければならない主な血管部位ではない、他の血管部位(misidentification)に該当し、一部が切れている(disconnection)など、エラーを示す場合が多い。したがって、臨床では、エンジンが主な血管を1次抽出した識別結果をそのまま使用されることができず、人力が投入されて血管造影映像を確認し前記エラーを人力が直接修正して活用して格納する。
課題を解決しようとする手段
一実施形態によれば、プロセッサによって行われる血管映像処理方法は、血管映像から対象血管を抽出するステップと、前記対象血管の抽出結果で前記対象血管に関する血管構造データ、前記対象血管の曲率情報、前記対象血管の直径情報、及び前記対象血管の明度情報のうち少なくとも1つに基づいてエラー部位を識別するステップと、前記対象血管の抽出結果でエラー部位が識別される場合に応答して、前記識別されたエラー部位を修正するステップとを含む。
一実施形態に係る血管映像処理方法の前記エラー部位を識別するステップは、前記血管映像から血管領域及び非血管領域を分割するステップと、前記血管領域で分枝地点を基準にして区分される血管分枝のそれぞれに対応するノードに基づいて前記血管構造データを生成するステップと、前記対象血管のノード間の連結データを前記血管構造データに比較するステップと、前記血管構造データのトポロジーに基づいてノード間の連結エラーが検出された場合、前記連結エラーに対応するノードを前記エラー部位として判断するステップとを含むことができる。
一実施形態に係る血管映像処理方法の前記エラー部位を修正するステップは、
前記エラー部位が断絶部位として識別される場合に応答して、前記断絶部位に対応する領域を連結するステップと、前記エラー部位が誤認部位として識別される場合に応答して、前記誤認部位に対応するノードを血管構造データにマッチングするノードに交替するステップとを含むことができる。
一実施形態に係る血管映像処理方法の前記エラー部位を識別するステップは、前記対象血管の抽出結果で互いに隣接する上位分枝と下位分枝との間の曲率が閾値曲率を超過する場合に応答して、前記下位分枝を誤認部位として識別するステップを含むことができる。
一実施形態に係る血管映像処理方法の前記エラー部位を修正するステップは、前記識別された誤認部位に対して、前記下位分枝を前記上位分枝との曲率が臨界曲率以下である他の候補下位分枝に交替するステップを含むことができる。
一実施形態に係る血管映像処理方法の前記エラー部位を識別するステップは、前記対象血管の抽出結果で互いに隣接する上位分枝と下位分枝との間の直径の差が閾値以上である場合に応答して、前記下位分枝をエラー部位として識別するステップを含むことができる。
一実施形態に係る血管映像処理方法の前記エラー部位を修正するステップは、前記識別されたエラー部位に対して、前記下位分枝を前記上位分枝との直径の差に基づいた交替、及び前記対象血管に対する新しい直径情報を含む新しい抽出結果の生成のうち少なくとも1つを行うことによって前記識別されたエラー部位を修正するステップを含むことができる。
一実施形態に係る血管映像処理方法の前記エラー部位を識別するステップは、前記対象血管の抽出結果で互いに隣接する上位分枝と下位分枝との間の明度差が閾値の明度差以上である場合に応答して、前記下位分枝を誤認部位として識別するステップを含むことができる。
一実施形態に係る血管映像処理方法の前記エラー部位を修正するステップは、前記識別されたエラー部位に対して、前記下位分枝を前記上位分枝との明度差が閾値の明度差未満である他の候補下位分枝に交替するステップを含むことができる。
一実施形態に係る血管映像処理方法の前記エラー部位を修正するステップは、エラー部位が識別された場合に応答して、前記対象血管に対する新しい抽出結果を生成するステップを含むことができる。
一実施形態に係る血管映像処理方法の前記エラー部位を修正するステップは、
前記エラー部位が誤認部位である場合に応答して、1つ以上の候補分枝をユーザに提供するステップと、前記ユーザから前記1つ以上の候補分枝のいずれかの分枝に対する選択入力を受信する場合に応答して、前記エラー部位に対応する分枝を前記選択された分枝に交替するステップとを含むことができる。
一実施形態に係る血管映像処理方法の前記エラー部位を修正するステップは、複数の地点でユーザ入力が検出される場合に応答して、前記ユーザ入力の開始地点に対応する血管分枝に対応する領域及び前記ユーザ入力の終了地点に対応する血管分枝に対応する領域を連結するステップを含むことができる。
一実施形態に係る血管映像処理装置は、血管映像を受信する映像受信部と、前記血管映像から対象血管を抽出し、前記対象血管の抽出結果で前記対象血管に関する血管構造データ、前記対象血管の曲率情報、前記対象血管の直径情報、及び前記対象血管の明度情報のうち少なくとも1つに基づいてエラー部位を識別し、前記対象血管の抽出結果でエラー部位が識別される場合に応答して、前記識別されたエラー部位を修正するプロセッサを含む。
一実施形態に係る血管映像処理装置によって行われる血管映像を処理する方法は、造影剤が挿入された血管映像から主な血管に対応する対象血管を抽出することができ、対象血管におけるエラー部位を人力が投入されて対象血管を確認する必要なく、自動にエラー部位を識別することができる。また、血管映像処理装置は、対象血管で識別されたエラー部位を自動に修正することができ、ユーザから入力を受信してエラー部位をユーザの意図の通りに修正することができる。
一実施形態に係る血管映像を処理する方法を示すフローチャートである。 一実施形態に係る対象血管が表示された血管映像を示す図である。 一実施形態に係る対象血管の抽出結果でエラー部位が発生した例示を示す。 血管分枝に対応するノードに基づいた血管構造データのトポロジーを用いて対象血管のエラーを識別する方法を示す。 一実施形態に係るエラー部位を修正する方法を示す。 一実施形態に係るエラー部位を修正する方法を示す。 対象血管の曲率情報を用いて対象血管のエラーを識別して修正する方法を示す。 対象血管の直径情報を用いて対象血管のエラーを識別して修正する方法を示す。 対象血管の明度情報を用いて対象血管のエラーを識別して修正する方法を示す。 一実施形態に係るエラー部位を修正する方法を示す。 一実施形態に係るエラー部位を修正する方法を示す。 一実施形態に係る血管映像処理装置を概略的に示すブロック図である。
実施形態に対する特定な構造的又は機能的な説明は単なる例示のための目的として開示されたものであって、様々な形態に変更されることができる。したがって、実施形態は特定な開示形態に限定されるものではなく、本明細書の範囲は技術的な思想に含まれる変更、均等物ないし代替物を含む。
第1又は第2などの用語を様々な構成要素を説明するために用いることがあるが、このような用語は1つの構成要素を他の構成要素から区別する目的としてのみ解釈されなければならない。例えば、第1構成要素は第2構成要素と命名することができ、同様に、第2構成要素は第1構成要素にも命名することができる。
いずれかの構成要素が他の構成要素に「連結」されていると言及されたときには、その他の構成要素に直接的に連結されているか又は接続されているが、中間に他の構成要素が存在し得るものと理解されなければならない。
単数の表現は、文脈上、明白に異なる意味をもたない限り複数の表現を含む。本明細書において、「含む」又は「有する」等の用語は、明細書上に記載した特徴、数字、ステップ、動作、構成要素、部品又はこれを組み合わせたものが存在することを示すものであって、1つ又はそれ以上の他の特徴や数字、ステップ、動作、構成要素、部品、又はこれを組み合わせたものなどの存在又は付加の可能性を予め排除しないものとして理解しなければならない。
異なるように定義さがれない限り、技術的又は科学的な用語を含んで、ここで用いる全ての用語は、本実施形態が属する技術分野で通常の知識を有する者によって一般的に理解されるものと同じ意味を有する。一般的に用いられる予め定義された用語は、関連技術の文脈上で有する意味と一致する意味を有するものと解釈されなければならず、本明細書で明白に定義しない限り、理想的又は過度に形式的な意味として解釈されることはない。以下、添付する図面を参照しながら実施形態を詳細に説明する。各図面に提示された同一の参照符号は同一の部材を示す。
図1は、一実施形態に係る血管映像を処理する方法を示すフローチャートである。
まず、ステップ110において、血管映像処理装置は、血管映像(blood vessel image)から対象血管(target vessel)を抽出する。対象血管は、主血管のように示してもよい。一実施形態によれば、血管映像処理装置の映像受信部は、血管映像撮影装置によって撮影された血管映像を受信する。血管映像は、生体(living body)の血管を撮影した映像であって、血管造影検査(coronary angiography、以下CAG)映像及び/又は磁気共鳴映像(Magnetic Resonance Imaging、以下、MRI)を用いて生成されてもよい。例示的に、血管映像は、造影剤の注入された生体に対してX線撮影することによって取得された映像であり得る。
一実施形態によれば、血管映像処理装置は、機械学習モデル(machine learning model)に基づいて血管映像から対象血管を抽出することができる。機械学習モデルは、血管映像の入力に応答して、血管映像から対象血管を抽出するように設計された機械学習構造を有する1つ以上のモデルとして、例えば、ニューラルネットワークを含んでもよい。血管映像処理装置は、受信された血管映像に対して上述した機械学習モデルによる演算を行うことで、対象血管の抽出結果を算出することができる。例えば、機械学習モデルの出力データは、血管映像の複数のピクセルで各ピクセルが対象血管を指示する可能性(例えば、確率)に対応するスコアを含む。血管映像処理装置は、出力データで閾値以上のスコアを有するピクセルを対象血管として決定することで、対象血管の抽出結果を生成することができる。異なる例として、機械学習モデルの出力データは、血管映像から分割(segment)された対象血管領域として、血管映像の複数のピクセルのうち、対象血管として抽出されたピクセルを含むことができる。対象血管の抽出結果は、例えば、血管映像のピクセルのうち対象血管として抽出されたピクセルの集合及び/又は血管映像から分割された対象血管領域に対応する映像(例えば、対象血管映像)であり得る。
参考として、ニューラルネットワークは、ディープニューラルネットワーク(DNN、deep neural network)を含む。DNNは、完全接続ネットワーク(fully connected network)、ディープコンボリューショナルネットワーク(deep convolutional network)及びリカレントニューラルネットワーク(recurrent neural network)などを含む。ニューラルネットワークは、ディープラーニングに基づいて非線形関係にある入力データ及び出力データを互いにマッピングすることで、オブジェクト分類、オブジェクト認識、及びレーダーイメージ認識などを行うことができる。ディープラーニングは、ビッグデータセットからオブジェクト認識のような問題を解決するための機械学習方式であって、教師有り学習又は教師なし学習を介して入力データ及び出力データを互いにマッピングすることができる。教師有り学習学習の場合、上述した機械学習モデルは、トレーニング入力(例えば、トレーニングのための血管映像)及び該当トレーニング入力にマッピングされたトレーニング出力(例えば、トレーニングのための血管映像に対して専門家などによって対象血管に分割された真の値(ground truth)映像)の対を含むトレーニングデータに基づいてトレーニングされ得る。例えば、機械学習モデルは、トレーニング入力からトレーニング出力を出力するようにトレーニングされてもよい。トレーニングのうちの機械学習モデル(以下、「臨時モデル」)は、トレーニング入力に応答して臨時出力を生成することができ、臨時出力及びトレーニング出力(例えば、真の値)間の損失が最小化されるようトレーニングされる。トレーニング過程の間に機械学習モデルのパラメータ(例えば、ニューラルネットワークでノード/レイヤ間の連結加重値)が損失に応じてアップデートされることができる。
但し、機械学習モデルが血管映像から対象血管を直ちに抽出する例示について説明したが、これに限定されることはない。例えば、機械学習モデルは、全体血管抽出モデル及び対象血管抽出モデルを含む。全体血管抽出モデルは、血管映像から全体血管領域を抽出するように設計されたモデルであってもよい。対象血管抽出モデルは、全体血管領域を指示する映像(例えば、全体血管映像)から対象血管領域を抽出するよう設計されたモデルであってもよい。また、全体血管抽出モデルの代わりに、血管映像処理装置は、血管映像でのピクセルの周辺ピクセルとのグレースケールレベル差に基づいて境界を検出することで、全体血管領域を抽出することができる。例示的に、血管映像処理装置は、任意のピクセルと周辺ピクセルのグレースケールレベルの勾配(gradient)値が閾値の勾配値よりも大きい場合、該当ピクセルを境界として検出することができる。したがって、血管映像処理装置は、グレースケールレベルが急激に変わる領域を境界として検出することができる。血管映像処理装置は、グレースケールレベルの勾配値に基づいて抽出された全体血管映像から対象血管抽出モデルを用いて対象血管映像を抽出してもよい。
さらに、血管映像処理装置は、血管の形態、種類及び/又は血管領域に応じて複数の機械学習モデルのうち、対象血管抽出に使用する機械学習モデルを選択的に使用してもよい。一実施形態によれば、血管映像処理装置は、血管の種類別(例えば、左主冠状動脈(Left Main Coronary Artery、LM)、左前下行動脈(Left Anterior Descending Artery、LAD)、左回旋動脈(Left Circumflex Artery、LCX)、右冠状動脈(Right Coronary Artery、RCA))及び/又は血管領域別(例えば、近位領域(proximal region)、中間領域(mid region)、及び遠位領域(distal region))の複数の機械学習モデルを格納することができる。参考として、血管領域は、カテーテル(catheter)の挿入される血管地点から離れた距離に応じて、近位地点(proximal portion)、中間領域(middle portion)、及び遠位領域(distal portion)に分類されるが、これに限定されない。血管領域は、血管映像を取得するために造影剤が血管挿入部に注入される地点から離れた距離、及び造影剤が注入され得る血管末端から離れた距離の比率に応じて分類してもよい。例えば、血管映像処理装置は、抽出しようとする血管の種類を選択し、識別された血管の種類に対応する機械学習モデルをロードすることができる。血管映像処理装置は、ロードされた機械学習モデルを用いて、血管映像から選択された血管の種類に対応する対象血管の抽出結果を生成することができる。例示的に、血管映像処理装置は、複数の心血管の種類(例えば、1個の右冠状動脈、2個の左冠状動脈)に対応する機械学習モデルを格納することができる。複数の心血管の種類に対応する機械学習モデルのそれぞれは、該当心血管の種類に対応するトレーニングデータに基づいてトレーニングされてもよい。心血管の種類ごとに機械学習モデルのトレーニングされたパラメータはそれぞれ異なってもよく、さらに、互いに機械学習構造(例えば、畳み込みニューラルネットワーク、U-net構造など)が異なってもよい。説明の便宜のために、血管の種類ごとにパラメータ及び/又は機械学習構造が区別される機械学習モデルについて説明したが、これに限定されることはない。血管映像処理装置は、血管の形態、種類、及び/又は血管領域ごとに互いに区別される複数の機械学習モデルを格納し、必要な血管モデルを選択的にロードして対象血管抽出に使用してもよい。
そして、ステップ120において、血管映像処理装置は、対象血管からエラー部位(error portion)を識別(identify)する。エラー部位は、断絶部位(discontinuity portion)及び/又は誤認部位(misidentification portion)を含む。断絶部位は、対象血管として抽出された対象血管領域及び/又は対象血管を指示するピクセルのうち、少なくとも1つの領域及び/又は少なくとも1つのピクセルが他の領域及び/又は他のピクセルから分離されたり離隔された部位を示す。誤認部位は、血管映像で実際に抽出しなければならない対象血管ではない他の血管に対応する領域及び/又はピクセルが対象血管であると間違って抽出された部位を示す。例示的に、血管映像処理装置は、対象血管の抽出結果から、対象血管に関する血管構造データ、対象血管の曲率情報、対象血管の直径情報、及び対象血管の明度情報のうち少なくとも1つに基づいてエラー部位を自動識別することができる。エラー部位識別に対する説明は後述する。
次に、ステップ130において、血管映像処理装置は、対象血管でエラー部位が識別される場合に応答して、識別されたエラー部位を修正する。血管映像処理装置は、識別されたエラー部位を自動に修正できるが、これに限定されることはない。血管映像処理装置は、ユーザから受信された入力に基づいてエラー部位を修正してもよい。エラー部位修正については後述する。
図2は、一実施形態に係る対象血管が表示された血管映像を示す図である。
血管映像撮影装置は、血管映像200を撮影することができる。血管映像処理装置は、血管映像撮影装置から血管映像200を受信する。血管映像処理装置は、受信された血管映像200から対象血管220を抽出する。血管映像処理装置は、機械学習モデルに基づいて対象血管220を抽出し、対象血管220に対する抽出結果からエラー部位を識別して修正できる。血管映像処理装置は、分枝地点(branch point)を基準にして、エラー部位を識別して修正することができる。分枝地点は、2以上の血管分枝が合わさる地点及び/又は複数の血管分枝に分かれる地点を示す。血管映像200にキャプチャーされた血管領域は、複数の分枝地点を有してもよい。
一実施形態に係る血管映像処理装置は、最終的に図2に示したような対象血管220を決定するために、対象血管220の抽出、抽出結果からエラー部位の識別及び修正を行うことができる。以下の図3では、説明の便宜のために分枝地点を含む領域211を中心に血管映像処理方法によるエラー部位識別及び修正について説明する。但し、血管映像処理方法は、いずれかの分枝地点に対してのみ適用されるものではなく、血管映像から抽出された対象血管を指示する対象血管領域内に存在する1つ以上の分枝地点(例えば、全ての分枝地点のそれぞれ)に対して適用されてもよい。また、血管映像処理装置は、分枝地点を含む領域211だけでなく、分枝地点を含んでいない領域212に対してもエラー部位の識別及び修正が可能である。
図3は、一実施形態に係る対象血管の抽出結果でエラー部位が発生した例示を示す。
一実施形態に係る血管映像処理装置は、血管映像から分枝地点及び分枝を識別することができる。例えば、血管映像処理装置は、全体血管領域を抽出し、抽出された全体血管領域から分枝地点及び分枝を識別することができる。血管映像処理装置は、分枝地点(branch point)を基準にして血管分枝を識別することができる。血管映像から全体血管領域を抽出することができる。一実施形態に係る血管映像処理装置は、全体血管抽出モデルに基づいて血管映像から全体血管領域を抽出することができる。全体血管抽出モデルは、上述したように、血管映像から血管領域及び残りの非血管領域が区分された結果を指示する出力データを生成するようにトレーニングされたモデルであり得る。
図3は、抽出された全体血管領域のうち、図2に図示された血管映像200の一部領域210に対応する部分300を示す。血管映像処理装置は、一部の領域210内で分枝地点を基準に血管分枝を識別することができる。参考として、本明細書において、上位血管分枝は、血流の進行方向又は造影剤の進行方向によって分枝地点以前の分枝を示し、下位血管分枝は、分枝地点以後の分枝を示し、母血管は、血管領域内の最上位分枝を示す。例えば、血管映像処理装置によって識別された結果として、第1分枝ないし第7分枝301~307及び分枝地点311、312が図示されている。第1分枝地点311を基準にして、第1分枝300は上位血管分枝であり、第2分枝302及び第3分枝303は下位血管分枝である。また、上位血管分枝及び下位血管分枝の関係は相対的なものであって、分枝地点ごとに変わり得る。例えば、第2分枝302は、第4分枝304及び第5分枝305については上位血管分枝である。血管映像処理装置は、血管分枝それぞれに対応するノードを生成してインデックスすることで、血管構造データを生成することができる。血管構造データは、下記の図4を参照して後述する。。
図3に示された例示で、血管映像処理装置は対象血管に属する血管分枝として、第1分枝301、第3分枝303、及び第5分枝305を抽出することができる。但し、図2に示された対象血管220に対応する血管経路320は、第1分枝301、第2分枝302、及び第5分枝305を通過する。したがって、例示的な対象血管の抽出結果において、第1分枝301から第3分枝303に向かう第1分枝地点311はエラー部位である。以下では、エラー部位の検出及び修正をについて説明する。
図4は、血管分枝に対応するノードに基づいた血管構造データのトポロジーを用いて対象血管のエラーを識別する方法を示す。
一実施形態によれば、血管映像処理装置は、対象血管に関する血管構造データ400に基づいてエラー部位(error portion)を識別することができる。血管構造データ400は、血管の構造を指示するトポロジーデータとして、例えば、血管分枝に対応するノードがインデクシングされたトリ構造データであってもよい。例えば、血管映像処理装置は、全体血管領域に対して血管分枝地点を基準にして血管分枝を識別し、識別された血管分枝のそれぞれに対応するノードを生成し、生成されたノードのインデクシングに基づいて血管構造データ400を生成することができる。
血管映像処理装置は、上位血管分枝から分枝地点を基準にして複数の下位血管分枝に分岐される場合、上位血管分枝に対応するノードと下位血管分枝に対応するノードを連結することができる。血管映像処理装置は、血管映像から識別される分枝地点及び血管分枝に対して上述したノード連結を繰り返すことで、全体血管に対する血管構造データ400を生成することができる。例えば、図4に示された血管構造データ400において、第1ノード~第7ノードN1~N7は、順次にそれぞれ図3に示された第1分枝~第7分枝301~307に対応する。第1ノードN1は最上位ノードであり、第2ノードN2及び第3ノードN3は第1ノードN1の下位ノードである。残りのノードも類似の上下関係を有する。
一実施形態に係る血管映像処理装置は、血管映像から対象血管として抽出された血管分枝に対応するノード間のノード連結データ410を算出することができる。例えば、血管映像処理装置は、図3に示された対象血管の抽出結果から第1ノードN1、第3ノードN3、及び第5ノードN5が連結されたノード連結データ410を算出することができる。
血管映像処理装置は、ノード連結データ410を血管構造データ400と比較することで対象血管のエラー部位を識別することができる。血管映像処理装置は、生成された血管構造データ400に基づいて対象血管のノード間の連結エラーを検出する。血管映像処理装置は、連結エラーが検出されたノードに対応する血管分枝及び/又は分枝地点をエラー部位として判断する。血管映像処理装置は、ノード連結データ410で血管構造データ400にマッチングしないノードに基づいてエラー部位を判断することができる。例示的に、図4で抽出されたノード連結データ410は、第1ノードN1、第3ノードN3、及び第5ノードN5が順に連結されたものを指示する。血管構造データ400において、第3ノードN3は第5ノードN5と連結されない。言い換えれば、血管映像処理装置は、ノード連結データ410において、第3ノードN3と第5ノードN5との間の連結が血管構造データ400にマッチングしないものと決定することができる。第3ノードN3と第5ノードN5との間の連結は、連結エラーである。血管映像処理装置は、連結エラーに対応するノードのうち少なくともいずれかのノードに対応する血管分枝をエラー部位として判断することができる。例えば、図4における血管映像処理装置は、第3ノードN3及び第5ノードN5のうち少なくとも1つに対応する血管分枝をエラー部位として判断することができる。上述した例示は、第3ノードN3又は第5ノードN5に対応する血管分枝が対象血管に属するものと間違って抽出されたため、誤認部位である。
異なる例として、血管映像処理装置は、血管構造データ400に基づいてノード連結データ420から脱落したノードに対応する血管分枝を断絶部位として検出することができる。血管映像処理装置は、第1ノードN1~第5ノードN5に直接連結されているノード連結データ420を算出することができる。血管映像処理装置は、血管構造データ400に基づいてノード連結データ420で第2ノードN2が脱落(missing)されたことを検出できる。この場合、血管映像処理装置は、第2ノードN2に対応する血管分枝を断絶部位として検出することができる。
図5及び図6は、一実施形態に係るエラー部位を修正する方法を示す。
図5は、図4においてノード連結データ410が算出された場合の修正を説明する。一実施形態によれば、血管映像処理装置は、誤認部位(misidentification portion)が識別される場合に応答して、誤認部位に対応するノードを血管構造データにマッチングするノードに交替できる。例示的に、血管映像処理装置は、ノード連結データ410で第5ノードN5が血管構造データにマッチングしないものと判断することができる。血管映像処理装置は、血管構造データにマッチングしない第5ノードN5をノード連結データから排除し、血管構造データにマッチングする第6ノードN6をノード連結データに追加することで、修正されたノード連結データ510を生成することができる。したがって、血管映像処理装置は、交替された第6ノードN6に対応する血管分枝(例えば、図3で第6分枝306)を対象血管として抽出することができる。
図6は、図4においてノード連結データ420が算出された場合の修正を説明する。一実施形態によれば、血管映像処理装置は、断絶部位が識別される場合に応答して、断絶部位に対応する領域を連結することができる。例えば、血管映像処理装置は、断絶部位に対応するノードの間に連結可能なノードを挿入することができる。例示的に、血管映像処理装置は、ノード連結データ420で第1ノードN1と第5ノードN5との間が脱落したものと判断する。血管映像処理装置は、ノードN1、N5の間に第2ノードN2を挿入することで、修正されたノード連結データ610を生成することができる。血管映像処理装置は、第2ノードN2に対応する血管分枝(例えば、図3で第2分枝302)を対象血管として抽出することができる。更なる例として、血管映像処理装置は、断絶部位に対応するノードに対応する領域を連結する修正も可能である。図4において、第1ノードN1、第2ノードN2及び第5ノードN5が順次連結されるものと示されているが、これに限定されることなく、第1ノードN1及び第5ノードN5が連結されてもよく、第1ノードN1に対応する領域及び第5ノードN5に対応する領域が断絶されてもよい。血管映像処理装置は、第1ノードN1に対応する領域及び第5ノードN5に対応する領域が断絶した場合、第1ノードN1に対応する領域と第5ノードN5に対応する領域が分離されないように自然に連結することができる。
上述したように、ノード交替及びノードさらによる誤認部位修正を介して、血管映像処理装置は、血管構造データのトポロジーに基づいたノード連結性エラーを除去することができる。
図7は、対象血管の曲率情報を用いて対象血管のエラーを識別して修正する方法を示す。
一実施形態によれば、血管映像処理装置は、対象血管で分枝地点を基準にして互いに隣接する上位分枝と下位分枝との間の曲率を用いて誤認部位を識別して修正することができる。例えば、血管映像処理装置は、対象血管で分枝地点を基準にして互いに隣接する上位分枝と下位分枝との間の曲率が閾値曲率を超過する場合に応答して、下位分枝を誤認部位として判断することができる。例示的に、血管映像処理装置は、血管分枝701,702を対象血管として1次抽出した場合、血管分枝701,702間の曲率を算出する。血管映像処理装置は、算出された曲率が閾値曲率を超過する場合に応答して、上述した血管分枝701,702のうち、下位血管分枝702を対象血管に対する抽出結果において誤認部位として判断することができる。血管映像で実際に抽出しなければならない対象血管に属する互いに隣接している上位血管分枝と下位血管分枝との間には比較的に小さな曲率を示しているところ、閾値曲率を超過する曲率を示す互いに隣接する上位血管分枝及び下位血管分枝は、エラー部位を含む可能性がある。
血管映像処理装置は、対象血管で閾値曲率を用いてエラー部位が識別される場合に応答して、識別されたエラー部位を修正することができる。血管映像処理装置は識別されたエラー部位に対して、対象血管のうちエラー部位が含む分枝地点を基準にして区分される上位分枝及び下位分枝のうち、下位分枝を前記分枝との曲率が閾値曲率以下である候補下位分枝に交替することができる。候補下位分枝は、任意の分枝地点を基準にして1つの血管分枝に連結された複数の下位血管分枝のうち、対象血管に抽出されていない残りの1つ以上の下位分枝を示す。
例示的に、血管映像処理装置は、対象血管から抽出されていない候補下位血管分枝703と上位血管分枝701との間の曲率が閾値曲率以下である場合、下位血管分枝702の代わりに、候補下位血管分枝703を対象血管として抽出することができる。血管映像処理装置は、対象血管の抽出結果から下位血管分枝702を候補下位血管分枝703に交替することができる。さらに、血管映像処理装置は、識別されたエラー部位を修正することにおいて、血管の種類(例えば、左主冠状動脈(Left Main Coronary Artery、LM)、左前下行動脈(Left Anterior Descending Artery、LAD)、左回旋動脈(Left Circumflex Artery、LCX)、右冠状動脈(Right Coronary Artery、RCA))又は、血管領域(例えば、近位領域、中間領域、及び遠位領域)に応じて互いに異なる閾値曲率を適用できる。
血管映像処理装置は、曲率差を分析するために血管分枝の主方向成分分析(例えば、主成分分析(PCA、principal component analysis)を行ってもよい。血管映像処理装置は、分枝地点を基準にして区分された血管分枝それぞれに対して、主方向成分を分析することができる。血管分枝の主方向成分は、該当血管分枝の方向を代表する成分を示す。例えば、血管映像処理装置は、各血管分枝に対する主方向成分分析を用いて、各血管分枝ごとに主方向ベクトルを取得することができる。主方向ベクトルは、分枝地点から血管分枝が延長する方向を代表する方向成分を有するベクトルであって、単位ベクトルの大きさを有する。例えば、任意の血管分枝に対する主方向ベクトルは、分枝地点から該当の血管分枝に対応する地点に向かうベクトル成分の主成分の方向を示す。例示的に、血管分枝701,702,703の主方向ベクトルは、それぞれベクトル721,722,723に対応する。血管映像処理装置は、対象血管で分枝地点を基準にして上位分枝と下位分枝との間の主方向ベクトルの内積値が閾値を超過する場合に応答して、下位分枝を誤認部位として判断することができる。
例示的、血管映像処理装置は、対象血管で上位分枝701及び下位分枝702の主方向ベクトル721,722の内積値が閾値を超過する場合、血管映像処理装置は、上位分枝701及び下位分枝702が誤認部位(misidentification portion)を含むものと判断する。各血管分枝の主方向ベクトルは、単位ベクトルの大きさを有するため、主方向ベクトルの内積値はベクトルが互いになしている角度に依存する。結局、上位分枝に対応する主方向ベクトル及び下位分枝に対応する主方向ベクトルは、互いに反対方向に向かっていることから2個の主方向ベクトルの内積は負数値を有するため、主方向ベクトルの内積値が閾値以上であることは、血管分枝間の曲率が大きいことを意味する。血管映像で実際に抽出しなければならない対象血管は1つの血管として、互いに隣接する上位血管分枝と下位血管分枝との間に大きい曲率の変化がないことが一般的であるため、互いに隣接する上位血管分枝と下位血管分枝との間の主方向ベクトルの内積値が閾値を超過することは、血管映像処理装置が実際に探すべき対象血管ではない他の血管を抽出したものを示す。さらに、識別されたエラー部位の修正において、血管映像処理装置は、血管の種類又は血管領域に応じて互いに異なる閾値を適用することができる。
血管映像処理装置は、抽出された対象血管に対して分枝地点を基準にして分割された血管分枝ごとに主方向ベクトルを求めることができるが、予め指定された長さで血管を分割し、予め指定された長さで分割された断片(piece)ごとに主方向ベクトルを求めてもよい。血管映像処理装置は、抽出された対象血管に対して、上位断片及び下位断片の主方向を比較することで、曲率が閾値以上に大きくなった部分を捜し出す方式で対象血管でエラー部位を識別してもよい。上位断片及び下位断片は、血流の進行方向又は造影剤の進行方向により定義される。言い換えれば、血管映像処理装置は、抽出された対象血管を閾値の長さ以下に分割することができる。血管映像処理装置は、予め指定された長さに分割された互いに隣接する上位断片及び下位断片の主方向ベクトルの内積値が閾値を超過する場合に応答して、互いに隣接する上位断片及び下位断片のうちの1つを誤認部位として判断することができる。血管映像処理装置は、対象血管を閾値の長さ以下に分割し、上位断片と下位断片との間の曲率を算出できるため、分枝地点を基準にして分割された血管分枝により対象血管のエラー部位を識別する場合よりもさらに正確に対象血管のエラー部位を識別することができる。
血管映像処理装置は、対象血管で主方向ベクトルを用いてエラー部位が識別される場合に応答して、識別されたエラー部位を修正することができる。血管映像処理装置は、識別されたエラー部位に対して、対象血管のうちエラー部位が含む分枝地点を基準にして区分される上位分枝及び下位分枝のうち、下位分枝を前記分枝との主方向ベクトルの内積値が閾値以下である候補下位分枝に交替することができる。
例示的に、血管映像処理装置は、対象血管から抽出されていない候補下位血管分枝703と上位血管分枝701との間の主方向ベクトル723,721の内積値が閾値以下である場合、下位血管分枝702の代わりに候補下位血管分枝703を対象血管として抽出することができる。血管映像処理装置は、対象血管の抽出結果において、下位血管分枝702を候補下位血管分枝703に交替することができる。
図8は、対象血管の直径情報を用いて対象血管のエラーを識別して修正する方法を示す。
一実施形態に係る血管映像処理装置は、対象血管で分枝地点を基準にして互いに隣接する上位分枝と下位分枝との間の直径の差を算出し、誤認部位を識別して修正することができる。一実施形態によれば、血管映像処理装置は、対象血管で分枝地点を基準にして互いに隣接する上位分枝と下位分枝との間の直径の差が閾値以上である場合に応答して、下位分枝を誤認部位として判断することができる。血管映像処理装置は、モデルアンサンブルに基づいて分枝地点を基準にして識別された血管分枝それぞれに対して直径情報を取得することができる。例えば、図8に示された血管分枝801,802,803の直径は、それぞれr1、r2、r3に対応する。例示的に、血管映像処理装置は、血管分枝801,802を対象血管として1次抽出した場合、血管分枝801,802間の直径の差を算出することができる。血管映像処理装置は、算出された直径の差が閾値直径の差以上である場合に応答して、上述した血管分枝801,802のうち少なくともいずれかの分枝を対象血管に対する抽出結果から誤認部位として判断する。血管映像において、実際に抽出しなければならない対象血管に属する互いに隣接する上位血管分枝と下位血管分枝との間では、比較的に小径の変化を示しており、閾値直径の差以上の直径の差を示す互いに隣接する上位血管分枝及び下位血管分枝は、エラー部位を含んでいる可能性がある。他の一実施形態によれば、血管映像処理装置は、対象血管で互いに隣接する上位分枝と下位分枝との間の直径の差が閾値直径の差以上である場合に応答して、互いに隣接する上位分枝及び下位分枝に対して取得した直径情報をエラーとして判断することができる。言い換えれば、血管映像処理装置は、対象血管で分枝地点を基準にして互いに隣接する上位分枝と下位分枝との間の直径の差が閾値直径の差以上である場合、互いに隣接する上位分枝及び下位分枝のうちの1つをエラー部位として識別し、上位分枝及び下位分枝のうち少なくとも1つの血管分枝に対する直径情報が誤ったと判断することができる。即ち、血管映像処理装置が互いに隣接する上位分枝及び下位分枝のうち少なくとも1つの血管分枝に対して、血管分枝の実際の直径よりも広いか狭くその直径を間違って予測した場合である。例示的に、血管映像処理装置は、血管分枝801,802を対象血管から1次抽出した場合、上述した血管分枝801,802のうち少なくともいずれかの分枝に対する直径を間違って予測した場合として判断することができる。
血管映像処理装置は、抽出された対象血管に対して隣接する上位血管分枝と下位血管分枝との間の直径の差を算出するが、血管分枝の直径情報を平滑化(smoothing)して直径の差を算出してもよい。血管映像において、実際に抽出しなければならない対象血管に属する互いに隣接する上位分枝及び下位分枝の間では、比較的に小径の変化を示すことが一般的である。しかし、抽出しなければならない対象血管が疾患のある血管分枝を含んでいる場合、疾患のある血管分枝と連結されている血管分枝間で直径の変化は大きく示される。したがって、血管映像処理装置は、複数の血管分枝をグルーピングして血管分枝の直径情報を平滑化し直径の差を算出することで、抽出された対象血管から誤認部位を識別することができる。
血管映像処理装置は、対象血管で直径情報を用いてエラー部位が識別される場合に応答して、識別されたエラー部位を修正することができる。血管映像処理装置は、識別されたエラー部位に対して下位分枝を上位分枝との直径の差に基づいた交替及び対象血管に対する新しい直径情報を含む新しい抽出結果の生成のうち少なくとも1つを行うことで、識別されたエラー部位を修正することができる。
一実施形態に係る血管映像処理装置は、識別された誤認部位に対して誤認部位が含む分枝地点を基準にして区分される下位分枝を上位分枝との直径の差が閾値未満である候補下位分枝に交替できる。血流の進行方向又は造影剤の進行方向に応じて血管の直径が減少することが一般的である。したがって、一実施形態に係る血管映像処理装置は、下位分枝を上位分枝よりも直径が小さく、上位分枝との直径の差が閾値直径の差未満である候補下位分枝に交替することができる。例示的に、血管映像処理装置は、対象血管から抽出されていない候補下位血管分枝803と上位血管分枝801との間の直径の差が閾値直径の差未満である場合、下位血管分枝802の代わりに、候補下位血管分枝803を対象血管として抽出することができる。血管映像処理装置は、対象血管の抽出結果で下位血管分枝802を候補下位血管分枝803に交替することができる。
異なる例として、血管映像処理装置が対象血管で互いに隣接する上位分枝と下位分枝との間の直径の差が閾値直径の差以上である場合に応答して、互いに隣接する上位分枝及び下位分枝のうち、少なくとも1つに対して取得した直径情報をエラーとして判断することができる。血管映像処理装置は、対象血管に対して取得した直径情報がエラーとして判断された場合、対象血管に対する新しい抽出結果を生成することができる。即ち、血管映像処理装置は、互いに隣接する上位分枝及び下位分枝のうち少なくとも1つの血管分枝に対して、血管分枝の実際の直径よりも広いか狭く予測された場合として判断することができる。血管映像処理装置は、対象血管に対する新しい直径情報を取得するために、以前の抽出結果で適用されたモデルアンサンブルとは区別されるモデルアンサンブルに基づいて、対象血管に対して新しい直径情報を含む新しい抽出結果を生成することができる。さらに、血管映像処理装置は、全体血管に対する新しい直径情報を取得することも可能であるが、直径情報がエラーとして判断された血管分枝に対する新しい直径情報のみを新しく取得する修正も可能である。血管映像処理装置は、識別されたエラー部位の修正において、血管の種類(例えば、LM、LAD、LCX、RCA)又は血管領域に応じて、互いに異なる閾値直径の差を適用することができる。
図9は、対象血管の明度情報を用いて対象血管のエラーを識別して修正する方法を示す。
一実施形態によれば、血管映像処理装置は、対象血管で分枝地点を基準にして互いに隣接する上位分枝と下位分枝との間の明度(brightness)差を算出し、誤認部位を識別して修正することができる。例えば、血管映像処理装置は、対象血管で分枝地点を基準にして互いに隣接する上位分枝と下位分枝との間の明度差が閾値の明度差以上である場合に応答して、下位分枝を誤認部位として判断することができる。血管映像処理装置は、モデルアンサンブルに基づいて分枝地点を基準にし識別された血管分枝それぞれに対して明度情報を取得することができる。例えば、血管分枝の明度は、血管映像を取得するために投与した造影剤の分布濃度を示す。例示的に、血管映像処理装置は、血管分枝901,902を対象血管として1次抽出した場合、血管分枝901,902間の明度差を算出することができる。血管映像処理装置は、算出された明度差が閾値の明度差以上である場合に応答して、上述した血管分枝901,902のうち、下位血管分枝902を対象血管に対する抽出結果において誤認部位として判断することができる。血管映像において、実際に抽出しなければならない対象血管に属する互いに隣接する上位血管分枝と下位血管分枝との間では、比較的に小さな明度の変化を示しているところ、閾値の明度差以上の明度差を示す互いに隣接する上位血管分枝及び下位血管分枝は、エラー部位を含む可能性がある。他の一実施形態によれば、血管映像処理装置は、対象血管で互いに隣接する上位分枝と下位分枝との間の明度差が閾値の明度差以上である場合に応答して、互いに隣接する上位分枝及び下位分枝に対して取得した明度情報をエラーとして判断することができる。言い換えれば、血管映像処理装置は、対象血管で分枝地点を基準にして互いに隣接する上位分枝と下位分枝との間の明度差が閾値の明度差以上である場合、互いに隣接する上位分枝及び下位分枝のうちの1つをエラー部位として識別し、上位分枝及び下位分枝のうち少なくとも1つの血管分枝に対する明度情報が誤ったものと判断する。即ち、血管映像処理装置が互いに隣接する上位分枝及び下位分枝のうち少なくとも1つの血管分枝に対して、血管分枝の実際の明度よりも明るいか暗く、その明度を間違って予測した場合である。例示的に、血管映像処理装置は、血管分枝901,902を対象血管から1次抽出した場合、上述した血管分枝901,902のうち少なくともいずれかの分枝に対する直径を間違って予測した場合として判断する。
血管映像処理装置は、対象血管で明度情報を用いてエラー部位が識別される場合に応答して、識別されたエラー部位を修正することができる。一実施形態によれば、血管映像処理装置は、対象血管で互いに隣接する上位血管分枝と下位血管分枝との間の明度差が閾値の明度差以上である場合に応答して、下位血管分枝を誤認部位として判断することができる。血管映像処理装置は、識別された誤認部位に対して、下位分枝を上位分枝との明度差が閾値の明度差未満である候補下位分枝に交替することができる。例示的に、血管映像処理装置は、対象血管から抽出されていない候補下位血管分枝903と上位血管分枝901との間の明度差が閾値の明度差未満である場合、下位血管分枝902の代わりに、候補下位血管分枝903を対象血管として抽出することができる。血管映像処理装置は、対象血管の抽出結果で、下位血管分枝902を候補下位血管分枝903に交替することができる。さらに、血管映像処理装置は、識別されたエラー部位の修正において、血管の種類(例えば、LM、LAD、LCX、RCA)又は血管領域に応じて互いに異なる閾値直径の差を適用できる。
他の一実施形態によれば、血管映像処理装置が対象血管で互いに隣接する上位分枝と下位分枝との間の明度差が閾値の明度差以上である場合に応答して、互いに隣接する上位分枝及び下位分枝のうち、少なくとも1つに対して取得した明度情報をエラーとして判断する。血管映像処理装置は、対象血管に対し取得した明度情報がエラーと判断された場合、対象血管に対する新しい抽出結果を生成することができる。即ち、血管映像処理装置は、互いに隣接する上位分枝及び下位分枝のうち少なくとも1つの血管分枝に対して、血管分枝の実際の明度よりも明るいか暗く予測した場合であると判断する。血管映像処理装置は、対象血管に対する新しい明度情報を取得するために、以前の抽出結果で適用されたモデルアンサンブルとは区別されるモデルアンサンブルに基づいて、対象血管に対して新しい直径情報を含む新しい抽出結果を生成することができる。さらに、血管映像処理装置は、全体血管に対する新しい明度情報を取得することも可能であるが、明度情報がエラーとして判断された血管分枝に対する新しい明度情報のみを新しく取得する修正も可能である。
一実施形態に係る血管映像処理装置は、分枝地点を含む領域だけでなく、分枝地点を含んでいない領域に対するエラー部位識別及び修正も可能である。例えば、血管映像処理装置は、1つの血管分枝内で発生するエラー部位に対しても修正が可能である。1つの血管分枝に対応する領域が空間的に連結された領域ではない、分離されたり離隔される断絶部位を示す。血管映像処理装置は、1つの血管分枝に対応する領域が断絶した場合、1つの血管分枝に対応する領域が互いに分離しないように自然に連結する修正が可能である。更なる例として、血管映像処理装置は、1つの血管分枝内で血管の直径変化が閾値直径以上である場合、血管分枝の直径情報をエラーとして判断して血管分枝に対する新しい直径情報を新しく取得する修正も可能である。同様に、血管映像処理装置は、1つの血管分枝内で血管の明度変化が閾値の明度以上である場合、血管分枝の明度情報をエラーとして判断し、血管分枝に対する新しい明度情報を新しく取得する修正も可能である。
一実施形態に係る血管映像処理装置は、エラー部位が識別される場合に応答して、以前の抽出結果で適用されたモデルアンサンブルとは区別されるモデルアンサンブルに基づいて、対象血管に対して新しい抽出結果を生成することができる。エラー部位が識別される場合、血管映像処理装置は、直接エラー部位を修正せずに、格納された複数の機械学習モデルのうち、以前の抽出結果で適用されたモデルアンサンブルとは区別される1つのモデルアンサンブルに基づいて、対象血管に対して新しい抽出結果を生成することができる。
また、血管映像処理装置は、格納された複数の機械学習モデルのうち、2以上のモデルアンサンブルに基づいて対象血管に対して抽出結果をそれぞれ生成することができる。血管映像処理装置は、モデルアンサンブルにより生成された対象血管の抽出結果ごとに対象血管に関する血管構造データ、対象血管の曲率情報、対象血管の直径情報、及び対象血管の明度情報のうち少なくとも1つに基づいてエラー部位を識別することができる。血管映像処理装置は、生成された複数の対象血管に対する抽出結果のうち、実際に抽出しなければならない対象血管と最も近い1つの抽出結果を対象血管として自動選択する。例えば、血管映像処理装置は、識別されたエラー部位が最も少ない対象血管に対する抽出結果を実際に抽出しなければならない対象血管と最も近い抽出結果として選択することができるが、これに限定されることはない。さらに、血管映像処理装置は、複数の対象血管に対する抽出結果に関する情報をユーザに提供することができる。血管映像処理装置は、ユーザから複数の対象血管に対する抽出結果のうちの1つが入力され、入力された抽出結果を対象血管として選択することができる。
一実施形態に係る血管映像処理装置は、対象血管の抽出結果でエラー部位が識別される場合に応答して、候補血管分枝それぞれに対する連結性スコア(connectivity score)に基づいて候補血管分枝のうち1つの血管分枝を対象血管に抽出することができる。血管映像処理装置は、対象血管の抽出結果からエラー部位が識別される場合、エラー部位を修正できる候補血管分枝のうち1つの血管分枝を選択する。血管映像処理装置は、候補血管分枝それぞれに対して連結性スコアを算出する。連結性スコアは、対象血管の抽出結果で任意の候補血管分枝と連結されている血管分枝と該当の候補血管分枝との間の連結性を指示するスコアとして、例えば、候補血管分枝及び候補血管分枝と連結された血管分枝間の血管構造データに対するマッチング程度、曲率差、直径の差、及び明度差などに基づいて算出され得る。連結性スコアは、候補血管分枝が対象血管に含まれている場合、修正された対象血管が実際に抽出しなければならない対象血管に類似の程度に対応する。
例えば、候補血管分枝が対象血管から抽出された場合、候補血管分枝と前記候補血管分枝と連結されている血管分枝の間で対象血管に関する血管構造データ、前記対象血管の曲率情報、前記対象血管の直径情報、及び前記対象血管の明度情報のうち少なくとも1つに基づいて連結性スコアを算出することができる。しかし、これに限定されることなく、様々な方法で候補血管分枝それぞれに対する連結性スコアを算出してもよい。血管映像処理装置は、連結性スコアに基づいて候補血管分枝のうち1つの血管分枝を選択することができる。例えば、候補血管分枝のうち連結性スコアが最も高い血管分枝として、エラーが最も小さいと判断される血管分枝を選択する。血管映像処理装置は、選択された血管分枝を対象血管に抽出することができる。
例示的に、図5において、血管映像処理装置が対象血管の抽出結果から第1ノードN1、第3ノードN3、及び第5ノードN5が連結されているノード連結データ410を算出した場合、血管映像処理装置は、ノード連結データ410で第3ノードN3及び第5ノードN5間の連結が血管構造データ400にマッチングしないものと決定することができる。血管映像処理装置は、連結エラーに対応するノードN3、N5のうち、少なくともいずれかのノードに対応する血管分枝をエラー部位として判断する。血管映像処理装置は、第3ノードN3に対応する血管分枝をエラーとして判断した場合、血管構造データにマッチングしない第3ノードN3をノード連結データから排除し、血管構造データにマッチングする第2ノードN2をノード連結データに追加する。また、血管映像処理装置は、第5ノードN5に対応する血管分枝をエラーとして判断した場合、血管構造データにマッチングしない第5ノードN3をノード連結データから排除し、血管構造データにマッチングする第6ノードN6又は第7ノードN7をノード連結データに追加する。結局、血管映像処理装置は、対象血管のノードになり得る第2ノード、第6ノード、第7ノードに対応する候補血管分枝それぞれに対する連結性スコアを算出することができる。血管映像処理装置は、候補血管分枝それぞれに対する連結性スコアに基づいて、エラーの最も小さいと判断される血管分枝を選択して対象血管に抽出することができる。
更なる例示として、図7において、血管映像処理装置が血管分枝701,702を対象血管にして1次抽出した場合、血管分枝701,702間の曲率が閾値曲率を超過してエラー部位として判断してもよい。血管映像処理装置は、互いに隣接する上位分枝及び下位分枝のうち、下位分枝を上位分枝との曲率が閾値曲率以下である候補下位分枝に交替してもよい。図7とは異なって、候補下位血管分枝が複数である場合、血管映像処理装置は、候補血管分枝それぞれに対して連結性スコアを算出してもよい。血管映像処理装置は、候補血管分枝それぞれに対する連結性スコアに基づいて、エラーが最も小さいと判断される血管分枝を選択し対象血管に抽出することができる。図8及び図9においても、血管映像処理装置は上述と同じ方法で、候補血管分枝のうちの1つを対象血管として選択して抽出することができる。
図10~図11は、一実施形態に係るエラー部位を修正する方法を示す。
図10は、識別されたエラー部位が誤認部位である場合、ユーザによる修正について説明する。一実施形態によれば、血管映像処理装置は、対象血管でエラー部位が識別される場合、ユーザの入力に基づいて識別されたエラー部位を修正する。血管映像処理装置は、対象血管で識別されたエラー部位が誤認部位である場合に応答して、1つ以上の候補分枝をユーザに提供する。血管映像処理装置は、ユーザから1つ以上の候補分枝のうち1つの血管分枝に対する選択入力を受信する場合に応答して、エラー部位に対応する分枝を選択された分枝に交替することができる。例示的に、血管映像処理装置が血管分枝1001,1002を対象血管として1次抽出した場合、対象血管に関する血管構造データ、対象血管の曲率情報、対象血管の直径情報、及び対象血管の明度情報のうち少なくとも1つに基づいて血管分枝1001,1002の少なくともいずれか1つの分枝を対象血管に対する抽出結果からエラー部位として判断することができる。血管映像処理装置は、識別されたエラー部位が誤認部位である場合、対象血管から抽出されていない1つ以上の候補血管分枝1010,1020,1030をユーザに提供する。血管映像処理装置の入出力インターフェースは、ディスプレイを介して候補血管分枝1010,1020,1030に対する選択をユーザに提供可能なグラフィックオブジェクト1071,1072,1073を出力する。グラフィックオブジェクトは、候補血管分枝に対応する領域を取り囲む形態であってもよい。図10において、グラフィックオブジェクト1071,1072,1073は長方状に図示されているが、形態がこれに限定されることはない。ユーザは、単純にディスプレイに出力される候補血管分枝に対応するグラフィックオブジェクトのうちの1つに対応する地点に対して、クリックしたりタッチして候補血管分枝を選択してもよい。即ち、血管映像処理装置は、グラフィックオブジェクトに対応する地点に対してユーザ入力を検出する場合、グラフィックオブジェクトに対応する候補血管分枝が選択されたものと判断し、エラー部位に対応する血管分枝を選択された候補血管分枝に交替することができる。例示的に、血管映像処理装置がユーザからグラフィックオブジェクト1071に対応する地点に対するユーザの入力を検出する場合、グラフィックオブジェクト1071に対応する候補血管分枝1002を対象血管として抽出することができる。即ち、血管映像処理装置は、対象血管の抽出結果で血管分枝1002を血管分枝1010に交替することができる。
図11は、識別された誤り部位が断絶部位である場合、ユーザによる修正について説明する。一実施形態によると、血管映像処理装置は、対象血管でエラー部位が識別される場合、ユーザの入力に基づいて識別されたエラー部位を修正することができる。血管映像処理装置は、対象血管で識別されたエラー部位が断絶部位である場合に応答し、断絶部位を含む血管領域をユーザに提供する。血管映像処理装置は、修正結果に対して複数の地点からユーザ入力が検出される場合に応答して、ユーザ入力の開始地点に対応する血管分枝に対応する領域及びユーザ入力の終了地点に対応する血管分枝に対応する領域を連結することができる。
一実施形態に係る血管映像処理装置は、ユーザから受信されるドラッグ入力に応答して、ドラッグ入力によって選択された断絶部位に対応する分枝を連結することができる。例示的に、血管映像処理装置が血管分枝1101,1102を対象血管に1次抽出した場合、血管映像処理装置は、抽出された対象血管の血管分枝1101,1102に対応する領域が互いに分離する断絶部位として判断することができる。血管映像処理装置の入出力インターフェースは、ディスプレイを介して血管分枝1101,1102に対応する領域間の連結に関するグラフィックオブジェクト1170を出力する。グラフィックオブジェクト1170は、断絶部位に対応する血管分枝1101,1102の少なくとも一部に対応する領域を取り囲む形態として、血管分枝1101,1102がユーザによって連結可能であることを示す。血管映像処理装置は、ユーザから受信されるドラッグ入力に応答して、ドラッグ入力によって選択された血管分枝1101,1102に対応する領域を連結し、断絶した血管分枝1101,1102間の血管分枝1103を対象血管として抽出することができる。ドラッグ入力によって選択された分枝は、ドラッグ開始地点に対応する血管分枝及びドラッグ終了地点に対応する血管分枝を示す。即ち、血管映像処理装置は、対象血管の抽出結果で血管分枝1103を含ませることができる。
他の一実施形態に係る血管映像処理装置は、ユーザから受信される2以上の地点に対応するポインティング(pointing)入力に応答して、ポインティング入力の開始地点と終了地点に対応する血管分枝の間の領域を連結することができる。ポインティング入力は、ディスプレイ上の一地点を指示する入力として、例えば、マウス操作によるクリック入力及び/又はタッチディスプレイに対するタッチ操作による入力(例えば、タッチ入力)を含んでもよい。但し、ポインティング入力がこれに限定されることなく、設計に応じて、様々なポインティングデバイス(例えば、トラックボールマウス、タッチパッド、及びトラックパッドなど)による位置指示入力を含んでもよい。即ち、血管映像処理装置は、ユーザから開始地点及び終了地点に対するクリック又はタッチ入力を受信することができる。血管映像処理装置は、ユーザから受信されるポインティング入力に応答して、ポインティング入力の開始地点に対応する血管分枝に対応する領域とポインティング入力の終了地点に対応する血管分枝に対応する領域とを連結することができる。
図12は、一実施形態に係る血管映像処理装置を概略的に示すブロック図である。
一実施形態に係る血管映像を処理するシステム1200は、血管映像処理装置1210及び血管映像撮影装置1220を含む。血管映像処理装置1210は、映像受信部1211、プロセッサ1212、及び入出力インターフェース1213を含む。映像受信部1211は、血管映像撮影装置1220によって撮影された血管映像を受信する。プロセッサ1212は、機械学習モデルを用いて映像受信部から受信した血管映像から対象血管を抽出し、エラー部位を識別し、識別されたエラー部位を修正することができる。プロセッサ1212の動作がこれに限定されることなく、プロセッサ1212は、図1~図11を参照して上述した動作を行うことができる。
入出力インターフェース1213は、ユーザから入力を受信してプロセッサに伝達する。例えば、入出力インターフェース1213は、マウス操作、及びタッチ操作などによる入力を受信してもよい。また、入出力インターフェース1213は、視覚的なフィードバックをユーザに提供することができる。例えば、入出力インターフェース1213は、対象血管の抽出結果をステップごとにディスプレイを通じて出力することができる。
以上で説明された実施形態は、ハードウェア構成要素、ソフトウェア構成要素、又はハードウェア構成要素及びソフトウェア構成要素の組み合せで具現される。例えば、本実施形態で説明した装置及び構成要素は、例えば、プロセッサ、コントローラ、ALU(arithmetic logic unit)、デジタル信号プロセッサ(digital signal processor)、マイクロコンピュータ、FPA(field programmable array)、PLU(programmable logic unit)、マイクロプロセッサー、又は命令(instruction)を実行して応答する異なる装置のように、1つ以上の汎用コンピュータ又は特殊目的コンピュータを用いて具現される。処理装置は、オペレーティングシステム(OS)及びオペレーティングシステム上で実行される1つ以上のソフトウェアアプリケーションを実行する。また、処理装置は、ソフトウェアの実行に応答してデータをアクセス、格納、操作、処理、及び生成する。理解の便宜のために、処理装置は1つが使用されるものとして説明する場合もあるが、当技術分野で通常の知識を有する者は、処理装置が複数の処理要素(processing element)及び/又は複数類型の処理要素を含むことが把握する。例えば、処理装置は、複数のプロセッサ又は1つのプロセッサ及び1つのコントローラを含む。また、並列プロセッサ(parallel processor)のような、他の処理構成も可能である。
ソフトウェアは、コンピュータプログラム、コード、命令、又はそのうちの一つ以上の組合せを含み、希望の通りに動作するよう処理装置を構成したり、独立的又は結合的に処理装置を命令することができる。ソフトウェア及び/又はデータは、処理装置によって解釈されたり処理装置に命令又はデータを提供するために、いずれかの類型の機械、構成要素、物理的装置、仮想装置、コンピュータ格納媒体又は装置、又は送信される信号波に永久的又は一時的に具体化することができる。ソフトウェアはネットワークに連結されたコンピュータシステム上に分散され、分散した方法で格納されたり実行され得る。ソフトウェア及びデータは一つ以上のコンピュータで読出し可能な記録媒体に格納され得る。
本実施形態による方法は、様々なコンピュータ手段を介して実施されるプログラム命令の形態で具現され、コンピュータ読み取り可能な記録媒体に記録される。記録媒体は、プログラム命令、データファイル、データ構造などを単独又は組み合せて含む。記録媒体及びプログラム命令は、本発明の目的のために特別に設計して構成されたものでもよく、コンピュータソフトウェア分野の技術を有する当業者にとって公知のものであり使用可能なものであってもよい。コンピュータ読み取り可能な記録媒体の例として、ハードディスク、フロッピー(登録商標)ディスク及び磁気テープのような磁気媒体、CD-ROM、DVDのような光記録媒体、フロプティカルディスクのような磁気-光媒体、及びROM、RAM、フラッシュメモリなどのようなプログラム命令を保存して実行するように特別に構成されたハードウェア装置を含む。プログラム命令の例としては、コンパイラによって生成されるような機械語コードだけでなく、インタプリタなどを用いてコンピュータによって実行される高級言語コードを含む。上記で説明したハードウェア装置は、本発明に示す動作を実行するために1つ以上のソフトウェアモジュールとして作動するように構成してもよく、その逆も同様である。
上述したように実施形態をたとえ限定された図面によって説明したが、当技術分野で通常の知識を有する者であれば、上記の説明に基づいて様々な技術的な修正及び変形を適用することができる。例えば、説明された技術が説明された方法と異なる順に実行され、及び/又は説明されたシステム、構造、装置、回路などの構成要素が説明された方法とは異なる形態に結合又は組み合わせられてもよく、他の構成要素又は均等物によって置き換え又は置換されたとしても適切な結果を達成することができる。

Claims (13)

  1. プロセッサによって実行される血管映像処理方法であって、
    血管映像から対象血管を抽出するステップと、
    前記対象血管の抽出結果で前記対象血管に関する血管構造データ、前記対象血管の曲率情報、前記対象血管の直径情報、及び前記対象血管の明度情報のうち少なくとも1つに基づいてエラー部位を識別するステップと、
    前記対象血管の抽出結果でエラー部位が識別される場合に応答して、前記識別されたエラー部位を修正するステップと、
    を含む血管映像処理方法。
  2. 前記エラー部位を識別するステップは、
    前記血管映像から血管領域及び非血管領域を分割するステップと、
    前記血管領域で分枝地点を基準にして区分される血管分枝のそれぞれに対応するノードに基づいて前記血管構造データを生成するステップと、
    前記対象血管のノード間の連結データを前記血管構造データに比較するステップと、
    前記血管構造データのトポロジーに基づいてノード間の連結エラーが検出された場合、前記連結エラーに対応するノードを前記エラー部位として判断するステップと、
    を含む、請求項1に記載の血管映像処理方法。
  3. 前記エラー部位を修正するステップは、
    前記エラー部位が断絶部位として識別される場合に応答して、前記断絶部位に対応する領域を連結するステップと、
    前記エラー部位が誤認部位として識別される場合に応答して、前記誤認部位に対応するノードを血管構造データにマッチングするノードに交替するステップと、
    を含む、請求項1に記載の血管映像処理方法。
  4. 前記エラー部位を識別するステップは、前記対象血管の抽出結果で互いに隣接する上位分枝と下位分枝との間の曲率が閾値曲率を超過する場合に応答して、前記下位分枝を誤認部位として識別するステップを含む、請求項1に記載の血管映像処理方法。
  5. 前記エラー部位を修正するステップは、前記識別された誤認部位に対して、前記下位分枝を前記上位分枝との曲率が臨界曲率以下である他の候補下位分枝に交替するステップを含む、請求項4に記載の血管映像処理方法。
  6. 前記エラー部位を識別するステップは、前記対象血管の抽出結果で互いに隣接する上位分枝と下位分枝との間の直径の差が閾値以上である場合に応答して、前記下位分枝をエラー部位として識別するステップを含む、請求項1に記載の血管映像処理方法。
  7. 前記エラー部位を修正するステップは、前記識別されたエラー部位に対して、前記下位分枝を前記上位分枝との直径の差に基づいた交替、及び前記対象血管に対する新しい直径情報を含む新しい抽出結果の生成のうち少なくとも1つを行うことによって前記識別されたエラー部位を修正するステップを含む、請求項6に記載の血管映像処理方法。
  8. 前記エラー部位を識別するステップは、前記対象血管の抽出結果で互いに隣接する上位分枝と下位分枝との間の明度差が閾値の明度差以上である場合に応答して、前記下位分枝を誤認部位として識別するステップを含む、請求項1に記載の血管映像処理方法。
  9. 前記エラー部位を修正するステップは、前記識別されたエラー部位に対して、前記下位分枝を前記上位分枝との明度差が閾値の明度差未満である他の候補下位分枝に交替するステップを含む、請求項8に記載の血管映像処理方法。
  10. 前記エラー部位を修正するステップは、エラー部位が識別された場合に応答して、前記対象血管に対する新しい抽出結果を生成するステップを含む、請求項1に記載の血管映像処理方法。
  11. 前記エラー部位を修正するステップは、
    前記エラー部位が誤認部位である場合に応答して、1つ以上の候補分枝をユーザに提供するステップと、
    前記ユーザから前記1つ以上の候補分枝のいずれかの分枝に対する選択入力を受信する場合に応答して、前記エラー部位に対応する分枝を前記選択された分枝に交替するステップを含む、請求項1に記載の血管映像処理方法。
  12. 前記エラー部位を修正するステップは、複数の地点でユーザ入力が検出される場合に応答して、前記ユーザ入力の開始地点に対応する血管分枝に対応する領域及び前記ユーザ入力の終了地点に対応する血管分枝に対応する領域を連結するステップを含む、請求項1に記載の血管映像処理方法。
  13. 血管映像処理装置であって、
    血管映像を受信する映像受信部と、
    前記血管映像から対象血管を抽出し、前記対象血管の抽出結果で前記対象血管に関する血管構造データ、前記対象血管の曲率情報、前記対象血管の直径情報、及び前記対象血管の明度情報のうち少なくとも1つに基づいてエラー部位を識別し、前記対象血管の抽出結果でエラー部位が識別される場合に応答して、前記識別されたエラー部位を修正するプロセッサを含む、血管映像処理装置。
JP2022561515A 2020-06-02 2021-02-05 血管映像を自動処理する方法及び装置 Active JP7427293B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
KR10-2020-0066267 2020-06-02
KR1020200066267A KR102395873B1 (ko) 2020-06-02 2020-06-02 자동으로 혈관 영상을 처리하는 방법 및 장치
PCT/KR2021/001538 WO2021246612A1 (ko) 2020-06-02 2021-02-05 자동으로 혈관 영상을 처리하는 방법 및 장치

Publications (2)

Publication Number Publication Date
JP2023521773A true JP2023521773A (ja) 2023-05-25
JP7427293B2 JP7427293B2 (ja) 2024-02-05

Family

ID=78831304

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022561515A Active JP7427293B2 (ja) 2020-06-02 2021-02-05 血管映像を自動処理する方法及び装置

Country Status (5)

Country Link
US (1) US12039729B2 (ja)
EP (1) EP4159131A1 (ja)
JP (1) JP7427293B2 (ja)
KR (2) KR102395873B1 (ja)
WO (1) WO2021246612A1 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102395873B1 (ko) 2020-06-02 2022-05-10 주식회사 메디픽셀 자동으로 혈관 영상을 처리하는 방법 및 장치
KR102521660B1 (ko) * 2020-11-30 2023-04-14 주식회사 메디픽셀 복수의 예측 결과를 활용한 혈관 영상 추출 방법 및 장치
KR102479645B1 (ko) * 2022-03-31 2022-12-23 재단법인 아산사회복지재단 병변 결정 방법 및 장치

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004283373A (ja) * 2003-03-20 2004-10-14 Toshiba Corp 管腔状構造体の解析処理装置
JP2008525126A (ja) * 2004-12-22 2008-07-17 バイオ−ツリー システムズ, インコーポレイテッド 病気の診断およびモニタリングのための医療画像化方法および装置ならびにその用途
JP2011212314A (ja) * 2010-03-31 2011-10-27 Fujifilm Corp 血管表示制御装置および方法、並びにプログラム
JP2013090799A (ja) * 2011-10-26 2013-05-16 Fujifilm Corp 画像処理装置および方法並びにプログラム
JP2014180492A (ja) * 2013-03-21 2014-09-29 Canon Inc 被検体情報取得装置およびその制御方法
US20140355858A1 (en) * 2013-06-03 2014-12-04 University Of Florida Research Foundation, Inc. Vascular anatomy modeling derived from 3-dimensional medical image processing
JP2016533815A (ja) * 2013-08-16 2016-11-04 ハートフロー, インコーポレイテッド 患者に特異的な解剖学的データからパーソナライズされた血管移植を特定するためのシステム及び方法
JP2019512295A (ja) * 2016-03-04 2019-05-16 4ディーエックス リミテッド 撮像方法およびシステム
JP2019122762A (ja) * 2018-01-11 2019-07-25 キヤノンメディカルシステムズ株式会社 医用情報処理装置、医用情報処理システム及び医用情報処理方法
JP2019521733A (ja) * 2016-05-16 2019-08-08 キャスワークス リミテッド 画像からの血管の選択方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101551862B (zh) 2009-05-13 2011-06-01 西安电子科技大学 基于边界距离场的血管中心路径提取方法
CN102404595B (zh) 2011-08-16 2013-11-20 上海交通大学 可提供3d节目拍摄指导的极线校正方法
KR101284865B1 (ko) * 2011-09-09 2013-07-09 삼성중공업 주식회사 혈관 인식 장치, 그 방법 및 기록 매체
US9814433B2 (en) 2012-10-24 2017-11-14 Cathworks Ltd. Creating a vascular tree model
JP6173751B2 (ja) * 2013-04-09 2017-08-02 東芝メディカルシステムズ株式会社 医用画像処理装置、x線診断装置及び医用画像処理プログラム
KR101625955B1 (ko) * 2014-11-03 2016-06-01 재단법인 아산사회복지재단 장기의 동맥 및 정맥의 구분 방법
US10037603B2 (en) * 2015-05-04 2018-07-31 Siemens Healthcare Gmbh Method and system for whole body bone removal and vascular visualization in medical image data
DE102017200489A1 (de) 2017-01-13 2018-07-19 Siemens Healthcare Gmbh Verfahren zur Reduzierung von aufgrund von Gefäßüberlagerungen auftretenden Artefakten in einem vierdimensionalen Angiographiedatensatz, Bildverarbeitungseinrichtung, Computerprogramm und elektronisch lesbarer Datenträger
CN108205807A (zh) * 2017-12-29 2018-06-26 四川和生视界医药技术开发有限公司 视网膜血管图像编辑方法及视网膜血管图像编辑装置
KR102050649B1 (ko) * 2018-01-08 2019-12-17 숭실대학교산학협력단 2차원 x-선 조영영상의 혈관 구조 추출 방법, 이를 수행하기 위한 기록매체 및 장치
KR102014385B1 (ko) * 2018-02-20 2019-08-26 (주)휴톰 수술영상 학습 및 학습 기반 수술동작 인식 방법 및 장치
KR102395873B1 (ko) 2020-06-02 2022-05-10 주식회사 메디픽셀 자동으로 혈관 영상을 처리하는 방법 및 장치

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004283373A (ja) * 2003-03-20 2004-10-14 Toshiba Corp 管腔状構造体の解析処理装置
JP2008525126A (ja) * 2004-12-22 2008-07-17 バイオ−ツリー システムズ, インコーポレイテッド 病気の診断およびモニタリングのための医療画像化方法および装置ならびにその用途
JP2011212314A (ja) * 2010-03-31 2011-10-27 Fujifilm Corp 血管表示制御装置および方法、並びにプログラム
JP2013090799A (ja) * 2011-10-26 2013-05-16 Fujifilm Corp 画像処理装置および方法並びにプログラム
JP2014180492A (ja) * 2013-03-21 2014-09-29 Canon Inc 被検体情報取得装置およびその制御方法
US20140355858A1 (en) * 2013-06-03 2014-12-04 University Of Florida Research Foundation, Inc. Vascular anatomy modeling derived from 3-dimensional medical image processing
JP2016533815A (ja) * 2013-08-16 2016-11-04 ハートフロー, インコーポレイテッド 患者に特異的な解剖学的データからパーソナライズされた血管移植を特定するためのシステム及び方法
JP2019512295A (ja) * 2016-03-04 2019-05-16 4ディーエックス リミテッド 撮像方法およびシステム
JP2019521733A (ja) * 2016-05-16 2019-08-08 キャスワークス リミテッド 画像からの血管の選択方法
JP2019122762A (ja) * 2018-01-11 2019-07-25 キヤノンメディカルシステムズ株式会社 医用情報処理装置、医用情報処理システム及び医用情報処理方法

Also Published As

Publication number Publication date
KR20220061940A (ko) 2022-05-13
KR102395873B9 (ko) 2022-09-30
KR102562317B1 (ko) 2023-08-03
US12039729B2 (en) 2024-07-16
KR102395873B1 (ko) 2022-05-10
WO2021246612A1 (ko) 2021-12-09
KR102562317B9 (ko) 2023-10-12
KR20210149350A (ko) 2021-12-09
JP7427293B2 (ja) 2024-02-05
EP4159131A1 (en) 2023-04-05
US20230196559A1 (en) 2023-06-22

Similar Documents

Publication Publication Date Title
JP7427293B2 (ja) 血管映像を自動処理する方法及び装置
Chen et al. Curve-like structure extraction using minimal path propagation with backtracking
US9471989B2 (en) Vascular anatomy modeling derived from 3-dimensional medical image processing
US11216950B2 (en) Method and system for automatically segmenting blood vessel in medical image by using machine learning and image processing algorithm
CN109377458B (zh) 一种冠脉分割断裂的修复方法及装置
KR102521660B1 (ko) 복수의 예측 결과를 활용한 혈관 영상 추출 방법 및 장치
US20230351550A1 (en) Method and device for processing blood vessel image on basis of user input
WO2008097520A2 (en) System and method for computer aided detection of pulmonary embolism in tobogganing in ct angiography
US20120002855A1 (en) Stent localization in 3d cardiac images
JP2011098195A (ja) 構造物検出装置および方法ならびにプログラム
KR102119057B1 (ko) 골절 진단모델의 학습 방법 및 장치
CN110176013B (zh) 血管分割方法、装置、电子设备及存储介质
CN111178420B (zh) 一种二维造影图像上冠脉段标注方法及***
CN114708259B (zh) 基于cta的头颈血管狭窄检测方法、装置、设备及介质
CN111507981B (zh) 图像处理方法和装置、电子设备、计算机可读存储介质
JP7479080B2 (ja) 血管分割に基づいて血管映像を処理する方法及び装置
Bacchuwar et al. Scale-space for empty catheter segmentation in PCI fluoroscopic images
KR102542972B1 (ko) 3차원 혈관 구조를 생성하는 방법 및 장치
Ghanaei et al. Automatic graph-based method for classification of retinal vascular bifurcations and crossovers
Qureshi et al. A probabilistic model for the optimal configuration of retinal junctions using theoretically proven features
CN115131508A (zh) 基于数据处理的dsa建模点云数据融合处理方法
Fang et al. Greedy Graph Searching for Vascular Tracking in Angiographic Image Sequences

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20221013

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230905

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231030

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20231030

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231204

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240117

R150 Certificate of patent or registration of utility model

Ref document number: 7427293

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150