(57)【要約】  Ｄ次元空間内の複数の点から成るドメインに対応する複数の輝度から成るディジタル画像中の管状構造を配向する方法を提供する。この方法は画像ドメイン内の１つの点を選択し（１０１）、選択された点の近傍において画像の勾配を計算し（１０２）、選択された点における基本構造を計算し（１０２）、選択された点の構造テンソルを求め（１０３）、構造テンソルの固有値を見つける（１０４）ステップを含む。最小固有値（１０５）に対応する固有ベクトルは管状構造と整列する。管状構造と整列する固有ベクトルによって画定される軸を中心とするカートホイール投影を計算（１０６）すればよい。
2007505711.tif 000002
【特許請求の範囲】
【請求項１】
画像がＤ次元空間内の点のドメインに対応する複数の輝度から成るディジタル画像中の管状構造を配向する方法において、
画像ドメイン内の１つの点を選択し；
前記選択された点の近傍において画像の勾配を計算し；
前記選択された点に対応の基本構造テンソルを計算し；
構造テンソルの固有値を見つけ；
前記固有値を分析することによって前記管状構造と整列する固有ベクトルを見出すステップを含むことを特徴とする前記方法。
【請求項２】
前記選択された点を中心とする前記近傍に亘って画像をガウス・カーネル導関数Ｇで畳込むことによって画像の勾配を測定することを特徴とする請求項１に記載の方法。
【請求項３】
前記ガウス・カーネルが画像のサイズよりもほぼ２桁小さい標準偏差σＧを有することを特徴とする請求項２に記載の方法。
【請求項４】
勾配とその転置行列との外積を、検出中のオブジェクトのサイズに対応する標準偏差σTのガウス・カーネルで畳込むことによって構造テンソルを求めることができることを特徴とする請求項１に記載の方法。
【請求項５】
Householder QL分解を行うことによって固有値を見出すことを特徴とする請求項１に記載の方法。
【請求項６】
管状構造と整列する固有ベクトルが最小固有値に対応することを特徴とする請求項１に記載の方法。
【請求項７】
前記選択された点の前記近傍を、前記近傍の２つの最小固有値を比較することによって、管状構造と比較し、前記近傍が局部的に管状構造と似ておれば前記２つの最小値が非近似であるとするステップをも含むことを特徴とする請求項１に記載の方法。
【請求項８】
前記管状構造と整列する固有ベクトルによって画定される軸を中心とするカートホイール投影を計算するステップをも含むことを特徴とする請求項１に記載の方法。
【請求項９】
Ｄ次元空間内の複数の点を含むドメインに相当する複数の輝度から成るディジタル画像中の管状構造を配向する方法において、
画像ドメイン内の１つの点を選択し；
前記選択された点の近傍において、画像のサイズより２桁小さい標準偏差σＧのガウス・カーネル導関数Ｇで画像を畳込むことによって画像の勾配を計算し；
勾配とその転置行列との外積を、検出中の対象物のサイズに対応する標準偏差σＴのガウス・カーネル導関数Ｇで畳込むことによって前記選択された点の構造テンソルを求め；
Householder QL分解を行うことによって構造テンソルの固有値を見つけ；
前記固有値のうちの最小固有値を見出すことによって前記管状構造と整列する対応の固有ベクトルを見つけ；
前記管状構造と整列する固有ベクトルによって確定される軸を中心とするカートホイール投影を計算するステップを含むことを特徴とする管状構造の配向方法。
【請求項１０】
前記選択された点の前記近傍を、前記近傍の２つの最小固有値を比較することによって管状構造と比較し、前記近傍が局部的に管状構造と似ておれば前記２つの最小値が非近似であるとするステップをも含むことを特徴とする請求項９に記載の方法。
【請求項１１】
Ｄ次元空間内の複数の点を含むドメインに相当する複数の輝度から成るディジタル画像中の管状構造を配向する方法ステップを実施するためコンピュータによって実行可能な命令プログラムを組み込まれたコンピュータ可読プログラム記憶装置において、
画像ドメイン内の１つの点を選択し；
前記選択された点の近傍において画像の勾配を計算し；
前記選択された点における基本構造を計算し；
前記選択された点の構造テンソルを求め；
構造テンソルの固有値を見つけ；
前記固有値を分析することによって前記管状構造と整列する固有ベクトルを見つけるステップを含むことを特徴とするコンピュータ可読プログラム記憶装置。
【請求項１２】
前記選択された点を中心とする、画像のサイズより２桁小さい標準偏差σTのガウス・カーネル導関数Ｇで画像を畳込むことによって勾配を計算するステップをも含むことを特徴とする請求項１１に記載のコンピュータ可読プログラム記憶装置。
【請求項１３】
勾配とその転置行列との外積を、検出中の対象物のサイズに対応する標準偏差σＴのガウス・カーネルで畳込むことによって構造テンソルを求めるステップをも含むことを特徴とする請求項１１に記載のコンピュータ可読プログラム記憶装置。
【請求項１４】
Householder QL分解を行うことによって構造テンソルの固有値を見つけるステップをも含むことを特徴とする請求項１１に記載のコンピュータ可読プログラム記憶装置。
【請求項１５】
前記管状構造と整列する、最小固有値と対応の固有ベクトルを見つけるステップをも含むことを特徴とする請求項１１に記載のコンピュータ可読プログラム記憶装置。
【請求項１６】
前記選択された点の前記近傍の２つの最小固有値を比較することによって前記近傍を管状構造と比較し、前記近傍が局部的に管状構造に似ておれば前記２つの最小固有値が互いに非類似であるとするステップをも含むことを特徴とする請求項１１に記載のコンピュータ可読プログラム記憶装置。
【請求項１７】
前記管状構造と整列する固有ベクトルによって確定される軸を中心とするカートホイールを計算するステップをも含むことを特徴とする請求項１１に記載のコンピュータ可読プログラム記憶装置。
【発明の詳細な説明】
【関連米国出願相互参照】
【０００１】
本願は、２００３年９月２２日付Pascal Cathierの米国仮出願第６０/５０４，９１０号“脈管構造の局部画像の自動的な配向”に基づく優先権を主張するものであり、この仮出願の内容を本願の一部として引用する。
【技術分野及び背景】
【０００２】
本発明はディジタル医療用画像における脈管構造の検出に係わる。
【０００３】
最新の撮像システムからのデータから得られる診断上有効な情報によって、潜在的に問題のある部位を早期且つ治療の可能性が高い段階で検出することができる。撮像システムから大量の詳細データが与えられたら、画像データを有効且つ正確に処理するため種々のアルゴリズムを開発しなければならない。コンピュータの助けにより、画像処理は通常、ディジタル画像またはディジタル化した画像を対象に行われる。
【０００４】
ディジタル画像は特定配列箇所で表される一連の生体構造部位点と対応する（グレースケールまたは磁場強さのような）特性を表す一連の数値から形成される。これら一連の生体構造部位点が画像のドメインを構成する。２-Ｄディジタル画像またはスライス断面において、不連続に配列された場所をピクセルと呼称する。３-Ｄディジタル画像は公知の種々の構成法によって積重ねられたスライス断面から形成することができる。３-Ｄ画像は２-Ｄ画像のピクセルから構成されるボクセルとも呼称される不連続な体積要素から成る。ピクセルまたはボクセルの特性を処理することによって、このようなピクセルまたはボクセルと関連する患者の生体構造部位に関する種々の特性を確認することができる。コンピュータ援用診断（“ＣＡＤ”）システムはディジタル画像データの分析及び可視化において重要な役割を演ずる。
【０００５】
コンピュータ断層（ＣＴ）撮影システム、磁気共鳴（ＭＲ）撮像システム及び３-ＤのＸ線（ＸＲ）撮影システムの重要な役割は、例えば、気道、管路、神経、血管等のような種々の蛇行性管状構造の分析を含む脈管分析のための３Ｄ画像データ・セットを作成することにある。このような３Ｄ画像データ・セットの作成は放射線技師にとって特に重要である。即ち、放射線技師は狭窄または動脈瘤パラメータの判定、長さ、断面寸法、角度、及び関連のパラメータの定量化を可能にする詳細な視覚的レポートの提供を要求される。特定の脈管部位における急性狭窄、特定の脈管部位における広範囲の動脈瘤、または脈管のねじれに関する情報は、医師が手術の段取りを決めるのに利用されるのが普通である。生産性の観点から、また、フィルム・コストを削減するためにも、３Ｄ画像データ・セットは少数の有意義な画像に的を絞るべきである。
【０００６】
脈管分析に役立つ情報を能率的に且つ容易に取得できるようにするため、従来の医療用撮像システムは３Ｄ可視化ソフトウェアを採用することがある。このようなソフトウェアは撮像システム自体に組み込まれるか、または分析ワークステーションに組み込まれ、長さ、角度または体積を測定し、例えば、断面、ナビゲータまたは体積レンダリングを利用する等の種々の方法で体積部を可視化する１組のツールを提供する。特に脈管の分析に関しては、ソフトウェアを利用することによって、特定脈管の複数の斜めスライスを得、これによって脈管を分析することができる。
【０００７】
気道、脈管、管路または神経のような蛇行性構造を分析することは、医療用撮像システムの重要な役割の１つである。今日では、このタスクは複数の斜めスライスを利用してこれらの構造の局部を分析することによって達成される。それぞれの斜めスライスはそれぞれのオブジェクトの短い断面の鮮明且つ歪みのない画像を提供するが、それぞれの全長に亘る画像を提供することは稀である。湾曲した矯正画像はこれらの管状対象物の全長を捕捉し、従って、この分析タスクに好適な合成画像を提供する。これらの画像から軸沿いの真の３Ｄ長さ測定値が得られ、測定値は多くの場合実際の生体構造と近似である。湾曲した矯正画像は湾曲部に沿う等間隔の点における値をサンプリングして線分を形成し、これらの線分をサンプリング・ベクトルによって変換して次の画像線を形成する。
【０００８】
従って、医療用撮像システム及び関連の３Ｄ可視化ソフトウェアが従来よりも効率的で、一貫性があり、リピート可能な、迅速且つオペレータの技能差に影響されない態様で有用な３Ｄ画像データ・セットを作成できるようにする新しい方法及び装置が開発されれば有益である。このような新しい方法及び装置が、管状構造及びその狭窄、動脈瘤及び歪みの分析及び撮像を含む脈管分析を容易にするなら、特に有益である。このような方法及び装置を撮像中にも撮像完了後の後処理にも利用できればさらに有益である。
【発明の概要】
【０００９】
本発明は画像がＤ次元空間内の点のドメインに対応する複数の輝度から成るディジタル画像中の管状構造を配向する方法を提供する。この方法は画像ドメイン内に点を選択し、選択された点の近傍において画像の勾配を計算し、選択された点において基本構造テンソルを計算し、選択された点に関する構造テンソルを求め、構造テンソルの固有値を見つけ、固有値を分析することによって管状構造と整列する固有ベクトルを見つけるステップを含む。
【００１０】
本発明の他の特徴として、選択された点を中心とする近傍に亘って画像をガウス・カーネル導関数Ｇで畳込むことによって画像の勾配を測定する。
【００１１】
本発明の他の特徴として、ガウス・カーネルは画像のサイズよりもほぼ２桁小さい標準偏差σＧを有する。
【００１２】
本発明の他の特徴として、画像勾配にその転置行列を乗算することによって基本構造テンソルを画定することができる。
【００１３】
本発明の他の特徴として、検出中の対象物のサイズに対応する標準偏差σＴのガウス・カーネルで基本構造テンソルを畳込むことによって構造テンソルを求めることができる。
【００１４】
本発明の他の特徴として、Householder QL分解を行うことによって固有値を見出すことができる。
【００１５】
本発明の他の特徴として、最小固有値に対応する固有ベクトルが管状構造と整列する。
【００１６】
本発明の他の特徴として、選択された点の近傍を、前記近傍の２つの最小固有値を比較することによって、管状構造と比較することができ、前記近傍が局部的に管状構造と似ておれば前記２つの最小値が非近似である。
【００１７】
本発明の他の特徴として、管状構造と整列する固有ベクトルによって画定される軸を中心にカートホイール投影を計算することができる。
【実施例】
【００１８】
本発明は脈管構造の軸心を自動的に検出し、軸心を利用してこのような構造の可視化を支援するＣＡＤ方法に係わる。
【００１９】
本発明の図示実施例を以下に説明する。煩雑になるのを避けるため、本明細書では実施の構成要件は一部その説明を省略する。尚、このような具体的な実施態様を開発するに当たって、開発者が特定の目標、例えば、システム及び業務に関連する固有の制約とのコンプライアンスを達成するには、幾つもの実施に特異な判断を下さねばならないことは言うまでもない。また、このような開発努力は複雑を極め、多大の時間を要するが、本発明の開示内容を活用しようとする当業者にとっては当然果すべき努力である。
【００２０】
本発明の実施に当たっては多様な形態が可能であるが、ここでは添付の図面を参照しながら特定の実施態様を詳細に説明する。但し、特定の実施態様に関する以下の説明は本発明をこれら特定の実施態様に制限することを意図するものではなく、添付の請求項に定義されている本発明の思想と範囲を逸脱しない限り、すべての変更案、等価案及び代案は本発明に包含される。
【００２１】
ここに使用する語“画像”は（例えば、２Ｄ画像ならピクセル、３Ｄ画像ならボクセルのような）不連続な画素から構成された多次元データを意味する。画像としては、例えば、コンピュータ断層画像、磁気共鳴断層画像、超音波画像のほか、当業者には公知の医療用画像が含まれる。画像は非医療撮像システム、例えば、遠隔感知システムや電子顕微鏡からも得られる。画像はＲ３からＲへの写像と考えることができるが、本発明の方法はこのような画像に制限されるものではなく、例えば、２Ｄ画面や３Ｄボリュームのような任意の次元の画像にも適用できる。本発明はコンピュータ・システム、例えば、本発明のアルゴリズムを実行するPentium(商標）級のパソコンにおいて、実施することが好ましい。コンピュータはプロセッサ、メモリ、及び種々の入/出力手段を含む。胸部の体積部を表す一連のディジタル画像をコンピュータに入力する。ここで使用する語“ディジタル”及び“ディジタライズド”は必要に応じてディジタル収集システムを介して、またはアナログ画像からの変換によって得られるディジタルまたはディジタル化形式の画像または体積部を指す。
【００２２】
脈管構造は医療用画像において広く見出される管状対象物の例である。医療用画像における管状対象物の他の例としては、血管、気管支、腸管、神経及び特殊な骨が含まれる。医療用画像における管状対象物の表現及び分析は医療担当者が患者の複雑な生体構造上の特徴の理解を助け、治療を容易にする。脈管構造の３Ｄ画像、例えば、ＣＴスキャンの精査に際して、医師は軸スライスを利用することによって（例えば、小結節または塞栓のような）異常構造を検出することができるが、構造の形状をより正確に分析するには、さらなる追加の画像があれば有用である。その可能性の１つとして、軸線を中心に投影面をカートホイールさせるカートホイール投影がある。これのよって医師は構造が丸いかどうかを容易に判断することができる。もう１つの可能性として、管状構造の軸線と直交する投影面の分析がある。これらの技術を駆使するには入力としての軸線が前提となる。この軸線は管状構造の軸線であることが好ましい。初期設定で恣意的な軸線を選ぶと、結果として好ましくない可視化となることがある。
【００２３】
典型的な分析状況において、医師は球状構造を探すために大きい画像、例えば、肺のＣＴ画像を精査する。３つの次元すべてにおいて画像が大きい。通常、医師は軸方向の画像、即ち、体積部のＸ-Ｙスライスだけを１枚ずつ往復して観察する。スライスは典型として５１２×５１２ピクセルから成るが、医師が観察する構造はせいぜい数ピクセルの広さである。従って、医師は画像の大部分を無視しがちであるが、時には医師が、より精密に構造を観察したいと考える場合がある。より精密に観察する場合、医師がＸ-Ｙスライスだけではなく、完全な３Ｄ情報を得たいと考えることが好ましい。この場合、３Ｄツールと共にカートホイール投影が有用であり、この３Ｄツールは独自に、またはカートホイール投影と併用することができる。
【００２４】
カートホイール投影を利用することによって球状構造と管状構造とを識別することができる。問題は、Ｘ-Ｙスライスだけを利用する場合、双方の構造はいずれも円形に見えることである。しかし、完全に３Ｄ回転させることによって、医師は脈管と小結節とを識別することができる。カートホイール投影は確かにこのような識別を可能にするが、原形において、カートホイールの回転軸は固定されており（通常はＸまたはＹ軸）、脈管軸線に対してはランダムである。従って、脈管と小結節とを識別するため、医師は脈管軸線がカートホイール面に含まれる時点での回転角を求めねばならない。ここに開示する方法は脈管軸線の検出を助け、この軸線を可視化技術に活用する方法である。
【００２５】
本発明の好ましい実施例では、局部構造テンソルを利用することにより、推測される局部主軸から始めて、脈管の主軸を位置検出することができる。主軸の位置はユーザーがマウスのような入力装置を使用してスライス中の点を選択することによって初期推測することができる。図１に示す構造テンソルの好ましい計算方法は基準軸に沿った画像偏導関数で形成される３Ｄベクトルである画像勾配を計算することから始まる：
2007505711.tif 000003
【００２６】
実際には、画像が不連続な点においてサンプリングされるだけであり、ノイズは免れない。好ましい実施例としては、画像をガウス導関数で畳込むことによって画像近傍の各点における勾配を推算することができる（ステップ１０１）、即ち：
2007505711.tif 000004
但し、Ｇは標準偏差σＧの不連続な、基準化されたＤ-次元ガウス・カーネルであり、
2007505711.tif 000005
演算子 * は畳込みである。通常、標準偏差は画像の全体サイズと比較して、例えば最大限でも３ボクセルと遥かに小さく、または画像の１つの次元におけるボクセル数より２桁も小さく、任意の点を中心とする小さい体積部に勾配計算を限定してノイズの影響を排除するように選択される。
【００２７】
構造テンソルはステップ１０３における畳込みによって得られる３×３行列であり、探索中の対象物に相当するサイズの空間フィルタを有する転置行列との外積である：
2007505711.tif 000006
【００２８】
上記式中のσは極めて大きく、対象物のサイズにほぼ比例する。他の畳込みカーネルを使用することもできるが、ガウス・カーネルが好ましい。構造テンソルの３つの固有値はステップ１０４において公知の適当な技術で計算することができる。その１つがHouseholder  QL分解である。
【００２９】
構造テンソルの固有値はすべて互いに直交関係にあり、ステップ１０５において、最小の固有値に対応する固有ベクトルが脈管構造の軸である可能性が最も高い。即ち、脈管の内側または近傍に起点を選択すれば、固有値が最小の固有ベクトルが脈管と整列する可能性が最も高い。脈管が湾曲しているか、または分岐している場合でも、最小固有値と連携する固有ベクトルを求めれば、少なくとも部分的に脈管と整列することになる。例えば、通常なら全く異なる筈の２つの最小固有値を比較することにより、画像が局部的に脈管構造に似ているか否かをチェックすることができる。
【００３０】
ステップ１０６において、クリック‐ポイントとこの固有ベクトルによって確定されるこの軸を中心にカートホイール投影を行う。図２は問題の対象物のカートホイール投影を示す。カートホイール投影が行われるそれぞれの軸方向画像スライスを、ＣＴ装置による肺２００の画像走査から得、問題の対象物、例えば、疑わしい部位構造にセンタリングする。図２に示すように、回転投影面２０１を問題の点２０２を中心に１８０°回転させ、（以下、カートホイール投影画像スライスと呼称する）軸方向画像スライスの一連の２Ｄ投影、例えば、カートホイール投影画像スライス２０３及び２０４を得る。回転角を所定のインターバル、例えば、５°に設定することができ、この場合、所与の軸方向画像スライスにおける入力される問題対象物毎に３６枚の個別カートホイール投影画像スライスが形成される。同様に、インターバルを１°に設定すると、入力される問題対象物毎に１８０枚の個別カートホイール投影画像スライスが形成される。
【００３１】
カートホイール投影を利用することによって球状構造と管状構造とを識別することができる。軸と直交する投影を利用することによって脈管と直交するスライスをも得ることができる。回転軸が脈管と整列すれば、最初の角度から最終の角度まですべての投影面が脈管の軸と直交するから、構造が脈管であることはカートホイール投影から明白である。
【００３２】
図３及び４はこのようなインターフェースの１例を示す。図３に十字線で示す画像中の局部に合わせてユーザーがクリックする。もしマウスが指摘する構造が脈管なら、その軸に沿ってカートホイール投影が行われ、所与の方向に自動的に表示される。図４は水平方向に表示された脈管を示す。
【００３３】
もし構造が小結節であれば、すべてのスライスがディスクの形状を呈するから、回転軸が多少ともランダムになる可能性がある。但し、この状況において、医師は通常、非ランダムな軸、典型的にはＸまたはＹ軸の使用を望む。従って、構造テンソルの計算に際して、もしすべての固有値が互いに近似であることが検知されたら、テンソル軸の使用から基準カートホイール軸の使用へシフトすればよい。
【００３４】
本発明は種々の形態のハードウェア、ソフトウェア、ファームウェア、特殊目的プロセス、またはこれらの組合せとして実施することができる。１例として、本発明をコンピュータ可読プログラム記憶装置に記憶されるアプリケーション・プログラムとして機能するソフトウェアの形に実施することができる。アプリケーション・プログラムは適当なアーキテクチャーを含むマシーンにアップロードされ、このマシーンによって実行される。
【００３５】
図４に示すように、本発明を実施するためのコンピュータ・システム４０１は、特に、中央処理装置（ＣＰＵ）４０２、メモリ４０３及び入/出力（Ｉ/Ｏ）インターフェース４０４を含む。コンピュータ・システム４０１は一般にＩ/Ｏインターフェース４０４を介してディスプレイ４０５及びマウスやキーボードのような種々の入力装置４０６と接続する。支援回路はキャッシュ、電源、クロック回路のような回路と通信用バスを含むことができる。メモリ４０３はランダムアクセスメモリ（ＲＡＭ）、読み出し専用メモリ（ＲＯＭ）、ディスク・ドライブ、テープ・ドライブ、等、またはこれらの組合せを含むことができる。本発明はメモリ４０３に記憶され、信号源４０８からの信号を処理するためＣＰＵによって実行されるルーチン４０７として実施することができる。従って、コンピュータ・システム４０１は本発明のルーチン４０７を実行する際に専用コンピュータ・システムとなる汎用コンピュータ・システムである。
【００３６】
コンピュータ・システム４０１はオペレーション・システム及びマイクロ命令コードをも含む。以上に述べた種々のプロセス及び機能はオペレーション・システムを介して実行されるマイクロ命令コードの一部またはアプリケーション・プログラム（または両者の組合せ）である。さらに、コンピュータ・プラットホームには、補足のデータ記憶装置やプリンターなどのような他の周辺装置を接続することができる。
【００３７】
添付図面に示すシステム構成成分及び方法ステップの幾つかはソフトウェアにおいて実行されるから、システム構成成分（または方法ステップ）間の実際の接続は本発明をプログラムする態様に応じて異なる。当業者ならば、本発明の教示内容に照らして、以上に述べた態様及び類似の態様で本発明を実施または構成できるであろう。
【００３８】
本発明には変更を加えることが可能であり、上記教示内容を利用し得る当業者ならば、種々の、但し等価の態様で実施することができるから、上記の具体的な実施例は飽くまでも本発明の説明がその目的である。また、本発明は図示の構成または設計のうち、下記の請求項に記載されている事項以外の細部に制限されるものではない。従って、上記特定の実施例に変更を加えることは可能であり、このような変更がすべて本発明の範囲及び思想に包含されることは言うまでもない。従って、ここに請求する特許権は下記請求項によって定義される。
【図面の簡単な説明】
【００３９】
【図１】本発明の好ましい方法のフローチャートである。
【図２】問題オブジェクトにおける画像スライスのカートホイール投影を示す図である。
【図３】分析したい構造をユーザーが選択する態様を示す。
【図４】本発明の方法によって発見される管状構造を示す。
【図５】本発明の方法の好ましい実施態様を実行するコンピュータの１例を示す。
【図１】
2007505711.tif 000007
【図２】
2007505711.tif 000008
【図３】
2007505711.tif 000009
【図４】
2007505711.tif 000010
【図５】
2007505711.tif 000011
【国際調査報告】
2007505711.tif 000012
2007505711.tif 000013
2007505711.tif 000014