(57)【特許請求の範囲】
【請求項１】
磁気共鳴イメージング装置が備える情報取得部、補間処理部、線形計算部、算出結果画像情報格納部、画像出力部の各部を当該磁気共鳴イメージング装置の制御手段が制御して、生体である測定対象の骨格画像を得る、磁気共鳴イメージング装置の画像処理方法であって、
前記情報取得部を制御して、測定対象に縦方向、横方向、高さ方向にそれぞれ所定の間隔で複数設けられた測定点それぞれについて磁気縦緩和測定及び磁気横緩和測定の２種類の測定方法により測定された結果得られた各磁気共鳴スペクトル強度値に基づき、各測定点における前記磁気共鳴スペクトル強度値の集合である磁気共鳴イメージング画像情報を、磁気縦緩和測定及び磁気横緩和測定についてそれぞれ取得し、
前記補間処理部を制御して、磁気縦緩和測定により得られた磁気共鳴イメージング画像と磁気横緩和測定により得られた磁気共鳴イメージング画像とを３次元的に位置合わせし、これら磁気縦緩和測定と磁気横緩和測定の一方の磁気共鳴イメージング画像から他方の磁気共鳴イメージング画像と同一測定点における磁気共鳴スペクトル強度値を補間することにより求め、
前記線形計算部を制御して、各測定点における磁気共鳴スペクトル強度値を、磁気縦緩和測定及び磁気横緩和測定の磁気共鳴スペクトル強度値に基づく線形計算により求め、
前記算出結果画像情報格納部を制御して、各測定点について前記線形計算の結果を算出結果画像情報として格納し、
前記画像出力部を制御して、前記算出結果画像情報に基づいて新たな画像情報として骨格画像を導出しこれを出力することを特徴とする磁気共鳴イメージング装置の画像情報処理方法。
【請求項２】
前記補間処理部の制御においては、磁気縦緩和測定により得られた磁気共鳴イメージング画像と磁気横緩和測定により得られた磁気共鳴イメージング画像とを３次元的に位置合わせし、磁気横緩和測定の磁気共鳴イメージング画像からこれら磁気縦緩和測定の磁気共鳴イメージング画像と同一測定点における磁気共鳴スペクトル強度値を補間することにより求める請求項１記載の磁気共鳴イメージング装置の画像処理方法。
【請求項３】
前記磁気共鳴スペクトル強度値が水素原子核の磁気共鳴スペクトル強度値であることを特徴とする請求項２記載の磁気共鳴イメージング装置の画像情報処理方法。
【請求項４】
前記線形計算部の制御においては、水分子中の水素原子核に由来する磁気共鳴スペクトルを消去する線形演算により各測定点における磁気共鳴スペクトル強度値を求める請求項３記載の磁気共鳴イメージング装置の画像処理方法。
【請求項５】
前記磁気共鳴スペクトル強度値は、各測定点における磁気共鳴イメージング画像をビット数で表したものであり、
前記線形計算部の制御においては、各測定点における磁気共鳴スペクトル強度値を、当該測定点における磁気縦緩和測定の磁気共鳴スペクトル強度値をビット反転した値から、前記補間により求められた磁気横緩和測定の磁気共鳴スペクトル強度値の定数倍値を減じることにより求め、当該定数として、水の磁気共鳴スペクトル強度値を０とする値を用いている請求項４記載の磁気共鳴イメージング装置の画像処理方法。
【請求項６】
生体である測定対象の骨格画像を得る磁気共鳴イメージング装置であって、
測定対象に縦方向、横方向、高さ方向にそれぞれ所定の間隔で複数設けられた測定点それぞれについて磁気縦緩和測定及び磁気横緩和測定の２種類の測定方法により測定された結果得られた各磁気共鳴スペクトル強度値に基づき、各測定点における前記磁気共鳴スペクトル強度値の集合である磁気共鳴イメージング画像情報を、磁気縦緩和測定及び磁気横緩和測定についてそれぞれ取得する情報取得部、
磁気縦緩和測定により得られた磁気共鳴イメージング画像と磁気横緩和測定により得られた磁気共鳴イメージング画像とを３次元的に位置合わせし、これら磁気縦緩和測定と磁気横緩和測定の一方の磁気共鳴イメージング画像から他方の磁気共鳴イメージング画像と同一測定点における磁気共鳴スペクトル強度値を補間する補間処理部、
各測定点における磁気共鳴スペクトル強度値を、磁気縦緩和測定及び磁気横緩和測定の磁気共鳴スペクトル強度値に基づく線形演算により求める線形計算部、
各測定点について前記線形計算部による計算結果を算出結果画像情報として格納する算出結果画像情報格納部、
前記算出結果画像情報格納部に格納した算出結果画像情報に基づいて新たな画像情報として骨格画像を導出しこれを出力する画像出力部、
を備えることを特徴とする磁気共鳴イメージング装置。
【請求項７】
前記補間処理部は、磁気縦緩和測定により得られた磁気共鳴イメージング画像と磁気横緩和測定により得られた磁気共鳴イメージング画像とを３次元的に位置合わせし、磁気横緩和測定の磁気共鳴イメージング画像からこれら磁気縦緩和測定の磁気共鳴イメージング画像と同一測定点における磁気共鳴スペクトル強度値を補間することにより求める処理を行う請求項６記載の磁気共鳴イメージング装置。
【請求項８】
前記磁気共鳴スペクトル強度値は、水素原子核の磁気共鳴スペクトル強度値である請求項７記載の磁気共鳴イメージング装置。
【請求項９】
前記線形計算部は、水分子中の水素原子核に由来する磁気共鳴スペクトルを消去する線形演算により各測定点における磁気共鳴スペクトル強度値を求める処理を行う請求項８記載の磁気共鳴イメージング装置。
【請求項１０】
前記磁気共鳴スペクトル強度値は、各測定点における磁気共鳴イメージング画像をビット数で表したものであり、
前記線形計算部は、各測定点における磁気共鳴スペクトル強度値を、当該測定点における磁気縦緩和測定の磁気共鳴スペクトル強度値をビット反転した値から、前記補間により求められた磁気横緩和測定の磁気共鳴スペクトル強度値の定数倍値を減じることにより求める処理を行い、当該定数として、水の磁気共鳴スペクトル強度値を０とする値を用いている請求項９記載の磁気共鳴イメージング装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、人体等の３次元物体内部を非破壊検査する際に用いられる磁気共鳴イメージング装置の画像情報処理方法並びに同装置に関する。
【０００２】
【従来の技術】
従来、人体の内部を検査する際には、Ｘ線撮影が広く用いられてきた。人体にＸ線を照射すると、Ｘ線は人体を透過するが、透過率は骨の部分では低いので、骨折箇所をＸ線写真により診断することができる。
【０００３】
しかし、Ｘ線を多量に人体に照射すると細胞が損傷するので、Ｘ線の照射は頻繁に行うことができないという問題がある。
【０００４】
この不具合を解消すべく、磁気共鳴現象を利用し、人体に磁場を照射して磁気共鳴スペクトル強度分布により人体の内部を検査する試みがなされている。（例えば、非特許文献１参照。）
【０００５】
【非特許文献１】
松尾政彦、濱田博文、藤川尚宏、二宮英彰、江田英雄、宮内哲「光による脳機能計測のための構造画像処理の研究」、第４１回日本エム・イー学会大会発表予稿集、ｐ．５５
【０００６】
【発明が解決しようとする課題】
ところが、骨格等、いくつかの部位については通常の磁気共鳴イメージング画像では像が現れず、このような部位を非破壊で検査する新たな方法が求められている。
【０００７】
これ以外でも、３次元物体内部を非破壊検査する際に、通常の磁気共鳴イメージング画像では求める画像情報が得られない場合があり、このような場合に対応する新たな方法が求められている。
【０００８】
本発明は以上に述べた課題を解決し、３次元物体内部を放射線等の有害電磁波を用いずに非破壊検査する新たな方法を提供するものである。
【０００９】
【課題を解決するための手段】
すなわち本発明に係る磁気共鳴イメージング装置の画像情報処理方法は、磁気共鳴イメージング装置が備える情報取得部、補間処理部、線形計算部、算出結果画像情報格納部、画像出力部の各部を当該磁気共鳴イメージング装置の制御手段が制御して、生体である測定対象の骨格画像を得る、磁気共鳴イメージング装置の画像処理方法であって、前記情報取得部を制御して、測定対象に縦方向、横方向、高さ方向にそれぞれ所定の間隔で複数設けられた測定点それぞれについて磁気縦緩和測定及び磁気横緩和測定の２種類の測定方法により測定された結果得られた各磁気共鳴スペクトル強度値に基づき、各測定点における前記磁気共鳴スペクトル強度値の集合である磁気共鳴イメージング画像情報を、磁気縦緩和測定及び磁気横緩和測定についてそれぞれ取得し、前記補間処理部を制御して、磁気縦緩和測定により得られた磁気共鳴イメージング画像と磁気横緩和測定により得られた磁気共鳴イメージング画像とを３次元的に位置合わせし、これら磁気縦緩和測定と磁気横緩和測定の一方の磁気共鳴イメージング画像から他方の磁気共鳴イメージング画像と同一測定点における磁気共鳴スペクトル強度値を補間することにより求め、前記線形計算部を制御して、各測定点における磁気共鳴スペクトル強度値を、磁気縦緩和測定及び磁気横緩和測定の磁気共鳴スペクトル強度値に基づく線形計算により求め、前記算出結果画像情報格納部を制御して、各測定点について前記線形計算の結果を算出結果画像情報として格納し、前記画像出力部を制御して、前記算出結果画像情報に基づいて新たな画像情報として骨格画像を導出しこれを出力することを特徴とする。
【００１０】
このような方法であれば、磁気縦緩和測定及び磁気横緩和測定の２種類の磁気共鳴イメージング測定から得られる画像情報に基づいた線形演算により、これまでは磁気共鳴イメージング画像では像が現れなかった生体の骨格の画像を得ることができ、骨格の画像を得る際に有害なＸ線等の放射線を用いる必要がなくなり、このような検査の安全性の向上が図られる。
【００１１】
特に、２種類の磁気共鳴スペクトル強度測定法として、磁気縦緩和測定及び磁気横緩和測定を用いているのは、水分子中の水素原子核の磁気縦緩和速度は低く、磁気横緩和速度は高いことに着目したためである。
【００１２】
すなわち、磁気縦緩和測定による磁気共鳴イメージング画像情報と、磁気横緩和測定による磁気共鳴イメージング画像情報とを同一の測定対象、測定点について取得し、水分子中の水素原子核に由来するスペクトルを消去することで、骨格の画像をより容易に導出することが可能となる。
【００１３】
補間処理部の制御においては、磁気縦緩和測定により得られた磁気共鳴イメージング画像と磁気横緩和測定により得られた磁気共鳴イメージング画像とを３次元的に位置合わせし、磁気横緩和測定の磁気共鳴イメージング画像からこれら磁気縦緩和測定の磁気共鳴イメージング画像と同一測定点における磁気共鳴スペクトル強度値を補間することにより求めるとよい。
【００１４】
ここで、磁気共鳴スペクトル強度値は、各測定点における磁気共鳴イメージング画像をビット数で表したものとすることができる。
【００１５】
この場合、線形計算部の制御においては、各測定点における磁気共鳴スペクトル強度値を、当該測定点における磁気縦緩和測定の磁気共鳴スペクトル強度値をビット反転した値から、前記補間により求められた磁気横緩和測定の磁気共鳴スペクトル強度値の定数倍値を減じることにより求め、当該定数として、水の磁気共鳴スペクトル強度値を０とする値を用いることで、具体的に骨格の画像情報が得られるようになる。
【００１６】
なお、本発明の方法を応用することで、骨格と同様にこれまで磁気共鳴イメージング画像としては得られなかった部分の画像を、放射線等の有害電磁波を用いずに得ることができるようになる。
【００１７】
このような画像情報処理方法を利用する本発明の磁気共鳴イメージング装置は、測定対象に縦方向、横方向、高さ方向にそれぞれ所定の間隔で複数設けられた測定点それぞれについて磁気縦緩和測定及び磁気横緩和測定の２種類の測定方法により測定された結果得られた各磁気共鳴スペクトル強度値に基づき、各測定点における前記磁気共鳴スペクトル強度値の集合である磁気共鳴イメージング画像情報を、磁気縦緩和測定及び磁気横緩和測定についてそれぞれ取得する情報取得部、磁気縦緩和測定により得られた磁気共鳴イメージング画像と磁気横緩和測定により得られた磁気共鳴イメージング画像とを３次元的に位置合わせし、これら磁気縦緩和測定と磁気横緩和測定の一方の磁気共鳴イメージング画像から他方の磁気共鳴イメージング画像と同一測定点における磁気共鳴スペクトル強度値を補間する補間処理部、各測定点における磁気共鳴スペクトル強度値を、磁気縦緩和測定及び磁気横緩和測定の磁気共鳴スペクトル強度値に基づく線形演算により求める線形計算部、各測定点について前記線形計算部による計算結果を算出結果画像情報として格納する算出結果画像情報格納部、前記算出結果画像情報格納部に格納した算出結果画像情報に基づいて新たな画像情報として骨格画像を導出しこれを出力する画像出力部、を備えることを特徴とするものである。
【００１８】
このような装置を用いることで、放射線を利用することなく安全に生体の骨格を画像として得ることができるようになる。
【００１９】
【発明の実施の形態】
以下に本発明の一実施形態について図面を参照して述べる。
【００２０】
本実施形態に係る磁気共鳴イメージング画像情報処理方法は、図１に示すような磁気共鳴イメージング装置（以下ＭＲＩ装置）Ｍを用いるものであり、このＭＲＩ装置Ｍにより、磁気縦緩和測定（以下Ｔ1測定）による磁気共鳴イメージング画像情報と、磁気横緩和測定（以下Ｔ2測定）による磁気共鳴イメージング画像情報を取得して行う。
【００２１】
前記ＭＲＩ装置Ｍは、医療検査用として広く用いられている周知のものと同様の構成を有し、同様に作用してＴ1測定による磁気共鳴イメージング画像情報及びＴ2測定による磁気共鳴イメージング画像情報を取得する。ここで、本実施形態では人体の左右方向にｘ軸、前後方向にｙ軸、上下方向にｚ軸を設定していて、ｘｙ平面をスライス画像平面とし、このスライス画像平面上にマトリクスをさらに設定している。測定点は、測定対象たる被験者の人体の左右方向及び前後方向、すなわち前記ｘ軸方向及びｙ軸方向には同一ピッチで同一マトリクス点数を設定している。
【００２２】
このＭＲＩ装置Ｍは、機能ブロック図を図２に示すように、Ｔ1測定及びＴ2測定の結果得られた各磁気共鳴スペクトル強度値に基づき、各測定点における前記磁気共鳴スペクトル強度値の集合である磁気共鳴イメージング画像情報を、Ｔ１測定とＴ2測定についてそれぞれ取得する情報取得部１、Ｔ1測定により取得した磁気共鳴イメージング画像情報を格納する第１取得画像情報格納部２、Ｔ2測定により取得した磁気共鳴イメージング画像情報を格納する第２取得画像情報格納部３、及び画像を出力する画像出力部８として機能する。前記第１取得画像情報格納部２及び前記第２取得画像情報格納部３は、それぞれこのＭＲＩ装置Ｍの内部メモリのメモリ空間に構成している。そして、取得した磁気共鳴イメージング画像情報は、前記第１取得画像情報格納部２及び前記第２取得画像情報格納部３に、図３に示すように格納している。具体的には、各位置でのスペクトル強度値は、最高を「６５５３５」、最低を「０」とする１６ビット（６５５３６段階）の数値で示していて、ｚ枚目のスライス画像平面のマトリクス点（ｘ、ｙ）のスペクトル強度をｚ枚目のｘ列ｙ段に格納している。ここで、水のスペクトル強度は、Ｔ1測定では低く、Ｔ2測定では高い。骨のスペクトル強度は、Ｔ1測定、Ｔ2測定のいずれにおいても低い。脳のスペクトル強度は、Ｔ1測定、Ｔ2測定のいずれにおいても中程度だが、Ｔ1測定においてやや高い値を示す。皮膚のスペクトル強度は、Ｔ1測定では中程度で、Ｔ2測定ではやや高い。以上に述べた傾向を図４に示す。また、Ｔ1測定により得られる画像の一例を図５、Ｔ2測定により得られる画像の一例を図６にそれぞれ示す。
【００２３】
しかして、本実施形態では、これら２つの磁気共鳴イメージング画像情報に対して、同一位置のスペクトル強度同士を線形演算して、頭部の骨格を示す画像情報を得るようにしている。
【００２４】
具体的には、図２に示すように、このＭＲＩ装置Ｍを補間処理部４、補間処理結果情報格納部５、線形計算部６、及び算出結果画像情報格納部７としても機能するようにしている。前記補間処理結果情報格納部５及び前記算出結果画像情報格納部７は、いずれもメモリ空間に構成している。前記補間処理部４は、Ｔ1測定により得られた磁気共鳴イメージング画像情報とＴ2測定により得られた磁気共鳴イメージング画像情報とを３次元的に位置合わせし、Ｔ2測定により得られた磁気共鳴イメージング画像情報からＴ1測定に用いた測定点と同一位置のスペクトル強度を補間により求めるようにしている。前記補間処理結果情報格納部５は、前記補間処理部４による計算結果を格納している。前記線形計算部６は、前記第１取得画像情報格納部２に格納されたスペクトル強度値をビット反転するとともに、このビット反転の結果と前記補間処理結果情報格納部５に格納されたスペクトル強度値の定数ａ倍との差を計算するようにしている。ここで、定数ａは水のスペクトル強度の計算結果が０になるように設定している。前記算出結果画像情報格納部７は、算出した結果の集合である算出結果画像情報を格納している。
【００２５】
さらに、具体的な処理の流れを、図７及び図８を参照して以下に示す。まず、情報取得部１がＴ1測定を行い、結果を第１取得画像情報格納部２に格納する（Ｓ１）。次いで、同じく情報取得部１がＴ2測定を行い、結果を第１取得画像情報格納部２に格納する（Ｓ２）。それから、Ｔ1測定により得られた磁気共鳴イメージング画像情報とＴ2測定により得られた磁気共鳴イメージング画像情報とを３次元的に位置合わせし、このＴ2測定より得られた磁気共鳴イメージング画像情報から、Ｔ1測定に用いた測定点と同一位置のスペクトル強度を補間処理部４が補間により求め、補間処理結果情報格納部５に格納する（Ｓ３）。具体的には、上述したようにマトリクスの大きさは同一であり、スライス幅を大きくしているので、両側最寄りのスライスの同一マス目のマトリクスのスペクトル強度値から求める測定点のスペクトル強度値を線形補間により求める。ここまでの処理の流れを図７に示す。それから、図８に示すように、マトリクスのｘ方向、すなわち図３の列を表す変数ｘ、マトリクスのｙ方向、すなわち図２の段を表す変数ｙ、スライス画像の番号を表す変数ｚを予め１にリセットし（Ｓ４）、ｚ枚目のｘ列ｙ段の信号強度について、第１取得画像情報格納部２のｚ枚目のｘ列ｙ段のスペクトル強度値をビット反転したものから補間処理結果情報格納部５のｚ枚目のｘ列ｙ段のスペクトル強度値のａ倍を減じた値を算出結果画像情報格納部７に格納する（Ｓ５、線形計算ステップ）。すなわち、求めるスペクトル強度値をｔ、Ｔ1測定により得られたスペクトル強度値をｔ1、Ｔ2測定により得られたスペクトル強度値をｔ2とおくと、
ｔ＝（６５５３５−ｔ1）−ａｔ2
で表される線形演算を行い、算出結果画像情報格納部７のｚ枚目のｘ列ｙ段の位置に格納する。そして、ｘがマトリクスのマス目数に一致しているか否かの判定を行い（Ｓ６）、一致していない場合にはｘに１を加算して（Ｓ７）線形計算ステップに戻る。一方、一致している場合には、次いでｙがマトリクスのマス目数に一致しているか否かの判定を行い（Ｓ８）、一致していない場合にはｙに１を加算するとともにｘを１にリセットして（Ｓ９）線形計算ステップに戻る。一方、一致している場合には、さらにｚがスライス画像の枚数に一致しているか否かの判定を行い（Ｓ１０）、一致していない場合にはｚに１を加算するとともにｘ及びｙを１にリセットして（Ｓ１１）線形計算ステップに戻る。そして、ｚがスライス画像の枚数に一致していた場合には、以上のステップで算出結果画像情報格納部７に格納した画像情報に基き画像を出力する（Ｓ１２）。
【００２６】
このように画像情報を得るようにすることで、通常のＭＲＩ測定検査では得られない骨格の形状データを、人体に有害な放射線を用いることなく取得することができる。特に、人間の頭部にこの磁気共鳴イメージング画像情報処理方法を用いると、頭蓋骨の形状を得ることができ、頭蓋骨の骨折箇所の判定等をより安全に行うことができるようになる。
【００２７】
なお、本発明は上述した実施の形態に限られない。
【００２８】
例えば、水素原子核密度測定による磁気共鳴イメージング画像情報を用いるようにしてもよい。さらに、この水素原子核密度測定による磁気共鳴イメージング画像情報と、磁気縦緩和測定による磁気共鳴イメージング画像情報と、磁気横緩和測定による磁気共鳴イメージング画像情報との３種類の画像情報の間で所定位置のスペクトル強度値間の線形演算を行い、新たな画像情報を得るようにしてもよい。この水素原子核密度測定による各部のスペクトル強度は、図４に示すように、磁気縦緩和測定及び磁気横緩和測定による各部のスペクトル強度とさらに異なる傾向を示すので、線形演算の変数にこの水素原子核密度測定による所定位置のスペクトル強度値をも採用すると、さらに新たな情報を得ることができ得る。
【００２９】
また、人体の骨格以外の部分、さらに人体以外に関する情報を得るために以上に述べたような磁気共鳴イメージング画像情報処理方法を用いてもよい。その際、線形演算のアルゴリズムは、求める情報の種類により任意に設定してもよい。
【００３０】
さらに、１種類の磁気共鳴イメージング画像情報の測定点に対応させて所定位置を決定するのでなく、ＭＲＩ装置が取得した磁気共鳴イメージング画像情報の測定点とは無関係に任意に所定位置を決定し、取得した磁気共鳴イメージング画像情報全てについて所定位置のスペクトル強度値を補間により求めてもよい。加えて、補間のアルゴリズムは、上述した実施形態で用いた線形補間でなく、他のアルゴリズムを採用してもよい。
【００３１】
そして、水素原子核だけでなく、炭素原子核や窒素原子核等、他の原子核の核磁気共鳴スペクトルを用いてもよい。
【００３２】
加えて、磁気イメージング画像情報により得られた画像情報とＸ線ＣＴにより得られる画像情報を比較してもよい。具体的には、磁気イメージング画像情報により得られた画像情報とＸ線ＣＴにより得られる画像情報とを同一画面上に同時に出力し、またはこれらの画像情報をともに用紙等の印刷媒体に出力して、これらの画像情報を同時に視認可能な状態にすることや、これらの画像情報の間で前記所定位置のスペクトル強度値間の線形演算を行い、さらに新たな画像情報を導出すること等を行うとよい。このようにすれば、Ｘ線ＣＴによる骨の位置を直接示す情報との比較を行い、骨の状態をより正確に得ることができる。
【００３３】
その他、本発明の趣旨を損ねない範囲で種々に変形してよい。
【００３４】
【発明の効果】
本発明は、複数の磁気共鳴イメージング画像情報に基づく線形演算により、骨格を示す画像情報を新たな画像情報として導出し出力するようにしているので、測定対象内部の非破壊検査を行う際に人体に有害なＸ線等の放射線を用いる必要がなくなり、このような検査の安全性の向上を図ることができる。
【図面の簡単な説明】
【図１】本発明の一実施形態に係る磁気共鳴イメージング画像情報処理方法に用いるＭＲＩ装置を示す概略図。
【図２】同実施形態に係るＭＲＩ装置の機能ブロック図。
【図３】同実施形態に係るＭＲＩ装置が取得した磁気共鳴イメージング画像情報の格納の態様を示す概略図。
【図４】同実施形態に係るＭＲＩ装置が取得した画像の各部の信号強度の傾向を示す図。
【図５】同実施形態に係るＭＲＩ装置が磁気縦緩和測定により取得した画像の一例を示す図。
【図６】同実施形態に係るＭＲＩ装置が磁気横緩和測定により取得した画像の一例を示す図。
【図７】同実施形態に係る磁気共鳴イメージング画像情報処理方法における処理の流れを示すフローチャート。
【図８】同実施形態に係る磁気共鳴イメージング画像情報処理方法における処理の流れを示すフローチャート。
【符号の説明】
Ｍ…磁気共鳴イメージング装置
１…情報取得部
２…第１取得画像情報格納部
３…第２取得画像情報格納部
４…補間処理部
６…線形計算部
７…算出結果画像情報格納部
８…画像出力部
【図１】
0003757276.tif 000002
【図２】
0003757276.tif 000003
【図３】
0003757276.tif 000004
【図４】
0003757276.tif 000005
【図５】
0003757276.tif 000006
【図６】
0003757276.tif 000007
【図７】
0003757276.tif 000008
【図８】
0003757276.tif 000009