(57)【要約】【目的】動脈瘤等の治療に関する医療上の多面的かつ有用な情報が得られることを課題とする。【構成】血管のＣＴ断層像に基づき血管部ＣＴ断面の中心点Ｏを求め、隣接する各血管部ＣＴ断面の中心点Ｏ１〜Ｏｎ間を滑らかな曲線で補間して血管部中心線Ｃを形成するステップと、前記中心線Ｃに垂直な血管部垂直断面を該中心線Ｃの微小区間毎に作成すると共に、該血管部垂直断面の輪郭形状を多角形で近似するステップと、隣接する多角形断面の各頂点間につき求めた微小距離を夫々累積加算することにより血管壁の長さを求めるステップとを備える。【選択図】    図１
2004201873.tif 000002
【特許請求の範囲】
【請求項１】
血管のＣＴ断層像に基づき血管部ＣＴ断面の中心点を求め、隣接する各血管部ＣＴ断面の中心点間を滑らかな曲線で補間して血管部中心線を形成するステップと、
前記中心線に垂直な血管部垂直断面を該中心線の微小区間毎に作成すると共に、該血管部垂直断面の輪郭形状を多角形で近似するステップと、
隣接する多角形断面の各頂点間につき求めた微小距離を夫々累積加算することにより血管壁の長さを求めるステップとを備えることを特徴とする画像処理方法。
【請求項２】
前記求めた血管壁の長さに基づき血管壁の最大長と最小長とを求めるステップを更に備えることを特徴とする請求項１記載の画像処理方法。
【請求項３】
血管部垂直断面につき求めた面積に、隣接する垂直断面間の距離を乗算して微小区間毎の血管部微小体積を求めると共に、該求めた各微小体積を加算して血管部全体の体積を求めるステップと、
前記ＣＴ断層像により作成した３次元画像につき所定のＣＴ閾値範囲内に含まれる血栓部のボクセルデータを抽出・計数して、該血栓部の体積を求めるステップと、
前記求めた血管部全体の体積と血栓部体積との比を求めるステップとを備えることを特徴とする請求項１記載の画像処理方法。
【請求項４】
血管のＣＴ断層像に基づき血管部ＣＴ断面の中心点を求め、隣接する各血管部ＣＴ断面の中心点間を滑らかな曲線で補間して血管部中心線を形成する中心線形成手段と、
前記中心線に垂直な血管部垂直断面を該中心線の微小区間毎に作成すると共に、該血管部垂直断面の輪郭形状を多角形で近似する垂直断面作成手段と、
隣接する多角形断面の各頂点間につき求めた微小距離を夫々累積加算することにより血管壁周囲の長さを求める血管長演算手段とを備えることを特徴とする画像処理装置。
【請求項５】
血管長演算手段は、求めた血管壁の長さに基づき血管壁の最大長と最小長とを求めることを特徴とする請求項４記載の画像処理装置。
【請求項６】
垂直断面作成手段は、隣接する垂直断面の一部が交差しない範囲内で中心線上の各微小区間を選択することを特徴とする請求項４記載の画像処理装置。
【請求項７】
血管部垂直断面につき求めた面積に、隣接する垂直断面間の距離を乗算して微小区間毎の血管部微小体積を求めると共に、該求めた各微小体積を加算して血管部全体の体積を求める血管部体積演算手段と、
前記ＣＴ断層像に基き血栓部の体積を求める血栓部体積演算手段と、
前記求めた血管部体積と血栓部体積との比を求める体積比演算手段とを備えることを特徴とする請求項６記載の画像処理装置。
【請求項８】
血管部両端面部の垂直断面につき第１，第２の血管径を求めると共に、中間の微小区間毎の各血管径を前記第１，第２の血管径の中間の血管径で滑らかに補間する血管径補間手段と、
前記補間された各血管径を有する仮想の血管部画像を作成する血管部作成手段とを備えることを特徴とする請求項６記載の画像処理装置。
【請求項９】
血管部垂直断面につき求めた面積に、隣接する垂直断面間の距離を乗算して微小区間毎の血管部微小体積を求めると共に、該求めた各微小体積を加算して血管部全体の体積を求める血管部体積演算手段と、
前記ＣＴ断層像に基き血栓部の体積を求める血栓部体積演算手段と、
前記求めた血管部体積と血栓部体積との比を求める体積比演算手段とを備えることを特徴とする請求項８記載の画像処理装置。
【請求項１０】
血栓部体積演算手段は、ＣＴ断層像により作成した３次元画像につき所定のＣＴ閾値範囲内に含まれる血栓部のボクセルデータを抽出・計数して、該血栓部の体積を求めることを特徴とする請求項７又は９に記載の画像処理装置。
【請求項１１】
血管部体積演算手段は、隣接する垂直断面間の距離を該隣接する多角形頂点間の最大の距離と最小の距離との和の１／２により求めることを特徴とする請求項７又は９に記載の画像処理装置。
【請求項１２】
コンピュータに請求項１乃至３の何れか１に記載の画像処理方法を実行させるためのプログラム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は画像処理方法及び画像処理装置並びにプログラムに関し、更に詳しくは、血管造影剤を使用して撮影した動脈瘤のＸ線ＣＴ画像解析処理に適用して好適なるものである。
【０００２】
【従来の技術】
従来より、動脈瘤の撮影目的で血管造影剤を用いたＸ線ＣＴ撮影が行われている。動脈瘤は、血管部と、この血管壁に付着し、ゼリー（jelly）状に血液が固まった血栓部分とからなっている。マルチディテクタを使用したＸ線ＣＴ撮影では被検体体軸方向の分解能を高めて血管造影を行うことができ、動脈瘤部の血管領域抽出を容易に行うことができる。
【０００３】
動脈瘤の定量的な評価は、正常な血管部の径と瘤部の径とを計測／比較することにより行われる。この場合に、血管部のＣＴ断層像は血中の造影剤により鮮明に得られるが、血管壁に付着した血栓部分には造影剤が流入しないため、鮮明な像は得られない。しかも、血栓部には濃度斑と呼ばれる明るさの濃淡が存在しており、これが画像処理技術による血栓部分の計測を妨げている。
【０００４】
この様な状況の下、従来は、造影領域を抽出するための閾値を適切に求める目的で、造影剤を注入しないで撮影したプレーンスカウト像と、造影剤を注入して撮影したアキシャル画像に基づき作成した造影スカウト像との差分をとって差分画像を求めると共に、該差分画像の造影領域の画素値に基づき造影領域抽出用の閾値を決定し、該閾値で造影アキシャル画像から造影領域を抽出し、鮮明な３Ｄ画像を作成するものが知られている（特許文献１）。
【０００５】
【特許文献１】
特開２０００−３２５３４０(要約）。
【０００６】
【発明が解決しようとする課題】
しかし、上記造影領域（即ち、血管領域）の３Ｄ画像が鮮明に観察できるだけでは、動脈瘤等の治療に関する医療上の多面的かつ有用な情報は得られない。
【０００７】
本発明は上記従来技術の問題点に鑑みなされたもので、その目的とする所は、動脈瘤等の治療に関する医療上の多面的かつ有用な情報が得られる画像処理方法及び画像処理装置並びにプログラムを提供することにある。
【０００８】
【課題を解決するための手段】
上記の課題は例えば図１の構成により解決される。即ち、本発明（１）の画像処理方法は、血管のＣＴ断層像に基づき血管部ＣＴ断面の中心点Ｏを求め、隣接する各血管部ＣＴ断面の中心点Ｏ１〜Ｏｎ間を滑らかな曲線で補間して血管部中心線を形成するステップと、前記中心線Ｃに垂直な血管部垂直断面を該中心線Ｃの微小区間毎に作成すると共に、該血管部垂直断面の輪郭形状を多角形で近似するステップと、隣接する多角形断面の各頂点間につき求めた微小距離を夫々累積加算することにより血管壁の長さを求めるステップとを備えるものである。
【０００９】
従って、様々な形状を有する動脈瘤（特に血管部）の形状、寸法を立体的、かつ忠実に把握できると共に、動脈瘤の多面的な診断が可能となる。なお、上記隣接する多角形断面の各頂点間につき求めた微小距離とは、好ましくは、隣接する多角形断面の各頂点間を滑らかな曲線で補間した場合の微小距離である。
【００１０】
本発明（２）では、上記本発明（１）において、前記求めた血管壁の長さに基づき血管壁の（長手方向の）最大長と最小長とを求めるステップを更に備える。従って、血管の曲がりのみならず、血管のねじれも含めたより実際的な形状、寸法を忠実に把握できる。
【００１１】
本発明（３）では、上記本発明（１）において、血管部垂直断面につき求めた面積に、隣接する垂直断面間の距離を乗算して微小区間毎の血管部微小体積を求めると共に、該求めた各微小体積を加算して血管部全体の体積を求めるステップと、前記ＣＴ断層像により作成した３次元画像につき所定のＣＴ閾値範囲内に含まれる血栓部のボクセルデータを抽出・計数して、該血栓部の体積を求めるステップと、前記求めた血管部全体の体積と血栓部体積との比を求めるステップとを備えるものである。従って、様々な形状を有する動脈瘤の立体的、定量的な診断が容易に行える。
【００１２】
また本発明（４）の画像処理装置は、血管のＣＴ断層像に基づき血管部ＣＴ断面の中心点を求め、隣接する各血管部ＣＴ断面の中心点間を滑らかな曲線で補間して血管部中心線を形成する中心線形成手段と、前記中心線に垂直な血管部垂直断面を該中心線の微小区間毎に作成すると共に、該血管部垂直断面の輪郭形状を多角形で近似する垂直断面作成手段と、隣接する多角形断面の各頂点間につき求めた微小距離を夫々累積加算することにより血管壁周囲の長さを求める血管長演算手段とを備えるものである。
【００１３】
本発明（５）では、上記本発明（４）において、血管長演算手段は、求めた血管壁の長さに基づき血管壁の（長手方向の）最大長と最小長とを求める。
【００１４】
本発明（６）では、上記本発明（４）において、垂直断面作成手段は、隣接する垂直断面の一部が交差しない範囲内で中心線上の各微小区間を選択する。従って、血栓部を除く血管部の全体を重複なく、効率よく微小体積区間に分割できる。
【００１５】
本発明（７）では、上記本発明（６）において、血管部垂直断面につき求めた面積に、隣接する垂直断面間の距離を乗算して微小区間毎の血管部微小体積を求めると共に、該求めた各微小体積を加算して血管部全体の体積を求める血管部体積演算手段と、前記ＣＴ断層像に基き血栓部の体積を求める血栓部体積演算手段と、前記求めた血管部体積と血栓部体積との比を求める体積比演算手段とを備える。
【００１６】
本発明（８）では、上記本発明（６）において、血管部両端面部の垂直断面につき第１，第２の血管径を求めると共に、中間の微小区間毎の各血管径を前記第１，第２の血管径の中間の血管径で滑らかに補間する血管径補間手段と、前記補間された各血管径を有する仮想の血管部画像を作成する血管部作成手段とを備えるものである。
【００１７】
本発明（８）によれば、計測された血管部両端面部の正常な血管径の情報をもとに、動脈瘤部も正常血管径であったと仮定した場合の正常血管を容易に作成(推定）できる。また、この正常血管の情報をもとに、ステントグラフト(stent graft)法等で用いる代替用血管部品を容易かつ適正に設計できる。更には、動脈瘤部の術前／術後の状態や、その後の経過を予測（シミュレーション）する上で多面的かつ極めて有用な医用情報が豊富に得られる。
【００１８】
本発明（９）では、上記本発明（８）において、血管部垂直断面につき求めた面積に、隣接する垂直断面間の距離を乗算して微小区間毎の血管部微小体積を求めると共に、該求めた各微小体積を加算して血管部全体の体積を求める血管部体積演算手段と、前記ＣＴ断層像に基き血栓部の体積を求める血栓部体積演算手段と、前記求めた血管部体積と血栓部体積との比を求める体積比演算手段とを備える。従って、正常であった場合の血管部と、血栓部との様々な関係を容易に見積もれる。
【００１９】
本発明（１０）では、上記本発明（７）又は（９）において、血栓部体積演算手段は、ＣＴ断層像により作成した３次元画像につき所定のＣＴ閾値範囲内に含まれる血栓部のボクセルデータを抽出・計数して、該血栓部の体積を求める。従って、もとのＣＴ断層像情報を利用できると共に、比較的容易な演算で血栓部の体積を計測できる。なお、この血栓部の体積に上記求めた血管部の体積も含めることで比較の対象を動脈瘤部の体積としてもよい。
【００２０】
本発明（１１）では、上記本発明（７）又は（９）において、血管部体積演算手段は、隣接する垂直断面間の距離を該隣接する多角形頂点間の最大の距離と最小の距離との和の１／２により求める。従って、簡単な演算により各微小体積をより正確に求められる。なお、上記隣接する多角形頂点間の最大の距離と最小の距離とは、好ましくは、隣接する多角形断面の各頂点間を滑らかな曲線で補間した場合の最大の距離と最小の距離である。
【００２１】
また本発明（１２）のプログラムは、コンピュータに上記本発明（１）乃至（３）の何れか１に記載の画像処理方法を実行させるためのプログラムである。この様なプログラムは、ＣＤ−ＲＯＭ等の記録媒体を介して、又は通信ネットワークを介するオンラインによりユーザに提供可能である。
【００２２】
【発明の実施の形態】
以下、添付図面に従って本発明に好適なる実施の形態を詳細に説明する。なお、全図を通して同一符号は同一又は相当部分を示すものとする。
【００２３】
図２は実施の形態によるＸ線ＣＴ装置の要部構成図であり、本装置は大きく分けて、Ｘ線ファンビームＸＬＦＢにより被検体１００のアキシャル／ヘリカルスキャン・読取等を行う走査ガントリ部３０と、被検体１００を搭載して体軸ＣＬｂの方向に移動させる撮影テーブル２０と、操作者の設定・操作に基づき前記走査ガントリ部３０及び撮影テーブル２０の遠隔制御を行うと共に、走査ガントリ部３０で検出された被検体の投影データ（ＲＡＷデータとも呼ばれる）に基づき該被検体のＣＴ断層像を再構成する操作コンソール部１０とを備える。
【００２４】
更に、走査ガントリ部３０において、４０はＸ線管、４０ＡはＸ線管制御部、５０はＸ線の体軸ＣＬｂ方向の曝射幅を制限するコリメータ、５０Ａはコリメータ制御部、９０はチャネルＣＨ方向に並ぶ多数（例えばｎ＝１０００程度）のＸ線検出素子が体軸ＣＬｂ方向の例えば４列Ｌ１〜Ｌ４に配列されている多列Ｘ線検出器（マルチディテクタとも呼ばれる）、９１はＸ線検出器９０の検出信号に基づき被検体１００の投影データｇ1（Ｘ，θ）〜ｇ4（Ｘ，θ）を生成し、収集するデータ収集部、３０Ａは上記Ｘ線撮影系を体軸ＣＬｂの回りに回転させる回転制御部である。なお、慣例に従い、被検体体軸ＣＬｂの方向は装置座標のｚ軸と一致している。
【００２５】
また、操作コンソール部１０において、１１はＸ線ＣＴ装置のスキャン制御並びにスキャン計画やＣＴ断層像の再構成処理等を行う中央処理装置、１１ａはそのＣＰＵ、１１ｂはＣＰＵ１１ａが使用するＲＡＭ，ＲＯＭ等からなる主メモリ（ＭＭ）、１２はキーボードやマウス等を含む指令やデータの入力装置、１３はスキャン計画情報や再構成されたＣＴ断層像等を表示するための表示装置（ＣＲＴ）、１４はＣＰＵ１１ａと走査ガントリ部３０及び撮影テーブル２０等との間で各種制御信号ＣＳやモニタ信号ＭＳのやり取りを行う制御インタフェース、１５はデータ収集部９１からの投影データを一時的に蓄積するデータ収集バッファ、１６は再構成されたＣＴ断層像データを蓄積・格納すると共に、Ｘ線ＣＴ装置の運用に必要な各種アプリケーションプログラムや各種演算／補正用のデータファイル等を格納している２次記憶装置（ハードディスク装置等）である。
【００２６】
そして、２００は病院内の１又は２以上の上記Ｘ線ＣＴ装置と遠隔の画像処理装置（ワークステーションとも呼ばれる）との間をオンラインで接続するためのローカルエリアネットワーク（ＬＡＮ）、１７はＬＡＮ２００に接続するためのＬＡＮインタフェース（ＬＩＦ）、５０は本発明による画像処理を行うワークステーション（ＷＳ）である。なお、このワークステーション５０は本発明の画像処理装置に相当する。
【００２７】
次に血管造影剤を使用したＸ線ＣＴ撮影処理を説明する。図３は実施の形態によるＸ線ＣＴ撮影処理のフローチャートであり、ＣＰＵ１１ａの制御下で実行される。図において、好ましくは、事前に被検体１００のスカウトスキャンを行った後、この処理に入力する。ステップＳ１１では、続く被検体１００のアキシャル／ヘリカルスキャンのためのスキャン計画画面を表示部１３に表示する。ステップＳ１２では、操作者がスキャン範囲、スキャン密度、Ｘ線強度等からなるスキャンパラメータを設定する。
【００２８】
次に、操作者が画面上のリコンタグをクリックすると、Ｘ線ＣＴ断層像の再構成パラメータを入力可能となる。ステップＳ１３では、操作者がＣＴ断層像のリコン範囲、リコン密度等からなるリコンパラメータを設定する。ステップＳ１４ではこれらの設定確認ボタン「ＣＯＮＦＩＲＭ」の入力を待ち、該ボタンが入されない場合は、ステップＳ１２に戻って各種パラメータ情報の設定変更が可能である。
【００２９】
やがて、ステップＳ１４で設定確認ボタン「ＣＯＮＦＩＲＭ」が入力されると、ステップＳ１５では上記設定されたスキャンパラメータに従って被検体１００のスキャン制御を行う。その際には、予め別途に設定されたシーケンスやタイミングに従って血管造影剤を投入する。ステップＳ１６では被検体１００の投影データｇ1（Ｘ，θ）〜ｇ4（Ｘ，θ）を収集・蓄積する。ステップＳ１７では所要撮影領域についての全スキャンを完了したか否かを判別し、完了でない場合はステップＳ１５に戻る。
【００３０】
こうして、やがて、全スキャンを完了すると、ステップＳ１８では上記設定されたリコンパラメータに従って被検体１００のＣＴ断層像を再構成する。ステップＳ１９では再構成されたＣＴ断層像をディスク１６に保存する。そして、ステップＳ２０では該作成したＣＴ断層像を表示装置１３に表示する。
【００３１】
図４は実施の形態による動脈瘤解析処理のフローチャートであり、ワークステーション５０のＣＰＵ（不図示）により実行される。また、図５〜図８は動脈瘤解析処理のイメージ図（１）〜（４）であり、以下、これらの図を参照して動脈瘤解析処理を詳細に説明する。図４において、造影剤を使用した動脈瘤のＣＴ画像データが入手可能になるとこの処理に入力する。
【００３２】
ステップＳ３１ではＸ線ＣＴ装置よりＬＡＮ２００を介して入力した動脈瘤を含む血管のＣＴ画像データに基き公知の処理により該血管部の３Ｄ画像を作成する。図５（Ａ）に血管部３Ｄ画像のイメージを示す。被検体体軸（ｚ軸）方向のＣＴ断層像を重ね合せて表示することにより血管部の３Ｄ画像が得られる。ここで、６１は血管部、６２は血栓部である。このままでも、動脈瘤の３Ｄ画像を観察できるが、診断や治療に有用な幾何学的情報は得られない。
【００３３】
ステップＳ３２では操作者の指示に従い注目すべき血管領域を抽出し、血管部６１の中心線を作成する。図５（Ｂ）に中心線Ｃの３Ｄ画像を示す。この場合に、上記図５（Ａ）の３Ｄ画像データを使用して血管６１の中心線Ｃを作成する処理としては、血管断面の面積重心を求める等の、公知の方法を採用できるが、ここでは断面中心を高速に検出可能な一例の方法を説明する。この処理はワークステーション５０のＣＰＵにより自動的に行われる。
【００３４】
図６に血管部中心線作成処理のイメージを示す。図６（Ａ）に血管部６１のあるＣＴスライス面の平面図を示す。一般に、血管部のスライス断面は円形又は楕円形をしているが、ここでは一般性を損なうことなく楕円形の場合を示す。まず直交線分Ｘ，Ｙを容易し、その中心Ｏが血管部６１の内側に来るようにセットする。線分中心Ｏが血管部６１の内側にある判定は、各線分長ＰＯ≠０，ＯＱ≠０，ＳＯ≠０，ＯＴ≠０であり、このとき、通常は、線分長ＰＯ≠ＯＱ，ＳＯ≠ＯＴである。なお、直交線分Ｘ，Ｙと血管部断面輪郭との交点Ｐ，Ｑ，Ｓ，Ｔは、血管部６１の造影撮影により、鮮明に得られる。
【００３５】
次に図６（Ｂ）では例えば線分長ＰＯ＝ＯＱになるまで直交線分Ｘ，Ｙをｘ軸方向（矢印ａ方向）にスライドする。このとき、線分長ＳＯ≠ＯＴのままである。そこで、次に図６（Ｃ）では線分長ＳＯ＝ＯＴになるまで直交線分Ｘ，Ｙをｙ軸方向（矢印ｂ方向）にスライドする。この時点で、もし両線分長ＰＯ，ＯＱが略等しい（許容範囲内である）状態に維持されていれば、その時の中心点Ｏの座標が血管６１の中心座標となる。特に、血管断面が円形の場合はこの時点で血管中心Ｏが求まる。しかし、図示の如く楕円形の場合は更に上記の処理を繰り返す必要がある。
【００３６】
即ち、この時点でまだ線分長ＰＯ≠ＯＱの場合は、更に図６（Ｄ）において、該線分長ＰＯ＝ＯＱになるまで直交線分Ｘ，Ｙをｘ軸方向（矢印ａ方向）にスライドする。こうして、血管部６１の中心を比較的に高速に検出できる。
【００３７】
そして、こうして求めた各血管部中心座標につき、隣接するＣＴ断面の中心点間を滑らかな曲線で補間することにより、血管部中心線Ｃを形成する。滑らかな曲線としてはＢスプライン曲線等を利用できる。
【００３８】
図４に戻り、ステップＳ３３では上記形成した中心線Ｃに沿って該中心線Ｃに垂直な断面を作成する。図７に血管部垂直断面の作成イメージを示す。ここで、図７（Ａ）に上記図５（Ｂ）の血管部ＣＴ断面を再度示す。また挿入図（ａ）に血管部垂直断面作成処理のイメージを示す。図（ａ）において、今、中心点Ｏｎの血管部断面に注目すると、該中心点Ｏｎのところで中心線Ｃを微分することにより接線の単位ベクトル<ｕ１>を得る。次に該ベクトル<ｕ１>に直交するベクトル<ｕ２>を求める。直交ベクトルを求める方法は公知であり、例えばグラムシュミットのベクトル直交化法を利用できる。そして、該ベクトル<ｕ２>に平行な多数の線分を想定すると共に、該線分と血管部周囲との交点座標をプロットすることにより、中心線Ｃに垂直な断面を形成できる。この場合に、血管壁部は周囲とのＣＴ値の相違によって鮮明に識別できる。図７（Ｂ）に血管部垂直断面のイメージを示す。
【００３９】
ステップＳ３４では、こうして得られた血管部垂直断面の幾何学的取り扱いが容易となるように、該断面の輪郭を例えば２０角形等からなる多角形で近似する。図８（Ａ）に多角形で近似された血管部垂直断面のイメージを示す。断面を多角形近似することによりコンピュータ処理（ＣＡＤによる処理）が容易となる。
【００４０】
ステップＳ３５では多角形断面の各頂点間を滑らかな曲線（Ｂスプライン曲線等）により補間する。このとき、各隣接する断面につき、どの頂点間を順次補間・接続するかについては幾分の自由度が存在する。今、ｚ軸方向に展開する血管部に曲がりはあるが、ねじれが無いとすると、隣接する多角形の各頂点はｚ軸方向に略揃っている。この場合は、各多角形断面のｚ軸方向に対応する（振れの少ない）頂点間をねじれ無く接続することにより、ねじれの無い血管部を再生できる。この処理は、ＣＰＵにより自動的に行える。一方、血管部にねじれが存在するような場合には、どの頂点間を接続するかの頂点の選択に人手を介入させてもよい。そして、ステップＳ３６ではこうして形成された血管部６１ａの情報に基き動脈瘤の診断や治療に有用な各種幾何学的情報を抽出する。次にこれらを説明する。
【００４１】
図９は実施の形態による各種医用情報の抽出処理を説明するイメージ図であり、動脈瘤をなす血管部全体の体積と血栓部体積との比を求める場合を示している。図９（Ａ）に血管部垂直断面の一部を示す。今、最初の垂直断面の中心点をＯ１とすると、その下側に第２，第３の各垂直断面に対応する中心点０２，０３等が続いている。最初の垂直断面と第２の垂直断面に着目すると、多角形の左端の頂点Ｐ１間の距離Ｄ１２ａが最大であり、その反対の右側の頂点Ｐｉ間の距離Ｄ１２ｂが最小となっている。この場合に、隣接する第１，第２の各垂直断面は、該断面の一部が交差しない範囲内で中心線Ｃ上の微小区間Ｏ１，Ｏ２が選択されている。更に、この例の中心線Ｃ（即ち、血管部）はｘｚ平面内で蛇行しており、よって多角形の頂点間の距離は図の左側が最大になったり、右側が最大になったりする。図示しないが、中心線Ｃがｙｚ平面内で蛇行する場合も同様に考えられる。更に、実際の血管部の曲がり方の方向は、もっと複雑かつ多様であり、これに応じて頂点間の距離が最大になる頂点の位置も変化する。そしてその反対側に頂点間の距離が最小となる頂点が表れる。
【００４２】
いずれにしても、隣接する多角形断面の各頂点間につき求めた微小距離を夫々累積加算することにより血管壁周囲の長さを求めることができる。例えば、Ｄ１２ａ＋Ｄ２３ａ＋Ｄ３４ｂ＋…の演算により血管部左面の長さが求まる。なお、図示の各微小距離Ｄ１２ａ，Ｄ２３ａ，Ｄ３４ｂ等は直線であるが、これに限らない。好ましくは、各微小距離Ｄ１２ａ，Ｄ２３ａ，Ｄ３４ｂ等は隣接する多角形断面の各頂点間を滑らかな曲線で補間した場合の微小距離とする。また、血管部にねじれが存在するような場合には、最大と最小の各微小距離の表れる位置が僅かな角度をずらしながら表れる。この場合はＤ１２ａ’＋Ｄ２３ａ’＋Ｄ３４ｂ’＋…の演算が行われる場合もある。この各微小距離Ｄ１２ａ’，Ｄ２３ａ’，Ｄ３４ｂ’等についても滑らかな曲線で補間したものを採用できる。また、他の血管壁の長さも同様にして求まる。かくして、動脈瘤をなす血管部の有用な幾何学的情報が得られる。
【００４３】
次に血管部体積の計算方法を説明する。図９（Ａ）において、血管部垂直断面につき求めた面積に、隣接する垂直断面間の距離を乗算して微小区間毎の血管部微小体積を求めると共に、該求めた各微小体積を中心線Ｃの方向に順次加算して血管部全体の体積を求める。この場合に、隣接する垂直断面間の距離は、隣接する多角形頂点間の最大の距離（例えばＤ１２ａ）と最小の距離（例えばＤ１２ｂ）との和の１／２により求める。従って、簡単な演算により各微小体積をより正確に求められる。
【００４４】
図９（Ｂ）に多角形の面積を求める方法のイメージを示す。今、中心点０から見た各頂点Ｐ１〜Ｐｍの位置ベクトルを<Ｐ１>〜<Ｐｍ>とすると、多角形の面積Ｓは次式の演算により求まる。
【００４５】
【数１】
2004201873.tif 000003
【００４６】
ここで、×はベクトルの外積、またｉ＝ｍの時のベクトル<Ｐi+1>＝<Ｐ１>である。なお、ＣＴ断層画像につき、造影剤部分のボクセルデータの個数をカウントする方法もあるが、そのカウント領域が複数のＣＴ断面を横切る垂直断面によって複雑な形状による制限を受けるため、上記数１による演算の方が簡単である。
【００４７】
一方、血栓部６２の体積については、ＣＴ断層像により作成した３次元画像につき血栓部組織のＣＴ閾値範囲内に含まれる血栓部のボクセルデータをＣＴ断層像毎に抽出し、累積計数することにより、該血栓部の体積を比較的容易に求めることが可能である。なお、隣接するＣＴ断層像間に距離がある場合は該距離分を加味してその微小体積分のボクセルデータを計数する。そして、前記求めた血管部全体の体積と血栓部体積との比を求めることにより、様々な形状を有する動脈瘤の立体的、定量的な診断情報が容易に得られる。
【００４８】
図１０は実施の形態によるステントグラフト(stent graft)医療への適用例を説明するイメージ図である。ステントグラフトは、柔らかい金属（ステント：stent）の回りに皮（グラフト：graft）をつけてできた代替用血管部品で、動脈瘤や奇形血管部の治療によく用いられる。ステントグラフト法では、カテーテルを使用して人工血管をはめ込むため、人工血管の寸法(血管径や長さ等）が合わないとうまくはまらない。従って、まず治療を受ける血管部に関する正確な寸法情報が必要である。この寸法情報は上記実施の形態により正確に得られた。これに対して、新たに置換すべき正常な血管部品（ステントグラフト）の詳細な寸法情報も必要である。
【００４９】
本実施の形態では、図１０に示す如く、正常な血管部両端面部の各垂直断面につき第１，第２の血管径を求めると共に、中間の微小区間毎の各血管径を前記第１，第２の血管径の中間の血管径で滑らかに補間することにより、該滑らかに変化する各血管径を有するような仮想の正常な血管部品６１ｂを設計する。好ましくは、正常な血管径は、正常部の複数垂直断面の血管径の平均により求める。
【００５０】
また、この例の血管部品６１ｂでは中央部の血管径が上記図８（Ｂ）のものに比べて太くなっており、よって置換後には血栓部６２による悪影響が除去される。なお、術後の血栓部は次第に消滅してゆくことが知られている。
【００５１】
この様に、本実施の形態では、計測された血管部両端面部の正常な血管径の情報をもとに、動脈瘤部も正常血管径であったと仮定した場合の正常血管を容易に作成(推定）できる。また、この正常血管の情報をもとに、ステントグラフト法等で用いる代替用血管部品を容易かつ適正に設計できる。必要なら、公知のＣＡＤ技術を使用して血管部品の様々な変更が行える。更には、動脈瘤部の術前／術後の状態や、手術計画、更には術後の経過等を予測（シミュレーション）する上で多面的かつ極めて有用な医用情報が豊富に得られる。
【００５２】
なお、上記実施の形態では、本発明の経皮的人工血管置換（即ち、カテーテルを用いた内皮的なステントグラフト置換）への応用例を述べたが、これに限らない。本発明は動脈瘤の外科的なバイパス手術等にも応用可能である。
【００５３】
また、上記実施の形態では本発明のＸ線ＣＴ画像への適用例を具体的に述べたが、これに限らない。本発明は核磁気共鳴（ＮＭＲ）法により被検体の空間的な核磁化分布を画像化するところの所謂ＭＲＩ（Magnetic Resonance Imaging）画像にも適用可能である。
【００５４】
また、上記本発明に好適なる実施の形態を述べたが、本発明思想を逸脱しない範囲内で各部の構成、制御、処理及びこれらの組み合わせの様々な変更が行えることは言うまでも無い。
【００５５】
【発明の効果】
以上述べた如く本発明によれば、動脈瘤等の治療に関する医療上の多面的かつ有用な情報が得られ、動脈瘤等の診断、治療に寄与するところが極めて大きい。
【図面の簡単な説明】
【図１】本発明の原理を説明する図である。
【図２】実施の形態によるＸ線ＣＴ装置の要部構成図である。
【図３】実施の形態によるＸ線ＣＴ撮影処理のフローチャートである。
【図４】実施の形態による動脈瘤解析処理のフローチャートである。
【図５】実施の形態による動脈瘤解析処理のイメージ図（１）である。
【図６】実施の形態による動脈瘤解析処理のイメージ図（２）である。
【図７】実施の形態による動脈瘤解析処理のイメージ図（３）である。
【図８】実施の形態による動脈瘤解析処理のイメージ図（４）である。
【図９】実施の形態による各種医用情報の抽出処理を説明するイメージ図である。
【図１０】実施の形態によるステントグラフト医療への適用例を説明するイメージ図である。
【符号の説明】
３０  走査ガントリ部
２０  撮影テーブル
１０  操作コンソール部
５０  ワークステーション（ＷＳ）
１００  被検体
２００  ローカルエリアネットワーク（ＬＡＮ）
【図１】
2004201873.tif 000004
【図２】
2004201873.tif 000005
【図３】
2004201873.tif 000006
【図４】
2004201873.tif 000007
【図５】
2004201873.tif 000008
【図６】
2004201873.tif 000009
【図７】
2004201873.tif 000010
【図８】
2004201873.tif 000011
【図９】
2004201873.tif 000012
【図１０】
2004201873.tif 000013