(57)【要約】  逆問題の解析を汎用的かつ高速に行い、物体内部の情報を簡便に映像化する。  散乱トモグラフィ方法は、波動を曲面上に配置された複数の送信アンテナ素子から物体に放射するステップと、散乱波を曲面上に配置された複数の受信アンテナ素子により受信するステップと、受信アンテナ素子により受信した散乱波を示す散乱波データから、物体の内部情報に関する画像を再構成するステップとを含み、画像を再構成するステップにおいて、物体の内部情報に関する画像を再構成するための関数φをあらかじめ設定し、関数の漸近式が満たす方程式を構築し、散乱波データから、方程式を解くことにより得られる映像化関数ρを導出し、映像化関数により、物体の内部情報に関する画像を再構成する。
2014125815.tif 000002
【特許請求の範囲】
【請求項１】
物体に放射した波動の散乱波を解析する散乱トモグラフィ方法であって、
前記波動を、曲面上に配置された複数の送信アンテナ素子から前記物体に放射するステップと、
前記散乱波を、曲面上に配置された複数の受信アンテナ素子により受信するステップと、
前記受信アンテナ素子により受信した散乱波を示す散乱波データを境界条件として用い、前記物体の内部情報に関する画像を再構成するステップとを含み、
前記画像を再構成するステップにおいて、
前記物体の内部情報に関する画像を再構成するための、（式Ａ）に定義される関数φをあらかじめ設定し、
前記関数φの漸近式が満たす、（式Ｂ）に定義される方程式を構築し、
測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、（式Ｃ）に定義される映像化関数ρを導出し、
前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する
散乱トモグラフィ方法。
【数６３】
2014125815.tif 000065
【数６４】
2014125815.tif 000066
【数６５】
2014125815.tif 000067
（式Ａ）から（式Ｃ）中において、ｘ、ｙ１、ｙ２、ｚ１、ｚ２は任意曲線上の座標、ξ、η、ζは領域中の任意の散乱点のｘ座標、ｙ座標、ｚ座標、ｃは伝播速度、ｋは波数であり、ω＝ｃｋである。
【請求項２】
物体に放射した波動の散乱波を解析する散乱トモグラフィ方法であって、
前記波動を、曲面上に配置された複数の送信アンテナ素子から前記物体に放射するステップと、
前記散乱波を、曲面上に配置された複数の受信アンテナ素子により受信するステップと、
前記受信アンテナ素子により受信した散乱波を示す散乱波データから、前記物体の内部情報に関する画像を再構成するステップとを含み、
前記画像を再構成するステップにおいて、
前記物体の内部情報に関する画像を再構成するための、（式Ｄ）に定義される関数φをあらかじめ設定し、
前記関数φの漸近式が満たす、（式Ｅ）に定義される方程式を構築し、
測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、（式Ｆ）に定義される映像化関数ρを導出し、
前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する
散乱トモグラフィ方法。
【数６６】
2014125815.tif 000068
【数６７】
2014125815.tif 000069
【数６８】
2014125815.tif 000070
（式Ｄ）から（式Ｆ）中において、ｘ１、ｘ２、ｙ１、ｙ２、ｚ１、ｚ２は任意曲面上の座標、ξ、η、ζは領域中の任意の散乱点のｘ座標、ｙ座標、ｚ座標、ｃは伝播速度、ｋは波数であり、ω＝ｃｋである。
【請求項３】
前記映像化関数は、高速フーリエ変換を用いて導出される
請求項１または２に記載の散乱トモグラフィ方法。
【請求項４】
前記波動は、電磁波または超音波である
請求項１〜３のいずれか１項に記載の散乱トモグラフィ方法。
【請求項５】
前記波動は、パルス波または所定の周波数を有する周期波である
請求項１〜４のいずれか１項に記載の散乱トモグラフィ方法。
【請求項６】
物体に放射した波動の散乱波を解析する散乱トモグラフィ装置であって、
曲面上に配置され、物体に波動を放射する複数の送信アンテナ素子と、
曲面上に配置され、前記放射された波動が前記物体において散乱した散乱波を受信する複数の受信アンテナ素子と、
前記受信した散乱波を示す散乱波データから、前記物体の内部情報に関する画像を再構成する画像再構成部とを備え、
前記画像再構成部は、
前記物体の内部情報に関する画像を再構成するための、（式Ａ）に定義される関数φをあらかじめ設定し、
前記関数φの漸近式が満たす、（式Ｂ）に定義される方程式を構築し、
測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、（式Ｃ）に定義される映像化関数ρを導出し、
前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する
散乱トモグラフィ装置。 
【数６９】
2014125815.tif 000071
【数７０】
2014125815.tif 000072
【数７１】
2014125815.tif 000073
（式Ａ）から（式Ｃ）中において、ｘ、ｙ１、ｙ２、ｚ１、ｚ２は任意曲線上の座標、ξ、η、ζは領域中の任意の散乱点のｘ座標、ｙ座標、ｚ座標、ｃは伝播速度、ｋは波数であり、ω＝ｃｋである。
【請求項７】
物体に放射した波動の散乱波を解析する散乱トモグラフィ装置であって、
曲面上に配置され、物体に波動を放射する複数の送信アンテナ素子と、
曲面上に配置され、前記放射された波動が前記物体において散乱した散乱波を受信する複数の受信アンテナ素子と、
前記受信した散乱波を示す散乱波データから、前記物体の内部情報に関する画像を再構成する画像再構成部とを備え、
前記画像再構成部は、
前記物体の内部情報に関する画像を再構成するための、（式Ｄ）に定義される関数φをあらかじめ設定し、
前記関数φの漸近式が満たす、（式Ｅ）に定義される方程式を構築し、
測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、（式Ｆ）に定義される映像化関数ρを導出し、
前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する
散乱トモグラフィ装置。 
【数７２】
2014125815.tif 000074
【数７３】
2014125815.tif 000075
【数７４】
2014125815.tif 000076
（式Ｄ）から（式Ｆ）中において、ｘ１、ｘ２、ｙ１、ｙ２、ｚ１、ｚ２は任意曲線上の座標、ξ、η、ζは領域中の任意の散乱点のｘ座標、ｙ座標、ｚ座標、ｃは伝播速度、ｋは波数であり、ω＝ｃｋである。
【請求項８】
前記画像再構成部は、高速フーリエ変換を用いて前記映像化関数を導出する
請求項６または７に記載の散乱トモグラフィ装置。
【請求項９】
前記波動は、電磁波または超音波である
請求項６〜８のいずれか１項に記載の散乱トモグラフィ装置。
【請求項１０】
前記波動は、パルス波または所定の周波数を有する周期波である
請求項６〜９のいずれか１項に記載の散乱トモグラフィ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、波動を利用して対象物の内部情報を取得し映像化（画像化）する技術に関し、特に、波動の散乱現象を利用した散乱トモグラフィ方法および散乱トモグラフィ装置に関する。
【背景技術】
【０００２】
従来、生体や建築物等の物体内部の情報を映像化する方法として、Ｘ線ＣＴ(Ｘ線トモグラフィ)、ＭＲＩ、ＰＥＴなどの方法が用いられている。具体的には、光、テラヘルツ、ミリ波、マイクロ波のような電磁波あるいは超音波、音波、弾性波などのフォノンのような波動を観測対象物である生体や物体、あるいはプラズマに照射し、その散乱波（反射波）を観測し解析することで、生体内部や固体内部、あるいは、プラズマの内部の情報を映像化している。また、最近では、波動に代えて磁場を用いて、生体や物体内部の情報の映像化も行われている。【０００３】
一般的に、これらの方法では、物体Ｏに電磁波や超音波のような波動ｕを照射し、物体Ｏの回りの多数の場所で物体Ｏから散乱される散乱波ｐを観測し、得られたデータを映像化するという技術が採用されている（例えば、特許文献１〜３、非特許文献１〜３参照）。【０００４】
特許文献１に記載の技術は、電波を用いて物体内部の情報の映像化を行うものである。円周上に配置されたセンサ素子で観測される散乱波のデータを、導電率や誘電率等のパラメータで修正しながら、データの取得を繰り返し、映像化するものである。【０００５】
特許文献２に記載の技術は、超音波フェーズドアレイ技術である。超音波センサ素子で受信した超音波のデータを平均ベクトルで補正し、映像化するものである。【０００６】
特許文献３に記載の技術は、物体に超音波を平面的に照射し、センサ素子で受信した超音波のデータを映像化するものである。曲面形状の物体内部の情報については、観測データ数を増加したり近似や補正を行ったりすることで映像化する。【０００７】
非特許文献１に記載の技術は、マルチパスリニアアレイレーダに関する技術であり、コンクリート内部の欠陥等の情報の映像化を行うものである。対象物の表面にセンサ素子を配置し、照射された波動の散乱波をセンサで観測し、観測データを解析し映像化するものである。【０００８】
非特許文献２に記載の技術は、超音波を用いて物体内部の情報の映像化を行うものである。物体に超音波を照射したときの散乱波を観測し、観測データをボルン近似（キルヒホフ近似）して映像化するものである。【０００９】
非特許文献３に記載の技術は、曲面上に設置したセンサ素子で得られる散乱波を観測し、センサ素子ごとに得られる時間波形データを、事前に取得したデータと比較し映像化するものである。
【先行技術文献】
【特許文献】
【００１０】
【特許文献１】特開２００７−１７７６５６号公報
【特許文献２】特開２００９−２８８１２９号公報
【特許文献３】特開２００４−５１２１１７号公報
【非特許文献】
【００１１】
【非特許文献１】トンネル覆工コンクリート検査用３次元映像化レーダを開発、三井造船技報、Ｎｏ．１８４、ｐ２４、２００５年２月
【非特許文献２】リニアアレイ探触子を用いた欠陥形状イメージング法の開発、日本機械学会、講演論文集、Ｐ６７９、２００５年９月
【非特許文献３】可撓性アレイプローブを用いた不規則面からの内部きずの超音波映像化、信学技報２０１２年６月
【発明の概要】
【発明が解決しようとする課題】
【００１２】
ここで、物体に波動を照射した際に生じる散乱現象は、演算子を用いて表すことができる。例えば、物体Ｏ、照射波ｕ、観測データｐで表される物理方程式を、演算子Ａを用いてｐ＝Ａｕ［Ｏ］と表すことができる。ここで、物体Ｏ、照射波ｕ、演算子（システム関数）Ａが既知である場合に、観測データｐを求める問題は、順方向問題と呼ばれている。順方向問題は数学的基盤が確かな方法であり、その解を求めることは普通の物理の教科書に書いてある方法で可能である。【００１３】
これに対して、医療や産業上で重要な課題は、照射波ｕ、システム関数Ａ、観測データｐが既知である場合に物体Ｏが何であるかを求める問題である。この問題は、物理現象の因果関係を逆方向に辿るという意味で、逆方向問題と呼ばれ、Ｏ＝Ａｕ−１［ｐ］と表すことができる。波動を物体に照射したときの散乱波を観測し解析して物体内部の情報を映像化するという方法（散乱トモグラフィ）は、この逆方向問題を利用するものである。【００１４】
逆方向問題は、順方向問題のように数学的基盤が確かな方法ではなく、現在まで確立した理論というものがないため、容易に解く事ができないという問題がある。したがって、逆方向問題を利用して物体内部の情報を映像化する方法では、例えば、物体の曲面形状等の条件が変わるたびに、理論や装置内部の構造等を変えてデータを再取得したり、取得したデータを補正したりする必要がある。そのため、逆方向問題を利用して物体内部の情報を映像化する方法は、汎用的な使用が難しい。また、データを再取得したり補正したりする必要があるため、計算速度の遅延やメモリの大量使用という問題がある。【００１５】
そこで、本発明は、逆問題の解析を汎用的かつ高速に行い、物体内部の情報を簡便に映像化することができる散乱トモグラフィ方法および散乱トモグラフィ装置を提供することを目的とする。
【課題を解決するための手段】
【００１６】
上記の課題を解決するため、本発明に係る散乱トモグラフィ方法は、物体に放射した波動の散乱波を解析する散乱トモグラフィ方法であって、前記波動を、曲面上に配置された複数の送信アンテナ素子から前記物体に放射するステップと、前記散乱波を、曲面上に配置された複数の受信アンテナ素子により受信するステップと、前記受信アンテナ素子により受信した散乱波を示す散乱波データを境界条件として用い、前記物体の内部情報に関する画像を再構成するステップとを含み、前記画像を再構成するステップにおいて、前記物体の内部情報に関する画像を再構成するための、後述する（式５）に定義される関数φをあらかじめ設定し、前記関数φの漸近式が満たす、後述する（式１０）に定義される方程式を構築し、測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、後述する（式２４）に定義される映像化関数ρを導出し、前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する。【００１７】
また、本発明に係る散乱トモグラフィ方法は、物体に放射した波動の散乱波を解析する散乱トモグラフィ方法であって、前記波動を、曲面上に配置された複数の送信アンテナ素子から前記物体に放射するステップと、前記散乱波を、曲面上に配置された複数の受信アンテナ素子により受信するステップと、前記受信アンテナ素子により受信した散乱波を示す散乱波データから、前記物体の内部情報に関する画像を再構成するステップとを含み、前記画像を再構成するステップにおいて、前記物体の内部情報に関する画像を再構成するための、後述する（式２７）に定義される関数φをあらかじめ設定し、前記関数φの漸近式が満たす、後述する（式３９）に定義される方程式を構築し、測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、後述する（式５２）に定義される映像化関数ρを導出し、前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する。【００１８】
これにより、画像を再構成するステップにおいて、関数φの設定は、３次元空間を対象に行われるので、曲率の大きい曲面を有する対象物の内部情報をより精度よく高速に映像化することができる。【００１９】
これにより、任意曲面上にセンサ素子を配置した解析モデルにおいて逆問題用の偏微分方程式を設定し、これを解くことで、曲率の大きい曲面を有する対象物の内部情報を汎用的かつ高速に映像化することができる。【００２０】
また、前記映像化関数は、高速フーリエ変換を用いて導出されてもよい。【００２１】
これにより、解析データの取得を高速に行うことができるので、対象物の内部情報を高速に映像化することができる。【００２２】
また、前記波動は、電磁波または超音波であってもよい。【００２３】
これにより、汎用的かつ簡便な方法で、曲率の大きい曲面を有する対象物の内部情報を映像化することができる。【００２４】
また、前記波動は、パルス波または所定の周波数を有する周期波であってもよい。【００２５】
これにより、汎用的かつ簡便な方法で、曲率の大きい曲面を有する対象物の内部情報を映像化することができる。【００２６】
上記の課題を解決するため、本発明に係る散乱トモグラフィ装置は、物体に放射した波動の散乱波を解析する散乱トモグラフィ装置であって、曲面上に配置され、物体に波動を放射する複数の送信アンテナ素子と、曲面上に配置され、前記放射された波動が前記物体において散乱した散乱波を受信する複数の受信アンテナ素子と、前記受信した散乱波を示す散乱波データから、前記物体の内部情報に関する画像を再構成する画像再構成部とを備え、前記画像再構成部は、前記物体の内部情報に関する画像を再構成するための、後述する（式５）に定義される関数φをあらかじめ設定し、前記関数φの漸近式が満たす、後述する（式１０）に定義される方程式を構築し、測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、後述する（式２４）に定義される映像化関数ρを導出し、前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する。【００２７】
これにより、任意曲面上にセンサ素子を配置した解析モデルにおいて逆問題用の偏微分方程式を設定し、これを解くことで、曲率の大きい曲面を有する対象物の内部情報を汎用的かつ高速に映像化することができる。【００２８】
また、本発明に係る散乱トモグラフィ装置は、物体に放射した波動の散乱波を解析する散乱トモグラフィ装置であって、曲面上に配置され、物体に波動を放射する複数の送信アンテナ素子と、曲面上に配置され、前記放射された波動が前記物体において散乱した散乱波を受信する複数の受信アンテナ素子と、前記受信した散乱波を示す散乱波データから、前記物体の内部情報に関する画像を再構成する画像再構成部とを備え、前記画像再構成部は、前記物体の内部情報に関する画像を再構成するための、後述する（式２７）に定義される関数φをあらかじめ設定し、前記関数φの漸近式が満たす、後述する（式３９）に定義される方程式を構築し、測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、後述する（式５２）に定義される映像化関数ρを導出し、前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する。【００２９】
これにより、画像を再構成するステップにおいて、関数φの設定は、３次元空間を対象に行われるので、曲率の大きい曲面を有する対象物の内部情報をより精度よく高速に映像化することができる。【００３０】
また、前記画像再構成部は、高速フーリエ変換を用いて前記映像化関数を導出してもよい。【００３１】
これにより、解析データの取得を高速に行うことができるので、対象物の内部情報を高速に映像化することができる。【００３２】
また、前記波動は、電磁波または超音波であってもよい。【００３３】
これにより、汎用的かつ簡便な方法で、曲率の大きい曲面を有する対象物の内部情報を映像化することができる。【００３４】
また、前記波動は、パルス波または所定の周波数を有する周期波であってもよい。【００３５】
これにより、汎用的かつ簡便な方法で、曲率の大きい曲面を有する対象物の内部情報を映像化することができる。
【発明の効果】
【００３６】
本発明により、逆問題の解析を汎用的かつ高速に行い、物体内部の情報を簡便に映像化することができる。
【図面の簡単な説明】
【００３７】
【図１】図１は、本発明の基礎となった技術について説明するための解析モデルである。
【図２】図２は、図１に示した解析モデルをさらに簡略化した図である。
【図３】図３は、実施の形態１に係るマルチパスアレイレーダの構成を示す概略図である。
【図４】図４は、実施の形態１に係るマルチパスアレイレーダの動作を示すフローチャートである。
【図５】図５は、実施の形態１に係る散乱トモグラフィ方法の原理を説明するための解析モデルである。
【図６】図６は、実施の形態１に係る散乱トモグラフィ方法の原理を説明するための解析モデルである。
【図７】図７は、実施の形態１に係る散乱トモグラフィ方法の原理を説明するための解析モデルである。
【図８】図８は、実施の形態１に係る散乱トモグラフィ方法の原理を説明するための解析モデルである。
【図９】図９は、実施の形態１に係る散乱トモグラフィ方法による解析のための解析モデルの例である。
【図１０Ａ】図１０Ａは、図９に示した解析モデルにおける散乱波データである。
【図１０Ｂ】図１０Ｂは、図９に示した解析モデルにおける散乱波データである。
【図１１Ａ】図１１Ａは、図９に示した解析モデルにおいてｃｕｒｖｆａｃｔｏｒ０のときのＭＰＬＡ  Ｒａｄａｒによる観測データである。
【図１１Ｂ】図１１Ｂは、図９に示した解析モデルにおいてｃｕｒｖｆａｃｔｏｒ０のときのＭＰＣＬＡ  Ｒａｄａｒによる観測データである。
【図１２Ａ】図１２Ａは、図９に示した解析モデルにおいてｃｕｒｖｆａｃｔｏｒ７．０２のときのＭＰＬＡ  Ｒａｄａｒによる観測データである。
【図１２Ｂ】図１２Ｂは、図９に示した解析モデルにおいてｃｕｒｖｆａｃｔｏｒ７．０２のときのＭＰＣＬＡ  Ｒａｄａｒによる観測データである。
【図１３Ａ】図１３Ａは、図９に示した解析モデルにおいてｃｕｒｖｆａｃｔｏｒ１４．０４のときのＭＰＬＡ  Ｒａｄａｒによる観測データである。
【図１３Ｂ】図１３Ｂは、図９に示した解析モデルにおいてｃｕｒｖｆａｃｔｏｒ１４．０４のときのＭＰＣＬＡ  Ｒａｄａｒによる観測データである。
【図１４】図１４は、実施の形態２に係る散乱トモグラフィ方法の原理を説明するための解析モデルである。
【図１５】図１５は、実施の形態２に係る散乱トモグラフィ方法の原理を説明するための解析モデルである。
【発明を実施するための形態】
【００３８】
（本発明の基礎となった知見）
本発明に係る実施の形態について説明する前に、本発明の基礎となった技術について説明する。
【００３９】
図１および図２は、本発明の基礎となった技術について説明をするための図である。
【００４０】
上記した、波動の散乱波を用いて非破壊で対象物の欠陥等を検出する方法（散乱トモグラフィ）として、Ｍｕｌｔｉ−Ｐａｔｈ  Ｌｉｎｅａｒ  Ａｒｒａｙ  Ｒａｄａｒ（ＭＰＬＡ  Ｒａｄａｒ）がある。この方法は、例えば、図１に示すように、センサとしてトンネル等の対象物にアンテナ素子１０を取り付け、当該アンテナ素子１０から電磁波を放射し、対象物において反射した反射波（散乱波）をアンテナ素子１０で検出し、放射した電磁波と反射波との関係から、非破壊で対象物の欠陥等を検出するものである。
【００４１】
ここで、ＭＰＬＡ  Ｒａｄａｒは、図２に示すように、曲率のある対象物であっても平面として近似して計算するものである。すなわち、ＭＰＬＡ  Ｒａｄａｒは、アンテナ素子１０が平面状に配列されていることを前提として理論が構築されているものであり、アンテナ素子１０を曲率の大きな曲面上に配置すると、電磁波の位相に影響を及ぼす誤差が増加し、計算により得られた像が焦点を結ばなくなるといった問題がある。
【００４２】
したがって、対象物が例えば新幹線のトンネル（曲率半径は約４．８ｍ）のように曲率が小さい（曲率半径が大きい）物体である場合には、平面として近似して計算しても誤差は小さく問題とならないが、対象物が在来線のトンネル（曲率半径は２．２〜２．８ｍ）のように曲率が大きい（曲率半径が小さい）物体である場合には、平面として近似して計算することは難しい。
【００４３】
平面として近似して計算した場合には、アンテナ上の送受素子の組み合わせによっては、位相に１波長相当以上の誤差が発生することもあり問題となる。例えば、図２の解析モデルにおいては、トンネルを構成する材料の比誘電率をε＝６．０、アレイサイズＬ＝１０００ｍｍ、電波の有効波長λ＝２０．４ｍｍ、電波の最高周波数３ＧＨｚと仮定すると、在来線のトンネルの半径Ｒ＝２．２ｍに対し、近似の際の誤差はΔｚ＝５７．６ｍｍとなり、電波の有効波長の約３倍の誤差が生じる。
【００４４】
そこで、以下、ＭＰＬＡ  Ｒａｄａｒを改良し、曲率が大きい対象物であっても誤差を低減して欠陥等の検出をすることができる散乱トモグラフィである、Ｍｕｌｔｉ−Ｐａｔｈ  ＣｕｒｖｉＬｉｎｅａｒ  Ａｒｒａｙ  Ｒａｄａｒ（ＭＰＣＬＡ  Ｒａｄａｒ）について説明する。ＭＰＣＬＡ  Ｒａｄａｒにより、逆問題の解析を汎用的かつ高速に行い、物体内部の情報を簡便に映像化することができる。
【００４５】
以下、図面を用いて、本発明に係る実施の形態について説明する。なお、図面において、同一の符号が付された構成要素は、同一または同種の構成要素を示す。
【００４６】
また、以下で説明する実施の形態は、本発明の好ましい一具体例を示す。以下の実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置および接続形態、ステップ、ステップの順序等は、一例であり、本発明を限定する主旨ではない。また、以下の実施の形態における構成要素のうち、本発明の最上位概念を示す独立請求項に記載されていない構成要素については、より望ましい形態を構成する任意の構成要素として説明される。
【００４７】
なお、一般的に、散乱トモグラフィとは、波動の散乱波を用いて非破壊で対象物の欠陥等を検出する方法のことをいうが、本件では、散乱トモグラフィが方法であることを明確にするために、以下、この方法を散乱トモグラフィ方法と呼ぶこととする。また、散乱トモグラフィ方法を行うための装置を散乱トモグラフィ装置と呼ぶこととする。
【００４８】
（実施の形態１）
＜センサアレイの構成＞
実施の形態１では、散乱トモグラフィ方法を行うためのセンサとして、一次元センサアレイを用いる場合について説明する。一次元センサアレイは、一次元的に配置された送信アンテナ素子および受信アンテナ素子（リニアアレイアンテナ）で構成されるセンサアレイである。
【００４９】
本実施の形態では、トンネル覆工コンクリートの内部情報を映像化する例を挙げるため、リニアアレイアンテナを構成する送信アンテナ素子および受信アンテナ素子は、トンネルの周方向に沿った一部にそれぞれ一列（一次元的）に配置され、トンネルの奥行き方向にリニアアレイアンテナを並進スキャンする場合について説明する。
【００５０】
はじめに、本実施の形態に係る散乱トモグラフィ装置であるマルチパスアレイレーダ２０の構成について説明する。マルチパスアレイレーダ２０は、送信用の複数のアンテナ素子と受信用の複数のアンテナ素子がそれぞれ一列に配置されたリニアアレイアンテナである。このリニアアレイアンテナにおいて、送信アンテナ素子列上の任意の素子で波動を送信し、受信アンテナ素子列上の任意の素子で反射波を受信する。例えば、送信アンテナ素子及び受信アンテナ素子がそれぞれｙ方向にｎｙ個配置された場合、ｎｙ２組の時系列データを得ることができる。
【００５１】
図３は、実施の形態１に係るマルチパスアレイレーダ２０の構成を示す概略図である。図３に示すように、マルチパスアレイレーダ２０は、センサ部３０と、画像再構成部４０とを備えている。
【００５２】
センサ部３０は、送信部３１と、受信部３２と、送信アンテナ素子３６と、受信アンテナ素子３７とを備えている。
【００５３】
送信アンテナ素子３６および受信アンテナ素子３７は、図３に示すように、それぞれ直線状に複数配置されてリニアアレイアンテナを構成している。送信アンテナ素子３６は、対象物に向かって波動を送信し、受信アンテナ素子３７は、対象物から反射（散乱）した反射波（散乱波）を受信する。ここでは、波動として電磁波を例として説明するが、電磁波に限らず超音波であってもよい。
【００５４】
送信部３１は、送信アンテナ素子３６から放射される電磁波の送信タイミングや回数、送信ゲインを調整する。
【００５５】
受信部３２は、受信アンテナ素子３７で受信した電磁波の散乱波のデータを、画像再構成部４０に搬送する。このとき、受信した散乱波のデータは、受信部３２により増幅されたりＡＤ変換等の信号処理をされてもよい。
【００５６】
画像再構成部４０は、受信部３２から搬送された散乱波のデータを解析し、後に説明する画像再構成アルゴリズムにより、散乱波のデータを映像化する。これにより、対象物の内部情報に応じた映像がモニタ５０に再生される。
【００５７】
なお、マルチパスアレイレーダ２０において、上記したリニアアレイアンテナがｘ軸方向へ移動することにより、ｎｘｎｙ２組の時系列データが得られる。時系列がｎｔ個の場合、得られるデータは総じてｎｘｎｙ２ｎｔ個となる。こうして得られたｎｘｎｙ２ｎｔ個のデータは、３次元の映像化に必要な量のデータｎｘｎｙｎｔ個に比べｎｙ倍の冗長度の情報量を持っている。したがって、このデータを、後に説明する画像再構成アルゴリズムにおける映像化関数に対応させることにより、対象物の内部情報に応じた３次元の画像を再生することができる。
【００５８】
また、送信アンテナ素子３６と受信アンテナ素子３７との間隔が自由に選択できるので、送信アンテナ素子３６および受信アンテナ素子３７の対に応じて波動のゲインを変える（経路依存可変ゲイン増幅機能）ことにより、対象物において検査可能な深さを変えることができる。
【００５９】
以下、散乱トモグラフィ方法による物体の内部情報の映像化の手順について説明する。
【００６０】
＜散乱トモグラフィ方法の手順＞
図４は、図３に係るマルチパスアレイレーダ２０の動作を示すフローチャート図である。
【００６１】
図４に示すように、散乱トモグラフィ方法による物体の内部情報を映像化（画像化）する方法は、以下のとおりである。
【００６２】
はじめに、送信アンテナ素子３６から対象物であるトンネル覆工コンクリートに向けて波動を放射する（Ｓ１１）。波動としては、例えば、超音波を用いる。送信部３１において超音波の波長、振幅等の調整を行い、送信アンテナ素子３６からコンクリートに放射する。
【００６３】
次に、コンクリート内部の欠陥で反射した反射波（散乱波）を受信アンテナ素子３７により受信する（Ｓ１２）。受信した散乱波は、受信部３２において、増幅されたりＡＤ変換等、画像再構成部４０における解析に適した形式に変換されてもよい。
【００６４】
次に、受信部３２から画像再構成部４０に、受信した散乱波を示す散乱波データが搬送される。画像再構成部４０では、搬送された散乱波データの解析が行われる（Ｓ１３）。ここでは、以下に説明する画像再構成アルゴリズムにより散乱波データの解析が行われる。具体的には、映像化関数が導出される。これにより、コンクリート内の欠陥に対応した映像（画像）が再構成される（Ｓ１４）。
【００６５】
さらに、再構成された映像のデータは、画像再構成部４０からモニタ５０に搬送され、モニタ５０に表示される。
【００６６】
以下、画像再構成部４０で行われる画像再構成アルゴリズムについて説明する。この画像再構成アルゴリズムは、本発明に係る散乱トモグラフィ方法の原理である。
【００６７】
＜画像再構成アルゴリズム＞
図５および図６は、本実施の形態に係る散乱トモグラフィ方法の原理を説明するための解析モデルである。以下、図５および図６に示すモデルを解析モデルとして、センサ部３０が一次元センサアレイである場合の映像化関数の導出について説明する。
【００６８】
以下に説明する画像再構成アルゴリズム（理論）は、曲線状に波動の放射点（送信アンテナ素子３６）と受信点（受信アンテナ素子３７）とを任意に設置し、放射点から放射された放射波を示す送信データと受信点で受信された散乱波を示す散乱波データを用いて、物体内部の情報を映像化するものである。
【００６９】
なお、図６に示すように、ｘ軸方向とは、トンネルの奥行き方向に対応し、ｘ軸方向に直交しｘ軸方向に沿って移動する２次元断面をｙｚ平面とする。つまり、この解析モデルでは、ｙｚ平面上に、曲線（円周の一部）状に送信アンテナ素子３６および受信アンテナ素子３７が設けられる。
【００７０】
画像再構成アルゴリズムについて数学的に簡単に説明すると、映像化に必要な漸近的な解（関数）をはじめに設定し、その解から方程式を構築し、送信データおよび受信データから、より厳密な解を求めるという、逆問題を解くものである。
【００７１】
すなわち、はじめに映像化に必要なグリーン関数を設定する。そして、このグリーン関数の解となる関数を導入する。そして、この関数が漸近的な解（関数）となる６次元空間（ｔ,ｘ,ｙ１,ｙ２,ｚ１,ｚ２）に関する偏微分方程式を構築する。そして、曲線状に配置された送信アンテナ素子３６から放射された送信データと曲線状に配置された受信アンテナ素子３７で受信された受信データとを境界条件としてこの偏微分方程式を解き、ｔ→０、ｘ→ｘ、ｚ２→ｚ１(＝ｚ)、ｙ１→ｙ２(＝ｙ)とおくことにより厳密な解（関数）を求める。これにより、厳密な解（関数）からなる映像化関数が求まり、コンクリート内部の欠陥等、物体内部の情報を汎用的かつ高速に映像化することができる。
【００７２】
具体的には、以下のとおりである。
【００７３】
１．散乱の逆問題とグリーン関数
図５により、ｒ１から出た波動が点ξで反射して点ｒ２へ戻ってくる状況を考える。周波数ωが一定という条件で、波動の送信点ｒ１と受信点ｒ２がある拘束条件を満たしながらｘ断面Ｄの内部を自由に動く。このとき得られるデータをＧ(ｒ１,ｒ２,ω)と書くと、この関数は領域内の反射点の分布に関係したものとなる。ここで、ωは角周波数であり、ω＝２πｆである。Ｇ(ｒ１,ｒ２,ω)は、全ての点ξからの反射信号の和であり、領域内には多くの反射点があるので、Ｇ(ｒ１,ｒ２,ω)は次の（式１）のように書くことができる。
【００７４】
【数１】
2014125815.tif 000003
【００７５】
ここで、
【数２】
2014125815.tif 000004
は、ｒ１から出て点ξで反射して点ｒ２へ戻る波動の信号強度を示す。
【００７６】
波動の送信点ｒ１と受信点ｒ２に課せられた拘束条件とは、ｒ１とｒ２のｘ座標が常に等しいという条件である。
【００７７】
この関数Ｇ(ｒ１,ｒ２,ω)を用いて散乱の逆問題の理論的構造を説明する。
【００７８】
３次元空間の部分領域をＤとし、その境界を∂Ｄとする。このとき、関数Ｇ(ｒ１,ｒ２,ω)は、領域Ｄの内部で次の（式２）のような微分方程式の解となっている。
【００７９】
【数３】
2014125815.tif 000005
【００８０】
さらに、Ｇ(ｒ１,ｒ２,ω)の境界∂Ｄでの値は、センサ部３０における測定値（送信データおよび受信データ）とする。上記方程式をこの境界条件のもとに解き、その結果からρ（ｒ）を次の（式３）のように定義する。
【００８１】
【数４】
2014125815.tif 000006
【００８２】
ここで、ρ（ｒ）とは、求めようとする領域Ｄ内の誘電率の勾配に関係した関数である。実際には、ここに現れる微分作用素L(∂/∂t,∂/∂r１,∂/∂r２)を求める必要がある。
【００８３】
２．任意曲面における散乱の逆問題L(∂/∂t,∂/∂r１,∂/∂r２)の導出
次に、上記した微分作用素を求める方法について説明する。図６〜図８は、微分作用素を求める方法について説明するための解析モデルである。
【００８４】
図６に示す一次元センサアレイにおいて、一次元センサアレイが回転スキャンしながらｘ軸方向へ移動する場合を考える。図７および図８は、図６に示すセンサアレイを部分的に拡大した図である。
【００８５】
図７および図８に示すように、点Ｐ１から放射した波動は、点Ｐで反射して点Ｐ２へ帰ってくる。また、計測点は断面曲線Ｓ上を動くと仮定する。
【００８６】
任意曲線上ではｒ１、ｒ２のｙｚ座標は、ｒ１＝（ｘ,ｙ１,ｚ１)、ｒ２＝（ｘ,ｙ２,ｚ２)と表される。このとき、関数Ｇは次の（式４）のように定義される。
【００８７】
【数５】
2014125815.tif 000007
【００８８】
次に、関数Ｇ(ｒ１,ｒ２,ω)を満たす方程式として、次の（式５）のような関数φを導入する。ここで、ｃは伝播速度、ｋは波数であり、波長をλとすると、ω＝ｃｋ、ｋ＝２π／λの関係がある。なお、（式５）に示すφは、本発明における物体の内部情報に関する画像を再構成するための関数（解）、すなわち、本発明における（式Ａ）に定義される関数に相当する。ξ、η、ζは、図６から図８に示す点Ｐ（ξ，η，ζ）のｘ座標、ｙ座標およびｚ座標である。なお、点Ｐ（ξ，η，ζ）は、領域中の任意の散乱点である。
【００８９】
【数６】
2014125815.tif 000008
【００９０】
ここでは、時間の因子はｅｘｐ（−ｉωｔ)に比例すると仮定し、波数をｋとしている。上式の被積分項にある関数は（式４）の
【数７】
2014125815.tif 000009
であり、以下の（式６）のように表される。
【００９１】
【数８】
2014125815.tif 000010
【００９２】
次に、（式５）が漸近的な解（関数）となる偏微分方程式を求める。そのために、微分した結果生ずる１／ρに関して高次の項を無視して計算する。なお、本発明において、漸近的な解である関数のことを漸近式と呼んでいる。以下、微分の略記法を以下の（式７）のように定義する。
【００９３】
【数９】
2014125815.tif 000011
【００９４】
このとき、φは次の（式８）に示す方程式を満たす。この（式８）が、本発明における漸近式に相当する。
【００９５】
【数１０】
2014125815.tif 000012
【００９６】
（式８）は定常状態の場合を仮定して導かれているが、（式８）を非定常状態の場合へ拡張することは容易である。そのために、次の（式９）のような変数の置換えをする。
【００９７】
【数１１】
2014125815.tif 000013
【００９８】
これにより、最終的に次の（式１０）が得られる。この（式１０）が、本発明における（式Ｂ）に定義される関数に相当する。
【００９９】
【数１２】
2014125815.tif 000014
【０１００】
次に、φの時間の因子をｅｘｐ(−ｉωｔ)に比例すると仮定して、（式１０）の解を求める。最初に、（式１０）におけるφをｔ,ｘ,ｙ１,ｙ２について多重フーリエ変換すると、以下の（式１１）のように表される。
【０１０１】
【数１３】
2014125815.tif 000015
【０１０２】
なお、高速フーリエ変換を用いることにより、解析データの取得を高速に行うことができる。したがって、解析時間を短縮することができる。
【０１０３】
ここで、ｚ１,ｚ２に関する偏微分をそれぞれＤｚ１,Ｄｚ２とおくと、以下の（式１２）が得られる。
【０１０４】
【数１４】
2014125815.tif 000016
【０１０５】
この方程式を解くのであるが、ｚ１,ｚ２の２変数がある。そのため、（式１２）を解くに当たり、固定された(ｘ,ｙ１,ｙ２)あるいは(ｋｘ,ｋｙ１,ｋｙ２)に関して、(ｚ１,ｚ２)空間の中の一次元の自由度を持った領域で与えられる境界条件が必要となる。
【０１０６】
ところが、レーダの計測により得られる境界条件は、(ｚ１,ｚ２)空間の中の１点(ｆ(ｙ１),ｆ(ｙ２))で与えられるのみである。
【０１０７】
したがって、この問題を解決するために、ｚ１＝ｚ，ｚ２＝ｚの場合の理論と本理論との整合性が要求される。すなわち、ｚ１,ｚ２を独立とする本理論が、ｚ１＝ｚ,ｚ２＝ｚという特別の解を含むことが要求される。
【０１０８】
そこで、（式１２）の解を、次の（式１３）のように仮定する。
【０１０９】
【数１５】
2014125815.tif 000017
【０１１０】
（式１３）は、ｚ１＝ｚ２＝ｚのとき、次の（式１４）のようになる。
【０１１１】
【数１６】
2014125815.tif 000018
【０１１２】
ここで、（式１３）を（式１２）へ代入すると、次の（式１５）が得られる。（式１５）が、１つの境界条件となる。
【０１１３】
【数１７】
2014125815.tif 000019
【０１１４】
さらに、境界条件としてもう１つの方程式が必要である。ここで、上記の整合性の要求から、次の（式１６）が導かれる。
【０１１５】
【数１８】
2014125815.tif 000020
【０１１６】
（式１５）と（式１６）から、ｓ１(ｋｘ,ｋｙ１,ｋｙ２)、ｓ２(ｋｘ,ｋｙ１,ｋｙ２)は次の（式１７）のように決定される。
【０１１７】
【数１９】
2014125815.tif 000021
【０１１８】
具体的な計算は後に行うことにし、このようなｓ１(ｋｘ,ｋｙ１,ｋｙ２)、ｓ２(ｋｘ,ｋｙ１,ｋｙ２)を用いて、（式１３）の方程式の解を次の（式１８）のように書くことができる。
【０１１９】
【数２０】
2014125815.tif 000022
【０１２０】
ここで、ｘを固定した断面曲線Ｓの方程式を、例えば次の（式１９）のように仮定する。本質的な仮定ではないが、例えば放物面を仮定する。
【０１２１】
【数２１】
2014125815.tif 000023
【０１２２】
このとき、断面曲線Ｓ上で与えられる境界条件は、次の（式２０）のようになる。
【０１２３】
【数２２】
2014125815.tif 000024
【０１２４】
この方程式は、ａ(ｋｘ,ｋｙ１,ｋｙ２)を決定するために用いられる。ここで、次の（式２１）に示す簡略記法を用いる。
【０１２５】
【数２３】
2014125815.tif 000025
【０１２６】
（式２１）に示す簡略記法により（式２０）は、次の（式２２）に示すａ(ｋｘ,ｋｙ１,ｋｙ２)に関する積分方程式となる。
【０１２７】
【数２４】
2014125815.tif 000026
【０１２８】
この式からａ(ｋ)を求めることができれば、映像化関数は（式１８）から次の（式２３）のように表される。
【０１２９】
【数２５】
2014125815.tif 000027
【０１３０】
（式２３）においてｚ１＝ｚ２＝ｚとおき、ｋに関してフーリエ変換すると、次の（式２４）のような関数が得られる。この（式２４）は、本発明における（式Ｃ）に定義される映像化関数に相当する。
【０１３１】
【数２６】
2014125815.tif 000028
【０１３２】
こうして、最終的な映像化関数ρ(ｒ)が得られる。
【０１３３】
３．積分方程式（式２４）の解
次に、上記した（式２４）の解を求める。すなわち、図８に示すように、断面曲面Ｓ上の点ＰＩ,ＰＪで測定されたデータからｘｙ面（平面ｚ＝０）上の点での近似データを求める。
【０１３４】
（式２４）におけるφについて、曲面上の点ＰＩ,ＰＪで測定されたデータφ(ｘ,ｙＩ,ｙＪ,ｚＩ,ｚＪ,ｔ)をフーリエ変換したΦ(ｋｘ,ｙＩ,ｙＪ,ｋ)を、次の（式２５）のようにおく。
【０１３５】
【数２７】
2014125815.tif 000029
【０１３６】
ここで、ｚＩ,ｚＪは、図８に示すように断面曲線上にあるので、次の（式２６）が成立する。
【０１３７】
【数２８】
2014125815.tif 000030
【０１３８】
さらに、（式１８）から次の（式２７）が得られる。
【０１３９】
【数２９】
2014125815.tif 000031
【０１４０】
（式２７）を、次の（式２８）のように解釈して書くことができる。
【０１４１】
【数３０】
2014125815.tif 000032
【０１４２】
（式２８）について両辺のフーリエ変換をとると、次の（式２９）が得られる。
【０１４３】
【数３１】
2014125815.tif 000033
【０１４４】
（式２９）の両辺を積分すると、
【数３２】
2014125815.tif 000034
【０１４５】
こうして、次の（式３１）に示すように、ａＩ,Ｊが得られる。
【０１４６】
【数３３】
2014125815.tif 000035
【０１４７】
全てのＩ,Ｊに関して和を求めると、次の（式３２）が得られる。
【０１４８】
【数３４】
2014125815.tif 000036
【０１４９】
こうして、平面ｚ＝０での境界条件へ変換することができる。
【０１５０】
ｚ＝０における境界条件から（式１０）の偏微分方程式の解を求めると、（式１８）から次の（式３３）のようになる。
【０１５１】
【数３５】
2014125815.tif 000037
【０１５２】
このとき、画像を再構成するための映像化関数は、次の（式３４）をｋについて積分することにより得られる。
【０１５３】
【数３６】
2014125815.tif 000038
【０１５４】
（式３４）をｋについて積分するに当たり、次のようにｋｚなる変数を導入する。（式３５）には、ｋをｋｚで表す式及びそれを微分した関数を示す。
【０１５５】
【数３７】
2014125815.tif 000039
【０１５６】
こうして、画像を再構成するための映像化関数ρ(ｘ,ｙ,ｚ)は、最終的に次の（式３６）のように表される。
【０１５７】
【数３８】
2014125815.tif 000040
【０１５８】
この映像化関数により、物体の内部情報に関する画像を再構成する。したがって、コンクリート内部の欠陥等、物体の内部情報に関する画像を汎用的かつ高速に映像化することができる。
【０１５９】
＜散乱トモグラフィ方法により得られる再構成画像＞
次に、上記した散乱トモグラフィ方法により得られる再構成画像について説明する。ここでは、図９に示すモデルを解析モデルとした再構成画像について説明する。
【０１６０】
図９は、本実施の形態に係る解析モデルを示す図である。図１０Ａおよび図１０Ｂは、生データをそのまま表示したＢモード画像であり、点状の反射体が双曲線状に大きく広がっている。図１１Ａ〜図１３Ｂは、従来のＭＰＬＡ  Ｒａｄａｒによる解析結果と本実施の形態に係るＭＰＣＬＡ  Ｒａｄａｒとを用いた解析結果である。
【０１６１】
図９に示す解析モデルは、トンネルを対象物とするものであり、トンネルの深さＺ＝Ｚ０の平面に点欠陥を３点有し、１次元のセンサアレイによって点欠陥を検出するものである。図９では、センサアレイ１０を放物線上に配置したと仮定したときのＣｕｒｖｆａｃｔｏｒをＣｕｒｖｆａｃｔｏｒ＝Ｎ＊Ｌ／（２Ｒ）と定義している。ここで、Ｎはｙ方向センサの数Ｎ＝ＮＹ、Ｌはセンサアレイの全長Ｌ＝Ｎ＊ｄｙ、Ｒはセンサアレイの曲率半径、ｄｘはｘ方向のスキャンピッチ、ｄｙはｙ方向のセンサ配置間隔を表している。
【０１６２】
図１０Ａ及び図１０Ｂに示す解析結果では、対象物であるトンネルのＣｕｒｖｆａｃｔｏｒを４．６８とし、Ｂモード、すなわち、受信アンテナ素子３７で受信した散乱波の強度をトンネルの断面図に対応する図として表示したものである。図１０Ａは、散乱体の深さを２００ｍｍ、すなわち、３つの欠陥が受信アンテナ素子３７が配置されたトンネル表面の位置より２００ｍｍ深い位置にある場合の散乱波を受信している。図１０Ｂでは、散乱体の深さを５００ｍｍ、すなわち、３つの欠陥が受信アンテナ素子３７が配置されたトンネル表面の位置より５００ｍｍ深い位置にある場合の散乱波を受信している。
【０１６３】
図１１Ａおよび図１１Ｂは、それぞれ、ｃｕｒｖｆａｃｔｏｒを０（平面）、散乱体の深さを５００ｍｍとしたときのＭＰＬＡ  ＲａｄａｒおよびＭＰＣＬＡ  Ｒａｄａｒによる解析結果である。図１１Ａおよび図１１Ｂによると、いずれの方法で解析した場合も３つの欠陥が計測され、ＭＰＬＡ  ＲａｄａｒおよびＭＰＣＬＡ  Ｒａｄａｒのいずれでもほとんど解析結果はほとんど変わらないことがわかる。
【０１６４】
また、図１２Ａおよび図１２Ｂは、それぞれ、ｃｕｒｖｆａｃｔｏｒを７．０２、散乱体の深さを５００ｍｍとしたときのＭＰＬＡ  ＲａｄａｒおよびＭＰＣＬＡ  Ｒａｄａｒによる解析結果である。図１２Ａおよび図１２Ｂによると、いずれの方法で解析した場合も３つの欠陥が計測されているが、図１２Ａに示すＭＰＬＡ  Ｒａｄａｒのほうでは像の焦点が合っておらずやや不明確な像となっている。一方、図１２Ｂに示すＭＰＣＬＡ  Ｒａｄａｒによる解析結果では、３つの欠陥が明確に計測されていることがわかる。
【０１６５】
また、図１３Ａおよび図１３Ｂは、それぞれ、ｃｕｒｖｆａｃｔｏｒを１４．０４、散乱体の深さを５００ｍｍとしたときのＭＰＬＡ  ＲａｄａｒおよびＭＰＣＬＡ  Ｒａｄａｒによる解析結果である。図１３Ａおよび図１３Ｂによると、いずれの方法で解析した場合も３つの欠陥が計測されているが、図１３Ａに示すＭＰＬＡ  Ｒａｄａｒのほうでは像の焦点が合っておらず不明確な像となっている。一方、図１３Ｂに示すＭＰＣＬＡ  Ｒａｄａｒによる解析結果では、３つの欠陥が明確に計測されていることがわかる。
【０１６６】
以上、本実施の形態に係る散乱トモグラフィ方法は、曲線状にセンサ素子を配置した解析モデルにおいて逆問題用の偏微分方程式を設定し、これを解くことで、映像化関数を得ている。これにより、物体に放射した波動の散乱波を解析する散乱トモグラフィ方法において、曲率の大きい曲面を有する対象物の内部情報を、汎用的かつ高速に映像化することができる。
【０１６７】
なお、マルチパスアレイレーダ２０において、上記したリニアアレイアンテナがｘ軸方向へ移動することにより、ｎｘｎｙ２組の時系列データが得られる。時系列がｎｔ個の場合、得られるデータは総じてｎｘｎｙ２ｎｔ個となる。こうして得られたｎｘｎｙ２ｎｔ個のデータは、３次元の映像化に必要な量のデータｎｘｎｙｎｔ個に比べｎｙ倍の冗長度の情報量を持っている。このように、計測データが大きな冗長性を持っているため、マルチパスアレイレーダ２０の出力はノイズに強い性質がある。
【０１６８】
なお、上記の演算式、および、演算式の導出手順は、一例であって、別の演算式、および、別の導出手順が用いられてもよい。
【０１６９】
また、本実施の形態では、波動として電磁波を用いたが、電磁波に限らず、超音波であってもよい。また、本実施の形態では、電磁波を用いているため、所定の周波数を有する周期波を用いているが、波動として周期波に限らずパルス波を用いてもよい。
【０１７０】
また、本実施の形態では、対象物としてトンネル覆工コンクリートを例に説明しているが、対象物はトンネル覆工コンクリートに限らず、他の物体であってもよい。
【０１７１】
（実施の形態２）
次に、本発明の第２の実施の形態について説明する。
【０１７２】
実施の形態２では、散乱トモグラフィ方法を行うためのセンサとして、二次元センサアレイを用いる場合について説明する。二次元センサアレイは、曲面状に二次元的に配置された送信アンテナ素子および受信アンテナ素子（マルチアレイアンテナ）で構成されるセンサアレイである。
【０１７３】
本実施の形態に係る散乱トモグラフィ装置であるマルチパスアレイレーダの構成は、実施の形態１にかかるマルチパスアレイレーダ２０とほぼ同様であるが、送信アンテナ素子および受信アンテナ素子が曲面状に二次元的に配置される点が実施の形態１にかかるマルチパスアレイレーダ２０と異なっている。したがって、本実施の形態に係るマルチパスアレイレーダの画像再構成部で行われる画像再構成アルゴリズムが、実施の形態１にかかるマルチパスアレイレーダ２０と異なっている。
【０１７４】
以下、本実施の形態に係る散乱トモグラフィ方法の画像再構成アルゴリズムについて説明する。
【０１７５】
＜画像再構成アルゴリズム＞
図１４および図１５は、本実施の形態に係る散乱トモグラフィ方法の原理を説明するための解析モデルである。以下、図１４および図１５に示すモデルを解析モデルとして、マルチパスアレイレーダのセンサ部が二次元センサアレイである場合の映像化関数の導出について説明する。
【０１７６】
以下に説明する画像再構成アルゴリズム（理論）は、曲面状に波動の放射点（一の送信アンテナ素子３６）と受信点（一の受信アンテナ素子３７）を任意に設置し、放射点から放射された放射波を示す送信データと受信点で受信された散乱波を示す散乱波データを用いて、物体内部の情報を映像化するものである。
【０１７７】
なお、図１５は、二次元センサアレイの模式図である。図１５に示すように、ｘ、ｙ座標は対象物の表面に置き、ｚ座標は対象物の表面の法線方向にとる。この解析モデルでは、ｘｙｚ空間に、曲面状に送信アンテナ素子３６および受信アンテナ素子３７が設けられる。
【０１７８】
画像再構成アルゴリズムについて数学的に簡単に説明すると、実施の形態１に係る散乱トモグラフィ方法と同様、映像化に必要な漸近的な解（関数）をはじめに設定し、その解から方程式を構築し、送信データおよび受信データから、より厳密な解を求めるという、逆問題を解くものである。
【０１７９】
すなわち、はじめに映像化に必要なグリーン関数を設定する。そして、このグリーン関数の解となる関数を導入する。ここで、次に、この関数が漸近的な解（関数）となる７次元空間（ｔ,ｘ１,ｘ２,ｙ１,ｙ２,ｚ１,ｚ２）に関する偏微分方程式を構築する。そして、曲面状に配置された送信アンテナ素子から放射された送信データと曲面状に配置された受信アンテナ素子で受信された受信データとを境界条件としてこの偏微分方程式を解き、ｔ→０、ｘ１→ｘ２(＝ｘ)、ｚ２→ｚ１(＝ｚ)、ｙ１→ｙ２(＝ｙ)とおくことにより厳密な解（関数）を求める。これにより、厳密な解（関数）からなる映像化関数が求まり、トンネル内部の欠陥等、物体内部の情報を汎用的かつ高速に映像化することができる。
【０１８０】
具体的には、以下のとおりである。
【０１８１】
１．散乱の逆問題とグリーン関数
図１４により、ｒ１から出た波動が点ξで反射して点ｒ２へ戻ってくる状況を考える。周波数ωが一定という条件で、波動の送信点ｒ１と受信点ｒ２がある拘束条件を満たしながら曲面上を自由に動く。このとき得られるデータをＧ(ｒ１,ｒ２,ω)と書くと、この関数は領域内の反射点の分布に関係したものとなる。このとき、ωは角周波数＝２πｆである。Ｇ(ｒ１,ｒ２,ω)は、全ての点ξからの反射信号の和であり、領域内に多くの反射点があるので、Ｇ(ｒ１,ｒ２,ω)は、実施の形態１に示した（式１）と同様に考えることができる。
【０１８２】
ただし、本実施の形態では、波動の送信点ｒ１と受信点ｒ２に課せられた拘束条件とは、ｒ１とｒ２がある曲面上を動くという条件である。
【０１８３】
この関数Ｇ(ｒ１,ｒ２,ω)を用いて散乱の逆問題の理論的構造を説明する。
【０１８４】
２．任意曲面における散乱の逆問題L(∂/∂t,∂/∂r１,∂/∂r２)の導出
以下、この微分作用素を求める方法について述べる。本実施の形態に係る解析モデルでは、実施の形態１の解析モデルとは異なり、送信点および受信点がある曲面上を動くということである。曲線上ではｒ１、ｒ２のｘ,ｙ,ｚ座標が必ずしも等しくない。具体的にはｒ１＝(ｘ１,ｙ１,ｚ１)、ｒ２＝(ｘ２,ｙ２,ｚ２)である。関数Ｇは実施の形態１に示した（式４）と同様に定義され、ｒ１＝(ｘ１,ｙ１,ｚ１)、ｒ２＝(ｘ２,ｙ２,ｚ２)から関数Ｇ(ｒ１,ｒ２,ω)の満たす方程式を求めると、以下の（式３７）のようになる。なお、（式３７）に示すφは、本発明における物体の内部情報に関する画像を再構成するための関数（解）、すなわち、本発明における（式Ｄ）に定義される関数に相当する。
【０１８５】
【数３９】
2014125815.tif 000041
【０１８６】
ここでは、時間の因子をｅｘｐ(−ｉωｔ)に比例すると仮定し、波数をｋとしている。ｘ１、ｘ２、ｙ１、ｙ２、ｚ１、ｚ２は任意曲線上の座標である。ｃは伝播速度であり、波長をλとすると、ω＝ｃｋ、ｋ＝２π／λの関係がある。ξ、η、ζは、任意の点の座標である。
【０１８７】
上式の被積分項にある関数は（式４）の
【数４０】
2014125815.tif 000042
であり、（式６）のように表される。
【０１８８】
ここで、（式３７）が漸近的な解（関数）となる偏微分方程式を求める。そのためには、微分した結果生ずる１／ρに関して高次の項を無視して計算する。なお、本発明において、漸近的な解である関数のことを漸近式と呼んでいる。以下、微分の略記法を以下の（式３８）ように定義する。
【０１８９】
【数４１】
2014125815.tif 000043
【０１９０】
このとき、実施の形態１と同様の計算により、φは次の（式３９）に示す方程式を満たす。この（式３９）が、本発明における漸近式に相当する。また、この（式３９）が、本発明における（式Ｅ）に定義される関数に相当する。
【０１９１】
【数４２】
2014125815.tif 000044
【０１９２】
次に、φの時間の因子をｅｘｐ(−ｉωｔ)に比例すると仮定して、（式３９）の解を求めることにする。まず、φをｔ,ｘ１,ｘ２,ｙ１,ｙ２について多重フーリエ変換すると、以下の（式４０）のように表される。
【０１９３】
【数４３】
2014125815.tif 000045
【０１９４】
なお、高速フーリエ変換を用いることにより、解析データの取得を高速に行うことができる。したがって、解析時間を短縮することができる。
【０１９５】
ここで、ｚ１,ｚ２に関する偏微分をそれぞれＤｚ１,Ｄｚ２とおくと、以下の（式４１）が得られる。
【０１９６】
【数４４】
2014125815.tif 000046
【０１９７】
この方程式を解くのであるが、ｚ１,ｚ２の２変数がある。そのため、（式４１）を解くに当たり、固定された(ｘ１,ｘ２,ｙ１,ｙ２)あるいは(ｋｘ１,ｋｘ２,ｋｙ１,ｋｙ２)に関して、(ｚ１,ｚ２)空間の中の一次元の自由度を持った領域で与えられる境界条件が必要となる。
【０１９８】
ところが、レーダの計測により得られる境界条件は(ｚ１,ｚ２)空間の中の１点{ｆ(ｘ１,ｙ１),ｆ(ｘ２,ｙ２)}で与えられるのみである。
【０１９９】
したがって、この問題を解決するために、ｚ１＝ｚ,ｚ２＝ｚの場合の理論と本理論との整合性が要求される。すなわち、ｚ１,ｚ２を独立とする本理論がｚ１＝ｚ,ｚ２＝ｚである特別の場合の解を含むことが要求される。
【０２００】
そこで、（式４１）の解を次の（式４２）のように仮定する。
【０２０１】
【数４５】
2014125815.tif 000047
【０２０２】
（式４２）は、ｚ１＝ｚ,ｚ２＝ｚのとき、次の（式４３）のようになる。
【０２０３】
【数４６】
2014125815.tif 000048
【０２０４】
一方、境界条件としてもう１つの方程式が必要である。ここで、上記の整合性の要求から、次の（式４４）が導かれる。
【０２０５】
【数４７】
2014125815.tif 000049
【０２０６】
（式４３）および（式４４）から、ｓ１,ｓ２は次の（式４５）のように決定される。
【０２０７】
【数４８】
2014125815.tif 000050
【０２０８】
このようなｓ１(ｋｘ,ｋｙ１,ｋｙ２),ｓ２(ｋｘ,ｋｙ１,ｋｙ２)を用いて方程式の解を次の（式４６）のように書くことができる。
【０２０９】
【数４９】
2014125815.tif 000051
【０２１０】
また、曲面Ｓの方程式を、例えば次の（式４７）のように仮定する。
【０２１１】
【数５０】
2014125815.tif 000052
【０２１２】
このとき、曲面Ｓ上で与えられる境界条件は、次の（式４８）のようになる。
【０２１３】
【数５１】
2014125815.tif 000053
【０２１４】
この方程式は、ａ(ｋｘ１,ｋｘ２,ｋｙ１,ｋｙ２)を決定するために用いられる。ここで、次の（式４９）に示す簡略記法を用いる。
【０２１５】
【数５２】
2014125815.tif 000054
【０２１６】
（式４９）に示す簡略記法により、（式４８）は次の（式５０）に示すａ(ｋｘ１,ｋｘ２,ｋｙ１,ｋｙ２)に関する積分方程式となる。
【０２１７】
【数５３】
2014125815.tif 000055
【０２１８】
この式からａ(ｋ)を求めることができれば、映像化は（式４６）から次の（式５１）のように表される。
【０２１９】
【数５４】
2014125815.tif 000056
【０２２０】
（式５１）において、ｚ１＝ｚ２＝ｚとおき、ｋに関してフーリエ変換すると、次の（式５２）のような関数が得られる。
【０２２１】
【数５５】
2014125815.tif 000057
【０２２２】
こうして最終的な映像化関数ρ(ｒ)が得られる。
【０２２３】
３．積分方程式（式５０）の解
次に、上記した（式５０）の解を求める。すなわち、曲面Ｓ上の点ＰＩ,ＰＪで測定されたデータから平面ｚ＝０上の任意の点での近似データを求める。
【０２２４】
（式５０）におけるφについて、曲面Ｓ上の点ＰＩ,ＰＪで測定されたデータφ(ｘＩ,ｘＪ,ｙＩ,ｙＪ,ｚＩ,ｚＪ,ｔ)をフーリエ変換したΦＩ，Ｊ(ｘＩ,ｘＪ,ｙＩ,ｙＪ,ｋ)を次の（式５３）のようにおく。
【０２２５】
【数５６】
2014125815.tif 000058
【０２２６】
さらに、（式４６）から次の（式５４）が得られる。
【０２２７】
【数５７】
2014125815.tif 000059
【０２２８】
（式５４）を、次の（式５５）のように解釈して書くことができる。
【０２２９】
【数５８】
2014125815.tif 000060
【０２３０】
（式５５）について両辺のフーリエ変換をとり、両辺を積分すると、次の（式５６）が得られる。
【０２３１】
【数５９】
2014125815.tif 000061
【０２３２】
こうして、次の（式５７）に示すように、ａＩ,Ｊが次のように得られる。
【０２３３】
【数６０】
2014125815.tif 000062
【０２３４】
全てのＩ,Ｊに関して和を求めると、次の（式５８）が得られる。
【０２３５】
【数６１】
2014125815.tif 000063
【０２３６】
こうして、平面ｚ＝０での境界条件へ変換することができる。
【０２３７】
ｚ＝０における境界条件から（式３９）の偏微分方程式の解を求めると、（式４６）に（式５８）を代入してから次の（式５９）のようになる。この（式５９）から上述した（式５２）が導かれ、これが本発明における（式Ｆ）に定義される映像化関数に相当する。
【０２３８】
【数６２】
2014125815.tif 000064
【０２３９】
この映像化関数により、物体の内部情報に関する画像を再構成する。したがって、物体内部の欠陥等、物体の内部情報に関する画像を汎用的かつ高速に映像化することができる。また、画像を再構成するステップにおいて、関数φの設定は、３次元空間を対象に行われるので、曲率の大きい曲面を有する対象物の内部情報をより精度よく高速に映像化することができる。
【０２４０】
以上、本実施の形態に係る散乱トモグラフィ方法は、曲面上に任意にセンサ素子を配置した解析モデルにおいて逆問題用の偏微分方程式を設定し、これを解くことで、映像化関数を得ている。これにより、物体に放射した波動の散乱波を解析する散乱トモグラフィ方法において、曲率の大きい曲面を有する対象物の内部情報を、汎用的かつ高速に映像化することができる。
【０２４１】
特に、本実施の形態では、マルチパスアレイレーダのセンサ部が二次元センサアレイである場合の散乱トモグラフィ方法における映像化関数を導出している。したがって、この映像化関数により、実施の形態１に示した、マルチパスアレイレーダのセンサ部が一次元センサアレイである場合よりも機械的な走査がないので高速に再構成画像を得ることができる。
【０２４２】
なお、上記の演算式、および、演算式の導出手順は、一例であって、別の演算式、および、別の導出手順が用いられてもよい。
【０２４３】
また、本実施の形態では、波動として電磁波を用いたが、電磁波に限らず、超音波であってもよい。また、本実施の形態では、電磁波を用いているため、所定の周波数を有する周期波を用いているが、波動として周期波に限らずパルス波を用いてもよい。
【０２４４】
また、本実施の形態では、対象物としてトンネル覆工コンクリートを例に説明しているが、対象物はトンネル覆工コンクリートに限らず、他の物体であってもよい。
【０２４５】
以上、本発明に係る散乱トモグラフィ方法および散乱トモグラフィ装置について、複数の実施の形態に基づいて説明したが、本発明は実施の形態に限定されるものではない。実施の形態に対して当業者が思いつく変形を施して得られる形態、および、複数の実施の形態における構成要素を任意に組み合わせて実現される別の形態も本発明に含まれる。
【０２４６】
例えば、散乱トモグラフィ装置において、特定の処理部が実行する処理を別の処理部が実行してもよい。また、散乱トモグラフィ方法において処理を実行する順番が変更されてもよいし、複数の処理が並行して実行されてもよい。
【０２４７】
また、本発明の散乱トモグラフィ方法におけるステップは、コンピュータによって実行されることとしてもよい。そして、本発明は、散乱トモグラフィ方法に含まれるステップを、コンピュータに実行させるためのプログラムとして実現できる。さらに、本発明は、そのプログラムを記録したＣＤ−ＲＯＭ等の非一時的なコンピュータ読み取り可能な記録媒体として実現できる。
【０２４８】
また、散乱トモグラフィ装置に含まれる複数の構成要素は、集積回路であるＬＳＩとして実現されてもよい。これらの構成要素は、個別に１チップ化されてもよいし、一部または全てを含むように１チップ化されてもよい。ここでは、ＬＳＩとしたが、集積度の違いにより、ＩＣ（Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）、システムＬＳＩ、スーパーＬＳＩまたはウルトラＬＳＩと呼称されることもある。
【０２４９】
また、集積回路化の手法はＬＳＩに限るものではなく、専用回路または汎用プロセッサで実現してもよい。プログラムすることが可能なＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）、または、ＬＳＩ内部の回路セルの接続および設定を再構成可能なリコンフィギュラブル・プロセッサを利用してもよい。
【０２５０】
さらには、半導体技術の進歩または派生する別技術によりＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて、散乱トモグラフィ装置に含まれる構成要素の集積回路化を行ってもよい。
【産業上の利用可能性】
【０２５１】
本発明に係る散乱トモグラフィ方法および散乱トモグラフィ装置は、曲率の大きい対象物に対する検査に有用であり、例えば、鉄道の在来線のトンネル内の欠陥の検査、道路のコンクリート覆工のレーダによる検査、コンクリート内部の鉄筋腐食の検査、被災地における鉄筋構造物の耐震検査、および、医療診断等に適用可能である。
【符号の説明】
【０２５２】
１０  アンテナ素子
２０  マルチパスアレイレーダ（散乱トモグラフィ装置）
３０  センサ部
３１  送信部
３２  受信部
３６  送信アンテナ素子
３７  受信アンテナ素子
４０  画像再構成部
５０  モニタ
【図１】
2014125815.tif 000077
【図２】
2014125815.tif 000078
【図３】
2014125815.tif 000079
【図４】
2014125815.tif 000080
【図５】
2014125815.tif 000081
【図６】
2014125815.tif 000082
【図７】
2014125815.tif 000083
【図８】
2014125815.tif 000084
【図９】
2014125815.tif 000085
【図１０Ａ】
2014125815.tif 000086
【図１０Ｂ】
2014125815.tif 000087
【図１１Ａ】
2014125815.tif 000088
【図１１Ｂ】
2014125815.tif 000089
【図１２Ａ】
2014125815.tif 000090
【図１２Ｂ】
2014125815.tif 000091
【図１３Ａ】
2014125815.tif 000092
【図１３Ｂ】
2014125815.tif 000093
【図１４】
2014125815.tif 000094
【図１５】
2014125815.tif 000095
【手続補正書】
【提出日】平成27年8月10日(2015.8.10)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００１０
【補正方法】変更
【補正の内容】
【００１０】
特開２００３−１７７６５６号公報
特開２００９−２８８１２９号公報
特開２００４−５１２１１７号公報
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００１６
【補正方法】変更
【補正の内容】
【００１６】
上記の課題を解決するため、本発明に係る散乱トモグラフィ方法は、物体に放射した波動の散乱波を解析する散乱トモグラフィ方法であって、前記波動を、曲面上に配置された複数の送信アンテナ素子から前記物体に放射するステップと、前記散乱波を、曲面上に配置された複数の受信アンテナ素子により受信するステップと、前記受信アンテナ素子により受信した散乱波を示す散乱波データを境界条件として用い、前記物体の内部情報に関する画像を再構成するステップとを含み、前記画像を再構成するステップにおいて、前記物体の内部情報に関する画像を再構成するための、後述する（式５）に定義される関数φをあらかじめ設定し、前記関数φが満たす、後述する（式１０）に定義される方程式を構築し、測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、後述する（式２４）に定義される映像化関数ρを導出し、前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００１７
【補正方法】変更
【補正の内容】
【００１７】
また、本発明に係る散乱トモグラフィ方法は、物体に放射した波動の散乱波を解析する散乱トモグラフィ方法であって、前記波動を、曲面上に配置された複数の送信アンテナ素子から前記物体に放射するステップと、前記散乱波を、曲面上に配置された複数の受信アンテナ素子により受信するステップと、前記受信アンテナ素子により受信した散乱波を示す散乱波データから、前記物体の内部情報に関する画像を再構成するステップとを含み、前記画像を再構成するステップにおいて、前記物体の内部情報に関する画像を再構成するための、後述する（式２７）に定義される関数φをあらかじめ設定し、前記関数φが満たす、後述する（式３９）に定義される方程式を構築し、測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、後述する（式５２）に定義される映像化関数ρを導出し、前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００２６
【補正方法】変更
【補正の内容】
【００２６】
上記の課題を解決するため、本発明に係る散乱トモグラフィ装置は、物体に放射した波動の散乱波を解析する散乱トモグラフィ装置であって、曲面上に配置され、物体に波動を放射する複数の送信アンテナ素子と、曲面上に配置され、前記放射された波動が前記物体において散乱した散乱波を受信する複数の受信アンテナ素子と、前記受信した散乱波を示す散乱波データから、前記物体の内部情報に関する画像を再構成する画像再構成部とを備え、前記画像再構成部は、前記物体の内部情報に関する画像を再構成するための、後述する（式５）に定義される関数φをあらかじめ設定し、前記関数φが満たす、後述する（式１０）に定義される方程式を構築し、測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、後述する（式２４）に定義される映像化関数ρを導出し、前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する。
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】００２８
【補正方法】変更
【補正の内容】
【００２８】
また、本発明に係る散乱トモグラフィ装置は、物体に放射した波動の散乱波を解析する散乱トモグラフィ装置であって、曲面上に配置され、物体に波動を放射する複数の送信アンテナ素子と、曲面上に配置され、前記放射された波動が前記物体において散乱した散乱波を受信する複数の受信アンテナ素子と、前記受信した散乱波を示す散乱波データから、前記物体の内部情報に関する画像を再構成する画像再構成部とを備え、前記画像再構成部は、前記物体の内部情報に関する画像を再構成するための、後述する（式２７）に定義される関数φをあらかじめ設定し、前記関数φが満たす、後述する（式３９）に定義される方程式を構築し、測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、後述する（式５２）に定義される映像化関数ρを導出し、前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する。
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】００７０
【補正方法】変更
【補正の内容】
【００７０】
画像再構成アルゴリズムについて数学的に簡単に説明すると、映像化に必要な基本解（関数）をはじめに設定し、その解から方程式を構築し、送信データおよび受信データから、より一般的な解を求めるという、逆問題を解くものである。
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】００７１
【補正方法】変更
【補正の内容】
【００７１】
すなわち、はじめに映像化に必要なグリーン関数を設定する。そして、この関数が厳密な解（関数）となる６次元空間（ｔ,ｘ,ｙ１,ｙ２,ｚ１,ｚ２）に関する偏微分方程式を構築する。そして、曲線状に配置された送信アンテナ素子３６から放射された送信データと曲線状に配置された受信アンテナ素子３７で受信された受信データとを境界条件としてこの偏微分方程式を解き、ｔ→０、ｘ→ｘ、ｚ２→ｚ１(＝ｚ)、ｙ１→ｙ２(＝ｙ)とおくことにより映像化関数を求める。これにより、映像化関数が求まり、コンクリート内部の欠陥等、物体内部の情報を汎用的かつ高速に映像化することができる。
【手続補正８】
【補正対象書類名】明細書
【補正対象項目名】００７８
【補正方法】変更
【補正の内容】
【００７８】
３次元空間の部分領域をＤとし、その境界を∂Ｄとする。このとき、関数Ｇ(ｒ１,ｒ２,ω)は、領域Ｄの内部で次の（式２）のような微分方程式の解となっている。この方程式を散乱場方程式と呼ぶ。
【手続補正９】
【補正対象書類名】明細書
【補正対象項目名】００８８
【補正方法】変更
【補正の内容】
【００８８】
次に、次の（式５）のような関数φを導入する。ここで、ｃは伝播速度、ｋは波数であり、波長をλとすると、ω＝ｃｋ、ｋ＝２π／λの関係がある。なお、（式５）に示すφは、本発明における物体の内部情報に関する画像を再構成するための関数（解）、すなわち、本発明における（式Ａ）に定義される関数に相当する。ξ、η、ζは、図６から図８に示す点Ｐ（ξ，η，ζ）のｘ座標、ｙ座標およびｚ座標である。なお、点Ｐ（ξ，η，ζ）は、領域中の任意の散乱点である。
【手続補正１０】
【補正対象書類名】明細書
【補正対象項目名】００９２
【補正方法】変更
【補正の内容】
【００９２】
次に、（式５）が漸近的な解（関数）となる偏微分方程式を求める。そのために、微分した結果生ずる１／ρに関して高次の項を無視して計算する。以下、微分の略記法を以下の（式７）のように定義する。
【手続補正１１】
【補正対象書類名】明細書
【補正対象項目名】００９４
【補正方法】変更
【補正の内容】
【００９４】
このとき、φは次の（式８）に示す方程式を満たす。この（式８）が、本発明における散乱場方程式に相当する。
【手続補正１２】
【補正対象書類名】明細書
【補正対象項目名】００９８
【補正方法】変更
【補正の内容】
【００９８】
これにより、最終的に次の（式１０）が得られる。この（式１０）が、本発明における（式Ｂ）に定義される方程式に相当する。
【手続補正１３】
【補正対象書類名】明細書
【補正対象項目名】０１０５
【補正方法】変更
【補正の内容】
【０１０５】
この方程式を解くのであるが、ｚ１,ｚ２の２変数がある。そのため、（式１２）を解くに当たり、固定された(ｘ,ｙ１,ｙ２)あるいは(ｋｘ,ｋｙ１,ｋｙ２)に関して、(ｚ１,ｚ２)空間の中の一次元の領域で与えられる境界条件が必要となる。
【手続補正１４】
【補正対象書類名】明細書
【補正対象項目名】０１２８
【補正方法】変更
【補正の内容】
【０１２８】
この式からａ(ｋ)を求めることができれば、（式８）の方程式の解は（式１８）から次の（式２３）のように表される。
【手続補正１５】
【補正対象書類名】明細書
【補正対象項目名】０１３３
【補正方法】変更
【補正の内容】
【０１３３】
３．積分方程式（式２２）の解
次に、上記した（式２２）の解を求める。すなわち、図８に示すように、断面曲線Ｓ上の点ＰＩ,ＰＪで測定されたデータからｘｙ面（平面ｚ＝０）上の点での近似データを求める。
【手続補正１６】
【補正対象書類名】明細書
【補正対象項目名】０１３４
【補正方法】変更
【補正の内容】
【０１３４】
（式２２）におけるφについて、曲線上の点ＰＩ,ＰＪで測定されたデータφ(ｘ,ｙＩ,ｙＪ,ｚＩ,ｚＪ,ｔ)をフーリエ変換したΦ(ｋｘ,ｙＩ,ｙＪ,ｋ)を、次の（式２５）のようにおく。
【手続補正１７】
【補正対象書類名】明細書
【補正対象項目名】０１７８
【補正方法】変更
【補正の内容】
【０１７８】
画像再構成アルゴリズムについて数学的に簡単に説明すると、実施の形態１に係る散乱トモグラフィ方法と同様、映像化に必要な基本解（関数）をはじめに設定し、その解から方程式を構築し、送信データおよび受信データから、より一般的な解を求めるという、逆問題を解くものである。
【手続補正１８】
【補正対象書類名】明細書
【補正対象項目名】０１７９
【補正方法】変更
【補正の内容】
【０１７９】
すなわち、はじめに映像化に必要なグリーン関数を設定する。そして、このグリーン関数の解となる関数を導入する。ここで、次に、この関数が厳密な解（関数）となる７次元空間（ｔ,ｘ１,ｘ２,ｙ１,ｙ２,ｚ１,ｚ２）に関する偏微分方程式を構築する。そして、曲面状に配置された送信アンテナ素子から放射された送信データと曲面状に配置された受信アンテナ素子で受信された受信データとを境界条件としてこの偏微分方程式を解き、ｔ→０、ｘ１→ｘ２(＝ｘ)、ｚ２→ｚ１(＝ｚ)、ｙ１→ｙ２(＝ｙ)とおくことに映像化関数を求める。これにより、映像化関数が求まり、トンネル内部の欠陥等、物体内部の情報を汎用的かつ高速に映像化することができる。
【手続補正１９】
【補正対象書類名】明細書
【補正対象項目名】０１８４
【補正方法】変更
【補正の内容】
【０１８４】
２．任意曲面における散乱の逆問題L(∂/∂t,∂/∂r１,∂/∂r２)の導出
以下、この微分作用素を求める方法について述べる。本実施の形態に係る解析モデルでは、実施の形態１の解析モデルとは異なり、送信点および受信点がある曲面上を動くということである。曲面上ではｒ１、ｒ２のｘ,ｙ,ｚ座標が必ずしも等しくない。具体的にはｒ１＝(ｘ１,ｙ１,ｚ１)、ｒ２＝(ｘ２,ｙ２,ｚ２)である。関数Ｇは実施の形態１に示した（式４）と同様に定義され、ｒ１＝(ｘ１,ｙ１,ｚ１)、ｒ２＝(ｘ２,ｙ２,ｚ２)から関数Ｇ(ｒ１,ｒ２,ω)の満たす方程式を求めると、以下の（式３７）のようになる。なお、（式３７）に示すφは、本発明における物体の内部情報に関する画像を再構成するための関数（解）、すなわち、本発明における（式Ｄ）に定義される関数に相当する。
【手続補正２０】
【補正対象書類名】明細書
【補正対象項目名】０１８６
【補正方法】変更
【補正の内容】
【０１８６】
ここでは、時間の因子をｅｘｐ(−ｉωｔ)に比例すると仮定し、波数をｋとしている。ｘ１、ｘ２、ｙ１、ｙ２、ｚ１、ｚ２は任意曲面上の座標である。ｃは伝播速度であり、波長をλとすると、ω＝ｃｋ、ｋ＝２π／λの関係がある。ξ、η、ζは、任意の点の座標である。
【手続補正２１】
【補正対象書類名】明細書
【補正対象項目名】０１８８
【補正方法】変更
【補正の内容】
【０１８８】
ここで、（式３７）が漸近的な解（関数）となる偏微分方程式を求める。そのためには、微分した結果生ずる１／ρに関して高次の項を無視して計算する。以下、微分の略記法を以下の（式３８）ように定義する。
【手続補正２２】
【補正対象書類名】明細書
【補正対象項目名】０１９０
【補正方法】変更
【補正の内容】
【０１９０】
このとき、実施の形態１と同様の計算により、φは次の（式３９）に示す方程式を満たす。この（式３９）が、本発明における散乱場方程式に相当する。また、この（式３９）が、本発明における（式Ｅ）に定義される関数に相当する。
【手続補正２３】
【補正対象書類名】明細書
【補正対象項目名】０１９７
【補正方法】変更
【補正の内容】
【０１９７】
この方程式を解くのであるが、ｚ１,ｚ２の２変数がある。そのため、（式４１）を解くに当たり、固定された(ｘ１,ｘ２,ｙ１,ｙ２)あるいは(ｋｘ１,ｋｘ２,ｋｙ１,ｋｙ２)に関して、(ｚ１,ｚ２)空間の中の一次元の領域で与えられる境界条件が必要となる。
【手続補正２４】
【補正対象書類名】明細書
【補正対象項目名】０２１８
【補正方法】変更
【補正の内容】
【０２１８】
この式からａ(ｋ)を求めることができれば、（式３９）の方程式の解は（式４６）から次の（式５１）のように表される。
【手続補正２５】
【補正対象書類名】明細書
【補正対象項目名】０２３９
【補正方法】変更
【補正の内容】
【０２３９】
この映像化関数により、物体の内部情報に関する画像を再構成する。したがって、物体内部の欠陥等、物体の内部情報に関する画像を汎用的かつ高速に映像化することができる。また、画像を再構成するステップにおいて、関数ρ（ｒ）の設定は、３次元空間を対象に行われるので、曲率の大きい曲面を有する対象物の内部情報をより精度よく高速に映像化することができる。
【手続補正２６】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
物体に放射した波動の散乱波を解析する散乱トモグラフィ方法であって、
前記波動を、曲面上に配置された複数の送信アンテナ素子から前記物体に放射するステップと、
前記散乱波を、曲面上に配置された複数の受信アンテナ素子により受信するステップと、
前記受信アンテナ素子により受信した散乱波を示す散乱波データを境界条件として用い、前記物体の内部情報に関する画像を再構成するステップとを含み、
前記画像を再構成するステップにおいて、
前記物体の内部情報に関する画像を再構成するための、（式Ａ）に定義される関数φをあらかじめ設定し、
前記関数φが満たす、（式Ｂ）に定義される方程式を構築し、
測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、（式Ｃ）に定義される映像化関数ρを導出し、
前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する
散乱トモグラフィ方法。
【数１】
2014125815.tif 000096
【数２】
2014125815.tif 000097
【数３】
2014125815.tif 000098
（式Ａ）から（式Ｃ）中において、ｘ、ｙ１、ｙ２、ｚ１、ｚ２は任意曲線上の座標、ξ、η、ζは領域中の任意の散乱点のｘ座標、ｙ座標、ｚ座標、ｃは伝播速度、ｋは波数であり、ω＝ｃｋである。
【請求項２】
物体に放射した波動の散乱波を解析する散乱トモグラフィ方法であって、
前記波動を、曲面上に配置された複数の送信アンテナ素子から前記物体に放射するステップと、
前記散乱波を、曲面上に配置された複数の受信アンテナ素子により受信するステップと、
前記受信アンテナ素子により受信した散乱波を示す散乱波データから、前記物体の内部情報に関する画像を再構成するステップとを含み、
前記画像を再構成するステップにおいて、
前記物体の内部情報に関する画像を再構成するための、（式Ｄ）に定義される関数φをあらかじめ設定し、
前記関数φが満たす、（式Ｅ）に定義される方程式を構築し、
測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、（式Ｆ）に定義される映像化関数ρを導出し、
前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する
散乱トモグラフィ方法。
【数４】
2014125815.tif 000099
【数５】
2014125815.tif 000100
【数６】
2014125815.tif 000101
（式Ｄ）から（式Ｆ）中において、ｘ１、ｘ２、ｙ１、ｙ２、ｚ１、ｚ２は任意曲面上の座標、ξ、η、ζは領域中の任意の散乱点のｘ座標、ｙ座標、ｚ座標、ｃは伝播速度、ｋは波数であり、ω＝ｃｋである。
【請求項３】
前記映像化関数は、高速フーリエ変換を用いて導出される
請求項１または２に記載の散乱トモグラフィ方法。
【請求項４】
前記波動は、電磁波または超音波である
請求項１〜３のいずれか１項に記載の散乱トモグラフィ方法。
【請求項５】
前記波動は、パルス波または所定の周波数を有する周期波である
請求項１〜４のいずれか１項に記載の散乱トモグラフィ方法。
【請求項６】
物体に放射した波動の散乱波を解析する散乱トモグラフィ装置であって、
曲面上に配置され、物体に波動を放射する複数の送信アンテナ素子と、
曲面上に配置され、前記放射された波動が前記物体において散乱した散乱波を受信する複数の受信アンテナ素子と、
前記受信した散乱波を示す散乱波データから、前記物体の内部情報に関する画像を再構成する画像再構成部とを備え、
前記画像再構成部は、
前記物体の内部情報に関する画像を再構成するための、（式Ａ）に定義される関数φをあらかじめ設定し、
前記関数φが満たす、（式Ｂ）に定義される方程式を構築し、
測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、（式Ｃ）に定義される映像化関数ρを導出し、
前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する
散乱トモグラフィ装置。
【数７】
2014125815.tif 000102
【数８】
2014125815.tif 000103
【数９】
2014125815.tif 000104
（式Ａ）から（式Ｃ）中において、ｘ、ｙ１、ｙ２、ｚ１、ｚ２は任意曲線上の座標、ξ、η、ζは領域中の任意の散乱点のｘ座標、ｙ座標、ｚ座標、ｃは伝播速度、ｋは波数であり、ω＝ｃｋである。
【請求項７】
物体に放射した波動の散乱波を解析する散乱トモグラフィ装置であって、
曲面上に配置され、物体に波動を放射する複数の送信アンテナ素子と、
曲面上に配置され、前記放射された波動が前記物体において散乱した散乱波を受信する複数の受信アンテナ素子と、
前記受信した散乱波を示す散乱波データから、前記物体の内部情報に関する画像を再構成する画像再構成部とを備え、
前記画像再構成部は、
前記物体の内部情報に関する画像を再構成するための、（式Ｄ）に定義される関数φをあらかじめ設定し、
前記関数φが満たす、（式Ｅ）に定義される方程式を構築し、
測定により得られた前記散乱波データから、前記方程式を解くことにより得られる、（式Ｆ）に定義される映像化関数ρを導出し、
前記映像化関数ρにより、前記物体の内部情報に関する画像を再構成する
散乱トモグラフィ装置。
【数１０】
2014125815.tif 000105
【数１１】
2014125815.tif 000106
【数１２】
2014125815.tif 000107
（式Ｄ）から（式Ｆ）中において、ｘ１、ｘ２、ｙ１、ｙ２、ｚ１、ｚ２は任意曲面上の座標、ξ、η、ζは領域中の任意の散乱点のｘ座標、ｙ座標、ｚ座標、ｃは伝播速度、ｋは波数であり、ω＝ｃｋである。
【請求項８】
前記画像再構成部は、高速フーリエ変換を用いて前記映像化関数を導出する
請求項６または７に記載の散乱トモグラフィ装置。
【請求項９】
前記波動は、電磁波または超音波である
請求項６〜８のいずれか１項に記載の散乱トモグラフィ装置。
【請求項１０】
前記波動は、パルス波または所定の周波数を有する周期波である
請求項６〜９のいずれか１項に記載の散乱トモグラフィ装置。
【国際調査報告】
2014125815.tif 000108
2014125815.tif 000109
2014125815.tif 000110
2014125815.tif 000111