JP2013536706A - 再帰的分離管状器官のツリーの定量的分析方法および装置 - Google Patents
再帰的分離管状器官のツリーの定量的分析方法および装置 Download PDFInfo
- Publication number
- JP2013536706A JP2013536706A JP2013526324A JP2013526324A JP2013536706A JP 2013536706 A JP2013536706 A JP 2013536706A JP 2013526324 A JP2013526324 A JP 2013526324A JP 2013526324 A JP2013526324 A JP 2013526324A JP 2013536706 A JP2013536706 A JP 2013536706A
- Authority
- JP
- Japan
- Prior art keywords
- tree
- centerline
- branch
- cross
- contour
- 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
- 210000000056 organ Anatomy 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 62
- 238000004445 quantitative analysis Methods 0.000 title claims abstract description 18
- 238000000926 separation method Methods 0.000 claims abstract description 74
- 230000033458 reproduction Effects 0.000 claims description 26
- 210000004204 blood vessel Anatomy 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 15
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 claims description 9
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 claims description 9
- 230000005855 radiation Effects 0.000 claims description 8
- 239000012530 fluid Substances 0.000 claims description 6
- 230000003902 lesion Effects 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 5
- 230000000007 visual effect Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 3
- 108020001568 subdomains Proteins 0.000 claims description 2
- 230000002093 peripheral effect Effects 0.000 claims 1
- 238000004590 computer program Methods 0.000 abstract 1
- 238000002955 isolation Methods 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 30
- 208000031481 Pathologic Constriction Diseases 0.000 description 28
- 230000036262 stenosis Effects 0.000 description 28
- 208000037804 stenosis Diseases 0.000 description 28
- 230000002792 vascular Effects 0.000 description 14
- 238000005259 measurement Methods 0.000 description 9
- 210000001367 artery Anatomy 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 6
- 238000013459 approach Methods 0.000 description 5
- 210000000709 aorta Anatomy 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 101100520142 Caenorhabditis elegans pin-2 gene Proteins 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000002059 diagnostic imaging Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 101150037009 pin1 gene Proteins 0.000 description 3
- 238000011002 quantification Methods 0.000 description 3
- 230000002966 stenotic effect Effects 0.000 description 3
- 238000002604 ultrasonography Methods 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000010339 dilation Effects 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000002583 angiography Methods 0.000 description 1
- 230000003143 atherosclerotic effect Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 210000001715 carotid artery Anatomy 0.000 description 1
- 238000005119 centrifugation Methods 0.000 description 1
- 210000001072 colon Anatomy 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 210000004351 coronary vessel Anatomy 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 230000008338 local blood flow Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B15/00—ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- 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/149—Segmentation; Edge detection involving deformable models, e.g. active contour models
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/04—Indexing scheme for image data processing or generation, in general involving 3D image data
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Medical Informatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Biophysics (AREA)
- Biotechnology (AREA)
- Evolutionary Biology (AREA)
- Software Systems (AREA)
- Chemical & Material Sciences (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Crystallography & Structural Chemistry (AREA)
- Physiology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Molecular Biology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
Abstract
【選択図】図2
Description
‐当該ツリーまたは当該ツリーの一部の3Dモデルを提供するステップであって、この3Dモデルが、前記ツリーまたは前記ツリーの一部を形成する前記管状器官のルーメン壁の表面を表示する、ステップと、
‐当該ツリーまたは前記ツリーの一部の3D中心線を定義するステップと、
‐前記ツリーの分岐を識別するステップと、
‐前記ツリーのまたは前記ツリーの一部のN-分岐を識別するステップであって、N-分岐が、近位の管状器官がN(N≧2)個の遠位管状器官に分岐している前記ツリーの一部である、ステップとを備える方法において、
‐更に、前記3Dモデルを得るために使用される前記モダリティーから独立して、各分岐を1個以上の領域に分割するステップであって、前記領域が、単一血管領域および分離領域と呼ばれる二種類のタイプが異なる領域であり、このような領域において異なる横断面表面が定義されている、ステップを備え、
‐当該分離領域が、分岐の前記近位の側と当該分岐の前記遠位側に存在することが可能であり、かつ、
‐各N-分岐が、分岐の前記遠位分離領域および当該分岐に直接遠位である前記N分岐の前記近位分離領域を備える。
‐前記方位を球状ドメインで定義し、これにより前記半径の次元を無視する、ステップ;
‐当該ドメインを両方の残りの角度の次元の各々において一サイクルに限定するステップ;
‐当該ドメインを矩形のサブドメイン内に分離し、かつ各サブドメイン毎に、当該サブドメインの前記中心値によって指定された前記方位に対する前記断面エリアを算出するステップ;
‐当該ドメインを、当該サブドメイン(より詳しくは、当該中心値に対する当該断面積の前記算出が、結果として、最小エリアとなるサブドメイン)の1つによって置換するステップ;
‐前記ドメイン・サイズが既定の閾値に達するまで、当該中心値の算出および当該サブドメインの置換を繰り返すステップ。
‐前記ルーメンの周囲の輪郭線を描く、対応する等値輪郭線を決定するために、その近位側で前記分岐を囲む前記周囲の輪郭線と前記遠位境界面を囲む前記輪郭線との間の前記レベル・セットを算出するステップと、
‐前記距離に位置する前記中心線ポイントに対応する前記等値に基づいて、前記等値輪郭線を前記中心線ポイントに連結させるステップであって、最初と最後の輪郭線の間の前記合計中心線距離の比率が、当該最初と最後の等値輪郭線の前記レベルの間で同じ比率を有する前記等値に等しい、ステップと、
‐各中心線ポイントに対し、当該中心線ポイントの前記連結された等値輪郭線によって決められる前記横断面の前記表面を算出するステップであって、当該表面が、湾曲していてかつ前記最小エネルギーを有する、ステップと、
‐各当該中心線に対し、当該表面の前記エリアを算出するステップと
によって、算出することができる。
‐前記近位境界面を囲む前記輪郭線とその遠位側で前記分岐を決める前記周囲輪郭線との間に前記レベル・セットを算出するステップと、
‐前記距離に位置する前記中心線ポイントに対応する前記等値に基づいて、当該輪郭線の間に最も長いルートを有する前記中心線の前記中心線ポイントに前記等値輪郭線を連結させるステップであって、前記最初と最後の等値輪郭線の間の前記合計中心線距離の前記比率が、当該最初と最後の等値輪郭線の前記レベルの間で同じ比率を有している前記等値に等しい、ステップと、
‐前記分岐の終端まで中心線の前記分離との間で測定される前記合計距離に比較される中心線ポイントの当該分離点への前記距離の前記比率に基づいて、すでに存在する連結を使用して、中心線の分離点の後の、前記残りの中心線の中心線ポイントに、前記等値輪郭線を連結するステップであって、当該すでに存在している連結が、前記分岐の当該分離と前記端との間の距離の等しい比率を有する、ステップと、
‐当該中心線ポイントの前記連結された等値輪郭線により決められる前記横断面の各中心線ポイントに対し、前記表面を算出するステップであって、当該表面が、湾曲していてかつ前記最小エネルギーを有している、ステップと、
‐各当該中心線ポイントに対して当該表面の前記エリアを算出するステップ
によって算出することができる。
‐前記近位境界面を囲む前記輪郭線と前記遠位境界面を囲む前記輪郭線との間に前記レベル・セットを算出するステップと、
‐各中心線ポイントをレベル・セット値の最小標準偏差を有する断面の輪郭線に連結するステップであって、これによって当該輪郭線が、当該中心線ポイントを含む横断面を囲む、ステップと、
‐各当該中心線ポイントに対し当該横断面のエリアを算出するステップ
によって算出することができる。
‐当該周囲輪郭線の各輪郭線頂点から当該周囲等輪郭線の前記重心への直線を定義し、かつ当該線を、mサンプルを含むnリングに分割するステップであって、当該周囲輪郭線の前記頂点が、前記第一のリングの前記サンプルであり、かつ各内側リングはその隣の外側リングより少ないサンプルを含み、かつ前記最初の2つ以外には、2個の連続したリングが同じ放射線をサンプリングしない、ステップ;
‐リングのポイントの間で表面メッシュを定義するステップ;
‐当該表面メッシュの合計エリアを算出するステップ;
‐当該第一リングの当該サンプルを固定し、かつメッシュ・エリアの変更および/またはサンプルの動きが、既定のしきい値以下になるまで、平均隣接方向にサンプルを互いに移動し続けるステップ。
‐輪郭線モデルとして、前記単一管状器官の健康な再現を作成するステップと、
‐前記N-分岐の健康な再現を作成するステップであって、当該健康な再現の各々が、N個の部分に分割され、当該部分の各々が、近位から遠位に行く最小ルーメン曲率を記述しかつ輪郭線モデルとして当該部分を結合する、ステップと、
‐全ての再現された単一管状器官およびN-分岐を1個の輪郭線モデルに結合するステップと、
‐前記ツリーまたは前記ツリーの一部を形成する前記健康な管状器官を形成する当該輪郭線モデルから健康なモデル表面メッシュを作成するステップと、
‐前記ツリーまたは前記ツリーの一部を形成する前記健康な管状器官に対し3D中心線を定義するステップと、
‐前記ツリーまたは前記ツリーの一部を形成する前記健康な管状器官のN-分岐を定義するステップと、
‐前記健康なモデル・ツリーの連結する中心線を前記ツリーの前記中心線に連結させるステップと、
‐断面を作成しかつN-分岐のエリア曲線を算出するステップと、
‐単一管状器官に対し断面を作成するステップと、
‐算出された参照に基づいて単一管状器官に対しエリア曲線を算出するステップと
を備える。
「近位の」は、血管または管または管状器官管の、流体が流れ込む側を意味する用語である。例えば、大動脈の近位側とは、大動脈の心臓に最も近い側のことである。
「遠位の」は、血管または管または管状器官の、流体が流れ出す側を意味する用語である。例えば、大動脈の遠位側とは、大動脈の心臓から最も離れた側のことである。
「ルーメン」は、流体が流れる血管または管または管状器官の内側または血管壁の外部境界、またはその間の中間位置を意味する用語である。
「健康である」は、狭窄または拡張が無いとみなせる、管状器官ディメンジョンを意味する用語である。
先行技術は、図1のブロック20により表される血管の3Dルーメン壁の定義を教示する。これは、流入面と流出面で開放されているルーメン壁表面メッシュであると仮定されている。X線の場合、例えば、3Dツリーをモデル化するためには、特許文献4に開示されているように、別々に再構成された二分岐を接続することにより、処理を行うことが出来る。MRAの別の具体例では、これは、非特許文献1のように実行することも出来る。超音波における別の具体例では、これは、特許文献5のように実行することもできる。CTAにおけるさらにもう一個の具体例では、これは、非特許文献2のようにも、達成することができる。
1) 3D輪郭線モデル(すなわち、(健康な)ルーメンを記述する中心線位置ごとの中心線および断面の輪郭線)が、各単一管状器官ごとに形成される(図1のブロック210を参照)。我々は、ツリーに対する3D狭窄境界と直径参照曲線を作成し(ユーザが定義しない場合、両方とも自動的に算出することができる)、かつこれから各単一管状器官の3Dの健康なモデルを作成する。これは、先行技術でありかつ例えば、特許文献4に記載される方法のサブセットにより達成することができる。
2b. 中心ポイントpcenterは、pin1とpout1を通る線の交叉ポイント並びにpin2とpout2を通る線の交叉ポイントとして見出される(図1のブロック240を参照)。
2c. 3D中心線は、pproxとpdistの間に生成される。角A(Aは、ベクトル長rproxとrdistを有するベクトルpprox-pcenterとpdist-pcenterの間の全ての現実の角に制限されている)の関数として、我々は、pcenterから距離rA(ここで、rAは、角Aに基づいてrproxとrdistの間に線形に内挿されている)で、中心線位置cAを定義する(図10、図1のブロック250を参照)。
2d. 各中心線位置cAに対し、断面の輪郭線は、中心cA、Aに従って直径pproxとpdistの間に線形に内挿された直径、および局所的な中心線方向に等しい法線を有する円である (図1のブロック260を参照)。
2c、2dを結合することは、N-分岐の部分iに対するステップ1におけるような輪郭線モデルを与える。
N-分岐の全ての部分の輪郭線モデルは、N-分岐の輪郭線モデルを形成するために結合される(図1のブロック270を参照)。
3Dの健康なモデルが定義されると、更なるステップを、元の3Dモデルに対するステップと同じものとすることが出来る。これは、図1において、ブロック290と30の連結によって示されている。もちろんまた、異なるステップを、実行することも可能である。
Claims (24)
- 再帰的分離管状器官のツリーのまたはツリーの一部の定量分析のための方法であって、
‐当該ツリーまたは当該ツリーの一部の3Dモデルを提供するステップであって、この3Dモデルが、前記ツリーまたは前記ツリーの一部を形成する前記管状器官のルーメン壁の表面を表示する、ステップと、
‐当該ツリーまたは前記ツリーの一部の3D中心線を定義するステップと、
‐前記ツリーの分岐を識別するステップと、
‐前記ツリーのまたは前記ツリーの一部のN-分岐を識別するステップであって、N-分岐が、近位の管状器官がN(N≧2)個の遠位管状器官に分岐している前記ツリーの一部である、ステップとを備える方法において、
‐更に、前記3Dモデルを得るために使用される前記モダリティーから独立して、各分岐を1個以上の領域に分割するステップであって、前記領域が、単一血管領域および分離領域と呼ばれる二種類のタイプが異なる領域であり、このような領域において異なる横断面表面が定義されている、ステップと、
を備え、
‐当該分離領域が、分岐の前記近位の側と当該分岐の前記遠位側に存在することが可能であり、かつ、
‐各N-分岐が、分岐の前記遠位分離領域および当該分岐に直接遠位である前記N分岐の前記近位分離領域を備える
ことを特徴とする方法。 - 分離領域と単一血管領域との間の境界面(EON, SON)が、各分岐に対して定義されていて、前記境界面が前記単一血管領域タイプの横断面である、請求項1に記載の方法。
- 前記分離領域の前記横断面が、湾曲表面であり、かつ前記単一血管領域の前記横断面が、平面である、請求項1または2に記載の方法。
- 前記分離領域の前記横断面が、前記分離方向のそれぞれに前記ツリーまたは前記ツリーの一部をかん流する流体の流れに対する垂直性を考慮するために、湾曲している、請求項3に記載の方法。
- 前記平面が、切断面(特に、前記局所分岐方位に垂直な切断面)に属している、請求項3または4に記載の方法。
- 前記中心線各位置での前記横断面エリアが、単一血管領域と分離領域に対して定義されている前記横断面表面のエリアを算出することにより決定される、請求項1から5の何れか1項に記載の方法。
- 分岐の単一血管領域の近位境界面(EON)が、当該分岐の前記中心線に渡って最も近位の平面横断面を連続して見出すことにより算出され、その輪郭線が完全に当該分岐の部分であり、当該中心線に前記対応する位置を含む単一面内にあり、かつそれが当該位置で最小エリアを有するように配置されていて、
分岐の単一血管領域の遠位境界面(SON)が、当該分岐の前記中心線に渡って最も近位の平面横断面を連続して見出すことにより算出され、その輪郭線が完全に当該分岐の部分であり、当該中心線に前記対応する位置を含む単一面内にあり、かつそれが当該位置で最小エリアを有するように配置されている、
請求項2から6の何れか1項に記載の方法。 - ‐前記方位を球状ドメインで定義し、これにより前記半径の次元を無視する、ステップと、
‐当該ドメインを両方の残りの角度の次元の各々において一サイクルに限定するステップと、
‐当該ドメインを矩形のサブドメイン内に分離し、かつ各サブドメイン毎に、当該サブドメインの前記中心値によって指定された前記方位に対する前記断面エリアを算出するステップと、
‐当該ドメインを、当該サブドメイン(より詳しくは、当該中心値に対する当該断面積の前記算出が、結果として、最小エリアとなるサブドメイン)の1個によって置換するステップと、
‐前記ドメイン・サイズが既定の閾値に達するまで、当該中心値の算出および当該サブドメインの置換を繰り返すステップと
により、当該横断面のエリアが最小となる当該横断面の方位を見出すことが、算出される、請求項7に記載の方法。 - 近位分離領域に対し、湾曲する表面横断面のエリアが、
‐前記ルーメンの周囲の輪郭線を描く、対応する等値輪郭線を決定するために、その近位側で前記分岐を囲む前記周囲の輪郭線と前記遠位境界面(EON)を囲む前記輪郭線との間の前記レベル・セットを算出するステップと、
‐前記距離に位置する前記中心線ポイントに対応する前記等値に基づいて、前記等値輪郭線を前記中心線ポイントに連結させるステップであって、最初と最後の輪郭線の間の前記合計中心線距離の比率が、当該最初と最後の等値輪郭線の前記レベルの間で同じ比率を有する前記等値に等しい、ステップと、
‐各中心線ポイントに対し、当該中心線ポイントの前記連結された等値輪郭線によって決められる前記横断面の前記表面を算出するステップであって、当該表面が、湾曲していてかつ前記最小エネルギーを有する、ステップと、
‐各当該中心線に対し、当該表面の前記エリアを算出するステップと
によって算出される、請求項1から8の何れか1項に記載の方法。 - 遠位分離領域に対して、前記湾曲した表面の横断面のエリアが、
‐前記近位境界面を囲む前記輪郭線とその遠位側で前記分岐を決める前記周囲輪郭線との間に前記レベル・セットを算出するステップと、
‐前記距離に位置する前記中心線ポイントに対応する前記等値に基づいて、当該輪郭線の間に最も長いルートを有する前記中心線の前記中心線ポイントに前記等値輪郭線を連結させるステップであって、前記最初と最後の等値輪郭線の間の前記合計中心線距離の前記比率が、当該最初と最後の等値輪郭線の前記レベルの間で同じ比率を有している前記等値に等しい、ステップと、
‐前記分岐の終端まで中心線の前記分離との間で測定される前記合計距離に比較される中心線ポイントの当該分離点への前記距離の前記比率に基づいて、すでに存在する連結を使用して、中心線の分離点の後の、前記残りの中心線の中心線ポイントに、前記等値輪郭線を連結するステップであって、当該すでに存在している連結が、前記分岐の当該分離と前記端との間の距離の等しい比率を有する、ステップと、
‐当該中心線ポイントの前記連結された等値輪郭線により決められる前記横断面の各中心線ポイントに対し、前記表面を算出するステップであって、当該表面が、湾曲していてかつ前記最小エネルギーを有している、ステップと、
‐各当該中心線ポイントに対して当該表面の前記エリアを算出するステップ
によって算出される、請求項1から9の何れか1項に記載の方法。 - 単一血管領域に対して、前記平面の横断面のエリアが、
‐前記近位境界面を囲む前記輪郭線と前記遠位境界面を囲む前記輪郭線との間に前記レベル・セットを算出するステップと、
‐各中心線ポイントをレベル・セット値の最小標準偏差を有する断面の輪郭線に連結するステップであって、これによって当該輪郭線が、当該中心線ポイントを含む横断面を囲む、ステップと、
‐各当該中心線ポイントに対し当該横断面のエリアを算出するステップ
によって算出される、請求項1から10の何れか1項に記載の方法。 - 分離領域内の周囲輪郭線の最小エネルギー表面が、
‐当該周囲輪郭線の各輪郭線頂点から当該周囲等輪郭線の前記重心への直線を定義し、かつ当該線を、mサンプルを含むnリングに分割するステップであって、当該周囲輪郭線の前記頂点が、前記第一のリングの前記サンプルであり、かつ各内側リングはその隣の外側リングより少ないサンプルを含み、かつ前記最初の2つ以外には、2個の連続したリングが同じ放射線をサンプリングしない、ステップと、
‐リングのポイントの間で表面メッシュを定義するステップと、
‐当該表面メッシュの合計エリアを算出するステップと、
‐当該第一リングの当該サンプルを固定し、かつメッシュ・エリアの変更および/またはサンプルの動きが、既定のしきい値以下になるまで、平均隣接方向にサンプルを互いに移動し続けるステップ
により、繰返し算出される、請求項1から11の何れか1項に記載の方法。 - 定量分析の結果が、典型的には、ツリーを参照して提示され、当該結果が、長さ、エリア、直径、分岐の間の角度、湾曲率を含む、請求項1から12の何れか1項に記載の方法。
- 最大の内接球体の前記半径を、単一血管領域内の各中心線位置での横断面の前記最小囲み円の前記半径によって分割することによって見出すことができる前記管状器官に渡って健全性を示す真円度グラフを決定するステップも、備える、請求項1から13の何れか1項に記載の方法。
- 健康な再現が、前記ツリーまたは前記ツリーの部分を形成する前記管状器官の前記健康なルーメン壁の前記表面を表示する3Dモデルとなるように、提供されおよび/または算出される、請求項1から14の何れか1項に記載の方法。
- ‐輪郭線モデルとして、前記単一管状器官の健康な再現を作成するステップと、
‐前記N-分岐の健康な再現を作成するステップであって、当該健康な再現の各々が、N個の部分に分割され、当該部分の各々が、近位から遠位に行く最小ルーメン曲率を記述しかつ輪郭線モデルとして当該部分を結合する、ステップと、
‐全ての再現された単一管状器官およびN-分岐を1個の輪郭線モデルに結合するステップと、
‐前記ツリーまたは前記ツリーの一部を形成する前記健康な管状器官を形成する当該輪郭線モデルから健康なモデル表面メッシュを作成するステップと、
‐前記ツリーまたは前記ツリーの一部を形成する前記健康な管状器官に対し3D中心線を定義するステップと、
‐前記ツリーまたは前記ツリーの一部を形成する前記健康な管状器官のN-分岐を定義するステップと、
‐前記健康なモデル・ツリーの連結する中心線を前記ツリーの前記中心線に連結させるステップと、
‐断面を作成しかつN-分岐のエリア曲線を算出するステップと、
‐単一管状器官に対し断面を作成するステップと、
‐算出された参照に基づいて単一管状器官に対しエリア曲線を算出するステップと
を備え、
断面及び/又はエリアが、形成されていて、及び/又は前記健康な3Dモデルに適用される前記請求項の何れかに記載の1個以上のステップにより定義されている、
請求項14または15に記載の方法。 - 前記中心線の連結が、前記真のツリーと同じトポロジーを有する前記健康なモデル・ツリーに基づいている、請求項16に記載の方法。
- 定量的分析の結果が、前記ツリーの再現との関係で前記健康なツリー再現を参照して提示される、請求項15から17の何れか1項に記載の方法。
- 健康な分離領域は、N-分岐の最初からそのN-分岐の終端の各々に行くルーメン壁の曲率を最小化するモデルによって、再現される、請求項15から18の何れか1項に記載の方法。
- 当該3Dモデルが、当該3Dモデルの一つ以上の狭窄された部分を除いて、健康とみなされ、これによって、当該狭窄された部分の各々は、当該ツリーのサブツリーの3Dモデルを反映し、かつ当該狭窄された部分は、自動的に算出される一つ以上またはユーザが提供した病変の位置に基づいて自動的に導出される、請求項15から19の何れか1項に記載の方法。
- 前記健康なモデルが、斑の形状および位置の視覚フィードバックを提供するために使用される、請求項15から20の何れか1項に記載の方法。
- 当該ツリーのN-分岐の何れかから生じる管状器官の一つ以上のサブツリーが、無視される、請求項1から21の何れか1項に記載の方法。
- コンピュータのメモリーに直接ロード可能であるコンピュータ製品であって、前記製品がコンピュータで駆動されると、請求項1から22の何れか1項に記載の方法を実行するためのソフトウエアコード部分を備えるコンピュータ製品。
- 再帰的分離管状器官のツリーまたはツリーの一部を定量的分析するための装置であって、前記ツリーまたは前記ツリーの一部を形成する前記管状器官の前記ルーメン壁の表面を表示する、当該ツリーの3Dモデルまたは当該ツリーの一部を受信しおよび/または算出する手段と、請求項1から22の何れか1項に記載の方法を実行するようにプログラムされている実行手段とを、備える装置。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/EP2010/062883 WO2012028190A1 (en) | 2010-09-02 | 2010-09-02 | Method and apparatus for quantitative analysis of a tree of recursively splitting tubular organs |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2013536706A true JP2013536706A (ja) | 2013-09-26 |
JP5750160B2 JP5750160B2 (ja) | 2015-07-15 |
Family
ID=43882678
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013526324A Active JP5750160B2 (ja) | 2010-09-02 | 2010-09-02 | 再帰的に管状器官が分離するツリーの定量的分析方法および装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8787641B2 (ja) |
EP (1) | EP2559001B1 (ja) |
JP (1) | JP5750160B2 (ja) |
WO (1) | WO2012028190A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7497471B2 (ja) | 2020-06-24 | 2024-06-10 | シャンハイ・パルス・メディカル・テクノロジー,インコーポレイテッド | 血管造影画像処理に用いる訓練方法及び装置、自動処理方法及び装置 |
Families Citing this family (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008107905A2 (en) | 2007-03-08 | 2008-09-12 | Sync-Rx, Ltd. | Imaging and tools for use with moving organs |
US10716528B2 (en) | 2007-03-08 | 2020-07-21 | Sync-Rx, Ltd. | Automatic display of previously-acquired endoluminal images |
US9629571B2 (en) | 2007-03-08 | 2017-04-25 | Sync-Rx, Ltd. | Co-use of endoluminal data and extraluminal imaging |
US8700130B2 (en) | 2007-03-08 | 2014-04-15 | Sync-Rx, Ltd. | Stepwise advancement of a medical tool |
US9968256B2 (en) | 2007-03-08 | 2018-05-15 | Sync-Rx Ltd. | Automatic identification of a tool |
US11197651B2 (en) | 2007-03-08 | 2021-12-14 | Sync-Rx, Ltd. | Identification and presentation of device-to-vessel relative motion |
US11064964B2 (en) | 2007-03-08 | 2021-07-20 | Sync-Rx, Ltd | Determining a characteristic of a lumen by measuring velocity of a contrast agent |
US8542900B2 (en) | 2007-03-08 | 2013-09-24 | Sync-Rx Ltd. | Automatic reduction of interfering elements from an image stream of a moving organ |
US9375164B2 (en) | 2007-03-08 | 2016-06-28 | Sync-Rx, Ltd. | Co-use of endoluminal data and extraluminal imaging |
US9101286B2 (en) | 2008-11-18 | 2015-08-11 | Sync-Rx, Ltd. | Apparatus and methods for determining a dimension of a portion of a stack of endoluminal data points |
US9144394B2 (en) | 2008-11-18 | 2015-09-29 | Sync-Rx, Ltd. | Apparatus and methods for determining a plurality of local calibration factors for an image |
US8855744B2 (en) | 2008-11-18 | 2014-10-07 | Sync-Rx, Ltd. | Displaying a device within an endoluminal image stack |
US9974509B2 (en) | 2008-11-18 | 2018-05-22 | Sync-Rx Ltd. | Image super enhancement |
US11064903B2 (en) | 2008-11-18 | 2021-07-20 | Sync-Rx, Ltd | Apparatus and methods for mapping a sequence of images to a roadmap image |
US10362962B2 (en) | 2008-11-18 | 2019-07-30 | Synx-Rx, Ltd. | Accounting for skipped imaging locations during movement of an endoluminal imaging probe |
US9095313B2 (en) | 2008-11-18 | 2015-08-04 | Sync-Rx, Ltd. | Accounting for non-uniform longitudinal motion during movement of an endoluminal imaging probe |
MY152058A (en) * | 2010-06-21 | 2014-08-15 | Univ Putra Malaysia | A method of constructing at least one 3 dimensional image |
EP2723231A4 (en) | 2011-06-23 | 2015-02-25 | Sync Rx Ltd | LUMINAL BACKGROUND CLEANING |
EP2570079B1 (en) | 2011-09-13 | 2017-06-14 | Pie Medical Imaging BV | Method and apparatus for determining optimal 3D reconstruction of an object |
CN103957832B (zh) * | 2011-10-26 | 2016-09-28 | 皇家飞利浦有限公司 | 血管树图像的内窥镜配准 |
JP6134789B2 (ja) | 2012-06-26 | 2017-05-24 | シンク−アールエックス,リミティド | 管腔器官における流れに関連する画像処理 |
JP2014100249A (ja) | 2012-11-19 | 2014-06-05 | Toshiba Corp | 血管解析装置、医用画像診断装置、血管解析方法、及び血管解析プログラム |
US9471989B2 (en) * | 2013-06-03 | 2016-10-18 | University Of Florida Research Foundation, Inc. | Vascular anatomy modeling derived from 3-dimensional medical image processing |
US9974442B2 (en) | 2013-06-24 | 2018-05-22 | Toshiba Medical Systems Corporation | Method of, and apparatus for, processing volumetric image data |
EP3656557B1 (en) * | 2013-10-11 | 2021-12-08 | Advanced Solutions Life Sciences, LLC | System and workstation for the design, fabrication and assembly of bio-material constructs |
EP3061015A2 (en) | 2013-10-24 | 2016-08-31 | Cathworks Ltd. | Vascular characteristic determination with correspondence modeling of a vascular tree |
US10424062B2 (en) * | 2014-08-06 | 2019-09-24 | Commonwealth Scientific And Industrial Research Organisation | Representing an interior of a volume |
WO2017019634A1 (en) | 2015-07-25 | 2017-02-02 | Lightlab Imaging, Inc. | Intravascular data visualization method |
EP4241694A3 (en) | 2016-05-16 | 2023-12-20 | Cathworks Ltd. | Selection of vascular paths from images |
EP3279865B2 (en) | 2016-08-01 | 2022-02-16 | 3mensio Medical Imaging B.V. | Method, device and system for simulating shadow images |
EP3403582B1 (en) | 2017-05-15 | 2023-06-07 | Pie Medical Imaging BV | Method and apparatus for determining blood velocity in x-ray angiography images |
US11080875B2 (en) * | 2018-03-22 | 2021-08-03 | Jvckenwood Corporation | Shape measuring apparatus, shape measuring method, non-transitory computer readable medium storing program |
US10861157B2 (en) | 2019-04-04 | 2020-12-08 | Medtronic Vascular, Inc. | System and methods for determining modified fractional flow reserve values |
US12039685B2 (en) | 2019-09-23 | 2024-07-16 | Cathworks Ltd. | Methods, apparatus, and system for synchronization between a three-dimensional vascular model and an imaging device |
KR102310740B1 (ko) * | 2019-11-12 | 2021-10-12 | 울산대학교 산학협력단 | 구조화된 혈관 분지 정보에 기초한 의료 시술 히스토리 검색 방법 |
DE102022204190A1 (de) * | 2022-04-28 | 2023-11-02 | Carl Zeiss GOM Metrology GmbH | Bestimmung einer geometrischen Begrenzung für einen Massenstrom, für einen Energiestrom oder für einen Kraftfluss durch einen Kanal mit einer Kanalwand |
WO2023237553A1 (en) | 2022-06-07 | 2023-12-14 | Pie Medical Imaging Bv | Method and system for assessing functionally significant vessel obstruction based on machine learning |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004535874A (ja) * | 2001-07-18 | 2004-12-02 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 磁気共鳴血管造影法及びその装置 |
JP2006075599A (ja) * | 2004-09-09 | 2006-03-23 | Siemens Ag | 管状の解剖学的構造の可視化方法 |
JP2007506531A (ja) * | 2003-09-25 | 2007-03-22 | パイエオン インコーポレイテッド | 管状臓器の3次元再構成を行うシステム及び方法 |
JP2007289687A (ja) * | 2006-04-21 | 2007-11-08 | Siemens Medical Solutions Usa Inc | 血管の3次元画像を作成する方法 |
JP2008529641A (ja) * | 2005-02-11 | 2008-08-07 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 3次元医療画像からの肺動脈樹の自動抽出の方法 |
WO2009102930A2 (en) * | 2008-02-13 | 2009-08-20 | Kitware, Inc. | Method and system for measuring tissue damage and disease risk |
JP2010042247A (ja) * | 2008-07-22 | 2010-02-25 | Pie Medical Imaging Bv | 複数の2d血管造影画像を用いた3dの定量的分岐解析のための方法、装置およびコンピュータプログラム |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6251072B1 (en) | 1999-02-19 | 2001-06-26 | Life Imaging Systems, Inc. | Semi-automated segmentation method for 3-dimensional ultrasound |
JP2004513673A (ja) | 2000-05-09 | 2004-05-13 | ペイエオン・インコーポレーテツド | 動脈の三次元的再構成のシステムと方法 |
US8086000B2 (en) * | 2007-08-27 | 2011-12-27 | Pie Medical Imaging B.V. | Method, apparatus and computer program for quantitative bifurcation analysis on angiographic images |
-
2010
- 2010-09-02 US US13/819,891 patent/US8787641B2/en active Active
- 2010-09-02 WO PCT/EP2010/062883 patent/WO2012028190A1/en active Application Filing
- 2010-09-02 EP EP10757058.2A patent/EP2559001B1/en active Active
- 2010-09-02 JP JP2013526324A patent/JP5750160B2/ja active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004535874A (ja) * | 2001-07-18 | 2004-12-02 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 磁気共鳴血管造影法及びその装置 |
JP2007506531A (ja) * | 2003-09-25 | 2007-03-22 | パイエオン インコーポレイテッド | 管状臓器の3次元再構成を行うシステム及び方法 |
JP2006075599A (ja) * | 2004-09-09 | 2006-03-23 | Siemens Ag | 管状の解剖学的構造の可視化方法 |
JP2008529641A (ja) * | 2005-02-11 | 2008-08-07 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 3次元医療画像からの肺動脈樹の自動抽出の方法 |
JP2007289687A (ja) * | 2006-04-21 | 2007-11-08 | Siemens Medical Solutions Usa Inc | 血管の3次元画像を作成する方法 |
WO2009102930A2 (en) * | 2008-02-13 | 2009-08-20 | Kitware, Inc. | Method and system for measuring tissue damage and disease risk |
JP2010042247A (ja) * | 2008-07-22 | 2010-02-25 | Pie Medical Imaging Bv | 複数の2d血管造影画像を用いた3dの定量的分岐解析のための方法、装置およびコンピュータプログラム |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7497471B2 (ja) | 2020-06-24 | 2024-06-10 | シャンハイ・パルス・メディカル・テクノロジー,インコーポレイテッド | 血管造影画像処理に用いる訓練方法及び装置、自動処理方法及び装置 |
Also Published As
Publication number | Publication date |
---|---|
US8787641B2 (en) | 2014-07-22 |
EP2559001A1 (en) | 2013-02-20 |
EP2559001B1 (en) | 2018-07-11 |
JP5750160B2 (ja) | 2015-07-15 |
WO2012028190A1 (en) | 2012-03-08 |
US20130158970A1 (en) | 2013-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5750160B2 (ja) | 再帰的に管状器官が分離するツリーの定量的分析方法および装置 | |
US11816837B2 (en) | Vascular characteristic determination with correspondence modeling of a vascular tree | |
Kreiser et al. | A survey of flattening‐based medical visualization techniques | |
CN100466974C (zh) | 使用n-MIP光线投射技术的自动冠状动脉分离 | |
Schaap et al. | Standardized evaluation methodology and reference database for evaluating coronary artery centerline extraction algorithms | |
JP2022169579A (ja) | リアルタイムで診断上有用な結果 | |
US7599535B2 (en) | System and method for tree-model visualization for pulmonary embolism detection | |
US7620225B2 (en) | Method for simple geometric visualization of tubular anatomical structures | |
JP4521271B2 (ja) | 変形可能な3次元対象の頭頂の動きに関する情報を表示する画像処理システム | |
US7583829B2 (en) | Method and apparatus for embolism analysis | |
KR20140141526A (ko) | 석회화된 혈관의 분할 | |
JP2008510499A (ja) | 解剖学的可視化/測定システム | |
RU2419882C2 (ru) | Способ визуализации секущих плоскостей для изогнутых продолговатых структур | |
JP2012511380A (ja) | 血管分析 | |
JP2010528750A (ja) | 管状構造の検査 | |
US20200118264A1 (en) | Shell-constrained localization of vasculature | |
Gotra et al. | Validation of a semiautomated liver segmentation method using CT for accurate volumetry | |
Luo et al. | IVUS validation of patient coronary artery lumen area obtained from CT images | |
Lidayova et al. | Skeleton-based fast, fully automated generation of vessel tree structure for clinical evaluation of blood vessel systems | |
KR101514795B1 (ko) | 혈관의 정량화 방법 | |
Kiraly et al. | Analysis of arterial subtrees affected by pulmonary emboli | |
Ghaffari et al. | Validation of parametric mesh generation for subject-specific cerebroarterial trees using modified Hausdorff distance metrics | |
JPWO2018042870A1 (ja) | Cpr画像生成装置、方法およびプログラム | |
Al Moussawi et al. | 3D Reconstruction of blood vessels | |
CN111798468B (zh) | 图像处理方法及装置、存储介质及电子终端 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130819 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20140122 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140128 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20140428 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20140508 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140527 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20141118 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150129 |
|
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: 20150507 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20150515 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5750160 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 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |