(57)【特許請求の範囲】
【請求項１】
診断部位の一連の動作のＸ線動画像が格納された動画像格納部と、
前記診断部位のＸ線静止画像が格納された静止画像格納部と、
前記動画像格納部から動画像を読み出し、該動画像に基づいて前記診断部位に関する機能を示す動態情報を生成する動画像処理部と、
前記静止画像格納部から静止画像を読み出し、該静止画像から前記診断部位に関する新たな画像を生成する静止画像処理部と、
前記動画像処理部により生成された動態情報と、前記動画像及び前記静止画像処理部により生成された新たな画像のうちの少なくとも一つの形態情報とを統合し、該統合した情報を画面に表示する統合部とを備えたことを特徴とするＸ線診断支援装置。
【請求項２】
請求項１に記載のＸ線診断支援装置において、
前記動画像格納部に格納されたＸ線動画像を、胸部の後面から撮影し、吸気から呼気までの一連の呼吸動作を示す胸部Ｘ線動画像とし、
前記静止画像格納部に格納されたＸ線静止画像を、胸部のスライス面を示すＣＴ画像とし、
前記動画像処理部は、胸部Ｘ線動画像の各フレームについて横隔膜の位置を計測し、該横隔膜の位置から最大吸気のフレーム及び最大呼気のフレームを特定し、予め設定された胸部のエリア毎に、前記最大吸気のフレームと最大呼気のフレームとの間のピクセル差分値を算出して換気情報とし、該換気情報を動態情報として生成し、
前記静止画像処理部は、ＣＴ画像に基づいて胸部の後面から見たレイサム画像、コロナル像、及びサジタル像のうちの少なくとも一つを新たな画像として生成することを特徴とするＸ線診断支援装置。
【請求項３】
請求項１に記載のＸ線診断支援装置において、
前記動画像格納部に格納されたＸ線動画像を、胸部の後面から撮影し、吸気から呼気までの一連の呼吸動作を示す胸部Ｘ線動画像とし、
前記静止画像格納部に格納されたＸ線静止画像を、胸部のスライス面を示すＣＴ画像とし、
前記動画像処理部は、胸部Ｘ線動画像の各フレームについて横隔膜の位置を計測し、予め設定された胸部のエリア毎に、胸部Ｘ線動画像の各フレームについて、時系列の前後におけるフレーム間のピクセル差分値を算出して換気情報とし、該換気情報を動態情報として生成し、
前記静止画像処理部は、ＣＴ画像に基づいて胸部の後面から見たレイサム画像、コロナル像、及びサジタル像のうちの少なくとも前記レイサム画像を含む新たな画像を生成し、該生成したレイサム画像について横隔膜の位置を計測し、
前記統合部は、静止画像処理部により計測された横隔膜の位置と、動画像処理部により各フレームについて計測された横隔膜の位置とを比較し、両位置が一致する動画像のフレームを特定し、動画像処理部により生成された動態情報と、前記動画像、前記特定した動画像のフレーム、静止画像処理部により生成されたレイサム画像、コロナル像、及びサジタル像のうちの少なくとも一つの形態情報とを統合し、該統合した情報を画面に表示することを特徴とするＸ線診断支援装置。
【請求項４】
請求項３に記載のＸ線診断支援装置において、
前記統合部は、特定した動画像のフレームと、静止画像処理部により生成されたレイサム画像とを、シフト及び回転させながら、各位置におけるピクセル差分値を算出し、該ピクセル差分値の絶対値の総和が最も小さいシフト及び回転位置を求め、該位置をマッチングした位置に決定し、該マッチング位置により、前記動態情報と形態情報とを統合して画面に表示することを特徴とするＸ線診断支援装置。
【請求項５】
診断部位の一連の動作のＸ線動画像が格納された動画像格納部と、前記診断部位のＸ線静止画像が格納された静止画像格納部とを備え、前記Ｘ線動画像及び静止画像を用いて診断を支援するＸ線診断支援装置が実行するプログラムであって、該Ｘ線診断支援装置を構成するコンピュータに、
前記動画像格納部から動画像を読み出し、該動画像に基づいて前記診断部位に関する機能を示す動態情報を生成する処理と、
前記静止画像格納部から静止画像を読み出し、該静止画像から前記診断部位に関する新たな画像を生成する処理と、
前記生成された動態情報と、動画像及び新たな画像のうちの少なくとも一つの形態情報とを統合し、該統合した情報を画面に表示する処理とを実行させるＸ線診断支援プログラム。
【請求項６】
請求項５に記載のＸ線診断支援プログラムを記録した記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＣＴ（Computed Tomography）画像及びＸ線動画像を用いてコンピュータ解析を行い、動態画像診断のための情報を得るＸ線診断支援技術に関する。
【背景技術】
【０００２】
従来、呼吸器系疾患の罹患率が世界的に上昇する傾向にあり、特に、慢性閉塞性肺疾患（ＣＯＰＤ）の罹患率が高い状況にある。この慢性閉塞性肺疾患は、大気汚染や喫煙が主な原因の肺気腫や慢性気管支炎等の疾患であり、高齢者に多く発症している。このため、呼吸器系疾患の罹患率の上昇に伴い、胸部の疾患に関する定量的な機能評価を簡便に繰り返して行う方法が必要とされている。【０００３】
一般に、呼吸器系疾患の機能評価は、肺機能検査により行われる。この肺機能検査は、体内に酸素を取り入れ、体外に二酸化炭素を排出する基本的な機能を評価するために、換気、ガスと血流の分布、拡散等が正常に行われているか否かを、定性的及び定量的に検査するものである。しかしながら、肺機能検査では、局所的にまたは左右の肺を別々に検査することができず、最大努力呼吸でないと正確な検査ができないという問題があった。また、受検者には、検査のために多大な体力が必要になるという負担もあった。【０００４】
このような肺機能検査の問題を解決し、受検者に与える負担を軽減するために、肺機能を局所的に定量評価するＸ線診断装置が開発されている（例えば、特許文献１を参照）。このＸ線診断装置は、胸部をＸ線撮影して得た静止画像を記憶手段に記憶し、この静止画像における胸部の任意の部位、濃度値の計算方法、及び結果の表示方法が設定されると、その設定された部位のピクセル値（濃度値）に基づいて濃度に関する解析処理を行い、その解析処理結果を画面に表示する。また、Ｘ線診断装置は、胸部をＸ線撮影して得た時系列のフレームの静止画像を記憶手段に記憶し、この静止画像における胸部の任意の部位、濃度値の計算方法、解析対象のフレーム、及び結果の表示方法が設定されると、その設定された部位のピクセル値に基づいて、フレーム間の濃度に関する解析処理を時系列に行い、その解析処理結果を画面に表示する。このような静止画像の解析処理、及び動画像の解析処理により、読影や診断のために必要な情報を表示するようにしたから、人間の目に頼ることのない診断を実現することができ、読影医によって診断結果がばらつく可能性が低くなる。また、Ｘ線診断により得た画像に関する様々な処理技術が開示されている（例えば、特許文献２〜４を参照。）【０００５】
ところで、従来、Ｘ線照射によって被検体の断層像（ＣＴ画像）を得るＸ線ＣＴ装置が利用されている。このＸ線ＣＴ装置は、被検体にＸ線を照射し、臓器、血液、灰白質等の人体組織のＸ線吸収率の差を検出器により検出し、これをコンピュータ処理して再構成することにより、検査部位の断層面（スライス面）の画像を得るものである。このようにして得られたＣＴ画像に基づいて、医師は患者の病状等を精度高く診断することができる。【０００６】
また、従来、被検体の胸部をスクリーニング検査する装置の開発が行われている。例えば、ＦＰＤ（Flat Panel Detector）は、Ｘ線を電気信号に変換することにより画像を得る機器であり、画像を直接デジタル化する平面検出器を使用した撮影機器である。このＦＰＤを使用することにより、画像の確認、濃度調整及び幾何学的調整を行い、安定した動画像を得ることができる。例えば、吸気から呼気までの一連の呼吸動作の動画像を得て、肺尖部から横隔膜までの距離を計測して動画像を解析処理し、その結果を読影や診断のために必要な情報として表示することができる。しかしながら、ＦＰＤを使用して得た動画像を読影するには、呼吸生理学の知識が必要である。【０００７】
〔特許文献１〕特開平７−１９４５８３号公報
〔特許文献２〕特開平５−７５７９号公報
〔特許文献３〕特開平４−１３４５６８号公報
〔特許文献４〕特開昭６４−１７６３１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
このような状況の下で、Ｘ線ＣＴ装置により得られたＣＴ画像や動画像の形態情報と、動画像を解析処理して得られた動態情報とを関連付け、例えば胸部に関する形態情報と呼吸性動態情報とを重ね合わせた新たな情報とし、この重ね合わせた情報を読影や診断のために利用することが望まれている。【０００９】
そこで、本発明は、このような状況を鑑みてなされたものであり、その目的は、被検体である肺等の機能について局所的な定量評価を行い、読影や診断のために有効に利用できる評価結果を得ることが可能なＸ線診断支援装置、プログラム及び記憶媒体を提供することにある。
【課題を解決するための手段】
【００１０】
本発明によるＸ線診断支援装置は、診断部位の一連の動作のＸ線動画像が格納された動画像格納部と、前記診断部位のＸ線静止画像が格納された静止画像格納部と、前記動画像格納部から動画像を読み出し、該動画像に基づいて前記診断部位に関する機能を示す動態情報を生成する動画像処理部と、前記静止画像格納部から静止画像を読み出し、該静止画像から前記診断部位に関する新たな画像を生成する静止画像処理部と、前記動画像処理部により生成された動態情報と、前記動画像及び前記静止画像処理部により生成された新たな画像のうちの少なくとも一つの形態情報とを統合し、該統合した情報を画面に表示する統合部とを備えたことを特徴とする。【００１１】
ここで、形態情報とは、例えば、血管や気管支の太さや分布、心臓の大きさや位置、横隔膜の位置や形、骨の浸潤や骨折、肺の性質、異常物（癌、水、異物、カビ、空気の塊等）の存在を認識できる静止画像または動画像に関する情報をいう。また、動態情報とは、動画像を解析処理して得られる生態の機能情報であり、例えば、血管や気管支の拡散または収縮、心臓の拍動、横隔膜の動きまたは可動領域、肺の換気、異常物の可動性に関する情報をいう。この動態情報には、呼吸性動態情報及び心拍動静動態情報等がある。【００１２】
また、本発明によるＸ線診断支援装置は、前記動画像格納部に格納されたＸ線動画像を、胸部の後面から撮影し、吸気から呼気までの一連の呼吸動作を示す胸部Ｘ線動画像とし、前記静止画像格納部に格納されたＸ線静止画像を、胸部のスライス面を示すＣＴ画像とし、前記動画像処理部が、胸部Ｘ線動画像の各フレームについて横隔膜の位置を計測し、該横隔膜の位置から最大吸気のフレーム及び最大呼気のフレームを特定し、予め設定された胸部のエリア毎に、前記最大吸気のフレームと最大呼気のフレームとの間のピクセル差分値を算出して換気情報とし、該換気情報を動態情報として生成し、前記静止画像処理部が、ＣＴ画像に基づいて胸部の後面から見たレイサム画像、コロナル像、及びサジタル像のうちの少なくとも一つを新たな画像として生成することを特徴とする。【００１３】
また、本発明によるＸ線診断支援装置は、前記動画像格納部に格納されたＸ線動画像を、胸部の後面から撮影し、吸気から呼気までの一連の呼吸動作を示す胸部Ｘ線動画像とし、前記静止画像格納部に格納されたＸ線静止画像を、胸部のスライス面を示すＣＴ画像とし、前記動画像処理部が、胸部Ｘ線動画像の各フレームについて横隔膜の位置を計測し、予め設定された胸部のエリア毎に、胸部Ｘ線動画像の各フレームについて、時系列の前後におけるフレーム間のピクセル差分値を算出して換気情報とし、該換気情報を動態情報として生成し、前記静止画像処理部が、ＣＴ画像に基づいて胸部の後面から見たレイサム画像、コロナル像、及びサジタル像のうちの少なくとも前記レイサム画像を含む新たな画像を生成し、該生成したレイサム画像について横隔膜の位置を計測し、前記統合部が、静止画像処理部により計測された横隔膜の位置と、動画像処理部により各フレームについて計測された横隔膜の位置とを比較し、両位置が一致する動画像のフレームを特定し、動画像処理部により生成された動態情報と、前記動画像、前記特定した動画像のフレーム、静止画像処理部により生成されたレイサム画像、コロナル像、及びサジタル像のうちの少なくとも一つの形態情報とを統合し、該統合した情報を画面に表示することを特徴とする。【００１４】
また、本発明によるＸ線診断支援装置は、前記統合部が、特定した動画像のフレームと、静止画像処理部により生成されたレイサム画像とを、シフト及び回転させながら、各位置におけるピクセル差分値を算出し、該ピクセル差分値の絶対値の総和が最も小さいシフト及び回転位置を求め、該位置をマッチングした位置に決定し、該マッチング位置により、前記動態情報と形態情報とを統合して画面に表示することを特徴とする。【００１５】
また、本発明によるＸ線診断支援プログラムは、診断部位の一連の動作のＸ線動画像が格納された動画像格納部と、前記診断部位のＸ線静止画像が格納された静止画像格納部とを備え、前記Ｘ線動画像及び静止画像を用いて診断を支援するＸ線診断支援装置が実行するプログラムであって、該Ｘ線診断支援装置を構成するコンピュータに、前記動画像格納部から動画像を読み出し、該動画像に基づいて前記診断部位に関する機能を示す動態情報を生成する処理と、前記静止画像格納部から静止画像を読み出し、該静止画像から前記診断部位に関する新たな画像を生成する処理と、前記生成された動態情報と、動画像及び新たな画像のうちの少なくとも一つの形態情報とを統合し、該統合した情報を画面に表示する処理とを実行させることを特徴とする。
【発明の効果】
【００１６】
本発明によれば、被検体である肺等の機能の局所的な定量評価を行い、例えば、Ｘ線ＣＴ装置により得られたＣＴ画像である形態情報と、動画像を解析処理して得られた換気に関する動態情報とを関連付けることができるから、胸部に関する形態情報と動態情報とを統合した新たな情報を、読影や診断のために有効に利用することが可能となる。
【図面の簡単な説明】
【００１７】
【図１】本発明の実施の形態によるＸ線診断支援装置を含む全体システムの構成を示す概略図である。
【図２】Ｘ線診断支援装置の機能構成を示すブロック図である。
【図３】Ｘ線診断支援装置の動作を示すフローチャート図である。
【図４】表示器に表示された画面の一例を示す図である。
【図５】吸気から呼気までの一連の呼吸動作における肺尖部から横隔膜までの距離を示すグラフである。
【図６】換気情報を重ね合わせたコロナル像の表示例を示す図である。
【図７】換気情報の表示例を示す図である。
【符号の説明】
【００１８】
１  Ｘ線診断支援装置
２  Ｘ線動画像装置
３  Ｘ線ＣＴ装置
１１  ＣＰＵ
１２  ＲＡＭ
１３  ＲＯＭ
１４  ＨＤ
１５，１６  Ｉ／Ｆ
１７  表示器
１８  マウス
１９  キーボード
２１  動画像格納部
２２  ＣＴ画像格納部
２３  動画像処理部
２４  ＣＴ画像処理部
２５  統合部
【発明を実施するための最良の形態】
【００１９】
以下、本発明の実施の形態について、図面を用いて詳細に説明する。
〔構成〕
図１は、本発明の実施の形態によるＸ線診断支援装置を含む全体システムの構成を示す概略図である。このシステムは、Ｘ線診断支援装置１、Ｘ線動画像装置２、及びＸ線ＣＴ装置３を備えており、Ｘ線診断支援装置１の内部構成はハードウェア資源を表している。このＸ線診断支援装置１は、プログラムに従って各処理を実行するＣＰＵ１１、プログラムやデータを一時的に格納するＲＡＭ１２、ＯＳ等のシステムプログラムやシステムデータが格納されているＲＯＭ１３、各処理を実行するためのプログラムやデータ、動画像及びＣＴ画像が格納されるＨＤ１４、Ｘ線動画像装置２及びＸ線ＣＴ装置３との間の情報の入出力を中継するＩ／Ｆ１５、動画像及びＣＴ画像等を画面に表示する表示器１７、オペレータの操作を入力するマウス１８、キーボード１９、及び、表示器１７等を中継するＩ／Ｆ１６を備えている。ＣＰＵ１１は、Ｉ／Ｆ１５を介してＸ線動画像装置２及びＸ線ＣＴ装置３からそれぞれ動画像及びＣＴ画像を入力し、ＨＤ１４に格納する。また、ＣＰＵ１１は、各処理を実行するためのプログラム及びデータをＨＤ１４から読み出してＲＡＭ１２に格納し、プログラムに従って、Ｉ／Ｆ１６を介してオペレータによるマウス１８及びキーボード１９の操作により、ＨＤ１４から動画像及びＣＴ画像を読み出し、各処理を実行し、実行結果をＩ／Ｆ１６を介して表示器１７に表示する。【００２０】
図２は、図１に示したＸ線診断支援装置１の機能構成を示すブロック図である。このＸ線診断支援装置１は、動画像格納部２１、ＣＴ画像格納部２２、動画像処理部２３、ＣＴ画像処理部２４、及び、統合部２５を備えている。動画像格納部２１には、前述したＦＰＤ等のＸ線動画像装置２により撮影された胸部Ｘ線動画像が格納されている。胸部Ｘ線動画像は、吸気から呼気までの一連の呼吸動作における動画像であり、受検者毎に格納されている。ＣＴ画像格納部２２には、Ｘ線ＣＴ装置３により撮影された胸部のＣＴ画像が格納されている。ＣＴ画像は、胸部の各部を水平にスライスした静止画像であり、受検者毎に格納されている。ここで、動画像格納部２１及びＣＴ画像格納部２２は、図１におけるＨＤ１４に相当する。【００２１】
動画像処理部２３は、オペレータの操作により、動画像格納部２１から動画像を読み出し、横隔膜の位置を計測して移動量を求め、最大吸気時及び最大呼気時の動画像のフレームを特定してこれらのフレームのピクセル値の差であるピクセル差分値や、フレーム間のピクセル差分値を用いて、分割した胸部エリア毎に相対的な換気情報（換気量）を求める。また、オペレータの操作により、動画像格納部２１から読み出した動画像、及び、計測した横隔膜の位置情報を画面に表示する。【００２２】
ＣＴ画像処理部２４は、オペレータの操作により、ＣＴ画像格納部２２からＣＴ画像を読み出し、読み出したＣＴ画像間で線形補間を行い、等方的なＣＴデータを作成し、そこから、表示のためのＭＰＲ（Multi Planer Reconstruction）画像（コロナル像（胸部を前面から見た断層像）やサジタル像（胸部を側面から見た断層像）、及びＸ線写真のようなＤＲＲ（Digital Reconstruction Radiographs：レイサム（Ray Sum）画像）を作成し、このレイサム画像から横隔膜の位置を計測する。また、オペレータの操作により、ＣＴ画像格納部２２から読み出したＣＴ画像、及び、作成したコロナル像及びサジタル像を画面に表示する。【００２３】
統合部２５は、ＣＴ画像と呼吸レベルが一致する動画像のフレームを特定し、その後、ＣＴ画像と呼吸レベルが一致する動画像のフレームと、ＣＴ画像から作成されたレイサム画像との間で位置合わせを行う。また、統合部２５は、動画像処理部２３から、胸部を分割したエリア情報を示す胸部分割情報及び胸部エリア毎の相対的な換気情報を入力し、当該エリア毎の換気情報をコロナル像及び動画像に重ね合わせる。また、統合部２５は、オペレータの操作により、換気情報を重ね合わせたコロナル像及び動画像を表示する。【００２４】
〔動作〕
次に、図１及び図２に示したＸ線診断支援装置１の動作について説明する。図３は、Ｘ線診断支援装置１の動作を示すフローチャート図である。図３に示すように、Ｘ線診断支援装置１は、胸部Ｘ線動画像及びＣＴ画像を読み出し、前処理を行い、動画像とＣＴ画像との間の位置合わせ、動画像の動態解析、重ね合わせの一連の処理を行う。以下、具体的に説明する。動画像格納部２１には、胸部Ｘ線動画像が、１３４４×１３４４ピクセルの画像サイズ（解像度）、０．３２ｍｍ／ピクセル、１２ビット／４０９６の階調にて、時系列に３０フレーム分格納されているものとする。ＣＴ画像格納部２２には、５１２×５１２ピクセルの画像サイズ、０．６４ｍｍ／ピクセル、１６ビット／６５５３６の階調にて、スライス厚５ｍｍのＣＴ画像が６６枚分格納されているものとする。【００２５】
動画像処理部２３は、動画像格納部２１から前記３０フレーム分の動画像を読み出す（ステップ３０１）。そして、ＣＴ画像の画像サイズに合わせるために、１３４４×１３４４ピクセルの画像サイズを５１２×５１２ピクセルの画像サイズに調整し、３０フレーム分の当該サイズの動画像を作成する（ステップ３０２）。【００２６】
ＣＴ画像処理部２４は、ＣＴ画像格納部２２から前記６６枚分のＣＴ画像を読み出す（ステップ３０３）。そして、５ｍｍのスライス厚の２枚のＣＴ画像に対して、線形補間により画像データのボクセル化を行う（ステップ３０４）。具体的には、前記２枚のＣＴ画像を用いて、１スライスのＣＴ画像を７．８１２５スライスのＣＴ画像に補間する。つまり、読み出した５１２ピクセル×５１２ピクセル×６６スライス（枚）のＣＴ画像をそれぞれ補間することにより、５１２ピクセル×５１２ピクセル×５１２スライス（枚）のＣＴ画像を生成する。この場合、全体で、３２７．６８ｍｍ×３２７．６８ｍｍ×３２７．６８ｍｍのサイズである。そして、ＣＴ画像処理部２４は、内挿補間したＣＴ画像を用いて、ＭＰＲ表示のためのコロナル像及びサジタル像を作成すると共に、ＤＲＲ（レイサム画像）を作成する（ステップ３０５）。この場合、コロナル像、サジタル像及びレイサム画像は、それぞれ５１２×５１２ピクセルの画像サイズ及び０．６４ｍｍ／ピクセルの画像となる。【００２７】
動画像処理部２３は、ステップ３０２により作成された全フレームの動画像について、肺野領域を認識する（ステップ３０６）。また、ＣＴ画像処理部２４は、閾値処理により、肺野領域を認識する（ステップ３０６）。具体的には、ＣＴ画像処理部２４は、空気のＣＴ値を閾値として設定し、各ピクセルにおけるＣＴ値と閾値とを比較し、ＣＴ値が閾値以下のピクセルを抽出し、その抽出したピクセルの領域を肺野領域として認識する。【００２８】
動画像処理部２３は、ステップ３０２により作成された全フレームの動画像について、肺尖部及び横隔膜の位置を検出し、肺尖部から横隔膜までの距離を算出する（ステップ３０７）。また、ＣＴ画像処理部２４は、ステップ３０５により作成されたレイサム画像について、肺尖部及び横隔膜の位置を検出し、肺尖部から横隔膜までの距離を算出する（ステップ３０７）。具体的には、動画像処理部２３は１フレーム目の動画像について、ＣＴ画像処理部２４はレイサム画像について、ピクセル値に基づき、その値が大きく変化する境界を検出し、肺野領域の上部において肺尖部の座標（位置）を決定し、肺野領域の下部において横隔膜の座標を決定する。そして、動画像処理部２３は、１フレーム目の動画像から決定した肺尖部の座標の周辺を関心領域（ＲＯＩ）として設定すると共に、同様に横隔膜の座標の周辺を関心領域として設定し、２フレーム目以降の動画像について、その関心領域に対して追跡し、ピクセル値が大きく変化する境界を検出し、肺尖部及び横隔膜の座標をそれぞれ決定する。【００２９】
統合部２５は、ステップ３０７において算出した肺尖部から横隔膜までの距離を動画像処理部２３及びＣＴ画像処理部２４から入力し、レイサム画像における当該距離と動画像の各フレームにおける距離とを比較し、距離が等しいフレームを特定する（ステップ３０８）。すなわち、ＣＴ画像の呼吸レベルと動画像の呼吸レベルとが一致するフレームを、複数のフレームから特定する。【００３０】
統合部２５は、レイサム画像をＣＴ画像処理部２４から入力し、特定した動画像のフレームとレイサム画像との間で位置合わせを行う（ステップ３０９）。具体的には、入力した特定フレームとレイサム画像とを少しずつシフト及び回転させながら、各ピクセルにおける両画像のピクセル値の差分を算出し、この差分値の絶対値の総和が最も小さいシフト及び回転位置を求める。この位置を、両画像がマッチングした位置に決定する。つまり、統合部２５は、以下の数式のＲを、両画像を少しずつシフト及び回転させてそれぞれ求め、Ｒが最小となるシフト及び回転位置をマッチング位置に決定する。
【数１】
0004797173.tif 000002
ここで、Ｒはピクセル差分値の絶対値の総和、Ｆ（ｘ，ｙ）は動画像の特定フレーム、Ｇ（ｘ，ｙ）はレイサム画像、Ｒｏｔａｔｉｏｎは回転関数（入力は画像及び回転角度、出力は回転された画像）、ｘは横方向の座標（０＜ｘ＜Ｍ）、ｙは縦方向の座標（０＜ｙ＜Ｎ）、Ｍは横方向の画像サイズ（単位はピクセル）、Ｎは縦方向の画像サイズ（単位はピクセル）、ｄｘは横方向のシフト量（０＜ｄｘ＜１０，単位はピクセル）、ｄｙは縦方向のシフト量（０＜ｄｙ＜１０，単位はピクセル）、ｄｒは回転量（０＜ｄｒ＜５，単位はｄｅｇ）とする。【００３１】
尚、統合部２５は、入力したフレーム及びレイサム画像の５１２×５１２ピクセルの画像サイズを、それぞれ１２８×１２８ピクセルの画像サイズに縮小し、この縮小したフレームとレイサム画像とを用いて、マッチング位置を決定するようにしてもよい。これにより、マッチング位置を決定するための時間を短縮することができる。【００３２】
動画像処理部２３は、ステップ３０８において特定した動画像のフレームについて、左肺野領域及び右肺野領域に対してｘ軸方向に水平なエリア毎に所定の数だけ分割する（ステップ３１０）。ここで、肺機能における生態上の追従精度の観点からすると、左右の肺野領域に対してそれぞれ８分割するのが好適である。【００３３】
動画像処理部２３は、ステップ３１０において分割したエリア毎に、ピクセル値の平均を算出し、時系列に隣り合う動画像のフレーム間のピクセル値の差分をそれぞれ算出する。また、最大吸気時のピクセル値と最大呼気時のピクセル値との差分も算出する（ステップ３１１）。具体的には、動画像処理部２３は、時系列の動画像フレームにおいて、第１のフレームと第２のフレームとの間のピクセル差分値（ピクセル値の平均の差分値）、第２フレームと第３フレームとの間のピクセル差分値、・・・をそれぞれエリア毎に算出する。また、動画像処理部２３は、ステップ３０７において算出した全ての動画像フレームの肺尖部から横隔膜までの距離から、その距離が最大となるフレームを特定すると共に、その距離が最小となるフレームを特定する。ここで、距離が最大のフレームは、最大吸気時の形態を示し、距離が最小のフレームは、最大呼気時の形態を示している。そして、動画像処理部２３は、最大吸気時のフレームのピクセル値と、最大呼気時のフレームのピクセル値との差分であるピクセル差分値をエリア毎に算出する。ここで、動画像処理部２３が算出したピクセル差分値は、換気情報となる。【００３４】
統合部２５は、胸部分割情報、エリア毎のフレーム間ピクセル差分値、及び最大吸気／呼気時ピクセル差分値を動画像処理部２３から入力し、エリア毎のフレーム間ピクセル差分値に応じた濃淡色を動画像に重ね合わせて表示し、また、エリア毎の最大吸気／呼気時ピクセル差分値に応じた濃淡色をコロナル像及び動画像に重ね合わせて表示する（ステップ３１２）。画面表示の詳細については後述する。【００３５】
〔画面表示〕
図４は、図１に示した表示器１７に表示される画面の一例を示す図である。この画面は、オペレータがマウス１８やキーボード１９を操作することにより表示される。図４において、画面左上の領域４０１にはＣＴ画像、画面中央上の領域４０２には動画像、画面左下の領域４０３にはコロナル像、その右隣の領域４０４にはサジタル像、その右隣の領域４０５には領域４０２と同じ動画像、その右隣の領域４０６には肺尖部から横隔膜までの距離を示すグラフ、画面右側にはオペレータによる操作領域４０７〜４０９が表示される。以下、オペレータの操作から画面表示までの動作について説明する。【００３６】
まず、オペレータが領域４０７のホルダーまたはファイルを選択し、領域４０９の「Open images」ボタンを押下すると、動画像がＨＤ１４（動画像格納部２１）から読み出され、領域４０２，４０５に表示される。この動画像は、吸気から呼気までの一連の呼吸動作の画像である。同様に、オペレータが領域４０７において同じ受検者のホルダーまたはファイルを選択し、領域４０９の「Open CT」ボタンを押下すると、ＣＴ画像がＨＤ１４（ＣＴ画像格納部２２）から読み出され、領域４０１に表示される。領域４０８には、ｘ軸方向の画像サイズ、ｙ軸方向の画像サイズ、ピクセルサイズ、スライス厚、階調ビット数が設定及び表示される。【００３７】
オペレータが領域４０９の「Lung area」ボタンを押下すると、動画像の肺野認識が行われるとともに、動画像における横隔膜の位置を計測して移動量が求められ、動画像の各フレームにおける肺尖部から横隔膜までの距離のグラフが領域４０６に表示される。この場合、オペレータが領域４０６の下部のスクロールバーを左右に操作すると、スクロールバーの位置に対応した動画像のフレームが領域４０２，４０５に表示される。すなわち、オペレータがこのスクロールバーを左端から右端まで移動させると、領域４０２，４０５には吸気から呼気までの一連の呼吸動作の動画像が表示される。【００３８】
図５は、図４に示した領域４０６を拡大したグラフであり、吸気から呼気までの一連の呼吸動作における肺尖部から横隔膜までの距離を示している。○（丸）印の特性は右の肺における前記距離であり、×（バツ）印は左の肺における前記距離である。【００３９】
図４に戻って、オペレータが領域４０９の「Isotropic data」ボタンを押下すると、ＣＴ画像のボクセル化が行われ、コロナル像、サジタル像及びレイサム画像が作成され、コロナル像が領域４０３に、サジタル像が領域４０４にそれぞれ表示される。尚、この段階では、領域４０３，４０５には、換気情報が重ね合わせ表示されていない。また、レイサム画像は画面には表示されない。この場合、オペレータが領域４０１の下部のスクロールバーを左右に操作すると、スクロールバーの位置に対応したサジタル像が領域４０４に表示される。また、オペレータが領域４０１の右部のスクロールバーを上下に操作すると、スクロールバーの位置に対応したコロナル像が領域４０３に表示される。また、オペレータが領域４０４の右部のスクロールバーを上下に操作すると、領域４０３，４０４においてスクロールバーの位置に対応したＣＴ画像（アキシャル画像）が領域４０１に表示される。また、レイサム画像における肺尖部から横隔膜までの距離が計測される。【００４０】
オペレータが領域４０９の「Registration」のボタンを押下すると、動画像の各フレームにおける肺尖部から横隔膜までの距離と、レイサム画像における肺尖部から横隔膜までの距離とを用いて呼吸レベルが一致する動画像のフレームが特定され、そのフレームが領域４０２，４０５に表示される。そして、該画像間で位置合わせが行なわれる。また、領域４０６のスクロールバーの位置は、その特定フレームに対応した位置となる。尚、この段階では、領域４０３，４０５には、換気情報が重ね合わせ表示されていない。【００４１】
オペレータが領域４０９の「Ventilation」ボタンを押下すると、胸部エリア毎の最大吸気／呼気時のピクセル差分値が換気情報として、領域４０３のコロナル像及び領域４０５の動画像の特定フレームに重ね合わせて表示される。また、胸部エリア毎のフレーム間ピクセル差分値が換気情報として、領域４０５の動画像に重ね合わせて表示される。この場合、オペレータが領域４０６の下部のスクロールバーを左右に操作すると、スクロールバーの位置に対応した動画像が領域４０５に表示されると共に、フレーム間ピクセル差分値の換気情報が重ね合わせて表示される。この領域４０５の動画像により、吸気及び呼気状態を把握することができる。【００４２】
図６は、図４に示した領域４０３を拡大した図であり、最大吸気／呼気時のピクセル差分値が換気情報としてコロナル像に重ね合わせて表示されている。図６に示すように、換気情報は、左右の胸部をそれぞれ８等分したエリア毎に表示され、ピクセル差分値が大きさに応じて濃淡で示され、ピクセル差分値が大きいほどエリアの色は濃くなっている。また、図７は、図４に示した領域４０５を拡大した図であり、最大吸気／呼気時のピクセル差分値が換気情報として動画像の特定フレームに重ね合わせて表示されている。図６と同様に、換気情報は、左右の胸部をそれぞれ８等分したエリア毎に表示され、ピクセル差分値が大きさに応じて濃淡で示され、ピクセル差分値が大きいほどエリアの色は濃くなっている。また、中央の山型の曲線は、中央の縦線を０とした各エリアのピクセル差分値を示しており、左側の曲線は左胸部の差分値であり、右側の曲線は右胸部の差分値である。胸部の下部に近いほど、ピクセル差分値が大きくなっているから、エリアの色は濃くなり、曲線は広がっている。【００４３】
図４及び図６を参照して、オペレータが領域４０３の下部のボタンを押下することにより、領域４０３において、オリジナルのコロナル像の表示、または、コロナル像及び最大吸気／呼気時のピクセル差分値（換気情報）の重ね合わせ表示を選択することができる。また、図４及び図７を参照して、領域４０５の下部のボタンを押下することにより、領域４０５において、動画像の表示、動画像の特定フレーム及び最大吸気／呼気時のピクセル差分値（換気情報）の重ね合わせ表示、または、動画像及びフレーム間ピクセル差分値の重ね合わせ表示を選択することができる。【００４４】
ここで、図には示してないが、領域４０５において、胸部エリア毎のフレーム間ピクセル差分値が換気情報として動画像に重ね合わせて表示された場合は、ピクセル差分値がプラスのときは呼気フェーズの色、ピクセル差分値がマイナスのときは吸気フェーズの色に分け、さらに、ピクセル差分値の大きさに応じたその色の濃淡で示される。【００４５】
また、領域４０９において、オペレータの操作により、画面のコントラスト等の設定や表示、ビットマップ等のデータ保存方式の設定やデータの保存等が行われる。【００４６】
以上のように、本発明の実施の形態によるＸ線診断支援装置１によれば、統合部２５が、分割された胸部エリア毎に、胸部の形態情報と換気情報である呼吸性動態情報とを重ね合わせて画面に表示するようにした。これにより、局所的な定量評価を実現し、読影や診断のために有効に利用できる評価結果を得ることが可能となる。【００４７】
以上、実施の形態を挙げて本発明を説明したが、本発明は上記実施の形態に限定されるものではなく、本発明の精神及び意図を逸脱しない限り、種々変形が可能である。例えば、前記実施の形態では、胸部におけるピクセル差分値を換気情報として形態情報に重ね合わせたが、部位は胸部に限るものではなく、心臓等であってもよく、また、換気情報に限定するものではなく、他の動態情報、例えば、血管の拡散及び縮小に関する情報であってもよい。【００４８】
また、胸部Ｘ線動画像の画像サイズ、階調数、フレーム数等の条件、及びＣＴ画像の画像サイズ、階調数、スライス厚等の条件は、前記実施の形態に示した条件に限定されるものではない。また、動画像処理部２３により調整される画像サイズ、及びＣＴ画像処理部２４によりボクセル化されるスライス枚数は、前記実施の形態に示した画像サイズ及びスライス枚数に限定されるものではない。また、図４，６，７では、換気情報が、左右の胸部をそれぞれ８等分したエリア毎に表示されているが、８等分のエリアに限定されるものではない。【００４９】
尚、Ｘ線診断支援装置１は、図１に示したように、ＣＰＵ１１、ＲＡＭ１２等の揮発性の記憶媒体、ＲＯＭ１３等の不揮発性の記憶媒体、マウス１８やキーボード１９、ポインティングデバイス等の入力装置、画像やデータを表示する表示器１７、及び外部の装置と通信をするためのインタフェースＩ／Ｆ１５を備えたコンピュータによって構成される。Ｘ線診断支援装置１に備えた動画像処理部２３、ＣＴ画像処理部２４及び統合部２５の各機能は、これらの機能を記述したプログラムをＣＰＵ１１に実行させることによりそれぞれ実現される。また、これらのプログラムは、磁気ディスク（フロッピィーディスク、ハードディスクＨＤ１４等）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤ等）、半導体メモリ等の記憶媒体に格納して頒布することもできる。
【図１】
0004797173.tif 000003
【図２】
0004797173.tif 000004
【図３】
0004797173.tif 000005
【図５】
0004797173.tif 000006
【図６】
0004797173.tif 000007
【図７】
0004797173.tif 000008
【図４】
0004797173.tif 000009