JPH07262370A - 画像のテンプレートマッチング方法 - Google Patents

Abstract

(57)【要約】 （修正有） 【目的】 テンプレートマッチング方法において、画像
サイズが大きくなっても計算量を大幅に増大させず、精
度良く高速にマッチングを行う。 【構成】 画像においてテンプレート41Ａを５画素ごと
に移動させてテンプレート41Ａと画像との相関値を求め
る粗サーチを行う。粗サーチにより得られた相関値が最
大となる領域を第１の基準領域51Ａ、最大値の次の値と
なる領域を第２の基準領域51Ｂと定め、基準領域51Ａ，
51を中心としてテンプレート41Ａを１画素だけ周囲にず
らして、テンプレート41Ａと画像との相関値を求める精
サーチを行う。精サーチにより得られた相関値が最大と
なる領域が基準領域51Ａ，51Ｂのいずれかであった場合
は、基準領域をテンプレート41Ａと画像43Ａとが一致す
る領域とし、精サーチにより得られた相関値が最大とな
る領域が基準領域51Ａ，51Ｂ以外の領域であった場合
は、最大となった領域を新たな基準領域として再度精サ
ーチを行い、基準領域が一致する領域となるまで精サー
チを繰り返す。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は画像のテンプレートマッ
チング方法、すなわち、ある画像の一部の領域に相当す
るテンプレートが他の画像の一部と一致する領域を探索
するテンプレートマッチング方法に関するものである。

【０００２】

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

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

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

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

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

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

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

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

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

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

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

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

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

【００１５】この方法によれば、この対応する領域の中
の点（対応点）をマーカーの代りに位置合せのために使
用することができるので、マーカー等を被写体とともに
記録することなく、迅速で精度の高い位置合せをするこ
とができる。

【００１６】このようなテンプレートマッチングを行う
方法としては、テンプレートを他の画像の全体に亘って
１画素ずつずらしながら、それぞれの位置において１画
素単位の解像度でテンプレートと画像との相関値を算出
し、この相関値が最大となる領域に含まれる対応点の位
置を求める方法が採用されていた。

【００１７】ところが、この方法ではテンプレートと画
像との相関値が最大となる対応領域の位置を求めるた
め、莫大な回数の画素間演算が必要であった。

【００１８】そこで、テンプレートマッチングを行う際
の演算量を低減させたテンプレートマッチング方法が提
案されている（特開昭63-211474 号公報）。この方法
は、複数の大きさのテンプレートを用意して、テンプレ
ートマッチングを複数段階に分割して行う階層化構造の
テンプレートマッチング方法である。すなわち、テンプ
レートの大きさに合わせて画像を間引いてテンプレート
マッチングを行い、結果として複数画素ずつテンプレー
トをずらしてテンプレートマッチングを行い、段階的に
テンプレートマッチングの解像度を高くして最終的な対
応点を求める方法である。この方法によれば、従来のテ
ンプレートマッチング方法と比較して大幅に少ない計算
量でテンプレートマッチングを行うことができる。

【００１９】

【発明が解決しようとする課題】しかしながら、上述し
た特開昭63- 211474号公報に開示された方法は、階層化
構造をしているため、テンプレートマッチングを行う画
像のサイズが大きい場合、精度良くマッチングを行うた
めには、層の数を大きくする必要があり、結果としてマ
ッチングを行うための計算量をそれ程削減できず、マッ
チングのための計算にやはり長時間を要するものとなっ
ていた。

【００２０】本発明は上記事情に鑑み、画像サイズが大
きくなっても計算量がそれ程増大せず、精度良く高速に
マッチングを行うことができるテンプレートマッチング
方法を提供することを目的とするものである。

【００２１】

【課題を解決するための手段】本発明によるテンプレー
トマッチング方法は、上記のようなテンプレートマッチ
ング方法において、先ずテンプレートを複数個の画素を
間引いて所定数の画素ごとに大きく移動させてテンプレ
ートと画像の一部との相関値を間引いた間隔の画素ごと
に求める粗サーチを行い、この粗サーチにより得られた
相関値が最大となる領域を基準領域と定め、次いでこの
基準領域を中心としてテンプレートを１画素だけその周
囲にずらして、テンプレートと画像の一部との相関値を
画素を間引かないで求める精サーチを行い、この精サー
チにより得られた相関値が最大となる領域が前記基準領
域である場合は、該基準領域を前記一致する領域とし、
該精サーチにより得られた相関値が最大となる領域が前
記基準領域以外の領域であった場合には、該最大となる
領域を新たな基準領域として再度前記精サーチを繰返し
て行い、前記基準領域が前記一致する領域となるまで前
記精サーチを繰り返すことを特徴とするものである。

【００２２】また、本発明のもう一つの方法によるテン
プレートマッチング方法は、上述したテンプレートマッ
チングを行う方法において、上述した粗サーチを行った
後に、該粗サーチにより得られた前記相関値が最大とな
る領域と次に大きい値となる領域を求めてこれを第１お
よび第２の基準領域と定め、該第１および第２の基準領
域を中心としてテンプレートを１画素だけ周囲にずらし
て、テンプレートと画像の一部との相関値を求める精サ
ーチを行い、該精サーチにより得られる相関値が最大と
なる領域が前記第１または第２の基準領域の少なくとも
一方である場合は、該第１または第２の基準領域を前記
一致する領域とし、該精サーチにより得られる相関値が
最大となる領域が前記第１または第２の基準領域以外の
領域である場合は、該最大となる領域が前記第１および
第２の基準領域からみて略同一方向に向かう側または互
いに略近付く方向に向かう側にあるか否かを検出し、該
最大となる領域が前記第１および第２の基準領域からみ
て略同一方向に向かう側または互いに略近付く方向に向
かう側にある場合は、該最大となる領域を新たな第１お
よび第２の基準領域と定め、前記第１または第２の基準
領域の少なくとも一方が前記一致する領域となるまで前
記精サーチを繰り返すことを特徴とするものである。

【００２３】なお、このテンプレートマッチング方法に
おいて、前記最大となる領域が前記第１および第２の基
準領域からみて略同一方向に向かう側または互いに略近
付く方向に向かう側にない場合は、粗サーチを行う際の
テンプレートを画像上において移動させる前記画素の間
隔を前記所定数よりも小さい数としてテンプレートマッ
チングをやり直すのが望ましい。

【００２４】また、これらのテンプレートマッチング方
法において、前記最大となる領域が前記第１および第２
の基準領域からみて互いに略近付く方向に向かう側にあ
る場合は、該第１または第２の基準領域の少なくともい
ずれか一方が前記最大となる領域となるまで前記精サー
チを繰り返すようにしてもよい。

【００２５】あるいは、この場合は、該第１および第２
の２つの基準領域が一致するまで前記精サーチを繰り返
すようにしてもよい。

【００２６】さらにこの場合は、該第１または第２の基
準領域のいずれか一方についてのみ前記最大となる領域
となるまで前記精サーチを繰り返すようにすることもで
きる。

【００２７】さらにこの場合は、該第１または第２の基
準領域のうち該同一方向に向かう側にある基準領域につ
いてのみ前記最大となる領域となるまで前記精サーチを
繰り返すようにしてもよい。

【００２８】

【作用】本発明によるテンプレートマッチング方法は、
まず粗サーチを行うことにより、１画素ずつずらしてテ
ンプレートマッチングを行う場合よりも少ない計算量で
テンプレートと他の画像との相関値が最大となる領域の
候補である基準領域を求めた後、この基準領域を中心と
してテンプレートを１画素だけ周囲にずらして相関値を
求める精サーチを行い、基準領域がテンプレートと他の
画像の一部とが一致する領域となるまで精サーチを繰り
返すようにしたものである。このように粗サーチ、精サ
ーチという２つのステップで相関値の最大値を求めるよ
うにしたため、テンプレートを画像上で１画素ずつずら
してテンプレートマッチングを行う場合よりも計算量は
少なくなり、さらに前述した階層構造のテンプレートマ
ッチング方法のように複数のテンプレートを用意して画
像の全面についてテンプレートと画像との相関値を求め
る必要もなくなる。

【００２９】また、、粗サーチのときに相関値の最大値
とともに最大値の次に大きい値となる領域を第１および
第２の基準領域とし、この２つの基準領域について前述
した精サーチを行うことにより、テンプレートと他の画
像の一部との相関値を求めるときの精度を向上させるこ
とができる。

【００３０】さらに、精サーチの際の相関値が最大値と
なる領域が第１および第２の基準領域からみて略同一方
向に向かう側または互いに略近付く方向に向かう側にな
い場合は、テンプレートと他の画像との相関値が最大と
なる領域を求めることができない。したがって、粗サー
チを行うときのテンプレートを移動させる画素の間隔を
小さくすることにより、第１および第２の基準領域をテ
ンプレートと他の画像との相関値が最大となる領域の候
補により近い領域として求めることができる。

【００３１】また、精サーチの際の相関値が最大値とな
る領域が第１および第２の基準領域からみて略同一方向
に向かう側または互いに略近付く方向に向かう側にある
場合は、第１または第２の基準領域のいずれか一方が最
大となる領域となるまで精サーチを繰り返す、あるいは
第１または第２の基準領域のいずれかについてのみ精サ
ーチを行うようにすれば、テンプレートマッチングのた
めの計算時間をより短縮できる。

【００３２】さらに、この場合において、第１および第
２の基準領域が一致するまで精サーチを行うことによ
り、テンプレートと他の画像の一部との相関値を求める
際の精度を向上させることができる。

【００３３】また、精サーチの際の相関値が最大値とな
る領域が第１および第２の基準領域からみて略同一方向
に向かう側にある場合は、第１または第２の基準領域の
いずれか一方が最大となる領域となるまで精サーチを繰
り返すようにすれば、テンプレートマッチングのための
計算時間をより短縮できる。

【００３４】

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

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

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

【００３７】この放射線源２から発せられた放射線３が
被写体４に照射される。被写体４を透過した放射線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上において所定範囲内で移動させてテンプレ
ートマッチングを行う。

【００５５】ここで、テンプレートマッチングの方法に
ついて説明する。

【００５６】まず、図５に示すテンプレート領域41Ａ〜
41Ｄを、図６に示す放射線画像4bの４等分された領域43
Ａ〜43Ｄ内でそれぞれ移動させる。この際各テンプレー
ト領域41Ａ〜41Ｄは５画素ごとに移動させる粗サーチを
行う。例えば、領域43Ａについては図７に示すように、
斜線部分の画素に、テンプレート領域41Ａの中心点を一
致させてテンプレート領域41Ａと放射線画像4bの一部の
領域43Ａとの相関値を求める。この相関値は前述した相
関法またはＳＳＤＡを用いて行う。

【００５７】このようにして粗サーチが行われた後、こ
の粗サーチにより得られた相関値が最大となる第１の基
準領域および最大の次に大きい値となった第２の基準領
域を求める。本実施例においては、図８に示す画素50Ａ
を中心とする第１の基準領域51Ａと、画素50Ｂを中心と
する第２の基準領域51Ｂとが求められる。

【００５８】第１および第２の基準領域51Ａ，51Ｂが求
められると、テンプレート領域41Ａの中心を各基準領域
51Ａ，51Ｂの中心の画素50Ａ，50Ｂの周囲に１画素だけ
ずらしてテンプレート領域41Ａと領域43Ａとの相関値を
求める精サーチを行う。例えば、第１の基準領域51Ａに
ついてみると、図９(a) に示すように、画素50Ａの周囲
をとりまく８つの画素52Ａ〜52Ｈのそれぞれとテンプレ
ート領域41Ａの中心とを一致させて８つの相関値を得
る。これを図９(b) に示すように第２の基準領域51Ｂの
周囲の画素53Ａ〜53Ｈについても同様に行い８つの相関
値を得る。

【００５９】このように各基準領域51Ａ，51Ｂについて
それぞれ８つの相関値が得られた後に、８つの相関値お
よびテンプレート領域41Ａの中心と各画素50Ａ，50Ｂと
が一致したときの相関値の比較を行い、この９つの相関
値の中から最大となる相関値を与える画素を求める。本
実施例においては図８のI-I 線断面図におけるテンプレ
ート41Ａと領域43Ａとの相関値は、図10に示すような分
布となっているため第１の基準領域51Ａにおいては画素
52Ｆとテンプレート領域41Ａの中心とが一致したときに
得られる相関値が最大となり、第２の基準領域51Ｂにお
いては画素53Ｂとテンプレート領域41Ａの中心とが一致
したときに得られる相関値が最大となる。

【００６０】このように、精サーチにより得られる相関
値が最大となる領域が基準領域51Ａ，51Ｂ以外の領域と
なる場合は、相関値が最大となる画素52Ｆ，53Ｂを中心
としたテンプレート領域41Ａと略同一の大きさの領域を
新たな第１および第２の基準領域として再度各画素52
Ｆ，53Ｂの周囲にテンプレート領域41Ａの中心を１画素
だけずらしてテンプレート領域41Ａと領域43Ａとの相関
値を求める精サーチを行う。

【００６１】そしてこのような精サーチを、第１および
第２の基準領域が一致するまで行うと、図11に示すよう
に、画素50Ａと画素50Ｂとの中間でやや画素50Ａ寄りに
ある画素54を中心としたテンプレート41Ａと同一の大き
さの領域が、第１および第２の基準領域が一致した領域
となる。

【００６２】このように、第１および第２の基準領域51
Ａ，51Ｂが一致したときに、各基準領域51Ａ，51Ｂの中
心となる画素54が本実施例において求めるべき対応点と
なる。

【００６３】このように、粗サーチ、精サーチという２
段階のステップでテンプレートマッチングを行うことに
より、従来のテンプレート領域を１画素ずつずらしてテ
ンプレートマッチングを行う場合、さらには階層化構造
のテンプレートマッチングを行う場合と比較して少ない
画素により相関値を求めることができるため、演算時間
を短縮することができ、しかも１画素ずつテンプレート
マッチングを行う場合と変わらない精度で対応点を求め
ることができる。

【００６４】以上説明したテンプレートマッチングを放
射線画像4bの４等分した領域43Ａ〜43Ｄについて行い、
図６に示すような４つの領域43Ａ〜43Ｄのそれぞれに４
つの対応点44Ａ〜44Ｄが求められる。

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

【００６６】

【数３】

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

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

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

【００７０】

【数４】

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

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

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

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

【００７５】なお、上記実施例においては、第１および
第２の基準領域51Ａ，51Ｂが一致するまで精サーチを行
うようにしているが、これに限定されるものではなく、
第１または第２の基準領域51Ａ，51Ｂのいずれか一方の
領域が、相関値が最大となった時点で精サーチを止め、
この相関値が最大となった領域の中心となる画素を前述
した対応点とするようにしてもよい。すなわち、図11に
示すように画素54が対応点となる場合、２つの画素51
Ａ，51Ｂを第１および第２の基準領域の中心として精サ
ーチを行うと、画素50Ａを中心とする第１の基準領域の
方が画素50Ｂを中心とする第２の基準領域と比較して早
い処理段階で領域の中心が画素54に一致する。すなわ
ち、第１の基準領域の中心は２回の精サーチにより画素
54と一致するが、第２の基準領域の中心は精サーチを３
回行わなければ画素54と一致しない。

【００７６】このような場合、第２の基準領域の中心を
画素54と一致させるため、３回目の精サーチは行わなく
とも対応点を設定することができるため、第１または第
２の基準領域51Ａ，51Ｂのいずれか一方の領域が、相関
値が最大となった時点で精サーチを止め、相関値が最大
となった領域の中心となる画素を対応点とすれば、演算
時間を短縮できるためより好ましい。

【００７７】また、図11に示すような場合は第１および
第２の基準領域の双方について精サーチを行う必要はな
く、相関値が最大となる領域が、第１および第２の基準
領域かにみて略近づく方向に向かう側にある場合は、第
１の基準領域または第２の基準領域のいずれか一方の領
域についてのみ精サーチを行えば、図11に示す画素54を
対応点として求めることができる。そしてこれにより、
テンプレートマッチングを行う際の計算量をさらに低減
させることができるためより好ましい。

【００７８】また、図14に示すように対応点となるべき
画素が、画素55である場合、第１および第２の基準領域
51Ａ，51Ｂは図示矢印Ａ，Ｂ方向へ移動し、各基準領域
51Ａ，51Ｂの中心が画素55と一致することとなるが、こ
の場合は前述した実施例と同様に、各基準領域51Ａ，51
Ｂが一致するまで精サーチを繰り返すようにしてもよ
い。また、１回目の精サーチにより相関値が最大となる
領域が第１および第２の基準領域からみて略同一方向に
向かう側にある場合は、演算時間短縮のために、図１に
示す矢印Ａ，Ｂ方向に向かう側にある第１の基準領域51
Ａのみを、相関値が最大となる領域となるまで精サーチ
を繰り返すようにすることが好ましい。

【００７９】さらに、上述した実施例においては、精サ
ーチを繰り返し行うことにより対応点を求めるようにし
ているが、粗サーチにより求められた基準領域の中心と
なる画素を中心として精サーチを行った結果、中心とな
る画素の周囲の画素とテンプレート領域との相関値より
も、粗サーチにより得られた相関値の方が大きい場合
は、粗サーチにより求められた基準領域の中心となる画
素が求めるべき対応点とになる。

【００８０】また、上述した実施例においては、精サー
チを繰り返すことにより第１および第２の基準領域が対
応点となるべき画素の方向へ向かって移動していくが、
例えば図15に示すように、粗サーチにより求められた第
１および第２の基準領域の中心となる画素50Ａ，50Ｂ
が、精サーチにより矢印Ｃ，Ｄに示すように全く異なる
方向（すなわち、同一方向でもなく、互いに近付く方向
でもない場合）に移動する場合がある。また、相関値が
図16に示すような分布をしている場合は、第１および第
２の基準領域の中心は略同一方向へ動くが、２つの基準
領域はそれぞれ異なる画素を中心としたときに相関値が
最大となるため、放射線画像4bの４つの領域中に唯一の
対応点を設定することができない。このような場合は、
粗サーチを行うときの画素の間隔を小さくして再度粗サ
ーチを行うようにすればよい。すなわち、図17の斜線に
示すように、テンプレート領域を３画素ごとに移動させ
ればよい。このように粗サーチの際の画素の間隔を小さ
くすることにより、図16に示すような相関値の分布を示
す場合であっても、図18に示すように相関値の最大値へ
と向かう画素を第１および第２の基準領域の中心として
求めることができることとなる。なお、再度粗サーチを
行っても相関値が最大となる唯一の画素が求められない
場合は、さらに粗サーチを行う際の画素の間隔を小さく
（例えば２画素）とすればよい。

【００８１】また、上記実施例においては、粗サーチの
結果得られた相関値が最大となる領域および最大の次に
大きい値となる領域の２つの領域を基準領域として求
め、その後精サーチを行うようにしているが、粗サーチ
により得られた相関値が最大となる１つの領域のみを基
準領域として精サーチを行うようにしてもよい。この場
合、基準領域を２つ設定して精サーチを行う場合よりも
精度は劣るが、演算時間を短縮することができる。

【００８２】また、上記実施例においては、図９に示す
ように粗サーチにより求められた基準領域の中心の画素
の周囲８画素についてテンプレートを移動させて相関値
を求めるようにしているが、計算量を低減させるために
基準領域の中心の画素の周囲４画素について精サーチを
行うようにしてもよい。すなわち、図９(a) に示す画素
52A ，52C ，52E ，52G 、図９(b) に示す画素53A ，53
C ，53E ，53G についてのみ精サーチを行うようにして
もよいものである。

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

【００８４】また、上述した実施例においては、放射線
画像上に４つの対応点を設定してアフィン変換を行うよ
うにしているが、この点の数は３点以上であればいくつ
でもよく、また対応点が多いほど位置合せの精度も向上
するものである。しかしながら、対応点の数が多いと演
算時間が長くなるため、位置合せ精度と演算時間との兼
ね合いから対応点の数を設定するのが好ましい。

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

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

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

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

【００８９】

【数５】

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

【００９１】

【発明の効果】以上詳細に説明したように、本発明によ
るテンプレートマッチング方法は、粗サーチおよび精サ
ーチという２つのステップによりテンプレートと画像の
一部との相関値を求めるようにしたため、テンプレート
を画像上で１画素ずつずらしてテンプレートマッチング
を行う場合よりも計算量は少なくなり、さらに前述した
階層構造のテンプレートマッチング方法のように複数の
テンプレートを用意して画像の全面についてテンプレー
トと画像との相関値を求める必要もなくなるため、少な
い計算量で高速かつ迅速に高精度のテンプレートマッチ
ングを行うことができる。

【図面の簡単な説明】

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

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

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

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

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

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

【図７】粗サーチを行う際の画像の一部の画素の間隔を
表す図

【図８】粗サーチにより相関値が最大となる第１の基準
領域および最大の次の値となる第２の基準領域が求めら
れた状態を表す図

【図９】精サーチにより相関値を求める際に基準領域の
中心となる画素の周囲の画素を表す図

【図１０】相関値の分布を表すグラフ

【図１１】第１の基準領域と第２の基準領域とが一致し
た状態を表す図

【図１２】本発明により位置合せがなされた結果を表す
図

【図１３】従来の方法により位置合せがなされた状態を
表す図

【図１４】第１および第２の基準領域が精サーチにより
移動する方向を表す図

【図１５】第１および第２の基準領域が精サーチにより
移動する方向を表す図

【図１６】相関値の分布を表すグラフ上で第１および第
２の基準領域が精サーチにより移動する方向を表す図

【図１７】粗サーチを行う際の画像の一部の画素の間隔
を小さくした状態を表す図

【図１８】相関値の分布を表す図

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

【符号の説明】

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

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

Claims (7)

【特許請求の範囲】
1. 【請求項１】 画像の一部に相当するテンプレートを他
   の画像の一部にマッチングさせることにより前記テンプ
   レートと前記他の画像の一部とが一致する領域を求める
   テンプレートマッチングを行う方法において、 前記他の画像上において前記テンプレートを所定数の画
   素ごとに移動させて前記テンプレートと前記他の画像の
   一部との相関値を求める粗サーチを行い、 該粗サーチにより得られた前記相関値が最大となる領域
   を基準領域と定め、 該基準領域を中心として前記テンプレートを１画素だけ
   周囲にずらして、前記テンプレートと前記他の画像の一
   部との相関値を求める精サーチを行い、 該精サーチにより得られる相関値が最大となる領域が前
   記基準領域である場合は、該基準領域を前記一致する領
   域とし、 該精サーチにより得られる相関値が最大となる領域が前
   記基準領域以外の領域である場合は、該最大となる領域
   を新たな基準領域として再度前記精サーチを行い、 前記基準領域が前記一致する領域となるまで前記精サー
   チを繰り返すことを特徴とするテンプレートマッチング
   方法。
2. 【請求項２】 画像の一部に相当するテンプレートを他
   の画像の一部にマッチングさせることにより前記テンプ
   レートと前記他の画像の一部とが一致する領域を求める
   テンプレートマッチングを行う方法において、 前記他の画像上において前記テンプレートを所定数の画
   素ごとに移動させて前記テンプレートと前記他の画像の
   一部との相関値を求める粗サーチを行い、 該粗サーチにより得られた前記相関値が最大および次に
   大きい値となる領域を第１および第２の基準領域と定
   め、 該第１および第２の基準領域を中心として前記テンプレ
   ートを１画素だけ周囲にずらして、前記テンプレートと
   前記他の画像の一部との相関値を求める精サーチを行
   い、 該精サーチにより得られる相関値が最大となる領域が前
   記第１または第２の基準領域の少なくとも一方である場
   合は、該第１または第２の基準領域を前記一致する領域
   とし、 該精サーチにより得られる相関値が最大となる領域が前
   記第１または第２の基準領域以外の領域である場合は、
   該最大となる領域が前記第１および第２の基準領域から
   みて略同一方向に向かう側または互いに略近付く方向に
   向かう側にあるか否かを検出し、 該最大となる領域が前記第１および第２の基準領域から
   みて略同一方向に向かう側または互いに略近付く方向に
   向かう側にある場合は、該最大となる領域を新たな第１
   および第２の基準領域と定め、 前記第１または第２の基準領域の少なくとも一方が前記
   一致する領域となるまで前記精サーチを繰り返すことを
   特徴とするテンプレートマッチング方法。
3. 【請求項３】 前記最大となる領域が前記第１および第
   ２の基準領域からみて略同一方向に向かう側または互い
   に略近付く方向に向かう側にない場合は、前記粗サーチ
   を行う際の前記テンプレートを前記他の画像上において
   移動させる前記画素の間隔を前記所定数よりも小さい数
   とすることを特徴とする請求項２記載のテンプレートマ
   ッチング方法。
4. 【請求項４】 前記最大となる領域が前記第１および第
   ２の基準領域からみて互いに略近付く方向に向かう側に
   ある場合は、該第１または第２の基準領域のいずれか一
   方が前記最大となる領域となるまで前記精サーチを繰り
   返すことを特徴とする請求項２または３記載のテンプレ
   ートマッチング方法。
5. 【請求項５】 前記最大となる領域が前記第１および第
   ２の基準領域からみて互いに略近付く方向に向かう側に
   ある場合は、該第１および第２の基準領域が一致するま
   で前記精サーチを繰り返すことを特徴とする請求項２ま
   たは３記載のテンプレートマッチング方法。
6. 【請求項６】 前記最大となる領域が前記第１および第
   ２の基準領域からみて互いに略近付く方向に向かう側に
   ある場合は、該第１または第２の基準領域のいずれか一
   方についてのみ前記最大となる領域となるまで前記精サ
   ーチを繰り返すことを特徴とする請求項２または３記載
   のテンプレートマッチング方法。
7. 【請求項７】 前記最大となる領域が前記第１および第
   ２の基準領域からみて略同一方向に向かう側にある場合
   は、該第１または第２の基準領域のうち該同一方向にに
   向かう側にある基準領域についてのみ前記最大となる領
   域となるまで前記精サーチを繰り返すことを特徴とする
   請求項２または３記載のテンプレートマッチング方法。