JP2023502409A - 血管中心線に基づき血管輪郭線を取得するステップと方法及び装置 - Google Patents
血管中心線に基づき血管輪郭線を取得するステップと方法及び装置 Download PDFInfo
- Publication number
- JP2023502409A JP2023502409A JP2022529118A JP2022529118A JP2023502409A JP 2023502409 A JP2023502409 A JP 2023502409A JP 2022529118 A JP2022529118 A JP 2022529118A JP 2022529118 A JP2022529118 A JP 2022529118A JP 2023502409 A JP2023502409 A JP 2023502409A
- Authority
- JP
- Japan
- Prior art keywords
- vessel
- blood vessel
- contour
- vascular
- centerline
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 132
- 238000000605 extraction Methods 0.000 claims abstract description 27
- 210000004351 coronary vessel Anatomy 0.000 claims abstract description 25
- 230000002792 vascular Effects 0.000 claims description 84
- 230000011218 segmentation Effects 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 12
- 238000010586 diagram Methods 0.000 claims description 12
- 230000037361 pathway Effects 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 230000003321 amplification Effects 0.000 claims description 5
- 238000002586 coronary angiography Methods 0.000 claims description 5
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 5
- 230000002708 enhancing effect Effects 0.000 claims 1
- 230000003313 weakening effect Effects 0.000 claims 1
- 230000006870 function Effects 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 201000000057 Coronary Stenosis Diseases 0.000 description 2
- 206010011089 Coronary artery stenosis Diseases 0.000 description 2
- 206010057469 Vascular stenosis Diseases 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000002107 myocardial effect Effects 0.000 description 2
- 208000031225 myocardial ischemia Diseases 0.000 description 2
- 230000002966 stenotic effect Effects 0.000 description 2
- 206010002383 Angina Pectoris Diseases 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013130 cardiovascular surgery Methods 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000000302 ischemic effect Effects 0.000 description 1
- 150000002632 lipids Chemical class 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 210000004165 myocardium Anatomy 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- 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
-
- 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
-
- 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
- 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/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/10—Segmentation; Edge detection
- G06T7/181—Segmentation; Edge detection involving edge growing; involving edge linking
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/66—Analysis of geometric attributes of image moments or centre of gravity
-
- 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/20021—Dividing image into blocks, subimages or windows
-
- 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/20112—Image segmentation details
- G06T2207/20156—Automatic seed setting
-
- 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/30048—Heart; Cardiac
-
- 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)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Vascular Medicine (AREA)
- Optics & Photonics (AREA)
- High Energy & Nuclear Physics (AREA)
- Physiology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Cardiology (AREA)
- Quality & Reliability (AREA)
- Geometry (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Image Analysis (AREA)
Abstract
Description
冠状動脈二次元造影画像に基づき血管中心線を抽出するステップと、
前記血管中心線に基づき、真っすぐに伸ばした血管画像を取得するステップと、
前記真っすぐに伸ばした血管画像上において、血管直径閾値D閾を設定するステップと、
前記D閾に基づき、血管中心直線の両側に血管プリセット輪郭線を生成するステップと、
前記血管プリセット輪郭線を前記血管中心直線に段階的に寄せ、真っすぐに伸ばした後の血管の輪郭線を取得するステップと、
前記真っすぐに伸ばした後の血管の輪郭線を前記血管中心線の画像上に投射し直し、血管輪郭線を取得するステップと、
を含む、
血管中心線に基づき血管輪郭線を取得する方法。
前記血管中心線に基づき、真っすぐに伸ばした血管画像を取得するステップは、
前記血管中心線を真っすぐに伸ばし、血管中心直線を取得するステップと、
前記始点から前記終点までの血管の延在方向に沿って、前記部分的血管区域図をx個のユニットに分けるステップ(但しxは正の整数である)と、
各前記ユニットの血管中心線を前記血管中心直線に沿って対応して設けるステップと、
対応して設けた後の画像が前記真っすぐに伸ばした血管画像になるステップと、
を含む、
上記の血管中心線に基づき血管輪郭線を取得する方法。
前記血管プリセット輪郭線を前記血管中心直線に段階的に寄せ、真っすぐに伸ばした後の血管の輪郭線を取得するステップは、
前記血管プリセット輪郭線をy個のユニットに分けるステップ(但しyは正の整数である)と、
各前記ユニットにおける各前記血管プリセット輪郭線上に位置するz個の点を取得するステップと、
前記血管中心直線に垂直な方向に沿って、z個点をそれぞれ前記血管中心直線に段階的に寄せ、z個の寄せ点を生成するステップ(但しzは正の整数である)と、
RGB差閾値をΔRGB閾として設定し、前記血管中心直線に垂直な方向に沿って、毎回の寄せは前記寄せ点のRGB値と前記血管中心直線上の点のRGB値を比較し、差がΔRGB閾以下になった場合に、前記寄せ点は前記血管中心直線に寄せるのを停止するステップと、
前記寄せ点を輪郭点として取得するステップと、
順に前記輪郭点を接続して形成された平滑な曲線が前記真っすぐに伸ばした後の血管の輪郭線になるステップと、
を含む、
上記の血管中心線に基づき血管輪郭線を取得する方法。
前記の冠状動脈二次元造影画像に基づき血管中心線を抽出するステップは、
冠状動脈二次元造影画像を読み取るステップと、
注目する血管段を取得するステップと、
前記注目する血管段の始点、シードポイント及び終点をピックアップするステップと、
始点、シードポイント、終点の隣接する二点間の二次元造影画像をそれぞれ分割し、少なくとも2つの部分的血管区域図を得るステップと、
各前記部分的血管区域図から少なくとも一本の血管部分経路線を抽出するステップと、
各前記部分的血管区域図上の対応する血管部分経路線を接続し、少なくとも一本の前記血管経路線を取得するステップと、
一本の前記血管経路線を前記血管中心線として選び取るステップと、
を含む、
上記の血管中心線に基づき血管輪郭線を取得する方法。
前記の各前記部分的血管区域図から少なくとも一本の血管部分経路線を抽出するステップは、
前記部分的血管区域図に対して画像増幅処理を行い、コントラストが強い粗血管図を得るステップと、
前記粗血管図に対して格子分割を行い、前記始点から前記終点方向に沿って、少なくとも一本の血管部分経路線を抽出するステップと、
を含む、
上記の血管中心線に基づき血管輪郭線を取得する方法。
任意選択的に、
前記の部分的血管区域図に対して画像増幅処理を行い、コントラストが強い粗血管図を得るステップは、
各前記部分的血管区域図において、前記注目する血管段を前景とし、他の区域を背景とし、前記前景を強くし、前記背景を弱くし、コントラストが強い前記粗血管図を得るステップと、
を含む、
上記の血管中心線に基づき血管輪郭線を取得する方法。
前記の前記粗血管図に対して格子分割を行い、前記始点から前記終点方向に沿って、少なくとも一本の血管部分経路線を抽出するステップは、
前記粗血管図に対して格子分割を行うステップと、
前記始点から前記終点までの血管の延在方向に沿って、前記始点と周辺のn個の格子上の交差点の間の最短時間経路を見つけ出して第二の点とし、前記第二の点と周辺のn個の格子上の交差点の間の最短時間経路を見つけ出して第三の点とし、前記第三の点以降も最短時間経路が終点に達するまで上記ステップを繰り返すステップ(但しnは1以上の正の整数である)と、
見つけ出す順に、前記始点から前記終点までの血管の延在方向に線を繋ぎ、少なくとも一本の血管部分経路線を得るステップと、
を含む、
上記の血管中心線に基づき血管輪郭線を取得する方法。
前記一本の前記血管経路線を前記血管中心線として選び取るステップは、
血管経路線が二本又は二本より多い場合、各血管経路線が前記始点から前記終点までに要する時間に対して和を求めるステップと、
用いる時に最も少ない前記血管経路線を前記血管中心線とするステップと、
を含む、
上記の血管中心線に基づき血管輪郭線を取得する方法。
順に接続された中心線抽出ユニットと、真っすぐに伸ばすユニットと、第一血管輪郭線ユニットと、第二血管輪郭線ユニットと、を備え、
前記中心線抽出ユニットは、冠状動脈二次元造影画像に基づき血管中心線を抽出するためのものであり、
前記真っすぐに伸ばすユニットは、前記の前記中心線抽出ユニットが抽出した血管中心線に基づき、真っすぐに伸ばした血管画像を取得するためのものであり、
前記第一血管輪郭線ユニットは、前記真っすぐに伸ばすユニットが送信した真っすぐに伸ばした血管画像上において、血管直径閾値D閾を設定するため、且つ、前記D閾に基づき、前記血管中心直線の両側に血管プリセット輪郭線を生成するため、且つ、前記血管プリセット輪郭線を前記血管中心直線に段階的に寄せ、真っすぐに伸ばした後の血管の輪郭線を取得するためのものであり、
前記第二血管輪郭線ユニットは、前記第一血管輪郭線ユニットが送信した前記真っすぐに伸ばした後の血管の輪郭線を前記血管中心線的画像上に投射し直し、血管輪郭線を取得するためのものである、
血管中心線に基づき血管輪郭線を取得する装置。
前記真っすぐに伸ばすユニットは、順に接続された中心線を真っすぐに伸ばすモジュールと、区域分割モジュールと、血管画像を真っすぐに伸ばすモジュールと、を備え、
前記中心線を真っすぐに伸ばすモジュールは、前記血管中心線を真っすぐに伸ばし、血管中心直線を取得するためのものであり、
前記区域分割モジュールは、前記始点から前記終点までの血管の延在方向に沿って、部分的血管区域図をx個のユニットに分けるためのものであり(但しxは正の整数である)、
前記血管画像を真っすぐに伸ばすモジュールは、対応して設けた後の画像が前記真っすぐに伸ばした血管画像になるように、各前記ユニットの血管中心線を前記血管中心直線に沿って対応して設けるためのものである、
上記の血管中心線に基づき血管輪郭線を取得する装置。
前記第一血管輪郭線ユニットはさらに、順に接続された輪郭線分割モジュールと、輪郭点抽出モジュールと、血管輪郭線を真っすぐに伸ばすモジュールと、を備え、
前記輪郭線分割モジュールは前記血管画像を真っすぐに伸ばすモジュールに接続され、真っすぐに伸ばした血管画像上において、血管直径閾値D閾を設定し、D閾に基づき、血管中心直線の両側に血管プリセット輪郭線を生成するため、且つ、前記血管プリセット輪郭線をy個のユニットに分けるためのものであり(但しyは正の整数である)、
前記輪郭点抽出モジュールは、前記輪郭線分割モジュールが送信したy個のユニットに基づき、各前記ユニットにおける各前記血管プリセット輪郭線上に位置するz個の点を取得するため、且つ、前記血管中心直線に垂直な方向に沿って、z個点をそれぞれ前記血管中心直線に段階的に寄せ、z個の寄せ点を生成するため(但しzは正の整数である)、且つ、RGB差閾値をΔRGB閾として設定し、前記血管中心直線に垂直な方向に沿って、毎回の寄せは前記寄せ点のRGB値と前記血管中心直線上の点のRGB値を比較し、差がΔRGB閾以下になった場合に、前記寄せ点は前記血管中心直線に寄せるのを停止するため、且つ、前記寄せ点を輪郭点として取得するためのものであり、
前記血管輪郭線を真っすぐに伸ばすモジュールは、形成された平滑な曲線が前記真っすぐに伸ばした後の血管の輪郭線になるように、前記輪郭点抽出モジュールが送信した輪郭点を順に接続するためのものである、
上記の血管中心線に基づき血管輪郭線を取得する装置。
本出願は、前記血管中心線に基づき、真っすぐに伸ばした血管画像を取得するステップと、前記真っすぐに伸ばした血管画像上において、血管直径閾値D閾を設定するステップと、前記D閾に基づき、前記血管中心直線の両側に血管プリセット輪郭線を生成するステップと、前記血管プリセット輪郭線を前記血管中心直線に段階的に寄せ、真っすぐに伸ばした後の血管の輪郭線を取得するステップと、前記真っすぐに伸ばした後の血管の輪郭線を前記血管中心線の画像上に投射し直し、血管輪郭線を取得するステップと、を含む冠状動脈二次元造影画像から血管中心線を抽出する方法を提供することにより、血管輪郭線の抽出が素早く、正確になる。
図1に示されるように、本出願は上記問題を解決するために、以下のステップを含む血管中心線に基づき血管輪郭線を取得する方法を提供する:
当該ステップは図2に示されるように以下のステップを含む:
S110 冠状動脈二次元造影画像を読み取るステップ、
S120 注目する血管段を取得するステップ、
S130 注目する血管段の始点、シードポイント及び終点をピックアップするステップ、
そのうち、シードポイントは演算の精度を高めるためであり、ピックアップが比較的正確な場合においては始点及び終点のみを採用してもよい、
S140 始点、シードポイント、終点の隣接する二点間の二次元造影画像をそれぞれ分割し、少なくとも2つの部分的血管区域図を得るステップ、
好ましくは、部分的血管区域図は矩形図である。各々の部分血管区域のみを計算するため、演算の複雑さ及び冗長程度を低減でき、演算速度を高められ、設計は科学的である、
S150 各部分的血管区域図から少なくとも一本の血管部分経路線を抽出するステップ、
当該ステップは図3に示されるように以下のステップを含む:
S151 部分的血管区域図に対して画像増幅処理を行い、コントラストが強い粗血管図を得るステップ、
該ステップは以下のステップを含む:
各部分的血管区域図において、注目する血管段を前景とし、他の区域を背景とし、前景を強くし、背景を弱くし、コントラストが強い粗血管図を得るステップ、
S152 粗血管図に対して格子分割を行い、始点から終点方向に沿って、少なくとも一本の血管部分経路線を抽出するステップ、
当該ステップは図4に示されるように以下のステップを含む:
S1521 粗血管図に対して格子分割するステップ、
S1522 始点から終点までの血管の延在方向に沿って、始点と周辺のn個の格子上の交差点の間の最短時間経路を見つけ出して第二の点とし、第二の点と周辺のn個の格子上の交差点の間の最短時間経路を見つけ出して第三の点とし、第三の点以降も最短時間経路が終点に達するまで上記ステップを繰り返すステップ(但しnは1以上の正の整数である)、
S1523 見つけ出す順に、始点から終点までの血管の延在方向に線を繋ぎ、少なくとも一本の血管部分経路線を得るステップ、
S160 各部分の血管区域図上の対応する血管部分経路線を接続し、少なくとも一本の血管経路線を取得するステップ、
S170 一本の血管経路線を血管中心線として選び取るステップ、
当該ステップは以下のステップを含む:
血管経路線が二本又は二本より多い場合、各血管経路線が始点から終点までに要する時間に対して和を求めるステップ、及び、用いる時に最も少ない血管経路線を血管中心線とするステップ、
S200 血管中心線に基づき、真っすぐに伸ばした血管画像を取得するステップ、
当該ステップは図5に示されるように以下のステップを含む:
S210 血管中心線真っすぐに伸ばし、血管中心直線を取得するステップ、
S220 始点から終点までの血管の延在方向に沿って、部分的血管区域図をx個のユニットに分けるステップ(但しxは正の整数である)、
S230 各ユニットの血管中心線を血管中心直線に沿って対応して設けるステップ、
S240 対応して設けた後の画像が真っすぐに伸ばした血管画像になるステップ、
S300 真っすぐに伸ばした血管画像上において、血管直径閾値D閾を設定するステップ、
S400 D閾に基づき、血管中心直線の両側に血管プリセット輪郭線を生成するステップ、
S500 血管プリセット輪郭線を血管中心直線に段階的に寄せ、真っすぐに伸ばした後の血管の輪郭線を取得するステップ、
当該ステップは図6に示されるように以下のステップを含む:
S510 血管プリセット輪郭線をy個のユニットに分けるステップ(但しyは正の整数である)、
S520 各ユニットにおける各血管プリセット輪郭線上に位置するz個点を取得するステップ、
S530 血管中心直線に垂直な方向に沿って、z個点をそれぞれ血管中心直線に段階的に寄せ、z個の寄せ点を生成するステップ(但しzは正の整数である)、
S540 RGB差閾値をΔRGB閾として設定し、血管中心直線に垂直な方向に沿って、毎回の寄せは寄せ点のRGB値と血管中心直線上の点のRGB値を比較し、差がΔRGB閾以下になった場合に、寄せ点は血管中心直線に寄せるのを停止するステップ、
S550 寄せ点を輪郭点として取得するステップ、
S560 順に輪郭点を接続して形成された平滑な曲線が真っすぐに伸ばした後の血管の輪郭線になるステップ、及び
S600 真っすぐに伸ばした後の血管の輪郭線を血管中心線の画像上に投射し直し、血管輪郭線を取得するステップ。
図7に示されるように、本出願は以下のユニットを備える血管中心線に基づき血管輪郭線を取得する装置を提供する:
順に接続された中心線抽出ユニット100、真っすぐに伸ばすユニット200、第一血管輪郭線ユニット300、第二血管輪郭線ユニット400、冠状動脈二次元造影画像に基づき血管中心線を抽出するための中心線抽出ユニット100、中心線抽出ユニット100が抽出した血管中心線に基づき、真っすぐに伸ばした血管画像を取得するための真っすぐに伸ばすユニット200、真っすぐに伸ばすユニット200が送信した真っすぐに伸ばした血管画像上において、血管直径閾値D閾を設定するため、且つ、D閾に基づき、血管中心直線の両側に血管プリセット輪郭線を生成するため、且つ、血管プリセット輪郭線を血管中心直線に段階的に寄せ、真っすぐに伸ばした後の血管の輪郭線を取得するための第一血管輪郭線ユニット300、及び第一輪郭線ユニット300が送信した真っすぐに伸ばした後の血管の輪郭線を血管中心線の画像上に投射し直し、血管輪郭線を取得するための第二血管輪郭線ユニット400。
中心線を真っすぐに伸ばすモジュール210は、血管中心線真っすぐに伸ばし、血管中心直線を取得するためのものであり、
区域分割モジュール220は、始点から終点までの血管の延在方向に沿って、部分的血管区域図をx個のユニットに分けるためのものであり(但しxは正の整数である)、
血管画像を真っすぐに伸ばすモジュール230は、対応して設けた後の画像が真っすぐに伸ばした血管画像になるように、各ユニットの血管中心線を血管中心直線に沿って対応して設けるためのものである。
輪郭線分割モジュール310は血管画像を真っすぐに伸ばすモジュール230に接続され、真っすぐに伸ばした血管画像上において、血管直径閾値D閾を設定し、D閾に基づき、血管中心直線の両側に血管プリセット輪郭線を生成し、前記血管プリセット輪郭線をy個のユニットに分けるためのものであり(但しyは正の整数である)、
輪郭点抽出モジュール320は、輪郭線分割モジュールが送信したy個のユニットに基づき、各ユニットにおける各血管プリセット輪郭線上に位置するz個点を取得するため、且つ、血管中心直線に垂直な方向に沿って、z個点をそれぞれ血管中心直線に段階的に寄せ、z個の寄せ点を生成するため(但しzは正の整数である)、かつ、RGB差閾値をΔRGB閾として設定し、血管中心直線に垂直な方向に沿って、毎回の寄せは寄せ点のRGB値と血管中心直線上の点のRGB値を比較し、差がΔRGB閾以下になった場合に、寄せ点は血管中心直線に寄せるのを停止するため、且つ、寄せ点を輪郭点として取得するためのものであり、
血管輪郭線を真っすぐに伸ばすモジュール330は、形成された平滑な曲線が真っすぐに伸ばした後の血管の輪郭線になるように、輪郭点抽出モジュールが送信した輪郭点順に接続するためのものである。
真っすぐに伸ばすユニット200
中心線を真っすぐに伸ばすモジュール210
区域分割モジュール220
血管画像を真っすぐに伸ばすモジュール230
第一血管輪郭線ユニット300
輪郭線分割モジュール310
輪郭点抽出モジュール320
血管輪郭線を真っすぐに伸ばすモジュール330
第二血管輪郭線ユニット400
Claims (13)
- 冠状動脈二次元造影画像に基づき血管中心線を抽出するステップと、
前記血管中心線に基づき、真っすぐに伸ばした血管画像を取得するステップと、
前記真っすぐに伸ばした血管画像上において、血管直径閾値D閾を設定するステップと、
前記D閾に基づき、血管中心直線の両側に血管プリセット輪郭線を生成するステップと、
前記血管プリセット輪郭線を前記血管中心直線に段階的に寄せ、真っすぐに伸ばした後の血管の輪郭線を取得するステップと、
前記真っすぐに伸ばした後の血管の輪郭線を前記血管中心線の画像上に投射し直し、血管輪郭線を取得するステップと、
を含む、
ことを特徴とする血管中心線に基づき血管輪郭線を取得する方法。 - 前記血管中心線に基づき、真っすぐに伸ばした血管画像を取得するステップは、
前記血管中心線を真っすぐに伸ばし、血管中心直線を取得するステップと、
始点から終点までの血管の延在方向に沿って、部分的血管区域図をx個のユニットに分けるステップ(但しxは正の整数である)と、
各前記ユニットの血管中心線を前記血管中心直線に沿って対応して設けるステップと、
対応して設けた後の画像が前記真っすぐに伸ばした血管画像になるステップと、
を含む、
ことを特徴とする請求項1に記載の血管中心線に基づき血管輪郭線を取得する方法。 - 前記血管プリセット輪郭線を前記血管中心直線に段階的に寄せ、真っすぐに伸ばした後の血管の輪郭線を取得するステップは、
前記血管プリセット輪郭線をy個のユニットに分けるステップ(但しyは正の整数である)と、
各前記ユニットにおける各前記血管プリセット輪郭線上に位置するz個の点を取得するステップと、
前記血管中心直線に垂直な方向に沿って、z個点をそれぞれ前記血管中心直線に段階的に寄せ、z個の寄せ点を生成するステップ(但しzは正の整数である)と、
RGB差閾値をΔRGB閾として設定し、前記血管中心直線に垂直な方向に沿って、毎回の寄せは前記寄せ点のRGB値と前記血管中心直線上の点のRGB値を比較し、差がΔRGB閾以下になった場合に、前記寄せ点は前記血管中心直線に寄せるのを停止するステップと、
前記寄せ点を輪郭点として取得するステップと、
順に前記輪郭点を接続して形成された平滑な曲線が前記真っすぐに伸ばした後の血管の輪郭線になるステップと、
を含む、
ことを特徴とする請求項2に記載の血管中心線に基づき血管輪郭線を取得する方法。 - 前記の冠状動脈二次元造影画像に基づき血管中心線を抽出するステップは、
冠状動脈二次元造影画像を読み取るステップと、
注目する血管段を取得するステップと、
前記注目する血管段の始点、シードポイント及び終点をピックアップするステップと、
始点、シードポイント、終点の隣接する二点間の二次元造影画像をそれぞれ分割し、少なくとも2つの部分的血管区域図を得るステップと、
各前記部分的血管区域図から少なくとも一本の血管部分経路線を抽出するステップと、
各前記部分的血管区域図上の対応する血管部分経路線を接続し、少なくとも一本の前記血管経路線を取得するステップと、
一本の前記血管経路線を前記血管中心線として選び取るステップと、
を含む、
ことを特徴とする請求項2に記載の血管中心線に基づき血管輪郭線を取得する方法。 - 前記の各前記部分的血管区域図から少なくとも一本の血管部分経路線を抽出するステップは、
前記部分的血管区域図に対して画像増幅処理を行い、コントラストが強い粗血管図を得るステップと、
前記粗血管図に対して格子分割を行い、前記始点から前記終点方向に沿って、少なくとも一本の血管部分経路線を抽出するステップと、
を含む、
ことを特徴とする請求項4に記載の血管中心線に基づき血管輪郭線を取得する方法。 - 前記部分的血管区域図に対して画像増幅処理を行い、コントラストが強い粗血管図を得るステップは、
各前記部分的血管区域図において、前記注目する血管段を前景とし、他の区域を背景とし、前記前景を強くし、前記背景を弱くし、コントラストが強い前記粗血管図を得るステップを含む、
ことを特徴とする請求項5に記載の血管中心線に基づき血管輪郭線を取得する方法。 - 前記粗血管図に対して格子分割を行い、前記始点から前記終点方向に沿って、少なくとも一本の血管部分経路線を抽出するステップは、
前記粗血管図に対して格子分割を行うステップと、
前記始点から前記終点までの血管の延在方向に沿って、前記始点と周辺のn個の格子上の交差点の間の最短時間経路を見つけ出して第二の点とし、前記第二の点と周辺のn個の格子上の交差点の間の最短時間経路を見つけ出して第三の点とし、前記第三の点以降も最短時間経路が終点に達するまで上記ステップを繰り返すステップ(但しnは1以上の正の整数である)と、
見つけ出す順に、前記始点から前記終点までの血管の延在方向に線を繋ぎ、少なくとも一本の血管部分経路線を得るステップと、
を含む、
ことを特徴とする請求項5に記載の血管中心線に基づき血管輪郭線を取得する方法。 - 一本の前記血管経路線を前記血管中心線として選び取るステップは、
血管経路線が二本又は二本より多い場合、各血管経路線が前記始点から前記終点までに要する時間に対して和を求めるステップと、
用いる時に最も少ない前記血管経路線を前記血管中心線とするステップと、
を含む、
ことを特徴とする請求項7に記載の血管中心線に基づき血管輪郭線を取得する方法。 - 順に接続された中心線抽出ユニットと、真っすぐに伸ばすユニットと、第一血管輪郭線ユニットと、第二血管輪郭線ユニットと、を備え、
前記中心線抽出ユニットは、冠状動脈二次元造影画像に基づき血管中心線を抽出するためのものであり、
前記真っすぐに伸ばすユニットは、前記中心線抽出ユニットが抽出した血管中心線に基づき、真っすぐに伸ばした血管画像を取得するためのものであり、
前記第一血管輪郭線ユニットは、前記真っすぐに伸ばすユニットが送信した真っすぐに伸ばした血管画像上において、血管直径閾値D閾を設定するため、且つ、前記D閾に基づき、前記血管中心直線の両側に血管プリセット輪郭線を生成するため、且つ、前記血管プリセット輪郭線を前記血管中心直線に段階的に寄せ、真っすぐに伸ばした後の血管の輪郭線を取得するためのものであり、
前記第二血管輪郭線ユニットは、前記第一血管輪郭線ユニットが送信した前記真っすぐに伸ばした後の血管の輪郭線を前記血管中心線の画像上に投射し直し、血管輪郭線を取得するためのものである、
ことを特徴とする請求項1~8のいずれかに記載の血管中心線に基づき血管輪郭線を取得する方法に用いられる血管中心線に基づき血管輪郭線を取得する装置。 - 前記真っすぐに伸ばすユニットは、順に接続された中心線を真っすぐに伸ばすモジュールと、区域分割モジュールと、血管画像を真っすぐに伸ばすモジュールと、を備え、
前記中心線を真っすぐに伸ばすモジュールは、前記血管中心線を真っすぐに伸ばし、血管中心直線を取得するためのものであり、
前記区域分割モジュールは、前記始点から前記終点までの血管の延在方向に沿って、部分的血管区域図をx個のユニットに分けるためのものであり(但しxは正の整数である)、
前記血管画像を真っすぐに伸ばすモジュールは、対応して設けた後の画像が前記真っすぐに伸ばした血管画像になるように、各前記ユニットの血管中心線を前記血管中心直線に沿って対応して設けるためものである、
ことを特徴とする請求項9に記載の血管中心線に基づき血管輪郭線を取得する装置。 - 前記第一血管輪郭線ユニットはさらに、順に接続された輪郭線分割モジュールと、輪郭点抽出モジュールと、血管輪郭線を真っすぐに伸ばすモジュールと、を備え、
前記輪郭線分割モジュールは前記血管画像を真っすぐに伸ばすモジュールに接続され、真っすぐに伸ばした血管画像上において、血管直径閾値D閾を設定し、D閾に基づき、血管中心直線の両側に血管プリセット輪郭線を生成するため、且つ、前記血管プリセット輪郭線をy個のユニットに分けるためのものであり(但しyは正の整数である)、
前記輪郭点抽出モジュールは、前記輪郭線分割モジュールが送信したy個のユニットに基づき、各前記ユニットにおける各前記血管プリセット輪郭線上に位置するz個の点を取得するため、且つ、前記血管中心直線に垂直な方向に沿って、z個点をそれぞれ前記血管中心直線に段階的に寄せ、z個の寄せ点を生成するため(但しzは正の整数である)、且つ、RGB差閾値をΔRGB閾として設定し、前記血管中心直線に垂直な方向に沿って、毎回の寄せは前記寄せ点のRGB値と前記血管中心直線上の点のRGB値を比較し、差がΔRGB閾以下になった場合に、前記寄せ点は前記血管中心直線に寄せるのを停止するため、且つ、前記寄せ点を輪郭点として取得するためのものであり、
前記血管輪郭線を真っすぐに伸ばすモジュールは、形成された平滑な曲線が前記真っすぐに伸ばした後の血管の輪郭線になるように、前記輪郭点抽出モジュールが送信した輪郭点を順に接続するためのものである、
ことを特徴とする請求項9に記載の血管中心線に基づき血管輪郭線を取得する装置。 - 請求項9~11のいずれかに記載の血管中心線に基づき血管輪郭線を取得する装置を含むことを特徴とする冠状動脈分析システム。
- コンピュータプログラムがプロセッサにより実行される時に請求項1~8のいずれかに記載の血管中心線に基づき血管輪郭線を取得する方法を実現することを特徴とするコンピュータ記憶媒体。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911135307.1A CN112116615B (zh) | 2019-11-19 | 2019-11-19 | 根据血管中心线获取血管轮廓线的方法和装置 |
CN201911135307.1 | 2019-11-19 | ||
PCT/CN2019/120333 WO2021097817A1 (zh) | 2019-11-19 | 2019-11-22 | 根据血管中心线获取血管轮廓线的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2023502409A true JP2023502409A (ja) | 2023-01-24 |
JP7409718B2 JP7409718B2 (ja) | 2024-01-09 |
Family
ID=73795593
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2022529118A Active JP7409718B2 (ja) | 2019-11-19 | 2019-11-22 | 血管中心線に基づき血管輪郭線を取得するステップと方法及び装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20220277447A1 (ja) |
EP (1) | EP4064181A4 (ja) |
JP (1) | JP7409718B2 (ja) |
CN (1) | CN112116615B (ja) |
WO (1) | WO2021097817A1 (ja) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10210956B2 (en) | 2012-10-24 | 2019-02-19 | Cathworks Ltd. | Diagnostically useful results in real time |
WO2014064702A2 (en) | 2012-10-24 | 2014-05-01 | Cathworks Ltd. | Automated measurement system and method for coronary artery disease scoring |
EP3061015A2 (en) | 2013-10-24 | 2016-08-31 | Cathworks Ltd. | Vascular characteristic determination with correspondence modeling of a vascular tree |
EP3461253B1 (en) | 2016-05-16 | 2023-08-09 | Cathworks Ltd. | Selection of vascular paths from images |
WO2017199245A1 (en) | 2016-05-16 | 2017-11-23 | Cathworks Ltd. | System for vascular assessment |
CN117058328B (zh) * | 2023-10-11 | 2024-01-09 | 杭州脉流科技有限公司 | 冠状动脉血管树分级方法、设备、存储介质和程序产品 |
CN117115150B (zh) * | 2023-10-20 | 2024-01-26 | 柏意慧心(杭州)网络科技有限公司 | 用于确定分支血管的方法、计算设备和介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150235360A1 (en) * | 2014-02-18 | 2015-08-20 | Siemens Aktiengesellschaft | Sparse Appearance Learning-based Segmentation |
CN109461138A (zh) * | 2018-09-29 | 2019-03-12 | 深圳睿心智能医疗科技有限公司 | 参数计算方法、系统、可读存储介质及计算机设备 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6047080A (en) * | 1996-06-19 | 2000-04-04 | Arch Development Corporation | Method and apparatus for three-dimensional reconstruction of coronary vessels from angiographic images |
JP4421203B2 (ja) * | 2003-03-20 | 2010-02-24 | 株式会社東芝 | 管腔状構造体の解析処理装置 |
CN102346803A (zh) * | 2010-07-28 | 2012-02-08 | 北京集翔多维信息技术有限公司 | 心血管造影图像分析系统 |
CN106127819B (zh) * | 2016-06-30 | 2019-10-08 | 上海联影医疗科技有限公司 | 医学图像中提取血管中心线的方法及其装置 |
CN106327487B (zh) * | 2016-08-18 | 2018-01-02 | 苏州润迈德医疗科技有限公司 | 基于x射线冠脉造影图像的冠状动脉血流储备分数计算方法 |
US10452813B2 (en) * | 2016-11-17 | 2019-10-22 | Terarecon, Inc. | Medical image identification and interpretation |
EP3503026A1 (en) * | 2017-12-20 | 2019-06-26 | Koninklijke Philips N.V. | Device, system and method for interacting with vessel images |
CN108335284B (zh) * | 2018-01-09 | 2022-06-28 | 北京理工大学 | 一种冠脉血管中心线匹配方法和系统 |
CN108564574B (zh) * | 2018-04-11 | 2021-04-20 | 上海联影医疗科技股份有限公司 | 确定血流储备分数的方法、计算机设备及计算机可读存储介质 |
CN108805815B (zh) * | 2018-06-14 | 2023-02-17 | 苏州润迈德医疗科技有限公司 | 基于x射线血管造影图像的血管拉直重建方法 |
CN109559326B (zh) * | 2018-11-05 | 2020-11-13 | 深圳睿心智能医疗科技有限公司 | 一种血流动力学参数计算方法、系统及电子设备 |
-
2019
- 2019-11-19 CN CN201911135307.1A patent/CN112116615B/zh active Active
- 2019-11-22 WO PCT/CN2019/120333 patent/WO2021097817A1/zh unknown
- 2019-11-22 EP EP19953666.5A patent/EP4064181A4/en active Pending
- 2019-11-22 JP JP2022529118A patent/JP7409718B2/ja active Active
-
2022
- 2022-05-17 US US17/746,087 patent/US20220277447A1/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150235360A1 (en) * | 2014-02-18 | 2015-08-20 | Siemens Aktiengesellschaft | Sparse Appearance Learning-based Segmentation |
CN109461138A (zh) * | 2018-09-29 | 2019-03-12 | 深圳睿心智能医疗科技有限公司 | 参数计算方法、系统、可读存储介质及计算机设备 |
Also Published As
Publication number | Publication date |
---|---|
CN112116615A (zh) | 2020-12-22 |
EP4064181A1 (en) | 2022-09-28 |
CN112116615B (zh) | 2023-12-05 |
EP4064181A4 (en) | 2023-11-29 |
US20220277447A1 (en) | 2022-09-01 |
WO2021097817A1 (zh) | 2021-05-27 |
JP7409718B2 (ja) | 2024-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7409718B2 (ja) | 血管中心線に基づき血管輪郭線を取得するステップと方法及び装置 | |
CN110889896B (zh) | 获取血管狭窄病变区间及三维合成方法、装置和系统 | |
JP2022169579A (ja) | リアルタイムで診断上有用な結果 | |
US11191490B2 (en) | Personalized assessment of patients with acute coronary syndrome | |
WO2021196536A1 (zh) | 精确提取血管中心线的方法、装置、分析系统和存储介质 | |
US11357573B2 (en) | Optimum treatment planning during coronary intervention by simultaneous simulation of a continuum of outcomes | |
CN112132949B (zh) | 具有狭窄病变区间的血管三维建模方法、装置和系统 | |
JP7162934B2 (ja) | 造影画像に基づき血管評定パラメータを取得する方法、装置及びシステム | |
CN112419484B (zh) | 三维血管合成方法、系统及冠状动脉分析系统和存储介质 | |
CN112419276B (zh) | 调节血管轮廓及中心线的方法及存储介质 | |
WO2021097821A1 (zh) | 从冠状动脉二维造影图像中提取血管中心线的方法和装置 | |
US20220261997A1 (en) | Method and apparatus for acquiring blood vessel evaluation parameter based on physiological parameter, and storage medium | |
JP2022554298A (ja) | 微小血管抵抗指数に基づき最大充血状態流速を調節する方法及び装置 | |
WO2021087961A1 (zh) | 测量舒张期血流速度的方法、装置、系统及存储介质 | |
US20220319116A1 (en) | Method and device for synthesizing mathematical model of blood vessel having stenotic lesion | |
CN112446867A (zh) | 血流参数的确定方法、装置、设备及存储介质 | |
WO2022109904A1 (zh) | 二维图像选取及三维血管合成的方法和存储介质 | |
JP7446645B2 (ja) | 深層学習に基づいて大動脈画像を取得するシステム | |
CN112419462B (zh) | 三维血管的渲染合成方法、系统及存储介质 | |
WO2021109116A1 (zh) | 用于流体力学分析的圆台血管数学模型的合成方法和装置 | |
CN115222665B (zh) | 一种斑块检测方法、装置、电子设备及可读存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220518 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20230509 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230802 |
|
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: 20231121 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20231214 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7409718 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |