(57)【要約】【課題】トモシンセシスにおいて、点広がり関数（ＰＳＦ）に基づいて構造（すなわち領域）を分類する能力を提供して、最終的な再構成画像でのアーティファクト及び／又はボケを低減する。【解決手段】データ集合において領域を分類するために点広がり関数型規則を用いる。他の観点では、データ集合において領域を分類するために点広がり関数型規則を用いることをコンピュータに命令するように構成されているプログラムを組み入れたコンピュータ読み取り可能な媒体が提供される。さらにもう一つの観点では、システムが、Ｘ線源と、該線源から放出されるＸ線を受光するように配置されているＸ線検出器と、線源及び検出器に動作に関して結合されているコンピュータとを含んでいる。コンピュータは、データ集合において構造を分類するために点広がり関数を用いるように構成されている。【選択図】        図１
2008043756.tif 000002
【特許請求の範囲】
【請求項１】
Ｘ線源（１０２）と、
該線源（１０２）から放出されるＸ線を受光するように配置されているＸ線検出器（１０８）と、
前記線源（１０２）及び前記検出器（１０８）に動作に関して結合されており、データ集合において領域を分類するために点広がり関数型規則を用いるように構成されているコンピュータと、
を備えたシステム（１００）。
【請求項２】
前記コンピュータは、トモシンセシス・データ集合において領域を分類するために前記点広がり関数型規則を用いるようにさらに構成されている、請求項１に記載のシステム（１００）。
【請求項３】
前記分類は、背景クラス、合焦点クラス、非合焦点クラス及び低周波数クラスを含むクラスへの分類である、請求項１に記載のシステム（１００）。
【請求項４】
前記分類は、背景クラス、合焦点クラス、非合焦点クラス及び低周波数クラスから成るクラスへの分類である、請求項１に記載のシステム（１００）。
【請求項５】
前記分類は、背景クラス、合焦点クラス、非合焦点クラス及び低周波数クラスを含むクラスへの分類である、請求項２に記載のシステム（１００）。
【請求項６】
前記コンピュータは、各々のピクセルに一つずつの複数の強度プロファイルを生成するようにさらに構成されている、請求項１に記載のシステム（１００）。
【請求項７】
前記コンピュータは、再構成画像に前記異なる分類の視覚的指標を設けるようにさらに構成されている、請求項１に記載のシステム（１００）。
【請求項８】
前記分類は、背景クラス、合焦点クラス、非合焦点クラス及び低周波数クラスを含むクラスへの分類である、請求項７に記載のシステム（１００）。
【請求項９】
前記コンピュータは、各々のピクセルに一つずつの複数の強度プロファイルを生成するようにさらに構成されている、請求項８に記載のシステム（１００）。
【請求項１０】
前記コンピュータは、トモシンセシス・データ集合において構造を分類するために前記点広がり関数を用いるようにさらに構成されている、請求項９に記載のシステム（１００）。
【請求項１１】
走査対象を設けるステップと、
データ集合を取得するために前記対象を走査するステップと、
前記データ集合において領域を分類するために点広がり関数型規則を用いるステップと、
を備えた方法。
【請求項１２】
前記データ集合は医用画像データ集合である、請求項１１に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は一般的には、Ｘ線の方法及び装置に関し、さらに具体的には、データ集合における構造の分類を提供する方法及び装置に関する。
【背景技術】
【０００２】
ディジタル・トモシンセシスは、静止型ディジタル検出器によって角度限定型のＸ線投影撮像から取得される対象の三次元（３Ｄ）再構成に広く用いられている。ディジタル・トモシンセシスは、１９３０年代以来公知である従来の幾何学的断層写真法の改良技術である。
【特許文献１】米国特許第５３７５１７５号
【特許文献２】米国特許第６０５８３２２号
【発明の開示】
【発明が解決しようとする課題】
【０００３】
幾何学的断層写真法と同様に、トモシンセシスでも関心平面外に位置する対象の残留ボケが起こる。重なった解剖学的構造によるこの断層写真のボケのため、関心平面内の細部が不明瞭になり、スライスのコントラスト強調が制限される。重なってボケた解剖学的構造を除去すると面内構造のコントラストが改善し、この除去は、画像のダイナミック・レンジを関心のある区画のダイナミック・レンジに制限して、かかる区画の幾つかの関心対象の周波数内容に類似した周波数内容を有し得る残留構造を除去することにより行なわれる。基本的レベルでは、トモシンセシス・データ集合の点広がり関数（ＰＳＦ）が撮像容積におけるボケの広がりの特徴を決定しており、この関数は本質的にシフト・バリアント（位置依存型）である。しかしながら、完全なボケの除去は瑣末なタスクではない。このタスクは計算集約的であり、またＰＳＦの範囲のため、ボケを排除することは容易ではない。
【課題を解決するための手段】
【０００４】
一観点では、方法が、データ集合において領域を分類するために点広がり関数型規則を用いるステップを含んでいる。
【０００５】
さらにもう一つの観点では、データ集合において領域を分類するために点広がり関数型規則を用いることをコンピュータに命令するように構成されているプログラムを組み入れたコンピュータ読み取り可能な媒体が提供される。
【０００６】
もう一つの観点では、システムが、Ｘ線源と、該線源から放出されるＸ線を受光するように配置されているＸ線検出器と、線源及び検出器に動作に関して結合されているコンピュータとを含んでいる。コンピュータは、データ集合において構造を分類するために点広がり関数を用いるように構成されている。
【０００７】
さらにもう一つの観点では、方法が、走査対象を設けるステップと、データ集合を取得するために対象を走査するステップと、データ集合において領域を分類するために点広がり関数型規則を用いるステップとを含んでいる。
【０００８】
さらにもう一つの観点では、方法が、データ集合にアクセスするステップと、データ集合において領域を分類するために点広がり関数型規則を用いるステップとを含んでいる。
【発明を実施するための最良の形態】
【０００９】
本書では、例えば限定しないがＸ線トモシンセシス・システムのようなイメージング・システムに有用な方法及び装置について記載する。これらの装置及び方法は図面に関して説明され、図面では類似の参照符号が全図面において同じ構成要素を指す。かかる図面は制限的ではなく説明的であることを意図しており、本発明の装置及び方法の実施形態の一例の説明を容易にするために本書に含まれている。Ｘ線トモシンセシス・システムの設定で記載されるが、本発明の利益はＸ線源を有する全てのシステムに齎されると思量される。この開示の一つの目的は、トモシンセシスのデータ集合を、ＰＳＦの性質に基づいて多数のカテゴリに分解する分類の枠組みを提供することにある。この分類は、本開示と同日に同時出願された米国特許出願第１１／４６４１０３号に記載されている処理経路の選択を誘導することができる。
【００１０】
図１は、例示的なＸ線トモシンセシス・イメージング・システム１００を示している。イメージング・システム１００は、被検構造１０６にＸ線フォトンを照射するＸ線源１０２を含んでいる。例として、Ｘ線源１０２はＸ線管であってよく、被検構造１０６は患者、試験用ファントム、及び／又は他の被験無生物体であってよい。
【００１１】
Ｘ線イメージング・システム１００はまた、処理回路に結合された検出器１０８を含んでいる。処理回路（例えばマイクロコントローラ、マイクロプロセッサ、又はカスタムＡＳＩＣ等）は、メモリ及び表示装置に結合されている。メモリ（例えばフレキシブル・ディスク、又はネットワーク若しくはインターネットのような他のディジタル・ソース等のコンピュータ読み取り可能な媒体から命令及び／又はデータを読み取るフレキシブル・ディスク・ドライブ、ＣＤ−ＲＯＭドライブ、ＤＶＤドライブ、光磁気ディスク（ＭＯＤ）装置、又はイーサネット装置等のネットワーク接続装置を含めたその他任意のディジタル装置、並びに開発中のディジタル手段等の１又は複数を含む）は、撮像データを記憶する。
【００１２】
メモリはまた、本書に記載される作用を具現化するために処理回路によって実行される命令を含んだコンピュータ・プログラムを記憶することができる。処理回路は、装置に表示するための画像を形成する。本書にさらに詳細に記載されるように、画像は様々な構造（例えば軟組織、骨）を表わしていてよい。検出器１０８は、例えばフラット・パネル型固体画像検出器であってよいが、メモリにディジタル形態で記憶された従来のフィルム画像を処理することもできる。一実施形態では、処理回路はファームウェア（図示されていない）に記憶されている命令を実行する。一般的には、プロセッサは後述する工程を実行するようにプログラムされる。
【００１３】
言うまでもなく、本書に記載される方法は、システム１００での実施に制限されてはおらず、イメージング・システムのその他多くの形式及び変形と関連して利用することができる。一実施形態では、処理回路は、本書に記載される作用を果たすようにプログラムされているコンピュータであり、従って、本書で用いられるコンピュータとの用語は当技術分野でコンピュータと呼ばれている集積回路のみに限らず、コンピュータ、プロセッサ、マイクロコントローラ、マイクロコンピュータ、プログラマブル論理コントローラ、特定応用向け集積回路、及び他のプログラム可能な回路を広範に指している。本書に記載される方法は、人間の患者の設定において記載されるが、本発明の利益は、小動物研究に典型的に用いられるシステムのような人間以外を対象とするイメージング・システムにも齎されると思量される。
【００１４】
本書で用いる場合には、単数形で記載されており単数不定冠詞を冠した要素またはステップとの用語は、排除を明記していない限りかかる要素又はステップを複数備えることを排除しないものと理解されたい。さらに、本発明の「一実施形態」に対する参照は、所載の特徴を同様に組み入れている他の実施形態の存在を排除すると解釈されることを意図していない。
【００１５】
また、本書で用いられる「画像を再構成する」との表現は、画像を表わすデータが生成されるが可視画像は形成されないような本発明の実施形態を排除するものではない。従って、本書で用いられる「画像」との用語は、可視画像及び可視画像を表わすデータの両方を広く指す。但し、多くの実施形態は少なくとも１枚の可視画像を形成する（か又は形成するように構成されている）。
【００１６】
本書に記載される方法は医療の設定において記載されるが、産業の設定、又は例えば限定しないが空港若しくは他の運輸拠点での手荷物走査システムのような運輸の設定で典型的に用いられるシステム等のような非医用イメージング・システムにおいても、本発明の利益が齎されると思量される。
【００１７】
再構成された複数のスライスにわたって生成される強度プロファイルに基づいて、トモシンセシス撮像容積を合焦点（in-focus）構造、非合焦点（out-of-focus）構造、背景構造及び低周波数構造に分類する新規の方法体系が開発された。
【００１８】
再構成のために、線形ディジタル・トモシンセシス幾何学的構成及び逆投影アルゴリズムについて検討した。但し、本発明の利益は、ＰＳＦがシフト・バリアントであるような全ての場合に齎されると思量され、円形の幾何学的構造にも、逆投影以外の再構成手法にも同等に適用可能である。撮像容積内の全ての点における三次元ＰＳＦは、その形状及び強度に関して決定されて特徴評価されている（Avinash等、２００６年３月２日、「CHARACTERIZATION OF POINT SPREAD FUNCTION IN LINEAR DIGITAL TOMOSYNTHESIS: A SIMULATION STUDY」、ＳＰＩＥ、第６１４２巻を参照されたい）。本分類手法では、ＰＳＦの強度プロファイルが重要な役割を果たす。これらの強度加重値は、有限の検出器分解能を考慮して、ＰＳＦの全てのスライスにおいて線形補間される。正確に行なわれれば、ＰＳＦの全てのスライスの加重値の積分は、（ＰＳＦを推定している）カレントのピクセルを通る投影射線の数（Ｎ）に等しくなる。撮像容積内の様々な位置でのシフト・バリアント型ＰＳＦを図２に示す。
【００１９】
容積内の全ての点について、ＰＳＦによって画定される線分に含まれる強度の積分値を各々のスライスにおいてプロットした。この方法を図３に図示しており、分類に用いられる強度プロファイルに対応するプロットも併せて示す。５枚のスライスを有するトモシンセシス再構成データ集合を考える。スライス４の任意の点について、ＰＳＦは図３（Ａ）に示すようなものである。同じ形状及び範囲を有する対応する強度マスクを図３（Ｂ）に示す。これらのマスクの単純乗算の結果を図３（Ｃ）に示す。Ｎ（スライス）が各々のスライスにおいてＰＳＦの範囲内にあるピクセル数を表わし、Ｗ（ピクセル，スライス）が各々のスライスでのピクセルの加重値を表わし、Ｉ（ピクセル，スライス）が当該ピクセルの強度を表わすとすると、この方法の数学的な公式化を下記のように与えることができる。
スライス１：
【００２０】
【数１】
2008043756.tif 000003
スライス２：
【００２１】
【数２】
2008043756.tif 000004
スライス５：
【００２２】
【数３】
2008043756.tif 000005
全てのスライスの強度を合計した後に、結果（Ａ、Ｂ、Ｃ、Ｄ、Ｅ）が図３（Ｄ）に示すようにプロットされる。プロット内の基準の点線は、当該ピクセルのカレントの強度に対応する。撮像容積におけるこれらのプロファイルについては、全ての構造が、別個であるが一貫した特徴を示した。全てのピクセルについて強度プロファイルの性質を解析することにより、これらの強度プロファイルを合焦点細部、非合焦点細部、背景細部及び低周波数細部に分類することができる。言うまでもなく、他の分類体系も存在する。例えば、ピクセル（及び関連するピクセルの領域）を、様々な順列で、合焦点高周波数及び合焦点低周波数等に分類することができる。
【００２３】
背景ピクセル及び低周波数ピクセルでは、強度プロットはそれぞれ比較的低強度及び比較的高強度を有する平坦な曲線を示した。合焦点ピクセルでは、プロットの最大強度は当該ピクセルのカレントの強度（スライス０）と一致しており、非合焦点ピクセルでは、プロットの最大強度（スライス２）は当該ピクセルのカレントの強度（スライス０）よりも大きかった。１回目の工程の後に、この分類工程を整合性検査のために全てのピクセルについて反復的に繰り返した。様々なカテゴリについての一連の強度プロットを図４に示す。
【００２４】
本書に記載される方法及び装置を、人間の観察者が分類した領域を「判断基準（gold standard）」として用いて検証した。２個の球体を有するシミュレーション・ファントムが図５に示すように画定され、ファントムの各領域を強度プロットを用いて分類した。次に、算出された分類を判断基準と比較して、この比較に基づいて本方法について感度及び特異性が画定された。
【００２５】
コンピュータでシミュレートされたファントムについて全てのピクセルをそれぞれのカテゴリに分類した後に、再構成されたデータ集合を視覚的表示のために色分けした。様々な濃淡レベル（グレイ・レベル）を用いて合焦点ピクセル、非合焦点ピクセル、背景ピクセル及び低周波数ピクセルを標識した。結果を図６に示す。この方法を、本開示と同時に出願された米国特許出願第１１／ＸＸＸ，ＸＸＸ号（出願人控番号第２９７４９４号）の添付開示に記載されているような適応型画像処理の入力として用いることができる。
【００２６】
一つの技術的効果は、本書に記載される方法及び装置が、ＰＳＦに基づいて構造（すなわち領域）を分類する能力を提供することであり、すると所望に応じて各々の異なるクラスについて異なるようにデータを選択的に処理することができる。構造／領域の様々なクラスを差異化するこの改良された能力によって、最終的な再構成画像でのアーティファクト及び／又はボケの低減が導かれる。
【００２７】
以上、実施形態の例について詳細に説明した。これらのアセンブリ及び方法は、本書に記載される特定の実施形態に限定されてはおらず、各々のアセンブリ及び／又は方法の構成要素を、本書に記載される他の構成要素と独立且つ別個に利用してよい。
【００２８】
本発明を様々な特定の実施形態について記載したが、当業者であれば、特許請求の範囲の要旨及び範囲内で改変を施して本発明を実施し得ることを認められよう。また、図面の符号に対応する特許請求の範囲中の符号は、単に本願発明の理解をより容易にするために用いられているものであり、本願発明の範囲を狭める意図で用いられたものではない。そして、本願の特許請求の範囲に記載した事項は、明細書に組み込まれ、明細書の記載事項の一部となる。
【図面の簡単な説明】
【００２９】
【図１】例示的なＸ線イメージング・システムを示す図である。
【図２ａ】再構成された容積内の様々な位置での三次元ＰＳＦを示す図である。
【図２ｂ】再構成された容積内の様々な位置での三次元ＰＳＦを示す図である。
【図２ｃ】再構成された容積内の様々な位置での三次元ＰＳＦを示す図である。
【図３ａ】複数のスライスにわたって強度プロファイルを用いたＰＳＦに基づく分類を示す図である。
【図３ｂ】複数のスライスにわたって強度プロファイルを用いたＰＳＦに基づく分類を示す図である。
【図３ｃ】複数のスライスにわたって強度プロファイルを用いたＰＳＦに基づく分類を示す図である。
【図３ｄ】複数のスライスにわたって強度プロファイルを用いたＰＳＦに基づく分類を示す図である。
【図４ａ】背景ピクセルの強度プロファイルを示す図である。
【図４ｂ】合焦点ピクセルの強度プロファイルを示す図である。
【図４ｃ】非合焦点ピクセルの強度プロファイルを示す図である。
【図４ｄ】低周波数ピクセルの強度プロファイルを示す図である。
【図５】コンピュータ・シミュレーションによる球体ファントムの模式的モデルの図である。
【図６ａ】逆投影再構成アルゴリズムを用いて再構成されたトモシンセシス・スライスを示す図である。
【図６ｂ】シミュレーション・ファントムのＰＳＦに基づく分類を示す図である。
【符号の説明】
【００３０】
１００  イメージング・システム
１０２  Ｘ線源
１０６  被検体
１０８  検出器
【図１】
2008043756.tif 000006
【図２ａ】
2008043756.tif 000007
【図２ｂ】
2008043756.tif 000008
【図２ｃ】
2008043756.tif 000009
【図３ａ】
2008043756.tif 000010
【図３ｂ】
2008043756.tif 000011
【図３ｃ】
2008043756.tif 000012
【図３ｄ】
2008043756.tif 000013
【図４ａ】
2008043756.tif 000014
【図４ｂ】
2008043756.tif 000015
【図４ｃ】
2008043756.tif 000016
【図４ｄ】
2008043756.tif 000017
【図５】
2008043756.tif 000018
【図６ａ】
2008043756.tif 000019
【図６ｂ】
2008043756.tif 000020