(57)【特許請求の範囲】
【請求項１】
第１の時間に立体を複数枚にスライスして得た第１の撮像群から、対象の指定を受け付け、
前記第１の撮像群の前記対象が指定されたスライスの画像から、前記指定を受け付けた対象を挟み込む複数の特徴点を抽出し、該抽出した複数の特徴点の各々が、第２の時間に前記立体を複数枚にスライスして得た第２の撮像群のいずれのスライスの位置にあるかを特定し、
前記特定した複数の特徴点を結ぶ線の各スライスの位置において、前記スライスの位置範囲内の所定の数の他のスライスの画像を用いて、前記特定した複数の特徴点の間の画像を生成して出力する、
ことを特徴とする画像処理装置。
【請求項２】
第１の時間に立体を複数枚にスライスして得た第１の撮像群から、対象の指定を受け付け、
前記第１の撮像群の前記対象が指定されたスライスの画像から、前記指定を受け付けた対象を挟み込む複数の特徴点を抽出し、該抽出した複数の特徴点の各々が、第２の時間に前記立体を複数枚にスライスして得た第２の撮像群のいずれのスライスの位置にあるかを特定し、
前記特定した複数の特徴点を結ぶ線の各スライスの位置において、前記スライスの位置範囲内の所定の数の他のスライスの画像を用いて、前記特定した複数の特徴点の間の画像を生成して出力する、
処理を、コンピュータに実行させることを特徴とする画像処理プログラム。
【請求項３】
コンピュータによる画像処理方法であって、該コンピュータが、
第１の時間に立体を複数枚にスライスして得た第１の撮像群から、対象の指定を受け付け、
前記第１の撮像群の前記対象が指定されたスライスの画像から、前記指定を受け付けた対象を挟み込む複数の特徴点を抽出し、該抽出した複数の特徴点の各々が、第２の時間に前記立体を複数枚にスライスして得た第２の撮像群のいずれのスライスの位置にあるかを特定し、
前記特定した複数の特徴点を結ぶ線の各スライスの位置において、前記スライスの位置範囲内の所定の数の他のスライスの画像を用いて、前記特定した複数の特徴点の間の画像を生成して出力する、
ことを特徴とする画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理装置、画像処理方法及び画像処理プログラムに関する。
【背景技術】
【０００２】
医療の現場では、医師が時間経過による患部の変化を確認することがある。確認に際して、例えばＣＴ（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）画像）を表示する場合に適用してもよい。
【０１４３】
なお、開示の技術では、以下に記載する付記のような形態が考えられる。
（付記１）
第１の時間に立体を複数枚にスライスして得た第１の撮像群から、対象の指定を受け付け、
前記第１の撮像群の前記対象が指定されたスライスの画像から、前記指定を受け付けた対象を挟み込む複数の特徴点を抽出し、該抽出した複数の特徴点の各々が、第２の時間に前記立体を複数枚にスライスして得た第２の撮像群のいずれのスライスの位置にあるかを特定し、
前記第２の撮像群から、前記特定した複数の特徴点のスライスの位置に基づいて特定したスライスの位置範囲に該当し、かつ、前記抽出した複数の特徴点の平面上の位置に基づいて特定した平面上の範囲に該当する画像を出力する、
ことを特徴とする画像処理装置。
（付記２）
前記第１の撮像群から、指定されたＸＹＺ軸で特定される一線上から前記複数の特徴点を抽出し、前記第２の撮像群の先に指定された２軸で特定される平面上から、前記抽出した複数の特徴点それぞれに閾値を超えて近似する箇所を特定し、残りの１軸の位置を特定することで、前記抽出した複数の特徴点の各々が、前記第２の撮像群のいずれのスライスの位置にあるかを特定することを特徴とする付記１に記載の画像処理装置。
（付記３）
前記特定した箇所を結ぶ線の各位置において、前記残りの１軸の方向にある所定の数のスライスの画像を用いて、前記特定した箇所の間の画像を再構成して出力することを特徴とする付記２に記載の画像処理装置。
（付記４）
再構成に用いられる前記スライスの画像の数は、前記抽出した複数の特徴点の平面上の距離に比例して設定されることを特徴とする付記３に画像処理装置。
（付記５）
前記特定した箇所の間の画像を再構成して得られた再構成画像を検索することで、前記対象に対応する箇所を抽出することを特徴とする付記３に記載の画像処理装置。
（付記６）
前記特定した箇所を結ぶ線の各位置において、前記残りの１軸の方向にある所定の数のスライスの画像を検索することで、前記対象に対応する箇所を抽出することを特徴とする付記３に記載の画像処理装置。
（付記７）
前記指定を受け付けた対象の移動に応じて、前記対象に対応する箇所が抽出されたスライスの画像を出力することを特徴とする付記６に記載の画像処理装置。
（付記８）
第１の時間に立体を複数枚にスライスして得た第１の撮像群から、対象の指定を受け付け、
前記第１の撮像群の前記対象が指定されたスライスの画像から、前記指定を受け付けた対象を挟み込む複数の特徴点を抽出し、該抽出した複数の特徴点の各々が、第２の時間に前記立体を複数枚にスライスして得た第２の撮像群のいずれのスライスの位置にあるかを特定し、
前記第２の撮像群から、前記特定した複数の特徴点のスライスの位置に基づいて特定したスライスの位置範囲に該当し、かつ、前記抽出した複数の特徴点の平面上の位置に基づいて特定した平面上の範囲に該当する画像を出力する、
処理を、コンピュータに実行させることを特徴とする画像処理プログラム。
（付記９）
コンピュータによる画像処理方法であって、該コンピュータが、
第１の時間に立体を複数枚にスライスして得た第１の撮像群から、対象の指定を受け付け、
前記第１の撮像群の前記対象が指定されたスライスの画像から、前記指定を受け付けた対象を挟み込む複数の特徴点を抽出し、該抽出した複数の特徴点の各々が、第２の時間に前記立体を複数枚にスライスして得た第２の撮像群のいずれのスライスの位置にあるかを特定し、
前記第２の撮像群から、前記特定した複数の特徴点のスライスの位置に基づいて特定したスライスの位置範囲に該当し、かつ、前記抽出した複数の特徴点の平面上の位置に基づいて特定した平面上の範囲に該当する画像を出力する、
ことを特徴とする画像処理方法。
【０１４４】
なお、上記実施形態に挙げた構成等に、その他の要素との組み合わせ等、ここで示した構成に本発明が限定されるものではない。これらの点に関しては、本発明の趣旨を逸脱しない範囲で変更することが可能であり、その応用形態に応じて適切に定めることができる。
【符号の説明】
【０１４５】
１００    ：ＣＴ画像撮影システム
１１０    ：ＣＴ装置
１２０    ：画像処理装置
１３０    ：画像ＤＢ
１４０    ：診断支援部
１４１    ：第１のレジストレーション部
１４２    ：第２のレジストレーション部
１４３    ：表示制御部
１５０    ：サンプリング点ＤＢ
３００    ：並列表示画面
４０１    ：拡大表示画面
４０２    ：拡大表示画面
６０１    ：サンプリング点抽出部
６０２    ：対応点検索部
６０３    ：ＤＢ生成部
１００１  ：サンプリング点取得部
１００２  ：検索範囲決定部
１００３  ：対応点検索部
１３０５  ：部分領域
１３０６  ：部分領域
１４００  ：再構成画像
１７００  ：部分領域群
２００２  ：拡大表示画面
【図１】
0006528386.tif 000002
【図２】
0006528386.tif 000003
【図３】
0006528386.tif 000004
【図４】
0006528386.tif 000005
【図５】
0006528386.tif 000006
【図６】
0006528386.tif 000007
【図７】
0006528386.tif 000008
【図８】
0006528386.tif 000009
【図９】
0006528386.tif 000010
【図１０】
0006528386.tif 000011
【図１１】
0006528386.tif 000012
【図１２】
0006528386.tif 000013
【図１３】
0006528386.tif 000014
【図１４】
0006528386.tif 000015
【図１５】
0006528386.tif 000016
【図１６】
0006528386.tif 000017
【図１７】
0006528386.tif 000018
【図１８】
0006528386.tif 000019
【図１９】
0006528386.tif 000020
【図２０】
0006528386.tif 000021
【図２１】
0006528386.tif 000022
【図２２】
0006528386.tif 000023
【図２３】
0006528386.tif 000024
【図２４】
0006528386.tif 000025