JPH0882880A - 放射線画像の位置合せ用対応点決定方法 - Google Patents

Abstract

(57)【要約】 （修正有） 【目的】 放射線画像の位置合せ用の対応点を決定する
方法において、位置合せに適した対応点を求めることが
できるようにする。 【構成】 放射線画像4a上に基準対応点40Ａ〜40Ｄを設
定し、この基準対応点を中心としたテンプレート領域41
Ａ〜41Ｄを設定する。このテンプレート領域41Ａ〜41Ｄ
を用いて放射線画像4b上に置いてテンプレートマッチン
グを行い。放射線画像4b上に複数の対応点を設定する。
この複数の対応点について、（１）テンプレートマッチ
ングにより得られた、テンプレート領域と他の放射線画
像の一部との相関値が所定値以上であること、（２）複
数の対応点近傍の相関値の分布が尾根状とならないこ
と、（３）任意の２基準対応点間の距離と、この２基準
対応点に対応する他の放射線画像上の２対応点の距離と
の比が一定値の範囲内にあること、の評価を行う。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は複数の放射線画像の位置
ずれを補正して画像の位置合せをする際に用いる対応点
を求める方法に関し、とくに詳細には放射線画像の重ね
合わせ処理またはサブトラクション処理のために放射線
画像の位置合せを行う際の対応点を求める方法に関する
ものである。

【０００２】

【従来の技術】蓄積性蛍光体を利用して、人体等の被写
体の放射線画像情報を一旦蓄積性蛍光体のシート（以
下、蓄積性蛍光体シートと称する）に記録し、これを励
起光で走査して輝尽発光させ、この輝尽発光光を光電的
に読み取って画像信号を得、この画像信号を処理して診
断適性の良い被写体の放射線画像を得る方法が知られて
いる。この最終的な画像はハードコピーとして再生した
り、あるいはＣＲＴ上に再生したりすることができる。

【０００３】一方、従来より放射線画像の重ね合せ処理
が公知となっている（例えば特開昭56-11399号参照）。
一般に、放射線画像は診断用その他の目的に使われる
が、その使用に当たっては被写体の微小な放射線吸収差
を良好に検出することが要求される。放射線画像におけ
るこの検出の程度をコントラスト検出能または単に検出
能と呼ぶが、この検出能の高いもの程診断性能も高く、
実用的価値が高い放射線画像であると言うことができ
る。したがって診断性能を高めるため、この検出能を高
くすることが望まれるが、その最も大きな障害要因は各
種ノイズである。重ね合せ処理は、その各種ノイズを大
幅に減少させ、被写体の僅かな放射線吸収差も最終画像
において明確に観察可能にして、検出能を大幅に向上さ
せる方法である。すなわち、複数枚重ねた蓄積性蛍光体
シートに放射線画像を撮影（蓄積記録）し、この複数枚
のシートを読取処理にかけて得た複数の画像信号を加算
処理し、このことにより、前述の各種ノイズを減少させ
るものである。

【０００４】従来、実際にこの重ね合せ処理を行うため
には、例えば、カセッテに蓄積性蛍光体シートを２枚重
ねて入れて被写体の撮影を行い、２枚の蓄積性蛍光体シ
ートに対して通常の読取処理を逐次行って２組の画像信
号を得、この２組の画像信号を加算処理する方法が用い
られている。

【０００５】また一方、従来より放射線画像のサブトラ
クション処理が公知となっている。この放射線画像のサ
ブトラクションとは、異なった条件で撮影した２つの放
射線画像を光電的に読み出してデジタル画像信号を得た
後、これらのデジタル画像信号を両画像の各画素を対応
させて減算処理し、放射線画像中の特定の構造物を抽出
させる差信号を得る方法であり、このようにして得た差
信号を用いれば、特定構造物のみが抽出された放射線画
像を再生することができる。

【０００６】このサブトラクション処理には、基本的に
次の２つの方法がある。すなわち、(1) 造影剤注入によ
り特定の構造物が強調された放射線画像の画像信号か
ら、造影剤が注入されていない放射線画像の画像信号を
引き算（サブトラクト）することによって特定の構造物
を抽出するいわゆる時間サブトラクション処理と、(2)
同一の被写体に対して相異なるエネルギー分布を有する
放射線を照射し、あるいは被写体透過後の放射線をエネ
ルギー分布状態を変えて２つの放射線検出手段に照射し
て、それにより特定の構造物が異なる画像を２つの放射
線画像間に存在せしめ、その後この２つの放射線画像の
画像信号間で適当な重み付けをした上で引き算（サブト
ラクト）を行って、特定の構造物の画像を抽出するいわ
ゆるエネルギーサブトラクション処理である。

【０００７】このサブトラクション処理は特に医療診断
上きわめて有効な方法であるため、近年大いに注目さ
れ、電子工学技術を駆使してその研究、開発が盛んに進
められている。

【０００８】しかしながら、上述したような蓄積性蛍光
体シートを用いた放射線画像の重ね合せ処理方法および
サブトラクション処理方法においては以下のような問題
が生じる。

【０００９】すなわち、蓄積性蛍光体シートを用いた前
記各処理方法においては、２枚（３枚以上の場合もあ
る）の蓄積性蛍光体シートを順次もしくは同時に撮影台
に挿入して重ね合わせまたはサブトラクションすべき放
射線画像を撮影し、その後に蓄積性蛍光体シートを個別
に読取装置に挿入し、その都度蓄積性蛍光体シートに励
起光を照射することにより発せられた輝尽発光光を検出
することにより放射線画像を読み出すが、この過程にお
いては、撮影および読取りに関わるすべての装置の機械
的精度を上昇させたとしても、重ね合わせまたはサブト
ラクションされるべき画像間で位置ズレおよび回転ズレ
が生じることとなる。この結果、重ね合せ処理において
は各種ノイズがこの処理により平均化されて減少するも
のの、画像中の構造物の縁の部分をはじめ画像全体にボ
ケが生じ、観察すべき画像が観察に適さなくなり、また
サブトラクション処理においては消去されるべき画像が
消去されなかったり、逆に抽出すべき画像が消去されて
偽画像が生じて正確なサブトラクション像を得ることが
できなくなる。このように前述した位置ズレおよび回転
ズレにより、診断上重大な支障が生じるということが見
出された。

【００１０】このようなズレが蓄積性蛍光体シートに蓄
積記録された放射線画像情報間に生じると放射線画像は
潜像として蓄積性蛍光体中に蓄積記録されているので、
Ｘ線画像を可視像としてとらえることのできるＸ線写真
フイルムの場合と異なって、目視によって２枚のＸ線写
真を合わせるといったことができず、ズレ補正は極めて
困難なものとなる。

【００１１】さらに、何らかの手段により２つの放射線
画像間に生じる位置ズレおよび回転ズレを検出しえたと
しても読み取られた放射線画像のデータを補正すべく従
来公知の演算処理を行うと、特に回転ズレの補正の際に
多大な時間が費やされ、実用上非常に大きな問題とな
る。

【００１２】そこで本出願人により特開昭58-163338 号
に、基準点または基準線を提供するような形状をもつマ
ーカーを用いた放射線画像のサブトラクション処理方法
を提案している。この方法は、マーカーを放射線画像に
対し固定した位置で２枚の蓄積性蛍光体シートに記録
し、この放射線画像の読取りの際にマーカーを検出し、
位置ズレおよび回転ズレを計算してサブトラクションす
べき放射線画像のいズレか一方をデジタルデータ上で回
転および／または移動し、この放射線画像の対応する各
画素間で画像データの引き算を行うものである。このマ
ーカーを用いた放射線画像のサブトラクション処理方法
における位置合せの工程は、上述した重ね合せ処理方法
にも適応することもできる。その場合、位置合せを行っ
た後に放射線画像の対応する各画素間で画像データの加
算処理を行えばよい。

【００１３】しかしながら、この方法においては放射線
画像の撮影の都度、上述したようなマーカーを被写体と
ともに蓄積性蛍光体シートに蓄積記録しなければならな
い。そして、この蓄積記録した放射線画像のマーカーの
位置と重なる部分からは被写体の画像情報を得ることが
できないという問題がある。

【００１４】そこで、本願出願人により位置合せのため
にマーカー等を用いることなく放射線画像の位置合せを
行う方法が提案されている（特願平4-318533号）。この
方法は、位置合せを行う複数の放射線画像のうちの１つ
の放射線画像にテンプレート領域を設定し、このテンプ
レート領域を用いて他の放射線画像についてテンプレー
トマッチングを行って各放射線画像に少なくとも２つの
マーカーの代りとなる対応点を求め、各放射線画像の対
応点が一致するように各対応点をアフィン変換して、各
放射線画像について回転移動補正、拡大または縮小率補
正および平行移動補正を行う方法である。

【００１５】この方法によれば、画像の中からマーカー
の代りとなる対応点を見つけ出して利用するものである
から、位置合せのためにマーカー等を被写体とともに記
録することなく、迅速で精度の高い位置合せをすること
ができる。

【００１６】

【発明が解決しようとする課題】しかしながら、上述し
た位置合せ方法においてはテンプレートマッチングによ
り各放射線画像の対応点を求めるようにしているが、例
えば、放射線画像中のノイズを対応点として求めてしま
う場合があり、このような対応点は各放射線画像の位置
的に共通した対応点とならず位置合せを精度良く行うこ
とができなかった。また、対応点が放射線画像のエッジ
上等の相関値の分布が尾根状となる位置に求められた場
合、位置合せのときにこの尾根に沿って画像がずれ、そ
の結果として位置合せを精度良く行うことができない場
合があった。

【００１７】本発明は上記事情に鑑み、精度良く位置合
せを行うために、位置合せに適した対応点を求めること
ができる放射線画像の位置合せ用の対応点を決定する方
法を提供することを目的とするものである。

【００１８】

【課題を解決するための手段】本発明による放射線画像
の位置合わせ用対応点決定方法は、複数の放射線画像の
位置合せを行うための対応点を決定する方法において、
前記放射線画像のうちの１つの放射線画像上にテンプレ
ート領域を設定し、該１つの放射線画像以外の他の放射
線画像上において前記テンプレート領域を移動させて該
テンプレート領域と前記他の放射線画像の一部との相関
値を求めるテンプレートマッチングを行うことにより、
前記各放射線画像の互いに対応する複数の対応点を求
め、前記各放射線画像のなかの１つの放射線画像の前記
各対応点を基準対応点とし、前記複数の対応点につい
て、（１）前記テンプレートマッチングにより得られ
た、前記テンプレート領域と前記他の放射線画像の一部
との前記相関値が所定値以上であること、（２）前記複
数の対応点近傍の前記相関値の分布が尾根状とならない
こと、（３）任意の２基準対応点間の距離と、該２基準
対応点に対応する他の放射線画像上の２対応点の距離と
の比が一定値の範囲内にあること、の条件のうち少なく
とも１つの条件を満たすか否かの評価を行い、該評価に
適合する対応点を前記位置合せのための対応点とするこ
とを特徴とするものである。

【００１９】また、前記評価に適合する対応点が所定数
未満である場合には、前記評価に適合しなかった対応点
から所定距離以上離れた放射線画像上において、前記テ
ンプレート領域の設定、前記テンプレートマッチング、
および前記評価を、該評価に適合する対応点が前記所定
数となるまで繰り返し行うことが好ましい。

【００２０】ここで、相関値の分布が尾根状になると
は、例えば放射線画像のエッジの部分を含むテンプレー
ト領域を設定した場合、このテンプレート領域と他の放
射線画像のエッジ付近におけるテンプレートマッチング
により得られる相関値は、この相関値の分布を３次元状
に見た場合、図８に示すように尾根状となる。このよう
な状態を相関値の分布が尾根状である状態と言うことと
する。

【００２１】またここで、テンプレートマッチングと
は、上述したように１の放射線画像上にテンプレート領
域を設定した場合、このテンプレート領域を他の放射線
画像上で動かすことにより最もマッチングのとれる場所
を探し出す処理で、その場所を表す点が対応点の座標を
与えるものである。

【００２２】このようなテンプレートマッチングにおい
て、そのマッチング度を表す評価尺度には、相関法およ
びＳＳＤＡ（Sequential Similarity Detection Algoli
thms）が挙げられる。

【００２３】この相関法とは、対応する各画素ごとに積
を算出し、その積の和を標準化した値（以下標準化値と
称する）を重ね合せの尺度とするものである。この標準
化は、それぞれの領域に於いて画素自身の積（２乗）の
和を算出してさらにそれぞれの和の積を算出し、この積
の平方根を対応する各画素ごとの積の和の分母とするこ
とにより行われる。重ね合せが完全な場合、雑音（ノイ
ズ）などによって分子の積が全て２乗の和とはならず、
このため標準化値は１にはならなくても１に最も近い最
大値になると考えられる。よって、テンプレート領域を
放射線画像上でいろいろ移動させて、上述した標準化値
が最大になる移動をもって重ね合せが達成されたと考え
られる。しかしながら、この標準化値が最大となる移動
は、全ての移動が終了しなければ判定することができな
い。この方法の詳細は、例えばSmith らの「Automated
cloud tracking using precisely aligned digital ATS
pictures 」ibid．、1972年 7月c-21巻、715-729 頁に
記載されている。

【００２４】また、ＳＳＤＡとは、各画素ごとに差の絶
対値の和（残差）を重ね合せの尺度とするものである。
重ね合せが完全な場合、雑音（ノイズ）などによって残
差は０にはならなくても最小にはなると考えられる。よ
って、テンプレート領域を画像上でいろいろ移動させ
て、残差が最小になる移動をもって重ね合せが達成され
たと考える。この際、重ね合せがズレていると、各画素
について順次に加算していくとき残差が急激に増大す
る。そこで加算の途中で残差があるしきい値を超えたら
早々に加算を打ち切り次の移動に移る方法がこのＳＳＤ
Ａである。用いる計算は加算だけであり、しかも多くの
場合途中で打ち切られるため、大幅に計算時間が短縮さ
れる。この方法の詳細は、例えばBarneaらの「A class
of algorithms for fast digital image registration
」IEEE．Trans ．、1972年 2月c-21巻、179-186 頁に
記載されている。

【００２５】

【作用】本発明による放射線画像の位置合せ用の対応点
を決定する方法は、前述した（１）の条件を満たすか否
かの評価により、放射線画像中のノイズのように位置合
せに適していない対応点は排除され、（２）の条件を満
たすか否かの評価により例えば放射線画像のエッジ部分
に求められた対応点のようにエッジに沿ってずれ易い対
応点は排除される。さらに、放射線画像が完全に一致す
るような場合は、各放射線画像における互いに対応する
２対応点間の距離は一定値となるため、（３）の条件を
満たすか否かの評価により各放射線画像に共通しない対
応点は排除されることとなる。このように、３つの条件
のうちのいずれか１の条件を満たすか否かの評価を行う
ことにより、位置合せに不適な対応点は排除され、より
位置合せに適した対応点のみが選択されることとなる。

【００２６】なお、複数の記録媒体を重ね合せて撮影を
行った場合に、放射線源に近い側にある放射線画像上に
ある任意の２対応点間の距離と、この２対応点に対応す
る放射線源に遠い側にある２対応点との距離とでは、後
者には放射線が広がって照射されるため後者の方が大き
いものとなる。したがって、（３）の評価の際に前者と
後者との比を一定値の範囲内とすることにより、複数の
記録媒体を重ね合せて撮影を行うことにより複数の放射
線画像を得た場合であっても、位置合せに適した対応点
を精度良く求めることができる。

【００２７】また、（１）、（２）および（３）の条件
のうちいずれの条件も満たさないため、対応点の数が少
な過ぎて、位置合せを正確に行うことができない場合が
生じることがある。このような場合、すなわち、上記評
価に適合する対応点が所定数未満である場合は、評価に
適合しなかった対応点から所定距離以上離れた放射線画
像上において、テンプレート領域の設定、テンプレート
マッチングおよび評価を、評価に適合する対応点が所定
数以上になるまで繰り返し行うことにより、位置合せの
ための対応点をより多く求めることができるため、この
対応点を用いることにより位置合せをより正確に行うこ
とができる。

【００２８】ここで、評価に適合しなかった対応点から
所定距離以上離れた放射線画像上において上記評価等を
繰り返し行うのは、例えば評価に適合しなかった対応点
において、相関値の分布が尾根状となるものであった場
合、この対応点の近傍で上記評価等を再度行うと、求め
られた対応点において相関値の分布が尾根状となってし
まい、上記（２）の条件に適合しないことがある。この
ため、評価に適合しなかった対応点から所定距離以上離
れた放射線画像上において、再度上記評価等を繰り返し
行うようにしたものである。

【００２９】

【実施例】以下、図面を参照して本発明の実施例につい
て詳細に説明する。

【００３０】図１は、本発明に用いられる放射線画像を
記録する装置の一実施例である放射線撮影装置の概略図
である。この撮影により得られる放射線画像は、エネル
ギーサブトラクション処理に用いられるものである。

【００３１】フィルタ６を間に挟み蓄積性蛍光体シート
５および７がシート７を下にして重ねられている。この
上には、被写体４を介して放射線３を発する放射線源２
が配置されている。以上のように、放射線撮影装置１が
構成されている。

【００３２】この放射線源２から発せられた放射線３が
被写体４に照射される。被写体４を透過した放射線3aは
第１の蓄積性蛍光体シート５に照射され、放射線3aのエ
ネルギーの一部が第１の蓄積性蛍光体シート５に記録さ
れ、これによりシート５に被写体４の放射線画像が蓄積
記録される。シート５を透過した放射線3bはさらにフィ
ルタ６を透過し、フィルタ６を透過した放射線3cが第２
の蓄積性蛍光体シート７に照射される。これによりシー
ト７にも被写体４の放射線画像が蓄積記録される。

【００３３】図２は、各蓄積性蛍光体シート５および７
に蓄積記録された放射線画像を模式的に表した図であ
る。各蓄積性蛍光体シート５および７の略全面に被写体
４の各放射線画像4a、4bが蓄積記録されている。すなわ
ち、放射線画像4aは上側の蓄積性蛍光体シート５から、
放射線画像4bは下側の蓄積性蛍光体シート７から得られ
た放射線画像となる。

【００３４】図３は、本発明に用いる放射線画像を読取
る読取ユニットの一実施例である放射線画像読取装置お
よび本発明の位置合せ方法を実施し、サブトラクション
処理を行う演算ユニットの一実施例である画像処理表示
装置の斜視図である。

【００３５】図１に示す放射線撮影装置１で撮影が行わ
れた後、第１および第２の蓄積性蛍光体シート５、７が
１枚ずつ放射線画像読取装置10の所定位置にセットされ
る。ここでは、第１の蓄積性蛍光体シート５に蓄積記録
された第１の放射線画像の読取りの場合について説明す
る。

【００３６】所定位置にセットされた、第１の放射線画
像が蓄積記録された蓄積性蛍光体シート５は、図示しな
い駆動手段により駆動されるエンドレスベルト等のシー
ト搬送手段15により、矢印Ｙ方向に搬送（副走査）され
る。一方、レーザ光源16から発せられた光ビーム17はモ
ータ18により駆動され矢印Ｚ方向に高速回転する回転多
面鏡19によって反射偏向され、ｆθレンズ等の集束レン
ズ20を透過した後、ミラー21により光路を変えて蓄積性
蛍光体シート５に入射し、副走査の方向（矢印Ｙ方向）
と略垂直な矢印Ｘ方向に主走査する。蓄積性蛍光体シー
ト５の光ビーム17が照射されたか所からは、蓄積記録さ
れている放射線画像情報に応じた光量の輝尽発光光22が
発せられ、この輝尽発光光22は光ガイド23によって導か
れ、フォトマルチプライヤ（光電子増倍管）24によって
光電的に検出される。光ガイド23はアクリル板等の導光
性材料を成形して作られたものであり、直線状をなす入
射端面23a が蓄積性蛍光体シート５上の主走査線にそっ
て延びるように配され、円環状に形成された射出端面23
b にフォトマルチプライヤ24の受光面が結合されてい
る。入射端面23a から光ガイド23内に入射した輝尽発光
光22は、該光ガイド23の内部を全反射を繰り返して進
み、射出端面23b から射出してフォトマルチプライヤ24
に受光され、放射線画像を表す輝尽発光光22がフォトマ
ルチプライヤ24によって電気信号に変換される。フォト
マルチプライヤ24から出力されたアナログ信号Ｓは、ロ
グアンプ25で対数的に増幅された後、Ａ／Ｄ変換器26に
入力され、サンプリングされてデジタルの画像信号ＳＯ
が得られる。この画像信号ＳＯは第１の蓄積性蛍光体シ
ート５に蓄積記録された第１の放射線画像を表すもので
あり、第１の画像信号ＳＯ₁ と呼ぶ。この第１の画像信
号ＳＯ₁ は画像処理表示装置30内の内部メモリーに一旦
記録される。

【００３７】この画像処理表示装置30は、種々の指示を
入力するキーボード31、指示のための補助情報や画像信
号に基づく可視画像を表示するＣＲＴディスプレイ32、
補助記憶媒体としてのフロッピーディスクが装填され駆
動されるフロッピィディスク駆動装置33、およびＣＰＵ
や内部メモリが内蔵された本体部34が備えられている。

【００３８】次に上記と同様にして、第２の蓄積性蛍光
体シート７に蓄積記録された第２の放射線画像を表す第
２の画像信号ＳＯ₂ が得られ、この第２の画像信号ＳＯ
₂ も画像処理表示装置30内の内部メモリに一旦記憶され
る。

【００３９】このようにしてサブトラクション演算を行
うべき２つの画像信号ＳＯ₁ ，ＳＯ₂ が内部メモリに記
憶されると、これら２つの画像信号ＳＯ₁ ，ＳＯ₂ が読
み出されて、これら２つの画像信号ＳＯ₁ ，ＳＯ₂ が担
持する各放射線画像の各画素間で対応したサブトラクシ
ョン演算が行われるように、画像の位置合せが行われ
る。

【００４０】ここで、本実施例における、画像信号ＳＯ
₁ ，ＳＯ₂ が表す２つの放射線画像の位置合せ方法につ
いて説明する。

【００４１】本発明による放射線画像の位置合せ方法に
おいては、図１における放射線源２に近い位置にある蓄
積性蛍光体シート５より得られた放射線画像4aを図４に
示すように４等分して４つの領域14Ａ，14Ｂ，14Ｃ，14
Ｄを定める。そして各領域について、クロス型のフィル
タを用いて画像をフィルタリング処理してフィルタの出
力最大値を与える点を検出する。

【００４２】ここで、フィルタのマトリクスをＡij（i=
1,2,……n 、j=1,2,……n ）としたとき、フィルタを表
す式は

【００４３】

【数１】

【００４４】となるが、クロス型フィルタとは、要素ａ
＝Ａi,j(i=j or i=n-j-1),ｂ＝Ａi,j(i ≠j and i ≠n-
j-1)としたときにａ≠ｂ（ａ＞ｂ）であるようなマトリ
クスをいうものであり、本実施例においては、例えば

【００４５】

【数２】

【００４６】というマトリクスを用いる。

【００４７】このようなクロス型のフィルタを用いるこ
とによって、図４に示すように、ろっ骨とろっ骨とが交
わるクロスエッジのような構造の複雑な部分、すなわ
ち、その周囲の領域からみて極端に濃度が変化する点を
検出することができ、本実施例においては４つの点40
Ａ，40Ｂ，40Ｃ，40Ｄを検出することができる。なお、
このようなクロスエッジは４つの領域14Ａ〜14Ｄ中に複
数存在するが、本実施例においては各領域においてフィ
ルタの出力が最大値となった点40Ａ，40Ｂ，40Ｃ，40Ｄ
を検出するものとする。

【００４８】このようにして４つの点40Ａ〜40Ｄが検出
されると、この４つの点を基準対応点とし、各基準対応
点を中心としたテンプレート領域を放射線画像4a上に設
定する。すなわち、図５に示すように４等分した各領域
14Ａ〜14Ｄに基準対応点40Ａ〜40Ｄを中心としたテンプ
レート領域41Ａ〜41Ｄを設定する。

【００４９】次いで各テンプレート領域41Ａ〜41Ｄを放
射線画像4b上において所定範囲内で移動させてテンプレ
ートマッチングを行う。ここで、テンプレートマッチン
グは上述した相関法またはＳＳＤＡを用いて行う。相関
法においては前述したように標準化値が最大となる点が
以下に記載する対応点の座標を与える。また、ＳＳＤＡ
においても上述したように残差の和が最小となる点が対
応点の座標を与える。

【００５０】このようにしてテンプレートマッチングを
行うことにより、図６に示すように放射線画像4bに複数
の対応点が求められる。

【００５１】次いで、このようにして求められた複数の
対応点について、以下の条件を満たすか否かの評価を行
う。すなわち、（１）テンプレートマッチングにより得
られた、テンプレート領域41Ａ〜41Ｄと放射線画像4bの
一部との相関値が所定値以上であること、（２）対応点
近傍の相関値の分布が尾根状とならないこと、（３）任
意の２基準対応点間の距離と、２基準対応点に対応する
放射線画像4b上の２対応点の距離との比が一定値の範囲
内にあること、の評価を行う。以下この評価について詳
細に説明する。

【００５２】まず、（１）の評価について説明する。

【００５３】図６に示す放射線画像4b上において求めら
れた複数の対応点について、各対応点におけるテンプレ
ートマッチングを行った際の相関値を求め、各対応点に
おける相関値の値が0.84未満となる対応点を排除する。
すなわち、各対応点における相関値をＣＯＲＲ（ｉ，
ｊ）（ｉ，ｊは対応点の座標）としたときに、 ＣＯＲＲ（ｉ，ｊ）＜0.84 …(3) となった対応点を図６に示す複数の対応点から排除す
る。これにより放射線画像4b中に求められた複数の対応
点のうち、ノイズのように位置合せに不適な対応点が排
除される。

【００５４】次いで、前述した評価（２）について説明
する。

【００５５】まず、図６に示す放射線画像4b上において
求められた複数の対応点について、各対応点を中心とし
て上下左右および45度斜め方向の８点における相関値の
傾きを求める。例えば、図７に示す対応点である画素50
を中心として８つの画素51Ａ〜51Ｈにおけるテンプレー
トマッチングの際の相関値の傾きを求める。

【００５６】すなわち、 Ａｋ＝（ＣＯＲＲ（ｉ，ｊ）−ＣＯＲＲｋ（ｉ，ｊ））／10 …(4) 但し、ＣＯＲＲ（ｉ，ｊ）：対応点における相関値 ＣＯＲＲｋ（ｉ，ｊ）：ｋ方向における10画素先の相関
値 Ａｋ：ｋ方向における傾き により傾きを求め、角傾きＡｋについて Ａｋ＜Ｔｈ（Ｔｈ＝0.00075 ） …(5) ならば、その画素を図６に示す複数の対応点から排除す
る。

【００５７】このように、ある方向について相関値の傾
きが小さい場合、その画素は図８に示すような相関値の
分布が尾根状となっている部分に存在するものと判断さ
れる。このように、相関値の分布が尾根状となっている
部分に対応点が存在すると、位置合せを行ったときこの
尾根に沿って画像がずれてしまうこととなる。したがっ
て、図８に示すようにある画素52の近傍の相関値の分布
が尾根状となり、対応点の位置が矢印Ａ方向にずれ易い
ような場合は、（２）の条件を満たすか否かの評価によ
りその対応点は位置合せのためには不適なものとして、
複数の対応点から排除される。

【００５８】次いで、前述した評価（３）について説明
する。

【００５９】まず、放射線画像4a上の４つの基準対応点
40Ａ〜40Ｄのうちの任意の２点間の距離と、この２点に
対応する放射線画像4b上の２対応点間の距離とを求め
る。例えば、図９に示すように、放射線画像4a上におけ
る基準対応点40Ａ，40Ｂ間の距離ＡＢと、放射線画像4b
上における対応点40Ａ′，40Ｂ′間の距離Ａ′Ｂ′を求
める。すなわち、基準対応点40Ａ，40Ｂの座標を（Ｕ
ｉ，Ｖｉ），（Ｕｊ，Ｖｊ）とし、対応点40Ａ′，40
Ｂ′の座標を（Ｘｉ，Ｙｉ），（Ｘｊ，Ｙｊ）としたと
きに

【００６０】

【数３】

【００６１】により距離ＡＢ，Ａ′Ｂ′を求める。

【００６２】次いで距離ＡＢ，Ａ′Ｂ′の比を求め、こ
の比が所定のしきい値の範囲内にあるか否かをみる。す
なわち Ｔｈ１＜Ａ′Ｂ′／ＡＢ＜Ｔｈ２ …(7) なる評価を行い、式(7) の条件を満たす対応点のみが位
置合せに適した対応点として用いられる。

【００６３】ここで、式(7) における比Ａ′Ｂ′／ＡＢ
は各点の位置が真に対応する場合、本来ならば１となる
はずであるが、図１に示す放射線画像の撮影の際には、
蓄積性蛍光体シート７に蓄積記録される放射線画像は、
蓄積性蛍光体シート５に蓄積記録される放射線画像と比
較して若干大きい画像となる。これは、蓄積性蛍光体シ
ート７は放射線源２からみてシート５よりも離れた位置
にあるため、放射線が広がって照射されることにより生
ずる現象である。したがって、式(7) におけるＡ′Ｂ′
／ＡＢの条件を１と定めずに、１より若干大きいしきい
値Ｔｈ１からＴｈ２の範囲にあるようにして、条件
（３）の評価を行うことが望ましい。

【００６４】このように、図６に示す複数の対応点の評
価を行い、最終的に位置合せに適した対応点として図10
に示すように最大４つの対応点44Ａ〜44Ｄが４つの領域
43Ａ〜43Ｄ内において選択される。

【００６５】次いで、放射線画像4a上における各基準対
応点40Ａ〜40Ｄの座標を（ｕi ，ｖi ）(i=1〜4)とし
て、アフィン変換

【００６６】

【数４】

【００６７】但し、ｍは拡大，縮小率、θは回転量、
Ａ，Ｂは回転移動補正および拡大または縮小率補正を示
す係数、Ｃ，Ｄは平行移動補正を示す係数にしたがっ
て、各対応点（ｘi ，ｙi ）の座標を変換することによ
り第１の放射線画像4aと第２の放射線画像4bとを重ね合
せる。ここで式(3) に基づく座標変換では、第２の放射
線画像全体をＸ方向とＹ方向とで互いに独立に拡大もし
くは縮小すること、第２の放射線画像全体を回転移動す
ること、および第２の放射線画像をＸ方向，Ｙ方向に平
行移動すること、のすべてが同時に行われる。

【００６８】ここで、式(8) に含まれる係数Ａ，Ｂ，
Ｃ，Ｄの求め方について説明する。

【００６９】本発明においては、最小二乗法により係数
Ａ，Ｂ，Ｃ，Ｄを決定するものであり、まず、基準対応
点（ｕi ，ｖi ）と対応点（ｘi ，ｙi ）との位置的な
誤差Ｅを、 Ｅ＝Σ（ｕ−ｕi ）² ＋Σ（ｖ−ｖi ）² ＝Σ（Ａ・ｘi −Ｂ・ｙi ＋Ｃ−ｕi ）² ＋Σ（Ｂ・ｘi ＋Ａ・ｙi ＋Ｄ−ｖi ）² …(9) とし、誤差Ｅの二乗誤差を最小にするために式(9) を解
くと、

【００７０】

【数５】

【００７１】となる。よって、式(10)から１次方程式を
解くと、係数Ａ，Ｂ，Ｃ，Ｄは以下のように求められ
る。

【００７２】 Ａ＝［ｄ・（ｅ＋ｆ）−ｂ・i −ｃ・j ］／Δ Ｂ＝［ｄ・（ｇ＋ｈ）＋ｃ・ｉ−ｂ・ｊ］／Δ Ｃ＝［−ｂ・（ｅ＋ｆ）＋ｃ・（−ｇ＋ｈ）＋ａ・ｉ］／Δ Ｄ＝［−ｃ・（ｅ＋ｆ）−ｂ・（−ｇ＋ｈ）＋ａ・ｊ］／Δ 但し、Δ＝ａ・ｄ−ｂ² −ｃ² ａ＝Σ（ｘi ² ＋ｙi ² ），ｂ＝Σｘi ，ｃ＝Σｙi ，ｄ＝Σ ｅ＝Σｘi ・ｕi ，ｆ＝Σｙi ・ｖi ， ｇ＝Σｙi ・ｕi ，ｈ＝Σｘi ・ｖi ，ｉ＝Σｕi ，ｊ＝Σｖi …(1 1) このようにして求められた係数、Ａ，Ｂ，Ｃ，Ｄを用い
て式(8) により座標変換を行うことにより、各対応点44
Ａ〜44Ｄと基準対応点40Ａ〜40Ｄとの誤差が最小となる
ように第２の放射線画像4bを第１の放射線画像4aに略等
しく合せることができる。すなわち、図11に示すよう
に、本発明の位置合せ方法により55組の放射線画像の位
置合せを行ったデータでみてみると、基準レベルを｜Ｒ
｜＝4.5 、max｜Ｒi ｜＝1.00とした場合、すべてのデ
ータが基準レベル以内、しかも｜Ｒ｜＝3.15、 max｜Ｒ
i ｜＝0.7 内に入っており、図12に示す従来の方法によ
る位置合せ結果と比較して、非常に高精度に位置合せが
行われていることがわかる。

【００７３】このようにしてアフィン変換が行われた
後、サブトラクション処理、すなわち画像信号ＳＯ₂ の
位置合せ後の画像信号をＳＯ₂ ′としたとき、 Ｓ１＝Ｗａ・ＳＯ₁ −Ｗｂ・ＳＯ₂ ′＋Ｃ …(12) 但し、Ｗａ，Ｗｂは重み付け係数、Ｃはバイアス分を表
す。

【００７４】により重み付け引き算が行われ、これによ
り２つの放射線画像の差の画像に対応する画像信号Ｓ１
が生成される。この画像信号Ｓ１は画像処理表示装置30
のＣＲＴディスプレイ32に送られ、この画像信号Ｓ１に
基づく可視画像（エネルギーサブトラクション画像）が
ＣＲＴディスプレイ32に再生表示される。なお、本体部
34で実行される上記サブトラクション処理を行う機能
（ハードウェアとソフトウェアとの組合せ）が本発明の
演算ユニットの一例と観念される。

【００７５】なお、上記実施例においては、図１におけ
る放射線源２からみて遠い位置にある蓄積性蛍光体シー
ト７より得られる放射線画像からテンプレートマッチン
グにより検出された対応点をアフィン変換するようにし
ているが、これは以下の理由によるものである。すなわ
ち、図１に示すようないわゆるワンショットで２枚の蓄
積性蛍光体シートに放射線画像を蓄積記録するような場
合は、２枚のシート５，７はフィルタ６を挟んで空間的
に互いに異なる位置に配置されることとなり、これによ
り２枚のシートの放射線源２および被写体４からの距離
が異なり、このため各シート５，７に記録される放射線
画像の寸法が異なり放射線源２からみて遠い位置にある
シート７に蓄積記録される放射線画像4bはシート５に蓄
積記録される放射線画像4aと比較してボケたり、散乱線
が多い画像となる。また、上述したアフィン変換を行っ
た場合、変換された画像については、各画素間は何らか
の補間をする必要があるため多少なりとも画質が劣化す
るものである。したがって、画質のよい放射線画像4aを
アフィン変換し画質を劣化させるよりも、放射線画像4a
と比較して画質が劣る放射線画像4bを変換した方が、結
果として得られるサブトラクション画像の画質が担保さ
れることとなる。したがって、放射線源２からみて遠い
位置にある蓄積性蛍光体シート７から得られた放射線画
像4bをアフィン変換することとしたものである。

【００７６】また、上述した実施例においては、前述し
た３つの条件（１），（２）および（３）のすべての条
件を満たした対応点の位置合せのために用いるようにし
ているが、３つの条件のうち少なくとも１つの条件を満
たす対応点を位置合せのために用いるようにしてもよ
い。但し、条件が多いほどより位置合せに適した対応点
が選択されるため、条件の数は求められる位置合せの精
度により決定すればよい。

【００７７】さらに、上述した実施例においては、条件
（１）における相関値のしきい値として0.84を用いてい
るが、このしきい値はいかなる値でもよく、さらに、図
13に示すように放射線画像の線量が大きい部分ほどしき
い値を大きくするようにしてもよい。このように線量に
応じてしきい値を変化させるのは、放射線量が低いとノ
イズが多くなりフィルタの出力値が小さくなって特徴点
を求め難くなるからである。

【００７８】また、上述した実施例においては、条件
（２）についての評価を行う際に、対応点である画素50
の上下左右および斜め45度方向の８つの画素51Ａ〜51Ｈ
における相関値の傾きを求めるようにしているが、画素
50の上下左右方向の４つの画素51Ａ，51Ｃ，51Ｅ，51Ｇ
における相関値の傾きを求めるようにしてもよいもので
ある。但し、この方向の数は多いほどより位置合せに適
した対応点が選択されるため、相関値の傾きを求める画
素の数は、求められる位置合せの精度により決定すれば
よい。

【００７９】なお、上述した条件（２）に関する評価を
行うことにより、例えば照射野絞りを用いて放射線画像
の撮影を行った場合のように、照射野エッジ上に対応点
が求められても、この対応点は排除されるため、照射野
エッジに沿って画像がずれて位置合せがなされることは
なくなる。

【００８０】また、上述した実施例においては、対応点
が選択された後に、アフィン変換により位置合せを行う
ようにしているが、位置合せの方法はアフィン変換に限
定されるものではなく、前述した特開昭58-163338 号に
開示されたマーカーを用いて位置合せを行う方法におい
て、マーカーの代わりに対応点を用いて位置ズレを補正
するようにしてもよい。

【００８１】さらに、上述した実施例において、評価を
行った結果、対応点が２個以下となってしまった場合、
すなわち、精度良く位置合せを行うには対応点が少なす
ぎる場合は、前述した評価における条件をゆるくして再
度評価を行って対応点を求めるようにすればよい。

【００８２】またこのような場合、上記（１）、（２）
および（３）の条件のうちいずれの条件をも満たさない
対応点から所定距離以上離れた位置、すなわち図14に示
すように評価に適合しない対応点を含む領域の、条件を
満たさない対応点から半径ｒ以上離れた位置において、
対応点が位置合せを正確に行うことができる所定数以上
となるまで、上述したテンプレート領域の設定、テンプ
レートマッチングおよび評価を繰り返し行うようにすれ
ばよい。

【００８３】すなわち、図14に示すように領域14A で求
められた基準対応点40A が評価に適合しない対応点であ
った場合、この対応点40A から半径ｒ離れた斜線の領域
において、放射線画像4aに対して上述したフィルタリン
グ処理を行ってテンプレート領域を設定し、この新たに
設定されたテンプレート領域を放射線画像4b上において
所定範囲内で移動させてテンプレートマッチングを行
う。そしてこのテンプレートマッチングにより新たな対
応点が求められた後、この対応点について上述した
（１）、（２）および（３）の条件のうち少なくとも１
つの条件を満たすか否かの評価を行う。そしてこれによ
り、位置合せを行う対応点の数が所定数以上（例えば３
個以上）となった場合は、この対応点を用いて位置合せ
を行う。一方、位置合せを行う対応点の数が所定数以上
とならない場合は、位置合せを行う対応点の数が所定数
以上となるまで、上記テンプレート領域の設定、テンプ
レートマッチングおよび評価を繰り返し行う。このよう
に、対応点の数が所定数以上となるまでテンプレート領
域の設定等を繰り返し行うことにより、位置合せを行う
対応点を確実に複数個求めることができ、このように複
数個の対応点を用いて位置合せを行うことにより、より
精度良く位置合せを行うことが可能となる。

【００８４】ここで、評価に適合しなかった対応点から
半径ｒ以上離れた放射線画像上において上記評価等を繰
り返し行うのは、例えば評価に適合しなかった対応点に
おいて、相関値の分布が尾根状となるものであった場
合、この対応点の近傍で上記評価等を再度行うと、求め
られた新たな対応点の近傍の領域において相関値の分布
がまた尾根状となってしまうことがあり、上記（２）の
条件に適合しないことがある。このため、評価に適合し
なかった対応点から半径ｒ以上離れた放射線画像上にお
いて、再度上記評価等を繰り返し行うようにしたもので
ある。なお、この半径ｒはこの尾根状となる部分を避け
ることができるように設定することが必要である。

【００８５】なお、このように評価を繰り返し行う場合
には、演算時間短縮のため繰り返し回数を制限すること
が好ましい。本出願人の胸部の放射線画像を対象とした
実験によれば、約500 例のデータベースにおいて位置合
せを行うための対応点の数を４個とした場合、平均的な
対応点の数は繰り返しを行わない場合は２．８個、繰り
返し回数を３回に制限した場合は３．８個であった。し
たがって、繰り返し回数を制限する場合は３回程度に制
限することが好ましい。

【００８６】また、上述した実施例においては、エネル
ギーサブトラクション処理を行うために２つの放射線画
像の位置合せを行うようにしているが、重ね合せ処理を
行うために２つの放射線画像の位置合せを行うようにし
てもよい。すなわち、図15に示すように図２と同様の放
射線画像撮影装置１においてフィルタ６を用いることな
く２枚の蓄積性蛍光体シート５′，７′に被写体４の放
射線画像を蓄積記録し、蓄積性蛍光体シート５′，７′
から図３に示す放射線画像読取装置により２つの放射線
画像を表す画像信号ＳＯ₁ ，ＳＯ₂ を得、２枚の蓄積性
蛍光体シート５′，７′から得られた放射線画像のう
ち、シート７′から得られた放射線画像を上述した実施
例と同様にアフィン変換を施して位置合せを行った後、
重ね合せ処理、すなわち画像信号ＳＯ₂ の位置合せ後の
画像信号ＳＯ₂ ′としたとき、 Ｓ２＝Ｗｃ・ＳＯ₁ ＋Ｗｄ・ＳＯ₂ ′ …(13) 但し、Ｗｃ，Ｗｄは重み付け係数 により、重み付け加算が行われ、これにより２つの放射
線画像の和の画像に対応する画像信号Ｓ２を生成するも
のである。

【００８７】また、上述した実施例においては、クロス
型フィルタにより放射線画像をフィルタリング処理し
て、フィルタ出力の最大値を与える点を検出するように
しているが、この点としては例えば、分割された各領域
の中の最大値を表す点、ある所定値以上の点等、いかな
る点を用いてもよいものである。

【００８８】さらに、上述した実施例においては、２つ
の放射線画像の位置合せについて説明しているが、位置
合せを行う放射線画像の数は３枚以上であってもよいも
のである。この場合、３以上の放射線画像のうち、１つ
の放射線画像に上述した基準対応点を設定し、他の放射
線画像についてテンプレートマッチングにより対応点を
設定し、各放射線画像の対応点をアフィン変換により基
準対応点に合わせるようにすればよい。

【００８９】また、上述した実施例においては式(8) に
示すアフィン変換により位置合せを行うようにしている
が、一般式である

【００９０】

【数６】

【００９１】（但し、ｕ，ｖは基準対応点の座標、ｘ，
ｙは変換される対応点の座標、ａ，ｂ，ｃ，ｄは回転移
動補正および拡大または縮小率補正を示す係数、ｅ，ｆ
は平行移動補正を示す係数）により表されるアフィン変
換を行うものであれば、いかなる係数により位置合せを
行うようにしてもよいものである。

【００９２】

【発明の効果】以上詳細に説明したように、本発明によ
る放射線画像の位置合せ用対応点決定方法は、テンプレ
ートマッチングにより求められた複数の対応点について
前述した（１）、（２）および（３）のうち少なくとも
１つの評価をするようにしたため、位置合せに適した対
応点をより精度良く求めることができ、このように求め
られた対応点を用いて位置合せを行えば、より精度良く
迅速に位置合せを行うことができる。

【００９３】また、上記評価に適合する対応点が所定数
未満であるときに、テンプレート領域の設定、テンプレ
ートマッチングおよび評価を、評価に適合する対応点が
所定数となるまで繰り返し行うことにより、位置合せの
ための対応点をより多く求めることができるため、この
対応点を用いることにより位置合せをより正確に行うこ
とができる。

【図面の簡単な説明】

【図１】本発明の実施例によるサブトラクションを行う
放射線画像を得る放射線画像記録装置の概略図

【図２】各蓄積性蛍光体シートに蓄積記録された放射線
画像を模式的に表す図

【図３】本発明による放射線画像位置合せ方法を実施す
るための装置を内包した画像読取装置の概略を表す図

【図４】放射線画像を分割して基準対応点を設定した状
態を表す図

【図５】放射線画像にテンプレート領域を設定した状態
を表す図

【図６】他の放射線画像に対応点を設定した状態を表す
図

【図７】評価（２）を説明するための図

【図８】相関値の分布が尾根状となる状態を表す図

【図９】評価（３）を説明するための図

【図１０】選択された対応点を表す図

【図１１】本発明による放射線画像位置合せ方法により
位置合せがなされた結果を表す図

【図１２】従来の放射線画像位置合せ方法により位置合
せがなされた状態を表す図

【図１３】放射線としきい値との関係を表す図

【図１４】評価に適合しなかった対応点から半径ｒ以上
離れた領域を斜線で表した図

【図１５】本発明の実施例による重ね合せを行う放射線
画像を得る放射線画像記録装置の概略図

【符号の説明】

１ 放射線撮影装置 ２ 放射線源 ３ 放射線 ４ 被写体 ５ 第１の蓄積性蛍光体シート ６ フィルタ ７ 第２の蓄積性蛍光体シート 10 放射線画像読取装置 15 シート搬送手段 16 レーザ光源 17 光ビーム 18 モータ 19 回転多面鏡 20 集束レンズ 21 ミラー 22 輝尽発光光 23 光ガイド 24 フォトマルチプライヤ 25 ログアンプ 26 Ａ／Ｄ変換器 30 画像処理表示装置 31 キーボード 32 ＣＲＴディスプレイ 33 フロッピィディスク駆動装置 34 本体部 40Ａ〜40Ｄ 基準対応点 41Ａ〜41Ｄ テンプレート領域 44Ａ〜44Ｄ 対応点

───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.⁶ 識別記号 庁内整理番号 ＦＩ 技術表示箇所 Ｇ０６Ｔ 7/00 Ｇ０６Ｆ 15/62 ３９０ Ａ 9061−5Ｈ 15/70 ３３０ Ｐ

Claims (2)

【特許請求の範囲】
1. 【請求項１】 複数の放射線画像の位置合せを行うため
   の対応点を決定する方法において、 前記放射線画像のうちの１つの放射線画像上にテンプレ
   ート領域を設定し、 該１つの放射線画像以外の他の放射線画像上において前
   記テンプレート領域を移動させて該テンプレート領域と
   前記他の放射線画像の一部との相関値を求めるテンプレ
   ートマッチングを行うことにより、前記各放射線画像の
   互いに対応する複数の対応点を求め、 前記各放射線画像のなかの１つの放射線画像の前記各対
   応点を基準対応点とし、 前記複数の対応点について、（１）前記テンプレートマ
   ッチングにより得られた、前記テンプレート領域と前記
   他の放射線画像の一部との前記相関値が所定値以上であ
   ること、（２）前記複数の対応点近傍の前記相関値の分
   布が尾根状とならないこと、（３）任意の２基準対応点
   間の距離と、該２基準対応点に対応する他の放射線画像
   上の２対応点の距離との比が一定値の範囲内にあるこ
   と、の条件のうち少なくとも１つの条件を満たすか否か
   の評価を行い、 該評価に適合する対応点を前記位置合せのための対応点
   とすることを特徴とする放射線画像の位置合せ用対応点
   決定方法。
2. 【請求項２】 前記評価に適合する対応点が所定数未満
   である場合は、前記評価に適合しなかった対応点から所
   定距離以上離れた放射線画像上において、前記テンプレ
   ート領域の設定、前記テンプレートマッチング、および
   前記評価を、該評価に適合する対応点が前記所定数とな
   るまで繰り返し行うことを特徴とする請求項１記載の放
   射線画像の位置合せ用対応点決定方法。