JPH11155862A - 超音波診断装置及び超音波画像処理装置 - Google Patents

Abstract

(57)【要約】 【課題】本発明の目的は、心筋等の収縮拡張機能を、客
観的、定量的、非侵襲的に、しかも簡易に高精度に提示
できる超音波診断装置及び超音波画像処理装置を提供す
ることにある。 【解決手段】本発明は、被検体を超音波で走査すること
により得られた超音波画像上で関心部位の領域を抽出
し、この領域の輪郭に応じた方向に従って組織の運動速
度を補正し、この補正された運動速度に基づいて関心部
位の収縮拡張機能を表す指標を演算するものである。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、心筋等の生体組織
の速度を推定し、推定した速度情報を処理して医学診断
に有効な情報を提供する超音波診断装置及び超音波画像
処理装置に関する。

【０００２】

【従来の技術】心臓の運動機能、つまり心機能を客観的
かつ定量的に評価することは、心臓病の診断にとって非
常に重要である。定量的評価法には、左室壁運動解析
法、ストレスエコー法などがあるが、虚血性心疾患に専
門の方法であり、汎用的ではない。

【０００３】また、左室壁運動解析法は、心機能を表す
指標として、心臓の収縮や拡張運動に伴う心筋の壁厚の
経時的変化をＢモード画像から計測するという方法であ
り、客観性に乏しい。また、左室壁運動解析法では、心
臓の長軸方向または心臓内膜あるいは外膜に沿う方向の
伸張や収縮を計測できない。

【０００４】このような問題を解決する方法として、セ
グメントショートニング法と呼ばれる方法があり、この
方法は、超音波の強反射体であるクリスタルを幾つか心
筋に直接埋め込み、強反射体間の距離変化から心機能を
比較的高精度に計測するというものであるが、この方法
は、現在のところ、動物実験等に止まり、その侵襲性の
高さから人体への適用は不可能と考えられている。つま
るところ、現状では、セグメントショートニング法と同
等の精度で、心機能に関する指標を、人体に対して非侵
襲的に高精度で計測する方法はないといえる。

【０００５】また、心機能評価の指標を計算する場合、
心臓の伸縮方向を見極めるのが重要である。この伸縮方
向は、心臓の部分部分によって異なっているので、心臓
全体にわたってこれらを計測するためには、心臓の各部
分について全て伸縮方向を個別に手動で設定しなければ
ならず、時間や手間の観点から実用的とはいえない。

【０００６】

【発明が解決しようとする課題】本発明の目的は、局所
的な心筋等の収縮や拡張機能を、客観的、定量的、非侵
襲的に、しかも簡易に高精度に提示できる超音波診断装
置及び超音波画像処理装置を提供することにある。

【０００７】

【課題を解決するための手段】本発明は、被検体を超音
波で走査することにより得られた超音波画像上で関心部
位の領域を抽出し、この領域の輪郭に応じた方向に従っ
て組織の運動速度を補正し、この補正された運動速度に
基づいて関心部位の収縮拡張機能を表す指標を演算する
ものである。

【０００８】

【発明の実施の形態】以下、本発明を好ましい実施形態
により図面を参照しながら説明する。図１は、本実施形
態による超音波診断装置の構成を示している。クロック
パルス発生器５から出力されるクロックパルスは、レー
トパルス発生器４で分周されてレートパルスとなり、送
信遅延回路３を通して送信パルス発生器２にトリガとし
て送られる。送信パルス発生器２は、アレイ振動子を有
するプローブ１に、ある中心周波数ｆ0 を持つパルス長
の短いパルス電圧をレートパルスに同期して繰り返し印
加する。なお、送受信の周期は、レート周波数の逆数で
与えられる。送信遅延回路３は、複数の振動素子に印加
される送信パルス電圧のタイミングを制御し、超音波ビ
ームの方向を変えたりビームを集束させたりする。これ
を、以下、ビームフォーミングという。

【０００９】プローブ１から体内に放射され、体内組織
の音響インピーダンスの境界で反射された超音波パルス
は同じプローブ１で受信される。この受信信号は、プリ
アンプ６、受信遅延回路７を通りミキサ９に入力され
る。受信遅延回路７は、送信遅延回路３と同様にビーム
フォーミングを行うためのものである。ミキサ９に入力
した受信信号はクロックパルス発生器５から分周され参
照信号発生器８から供給される参照信号と掛け合わされ
て、さらにローパスフィルタ１２を通って検波され、ア
ナログディジタル変換器（Ａ／Ｄ）１６によりディジタ
ル信号に変換される。アナログディジタル変換器１６の
出力信号は、イメージング回路１３に供給されて振幅情
報となり、さらにディジタルスキャンコンバータ（ＤＳ
Ｃ）１４でＢモード画像とＭモード画像の少なくとも一
方に生成されてＣＲＴディスプレイ１５上に表示され
る。

【００１０】また、アナログディジタル変換器１６の出
力信号は、速度演算部１７にも供給され、ここで組織ド
プライメージング法にしたがって、超音波ビームの走査
範囲内の複数のサンプル点（以下単に、“点”と称す
る）それぞれにおける心筋等の組織運動に関する速度が
演算され、この速度情報がディジタルスキャンコンバー
タ１４で２次元速度分布（組織ドプラ画像）に生成さ
れ、ＣＲＴデイスプレイ１５上にカラーで表示される。

【００１１】この画像上に、オペレータによる入力部１
１の操作によって、「速度補正方向」または「関心領
域」が設定される。この関心領域において、そして組織
ドプライメージングにより得られた速度に基づいて、指
標演算部１８により、心筋等の組織運動の機能の正常／
異常を判定する材料になり得る指標が演算され、その演
算された指標がＤＳＣ１４を経由してＣＲＴディスプレ
イ１５に表示される。

【００１２】ところで、入力部１１、指標演算部１８及
びＣＰＵ１０の一部の機能を除く主に組織ドプラ画像を
生成するための構成部分は、既に実施され実用化されて
いるものであり、その詳細は省略し、以下に本発明の特
徴部分について詳しく説明する。

【００１３】図２（ａ）は、心尖アプローチにより得ら
れた心臓の左室の長軸断層像（組織の形態画像）の模式
図である。この長軸断層像を超音波診断装置により撮影
すると同時に、組織ドプライメージング（ＴＤＩ）によ
り、この断層上の心筋組織の２次元速度分布（組織ドプ
ラ画像）が得られているものとする。

【００１４】ここで、求める指標としては、心臓の伸縮
運動機能の診断に有効であることが要求され、また、動
物実験等で有効性が確認されているセグメントショート
ニングと臨床診断上同等の意義をなしていることが望ま
れる。つまり、長軸断層像上で、心筋がある方向に収縮
または伸張していれば、求める指標としては、ゼロでな
い値となり、その指標値の正負は収縮又は伸張を表し、
指標値の大きさは収縮や伸張の度合い（程度）を反映し
ていることが好ましい。このような要求を満たす指標と
しては、例えば、速度差、速度勾配、速度微分、速度差
分などが考えられる。

【００１５】この指標を計算する上で重要なのは、計算
上扱う伸縮方向をどのように推定するかということにあ
る。なぜなら、この方向が、心筋の真の伸縮方向に近け
れば近いほど指標の精度としては向上するからである。
本実施形態では、この伸縮方向を、心臓左室の輪郭に基
づいて推定するものである。

【００１６】まず、心臓左室の長軸断層像上で、操作者
がトラックボール等の入力器具１１を使って心筋の輪郭
を手動でトレースする。この心筋の輪郭トレースは、内
膜と外膜をなぞるようにしてもよいし、またもっと簡易
な方法として、操作者が内膜と外膜上に離散的に数点だ
け指定すると、それら点を直線で結ぶ又は曲線近似で連
結するようにしてもよい。また、完全に自動化法とし
て、ＡＣＴ法などの手法を用いて輝度勾配に基づいて心
腔と心筋の境界を抽出するようにしても良い。このトレ
ースにより限定された領域は、関心領域１９として扱わ
れる（図２（ｂ））。

【００１７】次に、指標演算部１８では、関心領域に基
づいて、左室の内膜２０（又は外膜）に沿う方向、例え
ば接線方向を、内膜２０上の各点に設定する。そして、
指標演算部１８では、設定した方向に関する速度、つま
り各点における心筋の伸縮に関する真又はそれに近い速
度を求めるために、関心領域中の各点の速度データを、
次のように角度補正する。周知の通り、ドプラ法では、
ＴＤＩに限らず、計測される速度データは、真の速度の
超音波ビーム方向成分として与えられる。したがって、
図２（ｃ）に示すように、内膜線２０上の点２１の速度
データＶを、点２１での局所的な内膜曲線に接するベク
トル２２と、超音波ビームの向きを表すベクトル２３と
のなす角（θ）２４を用いて、“Ｖ／ｃｏｓθ”によ
り、角度補正する。

【００１８】なお、図２（ｄ）に示すように、内膜線２
０上にない関心領域（心筋部分）の内部の点２５に関す
る角度補正は、内膜線２０に向かって引いた垂線と、内
膜線２０との交点２６において設定した方向を適用し
て、その方向のベクトル２７と点２５での超音波ビーム
の向きを表すベクトル２８とのなす角２９を用いて、内
部点２５の速度を上記と同様に角度補正する。

【００１９】次に、角度補正と指標の演算方法について
説明する。図２（ｅ）に示すように、関心領域中のある
２点をそれぞれＡ、Ｂとする。それぞれの点の角度補正
前の速度データをＶＡ、ＶＢとし、Ａを原点としたＢの
位置ベクトルをベクトルＡＢとする。さらに、点Ａ、点
Ｂから超音波プローブヘ向かうベクトルとベクトルＡＢ
とのなす角をθＡ、θＢとする。指標演算部１８は、指
標として、例えば速度勾配を以下の式により、演算す
る。

【００２０】速度勾配＝（ＶＢ／ｃｏｓθＢ−ＶＡ／ｃ
ｏｓθＡ）／｜ＡＢ｜ 点Ａ，Ｂ間の速度差は、この式の分子の部分として与え
られる。周知の通り、速度ＶＡ、ＶＢは、その各点がプ
ローブヘ接近する方向に動いていれば、正で与えられ、
またプローブから離れる方向に動いていれば、負で与え
られる。したがって、速度勾配としては、２点ＡＢ間で
心筋が伸張していれば、正で与えられ、収縮していれば
逆に負で与えられ、伸縮率は値として与えられる。

【００２１】この演算を隣接した２点間で行えば、速度
差分として捉えることができ、また、隣接２点問の距離
が微少であれば、速度微分として捉えることができる。
また、速度勾配は２点問の全ての角度補正された速度デ
ータのプロファイルから最小自乗法で直線フィッティン
グして、その傾きとして求めても良い。このように求め
た指標の符号によって、組織が収縮しているか伸張して
いるかが分かり、更に、その指標値の大きさによって収
縮や伸張の度合いが分かる。

【００２２】なお、指標は、関心領域中の各点で求めて
も良いが、図３に示すように、関心領域を複数の微小領
域に分割して、分割した領域ごとに指標を１つずつ求め
ても良い。微小領域は互いに重なっても良い。関心領域
の分割の仕方は、図３（ａ）に示すように、内膜に対し
て複数の垂線を等間隔に引いてその垂線により関心領域
を分割してもよいし、図３（ｂ）に示すように、長軸に
対して複数の垂線を等間隔に引いてその垂線により関心
領域を分割してもよい。また、図３（ｃ）に示すよう
に、ＡＳＥ（Ａｍｅｒｉｃａｎ Ｓｏｃｉｅｔｙ ｏｆ
Ｅｃｈｏｃａｒｄｉｏｇｒａｐｈｙ）が推奨するよう
な領域（セグメント）に分割して、各領域毎に計算して
もよい。

【００２３】このように関心領域を分割する場合には、
指標としては、各分割領域内の複数点の速度の平均値等
から上述した式にしたがって求めることが考えられ、こ
の場合、演算結果のばらつきなどのノイズを除去できる
という効果が期待できるものである。

【００２４】なお、指標演算の方向としては、上述した
ように心筋内膜（又は外膜）の接線方向の他に、全ての
点で、心臓の長軸方向と平行に統一してもよい。いずれ
も、心尖アプローチまたは経食道アプローチでのＴＤＩ
による速度推定では、ドプラ角が小さく精度の高い計測
が期待できる。

【００２５】垂線３０上の角度補正された複数の速度デ
ータを平均値などをとって代表させたものを、内膜線上
の点２６の速度データとしてもよいし、また、関心領域
は心筋に限らず、内膜または外膜でもよくその場合の速
度データは、内膜線上または外膜線上かその近傍のデー
タを用いれば良い。

【００２６】演算処理部１８によって演算された指標の
表示法としては、図４（ａ）に示すように、指標の２次
元分布（指標マップ）を、その値に応じた色相でカラー
表示すると、心機能の空間的な変化を把握しやすい。カ
ラー表示の仕方は、例えば、心筋が内膜等に接する方向
等に収縮する場合には赤色系で、伸張する場合には青色
系でする。また、指標の絶対値に応じて色づけを行うな
ど様々な方法が考えられる。さらに、指標の値そのもの
を、直接、数値３８で表示することも重要である。カラ
ー表示や数値表示は、画像上の関心領域上に重ねて表示
しても良いが、別の部分に表示してももちろん良い。

【００２７】理想的には、断層画像上全ての部分につい
て、上記指標を充分な精度で計算するのが好ましいが、
指標の演算方向（速度の補正方向）と超音波ビーム方向
とのなす角が、垂直に近い所定角度以上になる部分で
は、角度補正の誤差が大きくなるので、そのような部分
を指標演算部１８で検知し、計測不能領域３９として指
標カラーとは別系統のカラーで表示するか、または指標
カラーを付けないことが望ましい。こうすることによ
り、大きな誤差を含む可能性のある部分から誤診する危
険性を回避できる。なお、図３に示したように、関心領
域を分割した場合に、ある分割領域に計測不能領域が存
在した場合には、その部分を除いた範囲の点の速度デー
タからその分割領域の指標を計算すれば良い。

【００２８】これまでは、心機能の２次元的な評価につ
いて説明してきた。心尖アプローチや経食道アプローチ
では、図５に示すように、プローブ１を例えば心臓左室
の長軸のまわりに回転させながら、複数の断層像４１を
収集すれば、心臓の全体にわたって形態情報や速度情報
を収集することが可能である。あるいは、断層面を扇型
にスキャンさせるか、平行移動させながら複数の断層像
を収集しても良い。収集した断層像上でそれぞれ指標を
計算し、各断層像のあいだを補間することにより、心臓
の全体にわたって局所的な心機能情報（指標）が得られ
ることになる。

【００２９】もちろん、速度データの補間後に指標を計
算しても良い。このようにして得られた心臓左室全体の
情報を、図６に示すように、通常の３次元画像として表
示し、さらに、任意の断面を切り出せるようにしておけ
ば、例えば梗塞部位がどの程度の拡がりを持つのかが３
次元的に容易に把握できる。

【００３０】また、核医学で用いられているようなプル
ズアイマップ、つまり図７（ａ）に示すように、心臓の
内部の視点から見たような展開図を作り表示すれば、心
臓の３次元的な情報を一度に表示することが可能であ
る。

【００３１】また、図７（ｂ）に示すように、分割領域
ごとの指標をブルズアイマップ表示すれば、各冠動脈の
支配領域ごとに心機能の判別がつく。心機能の低下して
いる境界部分等では、分割領域の設定の仕方により、分
割領域個々の指標値が影響されることが予想される。そ
こで、領域分割しないで得た指標マップと、領域分割し
た得た指標マップとを一括して表示するか、図７（ｃ）
のように、並列して同時に表示し、両方に領域分割線４
４をガイド表示させて、この分割線４４をリンクさせ、
一方を動かせば、それに応じて他方も動くようにすれ
ば、最適な領域分割を行うことができる。この時、分割
線４４の移動に応じて、各分割領域の指標をリアルタイ
ムに再計算して表示するのが望ましい。

【００３２】上述したように、心筋輪郭に基づいてその
心臓内膜や外膜に接する方向に関して指標を計算すると
説明したが、これは、関心領域を設定した後では、内膜
輪郭に接する接線方向は完全に自動化することができ、
操作者の手間を大幅に少なくすることができる。

【００３３】一方、心臓の短軸断層像の場合には、心筋
の各点から短軸断層像の重心へ向かう方向で計測速度を
角度補正して、壁厚変化に関連した指標を計算すること
が考えられる。以下、この方法について説明する。

【００３４】図８（ａ）は、心臓左室の短軸断層像の模
式図である。これが、ＣＲＴディスプレイ１５に表示さ
れているものとする。まず、操作者による入力器具１１
の操作により、左室短軸断層像の心筋の内膜４５がトレ
ースされる（図８（ｂ））。このトレースの方法として
は、内膜をなぞるようにしてもよいし、またもっと簡易
な方法として、操作者が内膜上に離散的に数点だけ指定
すると、それら点を直線で結ぶ又は曲線近似で連結する
ようにしてもよい。また、完全に自動化法として、ＡＣ
Ｔ法などの手法を用いて輝度勾配に基づいて心腔と心筋
の境界を抽出するようにしても良い。外膜４６も、内膜
と同様にトレースする。あるいは、内膜上の各点から一
定距離外側に離れた曲線を外膜としても良いし、内膜形
状の相似拡大としても良い。

【００３５】こうして内膜や外膜がトレースされると、
指標演算部１８では、内膜または外膜の輪郭座標から重
心位置を計算するか、あるいは内膜と外膜で囲まれる心
筋部分の面積重心などから重心４７を特定する（図８
（ｃ））。関心領域である心筋の各点で、重心に向かう
方向４８に指標演算方向を設定し、演算する（図８
（ｄ））。

【００３６】さらに、ＡＳＥで推奨するセグメントに関
心領域を分割するさいにも、重心点を中心として放射状
に分割線４９を引いて、関心領域を複数に分割分割線４
９を自動的に描画することが可能となり、操作者の手間
を省くことができる（図８（ｅ））。また、心筋の厚さ
方向の指標を計算する場合には、図８（ｆ）のように、
トレースした内膜の各部分で内膜曲線に垂直な方向５０
を求めれば良い。この方向設定後、速度差、速度勾配、
速度微分、速度差分を計算すればよい。この場合、これ
らの指標は、壁厚変化速度に関連したものとなる。

【００３７】以上のように、組織の収縮・伸張を表す指
標を演算する際に関心領域の輪郭情報を用いて、演算方
向の設定、関心領域の分割等を行えば、操作者による操
作は関心領域を設定するのみとなり、非常に簡便かつよ
り正確に指標を計測することが可能となる。また、指標
演算方向を関心領域中の任意の点で必要なときに表示
し、調節できるようにすれば、操作者は指標演算部１８
によって設定された方向を確認・変更することが可能に
なる。変更後、指標演算部１８は再び、変更された方向
に基づいて角度補正、指標の演算をする。

【００３８】なお、関心領域の設定、指標の演算や表示
などは、図９に示すように、図１の計測処理部の機能を
装備したＰＣやワークステーシヨンなどコンピュータ５
１により、超音波診断装置とは切り離して行っても良
い。ＤＳＣ１４から出力される断層画像、２次元速度分
布などのデータをコンピュータ５１に送り、演算した結
果をコンピュータ付属のモニタに表示しても良いし、診
断装置のＣＲＴ１５に表示しても良い。本発明は、上述
した実施形態に限定されることなく種々変形して実施可
能である。

【００３９】

【発明の効果】本発明によると、例えば心臓の運動機能
に関する指標を心臓左室全体にわたって簡便に、客観的
且つ定量的に評価可能であり、診断に有用な情報を提供
することができる。

【図面の簡単な説明】

【図１】本発明の好ましい実施形態に係る超音波診断装
置のブロック図。

【図２】（ａ）は心尖アプローチにより得た心臓の長軸
断層像を示す図、（ｂ）は手動又は自動の心筋トレース
により設定された関心領域を示す図、（ｃ）は心筋内膜
上の点に関する角度補正の説明図、（ｄ）は心筋内部の
点に関する角度補正の説明図、（ｅ）は２点間の速度勾
配等の指標の計算方法の説明図。

【図３】（ａ）は第１方法で分割された関心領域を示す
図、（ｂ）は第２方法で分割された関心領域を示す図、
（ｃ）は第３方法で分割された関心領域を示す図。

【図４】（ａ）は点毎に計算した指標マップの表示例、
（ｂ）は分割領域毎に計算した指標マップの表示例を示
す図。

【図５】心尖アプローチによる３次元走査の例を示す
図。

【図６】指標の３次元画像の表示例を示す図。

【図７】（ａ）は点毎に計算した指標の３次元画像のブ
ルズアイマップ表示例を示す図、（ｂ）は分割領域毎に
計算した指標の３次元画像のブルズアイマップ表示例を
示す図、（ｃ）は点毎と分割領域毎に計算した指標のブ
ルズアイマップ画像の同時表示例を示す図。

【図８】（ａ）は心臓左室の短軸断層像を示す図、
（ｂ）は心臓左室の内膜のトレースを示す図、（ｃ）は
心臓左室の外膜のトレースと短軸断面の重心を示す図、
（ｄ）は重心を中心にした伸縮指標の計算方向を示す
図、（ｅ）は重心から放射状に関心領域を分割する方法
を示す図、（ｆ）は心筋厚の方向を示す図。

【図９】本発明の変形例に関するプロック図。

【符号の説明】

１…プローブ、 ２…パルス発生器、 ３…送信遅延回路、 ４…レートパルス発生器、 ５…クロック発生器、 ６…プリアンプ、 ７…受信遅延回路、 ８…参照信号発生器、 ９…ミキサ、 １０…ＣＰＵ、 １１…入力器具、 １２…低域通過型フィルタ、 １３…イメージング回路、 １４…ディジタルスキャンコンバータ、 １５…ＣＲＴディスプレイ、 １６…アナログディジタルコンバータ、 １７…速度演算部、 １８…指標演算部、 ５１…コンピュータ。

Claims (13)

【特許請求の範囲】
1. 【請求項１】 被検体を超音波で走査して受信信号を得
   る送受信手段と、 前記受信信号に基づいて超音波画像を生成する手段と、 前記受信信号に基づいて前記超音波画像の複数点に関し
   て組織の運動速度を求める手段と、 前記超音波画像上で関心部位の領域を抽出する手段と、 前記抽出された関心部位の輪郭に応じた方向に従って前
   記運動速度を補正する手段と、 前記補正された運動速度に基づいて前記関心部位の収縮
   拡張機能を表す指標を演算する手段とを具備することを
   特徴とする超音波診断装置。
2. 【請求項２】 被検体の超音波画像データと、前記超音
   波画像の複数点に関する組織の運動速度データとを入力
   する手段と、 前記超音波画像上で関心部位の領域を抽出する手段と、 前記抽出された関心部位の輪郭に応じた方向に従って前
   記運動速度を補正する補正手段と、 前記補正された運動速度に基づいて前記関心部位の収縮
   拡張機能を表す指標を演算する指標演算手段とを具備す
   ることを特徴とする超音波画像処理装置。
3. 【請求項３】 前記補正手段は、前記関心部位の領域の
   輪郭に沿うように前記方向を設定することを特徴とする
   請求項１または請求項２記載記載の装置。
4. 【請求項４】 前記補正手段は、前記指標の絶対値が最
   大となるように前記方向を設定することを特徴とする請
   求項１または請求項２記載の装置。
5. 【請求項５】 前記補正手段は、前記関心部位の領域の
   輪郭に垂直に前記方向を設定することを特徴とする請求
   項１または請求項２記載の装置。
6. 【請求項６】 前記補正手段は、前記関心部位の領域の
   重心に向かうように前記方向を設定することを特徴とす
   る請求項１または請求項２記載の装置。
7. 【請求項７】 前記補正手段は、心臓の断層画像上で心
   臓内膜または心臓外膜に沿うように前記方向を設定する
   ことを特徴とする請求項１または請求項２記載の装置。
8. 【請求項８】 前記方向と超音波ビームのなす角度を表
   示する手段をさらに備えたことを特徴とする請求項１ま
   たは請求項２記載の装置。
9. 【請求項９】 前記方向と超音波ビームのなす角度を前
   記関心部位に重ねてカラーで表示する手段をさらに備え
   たことを特徴とする請求項１または請求項２記載の装
   置。
10. 【請求項１０】 前記指標演算手段で複数位置にわたっ
    て演算された指標を、２次元または３次元で表示する指
    標表示手段をさらに備えることを特徴とする請求項１ま
    たは請求項２記載の装置。
11. 【請求項１１】 前記指標表示手段は、前記指標を指標
    値に応じたカラーで表示することを特徴とする請求項１
    ０記載の装置。
12. 【請求項１２】 前記指標演算手段は、前記関心部位の
    領域を複数の分割領域に分割する手段を有し、前記分割
    領域ごとに前記指標を演算することを特徴とする請求項
    １または請求項２記載の装置。
13. 【請求項１３】 前記補正手段は、前記分割領域ごとに
    前記方向を１つずつ設定することを特徴とする請求項１
    ２記載の装置。