(57)【要約】【課題】ＣＴ画像の表示において、一方のＣＴ画像で指定した病巣部が時間経過により、どのように変化したかを容易に比較できるようにする。【解決手段】  画像表示装置は、異なる時期に撮影された撮影箇所が対応するＣＴ画像のうち、第１のＣＴ画像上で指定された箇所に対応する箇所を、第２の時期に撮影された第２のＣＴ画像上で、該第２のＣＴ画像の他の箇所とは異なる態様で表示する表示部を有することを特徴とする。【選択図】図１５
2015208539.tif 000002
【特許請求の範囲】
【請求項１】
異なる時期に撮影された撮影箇所が対応するＣＴ画像のうち、第１のＣＴ画像上で指定された箇所に対応する箇所を、第２の時期に撮影された第２のＣＴ画像上で、該第２のＣＴ画像の他の箇所とは異なる態様で表示する表示部、
を有することを特徴とする画像表示装置。
【請求項２】
前記第１のＣＴ画像上の箇所に対応する、前記第２のＣＴ画像上の箇所を特定する特定部を更に有し、
前記特定部は、
前記第１のＣＴ画像及び前記第２のＣＴ画像のそれぞれから、臓器の辺縁部分を特定し、
前記第１のＣＴ画像上の臓器の辺縁部分から、該第１のＣＴ画像で指定された箇所を通り対向する２点の辺縁部分を特定し、それらの辺縁部分の特徴量をもとに、前記第２のＣＴ画像の臓器の辺縁部分から前記第１のＣＴ画像で特定した辺縁部分に対応する辺縁部分を特定し、
前記第１のＣＴ画像で指定された箇所から、前記第１のＣＴ画像において特定した各辺縁部分までの距離を特定し、
前記第２のＣＴ画像から特定した前記第１のＣＴ画像の各々の辺縁部分に対応する辺縁部分から、特定した距離の比を満たす箇所を特定する、
ことを特徴とする請求項１に記載の画像表示装置。
【請求項３】
異なる時期に撮影された撮影箇所が対応するＣＴ画像のうち、第１のＣＴ画像上で指定された箇所に対応する箇所を、第２の時期に撮影された第２のＣＴ画像上で、該第２のＣＴ画像の他の箇所とは異なる態様で表示する、
処理を、コンピュータに実行させることを特徴とするプログラム。
【請求項４】
コンピュータによる画像表示方法であって、該コンピュータが、
異なる時期に撮影された撮影箇所が対応するＣＴ画像のうち、第１のＣＴ画像上で指定された箇所に対応する箇所を、第２の時期に撮影された第２のＣＴ画像上で、該第２のＣＴ画像の他の箇所とは異なる態様で表示する、
ことを特徴とする画像表示方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像表示装置、画像表示方法及びプログラムに関する。
【背景技術】
【０００２】
従来より、医療現場では、異なる時期に撮影したＣＴ（Computed Tomography）画像を並列に表示して、病巣部がどのように変化しているかを比較しながら診断が行われることがある。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００８−５１０２４７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、異なる時期に撮影されたＣＴ画像は、それぞれの撮影時における患者の体位の違いや、呼吸のタイミングの違い等に起因して、病巣部の位置や角度等にずれがある。このため、医師にとってはＣＴ画像間において対応箇所を把握するのが困難で、注目する病巣部の変化を比較しにくいという問題があった。
【０００５】
一つの側面では、一方のＣＴ画像で指定した病巣部が時間経過により、どのように変化したかを容易に比較できるようにすることを目的としている。
【課題を解決するための手段】
【０００６】
一態様によれば、画像表示装置は、異なる時期に撮影された撮影箇所が対応するＣＴ画像のうち、第１のＣＴ画像上で指定された箇所に対応する箇所を、第２の時期に撮影された第２のＣＴ画像上で、第２の画像の他の箇所と異なる態様で表示することを特徴とする。
【発明の効果】
【０００７】
一方のＣＴ画像で指定した病巣部が時間経過により、どのように変化したかを容易に比較できる。
【図面の簡単な説明】
【０００８】
【図１】ＣＴ画像撮影システムの一例を示す図である。
【図２】画像表示装置のハードウェア構成を示す図である。
【図３】画像表示装置における操作内容と表示内容との関係並びに診断支援プログラムの動作を示す図である。
【図４】異なる時期に撮影されたＣＴ画像の大域的位置合わせの概要を示した図である。
【図５】画像表示装置における局所的位置合わせ及び表示制御のための機能構成を示す図である。
【図６】特徴点・特徴量抽出部の機能の概要を示す図である。
【図７】特徴点・特徴量抽出部により実行される特徴点・特徴量抽出処理のフローチャートである。
【図８】特徴量ＤＢに格納される特徴量情報の一例を示す図である。
【図９】特徴点マッチング部の機能の概要を示す図である。
【図１０】特徴点マッチング部により実行される特徴点マッチング処理のフローチャートである。
【図１１】対応ベクトルＤＢに格納される対応点情報の一例を示す図である。
【図１２】対応点推定部の機能の概要を示す図である。
【図１３】対応点推定部により実行される対応点推定処理のフローチャートである。
【図１４】位置合わせ処理部により実行される位置合わせ処理のフローチャートである。
【図１５】異なる時期に撮影されたＣＴ画像の第１の表示例を示す図である。
【図１６】表示制御部により実行される第１の表示処理のフローチャートである。
【図１７】異なる時期に撮影されたＣＴ画像の第２の表示例を示す図である。
【図１８】表示制御部により実行される第２の表示処理のフローチャートである。
【図１９】表示制御部により実行される第３の表示処理のフローチャートである。
【発明を実施するための形態】
【０００９】
以下、実施形態について添付の図面を参照しながら説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複した説明を省く。
【００１０】
［第１の実施形態］
はじめに、第１の実施形態に係る画像表示装置を含むＣＴ（Computed Tomography）画像撮影システムについて説明する。図１は、ＣＴ画像撮影システムの一例を示す図である。
【００１１】
ＣＴ画像撮影システム１００は、ＣＴ装置１１０と画像表示装置１２０とを有する。ＣＴ装置１１０と画像表示装置１２０とは電気的に接続されており、両装置の間ではデータの送受信が行われる。
【００１２】
ＣＴ装置１１０は、放射線などを利用して患者の体内を走査し、コンピュータを用いて処理することで、患者の体内の断面画像を生成する（つまり、断面画像を撮影する）。ＣＴ装置１１０では、生成した断面画像（以下、ＣＴ画像と称す）を画像表示装置１２０に送信する。
【００１３】
画像表示装置１２０は、ＣＴ装置１１０において撮影されたＣＴ画像を格納するとともに、異なる時期に撮影されたＣＴ画像を位置合わせして（撮影箇所が互いに対応するようにして）並列に表示する。
【００１４】
画像表示装置１２０は、ＣＴ画像を格納する画像データベース（以下、データベースをＤＢと略す）１３０を有する。また、画像表示装置１２０は、撮影したＣＴ画像に基づいて、医師等の医療従事者が患者の診断を行う際に利用する診断支援プログラム１４０を有する。
【００１５】
診断支援プログラム１４０は、異なる時期に撮影されたＣＴ画像を並列に表示する際に、各ＣＴ画像間の位置ずれを線形変形により補正し、各ＣＴ画像間の位置合わせを行う大域的位置合わせプログラム１４１を有する。また、診断支援プログラム１４０は、大域的位置合わせが行われたＣＴ画像に残る位置ずれを非線形変形により補正し、各ＣＴ画像間の位置合わせを行う局所的位置合わせプログラム１４２を有する。更に、診断支援プログラム１４０は、局所的位置合わせが行われたＣＴ画像における注目点の指定や、注目点の拡大表示を行うための同期表示プログラム１４３を有する。
【００１６】
また、画像表示装置１２０は、局所的位置合わせプログラム１４２を実行する際に生成されるデータを格納する位置合わせＤＢ１５０を有する。
【００１７】
次に、画像表示装置１２０のハードウェア構成について説明する。図２は、画像表示装置のハードウェア構成を示す図である。図２に示すように、画像表示装置１２０は、ＣＰＵ２０１、ＲＯＭ（Read Only Memory）２０２、ＲＡＭ（Random Access Memory）２０３を備える。また、画像表示装置１２０は、記憶部２０４、接続部２０５、表示部２０６、操作部２０７を備える。なお、画像表示装置１２０の各部は、バス２０８を介して相互に接続されている。
【００１８】
ＣＰＵ２０１は、記憶部２０４に格納された各種プログラムを実行するコンピュータである。
【００１９】
ＲＯＭ２０２は不揮発性メモリである。ＲＯＭ２０２は、記憶部２０４に格納された各種プログラムをＣＰＵ２０１が実行するために必要な各種プログラム、データ等を格納する。具体的には、ＢＩＯＳ（Basic Input/Output System）やＥＦＩ（Extensible Firmware Interface）等のブートプログラムなどを格納する。
【００２０】
ＲＡＭ２０３は、ＤＲＡＭ（Dynamic Random Access Memory）やＳＲＡＭ（Static Random Access Memory）等の主記憶装置である。ＲＡＭ２０３は、記憶部２０４に各種プログラムがＣＰＵ２０１によって実行される際に展開される、作業領域として機能する。
【００２１】
記憶部２０４は、画像表示装置１２０にインストールされた各種プログラムや、プログラムを実行することで生成されるデータ等を格納する。接続部２０５は、ＣＴ装置１１０と接続され、ＣＴ装置１１０との間で、データの送受信を行う。表示部２０６は、ＣＴ装置１１０において撮影され画像ＤＢ１３０に格納されたＣＴ画像を表示する。操作部２０７は、医師等の医療従事者が画像表示装置１２０に対して行う各種操作を受け付ける。
【００２２】
次に、画像表示装置１２０において診断支援プログラム１４０を実行させた状態での、医師等の医療従事者による画像表示装置１２０に対する操作内容と表示部２０６に表示される表示内容との関係について説明する。図３は、画像表示装置における操作内容と表示内容との関係並びに診断支援プログラムの動作を示す図である。
【００２３】
画像表示装置１２０において診断支援プログラム１４０を実行させると、図３（ａ）に示すように、表示部２０６には、異なる時期に撮影されたＣＴ画像を並列して表示させるための並列表示画面が表示される。並列表示画面が表示された状態で、医師等の医療従事者により診断すべきＣＴ画像が選択されると、診断支援プログラム１４０では、画像ＤＢ１３０より、選択されたＣＴ画像を読み出す。
【００２４】
そして、診断支援プログラム１４０では、読み出したＣＴ画像を、例えば、投薬前のＣＴ画像として表示する（図３（ｂ）参照）。
【００２５】
投薬前のＣＴ画像が表示された状態で、医師等の医療従事者により投薬後のＣＴ画像が選択されると、診断支援プログラム１４０では、画像ＤＢ１３０より、選択されたＣＴ画像を読み出す。
【００２６】
このとき、診断支援プログラム１４０では、大域的位置合わせプログラム１４１を起動させ、読み出されたＣＴ画像に対して、大域的位置合わせを実行する。更に、局所的位置合わせプログラム１４２を起動させ、大域的位置合わせが実行されたＣＴ画像に対して、局所的位置合わせを実行する。
【００２７】
診断支援プログラム１４０では、局所的位置合わせが実行されたＣＴ画像を、例えば、投薬後のＣＴ画像として表示する（図３（ｃ）参照）。
【００２８】
投薬前のＣＴ画像と投薬後のＣＴ画像とが並列して表示された状態において、診断支援プログラム１４０は、同期表示プログラム１４３を実行させる。これにより、図３（ｄ）に示すように、医師等の医療従事者により一方のＣＴ画像上の任意の点が注目点として指定されると、他方のＣＴ画像上においても対応する点が注目点として指定される。また、一方のＣＴ画像において指定された注目点に対する拡大表示の指示が入力されると、他方のＣＴ画像においても指定された注目点が拡大表示される。
【００２９】
このように、画像表示装置１２０では、異なる時期に撮影されたＣＴ画像を並列して表示するにあたり、大域的位置合わせ及び局所的位置わせされたＣＴ画像を表示する。また、画像表示装置１２０では、異なる時期に撮影され並列に表示されたＣＴ画像のうち、一方のＣＴ画像に対して操作（例えば、注目点の指定）が行われると、同期して他方のＣＴ画像においても対応する位置にて同じ操作（注目点の指定）が自動的に行われる。このため、例えば、医師等の医療従事者は、ＣＴ画像間において対応箇所を容易に把握することができ、病巣部がどのように変化したのかを容易に比較することができる。
【００３０】
以下、診断支援プログラム１４０に含まれる大域的位置合わせプログラム１４１、局所的位置合わせプログラム１４２、同期表示プログラム１４３により実現される各機能の詳細について説明する。
【００３１】
はじめに、大域的位置合わせプログラム１４１が実行されることで実現される、異なる時期に撮影されたＣＴ画像の大域的位置合わせの概要について説明する。図４は、異なる時期に撮影されたＣＴ画像の大域的位置合わせの概要を示す図である。なお、図４に示す各ＣＴ画像は患者の胸部断面を患者の頭部側から見た様子を示しており、各ＣＴ画像には、患者の左側の肺４００ａ、右側の肺４００ｂ、心臓４００ｃ、脂肪４００ｄ（ＣＴ画像４０１参照）が描画されている。
【００３２】
このうち、ＣＴ画像４０１は、異なる時期に撮影されたＣＴ画像のうち、先に撮影されたＣＴ画像（例えば、術前や投薬前等に撮影されたＣＴ画像）であり、以下では、比較元ＣＴ画像と称する。
【００３３】
一方、ＣＴ画像４１１は、異なる時期に撮影されたＣＴ画像のうち、後に撮影されたＣＴ画像（例えば、術後や投薬後等に撮影されたＣＴ画像）であり、以下では、比較先ＣＴ画像と称する。
【００３４】
大域的位置合わせプログラム１４１が実行されることで、比較先ＣＴ画像４１１に対しては線形変形が施される。具体的には、比較元ＣＴ画像４０１に合わせて、比較先ＣＴ画像４１１に対して回転処理が行われ、比較先ＣＴ画像４１２が生成される。また、回転処理が行われた比較先ＣＴ画像４１２に対して平行移動処理が行われ、比較先ＣＴ画像４１３が生成される。なお、図４の例では明示していないが、必要に応じて、拡大または縮小処理が行われてもよい。
【００３５】
このように、比較先ＣＴ画像４１１を線形変形することで、比較先ＣＴ画像４１１と比較元ＣＴ画像４０１とを概ね合致させることができる。画像４２０は、比較元ＣＴ画像４０１と、大域的位置合わせが行われた比較先ＣＴ画像４１３とを重ね合わせた様子を示している。
【００３６】
画像４２０によれば、大域的位置合わせを行うことで、比較元ＣＴ画像４０１と比較先ＣＴ画像４１３とは概ね合致させることができることがわかる。一方で、大域的位置合わせだけでは、点線４２１や点線４２２に示すように、部分的な位置ずれが残る。これは比較元ＣＴ画像４０１を撮影した時と、比較先ＣＴ画像４１１を撮影した時とでは、患者の呼吸や心拍のタイミング等に違いがあるからである。
【００３７】
そこで、画像表示装置１２０では、局所的位置合わせプログラム１４２を実行させることで、大域的位置合わせを行った比較先ＣＴ画像４１３に対し、局所的位置合わせを施す。更に、局所的位置合わせの結果を利用して、同期表示プログラム１４３を実行させることで、比較元ＣＴ画像４０１と局所的位置合わせが行われた比較先ＣＴ画像４１３とを比較しやすくすべく表示制御を行う。以下、画像表示装置１２０による局所的位置合わせ及び表示制御の詳細について説明する。
【００３８】
はじめに、画像表示装置１２０による局所的位置合わせ及び表示制御のための機能構成について説明する。図５は、画像表示装置における局所的位置合わせ及び表示制御のための機能構成を示す図である。
【００３９】
画像表示装置１２０は、ＣＰＵ２０１が局所的位置合わせプログラム１４２及び同期表示プログラム１４３を実行することで、以下に説明する各部の機能を実現する。なお、本実施形態において、画像表示装置１２０の位置合わせＤＢ１５０には、特徴量ＤＢ５２１、対応ベクトルＤＢ５２２が含まれており、各ＤＢは、例えば記憶部２０４に設けられているものとする。
【００４０】
図５に示すように、画像表示装置１２０は、表示画像入力部５１１、特徴点・特徴量抽出部５１２、特徴点マッチング部５１３、対応点推定部５１４、位置合わせ処理部５１５を有する。また、表示制御部５１６を有する。
【００４１】
表示画像入力部５１１は、医師等の医療従事者による指定に応じて、比較元ＣＴ画像４０１と、大域的位置合わせ処理が施された比較先ＣＴ画像４１３とを入力する。
【００４２】
特徴点・特徴量抽出部５１２は、入力された比較元ＣＴ画像４０１に含まれる複数の特徴点について、特徴量を算出する。同様に、特徴点・特徴量抽出部５１２は、入力された比較先ＣＴ画像４１３に含まれる複数の特徴点について、特徴量を算出する。
【００４３】
特徴点マッチング部５１３は、比較元ＣＴ画像４０１の各特徴点の特徴量、比較先ＣＴ画像４１３の各特徴点の特徴量とに基づいて、比較元ＣＴ画像４０１の各特徴点に対応する比較先ＣＴ画像４１３の各特徴点（対応点）を決定する。また、特徴点マッチング部５１３は、決定した対応点を対応点情報として対応ベクトルＤＢ５２２に格納する。
【００４４】
対応点推定部５１４は、特徴点マッチング部５１３において対応点を決定することができなかった、比較元ＣＴ画像４０１の各特徴点について、比較先ＣＴ画像４１３における対応点を推定する。また、対応点推定部５１４は、推定した対応点を対応点情報に追加したうえで、対応点情報を対応ベクトルＤＢ５２２に格納する。
【００４５】
位置合わせ処理部５１５は、対応ベクトルＤＢ５２２に格納された対応点情報に基づいて、大域的位置合わせ処理が施された比較先ＣＴ画像４１３に、局所的位置合わせ処理を行う。
【００４６】
表示制御部５１６は、比較元ＣＴ画像４０１と、局所的位置合わせ処理が行われた比較先ＣＴ画像４１３とが並列に表示された状態で、ＣＴ画像内の一部が指定された場合に、他方のＣＴ画像内の対応する位置が自動的に指定されるように表示を制御する。更に、表示制御部５１６は、並列して表示されたＣＴ画像内の一部について、医師等の医療従事者が拡大表示の指示を入力した場合に、他方のＣＴ画像内の対応する位置が自動的に拡大表示されるように表示を制御する。
【００４７】
次に、画像表示装置１２０による局所的位置合わせのための各機能の詳細について説明する。画像表示装置１２０による局所的位置合わせのための各機能のうち、図６は、特徴点・特徴量抽出部の機能の概要を示す図である。
【００４８】
図６に示すように、特徴点・特徴量抽出部５１２は、比較元ＣＴ画像４０１を縦方向の複数の格子線及び横方向の複数の格子線により分け、各格子線の交点を特徴点として決定する。図６の例では、決定された特徴点を、それぞれα１、α２、α３・・・としている。
【００４９】
また、特徴点・特徴量抽出部５１２は、比較元ＣＴ画像４０１の各特徴点について、その周辺の所定画素分の領域（例えば、半径１０画素分の領域）について、その領域内に含まれる各画素の輝度値に基づいて、特徴量を算出する。
【００５０】
具体的には、特徴点・特徴量抽出部５１２では、特徴点α１周辺の領域６０１内に含まれる任意の２つの画素を抽出し、それぞれの画素の輝度値の差分を算出する。そして、算出結果が所定の閾値より大きければ"１"を、所定の閾値以下であれば"０"を導出し、特徴点α１の特徴量とする。
【００５１】
図６の例は、任意の２つの画素として、Ａ１点とＢ１点とが抽出され、両者の輝度差が閾値より大きかったため、特徴点α１の特徴量に"１"が設定されている。更に、任意の２つの画素として、Ａ２点とＢ２点とが抽出され、両者の輝度差が閾値より小さかったため、特徴点α１の特徴量に"０"が設定されている。なお、１つの特徴点の周辺の領域より２つの画素を抽出する処理は、所定回数（例えば、１２８回）繰り返すものとする。つまり、１つの特徴点の特徴量には、所定回数分の桁数（例えば、１２８桁）が含まれる。
【００５２】
同様に、特徴点・特徴量抽出部５１２は、比較先ＣＴ画像４１３を縦方向の複数の格子線及び横方向の複数の格子線により分け、各格子線の交点を特徴点として決定する。図６の例では、決定された特徴点を、それぞれβ１、β２、β３・・・としている。
【００５３】
また、特徴点・特徴量抽出部５１２は、比較先ＣＴ画像４１３の各特徴点について、その周辺の所定画素分の領域（例えば、半径１０画素分の領域）について、その領域内に含まれる各画素の輝度値に基づいて、特徴量を算出する。
【００５４】
具体的には、特徴点・特徴量抽出部５１２では、特徴点β１周辺の領域６１１内に含まれる任意の２つの画素を抽出し、それぞれの画素の輝度値の差分を算出する。そして、算出結果が所定の閾値より大きければ"１"を、所定の閾値以下であれば"０"を導出し、特徴点β１の特徴量とする。
【００５５】
図６の例は、任意の２つの画素として、Ａ１点とＢ１点とが抽出され、両者の輝度差が閾値以下であったため、特徴点β１の特徴量に"０"が設定されている。更に、任意の２つの画素として、Ａ２点とＢ２点とが抽出され、両者の輝度差が閾値より大きかったため、特徴点β１の特徴量に"１"が設定されている。なお、比較先ＣＴ画像４１３についても、１つの特徴点の周辺の領域より２つの画素を抽出する処理は、所定回数（例えば、１２８回）繰り返すものとする。つまり、１つの特徴点の特徴量には、所定回数分の桁数（例えば、１２８桁）が含まれる。
【００５６】
次に、特徴点・特徴量抽出部５１２により実行される特徴点・特徴量抽出処理の流れについて説明する。図７は、特徴点・特徴量抽出部により実行される特徴点・特徴量抽出処理のフローチャートである。なお、ここでは説明の簡略化のため、比較元ＣＴ画像４０１について特徴点・特徴量抽出処理を行う場合について説明するが、比較先ＣＴ画像４１３についても同様である。
【００５７】
図７に示すように、ステップＳ７０１において、特徴点・特徴量抽出部５１２は、比較元ＣＴ画像４０１を、複数の縦方向の格子線と複数の横方向の格子線とを用いて分ける。また、特徴点・特徴量抽出部５１２は、各格子線の交点を特徴点として決定する。
【００５８】
ステップＳ７０２において、特徴点・特徴量抽出部５１２は、特徴点をカウントするカウンタ値ｍに１を代入する。
【００５９】
ステップＳ７０３において、特徴点・特徴量抽出部５１２は、特徴量をカウントするカウンタ値ｎに１を代入する。
【００６０】
ステップＳ７０４において、特徴点・特徴量抽出部５１２は、特徴点αｍの周辺の領域６０１において、Ａｎ点とＢｎ点とを抽出し、抽出したＡｎ点の輝度値とＢｎ点の輝度値の差分（輝度差）を算出する。
【００６１】
ステップＳ７０５において、特徴点・特徴量抽出部５１２は、ステップＳ７０４において算出した輝度差が所定の閾値より大きいか否かを判定する。ステップＳ７０５において、所定の閾値より大きいと判定された場合には、ステップＳ７０６に進み、特徴点αｍの特徴量に"１"を追加する。一方、ステップＳ７０５において、所定の閾値以下であると判定された場合には、ステップＳ７０７に進み、特徴点αｍの特徴量に"０"を追加する。
【００６２】
ステップＳ７０８において、特徴点・特徴量抽出部５１２は、所定回数分（例えば、１２８回分）の輝度差を算出したか否かを判定する。ステップＳ７０８において、所定回数分の輝度差を算出していないと判定した場合には、ステップＳ７０９に進む。ステップＳ７０９において、特徴点・特徴量抽出部５１２は、特徴量をカウントするカウンタ値ｎをインクリメントした後、ステップＳ７０４に戻る。
【００６３】
一方、ステップＳ７０８において、所定回数分の輝度差を算出したと判定した場合には、ステップＳ７１０に進む。ステップＳ７１０において、特徴点・特徴量抽出部５１２は、特徴点αｍについて算出した特徴量に含まれる"１"の個数を算出する。算出の結果、特徴量に含まれる"１"の個数が所定の閾値以上あった場合、特徴点αｍが輝度差の大きい特徴点であると判定し、ステップＳ７１１に進む。
【００６４】
なお、ＣＴ画像において輝度差の大きい特徴点とは、例えば、患者の体表面や、右側及び左側の肺の辺縁部分、あるいは心臓の辺縁部分等に位置する特徴点である。
【００６５】
ステップＳ７１１において、特徴点・特徴量抽出部５１２は、特徴点αｍの特徴量を特徴量ＤＢ５２１に格納する。
【００６６】
一方、特徴量に含まれる"１"の個数が所定の閾値未満であった場合、特徴点・特徴量抽出部５１２では、特徴点αｍが輝度差の小さい特徴点であると判定し、ステップＳ７１２に進む。
【００６７】
なお、ＣＴ画像において輝度差の小さい特徴点とは、例えば、患者の右側及び左側の肺の内部や、心臓の内部等に位置する特徴点である。
【００６８】
ステップＳ７１２において、特徴点・特徴量抽出部５１２は、特徴点αｍの特徴量はなしと判定する。この場合、特徴点αｍの特徴量は特徴量ＤＢ５２１に格納されない。
【００６９】
ステップＳ７１３において、特徴点・特徴量抽出部５１２は、比較元ＣＴ画像４０１に含まれるすべての特徴点について、特徴量を算出したか否かを判定する。ステップＳ７１３において、特徴量を算出していない特徴点があると判定した場合には、ステップＳ７１４に進む。ステップＳ７１４において、特徴点・特徴量抽出部５１２は、特徴点をカウントするカウンタ値ｍをインクリメントした後、ステップＳ７０３に戻る。
【００７０】
一方、ステップＳ７１４において、比較元ＣＴ画像４０１に含まれるすべての特徴点について、特徴量を算出したと判定した場合には、特徴点・特徴量抽出処理を終了する。
【００７１】
次に、特徴量ＤＢ５２１に格納される特徴量情報について説明する。図８は特徴量ＤＢに格納される特徴量情報の一例を示す図である。
【００７２】
図８（ａ）、（ｂ）に示すように、特徴量ＤＢ５２１に格納される特徴量情報は、比較元ＣＴ画像４０１と比較先ＣＴ画像４１３とで分けて管理される。比較元ＣＴ画像４０１の特徴量情報８００は、情報の項目として、"特徴点"と"特徴量"とを含む。
【００７３】
"特徴点"には、比較元ＣＴ画像４０１より抽出された特徴点が格納される。"特徴量"には、各特徴点の特徴量が格納される。
【００７４】
同様に、比較先ＣＴ画像４１３の特徴量情報８１０は、情報の項目として、"特徴点"と"特徴量"とを含む。"特徴点"には、比較先ＣＴ画像４１３より抽出された特徴点が格納される。"特徴量"には、各特徴点の特徴量が格納される。
【００７５】
なお、特徴点・特徴量抽出処理において、特徴量なしと判定された特徴点（図８の例では、特徴点α３と特徴点βＭ−１）については、特徴量は格納されない。
【００７６】
次に、画像表示装置１２０による局所的位置合わせのための各機能のうち、特徴点マッチング部５１３の機能の概要について説明する。図９は、特徴点マッチング部の機能の概要を示す図である。
【００７７】
図９に示すように、特徴点マッチング部５１３では、比較元ＣＴ画像４０１の各特徴点に対応する、比較先ＣＴ画像４１３の特徴点を抽出するために、各特徴点の特徴量について排他的論理和を算出する。
【００７８】
具体的には、比較元ＣＴ画像４０１の特徴点α１の特徴量と、比較先ＣＴ画像４１３の特徴点β１の特徴量、特徴点β２の特徴量、特徴点β３の特徴量、・・・との排他的論理和をそれぞれ算出する。なお、図９の例では、比較元ＣＴ画像４０１の特徴点α１の特徴量との排他的論理和についてのみ示しているが、特徴点マッチング部５１３では、比較元ＣＴ画像４０１に含まれるすべての特徴点について、排他的論理和を算出する。
【００７９】
また、特徴点マッチング部５１３では、算出された排他的論理和に基づいて、各特徴点間の類似度を算出する。例えば、算出された排他的論理和に含まれる"０"の数に応じて類似度を算出する。
【００８０】
更に、特徴点マッチング部５１３では、各特徴量間において算出した類似度のうち、最大類似度の特徴点を「対応点」に決定する。例えば、比較元ＣＴ画像４０１の特徴点α１の特徴量と、比較先ＣＴ画像４１３の特徴点β１の特徴量との類似度が最大であった場合、比較元ＣＴ画像４０１の特徴点α１の対応点は、比較先ＣＴ画像４１３の特徴点β１に決定される。
【００８１】
次に、特徴点マッチング部５１３により実行される特徴点マッチング処理の流れについて説明する。図１０は、特徴点マッチング部により実行される特徴点マッチング処理のフローチャートである。
【００８２】
ステップＳ１００１において、特徴点マッチング部５１３は、比較元ＣＴ画像４０１の特徴点をカウントするカウンタ値ｍａに１を代入する。
【００８３】
ステップＳ１００２において、特徴点マッチング部５１３は、比較先ＣＴ画像４１３の特徴点をカウントするカウンタ値ｍｂに１を代入する。
【００８４】
ステップＳ１００３において、特徴点マッチング部５１３は、特徴量ＤＢ５２１より特徴点αｍａの特徴量を読み出す。ステップＳ１００４において、特徴点マッチング部５１３は、特徴点αｍａの特徴量があるか否かを判定する。
【００８５】
ステップＳ１００４において特徴点αｍａの特徴量がないと判定された場合には、ステップＳ１０１２に進み、特徴点αｍａの対応点はなしと判定する。一方、ステップＳ１００４において特徴点αｍａの特徴量があると判定された場合には、ステップＳ１００５に進む。
【００８６】
ステップＳ１００５において、特徴点マッチング部５１３は、特徴量ＤＢ５２１より特徴点βｍｂの特徴量を読み出す。ステップＳ１００６において、特徴点マッチング部５１３は、特徴点βｍｂの特徴量があるか否かを判定する。
【００８７】
ステップＳ１００６において特徴点βｍｂの特徴量がないと判定された場合には、ステップＳ１００９に進む。一方、ステップＳ１００６において特徴点βｍｂの特徴量があると判定された場合には、ステップＳ１００７に進む。
【００８８】
ステップＳ１００７において、特徴点マッチング部５１３は、特徴点αｍａの特徴量と特徴点βｍｂの特徴量との排他的論理和を算出する。ステップＳ１００８において、特徴点マッチング部５１３は、ステップＳ１００５における算出結果に基づいて、特徴点αｍａの特徴量と特徴点βｍｂの特徴量との類似度を算出する。
【００８９】
ステップＳ１００９において、特徴点マッチング部５１３は、比較先ＣＴ画像４１３のすべての特徴点について、特徴点αｍａとの排他的論理和の算出を行ったか否かを判定する。ステップＳ１００９において、比較先ＣＴ画像４１３の特徴点の中に、特徴点αｍａとの排他的論理和の算出を行っていない特徴点があると判定した場合には、ステップＳ１０１０に進む。ステップＳ１０１０において、特徴点マッチング部５１３は、カウンタ値ｍａをインクリメントし、ステップＳ１００５に戻る。
【００９０】
一方、ステップＳ１００９において、比較先ＣＴ画像４１３のすべての特徴点について、特徴点αｍａとの排他的論理和の算出を完了したと判定した場合には、ステップＳ１０１１に進む。
【００９１】
ステップＳ１０１１において、特徴点マッチング部５１３は、算出した全ての類似度が閾値未満か否かを判定する。ステップＳ１０１１において、算出した全ての類似度が閾値未満であると判定した場合には、ステップＳ１０１２に進む。ステップＳ１０１２において、特徴点マッチング部５１３は、比較元ＣＴ画像４０１の特徴点αｍａに対応する、比較先ＣＴ画像４１３の特徴点はなかったと判断し、特徴点αｍａの対応点はなしと判定する。
【００９２】
一方、ステップＳ１０１１において、算出した類似度の中に閾値以上の類似度があると判定した場合には、ステップＳ１０１３に進む。ステップＳ１０１３において、特徴点マッチング部５１３は、比較先ＣＴ画像の各特徴点のうち、最大類似度の特徴点を、比較元ＣＴ画像の特徴点αｍａの対応点に決定する。
【００９３】
ステップＳ１０１４において、特徴点マッチング部５１３は、ステップＳ１０１３において決定された対応点を、対応ベクトルＤＢ５２２に格納する。なお、ステップＳ１０１２において、対応点なしと判定された場合には、特徴点αｍａに対応する対応点はない旨の情報を、対応ベクトルＤＢ５２２に格納する。
【００９４】
ステップＳ１０１５において、特徴点マッチング部５１３は、比較元ＣＴ画像４０１のすべての特徴点について特徴点マッチング処理を実行したか否かを判定する。ステップＳ１０１５において、比較元ＣＴ画像４０１の特徴点の中に、特徴点マッチング処理を実行していない特徴点があると判定した場合には、ステップＳ１０１６に進む。ステップＳ１０１６において、特徴点マッチング部５１３は、カウンタ値ｍａをインクリメントした後、ステップＳ１００２に戻る。
【００９５】
一方、ステップＳ１０１５において、比較元ＣＴ画像４０１に含まれるすべての特徴点について特徴点マッチング処理を実行したと判定した場合には、特徴点マッチング処理を終了する。
【００９６】
このように、本実施形態における画像表示装置１２０では、特徴点マッチング処理において、特徴点αｍａの特徴量がない場合に、対応点を決定するための処理を実行することなく、対応点なしと判定する。また、特徴点βｍｂの特徴量がない場合には、類似度の算出を行わない。
【００９７】
この結果、画像表示装置１２０では、局所的位置合わせプログラム１４２を実行した際に、少ない計算量で特徴点マッチング処理を行うことができる。
【００９８】
次に、特徴点マッチング処理により対応ベクトルＤＢ５２２に格納された対応点情報について説明する。図１１は、対応ベクトルＤＢに格納される対応点情報の一例を示す図である。
【００９９】
図１１に示すように、対応点情報１１００は、情報の項目として"比較元ＣＴ画像の特徴点"と"比較先ＣＴ画像の特徴点"とが含まれる。"比較元ＣＴ画像の特徴点"には、比較元ＣＴ画像４０１の特徴点が格納される。"比較先ＣＴ画像の特徴点"には、比較先ＣＴ画像４１３の各特徴点のうち、比較元ＣＴ画像４０１の特徴点の対応点として決定された特徴点が格納される。なお、図１１の例では、比較元ＣＴ画像４０１の特徴点α３、αＭ−１は、比較先ＣＴ画像４１３内に対応する特徴点がなく、対応点なしと判定されたことを示している。
【０１００】
次に、画像表示装置１２０による局所的位置合わせのための各機能のうち、対応点推定部５１４の機能の概要について説明する。図１２は、対応点推定部の機能の概要を示す図である。
【０１０１】
対応点推定部５１４では、比較元ＣＴ画像４０１の特徴点のうち、特徴点マッチング部５１３において対応点がないと判定された特徴点に対応する特徴点を推定する。
【０１０２】
図１２において、拡大表示部１２０１、１２１１は、それぞれ、説明をわかりやすくするために比較元ＣＴ画像４０１及び比較先ＣＴ画像４１３の一部を便宜的に拡大して示した図である。拡大表示部１２０１において、特徴点αｘは、患者の肺の内部を示す点であり、特徴点マッチング部５１３により対応点なしと判定された特徴点である。
【０１０３】
境界線１２０２、１２１２は、患者の肺の内側と外側との境界を示す辺縁部分である。特徴点α１１は、境界線１２０２上の点であり、特徴点マッチング部５１３により、境界線１２１２上の特徴点β１１（拡大表示部１２１１参照）が対応点として決定されている点である。
【０１０４】
同様に、境界線１２０３、１２１３は、患者の肺の内側と外側の境界を示す辺縁部分である。特徴点α１２は、境界線１２０３上の点であり、特徴点マッチング部５１３により、境界線１２１３上の特徴点β１２（拡大表示部１２１１参照）が対応点として決定されている点である。
【０１０５】
ここで、特徴点α１１とα１２とは、特徴点αｘを通る直線１２０４と境界線１２０２または１２０３との交点として特定され、特徴点αｘを挟んで互いに対向した位置にある。また、特徴点αｘと特徴点α１１との間の距離はｒ１、特徴点αｘと特徴点α１２との間の距離はｒ２であるとする。
【０１０６】
対応点推定部５１４では、特徴点α１１、α１２と、特徴点αｘとの間のこれらの関係に基づいて、特徴点αｘの対応点である特徴点βｘを推定する。
【０１０７】
具体的には、特徴点β１１とβ１２とを通る直線１２１４を求め、直線１２１４上の点であって、特徴点β１１からの距離と、特徴点β１２からの距離との比が、ｒ１：ｒ２の関係となる特徴点βｘを、特徴点αｘの対応点として推定する。
【０１０８】
次に、対応点推定部５１４により実行される対応点推定処理の流れについて説明する。図１３は、対応点推定部により実行される対応点推定処理のフローチャートである。
【０１０９】
ステップＳ１３０１において、対応点推定部５１４は、対応ベクトルＤＢ５２２より、対応点なしと判定された比較元ＣＴ画像４０１の特徴点を処理対象の特徴点として、順次読み出す。
【０１１０】
ステップＳ１３０２において、対応点推定部５１４は、処理対象の特徴点（αｘ）を通る直線上の点であって、処理対象の特徴点を挟んで対向する最近傍の特徴点（対応点が決定されている辺縁部分の特徴点）を、参照点として対応ベクトルＤＢ５２２より読み出す。
【０１１１】
ステップＳ１３０３において、対応点推定部５１４は、処理対象の特徴点からステップＳ１３０２において読み出した各参照点（α１１、α１２）までの距離の比（ｒ１：ｒ２）を算出する。
【０１１２】
ステップＳ１３０４において、対応点推定部５１４は、各参照点（α１１、α１２）に対応する比較先ＣＴ画像４１３の特徴点を、対応ベクトルＤＢ５２２より探索し特定する。
【０１１３】
ステップＳ１３０５において、対応点推定部５１４は、ステップＳ１３０４において探索された、各参照点（α１１、α１２）に対応する特徴点（β１１、β１２）を通る直線を求める。また、対応点推定部５１４は、求めた直線上の点であって、各参照点に対応する特徴点からの距離の比が、ステップＳ１３０３において算出された距離の比（ｒ１：ｒ２）と等しくなる点を算出する。更に、対応点推定部５１４は、算出した点（βｘ）を、処理対象の特徴点の対応点に決定する。
【０１１４】
ステップＳ１３０６において、対応点推定部５１４は、ステップＳ１３０５において決定された対応点を、処理対象の特徴点の対応点として、対応ベクトルＤＢ５２２に格納する。
【０１１５】
ステップＳ１３０７において、対応点推定部５１４は、対応点なしと判定された比較元ＣＴ画像４０１の特徴点すべてについて、対応点推定処理を実行したか否かを判定する。ステップＳ１３０７において、対応点推定処理を実行していない特徴点があると判定した場合には、ステップＳ１３０１に進み、対応点なしと判定された次の特徴点を処理対象として、ステップＳ１３０１からステップＳ１３０６までの処理を行う。
【０１１６】
一方、ステップＳ１３０７において、対応点なしと判定されたすべての特徴点について、対応点推定処理を実行したと判定した場合には、対応点推定処理を終了する。
【０１１７】
このように、本実施形態における画像表示装置１２０では、２つの参照点からの距離の比に基づいて対応点を推定する。このため、特徴点マッチング処理により各特徴点の特徴量との類似度を算出したうえで対応点を決定する場合と比較して、少ない計算量で対応点を推定できる。
【０１１８】
つまり、本実施形態における画像表示装置１２０では、
・患者の体表面や、臓器（肺及び心臓）の辺縁部分のように、輝度差の大きい特徴点については、特徴点マッチング処理により高い精度で対応点を決定する。
・患者の肺及び心臓の内部のように、輝度差の小さい特徴点については、特徴点マッチング処理により高い精度で決定した対応点を用いて、距離の比に基づいて対応点を推定することで、計算量を削減する。
【０１１９】
これにより、本実施形態における画像表示装置１２０によれば、精度の高い位置合わせを少ない計算量で実現することができる。
【０１２０】
次に、画像表示装置１２０により局所的位置合わせのための各機能のうち、位置合わせ処理部５１５の機能の概要について説明する。図１４は、位置合わせ処理部５１５により実行される位置合わせ処理のフローチャートである。
【０１２１】
ステップＳ１４０１において、位置合わせ処理部５１５は、対応ベクトルＤＢ５２２に格納された比較元ＣＴ画像４０１の特徴点と、比較元ＣＴ画像４０１の特徴点に対応する、比較先ＣＴ画像４１３の対応点とを読み出す。
【０１２２】
ステップＳ１４０２において、位置合わせ処理部５１５は、ステップ１４０１において読み出された、特徴点とその対応点とにより規定される対応ベクトルを算出する。また、位置合わせ処理部５１５は、算出した対応ベクトルに基づいて、比較先ＣＴ画像４１３の各特徴点の移動先を決定する。このとき、位置合わせ処理部５１５では、特徴点を移動した後の比較先ＣＴ画像４１３が滑らかな画像になるように、補正処理を行う（つまり、対応ベクトルを補正する）。
【０１２３】
ステップＳ１４０３において、位置合わせ処理部５１５は、補正後の対応ベクトルに基づき、比較先ＣＴ画像４１３の各画素の位置合わせを行う。ステップＳ１４０４において、位置合わせ処理部５１５は、位置合わせした比較先ＣＴ画像を表示部２０６に表示する。
【０１２４】
次に、画像表示装置１２０による表示制御のための機能である、表示制御部５１６の機能の概要について説明する。表示制御部５１６では、並列して表示された比較元ＣＴ画像４０１と、位置合わせ処理部５１５により位置合わせが行われた比較先ＣＴ画像とを比較しやするするための表示制御を行う。
【０１２５】
図１５は、異なる時期に撮影されたＣＴ画像の第１の表示例であり、図１５（ａ）、（ｃ）は比較元ＣＴ画像４０１を、図１５（ｂ）、（ｄ）は、位置合わせ処理部５１５により位置合わせが行われた比較先ＣＴ画像１５０１をそれぞれ示している。
【０１２６】
ここで、図１５（ａ）は、表示された比較元ＣＴ画像４０１上において、座標（ｘ１，ｙ１）の点１５１１が、注目点として医師等の医療従事者により指示された様子を示している。表示制御部５１６では、表示された比較元ＣＴ画像４０１上において点１５１１が指示されたことを検知すると、指示された点１５１１の座標を識別する。更に、表示制御部５１６では、指示された点１５１１に対応する比較先ＣＴ画像１５０１上の点１５１２を、注目点として比較先ＣＴ画像１５０１上に自動的に表示する。
【０１２７】
なお、上述したとおり、位置合わせ処理部５１５による位置合わせ処理では、対応ベクトルの補正が行われている。このため、比較元ＣＴ画像４０１上の点１５１１に対応する比較先ＣＴ画像１５０１上の点１５１２は、比較元ＣＴ画像４０１上の点１５１１の座標と完全に同じ座標にはならない。このため、表示制御部５１６では、対応ベクトルの補正を加味した比較先ＣＴ画像１５０１上の座標（ｘ１'，ｙ１'）を指示することとなる。
【０１２８】
図１５（ｃ）は、医師等の医療従事者が、指示した点１５１１を含む領域１５２１を拡大表示する指示を入力した様子を示している。表示制御部５１６では、領域１５２１を拡大表示する指示が入力されたことを検知すると、領域１５２１に対応する比較先ＣＴ画像１５０１上の領域１５２２を、比較先ＣＴ画像１５０１上において自動的に拡大表示する。
【０１２９】
このように、表示制御部５１６では、異なる時期に撮影された比較元ＣＴ画像４０１と比較先ＣＴ画像１５０１とが並列に表示された状態で、一方のＣＴ画像上の点が注目点として指示されると他方のＣＴ画像上の点を注目点として自動的に指示する。また、一方のＣＴ画像上の点を含む領域について拡大表示の指示が入力されると、他方のＣＴ画像上の対応する点を含む領域を自動的に拡大表示する。このため、医師等の医療従事者は、注目点についての対応箇所を容易に把握することができ、病巣部がどのように変化しているかを容易に比較することができる。
【０１３０】
次に、表示制御部５１６により実行される表示処理の流れについて説明する。図１６は、表示制御部により実行される第１の表示処理のフローチャートである。なお、図１６に示すフローチャートの実行に際して、表示部２０６には、比較元ＣＴ画像４０１と比較先ＣＴ画像１５０１が並列に表示されているものとする。
【０１３１】
ステップＳ１６０１において、表示制御部５１６は、医師等の医療従事者により注目点の指定があるか否かを判定する。ステップＳ１６０１において、注目点の指定がないと判定された場合には、注目点の指定があると判定されるまで待機する。
【０１３２】
一方、ステップＳ１６０１において、注目点の指定があると判定された場合には、ステップＳ１６０２に進む。ステップＳ１６０２において、表示制御部５１６は、注目点の指定が、比較元ＣＴ画像上において行われたか否かを判定する。
【０１３３】
ステップＳ１６０２において、注目点の指定が、比較元ＣＴ画像上において行われたと判定された場合には、ステップＳ１６０３に進む。ステップＳ１６０３において、表示制御部５１６は、指定された注目点の座標を識別する。
【０１３４】
ステップＳ１６０４において、表示制御部５１６は、比較元ＣＴ画像４０１において識別した座標の表示態様を変更する。
【０１３５】
ステップＳ１６０５において、表示制御部５１６は、識別した座標の対応点を対応ベクトルＤＢ５２２より読み出し、読み出した対応点に対して位置合わせ処理部５１５による位置合わせ処理が行われた後の座標を算出する。更に、表示制御部５１６は、比較先ＣＴ画像１５０１において、算出した座標の表示態様を変更する。
【０１３６】
ステップＳ１６０６において、表示制御部５１６は、医師等の医療従事者により注目点の拡大表示の指示が入力されたか否かを判定する。ステップＳ１６０６において、注目点の拡大表示の指示が入力されなかったと判定した場合には、ステップＳ１６０１に戻る。
【０１３７】
一方、ステップＳ１６０６において、注目点の拡大表示の指示が入力されたと判定された場合には、ステップＳ１６０７に進む。ステップＳ１６０７において、表示制御部５１６は、比較元ＣＴ画像４０１において、識別した座標を含む所定領域を拡大表示する。
【０１３８】
ステップＳ１６０８において、表示制御部５１６は、比較先ＣＴ画像１５０１において、算出した座標を含む所定領域を拡大表示する。
【０１３９】
一方、ステップＳ１６０２において、注目点の指定が、比較先ＣＴ画像上において行われたと判定された場合には、ステップＳ１６０９に進む。ステップＳ１６０９において、表示制御部５１６は、指定された注目点の座標を識別する。
【０１４０】
ステップＳ１６１０において、表示制御部５１６は、比較先ＣＴ画像１５０１において識別した座標の表示態様を変更する。
【０１４１】
ステップＳ１６１１において、表示制御部５１６は、識別した座標についての位置合わせ処理部５１５による位置合わせ処理前の座標を算出し、算出した座標に対する比較元ＣＴ画像４０１上の対応点の座標を対応ベクトルＤＢ５２２より抽出する。更に、表示制御部５１６は、比較元ＣＴ画像４０１において、抽出した座標の表示態様を変更する。
【０１４２】
ステップＳ１６１２において、表示制御部５１６は、医師等の医療従事者により注目点の拡大表示の指示が入力されたか否かを判定する。ステップＳ１６１２において、注目点の拡大表示の指示が入力されなかったと判定した場合には、ステップＳ１６０１に戻る。
【０１４３】
一方、ステップＳ１６１２において、注目点の拡大表示の指示が入力されたと判定された場合には、ステップＳ１６１３に進む。ステップＳ１６１３において、表示制御部５１６は、比較先ＣＴ画像１５０１において、識別した座標を含む所定領域を拡大表示する。
【０１４４】
ステップＳ１６１４において、表示制御部５１６は、比較元ＣＴ画像４０１において、抽出した座標を含む所定領域を拡大表示する。
【０１４５】
ステップＳ１６１５において、表示制御部５１６は、表示処理の終了指示が入力されたか否かを判定し、入力されていないと判定した場合には、ステップＳ１６０１に戻る。一方、表示処理の終了指示が入力されたと判定した場合には、表示処理を終了する。
【０１４６】
以上の説明から明らかなように、本実施形態に係る画像表示装置１２０では、
・異なる時期に撮影したＣＴ画像を、比較元ＣＴ画像と比較先ＣＴ画像として並列して表示するにあたり、比較先ＣＴ画像に対して、大域的位置合わせ及び局所的位置合わせを行う。
・局所的位置合わせを行うにあたり、比較元ＣＴ画像の特徴点に対応する比較先ＣＴ画像の特徴点を対応点として算出する。
・比較元ＣＴ画像の特徴点の対応点がなかった場合、比較元ＣＴ画像の特徴点近傍の辺縁部分の特徴点であって、対向する位置にある特徴点を用いて、対応点を推定する。
・算出した対応点と推定した対応点とを用いて比較先ＣＴ画像の局所的位置合わせを行う。
・並列して表示した比較元ＣＴ画像または比較先ＣＴ画像の一方において、注目点が指定された場合、他方においても対応する位置の注目点を自動的に指定する。また、注目点を含む領域の拡大表示が指示されると他方においても対応する点を含む領域を自動的に拡大表示する。
【０１４７】
これにより、一方のＣＴ画像で指定した病巣部が時間経過により、どのように変化したかを容易に比較することができる。また、少ない計算量で局所的位置合わせを実現することができる。
【０１４８】
［第２の実施形態］
第２の実施形態において、表示制御部５１６は、比較元ＣＴ画像と比較先ＣＴ画像とを並列に表示するにあたり、比較先ＣＴ画像として、大域的位置合わせ及び局所的位置合わせを行っていないＣＴ画像を用いる。
【０１４９】
図１７は、異なる時期に撮影されたＣＴ画像の第２の表示例であり、図１７（ａ）、（ｃ）は比較元ＣＴ画像４０１を、図１７（ｂ）、（ｄ）は、大域的位置合わせ及び局所的位置合わせが行われる前の比較先ＣＴ画像４１１をそれぞれ示している。
【０１５０】
ここで、図１７（ａ）は、表示された比較元ＣＴ画像４０１上において、座標（ｘ１，ｙ１）の点１５１１が、医師等の医療従事者により指示された様子を示している。表示制御部５１６では、表示された比較元ＣＴ画像４０１上において、点１５１１が指示されたことを検知すると、指示された点１５１１の座標を識別する。更に、表示制御部５１６では、指示された点１５１１に対応する比較先ＣＴ画像４１１上の点１７１２を比較先ＣＴ画像４１１上に自動的に表示する。
【０１５１】
なお、比較先ＣＴ画像４１１は、大域的位置合わせ及び局所的位置合わせが行われていない画像であるため、表示制御部５１６では、対応ベクトルＤＢ５２２を参照することにより、点１５１１に対応する対応点を導出する。対応ベクトルＤＢ５２２を参照することにより導出される対応点は、大域的位置合わせが行われた比較先ＣＴ画像４１３上の点であるため、表示制御部５１６では、更に、大域的位置合わせが行われる前の点に変換する。このようにして求めた点１７１２（座標（ｘ２，ｙ２））を比較先ＣＴ画像４１１上に表示する。つまり、表示制御部５１６では、比較先ＣＴ画像４１１上の座標（ｘ２，ｙ２）を自動的に指示する。
【０１５２】
図１７（ｃ）は、医師等の医療従事者が、指示した点１５１１を含む領域１５２１を拡大表示する指示を入力した様子を示している。表示制御部５１６では、領域１５２１を拡大表示する指示が入力されたことを検知すると、領域１５２１に対応する比較先ＣＴ画像４１１上の領域１７２２を、比較先ＣＴ画像４１１上において自動的に拡大表示する。
【０１５３】
このように、表示制御部５１６では、異なる時期に撮影された比較元ＣＴ画像４０１と比較先ＣＴ画像４１１とが並列に表示された状態で、一方のＣＴ画像上の点が指示されると他方のＣＴ画像上の点を自動的に指示する。また、一方のＣＴ画像上の点を含む領域について拡大表示の指示が入力されると、他方のＣＴ画像上の点を含む領域を自動的に拡大表示する。このため、医師等の医療従事者は、注目点についての対応箇所を容易に把握することができ、病巣部がどのように変化しているかを容易に比較することができる。
【０１５４】
次に、表示制御部５１６により実行される表示処理の流れについて説明する。図１８及び図１９は、表示制御部により実行される第２及び第３の表示処理のフローチャートである。なお、図１８及び図１９に示すフローチャートの実行に際して、表示部２０６には、比較元ＣＴ画像４０１と比較先ＣＴ画像４１１が並列に表示されているものとする。また、比較先ＣＴ画像４１１に対しては、大域的位置合わせプログラム１４１及び局所的位置合わせプログラム１４２が実行され、対応ベクトルＤＢ５２２には対応点情報１１００が格納されているものとする。
【０１５５】
ステップＳ１８０１において、表示制御部５１６は、医師等の医療従事者により注目点の指定があるか否かを判定する。ステップＳ１８０１において、注目点の指定がないと判定された場合には、注目点の指定があると判定されるまで待機する。
【０１５６】
一方、ステップＳ１８０１において、注目点の指定があると判定された場合には、ステップＳ１８０２に進む。ステップＳ１８０２において、表示制御部５１６は、注目点の指定が、比較元ＣＴ画像４０１において行われたか否かを判定する。
【０１５７】
ステップＳ１８０２において、注目点の指定が、比較元ＣＴ画像４０１において行われたと判定された場合には、ステップＳ１８０３に進む。ステップＳ１８０３において、表示制御部５１６は、指定された注目点の座標を識別する。
【０１５８】
ステップＳ１８０４において、表示制御部５１６は、指定された比較元ＣＴ画像４０１において識別した座標の表示態様を変更する。
【０１５９】
ステップＳ１８０５において、表示制御部５１６は、識別した座標に対する比較先ＣＴ画像４１１上の対応点の座標を、対応ベクトルＤＢ５２２を参照することにより抽出する。また、表示制御部５１６は、抽出した対応点の座標を、大域的位置合わせ前の座標に変換する。
【０１６０】
ステップＳ１８０６において、表示制御部５１６は、比較先ＣＴ画像４１１上において、ステップＳ１８０５において変換された座標の表示態様を変更する。
【０１６１】
ステップＳ１８０７において、表示制御部５１６は、医師等の医療従事者により注目点の拡大表示の指示が入力されたか否かを判定する。ステップＳ１８０７において、注目点の拡大表示の指示が入力されなかったと判定した場合には、ステップＳ１８１０に進む。
【０１６２】
一方、ステップＳ１８０７において、注目点の拡大表示の指示が入力されたと判定された場合には、ステップＳ１８０８に進む。ステップＳ１８０８において、表示制御部５１６は、指定された比較元ＣＴ画像４０１において、識別した座標を含む所定領域を拡大表示する。
【０１６３】
ステップＳ１８０９において、表示制御部５１６は、比較先ＣＴ画像４１１において、変換された座標を含む所定領域を拡大表示する。
【０１６４】
ステップＳ１８１０において、表示制御部５１６は、表示処理の終了指示が入力されたか否かを判定し、入力されていないと判定した場合には、ステップＳ１８０１に戻る。
【０１６５】
一方、ステップＳ１８０２において、注目点の指定が、比較先ＣＴ画像４１１において行われたと判定された場合には、図１９のステップＳ１９０１に進む。
【０１６６】
ステップＳ１９０１において、表示制御部５１６は、指定された注目点の座標を識別する。ステップＳ１９０２において、表示制御部５１６は、比較先ＣＴ画像４１１において、識別した座標の表示態様を変更する。
【０１６７】
ステップＳ１９０３において、表示制御部５１６は、識別した座標に対して、大域的位置合わせを行った後の座標を算出する。また、表示制御部５１６は、大域的位置合わせを行った後の座標に対応する、比較元ＣＴ画像４０１上の対応点の座標を、対応ベクトルＤＢ５２２を参照することにより抽出する。
【０１６８】
ステップＳ１９０４において、表示制御部５１６は、ステップＳ１９０３において抽出された、比較元ＣＴ画像４０１上の座標の表示態様を変更する。
【０１６９】
ステップＳ１９０５において、表示制御部５１６は、医師等の医療従事者により注目点の拡大表示の指示が入力されたか否かを判定する。ステップＳ１９０５において、注目点の拡大表示の指示が入力されなかったと判定した場合には、ステップＳ１８１０に進む。
【０１７０】
一方、ステップＳ１９０５において、注目点の拡大表示の指示が入力されたと判定された場合には、ステップＳ１９０６に進む。ステップＳ１９０６において、表示制御部５１６は、指定された比較先ＣＴ画像４１１において、識別した座標を含む所定領域を拡大表示する。
【０１７１】
ステップＳ１９０７において、表示制御部５１６は、比較元ＣＴ画像４０１において、抽出された座標を含む所定領域を拡大表示する。その後、ステップＳ１８１０に進む。
【０１７２】
ステップＳ１８１０において、表示制御部５１６は、表示処理の終了指示が入力されたか否かを判定し、入力されたと判定した場合には、表示処理を終了する。
【０１７３】
以上のとおり、比較元ＣＴ画像４０１と比較先ＣＴ画像４１１とが並列に表示された場合であっても、本実施形態によれば、一方のＣＴ画像において注目点が指定されると他方のＣＴ画像においても対応する位置にて注目点を自動的に指定することができる。また、注目点を含む領域の拡大表示が指示されると他方のＣＴ画像においても対応する点を含む領域を自動的に拡大表示することができる。
【０１７４】
この結果、一方のＣＴ画像で指定した病巣部が時間経過により、どのように変化したかを容易に比較することができる。
【０１７５】
［第３の実施形態］
上記各実施形態では、患者の体軸方向の複数のＣＴ画像の突き合わせ方法について特に言及しなかったが、患者の体軸方向の複数のＣＴ画像の突き合わせ方法は、任意の方法が適用可能である。
【０１７６】
例えば、患者の頭から体軸方向に複数のＣＴ画像が撮影されていた場合、頭からの枚数が同じＣＴ画像どうしを突き合わせるようにしてもよい。
【０１７７】
また、上記各実施形態では、画像表示装置１２０が大域的位置合わせプログラム１４１、局所的位置合わせプログラム１４２及び同期表示プログラム１４３を有するものとしたが、画像表示装置１２０は、同期表示プログラム１４３のみを有していてもよい。
【０１７８】
この場合、その他のプログラムについては、例えば、別途設けたサーバ装置等にインストールされてもよい。
【０１７９】
また、上記各実施形態では、ＣＴ画像を表示する場合について説明したが、ＣＴ画像以外の医療画像（例えば、ＭＲＩ（Magnetic resonance imaging）画像）を表示する場合に適用してもよい。
【０１８０】
なお、開示の技術では、以下に記載する付記のような形態が考えられる。
（付記１）
異なる時期に撮影された撮影箇所が対応するＣＴ画像のうち、第１のＣＴ画像上で指定された箇所に対応する箇所を、第２の時期に撮影された第２のＣＴ画像上で、該第２のＣＴ画像の他の箇所とは異なる態様で表示する表示部、
を有することを特徴とする画像表示装置。
（付記２）
前記第２のＣＴ画像を前記第１のＣＴ画像に基づいて位置合わせする処理部を更に有し、
前記処理部は、
前記第１のＣＴ画像に含まれる複数の点と、前記第２のＣＴ画像に含まれる複数の点との対応関係を、それぞれの点を含む所定範囲内の各画素の輝度値に応じて算出される特徴量に基づいて決定するマッチング部を更に有することを特徴とする付記１に記載の画像表示装置。
（付記３）
前記処理部は、
前記マッチング部において、前記第１のＣＴ画像に含まれる点に対応する、前記第２のＣＴ画像に含まれる対応点が決定されなかった場合に、該第１のＣＴ画像に含まれる点に対応する、前記第２のＣＴ画像に含まれる対応点を推定する推定部を更に有することを特徴とする付記２に記載の画像表示装置。
（付記４）
前記推定部は、
前記第１のＣＴ画像及び前記第２のＣＴ画像のそれぞれから、臓器の辺縁部分を特定し、
前記第１のＣＴ画像上の臓器の辺縁部分から、該第１のＣＴ画像で指定された箇所を通り対向する２点の辺縁部分を特定し、それらの辺縁部分の特徴量をもとに、前記第２のＣＴ画像の臓器の辺縁部分から前記第１のＣＴ画像で特定した辺縁部分に対応する辺縁部分を特定し、
前記第１のＣＴ画像で指定された箇所から、前記第１のＣＴ画像において特定した各辺縁部分までの距離を特定し、
前記第２のＣＴ画像から特定した前記第１のＣＴ画像の各々の辺縁部分に対応する辺縁部分から、特定した距離の比を満たす箇所を特定する、
ことを特徴とする付記３に記載の画像表示装置。
（付記５）
前記表示部は、
前記第１のＣＴ画像上で指定された箇所に対応する箇所を、前記処理部により位置合わせされた前記第２のＣＴ画像上で、該位置合わせされた第２のＣＴ画像の他の箇所とは異なる態様で表示することを特徴とする付記２に記載の画像表示装置。
（付記６）
前記表示部は、
前記第１のＣＴ画像上で指定された箇所に対応する箇所を、前記マッチング部により決定された対応点または前記推定部により推定された対応点に基づいて算出し、該算出した箇所を、前記第２のＣＴ画像上で、該第２のＣＴ画像の他の箇所とは異なる態様で表示することを特徴とする付記４に記載の画像表示装置。
（付記７）
前記表示部は、
前記第１のＣＴ画像上で指定された箇所について拡大表示の指示があった場合、該第１のＣＴ画像上で指定された箇所に対応する箇所を、前記第２のＣＴ画像上で拡大表示することを特徴とする付記５または６に記載の画像表示装置。
（付記８）
異なる時期に撮影された撮影箇所が対応するＣＴ画像から臓器の辺縁部分を特定し、
一方のＣＴ画像上の臓器の辺縁部分から、一方のＣＴ画像で指定された位置を通り対向する２点の辺縁部分を特定し、それらの辺縁部分の特徴量をもとに、他方のＣＴ画像の臓器の辺縁部分から一方のＣＴ画像で特定した辺縁部分に対応する辺縁部分を特定し、
一方のＣＴ画像で指定された位置から、一方のＣＴ画像において特定した各辺縁部分までの距離を特定し、
他方のＣＴ画像から特定した一方のＣＴ画像の各々の辺縁部分に対応する辺縁部分から、特定した距離の比を満たす位置を特定する、
ことを特徴とする画像表示装置。
（付記９）
異なる時期に撮影された撮影箇所が対応するＣＴ画像のうち、第１のＣＴ画像上で指定された箇所に対応する箇所を、第２の時期に撮影された第２のＣＴ画像上で、該第２のＣＴ画像の他の箇所とは異なる態様で表示する、
処理を、コンピュータに実行させることを特徴とするプログラム。
（付記１０）
コンピュータによる画像表示方法であって、該コンピュータが、
異なる時期に撮影された撮影箇所が対応するＣＴ画像のうち、第１のＣＴ画像上で指定された箇所に対応する箇所を、第２の時期に撮影された第２のＣＴ画像上で、該第２のＣＴ画像の他の箇所とは異なる態様で表示する、
ことを特徴とする画像表示方法。
【０１８１】
なお、上記実施形態に挙げた構成等に、その他の要素との組み合わせなど、ここで示した構成に本発明が限定されるものではない。これらの点に関しては、本発明の趣旨を逸脱しない範囲で変更することが可能であり、その応用形態に応じて適切に定めることができる。
【符号の説明】
【０１８２】
１００    ：ＣＴ画像撮影システム
１１０    ：ＣＴ装置
１２０    ：画像表示装置
１３０    ：画像ＤＢ
１４１    ：大域的位置合わせプログラム
１４２    ：局所的位置合わせプログラム
１４３    ：同期表示プログラム
１５０    ：位置合わせＤＢ
５１１    ：表示画像入力部
５１２    ：特徴点・特徴量抽出部
５１３    ：特徴点マッチング部
５１４    ：対応点推定部
５１５    ：位置合わせ処理部
５１６    ：表示制御部
５２１    ：特徴量ＤＢ
５２２    ：対応ベクトルＤＢ
８００    ：特徴量情報
８１０    ：特徴量情報
１１００  ：対応点情報
【図１】
2015208539.tif 000003
【図２】
2015208539.tif 000004
【図３】
2015208539.tif 000005
【図４】
2015208539.tif 000006
【図５】
2015208539.tif 000007
【図６】
2015208539.tif 000008
【図７】
2015208539.tif 000009
【図８】
2015208539.tif 000010
【図９】
2015208539.tif 000011
【図１０】
2015208539.tif 000012
【図１１】
2015208539.tif 000013
【図１２】
2015208539.tif 000014
【図１３】
2015208539.tif 000015
【図１４】
2015208539.tif 000016
【図１５】
2015208539.tif 000017
【図１６】
2015208539.tif 000018
【図１７】
2015208539.tif 000019
【図１８】
2015208539.tif 000020
【図１９】
2015208539.tif 000021