(57)【要約】【課題】  本発明は、たとえば複数の障害が重なることにより生じる別種の障害を見つけることができる画像表示装置を提供する。【解決手段】  画像表示装置１００は、換気の流れを捉えた換気動画像データＤ１と、血流の流れを捉えた血流動画像データＤ２とを取得する。そして、換気動画像データＤ１から得られる換気障害を含む換気障害有無提示データＤｍと血流動画像データＤ２から得られる血流障害を含む血流障害有無提示データＤｎとを重ねることにより、異なる複数の障害状態を表す一の機能画像データＤｍｎを作成する。そして、当該機能画像データＤｍｎを用いて、異なる複数の障害状態を区別して表す機能画像を表示する。【選択図】  図１
2010057795.tif 000002
【特許請求の範囲】
【請求項１】
生体における第一の状態変化を捉えた第一の動画像データと、前記生体における第二の状態変化を捉えた第二の動画像データとを、少なくとも取得する動画像データ取得部と、
前記第一の動画像データから得られる、前記生体に関する第一の障害が存在する第一の障害領域を表示させることができる第一の画像データと、前記第二の動画像データから得られる、前記生体に関する第二の障害が存在する第二の障害領域を表示させることができる第二の画像データとを、少なくとも重ねることにより、異なる複数の障害状態を表す一の機能画像データを作成する機能画像データ作成部と、
前記機能画像データを用いて、前記異なる複数の障害状態を区別して表す機能画像を表示する表示部とを、備える、
ことを特徴とする画像表示装置。
【請求項２】
前記機能画像データ作成部は、
前記第一の障害領域と、前記第二の障害領域と、前記第一の障害領域と前記第二の障害領域との重複領域である障害重複領域とを、少なくとも表す前記機能画像データを作成する、
ことを特徴とする請求項１に記載の画像表示装置。
【請求項３】
前記動画像データ取得部は、
前記第一の障害を有すること若しくは前記第一の障害の度合いを示す前記第一の障害領域を表すことができる前記第一の画像データを、前記第一の動画像データから得られるピクセル差分値の分布状態を用いて作成し、
前記第二の障害を有すること若しくは前記第二の障害の度合いを示す前記第二の障害領域を表すことができる前記第二の画像データを、前記第二の動画像データから得られるピクセル差分値の分布状態を用いて作成し、
前記機能画像データ作成部は、
前記第一の画像データと前記第二の画像データとを用いて、前記機能画像データを作成する、
ことを特徴とする請求項２に記載の画像表示装置。
【請求項４】
前記機能画像データ作成部は、
前記第一の画像データにおける前記第一の障害および前記第二の障害を含む前記生体内の所定の部位の形状および位置と、前記第二の画像データにおける前記所定の部位の形状および位置とを合わせる画像変換処理を、前記第一の画像データおよび前記第二の画像データの少なくとも一方に施し、
前記画像変換処理後の、前記第一の画像データと前記第二の画像データとを重ね合わせることにより、前記機能画像データを作成する、
ことを特徴とする請求項３に記載の画像表示装置。
【請求項５】
前記第一の障害領域、前記第二の障害領域および前記障害重複領域と、障害の種別との関係を示すテーブルを用いて、前記機能画像データに含まれる前記第一の障害領域、前記第二の障害領域および前記障害重複領域の障害の種別を判別する障害判別部を、さらに備えている、
ことを特徴とする請求項２に記載の画像表示装置。
【請求項６】
前記生体に関するモデルデータを格納する記憶部を、さらに備えており、
前記機能画像データ作成部は、
前記モデルデータを用いて、前記機能画像データにおける前記第一の障害領域または前記第二の障害領域の範囲の補正を行う、
ことを特徴とする請求項２に記載の画像表示装置。
【請求項７】
前記第一の動画像データは、
前記肺野領域内における換気の流を捉えた動画像データであり、
前記第二の動画像データは、
肺野領域内における血液の流れを捉えた動画像データである、
ことを特徴とする請求項１に記載の画像表示装置。
【請求項８】
コンピュータによって実行されることにより、前記コンピュータを請求項１乃至請求項７の何れかに記載の画像表示装置として機能させる、
ことを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、画像表示装置およびプログラムに係る発明であり、特に、生体の障害の有無・種別を判別することができる画像表示装置およびプログラムに関する。
【背景技術】
【０００２】
医療現場では、Ｘ線等を用いて内蔵や骨格等に含まれる患部を撮影することにより、各種検査や診断が行われている。そして、近年では、デジタル技術の適用により、Ｘ線等を用いて患部の動きを捉えた動画像を比較的容易に取得することが可能となっている。
【０００３】
そして、患部の動きを捉えて診断することが有効な臓器としては、たとえば、呼吸によって臓器の形状が大きく変化する肺等が挙げられる。たとえば、肺は、疾病を伴う部分では拡大および収縮の動きが著しく低下する傾向を示す。このため、医師は、肺の挙動を動画像を通じて認識することで、診断を行うことが可能となる。
【０００４】
たとえば、３Ｈｅを投与し、肺における３Ｈｅからの磁気共鳴シグナルを検出し，シグナルの時間的変化から画像(換気画像、血流画像)を得て、換気および血流(灌流)の測定を可能としている技術が提案されている（特許文献１）。
【０００５】
【特許文献１】特表２００２−５１１３２９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
上記特許文献１では、換気画像からは換気障害のみを同定し、血流画像からは血流障害のみを同定している。このため、換気障害と血流障害とが重なることによって生じる別種の障害を見つけることができない。
【０００７】
そこで、本発明は、たとえば複数の障害が重なることにより生じる別種の障害を見つけることができる画像表示装置およびプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記の目的を達成するために、本発明に係る請求項１に記載の画像表示装置は、生体における第一の状態変化を捉えた第一の動画像データと、前記生体における第二の状態変化を捉えた第二の動画像データとを、少なくとも取得する動画像データ取得部と、前記第一の動画像データから得られる、前記生体に関する第一の障害が存在する第一の障害領域を表示させることができる第一の画像データと、前記第二の動画像データから得られる、前記生体に関する第二の障害が存在する第二の障害領域を表示させることができる第二の画像データとを、少なくとも重ねることにより、異なる複数の障害状態を表す一の機能画像データを作成する機能画像データ作成部と、前記機能画像データを用いて、前記異なる複数の障害状態を区別して表す機能画像を表示する表示部とを、備えている。
【０００９】
また、請求項２に記載の画像表示装置は、請求項１に記載の画像表示装置であって、前記機能画像データ作成部は、前記第一の障害領域と、前記第二の障害領域と、前記第一の障害領域と前記第二の障害領域との重複領域である障害重複領域とを、少なくとも表す前記機能画像データを作成する。
【００１０】
また、請求項３に記載の画像表示装置は、請求項２に記載の画像表示装置であって、前記動画像データ取得部は、前記第一の障害を有すること若しくは前記第一の障害の度合いを示す前記第一の障害領域を表すことができる前記第一の画像データを、前記第一の動画像データから得られるピクセル差分値の分布状態を用いて作成し、前記第二の障害を有すること若しくは前記第二の障害の度合いを示す前記第二の障害領域を表すことができる前記第二の画像データを、前記第二の動画像データから得られるピクセル差分値の分布状態を用いて作成し、前記機能画像データ作成部は、前記第一の画像データと前記第二の画像データとを用いて、前記機能画像データを作成する。
【００１１】
また、請求項４に記載の画像表示装置は、請求項３に記載の画像表示装置であって、前記機能画像データ作成部は、前記第一の画像データにおける前記第一の障害および前記第二の障害を含む前記生体内の所定の部位の形状および位置と、前記第二の画像データにおける前記所定の部位の形状および位置とを合わせる画像変換処理を、前記第一の画像データおよび前記第二の画像データの少なくとも一方に施し、前記画像変換処理後の、前記第一の画像データと前記第二の画像データとを重ね合わせることにより、前記機能画像データを作成する。
【００１２】
また、請求項５に記載の画像表示装置は、請求項２に記載の画像表示装置であって、前記第一の障害領域、前記第二の障害領域および前記障害重複領域と、障害の種別との関係を示すテーブルを用いて、前記機能画像データに含まれる前記第一の障害領域、前記第二の障害領域および前記障害重複領域の障害の種別を判別する障害判別部を、さらに備えている。
【００１３】
また、請求項６に記載の画像表示装置は、請求項２に記載の画像表示装置であって、前記生体に関するモデルデータを格納する記憶部を、さらに備えており、前記機能画像データ作成部は、前記モデルデータを用いて、前記機能画像データにおける前記第一の障害領域または前記第二の障害領域の範囲の補正を行う。
【００１４】
また、請求項７に記載の画像表示装置は、請求項１に記載の画像表示装置であって、前記第一の動画像データは、前記肺野領域内における換気の流れを捉えた動画像データであり、前記第二の動画像データは、肺野領域内における血液の流れを捉えた動画像データである。
【００１５】
また、請求項８に記載のプログラムは、コンピュータによって実行されることにより、前記コンピュータを請求項１乃至請求項７の何れかに記載の画像表示装置として機能させる。
【発明の効果】
【００１６】
本発明の請求項１に記載の画像表示装置または請求項８に記載のプログラムは、第一の画像データと第二の画像データとを重ねることにより、異なる複数の障害状態を表す一の機能画像データを作成し、当該機能画像データを用いて、異なる複数の障害状態を区別して表す機能画像を表示している。
【００１７】
したがって、複数の障害が重なることにより生じる別種の障害を見つけることができる画像表示装置を提供することができる。
【００１８】
また、請求項２に記載の画像表示装置は、第一の障害を示す第一の障害領域と第二の障害を示す第二の障害領域と上記障害重複領域とを、少なくとも表す機能画像データを作成する。
【００１９】
したがって、各障害領域を明確に区別された機能画像の表示が可能となる。
【００２０】
また、請求項３に記載の画像表示装置は、ピクセル差分値の分布状態（たとえばピクセル差分値の分散や当該差分値から得られる所定の絶対値量）を用いて、第一の画像データおよび第二の画像データを作成している。
【００２１】
したがって、第一の障害を含む第一の画像および第二の障害を含む第二の画像を、正確に作成することができる。
【００２２】
また、請求項４に記載の画像表示装置は、生体の所定の部位の形状・位置を合わせる画像変換処理を第一、二の画像データに施し、当該画像変換処理後の第一の画像データと第二の画像データとを重ね合わせることにより、機能画像データを作成している。
【００２３】
したがって、当該機能画像データを用いて機能画像を表示させることにより、ユーザは障害重複領域の位置を正確に把握することができる。
【００２４】
また、請求項５に記載の画像表示装置は、テーブルを用いて、機能画像データに含まれる各障害領域が如何なる障害の種別に当たるのかをを判別している。
【００２５】
したがって、第一の障害、第二の障害およびこれらの障害の重複により生じる障害の解析が可能となる。
【００２６】
また、請求項６に記載の画像表示装置は、モデルデータを用いて、機能画像データにおける第一の障害領域または第二の障害領域の範囲の補正を行う。
【００２７】
したがって、補正前の機能画像データには現れない隠れた障害部位を予測することができる。
【００２８】
また、請求項７に記載の画像表示装置では、第一の動画像データは、肺野領域内における換気の流れを捉えた動画像データであり、第二の動画像データは、肺野領域内における血液の流れを捉えた動画像データである。
【００２９】
したがって、肺野領域における換気障害、血流障害、これらの障害の重複により生じる障害を表示することができる。
【発明を実施するための最良の形態】
【００３０】
本明細書では、同一人物から取得された第一の動画像データおよび第二の動画像データから、異なる複数の障害状態を表す一の機能画像データを作成する。ここで、第一、二の動画像データは、同一生体における異なる状態の変化を各々捉えた動画像データであり、以下の説明では、第一の動画像データは、肺野領域内における換気の変化（第一の状態変化と把握できる）を捉えた動画像データとする。また、第二の動画像データは、肺野領域内における血液の流れの変化（第二の状態変化と把握できる）を捉えた動画像データであるとする。以下、第一の動画像データを換気動画像データと称し、第二の動画像データを血流動画像データと称する。
【００３１】
以下、この発明をその実施の形態を示す図面に基づいて具体的に説明する。
【００３２】
＜実施の形態１＞
図１は、本実施の形態に係る画像表示装置の構成を機能ブロック的に図示した図である。
【００３３】
図１に示す画像表示装置１００は、動画像データ取得部１０、機能画像データ作成部２０、障害判別部３０および表示部４０を、備えている。ここで、動画像データ取得部１０には、換気障害有無提示データ作成部１０Ａと血流障害有無提示データ作成部１０Ｂとが含まれている。
【００３４】
動画像データ取得部１０は、図１に示すように、換気動画像データＤ１および血流動画像データＤ２を外部から取得する。そして、動画像データ取得部１０内における換気障害有無提示データ作成部１０Ａは、取得した換気動画像データＤ１の情報から換気障害有無提示データ（第一の動画像データと把握できる）を作成する。これに対して、動画像データ取得部１０内における血流障害有無提示データ作成部１０Ｂは、取得した血流動画像データＤ２の情報から血流障害有無提示データ（第二の動画像データと把握できる）を作成する。
【００３５】
機能画像データ作成部２０では、上記換気障害有無提示データおよび上記血流障害有無提示データの情報を用いて、異なる複数の障害状態を表す一の機能画像データを作成する。
【００３６】
障害判別部３０では、上記機能画像データの情報を解析することにより、障害の種別を判断する。
【００３７】
表示部４０では、異なる複数の障害状態を表す機能画像を表示する。ここで、表示部４０は、前記障害の位置だけで無く、障害の種別情報をも付加した機能画像を表示することができる。
【００３８】
次に、各ブロック１０乃至４０の動作をより詳細に説明する。
【００３９】
｛動画像データ取得部１０の動作｝
Ｒｉｅ  ＴＡＮＡＫＡ，Ｓｈｉｇｅｒｕ  ＳＡＮＡＤＡ， ｅｔ  ａｌ，：Ｐｕｌｍｏｎａｒｙ  ｐｅｒｆｕｓｉｏｎ  ｉｍａｇｉｎｇ  ｗｉｔｈ  ａ  ｄｙｎａｍｉｃ  ｆｌａｔ−ｐａｎｅｌ  ｄｅｔｅｃｔｏｒ(ＦＰＤ)，Ｉｎｖｅｎｔ  Ｒａｄｉｏｌ．２００７（文献１）、若しくは、「第３５回日本磁気共鳴医学会大会  肺機能ワークショップ  田中理恵“動画対応フラットパネルＸ線ディテクタ装置による評価”」（文献２）に記されている手法にて、同一人物の肺野における、胸部動態Ｘ線換気動画像(換気動画像)および胸部動態Ｘ線血流動画像(血流動画像)を撮像し、動画像データ（換気動画像データおよび血流動画像データ）として獲得する。
【００４０】
｛換気障害有無提示データ作成部１０Ａの動作｝
上記換気動画像データおよび上記血流動画像データが、動画像データ取得部１０に入力されると、換気障害有無提示データ作成部１０Ａは、以下の手順により換気障害有無提示データを作成する。
【００４１】
換気障害有無提示データ作成部１０Ａでは、まず、取得した一連の換気動画像データの各々に対して、肺野形状（肺野部の輪郭）を合わせる処理を行う。図２は、当該肺野形状を合わせる処理を説明するための図である。図２において横軸は時間であり、図２の上段には、一連の連続する換気動画像データＤ１が図示されている。
【００４２】
たとえば、外部操作により、一連の換気動画像データＤ１の中から、１枚のデータを基準画像Ｂａとして選択する。当該選択は、肺野部の大きさがより大きなものを選択することがより望ましい。次に、当該選択された基準画像Ｂａ以外の画像Ｌａにおける肺野形状を、各々、基準画像Ｂａにおける肺野形状に合わせる。当該肺野形状を合わせる処理後の一連の換気動画像データＤ’１を、図２の下段に図示する。
【００４３】
ここで、肺野部の輪郭抽出の方法として、たとえば、Ｒｉｅ  ＴＡＮＡＫＡ，Ｓｈｉｇｅｒｕ  ＳＡＮＡＤＡ，ｅｔ ａｌ．：Ｃｏｍｐｕｔｅｒｉｚｘｅｄ  ｍｅｔｈｏｄ  ｆｏｒ  ｄｅｔｅｒｍｉｎｉｎｇ  ｒｅｓｐｉｒａｔｏｒｙ  ｐｈａｓｅ  ｏｎ  ｄｙｎａｍｉｃ  ｃｈｅｓｔ  ｒａｄｉｏｇｒａｐｈｓ  ｏｂｔａｉｎｅｄ  ｂｙ  ａ  ｄｙｎａｍｉｃ  ｆｌａｔ−ｐａｎｅｌ  ｄｅｔｅｃｔｏｒ(ＦＰＤ)  ｓｙｓｔｅｍ．Ｊｏｒｎａｌ  ｏｆ  Ｄｉｇｉｔａｌ  Ｉｍａｇｉｎｇ，１９(１)，４１−５１，２００６．または、Ｒｉｅ  ＴＡＮＡＫＡ，Ｓｈｉｇｅｒｕ  ＳＡＮＡＤＡ，ｅｔ  ａｌ．：Ｂｒｅａｔｈｉｎｇ  Ｃｈｅｓｔ  Ｒａｄｉｏｇｒａｐｈｙ  ｕｓｉｎｇ  ａ  Ｄｙｎａｍｉｃ  Ｆｌａｔ−Ｐａｎｅｌ  Ｄｅｔｅｃｔｏｒ(ＦＰＤ)  ｗｉｔｈ  Ｃｏｍｐｕｔｅｒ  Ａｎａｌｙｓｉｓ  ｆｏｒ  ａ  Ｓｃｒｅｅｎｉｎｇ  Ｅｘａｍｉｎａｔｉｏｎ．Ｍｅｄｉｃａｌ  Ｐｈｙｓｉｃｓ，３１(８)，２２５４−２２６２，２００４.を採用することができる。若しくは、肺野部の輪郭抽出の方法として、特開昭６３−２４０８３２号公報、または、平２−２５０１８０公報に開示される技術を採用しても良い。
【００４４】
また、上記肺野形状を合わせる方法として、モーフィング処理を用いることができる。モーフィング処理とは次のような処理である。図３は、モーフィング処理を説明するための図である。具体的に、図３は、点Ｘと点Ｘ’の対応関係の求め方を説明するための図である。
【００４５】
図３(ａ)では、点Ｐ'から点Ｑ'へのベクトル（Ｐ'Ｑ'ベクトル）、線分Ｐ'Ｑ'から点Ｘ'への垂線と合致する点Ｃ'から点Ｘ'へのベクトル（Ｃ'Ｘ'ベクトル）、および点Ｐ'から点Ｃ'へのベクトル（Ｐ'Ｃ'ベクトル）が示されている。一方、図３(ｂ)では、点Ｐから点Ｑへのベクトル（ＰＱベクトル）、線分ＰＱから点Ｘへの垂線と合致する点Ｃから点Ｘへのベクトル（ＣＸベクトル）、および点Ｐから点Ｃへのベクトル（ＰＣベクトル）が示されている。
【００４６】
ここで、２次元平面上で、Ｐ'Ｑ'ベクトルがＰＱベクトルに変換される場合に、Ｐ'Ｑ'ベクトルとＣ'Ｘ'ベクトルとＰ'Ｃ'ベクトルとの比と、ＰＱベクトルとＣＸベクトルとＰＣベクトルとの比が一定になるように、点Ｘ'が点Ｘに変換されるものとする。
【００４７】
このとき、点Ｘと点Ｘ'との間に下式(１)〜(３)の関係が成立する。
【００４８】
【数１】
2010057795.tif 000003
【００４９】
なお、ｕは、線分ＰＣの長さを線分ＰＱの長さで割った（正規化した）スカラーを示している。νは、直線ＰＱから点Ｘまでの距離を示している。Ｘ'は、点Ｘ'の座標を示している。（Ｘ−Ｐ）は、点Ｘの座標の値から点Ｐの座標の値を引いた座標の変化値、すなわち点Ｐから点Ｘへのベクトル（ＰＸベクトル）を示している。（Ｑ−Ｐ）は、点Ｑの座標の値から点Ｐの座標の値を引いた座標の変化値、すなわちＰＱベクトルを示している。Ｐｅｒ（Ｑ−Ｐ）は、ＰＱベクトルに対して垂直で且つ大きさが等しくＣＸベクトルと同じ方向を向いたベクトルを示している。Ｐｅｒ（Ｑ'−Ｐ'）は、Ｐ'Ｑ'ベクトルに対して垂直で且つ大きさが等しくＣ'Ｘ'ベクトルと同じ方向を向いたベクトルを示している。また、Ｐ'は、点Ｐ'の座標を示している。
【００５０】
上式(１)の右辺の分子は、ＰＸベクトルとＰＱベクトルとの内積を示している。上式(１)の右辺の分母は、ＰＱベクトルの大きさの２乗（線分ＰＱの長さの２乗）を示している。上式(２)の右辺の分子は、ＰＸベクトルとＰｅｒ（Ｑ−Ｐ）との内積を示している。上式（２）の右辺の分母は、ＰＱベクトルの大きさ（線分ＰＱの長さ）を示している。また、上式(３)の右辺第３項の分母は、Ｐ'Ｑ'ベクトルの大きさ（線分ＰＱの長さ）を示している。
【００５１】
ここで、換気動画像データＤ１の各画像データにおいて、肺野領域について、多数の特徴点のうちの任意の２点が点Ｐ，Ｑとされ、肺野形状を合わせる処理後の換気動画像データＤ’１の各画像データにおいて、この２点Ｐ，Ｑに対応する２点が点Ｐ'，Ｑ'と設定される。そして、換気動画像データＤ１の各画像データにおける各肺野領域内の任意の１点Ｘに対して、点Ｘと点Ｘ’との関係が多数求められる。
【００５２】
ただし、各点Ｐ，Ｑ（および点Ｐ'，Ｑ'）の設定に対して求められる点Ｘと点Ｘ'との関係は、完全には一致しない。そこで、加重平均を用いる。当該加重平均により、各点Ｘについて、下式(４)〜(６)で示すような点Ｘ'と点Ｘとの関係が求められる。なお、加重平均としては、たとえば、ＰＱベクトルの大きさや、点Ｐ，Ｑの検出精度や、ＰＱベクトルと点Ｘとの位置関係等を考慮した重み付けを用いるようにすれば良い。
【００５３】
【数２】
2010057795.tif 000004
【００５４】
ここでは、点Ｐ，Ｑをｎ通り設定されるものとし、上式(４)〜(６)では、ｎ通りの点Ｐ，Ｑの組合せのうち、ｉ番目の点Ｐ，Ｑがそれぞれ点Ｐｉ，Ｑｉと示されている。また、Ｘｉ’は、点ＰｉからＱｉへのベクトル（ＰｉＱｉベクトル）に対して、上式(３)により求められる座標である。さらに、ａ，ｂ，ｐは定数であり、ａ，ｂ，ｐの値が適宜変更されることで、換気動画像データＤ１における肺野領域の変形の仕方が調整される。また、上式(４)については、Ｘ’は点Ｘ’の座標を示し、Ｘは点Ｘの座標を示している。ωｉは上式(５)，(６)で規定される重みを示している。さらに、（Ｘｉ’−Ｘ）は、点Ｘから点Ｘi'へのベクトル（ＸＸｉ’ベクトル）を示している。
【００５５】
上式(５)で示すように、重みωｉは、点Ｐｉから点Ｑｉへのベクトル（ＰｉＱｉベクトル）のｂ乗を（ａ＋ＤＩＳＴｉ）で割って、ｐ乗することで求められる。なお、上式(６)で示す条件、ｕ＜０は、点Ｘが線分ＰＱのＰ側の延長線に対する垂線の上に位置する条件である。０≦ｕ≦１は、点Ｘが線分ＰＱの垂線の上に位置する条件である。さらに、ｕ＜１は、点Ｘが線分ＰＱのＱ側の延長線に対する垂線の上に位置する条件である。
【００５６】
上式(４)で示すように、点Ｘの座標に対して、右辺第２項のＸＸｉ’ベクトルの加重平均が加算されることで、点Ｘに対応する点Ｘ'の座標が求まる。
【００５７】
そして、上式(１)〜(６)で求まる点Ｘと点Ｘ'との関係が、換気動画像データＤ１の各画像データにおける肺野領域内の全画素について求められると、換気動画像データＤ１における肺野領域が、換気動画像データＤ’１における肺野領域に合致するように補正される。このような、２つの画像間のずれの補正は、画像領域の非線形の拡大および縮小を行う画像処理であり、たとえば、非線形のスケーリングを行う処理と称することができる。
【００５８】
なお、上記モーフィング処理は、肺野領域内の画像データに対してのみ行うことが望ましい。これは、肺野領域以外の画像データに対してモーフィング処理を行った場合、上記肺野形状を合わせる処理に長時間要する結果となるからである。
【００５９】
また、肺野形状を合わせる方法として、上記モーフィング処理以外に、特公昭６１−１４５５３号公報や特開昭６３−２７８１８３号公報等、他の公知の方法を採用することもできる。
【００６０】
さて、肺野形状を合わせる処理を行った後、換気障害有無提示データ作成部１０Ａでは、次に、当該肺野形状を合わせる処理を行った後の換気動画像データＤ’１を用いて、換気障害有無提示データを作成する。当該換気障害有無提示データは、一連の換気動画像データＤ１（またはＤ’１）から得られる。当該換気障害有無提示データを用いて画像表示させると、生体に関する換気障害（第一の障害と把握できる）が存在する第一の障害領域を表示させることができる。
【００６１】
換気障害有無提示データの作成方法を、図４を用いて具体的に説明する。図４の上段は、当該肺野形状を合わせる処理を行った後の換気動画像データＤ’１を図示している。図４の横方向は、時間軸である。また、図４の下段は、作成される換気障害有無提示データＤｍである。
【００６２】
まず、換気動画像データＤ’１の肺野領域内に存在する各ピクセルにおいて、ピクセル差分値を求める。ここで、ピクセル差分値の定義は、文献１，２に提示されている通りである（各画素の階調値）。つまり、時系列的に撮像された各換気動画像データＤ’１において、隣接する換気動画像データＤ’１同士において、対応する所定の画素間におけるピクセル値（階調値）差分を計算する。
【００６３】
具体的に、第一の換気動画像データＤ’１の第一の画素のピクセル値と、第二の換気動画像データＤ’１の第二の画素のピクセル値との差分を計算する。ここで、第二の換気動画像データＤ’１は、第一の換気動画像データＤ’１の撮像直後に撮像された動画像データである。また、第一の画素と第二の画素は、同じ座標位置に存する。当該ピクセル値の差分の計算を、換気動画像データＤ’１における肺野領域内の全画素に対して実施する。さらに、当該ピクセル値の差分の計算を、時系列的に配置されている全換気動画像データＤ’１に対して実施する。
【００６４】
図４の例では、換気動画像データｄ’１の画素ｍ１におけるピクセル差分値は、「８５」である。換気動画像データｄ’１の画素ｍ２におけるピクセル差分値は、「５」である。また、換気動画像データｄ’２の画素ｍ１におけるピクセル差分値は、「１０」である。換気動画像データｄ’２の画素ｍ２におけるピクセル差分値は、「２」である。さらに、換気動画像データｄ’ｎの画素ｍ１におけるピクセル差分値は、「−６０」である。換気動画像データｄ’ｎの画素ｍ２におけるピクセル差分値は、「−１」である。
【００６５】
次に、同じ座標位置に存する画素における上記ピクセル差分値から、分散(普遍分散)を算出する。当該分散の算出は、換気動画像データＤ’１における肺野領域内の全画素に対して実施する。そして、各画素毎に求めた当該分散の値と、予め換気障害有無提示データ作成部１０Ａが有している第一の閾値（本実施の形態では、「１００」）との、大小関係の比較（第一の閾値＜分散→正常、第一の閾値≧分散→異状）により、各画素毎に換気障害（第一の障害と把握できる）があるか否かを判断する。
【００６６】
定性的には、換気の流れが正常だと、階調が大きく変化し、分散が大きくなる。他方、換気の流れが良くないと、階調があまり変化せず、分散が小さくなる。したがって、上記第一の閾値（本実施の形態では「１００」）は、換気障害の正常・異状を区別するための閾値である。
【００６７】
たとえば、全換気動画像データｄ’１〜ｄ’ｎの画素ｍ１のピクセル差分値から、当該画素ｍ１における分散を求めた結果、分散が「５１３６．５９」であったとする。他方、全換気動画像データｄ’１〜ｄ’ｎの画素ｍ２のピクセル差分値から、当該画素ｍ２にける分散を求めた結果、分散が「５」であったとする。
【００６８】
画素ｍ１における分散の値が、「５１３６．５９」であり、第一の閾値１００より大きいので、換気障害有無提示データ作成部１０Ａは、当該画素ｍ１を正常（換気障害なし）と判断する。画素ｍ２における分散の値が、「５」であり、第一の閾値１００以下であるので、換気障害有無提示データ作成部１０Ａは、当該画素ｍ１を異状（換気障害あり）と判断する。
【００６９】
当該正常・異状の判断を肺野領域内の全画素に対して行った結果、図４で示した換気障害有無提示データ（第一の画像データと把握できる）Ｄｍが作成される。図４で示した当該図は、換気障害有無提示データＤｍを用いて画像表示させた場合に表示される映像である。当該換気障害有無提示データＤｍにおいて、砂地領域が正常（換気障害なし）と判断された領域であり、斜線領域が異状（換気障害あり）と判断された領域（第一の障害領域と把握できる）である。
【００７０】
なお、本実施の形態では、上記第一の閾値は「１００」と規定されている。しかしながら、当該第一の閾値は、対象となる人物が異なった場合、動画像データの撮像装置の仕様等に応じて変化するため、経験則によってまたは実験的に当該第一の閾値は調整されるべき値である。
【００７１】
｛血流障害有無提示データ作成部１０Ｂの動作｝
上記までは、換気障害有無提示データ作成部１０Ａにおける換気障害有無提示データＤｍ作成までの動作であるが、当該動作と同様の動作が血流障害有無提示データ作成部１０Ｂにおいても実施される。
【００７２】
具体的に、上記血流動画像データが動画像データ取得部１０に入力されると、血流障害有無提示データ作成部１０Ｂは、以下の手順により血流障害有無提示データを作成する。
【００７３】
血流障害有無提示データ作成部１０Ｂでは、まず、図２，３を用いて説明した方法と同様の方法により、肺野形状（肺野部の輪郭）を合わせる処理を行う。つまり、一連の血流動画像データＤ２から、一連の換気動画像データＤ’２を得る（図２参照）。その後、血流障害有無提示データ作成部１０Ｂでは、当該肺野形状を合わせる処理を行った後の血流動画像データＤ’２を用いて、血流障害有無提示データを作成する。
【００７４】
つまり、血流障害有無提示データは、血流動画像データＤ２（またはＤ’２）から得られる。当該血流障害有無提示データは、生体に関する血流障害（第二の障害と把握できる）が存在する第二の障害領域を表示させることができる。
【００７５】
血流障害有無提示データの作成方法を、図５を用いて具体的に説明する。
【００７６】
図５の上段は、当該肺野形状を合わせる処理を行った後の血流動画像データＤ’２を図示している。図５の横方向は、時間軸である。また、図５の下段は、作成される血流障害有無提示データＤｎである。
【００７７】
まず、血流動画像データＤ’２の肺野領域内に存在する各ピクセルにおいて、ピクセル差分値を求める。ここで、ピクセル差分値の定義は、文献１，２に提示されている通りである（各画素の階調値）。つまり、時系列的に撮像された各血流動画像データＤ’２において、隣接する血流動画像データＤ’２同士において、対応する所定の画素間におけるピクセル値（階調値）差分を計算する。
【００７８】
具体的に、第一の血流動画像データＤ’２の第一の画素のピクセル値と、第二の血流動画像データＤ’２の第二の画素のピクセル値との差分を計算する。ここで、第二の血流動画像データＤ’２は、第一の血流動画像データＤ’２の撮像直後に撮像された動画像データである。また、第一の画素と第二の画素は、同じ座標位置に存する。当該ピクセル値の差分の計算を、血流動画像データＤ’２における肺野領域内の全画素に対して実施する。さらに、当該ピクセル値の差分の計算を、時系列的に配置されている全血流動画像データＤ’２に対して実施する。
【００７９】
図５の例では、血流動画像データｄ’１１の画素ｍ３におけるピクセル差分値は、「８」である。血流動画像データｄ’１１の画素ｍ４におけるピクセル差分値は、「５」である。また、血流動画像データｄ’１２の画素ｍ３におけるピクセル差分値は、「１」である。血流動画像データｄ’１２の画素ｍ４におけるピクセル差分値は、「２」である。さらに、血流動画像データｄ’１ｎの画素ｍ３におけるピクセル差分値は、「−７」である。血流動画像データｄ’１ｎの画素ｍ４におけるピクセル差分値は、「−１」である。
【００８０】
次に、同じ座標位置に存する画素における上記ピクセル差分値から、分散(普遍分散)を算出する。当該分散の算出は、血流動画像データＤ’２における肺野領域内の全画素に対して実施する。そして、各画素毎に求めた当該分散の値と、予め血流障害有無提示データ作成部１０Ｂが有している第二の閾値（本実施の形態では、「１０」）との、大小関係の比較（第二の閾値＜分散→正常、第二の閾値≧分散→異状）により、各画素毎に血流障害（第二の障害と把握できる）があるか否かを判断する。
【００８１】
定性的には、血液の流れが正常だと、階調が大きく変化し、分散が大きくなる。他方、血液の流が良くないと、階調があまり変化せず、分散が小さくなる。したがって、上記第二の閾値（本実施の形態では「１０」）は、換気障害の正常・異状を区別するための閾値である。
【００８２】
たとえば、全血流換気動画像データｄ’１１〜ｄ’１ｎの画素ｍ３のピクセル差分値から、当該画素ｍ３における分散を求めた結果、分散が「５４．８４」であったとする。他方、全血流動画像データｄ’１１〜ｄ’１ｎの画素ｍ４のピクセル差分値から、当該画素ｍ４にける分散を求めた結果、分散が「５」であったとする。
【００８３】
画素ｍ３における分散の値が、「５４．８４」であり、第二の閾値１０より大きいので、血流障害有無提示データ作成部１０Ｂは、当該画素ｍ３を正常（血流障害なし）と判断する。画素ｍ４における分散の値が、「５」であり、第二の閾値１０以下であるので、血流障害有無提示データ作成部１０Ｂは、当該画素ｍ４を異状（血流障害あり）と判断する。
【００８４】
当該正常・異状の判断を肺野領域内の全画素に対して行った結果、図５で示した血流障害有無提示データ（第二の画像データと把握できる）Ｄｎが作成される。図５で示した当該図は、血流障害有無提示データＤｎを用いて画像表示させた場合に表示される映像である。当該血流障害有無提示データＤｎにおいて、砂地領域が正常（血流障害なし）と判断された領域であり、縦線領域が異状（血流障害あり）と判断された領域（第二の障害領域と把握できる）である。
【００８５】
なお、本実施の形態では、上記第二の閾値は「１０」と規定されている。しかしながら、当該第二の閾値は、対象となる人物が異なった場合、動画像データの撮像装置の仕様等に応じて変化するため、経験則によってまたは実験的に当該第二の閾値は調整されるべき値である。
【００８６】
なお、換気障害有無提示データＤｍおよび血流障害有無提示データＤｎの作成に際して、上記では、ピクセル差分値の分散を用いた。しかしながら、各ピクセルにおける差分値の最大値と最小値の差の絶対値（以下、単に絶対値と称する）を用いて、換気障害有無提示データＤｍおよび血流障害有無提示データＤｎを作成しても良い。
【００８７】
たとえば、一連の換気動画像データＤ’１（または一連の血流動画像データＤ’２）において、同じ座標に位置する画素Ｌに着目する。各動画像データ毎に画素Ｌのピクセル差分値を求めた結果、ピクセル差分値の最大値が「ＮＮ」で、ピクセル差分値の最小値が「ＭＭ」であったとする。この場合、上記絶対値は、｜ＮＮ−ＭＭ｜で求められる。この絶対値と、予め定められた第一の閾値（または第二の閾値）との大小関係を比較することにより、当該画素Ｌの正常・異状を判断する。
【００８８】
定性的には、換気（血液）の流れが正常だと、階調が大きく変化し、上記絶対値が大きくなる。他方、換気（血液）の流れが良くないと、階調があまり変化せず、上記絶対値が小さくなる。
【００８９】
上記分散値や上記絶対値を用いた方法から分かるように、換気障害有無提示データＤｍおよび血流障害有無提示データＤｎは、一連の換気動画像データＤ’１（または一連の血流動画像データＤ’２）における各画素のピクセル差分値の分布状態を用いて作成される。
【００９０】
｛機能画像データ作成部２０の動作｝
次に、機能画像データ作成部２０の動作、つまり異なる複数の障害状態を表す一の機能画像データを作成する方法について説明する。
【００９１】
上記で作成された、換気障害有無提示データＤｍおよび血流障害有無提示データＤｎを機能画像データ作成部２０が受信すると、当該機能画像データ作成部２０は、換気障害有無提示データＤｍから肺野領域を抽出し、さらに血流障害有無提示データＤｎから肺野領域を抽出する。
【００９２】
そして、前記で抽出した各肺野領域（生体の所定の部位と把握できる）の形状（輪郭）および座標位置を合わせる画像変換処理を実施する。ここで、換気障害有無提示データＤｍおよび血流障害有無提示データＤｎの両方若しくは何れか一方に対して当該画像変換処理を実施することにより、上記肺野領域の形状・位置合わせが可能となる。また、当該画像変換処理は、肺野領域内の各画素に対して実施される。
【００９３】
当該画像変換処理（つまり、肺野領域の形状・位置を合わせる処理）は、上記同様、モーフィング処理、特公昭６１−１４５５３号公報に係る技術や特開昭６３−２７８１８３号公報に係る技術等採用できる。なお、以下において肺野領域の形状等を合わせる処理について言及する場合には、前記で列記した各処理を採用されるものとする。
【００９４】
その後、機能画像データ作成部２０は、図６に示すように、上記肺野領域の形状・位置合わせ処理後の換気障害有無提示データＤｍおよび血流障害有無提示データＤｎ同士を、重ね合わせる。これにより、異なる複数の障害状態を表す一の機能画像データＤｍｎが作成される。ここで、肺野領域の輪郭が合致するように、両データＤｍ，Ｄｎの重ね合わせが実施される。より具体的には、肺野領域内の各画素に付与される各障害の有無を示す情報を、当該各画素毎に重ね合わせる。
【００９５】
機能画像データＤｍｎを用いた画像表示を行ったときには、図６に示すように、異なる複数の障害状態が現れる。具体的に、当該機能画像データＤｍｎにおいて、砂地領域は、換気障害も血流障害も無いと判断された領域であり、斜線領域は換気障害のみがあると判断された領域であり、縦線領域は血流障害のみがあると判断された領域であり、波線領域（障害重複領域と把握できる）は、換気障害も血流障害もあると判断された領域である。
【００９６】
｛障害判別部３０の動作｝
上記で作成された機能画像データＤｍｎを障害判別部３０が受信すると、当該障害判別部３０は、肺野領域における上記第一の障害領域と上記第二の障害領域の有無およびこれらの領域の重なり有無に応じて、当該機能画像データＤｍｎにおいて障害の種別の判別を行う。
【００９７】
たとえば、図６に示す機能画像データＤｍｎにおいて、砂地領域は、換気障害も血流障害も無いと判断され、正常として判断される。また、斜線領域は、換気障害のみがあると判断された領域であり、当該斜線領域はシャント状態であると判断される。また、縦線領域は、血流障害のみがあると判断された領域であり、当該縦線領域は肺胞死腔であると判断される。さらに、波線領域は、換気障害も血流障害もあると判断された領域であり、当該波線領域は（シャント状態＋肺胞死腔）と判断される。
【００９８】
なお、シャント状態とは、肺胞でガス交換が行われず、血流が心臓→肺動脈→肺静脈→心臓へと流れる状態を示す。肺胞死腔とは、肺胞において血流がなくガス交換できない状態を示す。
【００９９】
換気動画像データＤ１および血流動画像データＤ２が動画像データ取得部１０に入力された場合には、当該換気動画像データＤ１から作成される上記換気障害有無提示データＤｍに情報として含まれる第一の障害はシャント状態であること、および、当該血流動画像データＤ２から作成される上記血流障害有無提示データＤｎに情報として含まれる第二の障害は肺胞死腔であることが認識できるように、障害判別部３０は、設定・構成されている。
【０１００】
具体的には、図７に例示するようなテーブルが画像表示装置１００内に予め設定されている。図７は、機能画像データＤｍｎに含まれる各障害を示す領域と、障害の種別との関係が提示されている。具体的に、図７に示すテーブルでは、正常領域（図６の下段を用いて説明した砂地領域）、第一の障害領域（図６の下段を用いて説明した斜線領域）、第二の障害領域（図６の下段を用いて説明した縦線領域）および障害重複領域（図６の下段を用いて説明した波線領域）毎に、障害の各種別（シャント状態等）が列記されている。
【０１０１】
障害判別部３０は、当該図７に例示するテーブルを用いることにより、機能画像データＤｍｎにおける各障害を示す領域（第一、二の障害領域および障害重複領域）が、どの種別の障害（シャント状態、肺胞死腔、シャント状態＋肺胞死腔）に当たるかを判別している。
【０１０２】
｛表示部４０の動作｝
表示部４０が、上記機能画像データＤｍｎを受信すると、当該機能画像データＤｍｎを用いて、異なる複数の障害状態を区別して表す機能画像を表示する。具体的に、上記例では、機能画像データＤｍｎには、異なる複数の障害状態として、換気障害と血流障害と（換気障害＋血流障害）と正常状態とが含まれている。したがって、表示部４０が機能画像を表示する場合には、換気障害であることを示す第一の障害領域と血流障害であることを示す第二の障害領域と（換気障害＋血流障害）であることを示す障害重複領域と正常領域とを、明確に区別して視認できるように、当該機能画像を表示させる。
【０１０３】
たとえば、上記各領域を色分けして表示しても良く、図６に例示するように、複数の異なる模様を用いて各領域を区別しても良い。また、図７に示す表（各領域が何の障害であるかを判別した結果）を、異なる複数の障害状態を区別して表す機能画像と併せて表示しても良い。これにより、ユーザは、各領域がどのような障害を示しているのかを、一見して理解することができる。
【０１０４】
以上のように、本実施の形態に係る画像表示装置１００は、換気障害有無提示データ（第一の画像データ）Ｄｍと血流障害有無提示データ（第二の画像データ）Ｄｎとを重ねることにより、異なる複数の障害状態（シャント状態、肺胞死腔状態、シャント状態＋肺胞死腔状態）を表す一の機能画像データＤｍｎを作成し、当該機能画像データＤｍｎを用いて、異なる複数の障害状態を区別して表す機能画像を表示している。
【０１０５】
したがって、複数の障害（シャント状態、肺胞死腔）が重なることにより生じる別種の障害（シャント状態＋肺胞死腔）を見つけることができる画像表示装置１００を提供することができる。
【０１０６】
また、本実施の形態に係る画像表示装置１００は、換気障害を示す第一の障害領域と血流障害を示す第二の障害領域と上記障害重複領域とを、少なくとも表す機能画像データＤｍｎを作成する。
【０１０７】
したがって、各障害領域を明確に区別された機能画像の表示が可能となる。
【０１０８】
また、本実施の形態に係る画像表示装置１００は、ピクセル差分値の分布状態（たとえばピクセル差分値の分散や当該差分値から得られる所定の絶対値量）を用いて、換気障害有無提示データＤｍおよび血流障害有無提示データＤｎを作成している。
【０１０９】
したがって、換気障害有無提示データＤｍおよび血流障害有無提示データＤｎを、正確に作成することができる。
【０１１０】
また、本実施の形態に係る画像表示装置１００は、生体の肺野領域（所定の部位）の形状・位置を合わせる画像変換処理を換気有無提示データＤｍおよび血流障害有無提示データＤｎの少なくとも何れかに施し、当該画像変換処理後の換気有無提示データＤｍおよび血流障害有無提示データＤｎを重ね合わせることにより、機能画像データＤｍｎを作成している。
【０１１１】
したがって、当該機能画像データＤｍｎを用いて機能画像を表示させることにより、ユーザは障害重複領域の位置を正確に把握することができる。
【０１１２】
また、本実施の形態に係る画像表示装置１００は、図７に示すテーブルを用いて、機能画像データＤｍｎに含まれる各障害領域が如何なる障害の種別に当たるのかをを判別している。
【０１１３】
したがって、各障害領域がどのような障害に該当するかの解析が可能となる。
【０１１４】
＜実施の形態２＞
実施の形態１では、第一の閾値および第二の閾値と分散（または上記絶対値）とを用いて、換気障害と正常状態とを２値化で示す換気障害有無提示データＤｍと血流障害と正常状態とを２値化で示す血流障害有無提示データＤｎとを作成していた。
【０１１５】
これに対して、本実施の形態では、当該第一の閾値や当該第二の閾値や上記絶対値を用いずに、換気障害有無提示データおよび血流障害有無提示データを作成する場合について説明する。
【０１１６】
具体的に、実施の形態１と同様に、図４の上段に例示した一連の換気動画像データＤ’１を用いて、肺野領域における各画素についてピクセル差分値の分散を求める。次に、換気障害有無提示データ作成部１０Ａは、当該分散の値を各画素の情報として付与された画像データを、換気障害有無提示データとして作成する。図８は、当該作成された換気障害有無提示データＤｏを示す図である。
【０１１７】
図８に示すように、換気障害有無提示データ（第一の画像データと把握できる）Ｄｏを用いて画像を表示させた場合には、肺野領域において上記分散の値の大小がグラデーション的に表示される。つまり、換気障害（第一の障害と把握できる）の有無の度合いをグラデーション的に表示される（当該グラデーションの領域が第一の障害領域と把握できる）。
【０１１８】
図８の例では、分散の値が大きいほど階調が濃く表示され、分散の値が小さいほど階調が薄く表示されている。換気障害有無提示データＤｏから、階調がより濃い程、正常状態である可能性がより高く、階調がより薄い程、換気障害の可能性（若しくは換気障害の悪化のレベル）がより高いことが判別できる。
【０１１９】
また、図５の上段に例示した一連の血流動画像データＤ’２を用いて、肺野領域における各画素についてピクセル差分値の分散を求める。次に、血流障害有無提示データ作成部１０Ｂは、当該分散の値を各画素の情報として付与された画像データを、血流障害有無提示データとして作成する。図９は、当該作成された血流障害有無提示データＤｐを示す図である。
【０１２０】
図９に示すように、血流障害有無提示データ（第二の画像データと把握できる）Ｄｐを用いて画像を表示させた場合には、肺野領域において上記分散の値の大小がグラデーション的に表示されている。つまり、血流障害（第二の障害と把握できる）の有無の度合いをグラデーション的に表示される（当該グラデーションの領域が第二の障害領域と把握できる）。
【０１２１】
図９の例では、分散の値が大きいほど階調が濃く表示され、分散の値が小さいほど階調が薄く表示されている。血流障害有無提示データＤｐから、階調がより濃い程、正常状態である可能性がより高く、階調がより薄い程、血流障害の可能性（若しくは血流障害の悪化のレベル）がより高いことが判別できる。
【０１２２】
ここで、換気障害有無提示データＤｏのグラデーションと血流障害有無提示データＤｐのグラデーションとは、色分けされている。たとえば、換気障害有無提示データＤｏは青色のグラデーション情報として構成されており、血流障害有無提示データＤｐは赤色のグラデーション情報として構成されている。
【０１２３】
上記で作成された、換気障害有無提示データＤｏおよび血流障害有無提示データＤｐを機能画像データ作成部２０が受信する。すると、実施の形態１と同様に、当該機能画像データ作成部２０は、換気障害有無提示データＤｏから肺野領域を抽出し、さらに血流障害有無提示データＤｐから肺野領域を抽出する。
【０１２４】
そして、実施の形態１と同様に機能画像データ作成部２０は、前記で抽出した各肺野領域（生体の所定の部位と把握できる）の形状（輪郭）および座標位置を合わせる画像変換処理を実施する。ここで、当該画像変換処理は、換気障害有無提示データＤｏおよび血流障害有無提示データＤｐの両方若しくは何れか一方に対して、実施する。
【０１２５】
その後、機能画像データ作成部２０は、図１０に示すように、上記肺野領域の形状・位置合わせ処理後の換気障害有無提示データＤｏおよび血流障害有無提示データＤｐ同士を、重ね合わせる。これにより、異なる複数の障害状態を表す一の機能画像データＤｏｐが作成される。ここで、肺野領域の輪郭が合致するように、両データＤｏ，Ｄｐの重ね合わせが実施される。
【０１２６】
なお、上記の通り、換気障害有無提示データＤｏのグラデーションと血流障害有無提示データＤｐのグラデーションとは、色分けされている。したがって、機能画像データＤｏｐでは、異なる色情報に応じて異なる複数の障害状態が判別でき、グラデーション情報に応じて、障害の可能性（若しくは障害の悪化のレベル）が判別できる。
【０１２７】
その後、表示部４０は、作成された機能画像データＤｏｐに基づいて、たとえば色分けされたグラデーション表示にて、異なる複数の障害状態を区別して表す機能画像を表示する。当該機能画像が表示されることにより、この画像を基にユーザは障害状態をより詳細な部位まで判別できる。
【０１２８】
＜実施の形態３＞
本実施の形態は、生体に関するモデルデータを用いて、実施の形態１で説明した機能画像データに含まれる各障害領域（第一の障害領域、第二の障害領域および障害重複領域）の範囲を補正することができる技術に関するものである。
【０１２９】
具体的に、第一のモデルデータと機能画像データとを比較し、各障害のある部位を当該第一のモデルデータ上で特定する。そして、当該特定結果に応じて、複数の第二のモデルデータの中から所定の第二のモデルデータを抽出する。そして、当該抽出した第二のモデルデータを用いて、機能画像データにおける各障害の範囲を補正する。以下、図面を用いて具体例として説明する。
【０１３０】
本実施の形態では、図１に示した構成に加えて、上記第一，二のモデルデータ（生体に関するモデルデータと把握できる）を格納する記憶部をさらに備えている。
【０１３１】
本実施の形態では、記憶部には、第一のモデルデータとして図１１に示す「肺野血管・気管支区域モデル」が予め格納されている。さらに、当該記憶部には、第二のモデルデータとして図１２に示す「肺区画モデル」が予め格納されている。
【０１３２】
なお、図１１，１２では、説明容易化の観点で描写されている。したがって、本来のモデルと比較すると正確性には欠けている。より正確で詳細な「肺野血管・気管支区域モデル」は、新版 胸部単純Ｘ線診断  ｐ．１６に開示されている。また、より正確で詳細な「肺区画モデル」は、ｈｔｔｐ：／／ｗｗｗ．ｌａｂ．ｔｏｈｏ−ｕ．ａｃ．ｊｐ／ｍｅｄ／ｐｅｄｓ・ｔｅｘｔ／ｃｈｅｓｔｘ−ｒａｙ．ｂｍｐに開示されている。これら開示されている各モデルデータは、先見的知識により得られたものである。また、図１１では、説明簡略化のために、気管支のみを図示している、しかしながら、本来、肺野領域には、当該気管支だけが枝分かれして拡がっているので無く、動脈および静脈も枝分かれして拡がっている。
【０１３３】
図１１で示される各気管支は、Ｓ１〜Ｓ５の何れかで規定される。図１１においてＳ１で規定される気管支が枝分かれして、肺野領域の一部領域である第一の領域を占める。当該第一の領域が、図１２（Ａ）で示す斜線を付した領域Ｓ１である。同様に、図１１においてＳ２で規定される気管支が枝分かれして、肺野領域の一部領域である第二の領域を占める。当該第二の領域が、図１２（Ｂ）で示す斜線を付した領域Ｓ２である。他のＳ３〜Ｓ５においても、同様の関係が図１１と図１２との間で成立している（図１１および図１２（Ｃ），（Ｄ），（Ｅ）参照）。つまり、図１１においてＳｉで規定される気管支が枝分かれして、肺野領域の一部領域である所定の領域を占める。当該所定の領域が、図１２で示す斜線を付した領域Ｓｉである（ここで、ｉは１〜５である）。
【０１３４】
なお、上記の通り、図１１は気管支のみを図示しているが、動脈や静脈も同じＳ１〜Ｓ５で規定される。そして、たとえばＳで規定される動脈および静脈は、図１２で示す領域Ｓ１において網羅されている。つまり、Ｓｉで規定される静脈（または動脈）が枝分かれして、肺野領域の一部領域である所定の領域を占める。当該所定の領域が、図１２で示す斜線を付した領域Ｓｉである（ここで、ｉは１〜５である）。
【０１３５】
次に、本実施の形態に係る発明の動作を、例示を用いて具体的に説明する。
【０１３６】
実施の形態１で説明した方法により、機能画像データ作成部２０が機能画像データＤｍｎを作成したとする。当該作成した機能画像データＤｍｎを用いて画像表示を行った場合には、図１３に示すように、右肺野領域の正常領域（砂地領域）の内部に１部、換気障害領域（斜線領域）Ａｎが存しているものとする。なお、実施の形態１で説明したように、縦線領域は血管障害領域であり、波線領域は換気障害および血管障害が重複する領域である。
【０１３７】
次に、機能画像データ作成部２０は、記憶部から図１１に示した第一のモデルデータを読み出す。そして、機能画像データ作成部２０は、上記で作成した機能画像データＤｍｎ（図１３参照）に含まれる肺野領域および第一のモデルデータ（図１１参照）に含まれる肺野領域を抽出する。そして、当該抽出した各肺野領域の輪郭（形状）を合わせるために、機能画像データ作成部２０は、機能画像データＤｍｎに対して上記モーフィングを実施する。
【０１３８】
次に、機能画像データ作成部２０は、肺野形状の輪郭が同じとなった上記抽出した両肺野領域を対比する。具体的に、機能画像データ作成部２０は、たとえば斜線領域Ａｎが存する図１３の座標位置を特定し、図１１上において当該特定した座標位置と同じ座標位置には、Ｓ２で規定される気管支が存在すると判断する。
【０１３９】
すると、機能画像データ作成部２０は、記憶部に格納されている複数の第二のモデルデータの中から、図１２（Ｂ）に示すＳ２で規定される領域を表す第二のモデルデータを読み出す。
【０１４０】
次に、機能画像データ作成部２０は、当該読み出した第二のモデルデータから右側の肺野領域を抽出する。さらに、機能画像データ作成部２０は、当該抽出した肺野領域の輪郭と図１３に示す右側の肺野領域と輪郭とを合わせるために、当該抽出した第二のモデルデータに対して上記モーフィングを実施する。
【０１４１】
次に、機能画像データ作成部２０は、図１３に示した斜線領域Ａｎを、図１２（Ｂ）の右側肺野領域の一部領域を示す斜線領域Ｓ２と同じ範囲・位置・大きさの領域に置換する。つまり、図１２（Ｂ）に示す第二のモデルデータに合わせて、斜線領域（第一の領域であり、換気障害領域であると把握できる）Ａｓの面積を広げる補正を、機能画像データ作成部２０は行う。
【０１４２】
このような補正を行うのは、次のような考えに基づく。上記例では、機能画像データＤｍｎ内の右側肺野領域の上記領域Ａｎにおいて、Ｓ２で規定される気管支の一部に換気障害が存在することを示している。この場合には、当該気管支より先においてさらに枝分かれている細気管支（当該細気管支もＳ２で規定される）においても、換気障害のため換気が行われることはない。このことは結果として、図１２（Ｂ）の右側肺野領域のＳ２で規定される領域において、換気が正常に行われていないことを意味する。当該考えに基づき、機能画像データＤｍｎにおける上記領域Ａｎを、図１２（Ｂ）の右側肺野領域のＳ２で規定される領域へと補正する処理を行うのである。
【０１４３】
当該補正処理後の機能画像データＤ’ｍｎを用いて画像表示を行った場合には、図１４に示すように、右側の肺野領域の上部付近において、換気障害が拡大して表示される（斜線領域Ｂｎ参照）。
【０１４４】
なお、上記補正処理により、第一の障害領域（換気障害領域）が増加し、第二の障害領域（血流障害領域）と重複する領域が一部において生じたと場合には、機能画像データ作成部２０は、当該重複領域を重複障害領域（換気障害＋血流障害）に情報を変更（補正）する。
【０１４５】
また、障害判別部３０では、上記各補正処理後の機能画像データＤ’ｍｎを用いて、実施の形態１で説明した障害の種別を判別する。
【０１４６】
機能画像データＤｍｎ中において、肺野血管・気管支区域モデルの分岐次数が、たとえば２以下の気管支の領域で換気障害が現れる場合に、上記換気障害領域を変更（補正）する処理を実施するように、システムを設定・構成することが可能である。図１５は、当該気管支における分岐次数を定義している図である。なお、上記補正を行うか否かの判断基準である分岐次数の値は、対象となる人物が異なった場合には他の値が適切なこともある。したがって、分岐次数は２以下に限定するのでは無く、測定対象の人物に応じて当該分岐次数の値を調整することが望まれる。
【０１４７】
なお、上記では、換気障害の場合について換気障害領域（第一の障害領域）を補正する旨の説明を行った。しかしながら、今までの説明内容は、血流障害に関しても血流障害領域（第二の障害領域）を補正する場合にも同様に適用することができる。
【０１４８】
このように、本実施の形態では、モデルデータを用いて、機能画像データに含まれる第一，二の障害領域を拡げる補正を行うことができる。したがって、補正前の機能画像データ上では確認できなかった潜在的な障害を発見することができる。
【０１４９】
＜変形例＞
胸部動態Ｘ線画像を用いて撮像された換気動画像データおよび血流動画像データにより、換気障害有無提示データおよび血流障害有無提示データを取得する場合について説明した。しかしながら、３Ｈｅの空間的および時間的な分散を画像として表示することにより、換気および血流(灌流)の状態を示す換気障害有無提示データおよび血流障害有無提示データを取得することもできる（特表２００２−５１１３２９号公報参照）。または、ＣＴ画像より肺全体の換気分布状態を得、この換気分布状態より、換気障害有無提示データおよび血流障害有無提示データを取得することもできる（特開２００５−２８１２１号公報参照）。
【０１５０】
このように、換気障害有無提示データおよび血流障害有無提示データは様々な手法、装置によって得られる。したがって、撮影装置(モダリティ)によらずに、換気障害有無提示データおよび血流障害有無提示データを取得することが可能である。よって、複数のモダリティから、換気障害有無提示データおよび血流障害有無提示データを取得しても良い。
【０１５１】
また、換気障害有無提示データおよび血流障害有無提示データから、表示部４０で表示させる目的で、換気障害有無提示画像（たとえば図４の下段や図８に示す画像）や血流障害有無提示画像（たとえば図５の下段や図９に示す画像）を、機能画像データを作成する前等に、一旦作成しても良い。またはこれとは異なり、当該換気障害有無提示画像および血流障害有無提示画像を作成することなく、換気障害有無提示データおよび血流障害有無提示データから、直接、機能画像データを作成しても良い。
【０１５２】
また、図１の構成を、各実施の形態で説明した上記各動作を実現するために、各回路ブロックから成るハードウェア構成とすることができる。しかしながら、本発明は、ソフトウェア構成により実現することもできる。つまり、上記各動作・手順を規定するプログラムを作成し、当該プログラムを記録媒体に記憶し、当該記録媒体からコンピュータが当該プログラムを読み取り・実行する。これにより、当該コンピュータを上記各実施の形態に係る画像表示装置として機能させることができ、本願発明をソフトウェア構成として実現できる。
【図面の簡単な説明】
【０１５３】
【図１】実施の形態１に係る画像表示装置の構成を示すブロック図である。
【図２】一連の動画像データにおいて肺野領域の輪郭を統一させる動作を説明するための図である。
【図３】モーフィング処理を説明するための図である。
【図４】換気障害有無提示データの作成方法を説明するための図である。
【図５】血流障害有無提示データの作成方法を説明するための図である。
【図６】機能画像データの作成方法を説明するための図である。
【図７】各障害領域と障害の種別との関係を提示するテーブルを示す図である。
【図８】実施の形態２に係る画像表示装置の動作を説明するための図である。
【図９】実施の形態２に係る画像表示装置の動作を説明するための図である。
【図１０】実施の形態２に係る画像表示装置の動作を説明するための図である。
【図１１】第一のモデルデータを示す図である。
【図１２】第二のモデルデータを示す図である。
【図１３】第一の障害領域の補正処理前の機能画像データを示す図である。
【図１４】第一の障害領域の補正処理後の機能画像データを示す図である。
【図１５】気管支の分岐次数を説明するための図である。
【符号の説明】
【０１５４】
１０  動画像データ取得部
１０Ａ  換気障害有無提示データ作成部
１０Ｂ  血流障害有無提示データ作成部
２０  機能画像データ作成部
３０  障害判別部
４０  表示部
１００  画像表示装置
Ｄ１，Ｄ’１  換気動画像データ
Ｄ２，Ｄ’２  血流動画像データ
Ｄｍ，Ｄｏ  換気障害有無提示データ
Ｄｎ，Ｄｐ  血流障害有無提示データ
Ｄｍｎ，Ｄｏｐ  機能画像データ
Ｄ’ｍｎ  補正処理後の機能画像データ
【図１】
2010057795.tif 000005
【図２】
2010057795.tif 000006
【図３】
2010057795.tif 000007
【図４】
2010057795.tif 000008
【図５】
2010057795.tif 000009
【図６】
2010057795.tif 000010
【図７】
2010057795.tif 000011
【図８】
2010057795.tif 000012
【図９】
2010057795.tif 000013
【図１０】
2010057795.tif 000014
【図１１】
2010057795.tif 000015
【図１２】
2010057795.tif 000016
【図１３】
2010057795.tif 000017
【図１４】
2010057795.tif 000018
【図１５】
2010057795.tif 000019