(57)【要約】【課題】動態画像または動態画像から得られた解析結果画像と、他のモダリティで取得された機能情報画像とを容易に比較できるようにする。【解決手段】コンソール３の制御部３１によれば、放射線撮影装置１により被写体の胸部を動態撮影することにより取得された被写体の胸部の動態を示す複数のフレーム画像のそれぞれから肺野領域を特定するとともに、ＳＰＥＣＴ装置２により取得された前記被写体の肺の機能情報画像から肺野領域を特定する。そして、機能情報画像の肺野領域と複数のフレーム画像のそれぞれの肺野領域の位置合わせを行い、位置合わせされた複数のフレーム画像のそれぞれと機能情報画像とを統合し、統合画像を表示部３４に表示する。【選択図】図３
2018183493.tif 000002
【特許請求の範囲】
【請求項１】
放射線撮影装置により被写体の胸部を動態撮影することにより取得された動態画像の複数のフレーム画像のそれぞれから肺野領域を特定する第１の肺野領域特定手段と、
前記放射線撮影装置とは異なる他のモダリティにより取得された前記被写体の肺の機能情報画像から肺野領域を特定する第２の肺野領域特定手段と、
前記第２の肺野領域特定手段により特定された前記機能情報画像の肺野領域と前記第１の肺野領域特定手段により特定された前記複数のフレーム画像のそれぞれの肺野領域の位置合わせを行う変形手段と、
前記変形手段により位置合わせされた前記複数のフレーム画像のそれぞれと前記機能情報画像とを統合する統合手段と、
前記統合手段により統合された画像を表示する表示手段と、
を備える画像表示システム。
【請求項２】
放射線撮影装置により被写体の胸部を動態撮影することにより取得された動態画像の複数のフレーム画像のそれぞれから肺野領域を特定する第１の肺野領域特定手段と、
前記複数のフレーム画像間における信号値変化量又は所定の構造物の形態変化量を示す解析結果画像を生成する解析手段と、
前記放射線撮影装置とは異なる他のモダリティにより取得された前記被写体の肺の機能情報画像から肺野領域を特定する第２の肺野領域特定手段と、
前記第２の肺野領域特定手段により特定された前記機能情報画像の肺野領域と前記第１の肺野領域特定手段により特定された肺野領域に基づいて特定される前記解析結果画像の肺野領域との位置合わせを行う変形手段と、
前記変形手段により位置合わせされた前記解析結果画像と前記機能情報画像とを統合する統合手段と、
前記統合手段により統合された画像を表示する表示手段と、
を備える画像表示システム。
【請求項３】
前記解析結果画像と前記機能情報画像に基づいて異常部位を検出する検出手段を備え、
前記表示手段は、前記検出手段により検出された異常部位を前記統合手段により統合された画像、前記動態画像又は前記機能情報画像の対応する位置に表示する請求項２に記載の画像表示システム。
【請求項４】
前記機能情報画像は、肺換気シンチグラフィ又は肺血流シンチグラフィである請求項１〜３のいずれか一項に記載の画像表示システム。
【請求項５】
放射線撮影装置により被写体の胸部を動態撮影することにより取得された動態画像の複数のフレーム画像のそれぞれから肺野領域を特定する第１の肺野領域特定手段と、
前記放射線撮影装置とは異なる他のモダリティにより取得された前記被写体の肺の機能情報画像から肺野領域を特定する第２の肺野領域特定手段と、
前記第２の肺野領域特定手段により特定された前記機能情報画像の肺野領域と前記第１の肺野領域特定手段により特定された前記複数のフレーム画像のそれぞれの肺野領域の位置合わせを行う変形手段と、
前記変形手段により位置合わせされた前記複数のフレーム画像のそれぞれと前記機能情報画像とを統合する統合手段と、
を備える画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像表示システム及び画像処理装置に関する。
【背景技術】
【０００２】
従来、診断部位の動態を放射線撮影することにより得られた動態画像を解析することにより得られた動態情報（例えば、最大吸気のフレームと最大呼気のフレームとのピクセル差分値）を放射線撮影装置とは異なる他のモダリティで撮影された画像（例えば、ＣＴ装置で撮影されたコロナル像）に重ね合わせて表示する技術が提案されている（例えば、特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特許第４７９７１７３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、特許文献１に記載の技術では、他のモダリティで取得された診断部位の機能情報画像（例えば、肺シンチグラム）と動態画像を比較して位置の対応を取りながら比較読影をすることは困難である。また、他のモダリティで取得された診断対象部位の機能情報画像と動態画像を解析することにより得られる解析結果画像のミスマッチなどを比較するのはさらに困難である。
【０００５】
本発明の課題は、動態画像または動態画像から得られた解析結果画像と、他のモダリティで取得された機能情報画像とを容易に比較できるようにすることである。
【課題を解決するための手段】
【０００６】
上記課題を解決するために、請求項１に記載の発明の画像表示システムは、
放射線撮影装置により被写体の胸部を動態撮影することにより取得された動態画像の複数のフレーム画像のそれぞれから肺野領域を特定する第１の肺野領域特定手段と、
前記放射線撮影装置とは異なる他のモダリティにより取得された前記被写体の肺の機能情報画像から肺野領域を特定する第２の肺野領域特定手段と、
前記第２の肺野領域特定手段により特定された前記機能情報画像の肺野領域と前記第１の肺野領域特定手段により特定された前記複数のフレーム画像のそれぞれの肺野領域の位置合わせを行う変形手段と、
前記変形手段により位置合わせされた前記複数のフレーム画像のそれぞれと前記機能情報画像とを統合する統合手段と、
前記統合手段により統合された画像を表示する表示手段と、
を備える。
【０００７】
請求項２に記載の発明の画像表示システムは、
放射線撮影装置により被写体の胸部を動態撮影することにより取得された動態画像の複数のフレーム画像のそれぞれから肺野領域を特定する第１の肺野領域特定手段と、
前記複数のフレーム画像間における信号値変化量又は所定の構造物の形態変化量を示す解析結果画像を生成する解析手段と、
前記放射線撮影装置とは異なる他のモダリティにより取得された前記被写体の肺の機能情報画像から肺野領域を特定する第２の肺野領域特定手段と、
前記第２の肺野領域特定手段により特定された前記機能情報画像の肺野領域と前記第１の肺野領域特定手段により特定された肺野領域に基づいて特定される前記解析結果画像の肺野領域との位置合わせを行う変形手段と、
前記変形手段により位置合わせされた前記解析結果画像と前記機能情報画像とを統合する統合手段と、
前記統合手段により統合された画像を表示する表示手段と、
を備える。
【０００８】
請求項３に記載の発明は、請求項２に記載の発明において、
前記解析結果画像と前記機能情報画像に基づいて異常部位を検出する検出手段を備え、
前記表示手段は、前記検出手段により検出された異常部位を前記統合手段により統合された画像、前記動態画像又は前記機能情報画像の対応する位置に表示する。
【０００９】
請求項４に記載の発明は、請求項１〜３のいずれか一項に記載の発明において、
前記機能情報画像は、肺換気シンチグラム又は肺血流シンチグラムである。
【００１０】
請求項５に記載の発明の画像処理装置は、
放射線撮影装置により被写体の胸部を動態撮影することにより取得された動態画像の複数のフレーム画像のそれぞれから肺野領域を特定する第１の肺野領域特定手段と、
前記放射線撮影装置とは異なる他のモダリティにより取得された前記被写体の肺の機能情報画像から肺野領域を特定する第２の肺野領域特定手段と、
前記第２の肺野領域特定手段により特定された前記機能情報画像の肺野領域と前記第１の肺野領域特定手段により特定された前記複数のフレーム画像のそれぞれの肺野領域の位置合わせを行う変形手段と、
前記変形手段により位置合わせされた前記複数のフレーム画像のそれぞれと前記機能情報画像とを統合する統合手段と、
を備える。
【発明の効果】
【００１１】
本発明によれば、動態画像または動態画像から得られた解析結果画像と、他のモダリティで取得された機能情報画像とを容易に比較することが可能となる。
【図面の簡単な説明】
【００１２】
【図１】本実施形態に係る画像表示システムの全体構成を示す図である。
【図２】図１のコンソールの機能的構成を示すブロック図である。
【図３】第１の実施形態において図２の制御部により実行される画像統合処理Ａを示すフローチャートである。
【図４】動態画像のフレーム画像と機能情報画像の位置合わせを説明するための図である。
【図５】第２の実施形態において図２の制御部により実行される画像統合処理Ｂを示すフローチャートである。
【図６Ａ】肋骨の代表点における移動量の算出手法を模式的に示す図である。
【図６Ｂ】右肺野の各部における肋骨の形態変化量を示す解析結果画像の一例を示す図である。
【図７】第２の実施形態における統合画像の生成を模式的に示す図である。
【発明を実施するための形態】
【００１３】
以下、図面を参照して本発明に係る好適な実施形態について説明する。なお、本発明は、図示例に限定されるものではない。
【００１４】
＜第１の実施形態＞
（画像表示システム１００の構成）
まず、本発明に係る第１の実施形態の構成について説明する。
図１に、第１の実施形態における画像表示システム１００の全体構成例を示す。
図１に示すように、画像表示システム１００は、放射線撮影装置１と、ＳＰＥＣＴ（Single Photon Emission Computed Tomography）装置２と、コンソール３とを備えて構成されている。コンソール３は、ＬＡＮ（Local Area Network）等の通信ネットワークを介して放射線撮影装置１及びＳＰＥＣＴ装置と通信接続可能である。
【００１５】
（放射線撮影装置１）
放射線撮影装置１は、放射線源、曝射スイッチ、放射線検出部（ＦＰＤ：Flat Panel Detector）、通信部等を備えて構成され、被写体の胸部に動態撮影を行って胸部の動態画像を取得し、コンソール３に送信する装置である。
動態撮影とは、被写体に対し、Ｘ線等の放射線をパルス状にして所定時間間隔で繰り返し照射するか(パルス照射)、もしくは、低線量率にして途切れなく継続して照射する(連続照射)ことで、複数の画像を取得することをいう。動態撮影では、例えば、呼吸運動に伴う肺の膨張及び収縮の形態変化、心臓の拍動等の、周期性（サイクル）を持つ胸部の動態を撮影する。この連続撮影により得られた一連の画像を動態画像と呼ぶ。また、動態画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。
【００１６】
（ＳＰＥＣＴ装置２）
ＳＰＥＣＴ装置２は、放射性同位元素（ＲＩ=ラジオアイソトープ）を体内に吸入または静脈注射したときに体外へ放出される放射線を画像化してコンソール３に送信する装置である。本実施形態において、ＳＰＥＣＴ装置２は、肺換気の機能状態を示す機能情報画像として、キセノンガス（Xe）、クリプトンガス（Kr）又はテクネシウムガス（TC）等を被写体に吸い込ませたときに対外へ放出される放射線を画像化した肺野内のガスの分布を示す肺換気シンチグラムや、肺血流の機能状態を示す機能情報画像として、腕の静脈からテクネシウム・大凝集アルブミン（Tc-MAA）を注入したときに対外へ放出される放射線を画像化した肺野内の血流の分布を示す肺血流シンチグラムを取得する。肺換気シンチグラム及び肺血流シンチグラムの各画素は、ＲＩ値を示す。
【００１７】
（コンソール３）
コンソール３は、放射線撮影装置１による動態画像の取得動作及びＳＰＥＣＴ装置２による画像の取得動作を制御する。また、コンソール３は動態画像又はその解析結果画像と、ＳＰＥＣＴ装置２により取得された機能情報画像とを統合して統合画像を生成し、表示する画像処理装置である。
【００１８】
図２に、コンソール３の機能構成例を示す。図２に示すように、コンソール３は、制御部３１、記憶部３２、操作部３３、表示部３４、通信部３５等を備えて構成され、各部はバス３６により接続されている。
【００１９】
制御部３１は、ＣＰＵ（Central Processing Unit）、ＲＡＭ（Random Access Memory
）等により構成される。制御部３１のＣＰＵは、操作部３３の操作に応じて、記憶部３２に記憶されているシステムプログラムや各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って、コンソール３各部の動作や、放射線撮影装置１及びＳＰＥＣＴ装置２の動作を集中制御する。また、制御部３１は、展開されたプログラムに従って後述する画像統合処理Ａを始めとする各種処理を実行し、第１の肺野領域特定手段、第２の肺野領域特定手段、変形手段、統合手段として機能する。
【００２０】
記憶部３２は、不揮発性の半導体メモリーやハードディスク等により構成される。記憶部３２は、制御部３１で実行される各種プログラムやプログラムにより処理の実行に必要なパラメーター、或いは処理結果等のデータを記憶する。例えば、記憶部３２は、図３に示す画像統合処理Ａを実行するためのプログラムを記憶している。各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部３１は、当該プログラムコードに従った動作を逐次実行する。
【００２１】
また、記憶部３２には、放射線撮影装置１から送信された動態画像と、ＳＰＥＣＴ装置２から送信された機能情報画像とを被写体の患者情報、検査情報等に対応付けて記憶する。患者情報には、患者ＩＤ、患者の氏名、年齢、性別等が含まれる。検査情報には、検査ＩＤ、検査日時、検査部位（ここでは、胸部）又は検査対象の種類（例えば、肺換気、肺血流等）が含まれる。
【００２２】
操作部３３は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部３１に出力する。また、操作部３３は、表示部３４の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部３１に出力する。
【００２３】
表示部３４は、ＬＣＤ（Liquid Crystal Display）やＣＲＴ（Cathode Ray Tube）等のモニタにより構成され、制御部３１から入力される表示信号の指示に従って、操作部３３からの入力指示やデータ等を表示する。
【００２４】
通信部３５は、ＬＡＮアダプタ等を備え、ＬＡＮ等の通信ネットワークＮに接続された放射線撮影装置１、ＳＰＥＣＴ装置２を始めとする外部機器との間のデータ送受信を制御する。
【００２５】
（画像表示システム１００の動作）
次に、画像表示システム１００における動作について説明する。
図３に、操作部３３による指示に応じてコンソール３により実行される画像統合処理Ａの流れを示す。画像統合処理Ａは、制御部３１と記憶部３２に記憶されているプログラムとの協働により実行される。以下、図３を参照して画像統合処理Ａについて説明する。
【００２６】
まず、制御部３１は、処理対象の動態画像を取得する（ステップＳ１）。例えば、制御部３１は、表示部３４に画像検索画面等を表示し、記憶部３２に記憶された動態画像の中から操作部３３により入力された検索条件に合致する動態画像を記憶部３２から読み出して取得する。
【００２７】
次いで、制御部３１は、取得した動態画像の各フレーム画像から肺野領域を特定する（ステップＳ２）。
肺野の特定方法は、特に限定せず、公知の方法を適用することができる。例えば、特許第２９８７６３３号に開示されているように、Ｘ線画像では肺野領域は左右の肺部分の画像濃度が周辺より高濃度となる。よって、任意のフレーム画像の濃度ヒストグラムを作成し、その濃度ヒストグラムの形状や面積から肺野領域に該当する高濃度領域の画像部分を判断し、当該画像部分を肺野として特定すればよい。或いは、特開２００３−６６６１号公報に開示されているように、標準的な肺野領域の輪郭を定めたテンプレートを用いて任意のフレーム画像にテンプレートマッチングを行うことにより肺野領域を特定することもできる。
【００２８】
次いで、制御部３１は、取得した動態画像に対応する機能情報画像（肺換気シンチグラム及び／又は肺血流シンチグラム）を記憶部３２から読み出して取得する（ステップＳ３）。
例えば、制御部３１は、動態画像と同一患者に係る機能情報画像であって撮影日が同じ（撮影日が直近の）機能情報画像を動態画像に対応する機能情報画像として記憶部３２から読み出す。
【００２９】
次いで、制御部３１は、機能情報画像から肺野領域を特定する（ステップＳ４）。
例えば、制御部３１は、機能情報画像を予め定められた閾値で２値化して、ＲＩ値が閾値以上の領域を肺野領域として特定する。
【００３０】
次いで、制御部３１は、機能情報画像の肺野領域を変形して動態画像の各フレーム画像の肺野領域に位置合わせする（ステップＳ５）。
例えば、図４に示すように、動態画像の各フレーム画像の肺野領域の高さＹ１及び幅Ｘ１と、機能情報画像の肺野領域の高さＹ２及び幅Ｘ２とを求め、各フレーム画像のＹ１とＹ２、Ｘ１とＸ２とが等しくなるように、機能情報画像を変形させる。
【００３１】
次いで、制御部３１は、動態画像の各フレーム画像に機能情報画像を統合（合成）し、統合画像を生成する（ステップＳ６）。例えば、各フレーム画像と機能情報画像の対応する画素毎にアルファブレンディングを行って統合画像を生成する。
そして、制御部３１は、統合画像を表示部３４に表示させ（ステップＳ７）、画像統合処理Ａを終了する。統合画像は、動画表示してもよいし、フレーム画像を並べて表示してもよいし、操作部３３の操作に応じてフレーム画像を順次切り替えて表示してもよい。
【００３２】
第１の実施形態においては、動態画像の各フレーム画像に肺の機能情報画像（肺換気シンチグラム又は肺血流シンチグラム）を統合して表示する。したがって、ユーザーは、動態画像と、放射線撮影装置とは異なる他のモダリティ（ここではＳＰＥＣＴ装置２）で取得された機能情報画像とを容易に比較することが可能となる。
【００３３】
＜第２の実施形態＞
次に、本発明の第２の実施形態について説明する。
第２の実施形態では、動態画像に基づいて生成した解析結果画像と機能情報画像とを統合して表示する例について説明する。
【００３４】
第２の実施形態における構成は、コンソール３の記憶部３２に画像統合処理Ｂを実行するためのプログラムが記憶されている他は、第１の実施形態で説明したものと同様であるので説明を援用し、以下、コンソール３において実行される画像統合処理Ｂについて説明する。
【００３５】
図５は、第２の実施形態において、操作部３３による指示に応じてコンソール３により実行される画像統合処理Ｂを示すフローチャートである。画像統合処理Ｂは、制御部３１と記憶部３２に記憶されているプログラムとの協働により実行される。
【００３６】
まず、制御部３１は、処理対象の動態画像を取得する（ステップＳ１１）。ステップＳ１１の処理は、図３のステップＳ１の処理と同様であるので説明を援用する。
【００３７】
次いで、制御部３１は、取得した動態画像の各フレーム画像から肺野領域を特定する（ステップＳ１２）。ステップＳ１２の処理は、図３のステップＳ２の処理と同様であるので説明を援用する。
【００３８】
次いで、制御部３１は、動態画像を解析して解析結果画像を生成する（ステップＳ１３）。
【００３９】
例えば、制御部３１は、動態画像の肺野領域内において画素毎又は複数画素のブロック毎に信号値変化量を算出し、画素毎又はブロック毎の信号値変化量を示す解析結果画像を生成する。
例えば、検査対象の種類が肺換気である場合、まず、制御部３１は、動態画像の肺野領域内の画素毎の画素信号値（濃度値）の時間変化を時間方向のローパスフィルター（例えば、カットオフ周波数０．８５Ｈｚ）でフィルタリングする。あるいは、動態画像の各フレーム画像の肺野領域内を複数画素のブロックに分割し、各フレーム画像のブロック毎に画素信号値の代表値（例えば、平均値、中央値、最大値等）を算出し、算出した各ブロックの代表値の時間変化を時間方向のローパスフィルターでフィルタリングする。これにより、肺血流等による高周波成分を除去して肺換気の信号成分の時間変化を取得することができる。なお、ワーピング処理等によりフレーム画像間の肺野領域の位置合わせを行ってから上記処理を行うこととしてもよい。そして、各フレーム画像の画素毎又はブロック毎に、基準フレーム画像（例えば、最大呼気位のフレーム画像又は最大吸気位のフレーム画像）の対応する画素又はブロックとの画素信号値の差分値を算出し、各画素の画素信号値が基準フレーム画像からの画素信号値の差分値（つまり、信号変化量）を示す解析結果画像を生成する。
例えば、検査対象の種類が肺血流である場合、まず、制御部３１は、動態画像の肺野領域内の画素毎の画素信号値（濃度値）の時間変化を時間方向のハイパスフィルター（例えば、カットオフ周波数０．８０Ｈｚ）でフィルタリングする。あるいは、動態画像の各フレーム画像の肺野領域内を複数画素のブロックに分割し、各フレーム画像のブロック毎に画素信号値の代表値（例えば、平均値、中央値、最大値等）を算出し、算出した各ブロックの代表値の時間変化を時間方向のハイパスフィルターでフィルタリングする。これにより、肺換気等による低周波成分を除去して肺血流の信号成分の時間変化を取得することができる。なお、ワーピング処理等によりフレーム画像間の肺野領域の位置合わせを行ってから上記処理を行うこととしてもよい。そして、各フレーム画像の画素毎又はブロック毎に、基準フレーム画像（例えば、最大呼気位のフレーム画像又は最大吸気位のフレーム画像）の対応する画素又はブロックとの画素信号値の差分値を算出し、各画素の画素信号値が基準フレーム画像からの画素信号値の差分値（つまり、信号変化量）を示す解析結果画像を生成する。
【００４０】
また、例えば、制御部３１は、動態画像の肺野領域内において複数画素のブロック毎に肺野内の所定の構造物の形態変化量を算出し、ブロック毎の形態変化量を示す解析結果画像を生成する。
例えば、検査対象の種類が肺換気である場合、まず、制御部３１は、図６Ａに示すように、肋骨の輪郭を抽出し、各肋骨の代表点（例えば、縦隔側、肺野中央、外胸郭側）の各フレーム画像間の移動距離を算出する。肋骨の輪郭の抽出方法は、特に限定せず、公知の方法を適用することができる。例えば、特開平５−１７６９１９号公報に開示されているように、任意のフレーム画像の肺野領域において縦方向（頭−足の方向）にいくつもの輪郭線を定め、この輪郭線に対して予め定められたモデル関数をあてはめることで肋骨領域の輪郭部分を推定する。そして、この推定された輪郭部分に処理対象とする画像領域を複数定め、各画像領域の各画素についてＳｏｂｅｌオペレータによって勾配とその勾配に対応する方位とを求める。各画素の中で最も大きい勾配とその方位をその画像領域の勾配、方位とする。各画像領域の勾配と方位を座標として座標空間にプロットすると、座標空間では肋骨の上縁、下縁、肋骨の辺縁内部、肋骨以外といったような領域毎に集合が形成されるので、座標空間において肋骨の上縁、下縁等の辺縁に分類される画像領域を肋骨領域の輪郭として抽出することができる。
次いで、制御部３１は、図６Ｂの右肺野に示すように、肺野を上部、中部、下部に分割し、それぞれをさらに縦隔側、肺野中央、外胸郭側に分割し、分割された各領域に含まれる肋骨の代表点の平均移動量をその領域の形態変化量として算出し、肺野内の各画素の画素信号値がその画素の属する領域の形態変化量を示す解析結果画像を生成する。なお、解析結果画像は、画素信号値に応じた色を付した画像とする。
【００４１】
次いで、制御部３１は、取得した動態画像に対応する機能情報画像（肺換気シンチグラム及び／又は肺血流シンチグラム）を記憶部３２から読み出して取得する（ステップＳ１４）。
例えば、制御部３１は、動態画像と同一患者に係る機能情報画像であって撮影日が同じ（撮影日が直近の）機能情報画像を動態画像に対応する機能情報画像として記憶部３２から読み出す。
【００４２】
次いで、制御部３１は、機能情報画像から肺野領域を特定する（ステップＳ１５）。
例えば、制御部３１は、機能情報画像を予め定められた閾値で２値化して、ＲＩ値が閾値以上の領域を肺野領域として特定する。
【００４３】
次いで、制御部３１は、機能情報画像の肺野領域を変形して解析結果画像の各フレーム画像の肺野領域に位置合わせする（ステップＳ１６）。
例えば、図４に示すように、解析結果画像の各フレーム画像の肺野領域の高さＹ１及び幅Ｘ１と、機能情報画像の肺野領域の高さＹ２及び幅Ｘ２とを求め、各フレーム画像のＹ１とＹ２、Ｘ１とＸ２とが等しくなるように、機能情報画像を変形させる。なお、解析結果画像の各フレーム画像の肺野領域の位置（座標）は、ステップＳ１２において動態画像の各フレーム画像から特定された肺野領域により特定することができる。また、肺換気の解析結果画像には肺換気シンチグラムを、肺血流の解析結果画像には肺血流シンチグラムを位置合わせする。
【００４４】
次いで、制御部３１は、解析結果画像の各フレーム画像に機能情報画像を統合（合成）する（ステップＳ１７）。例えば、解析結果画像の各フレーム画像と機能情報画像の対応する画素毎にアルファブレンディングを行って統合画像を生成する。本実施形態では、肺換気の解析結果画像には肺換気シンチグラムを、肺血流の解析結果画像には肺血流シンチグラムを位置合わせする。
そして、統合画像を表示部３４に表示させ（ステップＳ１８）、画像統合処理Ｂを終了する。統合画像は、動画表示してもよいし、フレーム画像を並べて表示してもよいし、操作部３３の操作に応じてフレーム画像を順次切り替えて表示してもよい。
【００４５】
図７のＧ１は、放射線撮影装置１により取得された動態画像の任意のフレーム画像における肺野領域を囲んだ画像である。図７のＧ２は、Ｇ１のフレーム画像の解析結果画像（ここでは、肺換気の信号値変化量）である。図７のＧ３は、ＳＰＥＣＴ装置２により撮影された肺換気シンチグラムである。なお、Ｇ２とＧ３には、参考のため肺野領域の輪郭を示している。図７のＧ４は、Ｇ２とＧ３を統合した統合画像の一例を示す図である。
図７に示すように、第２の実施形態においては、動態画像の解析結果画像の各フレーム画像に肺換気シンチグラム（又は肺血流シンチグラム）を統合して表示する。したがって、ユーザーは、動態画像の解析結果画像と、放射線撮影装置とは異なる他のモダリティ（ここではＳＰＥＣＴ装置２）で取得された機能情報画像（肺換気シンチグラム、肺血流シンチグラム）とを容易に比較することが可能となる。例えば、解析結果画像と機能情報画像の表す特徴量間のミスマッチ等を容易に認識することが可能となる。
【００４６】
また、例えば、ＣＯＰＤの患者は気流制限のため呼気時間が長くなり、それに伴い肋骨の移動量にも差異が生じる可能性がある。そこで、動態画像の肋骨の形態変化量の解析結果画像と肺換気シンチグラムを統合して表示することで、診断性能を向上することができる。
【００４７】
なお、制御部３１は、検出手段として、解析結果画像と機能情報画像に基づいて異常部位を検出し、検出結果を表示部３４に表示させることとしてもよい。
【００４８】
例えば、肺換気に係る信号値変化量を示す解析結果画像の信号値変化量が所定の閾値より大きい領域（すなわち、肺胞の密度変化が大きい領域）において、肺換気シンチグラフィのＲＩ値が所定の閾値より小さい領域がある場合、その領域は肺胞の密度変化は大きいが換気量が小さいというミスマッチが起きているので、異常部位として検出する。そして、例えば、統合画像、動態画像又は機能情報画像の異常部位の領域にアノテーションを付して表示部３４に表示する。これにより、ユーザーがより容易に異常部位を認識することが可能となる。
【００４９】
また、例えば、所定の構造物（例えば、肋骨）の形態変化量を示す解析結果画像の信号値変化量が所定の閾値より大きい領域において、肺換気シンチグラフィのＲＩ値が所定の閾値より小さい領域がある場合、その領域は肋骨の動きは大きいが換気量が小さいというミスマッチが起きているので、異常部位として検出する。そして、例えば、統合画像、動態画像又は機能情報画像の異常部位の領域にアノテーションを付して表示部３４に表示する。これにより、ユーザーがより容易に異常部位を認識することが可能となる。
【００５０】
以上、本発明の第１〜第２の実施形態について説明したが、上記実施形態及び変形例における記述内容は、本発明に係る画像表示システムの好適な一例であり、これに限定されるものではない。
【００５１】
例えば、上記実施形態においては、放射線撮影装置１とは異なる他のモダリティがＳＰＥＣＴ装置２であり、機能情報画像が肺換気シンチグラム又は肺血流シンチグラムである場合を例にとり説明したが、これに限定されない。例えば、酵素ガスを造影剤として用いてＭＲＩ（Magnetic Resonance Imaging）装置で画像化することにより得られる肺の換気画像（例えば、岩澤多恵「ＭＲ換気画像の現状と将来展望」，INNERVISION（28・10）2013参照）を機能情報画像としてもよい。また、例えば、キセノン吸入２重エネルギーＣＴ画像（例えば、青木耕平ら「キセノン吸入および造影剤投与後の２重エネルギーＣＴを用いた肺切除術前後における肺形態ならびに肺換気血流比率の評価」，埼玉医科大学雑誌第４２巻第１号平成２７年８月参照）を機能情報画像としてもよい。
【００５２】
また、上記第２の実施形態においては、解析結果画像と機能情報画像とを統合して表示することとして説明したが、動態画像の各フレーム画像上に、解析結果画像と機能情報画像とを統合して表示することとしてもよい。
また、予め用意されたシェーマ上に、動態解析画像と機能情報画像とを統合して表示することとしてもよい。
【００５３】
また、上記の説明では、本発明に係るプログラムのコンピューター読み取り可能な媒体としてハードディスクや半導体の不揮発性メモリー等を使用した例を開示したが、この例に限定されない。その他のコンピューター読み取り可能な媒体として、ＣＤ−ＲＯＭ等の可搬型記録媒体を適用することが可能である。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウエーブ（搬送波）も適用される。
【００５４】
その他、画像表示システムを構成する各装置の細部構成及び細部動作に関しても、発明の趣旨を逸脱することのない範囲で適宜変更可能である。
【符号の説明】
【００５５】
１００  画像表示システム
１  放射線撮影装置
２  ＳＰＥＣＴ装置
３  コンソール
３１  制御部
３２  記憶部
３３  操作部
３４  表示部
３５  通信部
３６  バス
【図１】
2018183493.tif 000003
【図２】
2018183493.tif 000004
【図３】
2018183493.tif 000005
【図４】
2018183493.tif 000006
【図５】
2018183493.tif 000007
【図６Ａ】
2018183493.tif 000008
【図６Ｂ】
2018183493.tif 000009
【図７】
2018183493.tif 000010