(57)【要約】【課題】過去の診断情報で抽出された関心領域に基づいてマッチング・フィッティングのパラメータを変化させて、経過観察における差分画像精度の向上を図る。【解決手段】Ｓ１１で現在画像及び過去画像に複数のマッチングＲＯＩを設定する。Ｓ１２で各対応するマッチングＲＯＩ間でマッチング処理が行われ、現在画像上のマッチングＲＯＩの位置に対する過去画像上の対応点が計算される。Ｓ１３で求められた各マッチングＲＯＩにおけるシフトベクトルに対して二次元のフィッティングを行う。現在画像のマッチングＲＯＩが先に入力された関心領域に含まれるか、その一部がオーバーラップしている場合には、そのシフトベクトルに対しては他のシフトベクトルに対して小さな重み付けが与えられる。得られたフィッティング係数に基づいて、Ｓ１４において現在画像が変形される。変形された現在画像は、Ｓ１５において過去画像との差分処理が行われ、差分画像が表示される。【選択図】図４
2005270635.tif 000002
【特許請求の範囲】
【請求項１】
第１の画像データ及び第２の画像データの少なくとも何れか一方に関する領域情報に基づいて経時差分処理を行う画像処理方法であって、前記領域情報を用いて前記第１の画像データと第２の画像データとの位置合わせを行う位置合わせ工程と、位置合わせ後の前記第１の画像データと第２の画像データとの差分処理を行う差分処理工程とを備えることを特徴とする画像処理方法。
【請求項２】
前記位置合わせ工程は、前記第１の画像データ上の複数の位置を示す代表点Ｓと、該代表点Ｓに対応する前記第２の画像データ上の対応点を算出する工程と、前記代表点Ｓと前記対応点Ｓとの位置の差をシフトベクトルとして算出する工程と、前記シフトベクトルを重み係数で重み付けし変更係数を算出する工程と、該変更係数に従い、前記画像データの少なくとも一方の画像データの座標を変更する工程とを備え、前記重み付け係数は前記領域情報に基づき前記重み係数の値が定められることを特徴とする請求項１に記載の画像処理方法。
【請求項３】
前記第１の画像データのｎ個の代表点の座標を（Ｘｉ，Ｙｉ）、前記変更係数をＣｎ、前記対応点の座標とのＸ座標での差をΔｘｉ、重み係数をｗｉとするとき、
【数１】
2005270635.tif 000006
を最小とするように前記変更係数Ｃｎを算出し、前記重み係数ｗｉは前記領域情報に基づいて定めることを特徴とする請求項２に記載の画像処理方法。
【請求項４】
前記領域情報に基づく関心領域内の前記重み付け係数の値を、前記関心領域外の重み付け係数の値に対して小さくすることを特徴とする請求項１又は２に記載の画像処理方法。
【請求項５】
前記領域情報は病変領域の座標であることを特徴とする請求項１〜４の何れか１つの請求項に記載の画像処理方法。
【請求項６】
前記差分処理結果と前記領域情報とを表示する表示工程を備え、該表示工程は前記診断情報に基づいて表示内容を変更することを特徴とする請求項１に記載の画像処理方法。
【請求項７】
前記表示工程は前記差分処理結果に対し、前記過去の診断情報に基づいて決定される関心領域を表すマーカを追加表示することを特徴とする請求項６に記載の画像処理方法。
【請求項８】
第１の画像データ及び第２の画像データの少なくとも何れか一方に関する領域情報に基づいて経時差分処理を行う画像処理方法をコンピュータで実行するためのプログラムであって、前記領域情報を用いて前記第１の画像データと第２の画像データとの位置合わせを行う位置合わせ工程と、位置合わせ後の前記第１の画像データと第２の画像データとの差分処理を行う差分処理工程とをコンピュータに実行させるプログラム。
【請求項９】
第１の画像データ及び第２の画像データの少なくとも何れか一方に関する領域情報に基づいて経時差分処理を行う画像処理方法をコンピュータで実行するためのプログラムを記録したコンピュータ読み取り可能なコンピュータ可読記憶媒体であって、前記領域情報を用いて前記第１の画像データと第２の画像データとの位置合わせを行う位置合わせ工程と、位置合わせ後の前記第１の画像データと第２の画像データとの差分処理を行う差分処理工程とをコンピュータに実行させるためのプログラムを記録しコンピュータ読み取りを可能とするコンピュータ可読記憶媒体。
【請求項１０】
第１の画像データ及び第２の画像データの少なくとも何れか一方に関する領域情報に基づいて経時差分処理を行う画像処理装置であって、前記領域情報を用いて前記第１の画像データと第２の画像データとの位置合わせを行う位置合わせ手段と、位置合わせ後の前記第１の画像データと第２の画像データとの差分処理を行う差分処理手段とを備えたことを特徴とする画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、単純Ｘ線画像等の画像を基に被写体の経時変化を抽出する画像処理方法及び処理装置に関するものである。
【背景技術】
【０００２】
近年の医療画像分野においては、ＣＡＤ（コンピュータ支援診断：Computer Aided Diagnosis）の研究が活発に行われており、単純Ｘ線、ＣＴ（Computer Tomography）などの画像をコンピュータにより解析し、病変部分の疑いがある部分を検出することで、疾病の早期発見に貢献することが期待されている。
【０００３】
ＣＡＤには、対象とする画像の種類や検出対象により様々な種類が考えられているが、その中に同一部位を撮影した画像の差分を求めることにより、経時変化を強調した画像を得る経時差分方式によるものが注目されている。
【０００４】
この経時差分ＣＡＤにおいては、例えば時間的に異なる時点で撮影された１組の胸部単純Ｘ線画像を入力し、画像解析を行ってそれぞれの画像において解剖学的に同一となる位置を求め、現在又は過去の何れか一方の画像を変形し、画素毎の差分処理を行う。経時差分ＣＡＤは時間間隔をおいて撮影された画像間の変化を、正常な解剖学的構造を除去した差分画像として表示できるため、病変部の変化の様子を観察することに適している。
【０００５】
このような経時差分ＣＡＤの手法に関しては、例えば特許文献１等に開示されており、図７はこの特許文献１に開示されている経時差分の処理工程を示している。
【０００６】
【特許文献１】米国特許５，９８２，９１５号公報
【０００７】
図７では、２０４８×２０４８の解像度を持つ現在及び過去画像を工程Ｐ１において所定の解像度変換手法により、例えば５１２×５１２及び１２８×１２８の大きさに画素数を減じ、更に工程Ｐ２において１２８×１２８の画像に対して、ガウシアンフィルタにより高周波成分をカットすることで、画像の概要を表す縮小画像を生成する。
【０００８】
１２８×１２８の２つの縮小画像に対して、画像間の全体的な位置ずれを表すグローバルシフトベクトルを求める。更に、工程Ｐ４においてグローバルシフトベクトルが表す２つの画像間の位置関係に基づいて、５１２×５１２の解像度を持つ各画像に関心領域（ＲＯＩ：Region Of Interest）を設定して、各画像において対応する関心領域毎にマッチングを行って、更に詳細なシフトベクトルを求め、これを二次元多項式で補間して、任意の位置のシフトベクトルが計算可能なパラメータを求め、このパラメータを用いて何れか一方の画像を変形して差分を求め差分画像を出力する。
【０００９】
このようにして、被写体の撮像時における姿勢の違い等に起因する画像上の位置ずれを補正し、次の差分処理を行う際に縮小画像を用いて全体的なずれを求めるという２段階補正を行うことで、より精度の高い位置補正を行っている。
【００１０】
しかし、原理的に２つの画像の差分処理を行っているために、元々病変部が存在しており、それが変化していない場合には、差分画像上には全くこの病変部による差信号は得られない。従って、何らかの方法により既に病変の存在が確認されている場合において、その情報が経時差分処理において用いられないことは、誤判断に結びつく可能性がある。
【００１１】
また、既に病変部の存在が確認されており、その後の経過観察に差分画像が用いられている場合に、前述した差分処理において病変部分の変化がより詳細に観察可能な精度の高い差分処理画像が得られることが望ましい。
【発明の開示】
【発明が解決しようとする課題】
【００１２】
従来においては、患者に関する過去の診断情報を有していても、経時差分処理に用いられなかったため、診断精度を落とす場合があり、それに対する種々の対策が望まれている。
【００１３】
本発明の目的は、前述した問題点に鑑み、効果的な経時差分処理を行う画像処理方法及び処理装置を提供することにある。
【課題を解決するための手段】
【００１４】
上記目的を達成するための本発明に係る画像処理方法は、第１の画像データ及び第２の画像データの少なくとも何れか一方に関する領域情報に基づいて経時差分処理を行う画像処理方法であって、前記領域情報を用いて前記第１の画像データと第２の画像データとの位置合わせを行う位置合わせ工程と、位置合わせ後の前記第１の画像データと第２の画像データとの差分処理を行う差分処理工程とを備えることを特徴とする。
【００１５】
また、本発明に係るプログラムは、第１の画像データ及び第２の画像データの少なくとも何れか一方に関する領域情報に基づいて経時差分処理を行う画像処理方法をコンピュータで実行するためのプログラムであって、前記領域情報を用いて前記第１の画像データと第２の画像データとの位置合わせを行う位置合わせ工程と、位置合わせ後の前記第１の画像データと第２の画像データとの差分処理を行う差分処理工程とをコンピュータに実行させる。
【００１６】
更に、本発明に係るコンピュータ可読記憶媒体は、第１の画像データ及び第２の画像データの少なくとも何れか一方に関する領域情報に基づいて経時差分処理を行う画像処理方法をコンピュータで実行するためのプログラムを記録したコンピュータ読み取り可能なコンピュータ可読記憶媒体であって、前記領域情報を用いて前記第１の画像データと第２の画像データとの位置合わせを行う位置合わせ工程と、位置合わせ後の前記第１の画像データと第２の画像データとの差分処理を行う差分処理工程とをコンピュータに実行させるためのプログラムを記録しコンピュータ読み取りを可能とする。
【００１７】
本発明に係る画像処理装置は、第１の画像データ及び第２の画像データの少なくとも何れか一方に関する領域情報に基づいて経時差分処理を行う画像処理装置であって、前記領域情報を用いて前記第１の画像データと第２の画像データとの位置合わせを行う位置合わせ手段と、位置合わせ後の前記第１の画像データと第２の画像データとの差分処理を行う差分処理手段とを備えたことを特徴とする。
【発明の効果】
【００１８】
本発明に係る画像処理方法及び処理装置によれば、画像効果的な経時差分処理を行うことができる。
【発明を実施するための最良の形態】
【００１９】
本発明を図１〜図６に図示の実施例に基づいて詳細に説明する。
【実施例１】
【００２０】
図１は実施例１の画像処理システムのブロック回路構成図である。医療画像を取得するための撮像部１の出力は、画像を記憶する例えばハードディスク、光磁気ディスク等から成る蓄積型の記憶装置である保存部２、及び画像入力部３に接続され、画像入力部３の出力は差分処理部４、表示制御部５、及びＣＲＴモニタ、液晶ディスプレイ、或いはプリンタ等から成る表示部６に順次に接続されている。また、保存部２の出力は画像入力部３、診断情報入力部７に接続され、診断情報入力部７の出力は差分処理部４、表示制御部５に接続されている。
【００２１】
撮像部１は診断に用いられる画像を新たに生成する撮像装置であり、例えばＸ線撮像装置等のモダリティである。本実施例においては、撮像部１は胸部を対象とした単純Ｘ線撮像装置とするが、他の撮像装置であっても本発明の適用が可能である。
【００２２】
撮像部１で生成された現在画像Ｉｃは、保存部２及び画像入力部３に出力される。ただし、必ずしも撮像装置に直結されている必要はなく、ネットワークを介して画像の伝送が可能な画像サーバであってもよい。また、撮像部１から出力されるデータは現在画像Ｉｃだけではなく、撮影の日付、患者氏名、患者のＩＤ番号、生年月日、所見等の付帯情報Ｄｃも同時に出力される。
【００２３】
画像入力部３は撮像部１から撮影された画像を入力する一方で、保存部２から同一患者の過去に撮影された過去画像Ｉｐ及び付帯情報Ｄｐを入力する。本実施例においては、同一患者の過去に撮影された過去画像Ｉｐが、保存部２に保存されているものとして説明する。ただし、現在画像Ｉｃ及び付帯情報Ｄｃが撮像部１から入力されることに限定されるものではなく、現在画像Ｉｃ及び付帯情報Ｄｃも保存部２から入力されてもよい。
【００２４】
また、保存部２から出力された過去画像Ｉｐの付帯情報Ｄｐは診断情報入力部７にも出力される。ここで、付帯情報Ｄｐに含まれる所見が有所見の場合に、本実施例においては付帯情報Ｄｐに関心領域が含まれているものとする。ここで、関心領域ＲＯＩは過去画像Ｉｐにおける病変部分の位置を現す座標値であり、例えば図２に示すように、病変部分を含む矩形状のｘ１、ｘ２、ｙ１、ｙ２として定義されているものとする。診断情報入力部７はこの関心領域ＲＯＩの座標値を差分処理部４に出力する。
【００２５】
一方、画像入力部３は現在画像Ｉｃ及び過去画像Ｉｐを差分処理部４に出力する。差分処理部４は従来技術の項において説明したように、基本的に図７に示す処理により、２つの画像の差分画像を生成して出力するが、本実施例では図７における工程Ｐ４に対応する差分処理において、関心領域の情報が用いられる。
【００２６】
図３は現在画像Ｉｃに対して胸郭を関心領域（ＲＯＩ：Region Of Interest）として設定した画像の説明図であり、図中の各点は関心領域に対して均等な間隔で置かれた代表点を示している。
【００２７】
図４は差分処理の基本フローチャート図である。ステップＳ１１において、現在画像Ｉｃの代表点ごとに複数のマッチングＲＯＩｍｃを設定する。ここで、ＲＯＩｍｃは各代表点を中心とした一定範囲の領域であり、例えば画素２１×２１の方形領域とする。なお、これらのマッチングＲＯＩは、図７における工程Ｐ３において求められたグローバルシフトベクトルを考慮し、胸郭である各ＲＯＩが２つの画像において略一致している状態で設定される。
【００２８】
次に、ステップＳ１２において各対応するマッチングＲＯＩｍｃ毎にマッチング処理が行われ、現在画像Ｉｃ上の各代表点の位置に対する過去画像Ｉｐ上の対応点が計算される。このときのマッチング処理としては、相互相関係数が最大値となる位置を求めるなどの公知の方法が用いられる。
【００２９】
具体的には、代表点を中心とする方形領域と最も相互相関係数が最大となる領域を算出し、算出した領域の中心点を対応点とするものである。これにより、現在画像Ｉｃ上の代表点（ｘｉ，ｙｉ）を過去画像Ｉｐ上の対応点（ｘｉ'，ｙｉ'）に変換するシフトベクトル（Δｘｉ，Δｙｉ）が求められる。
【００３０】
続いて、ステップＳ１３においては、ステップＳ１２において求められた代表点、対応点間におけるシフトベクトル（Δｘｉ，Δｙｉ）に対して二次元のフィッティングを行うが、このフィッティングには例えば二次元の多項式近似が用いられる。具体的には、式（１）のＳｘを最小とする係数Ｃｎを算出するものである。最小二乗法により多項式の係数が計算されるが、ここで、計算の際に各シフトベクトルに対して重み係数ｗｉを用いて重み付けを行い、フィッティングにおける寄与の度合いを調整する。
【００３１】
【数１】
2005270635.tif 000003
【００３２】
即ち、現在画像ＩｃのマッチングＲＯＩｍｃが先に診断情報入力部から差分処理部４に入力された関心領域に含まれるか、その一部がオーバーラップしている場合には、そのシフトベクトルに対しては他のシフトベクトルに対して小さな重み付けが与えられる。例えば、関心領域にかかるシフトベクトルに対しては重み係数として０．５を与え、その他のシフトベクトルに対しては、１等の値を与えればよい。
【００３３】
このようにして得られたフィッティング係数に基づいて、ステップＳ１４において現在画像Ｉｃが変形される。即ち、現在画像ＩｃのＸ座標（ｘｉ）は、過去画像ＩｐのＸ座標（ｘｉ'）に移動する。
【００３４】
【数２】
2005270635.tif 000004
【００３５】
同様に、各座標点のＹ方向のシフトベクトルも計算される。
【００３６】
このように、二次元多項式補間は最小自乗誤差法などにより、入力された二次元データを利用し、各多項式項の係数を求めて、入力された二次元データを近似する。二次元多項式による補間は、対応点全体領域を渡って係数を求めるので、入力された二次元データを滑らかに補間でき、データに潜んでいるノイズも或る程度抑制できる。そのため、対応点シフトベクトルを補間する際に、二次元多項式補間を用いることが多い。
【００３７】
次に、変形された現在画像Ｉｃは、ステップＳ１５において過去画像Ｉｐとの差分処理が行われ、差分画像が出力され表示制御部５に対して出力される。表示制御部５は差分画像、現在画像Ｉｃ、過去画像Ｉｐ及びそれらの付帯情報を所定の配置に構成して表示部６に出力する。
【００３８】
このように、本実施例では差分処理部４においてシフトベクトルのフィッティングを行う際に、別途に保存部２から入力した過去の診断情報における関心領域に基づいて、この関心領域のシフトベクトルをその他のシフトベクトルに対して、より小さな重み付けを行って画像を変形し差分画像を生成する。これによって、正常な構造の信号成分はより減弱されることなり、過去の診断により決定される関心領域に対しては、変化がより強調された差分画像が生成される。
【００３９】
更に、これに加えて画素値の平均値の大きさに応じて重み付けを変更することもできる。胸部Ｘ線画像においては、縦隔及び腹腔部分の輝度レベルは平均的に高い値を持ち、コントラストが低いために、十分なマッチング精度が得られないことがある。従って、関心領域には含まれず、かつ高輝度レベルの領域に属するシフトベクトルに対しては、低い重み付けを与えるようにすればよい。
【００４０】
一方、縦隔、腹腔の領域であっても、それが関心領域に含まれるシフトベクトルに関しては、低い重み付けを与えないようにすれば、診断上重要な部分に対する必要以上の画質低下を防止することができる。
【実施例２】
【００４１】
上述の実施例１においては、過去の診断情報から得られた関心領域に基づいてシフトベクトルを求めるマッチング時の重み付けを制御したが、マッチングを行う際の処理方法を変更してもよい。
【００４２】
実施例２においては、現在画像Ｉｃ及び過去画像Ｉｐに複数のマッチングを設定し、マッチングを行う際の画素間隔を変更する。即ち、マッチング尺度として相互相関係数を用いた場合に、位置（ａ，ｂ）における相互相関係数Ｃ（ａ，ｂ）は次式で計算される。ただし、Ｉｐ（ａ，ｂ）は過去画像Ｉｐの位置（ａ，ｂ）に設定されたマッチングＲＯＩｍｐの画像データ、Ｉｃは現在画像Ｉｃに設定された対応領域の画像データである。
【００４３】
【数３】
2005270635.tif 000005
【００４４】
ここで、現在画像Ｉｃは画像処理のマッチング技術においては屡々テンプレートと呼ばれる。実施例１においては、式（３）で相互相関係数Ｃ（ａ，ｂ）を計算する際に、このテンプレートを１画素ずつ移動して相互相関係数Ｃ（ａ，ｂ）を計算する。
【００４５】
本実施例２においては、診断情報入力部７から得られた関心領域に基づいて、マッチングの際の座標が関心領域の内部に位置しない場合は、１画素ずつテンプレートを移動する。しかし、座標が関心領域の内部に位置する場合は、１画素よりも大きい値、例えば２画素単位でテンプレートを移動する。
【００４６】
これにより、関心領域では若干マッチングの精度が低下する可能性があるが、全体的には処理の速度を向上させ、関心領域以外に関してはより精度の高いマッチングを行うことができる。
【実施例３】
【００４７】
実施例１、２においては、過去の診断情報の関心領域に基づいて適応的な処理を行ったが、他の処理において用いることも可能である。
【００４８】
図５は実施例３における差分処理のフローチャート図である。このフローチャート図において、ステップＳ１５の差分計算までは図４と同様であるが、ステップＳ１６以下において生成された差分画像の画質を次のように評価する。
【００４９】
即ち、ステップＳ１６において、差分画像中の肺野領域に対応する部分に対してサンプル領域を設定する。図６（ａ）は差分画像に対して２つのサンプル領域ＲＯＩｓｒ及びＲＯＩｓｌの設定例である。これらのサンプル領域は過去画像Ｉｐから抽出された肺野領域の中央に設定されるが、肺野抽出の方法としては公知の方法によることができるため、詳細な説明は省略する。
【００５０】
なお図６においては、２個所のサンプル領域が設定されているが、必ずしも２個所である必要はなく、肺野内の診断上重要な位置に任意の数が設定される。ただし、本実施例３において、サンプル領域は図６（ｂ）に示すように、診断情報入力部７から得られた関心領域に重ならない位置に設定されるものとする。
【００５１】
次に、ステップＳ１７において各サンプル領域内の差分画像の画素分布が解析され、分布が所定の閾値以上の場合にステップＳ１８に進み、そうでない場合には差分画像を出力して処理を終了する。画素分布の解析方法としては、画素値のヒストグラム幅や分散を用いることができる。
【００５２】
ステップＳ１８においては、先にステップＳ１３においてフィッティング係数を計算する際に用いたシフトベクトルを新たに選択し、再びステップＳ１３に戻りフィッティング係数が計算される。
【００５３】
このとき、先に画素分布が所定の大きさを超えたサンプル領域に属するシフトベクトルは、再びステップＳ１３においてフィッティング係数が計算される際に用いられないように、ステップＳ１８においてシフトベクトルが選択されている。
【００５４】
これにより、局所的な外乱要因で画像が歪み、当該部分のシフトベクトルに誤差が多く含まれる場合にこれを除外し、全体的な位置合わせに及ぼす影響を抑制することができるが、関心領域に関してはサンプル領域に含めないことにより、常にフィッティングの対象とすることができる。
【００５５】
従って、経過観察等において関心領域の部分に局所的な変化が生じている場合であっても、この部分についてはできるだけ正確な差分画像を生成することが可能となる。
【実施例４】
【００５６】
上述の実施例１〜３においては、差分画像処理の工程において過去の診断情報に基づいた処理が行われたが、生成された差分画像を表示する際に表示形態を切換えてもよい。
【００５７】
実施例４においては、図１において診断情報入力部７からの診断情報は表示制御部５に出力される。表示制御部５は差分処理部４から入力した差分画像を表示部６に出力する際に、診断情報入力部７から入力した関心領域の位置にマーカを重ね合わせて出力する。
【００５８】
差分処理においては、処理対象となる画像の対応する部分間で変化がないと、それは差分画像上には現れない。従って、何らかの病変が過去から引き続き存在したとしても、十分に注意を払わなければ誤判断に繋がる可能性がある。しかし、別に入力した過去の診断情報から関心領域を抽出して重ね合わせることにより、このような誤判断を防止することが可能となる。
【実施例５】
【００５９】
実施例４においては、差分画像の関心領域にマーカを追加表示したが、過去の診断情報において無所見であった場合には、それを示すようにしてもよい。即ち、表示制御部５は診断情報入力部７から関心領域に関する情報が入力されなかった場合は、差分画像に対してマーカを重ね合わせる代りに、無所見であることを示す表示を表示部６に出力する。表示の方法としては、ウィンドウシステムをベースとした場合には、小ウィンドウの追加表示や、ウィンドウ内に過去の所見の有無を示すランプ状の表示を行ってもよい。
【実施例６】
【００６０】
上述の実施例においては、対象患者の付帯情報Ｄｐは過去画像Ｉｐに対応するものとして説明した。しかし、これは必ずしも対応しなくともよく、例えば対象患者の過去のＣＴ撮影結果に対する付帯情報Ｄｐを使用して、過去画像Ｉｐにおいて関心領域ＲＯＩを決定することもできる。
【００６１】
ＣＴ画像の付帯情報Ｄｐから単純Ｘ線画像上での関心領域ＲＯＩを決定する方法を説明すると、近年のマルチスライスＣＴによれば、等方性画像（Isotropic）を得ることができ、この等方性画像からはＲＡＹＳＵＭ処理により単純Ｘ線画像と同等の得られることは周知である。ＲＡＹＳＵＭ処理は、等方性画像が形成する三次元ボリュームデータを空間上の一点から投影する処理である。
【００６２】
その際に、三次元ボリュームデータ中にマークされた関心領域ＲＯＩも同時に投影することが可能である。ＣＴ画像において、関心領域ＲＯＩを示すマークはスライス画像毎に設定されるからである。このようにして作成されたＲＡＹＳＵＭ画像と関心領域ＲＯＩの関係は図２と同様である。
【００６３】
ＣＴ画像と単純Ｘ線画像を例にとって説明したが、ＭＲＩ画像やＰＥＴ画像と単純Ｘ線画像であっても同様の手順で、異なるモダリティによる付帯情報Ｄｐを利用することができる。ＣＴ画像の付帯情報Ｄｐを利用して、単純Ｘ線画像における経時差分をすることの利点は被曝線量の低減にある。ＣＴ撮影は撮影線量が多く、頻回の経過観察には適さないため、被曝線量の小さい単純Ｘ線画像を使用して経過観察することが、対被曝効果の面で有利である。
【００６４】
なお、上述の実施例１〜６において、保存部２に画像と共に保存されている診断情報は、次の何れかの方法によることができる。即ち、医師が画像を観察して入力を行った読影レポートを基に関心領域ＲＯＩを決定することができる。この場合に、関心領域の位置は読影レポート内に手動で設定され、所定のタグと共に保存される。一方、画像内の関心領域がコンピュータによる画像解析手法により、自動的に検出されたものであってもよい。
【００６５】
更に、上述の画像解析手法に関しては、１枚の画像から病変部分を自動的に検出するＣＡＤによるものであってもよいし、１組の画像から得られた差分画像に対して画像解析により病変部分を検出するものであってもよい。
【００６６】
更に、本発明は上記の実施例を実現するためのシステム、装置及び方法のみに限定されるものではなく、上記システム又は装置内のＣ工程ＰＵやＭ工程ＰＵのコンピュータに、実施例を実現するためのソフトウエアのプログラムコードを供給し、このプログラムコードに従って上記システム或いは装置のコンピュータが各種デバイスを動作させることにより、上述の実施例を実現する場合も本発明の範疇に含まれる。
【００６７】
また、この場合にソフトウエアのプログラムコード自体が実施例の機能を実現することになり、そのプログラムコード自体、及びそのプログラムコードをコンピュータに供給するための手段、具体的にはプログラムコードを格納した記憶媒体は本発明の範疇に含まれる。
【００６８】
このようなプログラムコードを格納する記憶媒体としては、例えばフロッピディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。
【００６９】
またコンピュータが供給されたプログラムコードのみに従って各種デバイスを制御することにより、実施例の機能が実現される場合だけではなく、プログラムコードがコンピュータ上で稼働しているＯＳ（オペレーティングシステム）、或いは他のアプリケーションソフト等と共働して実施例を実行する場合にも、このようなプログラムコードは本発明の範疇に含まれる。
【００７０】
更に、この供給されたプログラムコードが、コンピュータの機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに格納された後に、そのプログラムコードの指示に基づいてその機能拡張ボードや機能格納ユニットに備わるＣ工程ＰＵ等が実際の処理の一部又は全部を行い、その処理によって実施例が実現される場合も本発明の範疇に含まれる。
【図面の簡単な説明】
【００７１】
【図１】実施例のブロック回路構成図である。
【図２】関心領域の説明図である。
【図３】胸郭を関心領域として設定した画像の説明図である。
【図４】差分処理の基本フローチャート図である。
【図５】差分処理のフローチャート図である。
【図６】サンプル領域の設定例の説明図である。
【図７】従来例の差分処理の説明図である。
【符号の説明】
【００７２】
１  撮像部
２  保存部
３  画像入力部
４  差分処理部
５  表示制御部
６  表示部
７  診断情報入力部
【図１】
2005270635.tif 000007
【図２】
2005270635.tif 000008
【図３】
2005270635.tif 000009
【図４】
2005270635.tif 000010
【図５】
2005270635.tif 000011
【図６】
2005270635.tif 000012
【図７】
2005270635.tif 000013