(57)【要約】【課題】  診断対象となる診断部位を含む診断用画像を精度良く、かつ高速に判定して、出力することができる画像処理装置及びその方法、プログラムを提供する。【解決手段】  放射線画像から、前記被検体のアキシャル画像以外の第１断面画像を作成する。第１断面画像中の陰影を検出する。検出した陰影を含む第２断面画像を、放射線画像から作成する。第２断面画像の特徴量を算出する。算出された特徴量に基づいて、検出した陰影が診断対象部位であるか否かを判定する。その判定結果に基づいて、検出した陰影を含む診断用画像を出力する。【選択図】  図１
2005296348.tif 000002
【特許請求の範囲】
【請求項１】
放射線を利用して得られる被検体の放射線画像を処理して、診断用画像を出力する画像処理装置であって、
前記放射線画像から、前記被検体のアキシャル画像以外の第１断面画像を作成する第１作成手段と、
前記第１断面画像中の陰影を検出する検出手段と、
前記検出手段で検出した陰影を含む第２断面画像を、前記放射線画像から作成する第２作成手段と、
前記第２断面画像の特徴量を算出する算出手段と、
前記算出手段で算出された特徴量に基づいて、前記検出手段で検出した陰影が診断対象部位であるか否かを判定する判定手段と、
前記判定手段の判定結果に基づいて、前記前記検出手段で検出した陰影を含む診断用画像を出力する出力手段と
を備えることを特徴とする画像処理装置。
【請求項２】
前記第１断面画像は、前記被検体のコロナル画像、サジタル画像及びそれらのＲＡＹＳＵＭ画像のいずれかである
ことを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記第２断面画像は、前記被検体のアキシャル画像である
ことを特徴とする請求項１に記載の画像処理装置。
【請求項４】
前記第１断面画像が前記被検体のコロナル画像である場合、前記第２断面画像は、前記被検体のアキシャル画像あるいはサジタル画像である
ことを特徴とする請求項１に記載の画像処理装置。
【請求項５】
前記放射線画像は、等方性画素で構成されている
ことを特徴とする請求項１に記載の画像処理装置。
【請求項６】
前記算出手段は、前記第２断面画像中の陰影の円形度を特徴量として算出する
ことを特徴とする請求項１に記載の画像処理装置。
【請求項７】
放射線を被検体に照射して、該被検体の放射線画像を撮像して出力する撮像手段を更に備える
ことを特徴とする請求項１に記載の画像処理装置。
【請求項８】
放射線を利用して得られる被検体の放射線画像を処理して、診断用画像を出力する画像処理方法であって、
前記放射線画像から、前記被検体のアキシャル画像以外の第１断面画像を作成する第１作成工程と、
前記第１断面画像中の陰影を検出する検出工程と、
前記検出工程で検出した陰影を含む第２断面画像を、前記放射線画像から作成する第２作成工程と、
前記第２断面画像の特徴量を算出する算出工程と、
前記算出工程で算出された特徴量に基づいて、前記検出工程で検出した陰影が診断対象部位であるか否かを判定する判定工程と、
前記判定工程の判定結果に基づいて、前記前記検出工程で検出した陰影を含む診断用画像を出力する出力工程と
を備えることを特徴とする画像処理方法。
【請求項９】
放射線を利用して得られる被検体の放射線画像を処理して、診断用画像を出力する画像処理を実現するプログラムであって、
前記放射線画像から、前記被検体のアキシャル画像以外の第１断面画像を作成する第１作成工程のプログラムコードと、
前記第１断面画像中の陰影を検出する検出工程のプログラムコードと、
前記検出工程で検出した陰影を含む第２断面画像を、前記放射線画像から作成する第２作成工程のプログラムコードと、
前記第２断面画像の特徴量を算出する算出工程のプログラムコードと、
前記算出工程で算出された特徴量に基づいて、前記検出工程で検出した陰影が診断対象部位であるか否かを判定する判定工程のプログラムコードと、
前記判定工程の判定結果に基づいて、前記前記検出工程で検出した陰影を含む診断用画像を出力する出力工程のプログラムコードと
を備えることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、放射線を利用して得られる被検体の放射線画像を処理して、診断用画像を出力する画像処理装置及びその方法、プログラムに関するものである。
【背景技術】
【０００２】
従来、被検体に対してＸ線を曝射し、該被検体を透過或いは被検体で散乱したＸ線をＸ線検出器で検出し、このＸ線検出出力（Ｘ線のフォトン数）に基づいて被検体の透視画像、断層像或いは三次元画像を撮像するＸ線ＣＴ装置が知られている。
【０００３】
このＸ線ＣＴ装置として、コーンビームＣＴ装置が開発されている。通常のＸ線ＣＴ装置では、Ｘ線ビームはＺ方向に薄く切り出されており、ファンビームと呼ばれるが、コーンビームＣＴ（ＣＢＣＴ）では、Ｚ方向にも広がったＸ線ビームを用い、このＸ線ビームはコーンビームと呼ばれる。
【０００４】
ところで、マルチスライスＣＴやＣＢＣＴでは、その撮影結果として、大量の画像が発生するので、その大量の画像から診断用に適した画像を選別して提示するための、医師に対する診断支援装置が重要になっている。特には、比較的、病理を見つけやすい胸部ＣＴ画像から石灰化（ノジュール）を見つけることが、ＣＴを検診目的に使用する場合の要求にかなっている。従来より、アキシャル画像中のノジュールの検出を行うことが良く行われている。
【０００５】
例えば、まず、２次元画像（アキシャル）情報から腫瘍か否かを判定すべき特定の領域（陰影）抽出し、当該領域の面積Ｓを算出し、この抽出した領域の重心を算出する。次に、抽出領域と同一面積を有する円の半径ｒ（以下、抽出領域の等価半径と称する）を次式により算出する。
【０００６】
ｒ＝（Ｓ／π）１／２                                            （１）
続いて、抽出領域の円形度を算出する。円形度は、算出した重心を中心として、算出した等価半径ｒの円に含まれる抽出領域の面積の割合により、次式のように定義する。
【０００７】
円形度＝円に含まれる抽出領域の面積／Ｓ                          （２）
そして、得られる等価半径ｒと円形度とで決まる特徴空間において、抽出領域の等価半径ｒと円形度が一定の範囲内に含まれているときは、抽出領域が腫瘍であるとみなす。
【０００８】
特許文献１では、更に、１枚の画像のみならず、３次元的に分布する多数枚の２次元画像情報に基づいて、それぞれの２次元画像内で領域を抽出し、抽出領域の等価半径と円形度を算出して、それによって各画像毎に腫瘍を検出する装置が提案されている。
【０００９】
特許文献２には、再構成して得た断層像の濃度データ値を、濃度データ値の低い所望の第１閾値を以て切り出し、注目画像として得る第１処理機能と、切り出した注目画像についてその形状等の特徴量から陰影か否かを判定する第１判定機能と、断層像の濃度データ値を、濃度データ値の高い所望の第２閾値を以て切り出し、注目画像として得る第２処理機能と、この切り出した注目画像についてその形状等の特徴量から陰影の候補となるか否かを判定する候補判定機能と、候補と判定された注目画像について、投影データからズーミングにより再構成して所定の倍率の拡大断層像を取得し、これを第１閾値を以て切り出し、注目画像として得る第２処理機能と、この第２処理機能により切り出した注目画像についてその形状等の特徴量から陰影か否かを判定する第２判定機能とを有する診断支援システムが提案されている。
【００１０】
特許文献３には、被検者の３次元領域を投影して得られた投影データから複数の断層像を算出し、この複数の断層像に基づいて当該被検者の診断を支援するための情報を出力デバイスを介して提供する医用診断支援システムにおいて、複数の断層像から病変部の疑いのある断層像を注目画像として決定する断層像解析手段と、注目画像に対応する３次元領域の複数の第１の２次元画像を残りの３次元領域の複数の第２の２次元画像から強調して出力デバイスに表示させる表示制御手段とを備える医用診断支援システムが提案されている。
【００１１】
特許文献４には、マルチスライスの第１の断層画像と、第１の断層画像と同じ被検者の同じ部位に関し、且つ第１の断層画像と撮影時期の異なるマルチスライスの第２の断層画像との各々から特定臓器の領域を抽出する領域抽出部と、第１の断層画像から抽出した特定臓器の領域と第２の断層画像から抽出した特定臓器の領域とのサイズ及び形態に関わるインデックスに基づいて、第１の断層画像と第２の断層画像との間で体軸方向に関する解剖学的な位置を揃えるスライスマッチング部とを備えるコンピュータ支援診断システムが提案されている。
【００１２】
特許文献５には、３次元デジタル画像データの診断支援のための計算機支援診断方法において、３次元デジタル画像データ内の３次元対象物を識別するステップを有し、１つの所定の３次元対象物に対して、所定の対象物に対してローカル回転平面を求めるステップを有し、ローカル回転平面は、所定対象物のセントロイド及びローカル回転軸にセンタリングされており、ローカル回転平面を、３６０°のうちの少なくとも部分的な角度にわたり回転させるステップを有し、前もって規定された回転のインクリメントにおいて所定対象物の観察画像を生成することによって、所定の対象物の複数の観察画像を得るステップを有する、３次元デジタル画像データの診断支援のための計算機支援診断方法が提案されている。
【００１３】
以上の特許文献１〜５は、処理対象画像が、全て、アキシャル画像をベースにした技術である。従来、アキシャル画像をベースに概略円形の陰影を検出する一次検出を行い、隣接する他のアキシャル画像での候補の連結性をもとに判定される。
【００１４】
アキシャル画像以外のコロナル画像を使用する診断支援技術としては、例えば、特許文献６として、複数のアキシャル断層像を利用する診断支援装置において、複数のアキシャル断層像の少なくとも１枚から関心領域を抽出する手段と、複数のアキシャル断層像のうちの抽出手段により抽出された関心領域内の画像を断面変換して複数のコロナル断層像を作成する手段を具備し、複数のコロナル断層像が診断に用いられることを特徴とする診断支援装置が提案されている。
【００１５】
この特許文献６によれば、アキシャル断層像を関心領域のみ断面変換してコロナル断層像を算出し、コロナル断層像により診断が行なわれることにより、アキシャル断層像を読影する場合に比べて医師が読影する画像枚数を少なくすることができる診断支援装置を提案している。
【特許文献１】特開平０７−２３６６３４号公報
【特許文献２】特開平０８−０６６３９２号公報
【特許文献３】特開平０８−１６６９９５号公報
【特許文献４】特開２００１−１３７２３０号公報
【特許文献５】特開２００２−０７８７０６号公報
【特許文献６】特開平７−２６５３００号公報
【発明の開示】
【発明が解決しようとする課題】
【００１６】
ところで、従来、アキシャル画像をベースにＣＡＤ処理が行われているのは、スパイラルＣＴ以前は体軸方向の連続性が保証されていなかったからで、断面変換を行うと体軸方向に接続不良による偽像が発生するからである。
【００１７】
近年のＭＤスパイラルＣＴにおいては、概略等方性の画像が構成可能になったが、スパイラルを行っている関係で完全ではない。これに対して、肺野などの器官を一回転でスキャンするＣＢＣＴにおいては、完全な等方性が保証されるので完全なコロナル画像が作成できる。
【００１８】
他方、肺野における血管枝の走行が上下方向に多いことを考えると、概略円形な陰影をアキシャル画像から検出することは血管と誤り易く利益がない。
【００１９】
本発明は上記の課題を解決するために、被検体画像から、診断対象となる診断部位を含む診断用画像を精度良く、かつ高速に判定して、出力することができる画像処理装置及びその方法、プログラムを提供することを目的とする。
【課題を解決するための手段】
【００２０】
上記の目的を達成するための本発明による画像処理装置は以下の構成を備える。即ち、  放射線を利用して得られる被検体の放射線画像を処理して、診断用画像を出力する画像処理装置であって、
前記放射線画像から、前記被検体のアキシャル画像以外の第１断面画像を作成する第１作成手段と、
前記第１断面画像中の陰影を検出する検出手段と、
前記検出手段で検出した陰影を含む第２断面画像を、前記放射線画像から作成する第２作成手段と、
前記第２断面画像の特徴量を算出する算出手段と、
前記算出手段で算出された特徴量に基づいて、前記検出手段で検出した陰影が診断対象部位であるか否かを判定する判定手段と、
前記判定手段の判定結果に基づいて、前記前記検出手段で検出した陰影を含む診断用画像を出力する出力手段と
を備える。
【００２１】
また、好ましくは、前記第１断面画像は、前記被検体のコロナル画像、サジタル画像及びそれらのＲＡＹＳＵＭ画像のいずれかである。
【００２２】
また、好ましくは、前記第２断面画像は、前記被検体のアキシャル画像である。
【００２３】
また、好ましくは、前記第１断面画像が前記被検体のコロナル画像である場合、前記第２断面画像は、前記被検体のアキシャル画像あるいはサジタル画像である。
【００２４】
また、好ましくは、前記放射線画像は、等方性画素で構成されている。
【００２５】
また、好ましくは、前記算出手段は、前記第２断面画像中の陰影の円形度を特徴量として算出する。
【００２６】
また、好ましくは、放射線を被検体に照射して、該被検体の放射線画像を撮像して出力する撮像手段を更に備える。
【００２７】
上記の目的を達成するための本発明による画像処理方法は以下の構成を備える。即ち、
放射線を利用して得られる被検体の放射線画像を処理して、診断用画像を出力する画像処理方法であって、
前記放射線画像から、前記被検体のアキシャル画像以外の第１断面画像を作成する第１作成工程と、
前記第１断面画像中の陰影を検出する検出工程と、
前記検出工程で検出した陰影を含む第２断面画像を、前記放射線画像から作成する第２作成工程と、
前記第２断面画像の特徴量を算出する算出工程と、
前記算出工程で算出された特徴量に基づいて、前記検出工程で検出した陰影が診断対象部位であるか否かを判定する判定工程と、
前記判定工程の判定結果に基づいて、前記前記検出工程で検出した陰影を含む診断用画像を出力する出力工程と
を備える。
【００２８】
上記の目的を達成するための本発明によるプログラムは以下の構成を備える。即ち、
放射線を利用して得られる被検体の放射線画像を処理して、診断用画像を出力する画像処理を実現するプログラムであって、
前記放射線画像から、前記被検体のアキシャル画像以外の第１断面画像を作成する第１作成工程のプログラムコードと、
前記第１断面画像中の陰影を検出する検出工程のプログラムコードと、
前記検出工程で検出した陰影を含む第２断面画像を、前記放射線画像から作成する第２作成工程のプログラムコードと、
前記第２断面画像の特徴量を算出する算出工程のプログラムコードと、
前記算出工程で算出された特徴量に基づいて、前記検出工程で検出した陰影が診断対象部位であるか否かを判定する判定工程のプログラムコードと、
前記判定工程の判定結果に基づいて、前記前記検出工程で検出した陰影を含む診断用画像を出力する出力工程のプログラムコードと
を備える。
【発明の効果】
【００２９】
本発明によれば、診断対象となる診断部位を含む診断用画像を精度良く、かつ高速に判定して、出力することができる画像処理装置及びその方法、プログラムを提供できる。
【発明を実施するための最良の形態】
【００３０】
以下、本発明の実施の形態について図面を用いて詳細に説明する。
【００３１】
本発明は、例えば、Ｘ線等の放射線を利用して画像撮影を行うＸ線ＣＴ装置等の放射線一般を使用して被検体内の放射線特性分布を画像化する放射線撮像技術に関するものである。特には、コーンビームＣＴ装置において、コンピュータ診断支援機能を有するもの、あるいはコーンビームＣＴ画像からのコンピュータ支援診断処理機能を有する画像処理技術に関するものである。
【００３２】
図１は本発明の画像処理システムの全体像を説明するための概念図である。
【００３３】
例えば、ＣＢＣＴ撮影装置（不図示）から得られる被写体の３６０度（あるいは１８０度＋ファン角）からの複数の２次元投影データ１に、逆ラドン変換２を施すことにより、被検体の３次元画像である３次元ボクセル３を再構成する。
【００３４】
尚、図１では、ＣＢＣＴ装置からの画像を処理対象としているが、これに限定されず、完全な等方性が保証される画素（等方性画素）で構成される画像を出力可能なＣＴ装置であれば、これに限定されるものではない。
【００３５】
次に、３次元ボクセル３から断面変換によりコロナル画像４を作成する。コロナル画像４は、一枚だけでなく被検体の背中から正面、あるいは正面から背中へと複数枚のコロナル画像を順次作成する。ここで、コロナル画像４を作成する際の厚みは、検出する陰影の大きさに依存するように予め設定される。
【００３６】
次に、コロナル画像４に対して、陰影検出処理５を適用する。陰影検出処理５で検出した陰影候補を含むアキシャル画像を３次元ボクセル３から切り出すために、切出位置指定６を実行する。切出処理７によって、指定された切出位置で３次元ボクセル３から切り出されたアキシャル画像８に対して、対応する陰影の特徴量を算出し、その算出結果に基づいて、陰影が診断対象部位（例えば、病気、腫瘍等）であるか否かの真偽判定９を行う。
【００３７】
次に、本実施形態の画像処理システムの構成例について、図２を用いて説明する。
【００３８】
図２は本発明の実施形態の画像処理システムの構成例を示す図である。
【００３９】
尚、図２（ａ）は、画像処理システム構成のトップビュー（上面図）であり、図２（ｂ）は、画像処理システム構成のサイドビュー（側面図）である。
【００４０】
Ｘ線発生部（Ｘ線焦点）１１から発射されたＸ線は、被検体（この場合、人体）１６を透過し、胸当て板１３及び散乱線除去グリッド（不図示）を通過した後に２次元検出器１２に到達する。
【００４１】
２次元検出器１２は、半導体センサから構成され、１画素が、例えば、２５０ｘ２５０μｍ、センサ外形が４３ｘ４３ｃｍである。また、画素数は、例えば、１７２０ｘ１７２０画素である。２次元検出器１２で取得されたデータ（被検体１６の３次元領域を投影して得られた投影データ）は、再構成部１４に転送されて画像再構成が実行される。ここで、Ｘ線発生部１１と２次元検出器１２の幾何学的配置によりファン角φ、コーン角が決定される。
【００４２】
尚、本実施形態では、検出面が正方形の２次元検出器を使用しているので、ファン角とコーン角は同じである。
【００４３】
次に、本実施形態の画像処理システムの詳細構成について、図３を用いて説明する。
【００４４】
図３は本発明の実施形態の画像処理システムの詳細構成を示すブロック図である。
【００４５】
システム全体は、コンピュータシステムにより構成される。ＢＵＳ２４は、コンピュータの内部バスと考えることができ、このＢＵＳ２４を介して制御信号やデータの送受信が行われる。
【００４６】
制御部１８は、コンピュータのＣＰＵである。メモリ２５は、例えば、ＲＯＭ、ＲＡＭ、大容量記憶装置（例えば、ハードディスク等）で構成される。ＲＯＭには、当該システムの各種処理を実行するためのプログラムや設定データが記憶されている。ＲＡＭは、データの作業領域や一時退避領域として機能する。大容量記憶装置は、投影データや画像データ等の各種データを記憶する。
【００４７】
インタフェース部２１は、入力部と出力部から構成されており、入力部は、例えば、キーボードやポインティングデバイス（マウスや入力ペン）であり、出力部は、例えば、ＣＲＴやＬＣＤ等のディスプレイである。このインタフェース部２１の入力部を介して、撮影開始指示が入力される。
【００４８】
撮影開始指示が入力されると、制御部１８からの指示により、被検体１６が固定された回転テーブル１５が回転を開始する。制御部１８は、回転テーブル１５から発生されるエンコーダ信号（不図示）を監視し、所定の一定速度及び角度に到達したか否かを確認する。所定の一定速度及び角度に到達した時点で、Ｘ線発生部１１にイネーブル信号を送信し、Ｘ線曝射を開始する。尚、このエンコーダ信号は、データの積分タイミング決定にも使用される。
【００４９】
例えば、テーブル一回転あたり２５０００パルスを発生させるエンコーダを使用し、一回転に対して１０００ビューの投影データを収集する場合、エンコーダ信号２５パルス毎に２次元検出器１２から投影データが収集されることになる。制御部１８では、このエンコーダ信号のパルスをカウントして、２５パルス毎に積分信号を発生させて、２次元検出器１２に到達したＸ線量をカウントする。
【００５０】
尚、本実施形態においては、Ｘ線は連続に発生されることを想定しているが、これに限定されるものでなく、エンコーダ信号をもとに２次元検出器１２の積分区間に合わせてパルス状のＸ線を発生させてもよい。２次元検出器１２からの投影データは、ＢＵＳ２４を介して逐次的に再構成部１４に転送される。データの転送は、回転テーブル１５が所定の回転角度を回転し、所定のビュー（画像）数が収集されるまで継続する。Ｘ線曝射が完了した直後に、最後の投影データが収集される。
【００５１】
収集された投影データは、再構成部１４に転送され、３次元ボクセルの再構成が実行される。再構成部１４は、前処理、フィルタ処理、逆投影処理から構成される。前処理は、例えば、オフセット処理、ＬＯＧ変換、ゲイン補正、欠陥補正から構成される。
【００５２】
フィルタ処理では、ラマチャンドラン関数あるいはシェップローガン関数が一般的であり、本実施形態でもこれらを使用する。フィルタ処理されたデータは、逆投影処理によって逆投影される。これらフィルタ処理から逆投影処理までのアルゴリズムは、例えば、フェルドカンプのアルゴリズムを使用している。逆投影処理が完了して、被検体の断面画像が再構成されると、その断面画像は、画像表示部１９に表示される。
【００５３】
尚、再構成部１４で実行する処理アルゴリズムは、フェルドカンプのアルゴリズムを使用するが、これに限定されるものではない。このフェルドカンプのアルゴリズムに関する参考文献には、フェルドカンプ（Ｆｅｌｄｋａｍｐ）とデイビス（Ｄａｖｉｓ）およびクレス（Ｋｒｅｓｓ）が記載した方法（「実用コーンビームアルゴリズム」（「Ｐｒａｃｔｉｃａｌ  Ｃｏｎｅ−Ｂｅａｍ  Ａｌｇｏｒｉｔｈｍ」），Ｊ．Ｏｐｔ．Ｓｏｃ．Ａｍ．Ａ１，６１２〜６１９，１９８４がある。
【００５４】
再構成された３次元ボクセルは、３次元ボクセル保存部２０に保存される。この３次元ボクセル保存部２０は、例えば、ハードディスク等の大容量記憶装置で実現される。断面変換部２２は、制御部１８に指示により、３次元ボクセル保存部２０に保存された３次元ボクセルから、所定枚数のコロナル画像を作成して、陰影検出部１７に転送する。
【００５５】
陰影検出部１７は、コロナル画像中から陰影候補を検出する。陰影候補は複数の場合もある。また、陰影検出部１７は、３次元ボクセルから陰影候補を含むアキシャル画像を切り出すために、その切出位置を算出する。断面変換部２２は、３次元ボクセル保存部２０に保存されている３次元ボクセルから、算出された切出位置でアキシャル画像を切り出して、真偽判定部２３に転送する。
【００５６】
真偽判定部２３は、当該陰影候補を含むアキシャル画像上での特徴量をもとに、陰影が診断用画像であるか否かの真偽判定を行う。診断用画像と判定された陰影は、体積をもって３次元ボクセル保存部２０上へ記録され、画像表示部１９上に表示される。
【００５７】
尚、図３では、Ｘ線発生部１１と２次元検出器１２とで、被検体画像を撮像する撮像部として機能することになるが、この撮像部が別体で構成されたシステムであっても良い。この場合、この撮像部から得られる２次元投影データを用いて、再構成部１４が３次元ボクセルを再構成することになる。
【００５８】
次に、本実施形態の画像処理システムが実行する処理について、図４を用いて説明する。
【００５９】
図４は本発明の実施形態の画像処理システムが実行する処理を示すフローチャートである。
【００６０】
インタフェース部２１の入力を介して、撮影開始指示を入力する（ステップＳ１０１）。撮影開始指示が入力されると、インタフェース部２１の出力部に、撮影準備完了画面（不図示）が表示される。制御部１８からの指示により、回転テーブル１５が回転を開始する（ステップＳ１０２）。
【００６１】
制御部１８は、回転テーブル１５から発生されるエンコーダ信号（不図示）を監視し、所定の一定速度及び角度に到達したか否かを確認する。所定の一定速度及び角度に到達した時点で、Ｘ線発生部１１にイネーブル信号を送信し、Ｘ線曝射を開始するとともに、２次元検出器１２で投影データを検出する。回転テーブル１５が所定の回転角度回転し、所定のビュー数（プロジェクション数）に達すれば、制御部１８はＸ線発生部１１に指令して、Ｘ線の曝射を停止する。その後、回転テーブル１５を減速させながら停止まで制御する（ステップＳ１０３）。
【００６２】
収集されたデータは再構成部１４に転送されて、再構成部１４によって、３次元ボクセルを再構成する（ステップＳ１０４）。そして、再構成された３次元ボクセルを３次元ボクセル保存部２０に保存する。
【００６３】
次に、３次元ボクセル保存部２０に保存された３次元ボクセルから、診断支援処理対象のコロナル画像として、所定枚数のコロナル画像の陰影検出が終了したか否かを判定する（ステップＳ１０５）。所定枚数のコロナル画像の陰影検出が終了していない場合（ステップＳ１０５でＮＯ）、ステップＳ１０６に進む。一方、所定枚数のコロナル画像の陰影検出が終了した場合（ステップＳ１０５でＹＥＳ）、ステップＳ１０９に進む。
【００６４】
尚、コロナル画像の作成方法は、インタフェース部２１の入力部により設定される、コロナル画像のスライス厚及びスライスピッチに従って決定される。また、スライス厚及びスライスピッチは、メモリ２５に予め記憶されていても良い。
【００６５】
ステップＳ１０５において、所定枚数のコロナル画像の陰影検出が終了していない場合、指定されたスライス厚及びスライスピッチにより、コロナル画像を作成する（ステップＳ１０６）。
【００６６】
次に、陰影検出部１７によって、コロナル画像に対して陰影検出を実行する（ステップＳ１０７）。
【００６７】
陰影検出方法には、種々の方法が提案されているが、本実施形態では、例えば、図５に示すＤＩＳＫフィルタ３０とＲＩＮＧフィルタ３１を用いた陰影検出方法（非特許文献１：肺がん検診用ＣＴ（ＬＳＣＴ）の診断支援システム、コンピュータ支援画像診断学会論文誌ＶＯＬ２、ＮＯ３、Ｊｕｌｙ１９９８）を使用する。
【００６８】
ここで、ＤＩＳＫフィルタ３０とＲＩＮＧフィルタ３１を用いた陰影検出処理の詳細にについて、図６を用いて説明する。
【００６９】
図６は本発明の実施形態の陰影検出処理の詳細を示すフローチャートである。
【００７０】
まず、処理対象のコロナル画像に対して、ＤＩＳＫフィルタ３０とＲＩＮＧフィルタ３１を適用したＤｉｌａｔｉｏｎ画像（ステップＳ２０１及びステップＳ２０２）を作成し、差分処理（ステップＳ２０３）を実行すると、ＲＩＮＧフィルタ３１の半径Ｒを持つ概略円形の陰影が強調される。
【００７１】
尚、Ｄｉｌａｔｉｏｎ画像とは、モフォロジカル処理の一種である（非特許文献２：Ｆｕｎｄａｍｅｎｔａｌｓ  Ｏｆ  Ｄｉｇｉｔａｌ  Ｉｍａｇｅ  Ｐｒｏｃｅｓｓｉｎｇ、Ａｎｉｌ  Ｋ．  Ｊａｉｎ著、Ｐｒｅｎｔｉｃｅ  Ｈａｌｌ社、１９８９、３８４ページ）。
【００７２】
ＤＩＳＫフィルタ３０とＲＩＮＧフィルタ３１を使用するメリットは、コロナル画像中の上下方向に走行する血管画像を効率的に排除できる点である。従来は、非特許文献１で示されるように、ＤＩＳＫフィルタとＲＩＮＧフィルタをアキシャル画像に適用しているが、アキシャル画像中の血管は陰影と同じく概略円形をしており、上記フィルタ処理では血管の殆どを検出してしまう欠点があった。
【００７３】
ＤＩＳＫフィルタ３０とＲＩＮＧフィルタ３１を適用したＤｉｌａｔｉｏｎ画像の差分画像に対し、２値化処理（ステップＳ２０４）を実行する。そして、得られる２値化画像に対して、ラベリング処理（ステップＳ２０５）を実行して、２値化画像中の小信号（ノイズ）を削除する。
【００７４】
ここで、ラベリング処理後の画像中にも血管画像に依存するノイズ、及び胸郭の輪郭画像から発生するノイズが含まれる。そこで、これらのノイズを削除するために、連結性判定処理（ステップＳ２０６）を実行する。
【００７５】
連結性判定処理は、隣り合うコロナル画像間での陰影面積と陰影重心の変動率で、腫瘍となる陰影の連結性を判定する。陰影面積と陰影重心の変動率が所定閾値よりも大きい場合は、概略円形でないものを検出している蓋然性が高いので陰影候補から除外する。
【００７６】
図４の説明に戻る。
【００７７】
図６のステップＳ２０６の連結性判定処理の結果、腫瘍となる陰影の可能性が高いと判定される陰影を陰影候補として、メモリ２５に保存する（ステップＳ１０８）。その後、ステップＳ１０５に戻る。
【００７８】
ステップＳ１０５において、所定枚数のコロナル画像からの陰影候補の検出が終了した場合（ステップＳ１０５でＹＥＳ）、陰影候補全てに対応するアキシャル画像を３次元ボクセルから切り出し、そのアキシャル画像を用いて、陰影候補が診断部位画像であるか否かの真偽判定を行うことになる。
【００７９】
そして、すべての陰影候補の真偽判定が終了したか否かを判定する（ステップＳ１０９）。すべての陰影候補の真偽判定が終了した場合（ステップＳ１０９でＹＥＳ）、ステップＳ１１３に進む。一方、すべての陰影候補の真偽判定が終了していない場合（ステップＳ１０９でＮＯ）、ステップＳ１１０に進む。
【００８０】
アキシャル画像の作成時にも、コロナル画像作成時のスライス厚及びスライスピッチが使用される。陰影候補の立体的形状は、複数のコロナル画像間での形状から把握されているので、ひとつの陰影形状を含むようにアキシャル画像を作成する（ステップＳ１１０）。尚、ここで作成するアキシャル画像は、被検体のアキシャル断面の全面に相当するアキシャル画像である必要はなく、陰影候補を含む部分アキシャル画像であればよい。
【００８１】
次に、アキシャル画像中の陰影候補の特徴量を算出して、その算出結果に基づいて、陰影候補が診断用画像であるか否かの真偽判定を行う。
【００８２】
ここで、アキシャル画像中の陰影候補の真偽判定処理の詳細について、図７を用いて説明する。
【００８３】
図７は本発明の実施形態の真偽判定処理の詳細を示すフローチャートである。
【００８４】
まず、アキシャル画像中の陰影候補の特徴量を算出する（ステップＳ２０７）。本実施形態では、アキシャル画像中の陰影候補の特徴量には、例えば、円形度を使用するが、これに限定されない。
【００８５】
但し、図６のステップＳ２０１からステップＳ２０６の処理において、陰影画像と血管画像を誤って検出する条件は、血管が肺野正面から背面に向かって走行する場合である。このような状況の場合、アキシャル画像では、血管に相当する画像は、線状の陰影として検出されるので、円形度を特徴量として採用すると良好な真偽判定を実現することができる。尚、円形度は、従来例の欄で説明したように、式（１）と式（２）で計算できる。
【００８６】
そして、算出した円形度が所定値よりも小さい場合は、腫瘍の陰影の蓋然性が低い。また、隣接するアキシャル画像間で円形度の変動率を計算することによって、変動率が所定値よりも大きい場合は、腫瘍の陰影の蓋然性が低い。このような判定方法で、ひとつの陰影候補に対する真偽判定処理（ステップＳ２１１）を実行し、真偽結果をメモリ２５に保存する。
【００８７】
図４の説明に戻る。
【００８８】
真偽判定処理を実行したら、その判定結果（真偽結果）をメモリ２５に保存する。そして、陰影候補の全てに対して、ステップＳ１１０からステップＳ１１２の処理を実行した場合、つまり、すべての陰影候補の真偽判定が終了した場合（ステップＳ１０９でＹＥＳ）、メモリ２５に記憶されている真偽結果を画像表示部１９に表示する（ステップＳ１１３）。
【００８９】
ここで、真偽結果は、ステップＳ１０８で保存した各陰影候補に対して、それぞれの陰影が腫瘍となる診断用画像であるか否かを示す真偽属性を付加したものである。真と判定された陰影領域は、３次元ボクセル上で特殊な値に置換すると、断面変換を実行した際に腫瘍に相当する画像を自動的に表示できるので、診断用画像に適している。
【００９０】
以上説明したように本発明によれば、コロナル画像は、アキシャル画像に比較して略円形状の陰影が少ないので、処理対象画像の陰影検出を行うための処理時間を短縮することができ、また、病理（腫瘍）に対する疑陽性の数を低減することができる。
【００９１】
尚、上記実施形態では、陰影検出対象の画像として、コロナル画像を用いる例を挙げて説明したが、これ以外にも、サジタル画像、オブリーク画像、それらのＲＡＹＳＵＭ画像（３次元上の一点からの厚みをもった投影画像であり、同時に医師が診断する場合に適している）等の、３次元ボクセルから得られる他の断面画像を用いることが可能である。但し、アキシャル画像は、陰影検出には利用しない。これは、上述のように、腫瘍と血管との判定が誤り易いためであり、そして、このアキシャル画像以外の他の断面画像を陰影検出に利用する点が、本願発明の特徴の１つとなっている。
【００９２】
つまり、本願発明では、陰影検出には、アキシャル画像（被検体の体軸に直交する断面画像）以外の断面画像（第１断面画像）を用いて陰影検出を行い、その陰影検出によって検出した陰影の真偽判定には、第１断面画像とは異なる第２断面画像を用いる。
【００９３】
この第１断面画像と第２断面画像の組み合わせは、種々の組み合わせが考えられるが、より実用的な組み合わせとしては、以下のものが考えられる。
【００９４】
第１断面画像      ：      第２断面画像
コロナル画像      ：      アキシャル画像あるいはサジタル画像
サジタル画像      ：      アキシャル画像あるいはコロナル画像
この組み合わせは、陰影の形状を３次元的に判定することができる組み合わせとなっており、この組み合わせを用いることで、より精度の高い陰影検出及び真偽判定を実現することができる。
【００９５】
以上、実施形態例を詳述したが、本発明は、例えば、システム、装置、方法、プログラムもしくは記憶媒体等としての実施態様をとることが可能であり、具体的には、複数の機器から構成されるシステムに適用しても良いし、また、一つの機器からなる装置に適用しても良い。
【００９６】
尚、本発明は、前述した実施形態の機能を実現するソフトウェアのプログラム（実施形態では図に示すフローチャートに対応したプログラム）を、システムあるいは装置に直接あるいは遠隔から供給し、そのシステムあるいは装置のコンピュータが該供給されたプログラムコードを読み出して実行することによっても達成される場合を含む。
【００９７】
従って、本発明の機能処理をコンピュータで実現するために、該コンピュータにインストールされるプログラムコード自体も本発明を実現するものである。つまり、本発明は、本発明の機能処理を実現するためのコンピュータプログラム自体も含まれる。
【００９８】
その場合、プログラムの機能を有していれば、オブジェクトコード、インタプリタにより実行されるプログラム、ＯＳに供給するスクリプトデータ等の形態であっても良い。
【００９９】
プログラムを供給するための記録媒体としては、例えば、フロッピー（登録商標）ディスク、ハードディスク、光ディスク、光磁気ディスク、ＭＯ、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、ＣＤ−ＲＷ、磁気テープ、不揮発性のメモリカード、ＲＯＭ、ＤＶＤ（ＤＶＤ−ＲＯＭ，ＤＶＤ−Ｒ）などがある。
【０１００】
その他、プログラムの供給方法としては、クライアントコンピュータのブラウザを用いてインターネットのホームページに接続し、該ホームページから本発明のコンピュータプログラムそのもの、もしくは圧縮され自動インストール機能を含むファイルをハードディスク等の記録媒体にダウンロードすることによっても供給できる。また、本発明のプログラムを構成するプログラムコードを複数のファイルに分割し、それぞれのファイルを異なるホームページからダウンロードすることによっても実現可能である。つまり、本発明の機能処理をコンピュータで実現するためのプログラムファイルを複数のユーザに対してダウンロードさせるＷＷＷサーバも、本発明に含まれるものである。
【０１０１】
また、本発明のプログラムを暗号化してＣＤ−ＲＯＭ等の記憶媒体に格納してユーザに配布し、所定の条件をクリアしたユーザに対し、インターネットを介してホームページから暗号化を解く鍵情報をダウンロードさせ、その鍵情報を使用することにより暗号化されたプログラムを実行してコンピュータにインストールさせて実現することも可能である。
【０１０２】
また、コンピュータが、読み出したプログラムを実行することによって、前述した実施形態の機能が実現される他、そのプログラムの指示に基づき、コンピュータ上で稼動しているＯＳなどが、実際の処理の一部または全部を行ない、その処理によっても前述した実施形態の機能が実現され得る。
【０１０３】
さらに、記録媒体から読み出されたプログラムが、コンピュータに挿入された機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに書き込まれた後、そのプログラムの指示に基づき、その機能拡張ボードや機能拡張ユニットに備わるＣＰＵなどが実際の処理の一部または全部を行ない、その処理によっても前述した実施形態の機能が実現される。
【図面の簡単な説明】
【０１０４】
【図１】本発明の画像処理システムの全体像を説明するための概念図である。
【図２】本発明の実施形態の画像処理システムの構成例を示す図である。
【図３】本発明の実施形態の画像処理システムの詳細構成を示すブロック図である。
【図４】本発明の実施形態の画像処理システムが実行する処理を示すフローチャートである。
【図５】本発明の実施形態のＤＩＳＫフィルタとＲＩＮＧフィルタの一例を示す図である。
【図６】本発明の実施形態の陰影検出処理の詳細を示すフローチャートである。
【図７】本発明の実施形態の真偽判定処理の詳細を示すフローチャートである。
【符号の説明】
【０１０５】
１１  Ｘ線発生部
１２  ２次元堅検出器
１３  胸当て板
１４  再構成部
１５  回転テーブル
１６  人体
１７  陰影検出部
１８  制御部
１９  画像表示部
２１  インタフェース部
２３  真偽判定部
２４  ＢＵＳ
３０  ＤＩＳＫフィルタ
３１  ＲＩＮＧフィルタ
【図１】
2005296348.tif 000003
【図２】
2005296348.tif 000004
【図３】
2005296348.tif 000005
【図４】
2005296348.tif 000006
【図５】
2005296348.tif 000007
【図６】
2005296348.tif 000008
【図７】
2005296348.tif 000009