(57)【要約】  画像化方法は、波動を対象物（１０）に放射するステップと、波動が対象物（１０）において散乱した散乱波を受信するステップと、散乱波を示す散乱波データから、対象物（１０）の内部情報に関する画像を再構成するステップとを含む。画像を再構成するステップにおいては、散乱波データと形状を示す解析モデルとを用いて偏微分方程式を解くことにより再構成関数を導出し、再構成関数により、対象物（１０）の内部情報に関する画像を再構成する。ここで偏微分方程式は、対象物（１０）の内部情報に関する画像を再構成するための再構成関数が満たす方程式である。
2017057524.tif 000002
【特許請求の範囲】
【請求項１】
波動を対象物に放射するステップと、
前記波動が前記対象物において散乱した散乱波を受信するステップと、
前記散乱波を示す散乱波データから、前記対象物の内部情報に関する画像を再構成するステップとを含み、
前記画像を再構成するステップにおいて、
前記散乱波データと形状を示す解析モデルとを用いて偏微分方程式を解くことにより、前記対象物の内部情報に関する画像を再構成するための再構成関数を導出し、
前記再構成関数を用いて、前記対象物の内部情報に関する画像を再構成し、
前記偏微分方程式は、前記再構成関数が満たす方程式である
画像化方法。
【請求項２】
前記偏微分方程式は、
前記対象物内の送信点および受信点の位置を示す独立変数を用いて表され、
前記独立変数の数と同一の次元を持つ空間における各点において、前記散乱波が生じている場の関数である散乱場関数が解となる線形偏微分方程式であり、
前記画像を再構成するステップでは、前記再構成関数の時間変数の極限値である映像化関数を導出し、前記映像化関数を用いて前記対象物の内部情報に関する画像を再構成する請求項１に記載の画像化方法。
【請求項３】
前記解析モデルが示す形状は円錐であり、
前記映像化関数は、以下の（式Ａ）により表され、
前記画像を再構成するステップでは、前記映像化関数をθで積分することにより得られる以下の（式Ｂ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成する請求項２に記載の画像化方法。
【数１】
2017057524.tif 000090
..............（式Ａ）
【数２】
2017057524.tif 000091
...........（式Ｂ）
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記円錐の頂点を原点とし、前記原点から前記円錐の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、φＲはφ＝φＲδ（ｘ）が成り立つ関数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記円錐の軸を回転軸とする回転角度である。
【請求項４】
前記解析モデルが示す形状は曲線状の母線を有する略円錐であり、
前記映像化関数は、以下の（式Ｃ）により表され、
前記画像を再構成するステップでは、前記映像化関数をθで積分することにより得られる以下の（式Ｄ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成する請求項２に記載の画像化方法。
【数３】
2017057524.tif 000092
..............（式Ｃ）
【数４】
2017057524.tif 000093
...........（式Ｄ）
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記略円錐の頂点を原点とし、前記原点から前記略円錐の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、ａθは係数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記略円錐の軸を回転軸とする回転角度である。
【請求項５】
前記解析モデルが示す形状は非対称の錐体であり、
前記映像化関数は、ξ、η、およびζがそれぞれ以下の（式Ｅ）を満たす変数であるとき以下の（式Ｆ）により表され、
前記画像を再構成するステップでは、前記映像化関数をθで積分することにより得られる以下の（式Ｇ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成する請求項２に記載の画像化方法。
【数５】
2017057524.tif 000094
..............（式Ｅ）
【数６】
2017057524.tif 000095
..............（式Ｆ）
【数７】
2017057524.tif 000096
...........（式Ｇ）
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記錐体の頂点を原点とし、前記原点から前記錐体の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、φＲはφ＝φＲδ（ｘ）が成り立つ関数であり、ａθは係数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記錐体の軸を回転軸とする回転角度であり、αは前記Ｚ方向を基準とした傾き角度である。
【請求項６】
前記放射するステップでは、前記散乱波を受信する受信アンテナ素子と前記対象物に前記波動を放射する送信アンテナ素子とが一体に設けられ、前記受信アンテナ素子と前記送信アンテナ素子との間に電波吸収手段が位置しているプローブを用いて前記波動を前記対象物に放射し、
前記受信するステップでは、前記プローブを用いて前記散乱波を受信する請求項１または２に記載の画像化方法。
【請求項７】
前記解析モデルが示す形状は円錐であり、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち円錐の母線に相当する線に沿って移動される
請求項６に記載の画像化方法。
【請求項８】
前記解析モデルが示す形状は、曲線状の母線を有する略円錐であり、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち略円錐の曲線状の母線に相当する線に沿って移動される
請求項６に記載の画像化方法。
【請求項９】
前記解析モデルが示す形状は非対称の錐体であり、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち非対称の錐体の母線に相当する線に沿って移動される
請求項６に記載の画像化方法。
【請求項１０】
前記波動は、マイクロ波である
請求項１〜９のいずれか１項に記載の画像化方法。
【請求項１１】
前記波動は、パルス波または所定の周波数を有する周期波である
請求項１〜１０のいずれか１項に記載の画像化方法。
【請求項１２】
波動を対象物に放射する送信アンテナ素子と、
前記送信アンテナ素子から放射された前記波動が前記対象物において散乱した散乱波を受信する受信アンテナ素子と、
前記受信アンテナ素子により受信した散乱波を示す散乱波データから、前記対象物の内部情報に関する画像を再構成する画像再構成部とを備え、
前記画像再構成部は、
前記散乱波データと形状を示す解析モデルとを用いて偏微分方程式を解くことにより、前記対象物の内部情報に関する画像を再構成するための再構成関数を導出し、
前記再構成関数を用いて、前記対象物の内部情報に関する画像を再構成し、
前記偏微分方程式は、前記再構成関数が満たす方程式である
画像化装置。
【請求項１３】
前記偏微分方程式は、
前記対象物内の送信点および受信点の位置を示す独立変数を用いて表され、
前記独立変数の数と同一の次元を持つ空間における各点において、前記散乱波が生じている場の関数である散乱場関数が解となる線形偏微分方程式であり、
前記画像再構成部は、前記再構成関数の時間変数の極限値である映像化関数を導出し、前記映像化関数を用いて前記対象物の内部情報に関する画像を再構成する請求項１２に記載の画像化装置。
【請求項１４】
前記解析モデルが示す形状は円錐であり、
前記映像化関数は、以下の（式Ａ）により表され、
前記画像再構成部は、前記映像化関数をθで積分することにより得られる以下の（式Ｂ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成する請求項１３に記載の画像化装置。
【数８】
2017057524.tif 000097
..............（式Ａ）
【数９】
2017057524.tif 000098
...........（式Ｂ）
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記円錐の頂点を原点とし、前記原点から前記円錐の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、φＲはφ＝φＲδ（ｘ）が成り立つ関数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記円錐の軸を回転軸とする回転角度である。
【請求項１５】
前記解析モデルが示す形状は曲線状の母線を有する略円錐であり、
前記映像化関数は、以下の（式Ｃ）により表され、
前記画像再構成部は、前記映像化関数をθで積分することにより得られる以下の（式Ｄ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成する請求項１３に記載の画像化装置。
【数１０】
2017057524.tif 000099
..............（式Ｃ）
【数１１】
2017057524.tif 000100
...........（式Ｄ）
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記略円錐の頂点を原点とし、前記原点から前記略円錐の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、ａθは係数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記略円錐の軸を回転軸とする回転角度である。
【請求項１６】
前記解析モデルが示す形状は非対称の錐体であり、
前記映像化関数は、ξ、η、およびζがそれぞれ以下の（式Ｅ）を満たす変数であるとき以下の（式Ｆ）により表され、
前記画像再構成部は、前記映像化関数をθで積分することにより得られる以下の（式Ｇ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成する請求項１３に記載の画像化装置。
【数１２】
2017057524.tif 000101
..............（式Ｅ）
【数１３】
2017057524.tif 000102
..............（式Ｆ）
【数１４】
2017057524.tif 000103
...........（式Ｇ）
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記錐体の頂点を原点とし、前記原点から前記錐体の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、φＲはφ＝φＲδ（ｘ）が成り立つ関数であり、ａθは係数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記錐体の軸を回転軸とする回転角度であり、αは前記Ｚ方向を基準とした傾き角度である。
【請求項１７】
前記受信アンテナ素子と前記送信アンテナ素子とが一体に設けられたプローブをさらに備え、
前記受信アンテナ素子と前記送信アンテナ素子との間には電波吸収手段が位置する請求項１２または１３に記載の画像化装置。
【請求項１８】
前記解析モデルが示す形状は円錐であり、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち円錐の母線に相当する線に沿って移動される
請求項１７に記載の画像化装置。
【請求項１９】
前記解析モデルが示す形状は、曲線状の母線を有する略円錐であり、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち略円錐の曲線状の母線に相当する線に沿って移動される
請求項１７に記載の画像化装置。
【請求項２０】
前記解析モデルが示す形状は非対称の錐体であり、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち非対称の錐体の母線に相当する線に沿って移動される
請求項１７に記載の画像化装置。
【請求項２１】
前記波動は、マイクロ波である
請求項１２〜２０のいずれか１項に記載の画像化装置。
【請求項２２】
前記波動は、パルス波または所定の周波数を有する周期波である
請求項１２〜２１のいずれか１項に記載の画像化装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、波動を利用して物体の内部情報を取得し映像化（画像化）する技術に関し、特に、波動の散乱現象を利用して物体の内部情報を取得する画像化方法および画像化装置に関する。
【背景技術】
【０００２】
従来、生体や建築物等の物体内部の情報を映像化する方法として、Ｘ線ＣＴ（Ｘ線トモグラフィ）、ＭＲＩ（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ：核磁気共鳴画像法）、ＰＥＴ（Ｐｏｓｉｔｒｏｎ  Ｅｍｉｓｓｉｏｎ  Ｔｏｍｏｇｒａｐｈｙ：ポジトロン断層法）などの方法が用いられている。具体的には、光、テラヘルツ、ミリ波、マイクロ波のような電磁波あるいは超音波などの波動を観測対象物である生体や物体、あるいは、プラズマに放射し、その散乱波（反射波）を観測し解析することで、生体内部や固体内部、あるいは、プラズマの内部の情報を映像化している。また、最近では、波動に代えて磁場を用いて、生体や物体内部の情報の映像化も行われている。【０００３】
一般的に、これらの方法では、物体Ｏに電磁波や超音波のような波動ｕを放射し、物体Ｏの回りの多数の場所で物体Ｏから散乱される散乱波ｐを観測し、得られたデータを映像化するという技術が採用されている（例えば、特許文献１、非特許文献１参照）。【０００４】
特許文献１に記載の技術は、電波を用いて物体内部の情報の映像化を行うものである。円周上に配置されたセンサ素子で観測される散乱波のデータを、導電率や誘電率等のパラメータで修正しながら、データの取得を繰り返し、映像化するものである。【０００５】
非特許文献１に記載の技術は、マルチパスリニアアレイレーダに関する技術であり、コンクリート内部の欠陥等の情報の映像化を行うものである。対象物の表面に直線状または曲線状に配置されたセンサ素子（直線状のマルチアレイアンテナ）を配置し、放射された波動の散乱波をセンサで観測し、観測データを解析し映像化するものである。【０００６】
また、医療の分野においては、例えば、乳がんの検出を行うためのマンモグラフィ装置等に、このような波動を使った観測方法を利用することが想定される。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００３−１７７６５６号公報
【非特許文献】
【０００８】
【非特許文献１】トンネル覆工コンクリート検査用３次元映像化レーダを開発、三井造船技報、Ｎｏ．１８４、ｐ２４、２００５年２月
【発明の概要】
【発明が解決しようとする課題】
【０００９】
特許文献１および非特許文献１に開示された技術では、例えば、物体の曲面形状等の条件が変わるたびに、理論や装置内部の構造等を変えてデータを再取得したり、取得したデータを補正したりする必要があるため、汎用的な使用が難しい。特に、生体等のフレキシブルな形状を有する対象物については、対象物の形状や外形の曲面の曲率が一定でない場合が多いため、従来の直線状マルチアレイアンテナや一定の形状の対象物をモデルとしたアレイアンテナを一律に用いることは困難である。また、データを再取得したり補正したりする必要があるため、計算速度の遅延やメモリの大量使用という問題がある。【００１０】
また、対象物を生体とした場合の従来技術として、マンモグラフィ装置などＸ線を利用した装置がある。しかし、Ｘ線を使用したマンモグラフィ装置では、特にアジア人の大半を占める高濃度乳房に対しては、炭酸カルシウムを撮像することによって間接的にがんの診断を行うことがあるため、直接微細が乳がん組織を他の組織と識別することは困難である。【００１１】
さらに、近年では、高輝度Ｘ線の散乱を利用した手法が提案されているが、透過光に対して数百万分の１の強度の散乱光を測定するものであるため、正確な診断を行うにはＸ線の強度を高くする必要がある。この場合、生体に対するＸ線の影響が懸念される。【００１２】
そこで、本発明は、汎用的で、物体内部の情報を簡便かつ高速に映像化することができる検査方法および検査装置を提供することを目的とする。
【課題を解決するための手段】
【００１３】
上記の課題を解決するため、本発明に係る画像化方法は、
波動を対象物に放射するステップと、
前記波動が前記対象物において散乱した散乱波を受信するステップと、
前記散乱波を示す散乱波データから、前記対象物の内部情報に関する画像を再構成するステップとを含み、
前記画像を再構成するステップにおいて、
前記散乱波データと形状を示す解析モデルとを用いて偏微分方程式を解くことにより、前記対象物の内部情報に関する画像を再構成するための再構成関数を導出し、
前記再構成関数を用いて、前記対象物の内部情報に関する画像を再構成し、
前記偏微分方程式は、前記再構成関数が満たす方程式である。【００１４】
これにより、送信アンテナ素子と受信アンテナ素子とを有するセンサ素子を対象物上で自由に移動させながら対象物の内部情報の観測（検査）を行う解析モデルにおいて、逆問題用の偏微分方程式を設定し、これを解くことで、対象物の内部情報を汎用的かつ高速に映像化することができる。また、画像を再構成するステップにおいて、再構成関数φの設定は、３次元空間を対象に行われるので、対象物の内部情報を３次元で精度よく高速に映像化することができる。また、対象物を、他の物体についての観測結果を介して観測するのではなく、直接的に観測することができる。【００１５】
さらに、複数の送信部と受信部とを備えたマルチスタティック法による散乱トモグラフィを利用した検査装置と比較すると、この検査方法は、モノスタティック法による散乱トモグラフィを利用した検査方法であり、形状がフレキシブルな検査対象であっても精度よく計測することができる。【００１６】
また、前記偏微分方程式は、
前記対象物内の送信点および受信点の位置を示す独立変数を用いて表され、
前記独立変数の数と同一の次元を持つ空間における各点において、前記散乱波が生じている場の関数である散乱場関数が解となる線形偏微分方程式であり、
前記画像を再構成するステップでは、前記再構成関数の時間変数の極限値である映像化関数を導出し、前記映像化関数を用いて前記対象物の内部情報に関する画像を再構成してもよい。【００１７】
これにより、対象物の内部情報をより精度よく高速に映像化することができる。【００１８】
また、前記解析モデルが示す形状は円錐であり、
前記映像化関数は、以下の（式Ａ）により表され、
前記画像を再構成するステップでは、前記映像化関数をθで積分することにより得られる以下の（式Ｂ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成してもよい。
【数１】
2017057524.tif 000003
..............（式Ａ）
【数２】
2017057524.tif 000004
...........（式Ｂ）
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記円錐の頂点を原点とし、前記原点から前記円錐の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、φＲはφ＝φＲδ（ｘ）が成り立つ関数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記円錐の軸を回転軸とする回転角度である。【００１９】
これにより、およそ円錐状の形状を有する対象物の内部情報をより精度よく高速に映像化することができる。【００２０】
また、前記解析モデルが示す形状は曲線状の母線を有する略円錐であり、
前記映像化関数は、以下の（式Ｃ）により表され、
前記画像を再構成するステップでは、前記映像化関数をθで積分することにより得られる以下の（式Ｄ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成してもよい。
【数３】
2017057524.tif 000005
..............（式Ｃ）
【数４】
2017057524.tif 000006
...........（式Ｄ）
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記略円錐の頂点を原点とし、前記原点から前記略円錐の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、ａθは係数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記略円錐の軸を回転軸とする回転角度である。【００２１】
これにより、およそ略円錐の形状を有する対象物の内部情報をより精度よく高速に映像化することができる。【００２２】
また、前記解析モデルが示す形状は非対称の錐体であり、
前記映像化関数は、ξ、η、およびζがそれぞれ以下の（式Ｅ）を満たす変数であるとき以下の（式Ｆ）により表され、
前記画像を再構成するステップでは、前記映像化関数をθで積分することにより得られる以下の（式Ｇ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成してもよい。
【数５】
2017057524.tif 000007
..............（式Ｅ）
【数６】
2017057524.tif 000008
..............（式Ｆ）
【数７】
2017057524.tif 000009
...........（式Ｇ）
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記錐体の頂点を原点とし、前記原点から前記錐体の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、φＲはφ＝φＲδ（ｘ）が成り立つ関数であり、ａθは係数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記錐体の軸を回転軸とする回転角度であり、αは前記Ｚ方向を基準とした傾き角度である。【００２３】
これにより、およそ非対称の錐体の形状を有する対象物の内部情報をより精度よく高速に映像化することができる。【００２４】
また、前記放射するステップでは、前記散乱波を受信する受信アンテナ素子と前記対象物に前記波動を放射する送信アンテナ素子とが一体に設けられ、前記受信アンテナ素子と前記送信アンテナ素子との間に電波吸収手段が位置しているプローブを用いて前記波動を前記対象物に放射し、
前記受信するステップでは、前記プローブを用いて前記散乱波を受信してもよい。【００２５】
これにより、送信アンテナ素子から発信された波動のうち、直接受信アンテナ素子に受信される波動の大きさを小さくすることができ、その直接受信される波動の前記大きさに対する、受信アンテナ素子に受信される対象物からの散乱波の大きさの比を高めることができる。したがって、より鮮明な対象物内の画像を得ることができる。【００２６】
また、前記解析モデルが示す形状は円錐であり、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち円錐の母線に相当する線に沿って移動されてもよい。【００２７】
これにより、円錐の母線上をセンサ素子が移動する解析モデルにおいて逆問題用の偏微分方程式を設定し、これを解くことで、対象物の内部情報を汎用的かつ高速に映像化することができる。特に、円錐状の形状を有する対象物に対して、内部情報をより精度よく高速に映像化することができる。【００２８】
また、解析モデルが示す形状は、曲線状の母線を有する略円錐であり、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち略円錐の曲線状の母線に相当する線に沿って移動されてもよい。【００２９】
これにより、曲線状の母線上をセンサ素子が移動する解析モデルにおいて逆問題用の偏微分方程式を設定し、これを解くことで、対象物の内部情報を汎用的かつ高速に映像化することができる。特に、略円錐状、例えば、半球またはドーム状の形状を有する対象物に対して、内部情報をより精度よく高速に映像化することができる。【００３０】
また、前記解析モデルが示す形状は非対称の錐体であり、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち非対称の錐体の母線に相当する線に沿って移動されてもよい。【００３１】
これにより、非対称の錐体の形状の母線上をセンサ素子が移動する解析モデルにおいて逆問題用の偏微分方程式を設定し、これを解くことで、対象物の内部情報を汎用的かつ高速に映像化することができる。特に、円錐状または略円錐状の一部が欠けた形状を有する対象物やフレキシブルな形状を有する対象物に対して、内部情報をより精度よく高速に映像化することができる。【００３２】
また、前記波動は、マイクロ波であってもよい。【００３３】
これにより、水分が多いか少ないかにかかわらず、汎用的かつ簡便な方法で、対象物の内部情報を映像化することができる。【００３４】
また、前記波動は、パルス波または所定の周波数を有する周期波であってもよい。【００３５】
これにより、汎用的かつ簡便な方法で、フレキシブルな形状を有する対象物の内部情報を映像化することができる。【００３６】
上記の課題を解決するため、本発明に係る画像化装置は、波動を対象物に放射する送信アンテナ素子と、
前記送信アンテナ素子から放射された前記波動が前記対象物において散乱した散乱波を受信する受信アンテナ素子と、
前記受信アンテナ素子により受信した散乱波を示す散乱波データから、前記対象物の内部情報に関する画像を再構成する画像再構成部とを備え、
前記画像再構成部は、
前記散乱波データと形状を示す解析モデルとを用いて偏微分方程式を解くことにより、前記対象物の内部情報に関する画像を再構成するための再構成関数を導出し、
前記再構成関数により、前記対象物の内部情報に関する画像を再構成し、
前記偏微分方程式は、前記再構成関数が満たす方程式である。【００３７】
これにより、送信アンテナ素子と受信アンテナ素子とを有するセンサ素子を対象物上で自由に移動させながら対象物の観測（検査）を行う解析モデルにおいて、逆問題用の偏微分方程式を設定し、これを解くことで、対象物の内部情報を汎用的かつ高速に映像化することができる。また、画像を再構成する際に、再構成関数φの設定は３次元空間を対象に行われるので、対象物の内部情報を３次元で精度よく高速に映像化することができる。また、対象物を、他の物体の観測結果を介して観測するのではなく、直接的に観測することができる。【００３８】
さらに、複数の送信部と受信部とを備えたマルチスタティック法による散乱トモグラフィを利用した検査装置と比較すると、この検査装置は、モノスタティック法による散乱トモグラフィを利用した検査装置であり、形状がフレキシブルな検査対象であっても精度よく計測することができる。【００３９】
また、前記偏微分方程式は、
前記対象物内の送信点および受信点の位置を示す独立変数を用いて表され、
前記独立変数の数と同一の次元を持つ空間における各点において、前記散乱波が生じている場の関数である散乱場関数が解となる線形偏微分方程式であり、
前記画像再構成部は、前記再構成関数の時間変数の極限値である映像化関数を導出し、前記映像化関数を用いて前記対象物の内部情報に関する画像を再構成してもよい。【００４０】
これにより、対象物の内部情報をより精度よく高速に映像化することができる。【００４１】
また、前記解析モデルが示す形状は円錐であり、
前記映像化関数は、以下の（式Ａ）により表され、
前記画像再構成部は、前記映像化関数をθで積分することにより得られる以下の（式Ｂ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成してもよい。
【数８】
2017057524.tif 000010
..............（式Ａ）
【数９】
2017057524.tif 000011
...........（式Ｂ）
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記円錐の頂点を原点とし、前記原点から前記円錐の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、φＲはφ＝φＲδ（ｘ）が成り立つ関数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記円錐の軸を回転軸とする回転角度である。【００４２】
これにより、およそ円錐状の形状を有する対象物の内部情報をより精度よく高速に映像化することができる。【００４３】
また、前記解析モデルが示す形状は曲線状の母線を有する略円錐であり、
前記映像化関数は、以下の（式Ｃ）により表され、
前記画像再構成部は、前記映像化関数をθで積分することにより得られる以下の（式Ｄ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成してもよい。
【数１０】
2017057524.tif 000012
..............（式Ｃ）
【数１１】
2017057524.tif 000013
...........（式Ｄ）
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記略円錐の頂点を原点とし、前記原点から前記略円錐の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、ａθは係数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記略円錐の軸を回転軸とする回転角度である。【００４４】
これにより、およそ略円錐の形状を有する対象物の内部情報をより精度よく高速に映像化することができる。【００４５】
また、前記解析モデルが示す形状は非対称の錐体であり、
前記映像化関数は、ξ、η、およびζがそれぞれ以下の（式Ｅ）を満たす変数であるとき以下の（式Ｆ）により表され、
前記画像再構成部は、前記映像化関数をθで積分することにより得られる以下の（式Ｇ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成してもよい。
【数１２】
2017057524.tif 000014
..............（式Ｅ）
【数１３】
2017057524.tif 000015
..............（式Ｆ）
【数１４】
2017057524.tif 000016
...........（式Ｇ）
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記錐体の頂点を原点とし、前記原点から前記錐体の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、φＲはφ＝φＲδ（ｘ）が成り立つ関数であり、ａθは係数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記錐体の軸を回転軸とする回転角度であり、αは前記Ｚ方向を基準とした傾き角度である。【００４６】
これにより、およそ非対称の錐体の形状を有する対象物の内部情報をより精度よく高速に映像化することができる。【００４７】
また、画像化装置は、前記受信アンテナ素子と前記送信アンテナ素子とが一体に設けられたプローブをさらに備え、
前記受信アンテナ素子と前記送信アンテナ素子との間には電波吸収手段が位置してもよい。【００４８】
これにより、送信アンテナ素子から発信された波動のうち、直接受信アンテナ素子に受信される波動の大きさを小さくすることができ、その直接受信される波動の前記大きさに対する、受信アンテナ素子に受信される対象物からの散乱波の大きさの比を高めることができる。したがって、より鮮明な対象物内の画像を得ることができる。【００４９】
また、前記解析モデルが示す形状は円錐であり、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち円錐の母線に相当する線に沿って移動されてもよい。【００５０】
これにより、円錐の母線上をセンサ素子が移動する解析モデルにおいて逆問題用の偏微分方程式を設定し、これを解くことで、対象物の内部情報を汎用的かつ高速に映像化することができる。特に、円錐状の形状を有する対象物に対して、内部情報をより精度よく高速に映像化することができる。【００５１】
また、前記解析モデルが示す形状は、曲線状の母線を有する略円錐であり、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち略円錐の曲線状の母線に相当する線に沿って移動されてもよい。【００５２】
これにより、曲線状の母線を有する略円錐の母線上をセンサ素子が移動する解析モデルにおいて逆問題用の偏微分方程式を設定し、これを解くことで、対象物の内部情報を汎用的かつ高速に映像化することができる。特に、略円錐状、例えば、半球またはドーム状の形状を有する対象物に対して、内部情報をより精度よく高速に映像化することができる。【００５３】
また、前記解析モデルが示す形状は非対称の錐体であり、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち非対称の錐体の母線に相当する線に沿って移動されてもよい。【００５４】
これにより、非対称の錐体の形状の母線上をセンサ素子が移動する解析モデルにおいて逆問題用の偏微分方程式を設定し、これを解くことで、対象物の内部情報を汎用的かつ高速に映像化することができる。特に、円錐状または略円錐状の一部が欠けた形状を有する対象物や様々な曲面形状を有する対象物に対して、内部情報をより精度よく高速に映像化することができる。【００５５】
また、前記波動は、マイクロ波であってもよい。【００５６】
これにより、対象物の水分が多いか少ないかにかかわらず、汎用的かつ簡便な方法で、対象物の内部情報を映像化することができる。【００５７】
また、前記波動は、パルス波または所定の周波数を有する周期波であってもよい。【００５８】
これにより、汎用的かつ簡便な方法で、様々な曲面形状を有する対象物の内部情報を映像化することができる。
【発明の効果】
【００５９】
本発明により、逆問題の解析を汎用的かつ高速に行い、様々な曲面形状を有する物体内部の情報を簡便に映像化することができる。様々な曲面形状を有する対象物の内部情報をより精度よく高速に映像化することができる。
【図面の簡単な説明】
【００６０】
上述した目的、およびその他の目的、特徴および利点は、以下に述べる好適な実施の形態、およびそれに付随する以下の図面によってさらに明らかになる。【００６１】
【図１】実施の形態１にかかるマンモグラフィ装置の構成を示す概略図
【図２】実施の形態１にかかるマンモグラフィ装置におけるプローブの構成を示す概略図
【図３】図２に示したプローブの使用方法の一例を示す図
【図４】実施の形態１にかかるマンモグラフィ装置の動作を示すフローチャート
【図５】実施の形態１にかかるマンモグラフィ方法の原理を説明するための解析モデル
【図６】実施の形態１にかかるマンモグラフィ方法の原理を説明するための解析モデル
【図７】実施の形態１にかかるマンモグラフィ方法の原理を説明するための解析モデル
【図８】実施の形態１にかかるマンモグラフィ装置を用いた計測実験のための実験モデルを示す概略図
【図９】図８に示した実験モデルをマンモグラフィ装置で計測した計測結果を示す図
【図１０】実施の形態２にかかるマンモグラフィ方法の原理を説明するための解析モデル
【図１１】実施の形態２にかかるマンモグラフィ方法の原理を説明するための解析モデル
【図１２】実施の形態３にかかるマンモグラフィ方法の原理を説明するための解析モデル
【発明を実施するための形態】
【００６２】
（本発明の基礎となった知見）
本発明に係る実施の形態について説明する前に、本発明の基礎となった技術について説明する。
【００６３】
本発明にかかる検査方法および検査装置は、散乱トモグラフィを利用した検査技術である。
【００６４】
散乱トモグラフィとは、波動を物体に放射したときの散乱波を観測し解析して物体内部の情報を映像化する技術である。すなわち、本発明にかかる検査方法および検査装置は、物体に波動を放射して生じた散乱波を解析することで、非破壊で対象物の欠陥等を検出する技術である。散乱トモグラフィの詳細な原理は、以下のとおりである。
【００６５】
一般に、物体に波動を放射した際に生じる散乱現象は、演算子を用いて表すことができる。例えば、物体Ｏ、放射波ｕ、観測データｐで表される物理方程式を、演算子Ａを用いてｐ＝Ａｕ［Ｏ］と表すことができる。ここで、物体Ｏ、放射波ｕ、演算子（システム関数）Ａが既知である場合に、観測データｐを求める問題は、順方向問題と呼ばれている。
【００６６】
これに対して、医療や産業上で重要な課題は、放射波ｕ、システム関数Ａ、観測データｐが既知である場合に物体Ｏが何であるかを求める問題である。この問題は、物理現象の因果関係を逆方向に辿るという意味で、逆方向問題と呼ばれ、Ｏ＝Ａｕ−１［ｐ］と表すことができる。散乱トモグラフィは、この逆方向問題を利用するものである。
【００６７】
散乱トモグラフィを利用した一例として、Ｍｕｌｔｉ−Ｐａｔｈ  Ｌｉｎｅａｒ  Ａｒｒａｙ  Ｒａｄａｒ（ＭＰＬＡ  Ｒａｄａｒ）がある。この方法は、例えば、センサ素子として対象物にアンテナ素子を取り付け、曲率のある対象物であっても平面として近似して、アンテナ素子から放射した電磁波と対象物において反射した反射波（散乱波）との関係から、非破壊で対象物の欠陥等を検出するものである。
【００６８】
ＭＰＬＡ  Ｒａｄａｒに使用されるアンテナ素子は、対象物に波動を放射する複数の送信アンテナ素子と、対象物で生じた散乱波を受信する複数の受信アンテナ素子とで構成されている。各アンテナ素子は、事前に想定された解析モデルにおいて、配置位置が決定されている。そして、ある位置に配置された送信アンテナ素子から放射された波動の散乱波を、他の位置に配置された受信アンテナ素子により受信する。これにより、ＭＰＬＡ  Ｒａｄａｒでは、受信アンテナ素子で受信された散乱波を用いてマルチスタティックな解析を行い、物体内部の情報を映像化する。
【００６９】
ここで、解析モデルとして、生体モデルを想定することで、生体の内部情報を映像化することができる。例えば、解析モデルとして乳房を想定すると、乳房の内部における癌組織の位置および大きさを検出するマンモグラフィとして利用することができる。
【００７０】
生体等を対象物としてその内部情報を映像化する場合、対象物の形状が、想定された解析モデルと適合しないことがある。特に、対象物がフレキシブルな形状を有する場合には、対象物の形状が配置位置や向きにより変わるため、想定された解析モデルにおける受信アンテナ素子では、散乱波が正常に受信できない等の問題が生じる。したがって、アンテナ素子の位置と受信信号の強度が正確に把握できず、解析を行った場合に内部情報が正確に把握できないことが考えられる。
【００７１】
そこで、本発明にかかる検査方法および検査装置は、一対の送信アンテナ素子と受信アンテナ素子とで構成されるアンテナ素子を、事前に想定された対象物の解析モデルにおいて、対象物上を自由に移動させながら解析を行うものである。すなわち、本発明における検査方法および検査装置では、対になった送信アンテナ素子と受信アンテナ素子のうちの送信アンテナ素子から対象物に波動を放射し、その散乱波を、波動を放射した送信アンテナと対になった受信アンテナ素子により受信する。これにより、当該マンモグラフィでは、受信アンテナ素子で受信された散乱波を用いてモノスタティックな解析を行い、物体内部の情報を映像化する。
【００７２】
以下、図面を用いて、本発明に係る実施の形態について説明する。以下の実施の形態では、乳房を検査の対象物として乳房の内部における癌組織の位置および大きさを検出するマンモグラフィを例として説明する。なお、図面において、同一の符号が付された構成要素は、同一または同種の構成要素を示す。
【００７３】
また、以下で説明する実施の形態は、本発明の好ましい一具体例を示す。以下の実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置、接続形態、ステップおよびステップの順序等は、一例であり、本発明を限定する主旨ではない。また、以下の実施の形態における構成要素のうち、本発明の最上位概念を示す独立請求項に記載されていない構成要素については、より望ましい形態を構成する任意の構成要素として説明される。
【００７４】
（実施の形態１）
＜検査装置の構成＞
以下、実施の形態１にかかる検査装置の構成について、図１〜３を用いて説明する。本実施の形態では、生体の内部情報、特に、欠陥組織の位置情報を映像化する例として、直線状の母線を持つ円錐体をモデルとする。具体的には、検査装置としてマンモグラフィ装置を、略円錐状の形状を有する対象物として乳房を、欠陥組織として乳房における癌組織を例として説明する。また、この装置において、センサ素子であるプローブを構成する送信アンテナ素子および受信アンテナ素子は、円錐体の母線に沿って移動されるものである。すなわち、送信アンテナ素子および受信アンテナ素子は、一体となって対象物のうち円錐の母線に相当する線に沿って移動される。ここで、円錐の母線に相当する線とは、対象物１０の形状を円錐とみなしたとき、その円錐の母線に相当する線である。
【００７５】
図１は、本実施の形態にかかるマンモグラフィ装置の構成を示す概略図である。図２は、本実施の形態にかかるマンモグラフィ装置におけるプローブの構成を示す概略図である。図３は、図２に示したプローブの使用方法の一例を示す図である。
【００７６】
はじめに、本実施の形態に係るマンモグラフィ装置１の構成について説明する。
【００７７】
図１に示すように、マンモグラフィ装置１は、プローブ３６と、検出部２０とを備えている。
【００７８】
プローブ３６は、対象物１０の内部情報を検出するセンサ素子であり、図２に示すように、一対の送信アンテナ素子３６ａと受信アンテナ素子３６ｂとで構成されている。送信アンテナ素子３６ａおよび受信アンテナ素子３６ｂはそれぞれ、複数の送信点および受信点を備えるアンテナアレイであっても良い。
【００７９】
プローブ３６は、図２に示すように、送信アンテナ素子３６ａと受信アンテナ素子３６ｂとが波動吸収手段（電波吸収手段）３８を介して接合された構成をしている。また、送信アンテナ素子３６ａと受信アンテナ素子３６ｂとは、それぞれ同軸ケーブル３４ａおよび３４ｂに接続されている。波動吸収手段は送信アンテナ素子３６ａで受信される波動を吸収する手段であり、波動がマイクロ波等の電波であるとき、電波吸収手段である。以下では、波動吸収手段が電波吸収手段である例について説明するが、この例に限定されない。
【００８０】
プローブ３６において送信アンテナ素子３６ａと受信アンテナ素子３６ｂとの相対位置は固定されている。ただし、プローブ３６の構成はこの構成に限定されない。送信アンテナ素子３６ａと受信アンテナ素子３６ｂとは、電波吸収手段３８を介して一体に設けられていることが好ましい。すなわち、プローブ３６において受信アンテナ素子３６ｂと送信アンテナ素子３６ａとは一体に設けられており、受信アンテナ素子３６ｂと送信アンテナ素子３６ａとの間には電波吸収手段３８が位置していることが好ましい。
【００８１】
「一体に設けられ」とは、プローブ３６に送信アンテナ素子３６ａと受信アンテナ素子３６ｂとが設けられ、送信アンテナ素子３６ａと受信アンテナ素子３６ｂとが一緒に移動できる状態をいう。送信アンテナ素子３６ａと受信アンテナ素子３６ｂとが一体に設けられている場合には、送信アンテナ素子３６ａと受信アンテナ素子３６ｂとが、電波吸収手段３８を介して、密着して、配置されている場合だけでなく、アンテナ素子と電波吸収手段３８との間に隙間がある場合が含まれる。アンテナ素子と電波吸収手段３８との間に隙間がある場合とは、送信アンテナ素子３６ａと電波吸収手段３８との間、および受信アンテナ素子３６ｂと電波吸収手段３８との間の少なくとも一方に隙間を設けた状態で、電波吸収手段が配置されている場合である。電波吸収手段３８としては、電波吸収体（波動吸収体）を送信アンテナ素子３６ａと受信アンテナ素子３６ｂとの間に配置して使用できる。
【００８２】
「電波吸収手段３８を介して」とは、送信アンテナ素子３６ａの受信アンテナ素子３６ｂに向いている一面、又は、受信アンテナ素子３６ｂの送信アンテナ素子３６ａに向いている一面に、電波吸収材料（波動吸収材料）を塗布等の手段で膜付けした場合や、電波吸収材料で製作した板状の部品を前記面上に接着して固着を行った場合も含むものとする。電波吸収手段３８は、送信アンテナ素子３６ａから発信されるマイクロ波（波動）のうち、直接、受信アンテナ素子３６ｂに受信されるマイクロ波（波動）の大きさを低減できればよい。電波吸収材料は導電性材料であり、具体的には、金属、カーボン等である。このような構成とすることで、送信アンテナ素子３６ａから発信されたマイクロ波（波動）のうち、直接受信アンテナ素子３６ｂに受信されるマイクロ波（波動）の大きさを小さくすることができ、その直接受信されるマイクロ波（波動）の前記大きさに対する、受信アンテナ素子３６ｂに受信される対象物１０からの反射波（散乱波）の大きさの比を高めることができる。すなわち、電波吸収手段３８を設けることにより、より鮮明な対象物内の画像を得ることができる。
【００８３】
符号発生器３２ａで生成された符号は、同軸ケーブル３４ａを介して、送信アンテナ素子３６ａに伝搬される。
【００８４】
アンテナ素子３６ｂで受信された信号は、同軸ケーブル３４ｂを介して、受信回路３２ｂに伝搬される。これにより、受信回路３２ｂは受信信号を検波する。
【００８５】
プローブ３６は、図３に示すように、ユーザの指示（動作）によって対象物１０の上を移動しながら対象物１０の内部情報の計測を行う。プローブ３６は、対象物１０に対して、自由に動かせるよう構成されている。送信アンテナ素子３６ａは、対象物１０に向かって波動を送信し、受信アンテナ素子３６ｂは、対象物１０の内部の散乱体１０ａにおいて散乱した散乱波を受信する。ここでは、波動としてマイクロ波を例として説明するが、マイクロ波に限らず他の帯域の電磁波または超音波であってもよい。
【００８６】
また、検出部２０は、図１に示すように、制御装置３０と、画像再構成部４０と、符号発生器３２ａと、受信回路３２ｂとを有している。
【００８７】
制御装置３０は、符号発生器３２ａから対象物１０への波動の送信、受信回路３２ｂによる散乱波の受信、および、画像再構成部４０による受信した散乱波の解析について制御を行う。制御装置３０は、符号発生器３２ａに対し、トリガ信号を伝搬する。
【００８８】
制御装置３０による制御によって、マンモグラフィ装置１は、送信アンテナ素子３６ａから対象物へ波動を送信し、受信アンテナ素子３６ｂで散乱体１０ａからの散乱波を受信し、さらに、画像再構成４０で散乱波を解析して画像を生成する。
【００８９】
符号発生器３２ａは、制御装置３０からのトリガ信号に基づいて、送信アンテナ素子３４ａから放射される波動を符号として生成し、送信アンテナ素子３４ａに搬送する。符号発生器３２ａは、符号として、例えば、送信アンテナ素子３６ａから放射される波動の送信タイミングや回数、送信ゲインを生成する。なお、送信アンテナ素子３６ａから対象物１０に放射される波動は、例えばマイクロ波である。
【００９０】
受信回路３２ｂは、受信アンテナ素子３６ｂで受信したマイクロ波の散乱波のデータを、画像再構成部４０に搬送する。このとき、受信した散乱波のデータは、受信回路３２ｂにより増幅されたり、制御装置３０においてＡＤ変換等の信号処理をされたりしてもよい。
【００９１】
画像再構成部４０は、受信回路３２ｂから搬送された散乱波のデータを解析し、後に説明する画像再構成アルゴリズムにより、散乱波のデータを映像化する。これにより、対象物１０の内部情報に応じた映像が画像表示器５０に再生される。
【００９２】
画像表示器５０は、モニタ画面であり、画像再構成部４０で演算されたデータを、映像として出力するものである。
【００９３】
なお、画像再構成アルゴリズムで使用される解析モデルとしては、円錐状、または、略円錐状の形状を有する解析モデルを想定してもよい。例えば、本実施の形態においては、一例として、底面の円の直径が３０ｃｍの円錐状の解析モデルを想定する。当該解析モデルにおいて、プローブ３６は、円錐の頂点から底面に向かう母線に沿ってライン状に移動されるとしてもよい。すなわち、送信アンテナ素子３６ａおよび受信アンテナ素子３６ｂは、少なくとも一方向から平面視したときに円錐の回転対象の軸と同一方向にライン状となるように移動されてもよい。ここで、円錐の回転対象の軸は、円錐の頂点と円錐の底面の中心とを結ぶ直線を指す。なお、画像再構成アルゴリズムについては、後に詳述する。
【００９４】
また、送信アンテナ素子３６ａと受信アンテナ素子３６ｂとは対をなして移動する。
【００９５】
以下、本実施の形態にかかるマンモグラフィ装置１による生体の内部情報、すなわち、乳房における癌組織の位置情報の観測（映像化）の手順について、説明する。
【００９６】
＜生体の内部情報の観測手順＞
図４は、図１に示した本実施の形態にかかるマンモグラフィ装置１の動作を示すフローチャート図である。
【００９７】
図４に示すように、マンモグラフィ装置１により、対象物１０である乳房の内部情報を映像化（画像化）する方法は、以下のとおりである。
【００９８】
はじめに、対象物１０の上にプローブ３６を配置する（Ｓ１０）。そして、送信アンテナ素子３６ａから対象物１０に向けて波動を放射する（Ｓ１１）。波動としては、例えば、マイクロ波を用いる。符号発生器３２ａにおいてマイクロ波の波長、振幅等の調整を行い、送信アンテナ素子３６ａから対象物１０に放射する。
【００９９】
次に、対象物１０の内部の正常組織および癌組織で反射した散乱波（散乱波）を受信アンテナ素子３６ｂにより受信する（Ｓ１２）。なお、正常組織および癌組織は、図３に示した散乱体１０ａに相当する。波動がマイクロ波等の電磁波である場合、正常組織と癌組織とでは、誘電率が異なるため、散乱波の強度が異なる。受信した散乱波は、受信回路３２ｂにおいて、増幅されたり、ＡＤ変換等の画像再構成部４０における解析に適した形式に変換されてもよい。
【０１００】
次に、受信回路３２ｂから画像再構成部４０に、受信した散乱波を示す散乱波データが搬送される。画像再構成部４０では、搬送された散乱波データの解析が行われる（Ｓ１３）。ここでは、以下に説明する画像再構成アルゴリズムにより散乱波データの解析が行われる。具体的には、映像化関数（画像化関数）が導出される。これにより、対象物１０の内部の正常組織と癌組織とに対応した映像（画像）が再構成される（Ｓ１４）。そして、プローブ３６を移動して、対象物１０の他の観測位置にプローブ３６を配置し（Ｓ１０）、上述した手順を繰り返す。なお、映像化関数（画像化関数）を用いた再構成で得られるデータは動画である必要は無く、画像、すなわち静止画であってよい。
【０１０１】
さらに、再構成された映像のデータは、画像再構成部４０からモニタ５０に搬送され、モニタ５０に表示される。これにより、対象物１０の内部の癌組織の存在およびその位置、形状および大きさが確認できる。
【０１０２】
上記した通り、画像再構成部４０は受信アンテナ素子３６ｂにより受信した散乱波を示す散乱波データから、対象物１０の内部情報に関する画像を再構成する。詳しくは画像再構成部４０は、散乱波データと形状を示す解析モデルとを用いて偏微分方程式を解くことにより再構成関数を導出する。そして、再構成関数により、対象物１０の内部情報に関する画像を再構成する。ここで、偏微分方程式は、対象物の内部情報に関する画像を再構成するための再構成関数が満たす方程式である。
【０１０３】
具体的には、偏微分方程式は、対象物内の送信点および受信点の位置を示す独立変数を用いて表される式である。また、偏微分方程式は、独立変数の数と同一の次元を持つ空間における各点において、散乱波が生じている場の関数である散乱場関数が解となる線形偏微分方程式である。そして、画像再構成部は、再構成関数の時間変数の極限値である映像化関数を導出し、映像化関数を用いて対象物の内部情報に関する画像を再構成する。
【０１０４】
以下、画像再構成部４０で行われる画像再構成アルゴリズムについて説明する。この画像再構成アルゴリズムは、直線状の母線を持つ錐体を解析モデルとする、本実施の形態に係るマンモグラフィ装置１の画像再構成原理である。
【０１０５】
＜画像再構成アルゴリズム＞
図５〜図７は、本実施の形態に係る散乱トモグラフィ方法の原理を説明するための解析モデルである。
【０１０６】
以下に説明する画像再構成アルゴリズム（理論）では、映像化対象が直線状の母線を持つ円錐体であると仮定する。この円錐体を解析モデルとして、当該母線上に波動の放射点（送信アンテナ素子３６ａ）と受信点（受信アンテナ素子３６ｂ）とを任意に設置する。すなわち、プローブ３６が円錐体の母線上を移動するという解析モデルとする。なお、この画像再構成アルゴリズムは、放射点と受信点とが同じ位置であることを想定したモノスタティックな解析を行うものである。そして、放射点から放射された放射波を示す送信データと受信点で受信された散乱波を示す散乱波データとを用いて、映像化対象物の内部の情報を映像化する。
【０１０７】
画像再構成アルゴリズムについて数学的に簡単に説明すると、映像化に必要な解（関数）をはじめに設定し、その解から方程式を構築し、送信データおよび受信データから、より厳密な解を求めるという、逆問題を解くものである。
【０１０８】
すなわち、はじめに映像化に必要なグリーン関数（再構成関数）を設定する。そして、この関数が解（関数）となる（ｔ,ｘ,ｙ,ｚ）の４つの変数からなる３次元空間に関する偏微分方程式を構築する。そして、送信アンテナ素子３６ａから放射された送信データと受信アンテナ素子３６ｂで受信された受信データ（散乱データ）とを境界条件としてこの偏微分方程式を解く。これにより、映像化関数が求まり、物体内部の情報を高品質かつ高速に映像化することができる。
【０１０９】
具体的には、以下のとおりである。
【０１１０】
１．散乱の逆問題とグリーン関数
図５により、点ｒ１から出た波動が点ξで反射して点ｒ２へ戻ってくる状況を考える。ここで、本実施の形態にかかるマンモグラフィ装置１では、点ｒ１と点ｒ２とは同一の位置である。周波数ωが一定という条件で、波動の送信点ｒ１と受信点ｒ２がｘ断面Ｄ（円錐体６０ａの側面）の内部を自由に動く。このとき得られるデータをＧ（ｒ１,ｒ２,ω）と書くと、この関数は領域内の反射点の分布に関係したものとなる。ここで、ωは角周波数であり、ω＝２πｆである。Ｇ（ｒ１,ｒ２,ω）は、全ての点ξからの反射信号の和であり、領域内には多くの反射点があるので、Ｇ(ｒ１,ｒ２,ω)は次の（式１）のように書くことができる。
【０１１１】
【数１５】
2017057524.tif 000017
....（式１）
【０１１２】
ここで、
【数１６】
2017057524.tif 000018
は、点ｒ１から出て点ξで反射して点ｒ２へ戻る波動の信号強度を示す。なお、式中のパラメータｒ１およびｒ２はベクトルであり、式中では太字で示している。
【０１１３】
ここで、波動の送信点ｒ１と受信点ｒ２の座標は常に等しい。
【０１１４】
この関数Ｇ(ｒ,ω)≡Ｇ(ｒ１,ｒ２,ω)を用いて散乱の逆問題の理論的構造を説明する。
【０１１５】
３次元空間の部分領域をＤとし、その境界を∂Ｄとする。このとき、関数Ｇ(ｒ,ω)は、領域Ｄの内部で次の（式２）のような微分方程式の解となっている。
【０１１６】
【数１７】
2017057524.tif 000019
.......（式２）
【０１１７】
さらに、Ｇ(ｒ,ω)の境界∂Ｄでの値は、プローブ３６における測定値（送信データおよび受信データ）とする。上記方程式をこの境界条件のもとに解き、その結果から求めようとする領域Ｄ内の誘電率の勾配に関係した関数ρ（ｒ）を定義する。実際には、ここに現れる微分作用素Ｌ（∂／∂ｔ,∂／∂ｒ１,∂／∂ｒ２）を求める必要がある。
【０１１８】
２．円錐体を解析モデルとする、回転対称曲面上のマルチパス逆散乱理論
次に、上記した微分作用素を求める方法について説明する。図６および図７は、微分作用素を求める方法について説明するための解析モデルである。
【０１１９】
本実施の形態にかかる解析モデルは、図６に示すように、映像化対象が直線状の母線を持つ円錐体６０ａであると仮定する。プローブ３６は、円錐体６０ａの母線に沿ってライン状に移動する。また、プローブ３６は、円錐体６０ａの頂点と円錐体６０ａの底面の中心とを通る軸を中心として回転移動もしてもよい。すなわち、送信アンテナ素子３６ａおよび受信アンテナ素子３６ｂは、少なくとも一方向から平面視したときに解析モデルとする円錐体６０ａの回転対象の軸と同一方向にライン状に移動しながら計測を行う。ここで、円錐体の回転対象の軸は、円錐体の頂点と円錐体の底面の中心とを結ぶ直線を指す。
【０１２０】
詳細には、プローブ３６の位置での円錐体６０ａの接平面８０ａを想定し、仮想的な接平面８０ａが回転し、全ての回転角度θに対して散乱波データが得られるものとする。このような円錐体６０ａの表面におけるモノチスタティックな時系列データから、円錐体６０ａの内部の３次元的な構造を再構成する。その理論を以下に説明する。
【０１２１】
なお、本解析モデルでは、図６に示すように、円錐体６０ａの頂点を点Ｏ、円錐体６０ａの底面の円の中心を点Ｏ'とする。点Ｏから点Ｏ'に向かう方向をＺ方向（Ｚ軸方向）とする。また、点Ｏを通り円錐体６０ａの底面と平行な面を基準平面７０ａとする。基準平面７０ａにおける任意の一方向をＸ方向（Ｘ軸方向）、基準平面７０ａにおけるＸ方向と垂直な方向をＹ方向（Ｙ軸方向）とする。
【０１２２】
また、円錐体６０ａの底面において、Ｘ方向と平行な方向をＸ'方向（Ｘ'軸方向）、Ｙ方向と平行な方向をＹ'方向（Ｙ'軸方向）とする。
【０１２３】
また、円錐体６０ａの側面において、プローブ３６の位置で円錐体６０ａに接する平面を、接平面８０ａとする。この接平面８０ａと基準平面７０ａとが交差する交差線の一方向をｘ方向（ｘ軸方向）とする。接平面８０ａにおいて、点Ｏから円錐体６０ａの底面に向かう方向をｙ方向（ｙ軸方向）とする。また、ｘ方向およびｙ方向に垂直な方向をｚ方向（ｚ軸方向）とする。
【０１２４】
また、基準平面７０ａにおいて、ｘ軸はＺ軸を中心に回転しているものとし、Ｘ軸とｘ軸との間の角をθとする。また、Ｚ軸とｚ軸との間の角をαとする。
【０１２５】
図７は、図６に示す円錐体６０ａの母線上を移動するプローブ３６の軌跡を部分的に拡大した模式図である。
【０１２６】
図７において、ｙ軸上にモノスタティック座標ｙ１およびｙ２（ｙ１＝ｙ２＝ｙ）を定義すると、以下のような再構成理論が成立する。なお、以下の理論では、接平面８０ａ上において、ｘ座標に関しては１点（ｘ＝０,ｙ１,ｙ２）のデータが得られる。
【０１２７】
図７に示すように、プローブ３６の点Ｐ１（ｘ,ｙ１,ｚ）＝（ｘ,ｙ,ｚ）における送信アンテナ３６ａから放射した波動は、点Ｐで反射して、リニアアレイアンテナ３６の点Ｐ２（ｘ,ｙ２,ｚ）＝Ｐ１（ｘ,ｙ,ｚ）における受信アンテナ３６ｂで受信される。また、計測点Ｐ１、Ｐ２は、リニアアレイアンテナ３６が配置された任意直線上を動くと仮定する。
【０１２８】
任意直線上ではｒ１、ｒ２のｘ、ｙ、ｚ座標は、r１＝（ｘ,ｙ１,ｚ）、ｒ２＝（ｘ,ｙ２,ｚ）と表される。このとき、関数Ｇは次の（式３）のように定義される。
【０１２９】
【数１８】
2017057524.tif 000020
................（式３）
【０１３０】
次に、関数Ｇ（ｒ,ω）≡Ｇ（ｒ１,ｒ２,ω）が満たす方程式として、次の（式４）のような関数φを導入する。すなわち、点Ｐ１が領域全体を動くときＰ１で受信される信号は、次の（式４）ように表すことができる。ここで、ｃは伝播速度、ｋは波数であり、波長をλとすると、ω＝ｃｋ、ｋ＝２π／λの関係がある。なお、（式４）に示すφは、本発明における物体の内部情報に関する画像を再構成するための再構成関数（解）に相当する。ξ、η、ζは、図６および図７に示す点Ｐ（ξ，η，ζ）のｘ座標、ｙ座標およびｚ座標である。なお、点Ｐ（ξ，η，ζ）は、領域中の任意の散乱点である。また、ベクトルｒ１およびベクトルｒ２について、ｒ１＝ｒ２＝ｒ（ｘ，ｙ，ｚ）が成り立つ。
【０１３１】
【数１９】
2017057524.tif 000021
...........（式４）
【０１３２】
ここで、時間の因子をｅｘｐ（−ｉωｔ）に比例すると仮定する。上式の被積分項にあるカーネル関数をΦと置く。
【０１３３】
【数２０】
2017057524.tif 000022
【０１３４】
この式、及び、この式をξ、η、ζなどで微分および積分したものが解となる偏微分方程式を求める。そのためには、微分した結果生ずる１／ρに関して高次の項を無視して計算すればよい。
【０１３５】
ここで、微分の略記法を、以下の（式５）ように定義する。
【０１３６】
【数２１】
2017057524.tif 000023
............（式５）
【０１３７】
この（式５）を用いてΦの各次数の微分を求めると、以下の（式６）のようになる。
【０１３８】
【数２２】
2017057524.tif 000024
........（式６）
【０１３９】
以下の式では、繁雑なｏ（＊）の項は省略する。なお、＊は任意の変数を意味する。２階微分に関する４つの式の和をとると、以下の（式７）のようになる。
【０１４０】
【数２３】
2017057524.tif 000025
............（式７）
【０１４１】
（式７）を整理すると、（式８）のようになる。
【０１４２】
【数２４】
2017057524.tif 000026
............（式８）
【０１４３】
この（式８）は定常状態を仮定して導かれたが、非定常の場合へ拡張することは容易である。（式８）を非定常の場合へ拡張するために、（式８）について次の（式９）のような変数の置換えをする。
【０１４４】
【数２５】
2017057524.tif 000027
..........（式９）
【０１４５】
こうして、最終的に次の（式１０）に示す方程式を得る。なお、（式１０）は、本発明における偏微分方程式に相当する。
【０１４６】
【数２６】
2017057524.tif 000028
...................（式１０）
【０１４７】
（式１０）において微分を積分核へ適用することにより、φもまた上記偏微分方程式を満たす。この方程式は、変数がｔ、ｘ、ｙ、ｚの４個からなる３次元の擬似波動方程式である。
【０１４８】
この（式１０）を、フーリエ変換を用いて解く。まず最初に、φをｔ、ｘ、ｙについて多重フーリエ変換すると、以下の（式１１）が得られる。
【０１４９】
【数２７】
2017057524.tif 000029
..................（式１１）
【０１５０】
（式１１）において、ｚに関する微分をＤｚと書くと、以下の（式１２）を得る。
【０１５１】
【数２８】
2017057524.tif 000030
.........（式１２）
【０１５２】
ここで、ω＝ｃｋの関係を用いると、この方程式（式１２）の２つの基本解は次のようになる。
【０１５３】
【数２９】
2017057524.tif 000031
................（式１３）
【０１５４】
（式１３）において、時間のファクターがｅ−ｉωｔであること、および放射電波のパスで位相が加算されること、また、物体で反射した電波が測定面へ向かって跳ね返っていくことを考慮すると、Ｅ１が唯一の意味のある解となる。従って、次の（式１４）が得られる。
【０１５５】
【数３０】
2017057524.tif 000032
.......（式１４）
【０１５６】
これをフーリエ変換して、波動方程式（式１０）の解φが次の（式１５）のように求まる。
【０１５７】
【数３１】
2017057524.tif 000033
.............（式１５）
【０１５８】
計測されるデータは、ｘ＝０にしかないから、次の（式１６）が成り立つ。
【０１５９】
【数３２】
2017057524.tif 000034
...............（式１６）
【０１６０】
これを、上式でｚ＝０とおいたものへ適用すると、次の（式１７）のようになる。
【０１６１】
【数３３】
2017057524.tif 000035
.............（式１７）
【０１６２】
この（式１７）の両辺を（ｘ,ｙ）でフーリエ変換すると、ａ（ｋｘ,ｋｙ,ｋ）が次のように求まる。
【０１６３】
【数３４】
2017057524.tif 000036
............（式１８）
【０１６４】
こうして、偏微分方程式の解φ（ｘ,ｙ,ｚ,ｋ）が得られる。
【０１６５】
【数３５】
2017057524.tif 000037
..........（式１９）
【０１６６】
ｋに関する積分を遂行するために、次の（式２０）のような変数変換を行う。
【０１６７】
【数３６】
2017057524.tif 000038
........（式２０）
【０１６８】
次に、（式２０）をｋに関してフーリエ変換し、ｔ＝０を適用すると、（式２１）に示すように、角度θにおける局所座標系での映像化関数ρ（ｒ,θ）が求まる。
【０１６９】
ここで、
【数３７】
2017057524.tif 000039
はθに依存するため、θ依存性を明示するために、
【数３８】
2017057524.tif 000040
のように書くことにする。
【０１７０】
【数３９】
2017057524.tif 000041
..............（式２１）
【０１７１】
さらに、各θで得られた結果を積分することにより、（式２２）に示す３次元再構成画像が得られる。
【０１７２】
【数４０】
2017057524.tif 000042
...........（式２２）
【０１７３】
次に、接空間で計算された結果を全体座標（Ｘ,Ｙ,Ｚ）へ変換する。ｙ軸の（Ｘ,Ｙ）平面への射影をｙ'とすると、次の（式２３）が成立する。
【０１７４】
【数４１】
2017057524.tif 000043
............（式２３）
【０１７５】
さらに（ｘ,ｙ'）から（Ｘ,Ｙ）へ変換する式は、次の（式２４）のようになる。
【０１７６】
【数４２】
2017057524.tif 000044
.............（式２４）
【０１７７】
以上をまとめると、ｘ、ｙ、ｚは、次の（式２５）のようになる。
【０１７８】
【数４３】
2017057524.tif 000045
.............（式２５）
【０１７９】
この変換公式を、次の（式２６）へ適用していく。
【０１８０】
【数４４】
2017057524.tif 000046
............（式２６）
【０１８１】
さらに、スペクトラム領域で、次の（式２７）に示す変数変換を行う。
【０１８２】
なお、以下の（式２７）および（式２８）で使用する変数ξ、η、ζは、図７および（式２２）等で使用したＰ点の座標（ξ，η，ζ）ではなく、変数変換のために新たに用いた変数である。
【０１８３】
【数４５】
2017057524.tif 000047
.............（式２７）
【０１８４】
この逆変換は、次の（式２８）で与えられる。
【０１８５】
【数４６】
2017057524.tif 000048
.............（式２８）
【０１８６】
角度θにおける再構成関数は、次の（式２９）のようになる。
【０１８７】
【数４７】
2017057524.tif 000049
.............（式２９）
【０１８８】
ここで、ｋｘ、ｋｙ、ｋｚなどは、（式２６）、（式２８）のようにξ、η、ζの関数であるから、フーリエ変換だけを用いて、角度θの局所座標におけるデータを全体座標へ変換できる。
【０１８９】
最後に角度θで積分して、（式３０）に示すような再構成画像、すなわち、映像化関数が得られる。
【０１９０】
【数４８】
2017057524.tif 000050
.............（式３０）
【０１９１】
この映像化関数により、物体の内部情報に関する画像を再構成することができる。
【０１９２】
＜検査装置により得られる再構成画像＞
上記したマンモグラフィ装置１により再構成画像を得ることの効果について、説明する。
【０１９３】
本実施の形態にかかるマンモグラフィ装置１は、特に、乳房中の欠陥組織の観測に有効である。
【０１９４】
乳房における欠陥組織を観測する手段としては、従来から、Ｘ線、超音波またはＭＲＩによる観測方法があるが、Ｘ線による観測では、Ｘ線は、有機組織の種別の判別がむずかしく、乳房中の欠陥組織の観測には適していないという問題がある。また、超音波による観測では、領域の大分部を占める脂質中での減衰が大きいため構成画像のＳ／Ｎの向上が難しく、乳房中の欠陥組織の観測には適していないという問題がある。また、ＭＲＩによる観測では、造影剤が必要であること、鮮明な構成画像が得られるが、構成画像のうちのどれが欠陥組織であるかの判断が難しいという点から、乳房中の欠陥組織の観測には適していないという問題がある。加えて、強磁場の発生は、超伝導コイルとその冷却システムとが必要であり、大型で高コストである。
【０１９５】
これに対し、本実施の形態にかかるマンモグラフィ装置１による観測では１〜１０ＧＨｚのＵＷＢ（ウルトラワイドバンド）マイクロ波を用いる。マイクロ波は、生体、特に、脂質等に対してマイクロ波の減衰が極めて小さいため、乳房中の欠陥組織の観測に有効である。また、装置構成がＭＲＩに比べると簡単であり、造影剤など他の物質も必要ないため、汎用的に用いることができる。さらに、本実施の形態にかかるマンモグラフィ装置１では、上記した逆散乱理論を用い、マイクロ波の多重経路（マルチスタティック）散乱波データを基にして、乳房などの生体内部の３次元的構成要素を映像化するため、生体の内部情報（構成）を汎用的かつ高速に映像化することができる。
【０１９６】
なお、上記したマンモグラフィ装置１の確認実験の例としては、これまでのところ、例えば、ケースにアルミナボールを充填した生体モデルにおける観測、および、動物の乳房の観測を行っている。以下、観測結果について説明する。
【０１９７】
図８は、実施の形態１にかかるマンモグラフィ装置１を用いた計測実験のための実験モデルを示す概略図である。図９は、図８に示した実験モデルをマンモグラフィ装置で計測した計測結果を示す図である。
【０１９８】
実験モデル１００は、豚の乳房を想定し、図８の（ａ）および（ｂ）に示すように、内径３０ｃｍ、高さが１０ｃｍ程度のプラスチック容器１００ａ内に、直径が１ｍｍで比誘電率が６のアルミナボール１００ｂを充填し、さらに、アルミナボール１００ｂ内に埋め込むように、断面直径が５ｍｍの水入りゴムホース１００ｃを５本配置したものである。詳細には、図８の（ｂ）における高さｄはｄ＝３４ｍｍ、高さｈはｈ＝５４．６ｍｍ、角度αはα＝２０°である。
【０１９９】
測定条件の一例として、送信アンテナから送信するマイクロ波には、１０ＧＨｚ以下のマイクロ波を用いた。この場合、送信アンテナ３６ａおよび受信アンテナ３６ｂは、実験モデル１００に当接する部分の大きさが１５ｍｍ×１５ｍｍ程度の大きさのものを使用した。この条件の下で、本実施の形態にかかるマンモグラフィ装置１により実験モデル１００の観測を行ったところ、図９に示すように、実験モデル１００に配置した５本の水入りゴムホース１００ｃを観測することができた。また、図示を省略しているが、マイクロ波の周波数を変更することにより、深さ方向の観測位置を変更しながら観測を行うことができた。図９に示す観測結果では、表面からの深さが５４ｍｍ程度のところで水入りゴムホース１００ｃが観測された。図８のモデルに対してこの観測結果は妥当である。さらに、この観測結果から、上記したＸ線、超音波およびＭＲＩによる観測で得られた画像よりも精度のよい再構成画像が得られることが確認できた。
【０２００】
以上より、本実施の形態にかかるマンモグラフィ装置１は、特に、乳房中の欠陥組織の観測に有効である。
【０２０１】
以上、本実施の形態に係るマンモグラフィ装置１によると、円錐体６０ａの母線に沿ってライン状にプローブ３６を移動する解析モデルにおいて、逆問題用の偏微分方程式を設定し、これを解くことで映像化関数を得ている。これにより、本実施の形態にかかるマンモグラフィ装置１において、対象物１０の内部情報を汎用的かつ高速に映像化することができる。
【０２０２】
具体的な効果としては、本実施の形態にかかるマンモグラフィ装置１は、波動としてマイクロ波を用いているので、Ｘ線を用いた観測方法と比較して、生体に放射したときに被爆の可能性がなく安全である。また、Ｘ線を用いた観測方法と比較して、画像化したときに癌組織と正常組織とのコントラスト比が大きいため、高感度かつ短時間で観測することができる。また、プローブを検査対象箇所に当接して観測を行うことができるので、装置の小型化および低コスト化を達成することができる。
【０２０３】
また、ユーザすなわち医師等の診断者の主観や技量に左右されない画像解析を行うことができるため、診断者が異なる場合であっても正確に診断することができる。
【０２０４】
また、検査対象箇所を圧迫することなくプローブを検査対象箇所に当接して観測を行うことができるので、患者の苦痛を伴うことなく、かつ、着衣のままで簡便に観測を行うことができる。
【０２０５】
なお、上記の演算式、および、演算式の導出手順は、一例であって、別の演算式、および、別の導出手順が用いられてもよい。
【０２０６】
また、本実施の形態では、波動としてマイクロ波を用いたが、マイクロ波に限らず、他の周波数帯域の電磁波であってもよいし、超音波であってもよい。また、本実施の形態では、マイクロ波を用いているため、所定の周波数を有する周期波を用いているが、波動として周期波に限らずパルス波を用いてもよい。
【０２０７】
また、本実施の形態では、対象物として乳房を例に説明しているが、対象物は乳房に限らず、他の生体や物体、例えば、円錐状のコンクリート支柱などであってもよい。
【０２０８】
（実施の形態２）
次に、本発明の実施の形態２について説明する。
【０２０９】
本実施の形態に係るマンモグラフィ装置１の構成は、実施の形態１にかかるマンモグラフィ装置１とほぼ同様であるが、実施の形態１に示した円錐体の解析モデルにおける円錐体の母線が曲線状であり、この曲線状の母線に沿って、送信アンテナ素子３６ａおよび受信アンテナ素子３６ｂが移動する点が実施の形態１に示したマンモグラフィ装置１と異なっている。したがって、本実施の形態に係るマンモグラフィ装置１の画像再構成部で行われる画像再構成アルゴリズムも、実施の形態１にかかるマンモグラフィ装置１と異なっている。
【０２１０】
実施の形態２に係るマンモグラフィ装置１では、送信アンテナ素子３６ａおよび受信アンテナ素子３６ｂは、一体となって対象物１０のうち略円錐の曲線状の母線に相当する線に沿って移動される。ここで、略円錐とは、曲線状の母線を持つ円錐であり、略円錐の母線に相当する線とは、対象物１０の形状を曲線状の母線を有する略円錐とみなしたとき、その略円錐の曲線状の母線に相当する線である。
【０２１１】
以下、本実施の形態に係るマンモグラフィ装置１の画像再構成アルゴリズムについて説明する。
【０２１２】
＜画像再構成アルゴリズム＞
図１０および図１１は、本実施の形態に係る散乱トモグラフィ方法の原理を説明するための解析モデルである。以下、図１０および図１１に示すモデルを解析モデルとして、本実施の形態にかかるマンモグラフィ装置の映像化関数の導出について説明する。
【０２１３】
以下に説明する画像再構成アルゴリズム（理論）では、映像化対象（対象物１０）が曲線状の母線を持つ円錐体、すなわち半球状またはドーム状の円錐体であると仮定する。この円錐体を解析モデルとして、当該母線上に波動の放射点（一の送信アンテナ素子３６ａ）と受信点（一の受信アンテナ素子３６ｂ）を任意に設置する。すなわち、プローブ３６が曲線状の母線を有する円錐体の母線上を移動するという解析モデルとする。なお、この画像再構成アルゴリズムは、放射点と受信点とが同じ位置であることを想定したモノスタティックな解析を行うものである。そして、放射点から放射された放射波を示す送信データと受信点で受信された散乱波を示す散乱波データとを用いて、映像化対象物の内部の情報を映像化する。
【０２１４】
具体的には、本実施の形態にかかる解析モデルは、図１０に示すように、映像化対象が曲線状の母線を持つ円錐体（略円錐体）６０ｂであると仮定する。プローブ３６は、円錐体６０ｂの母線に沿って移動する。また、プローブ３６は、円錐体６０ｂの頂点と円錐体６０ｂの底面の中心とを通る軸を中心として回転移動してもよい。すなわち、複数の送信アンテナ素子３６ａと複数の受信アンテナ素子３６ｂとで構成されるプローブ３６は、少なくとも一方向から平面視したときに解析モデルとする円錐体６０ｂの回転対象の軸と同一方向にライン状に移動しながら計測を行う。ここで、円錐体の回転対象の軸は、円錐体の頂点と円錐体の底面の中心とを結ぶ直線を指す。
【０２１５】
詳細には、プローブ３６の位置での円錐体６０ｂの接平面８０ｂを想定し、仮想的な接平面８０ｂが円錐体６０ｂの頂点と円錐体６０ｂの底面の中心とを通る軸を中心として回転し、全ての回転角度θに対して散乱波データが得られるものとする。このような円錐体６０ｂの表面におけるモノスタティックな時系列データから、円錐体６０ｂの内部の３次元的な構造を再構成する。その理論を以下に説明する。
【０２１６】
なお、本解析モデルでは、図１０に示すように、円錐体６０ｂの頂点を点Ｏ、円錐体６０ｂの底面の円の中心を点Ｏ'とする。点Ｏから点Ｏ'に向かう方向をＺ方向（Ｚ軸方向）とする。また、円錐体６０ｂの側面において、プローブ３６の位置で円錐体６０ｂに接する平面を、接平面８０ｂとする。この接平面８０ｂがＺ軸と交差する点を含み、円錐体６０ｂの底面と平行な平面を基準面７０ｂとする。基準平面７０ｂにおける任意の一方向をＸ方向（Ｘ軸方向）、基準平面７０ｂにおけるＸ方向と垂直な方向をＹ方向（Ｙ軸方向）とする。
【０２１７】
また、円錐体６０ｂの底面において、Ｘ方向と平行な方向をＸ'方向（Ｘ'軸方向）、Ｙ方向と平行な方向をＹ'方向（Ｙ'軸方向）とする。
【０２１８】
この接平面８０ｂと基準平面７０ｂとが交差する交差線の一方向をｘ方向（ｘ軸方向）とする。接平面８０ｂにおいて、点Ｏから円錐体６０ｂの底面に向かう方向をｙ方向（ｙ軸方向）とする。また、ｘ方向およびｙ方向に垂直な方向をｚ方向（ｚ軸方向）とする。
【０２１９】
また、基準平面７０ｂにおいて、ｘ軸はＺ軸を中心に回転しているものとし、Ｘ軸とｘ軸との間の角をθとする。また、Ｚ軸とｚ軸との間の角をαとする。
【０２２０】
図１１は、図１０に示す円錐体６０ｂの母線上を移動するプローブ３６の軌跡を部分的に拡大した模式図である。図１１に示すように、ｘ、ｙ座標は対象物の表面に置き、ｚ座標は対象物の表面の法線方向にとる。この解析モデルでは、プローブ３６は、ｘｙｚ空間において曲面状に移動する。
【０２２１】
本実施の形態にかかる画像再構成アルゴリズムにおいては、実施の形態１に係るマンモグラフィ装置１と同様、はじめに映像化に必要なグリーン関数（再構成関数）を設定する。そして、この関数が解（関数）となる（ｔ，ｘ，ｙ，ｚ）４つの変数からなる３次元空間に関する偏微分方程式を構築する。そして、送信アンテナ素子３６ａから放射された送信データと受信アンテナ素子３６ｂで受信された受信データ（散乱データ）とを境界条件としてこの偏微分方程式の厳密な映像化関数を求める。これにより、物体内部の情報を高品質かつ高速に映像化することができる。
【０２２２】
具体的には、以下のとおりである。
【０２２３】
１．散乱の逆問題とグリーン関数
図１０により、点ｒ１から出た波動が点ξで反射して点ｒ２へ戻ってくる状況を考える。ここで、本実施の形態にかかるマンモグラフィ装置１では、点ｒ１と点ｒ２とは同一の位置である。周波数ωが一定という条件で、波動の送信点ｒ１と受信点ｒ２がある拘束条件を満たしながら曲面（円錐体６０ｂの側面）上を自由に動く。このとき得られるデータをＧ（ｒ１，ｒ２，ω）と書くと、この関数は領域内の反射点の分布に関係したものとなる。このとき、ωは角周波数＝２πｆである。Ｇ（ｒ１，ｒ２，ω）は、全ての点ξからの反射信号の和であり、領域内に多くの反射点があるので、Ｇ（ｒ１，ｒ２，ω）は、実施の形態１に示した（式１）と同様に考えることができる。
【０２２４】
ただし、本実施の形態では、波動の送信点ｒ１と受信点ｒ２に課せられた拘束条件とは、ｒ１とｒ２がある曲面上を動くものであり、かつ、ｒ１とｒ２が略同一の位置（座標）であるという条件である。
【０２２５】
この関数Ｇ（ｒ１,ｒ２,ω）を用いて散乱の逆問題の理論的構造を説明する。
【０２２６】
２．曲線状の母線を持つ円錐体を解析モデルとする、回転対称曲面上のマルチパス逆散乱理論
以下、この微分作用素を求める方法について述べる。本実施の形態に係る解析モデルでは、図１０および図１１に示したように、実施の形態１の解析モデルにおける円錐体の母線が曲線状であり、この曲線状の母線に沿って、送信アンテナおよび受信アンテナが移動する。そして、アンテナアレイが、Ｚ軸を中心に回転移動もしてもよい。曲線状の母線の上ではｒ１、ｒ２のｘ、ｙ、ｚ座標は、ｒ１＝（ｘ,ｙ１,ｚ１）（＝（ｘ,ｙ,ｚ））、ｒ２＝（ｘ,ｙ２,ｚ２）（＝（ｘ,ｙ,ｚ））と表される。
【０２２７】
ここで、関数Ｇは、実施の形態１に示した（式３）と同様に定義され、ｒ１＝（ｘ,ｙ１,ｚ１）、ｒ２＝（ｘ,ｙ２,ｚ２）から関数Ｇ（ｒ１,ｒ２,ω）の満たす方程式を求める。関数Ｇ（ｒ１,ｒ２,ω）を満たす方程式として、次の（式３１）のような関数φを導入する。なお、（式３４）に示すφは、本発明における物体の内部情報に関する画像を再構成するための再構成関数（解）に相当する。
【０２２８】
図１１において、断面曲線Ｓ上の点ＰＩで放射され、点Ｐで反射した後、点ＰＩで受信される波動を考えると、関数φは、次の（式３１）ように書くことができる。なお、ベクトルｒ１およびベクトルｒ２について、ｒ１＝ｒ２＝ｒ（ｘ，ｙ，ｚ）が成り立つ。
【０２２９】
【数４９】
2017057524.tif 000051
...............（式３１）
【０２３０】
ここで、時間の因子をｅｘｐ−ｉωｔに比例すると仮定し、波数をｋとおく。また、上式の被積分項にあるカーネル関数をφ置く。
【０２３１】
【数５０】
2017057524.tif 000052
..................（式３２）
【０２３２】
この式、及び、この式をξ、η、ζなどで微分および積分したものが解となる偏微分方程式を求める。そのためには、微分した結果生ずる１／ρに関して高次の項を無視して計算すればよい。
【０２３３】
ここで、微分の略記法を、以下の（式３３）のように定義する。
【０２３４】
【数５１】
2017057524.tif 000053
.............（式３３）
【０２３５】
この（式３３）を用いてΦの各次数の微分を求めると、以下の（式３４）のようになる。
【０２３６】
【数５２】
2017057524.tif 000054
...........（式３４）
【０２３７】
以下の式では、繁雑なｏ（＊）の項は省略する。なお、＊は任意の変数を意味する。２階微分に関する４つの式の和をとると、以下の（式３５）のようになる。
【０２３８】
【数５３】
2017057524.tif 000055
..............（式３５）
【０２３９】
従って、
【数５４】
2017057524.tif 000056
..............（式３６）
【０２４０】
さらに、式を整理することにより最終的に次の（式３７）が得られる。なお、（式３７）は、本発明における偏微分方程式に相当する。
【０２４１】
【数５５】
2017057524.tif 000057
.............（式３７）
【０２４２】
φの時間の因子をｅｘｐ−ｉωｔに比例すると仮定して、（式３７）の解を求めることにする。まず最初に、φをｔ、ｘ、ｙについて多重フーリエ変換すると、以下の（式３８）が得られる。
【０２４３】
【数５６】
2017057524.tif 000058
..............（式３８）
【０２４４】
ｚ１、ｚ２に関する偏微分をそれぞれＤｚ１、Ｄｚ２と書くと以下の式を得る。
【０２４５】
【数５７】
2017057524.tif 000059
.......（式３９）
【０２４６】
この方程式を解くのであるが、ω＝ｃｋの関係を用いると、上記方程式の２つの基本解は以下の（式４０）のようになる。
【０２４７】
【数５８】
2017057524.tif 000060
.........（式４０）
【０２４８】
ここで、時間のファクターがｅ−ｉωｔであること、放射電波のパスで移送が加算されること、および、物体で反射した電波が測定面へ向かって跳ね返っていくことを考慮すると、Ｅ１が唯一の意味のある解となる。したがって、（式３７）の解は、逆フーリエ変換により、次の（式４１）のように書くことができる。
【０２４９】
【数５９】
2017057524.tif 000061
.........（式４１）
【０２５０】
ｘを固定した断面曲線Ｓの方程式を、次の（式４２）のように仮定する。
【０２５１】
【数６０】
2017057524.tif 000062
..（式４２）
【０２５２】
断面曲線Ｓ上で与えられる境界条件は、次の（式４３）のようになる。
【０２５３】
【数６１】
2017057524.tif 000063
....（式４３）
【０２５４】
ここで、積分方程式（式５６）の解を求める。曲面上の点ＰＩで測定された時系列データφ（ｘ,ｙＩ,ｚＩ,ｔ)をフーリエ変換した関数Φを次の（式４４）ように書くことができる。
【０２５５】
【数６２】
2017057524.tif 000064
..............（式４４）
【０２５６】
ここでｚＩ,ｚＪは次の式を満たす。
【０２５７】
【数６３】
2017057524.tif 000065
..............（式４５）
【０２５８】
さらに、ｘ＝０以外ではデータが無いので、φは次のように書ける。
【０２５９】
【数６４】
2017057524.tif 000066
..............（式４６）
【０２６０】
すると（式４４）は次のようになる。
【０２６１】
【数６５】
2017057524.tif 000067
..............（式４７）
【０２６２】
この（式４７）の右辺には、ｋｘが含まれていない。しかし、ある角度θでのデータであることを考慮して、右辺を次の（式４８）のような関数で表す。
【０２６３】
【数６６】
2017057524.tif 000068
..............（式４８）
【０２６４】
この記号を用いると、次の（式４９）が得られる。
【０２６５】
【数６７】
2017057524.tif 000069
.............（式４９）
【０２６６】
（式４９）を次の（式５０）のように書き直す。
【０２６７】
【数６８】
2017057524.tif 000070
..............（式５０）
【０２６８】
両辺のフーリエ変換をとると、次の（式５１）が得られる。
【０２６９】
【数６９】
2017057524.tif 000071
.......（式５１）
【０２７０】
これを積分すると、以下の様になる。
【数７０】
2017057524.tif 000072
..........（式５２）
【０２７１】
ここで、ａ（ｋｘ，ｋｙ，ｋ）は、次の（式５３）のように得られる。
【０２７２】
【数７１】
2017057524.tif 000073
............（式５３）
【０２７３】
さらに、全てのＩに関して和を求めると、次の（式５４）が得られる。
【０２７４】
【数７２】
2017057524.tif 000074
......（式５４）
【０２７５】
偏微分方程式（式３７）の解は、次の（式５５）のようになる。
【０２７６】
【数７３】
2017057524.tif 000075
.........（式５５）
【０２７７】
これをさらに変形して、次の（式５６）が得られる。
【０２７８】
【数７４】
2017057524.tif 000076
..（式５６）
【０２７９】
各θで得られた結果を積分することにより、次の（式５７）で示される３次元再構成画像が得られる。
【０２８０】
【数７５】
2017057524.tif 000077
....（式５７）
【０２８１】
ここで、接空間で計算された結果を全体座標（Ｘ，Ｙ，Ｚ）へ変換する。ｙ軸の（Ｘ，Ｙ）平面への射影をｙ'とすると、次の（式５８）が成立する。
【０２８２】
【数７６】
2017057524.tif 000078
.........（式５８）
【０２８３】
さらに、（ｘ，ｙ'）から（Ｘ，Ｙ）へ変換する式は、次の（式５９）のようになる。
【０２８４】
【数７７】
2017057524.tif 000079
.............（式５９）
【０２８５】
この変換公式を、次の（式６０）へ適用していく。
【０２８６】
【数７８】
2017057524.tif 000080
...........（式６０）
【０２８７】
さらに、上述した式を整理すると、上の（式６０）は、次の（式６１）のようになる。
【０２８８】
【数７９】
2017057524.tif 000081
..............（式６１）
【０２８９】
さらに、スペクトラム領域で次の（式６２）に示す変数変換（逆変換）を行う。
【０２９０】
なお、以下の（式６２）および（式６３）で使用する変数ξ、η、ζは、図１０および（式３１）等で使用したＰ点の座標（ξ，η，ζ）ではなく、変数変換のために新たに用いた変数である。
【０２９１】
【数８０】
2017057524.tif 000082
..............（式６２）
【０２９２】
この逆変換は、次の（式６３）で与えられる。
【０２９３】
【数８１】
2017057524.tif 000083
..............（式６３）
【０２９４】
角度θにおける再構成関数は、次の（式６４）のようになる。
【０２９５】
【数８２】
2017057524.tif 000084
..............（式６４）
【０２９６】
ここで、被積分項に現れる関数は次の（式６５）のようになる。
【０２９７】
【数８３】
2017057524.tif 000085
..............（式６５）
【０２９８】
ここで、ｋｘ、ｋｙ、ｋｚは（式２８）のようにξ、η、ζの関数であるから、フーリエ変換だけを用いて、角度θの局所座標におけるデータを全体座標へ変換できる。全体座標へ変換するには次の（式６６）を用いる。つまり、最後に角度θで積分して、映像化関数が得られる。
【０２９９】
【数８４】
2017057524.tif 000086
.............（式６６）
【０３００】
この映像化関数により、対象物１０の内部情報に関する画像を再構成する。したがって、対象物１０の内部の欠陥等、対象物１０の内部情報に関する画像を汎用的かつ高速に映像化することができる。また、画像を再構成するステップにおいて、関数φの設定は、３次元空間を対象に行われるので、曲率の大きい曲面を有する対象物１０の内部情報をより精度よく高速に映像化することができる。
【０３０１】
以上、本実施の形態に係るマンモグラフィ装置は、曲面状にプローブ３６を移動する解析モデルにおいて逆問題用の偏微分方程式を設定し、これを解くことで、映像化関数を得ている。これにより、対象物である生体に放射した波動の散乱波を解析するマンモグラフィ装置において、曲率の大きい曲面を有する対象物の内部情報を、汎用的かつ高速に映像化することができる。
【０３０２】
（実施の形態３）
次に、実施の形態３について説明する。本実施の形態にかかるマンモグラフィ装置が実施の形態１に係るマンモグラフィ装置１と異なる点は、非対称の円錐状の形状を有する対象物を解析モデルとしている点である。
【０３０３】
実施の形態３に係るマンモグラフィ装置１では、送信アンテナ素子３６ａおよび受信アンテナ素子３６ｂは、一体となって対象物１０のうち非対称の錐体の母線に相当する線に沿って移動される。ここで、非対称の錐体の母線に相当する線とは、対象物１０の形状を非対称の錐体とみなしたとき、その錐体の母線に相当する線である。
【０３０４】
生体を観測対象とする場合、生体の形状は、実施の形態１および２に示した解析モデルのように、回転対称の形状を有する円錐体または略円錐体の形状に一部は当てはまるものの完全には当てはまらない場合がある。このような場合に、本実施の形態にかかるマンモグラフィ装置１における解析モデルを用いてもよい。
【０３０５】
以下、回転非対称の錐体を解析モデルとするマンモグラフィ装置について説明する。図１２は、本実施の形態に係る散乱トモグラフィ方法の原理を説明するための解析モデルである。
【０３０６】
本実施の形態にかかるマンモグラフィ装置の画像再構成アルゴリズムは、実施の形態１にかかるマンモグラフィ装置の画像再構成アルゴリズムとほぼ同様であるが、図２に示す解析モデルにおける角αが母線の回転角θに依存する点が実施の形態１にかかるマンモグラフィ装置１と異なっている。つまり、対象物１０は、回転非対称の錐体の形状を有しており、かつ、一部が欠けた形状をしている。
【０３０７】
以下、画像再構成アルゴリズムにおいて実施の形態１にかかるマンモグラフィ装置１と異なる点のみ説明する。
【０３０８】
実施の形態１に示した（式２６）について、スペクトラム領域で変数変換を行う際の変数変換公式とその逆変換公式は、次の（式６７）で表される。
【０３０９】
【数８５】
2017057524.tif 000087
.............（式６７）
【０３１０】
角度θにおける再構成関数は、次の（式６８）のようになる。
【０３１１】
【数８６】
2017057524.tif 000088
.............（式６８）
【０３１２】
ここで、変数ξ、η、ζの関数であるため、角度θの局所座標におけるデータを全体座標へフーリエ変換だけを用いて変換することができる。最後に、角度θで積分することにより、以下の（式６９）に示す再構成画像が得られる。
【０３１３】
【数８７】
2017057524.tif 000089
.............（式６９）
【０３１４】
以上、本実施の形態にかかるマンモグラフィ装置によると、観測対象である生体の形状が回転対称の形状を有する円錐体または略円錐体の形状に完全には当てはまらない場合であっても、逆問題用の偏微分方程式を設定し、これを解くことで、映像化関数を得ることができる。これにより、物体に放射した波動の散乱波を解析するマンモグラフィ装置において、対象物がフレキシブルであったり、完全な円錐体または略円錐体でない場合であっても、対象物の内部情報を、汎用的かつ高速に映像化することができる。
【０３１５】
なお、上記した実施の形態では、実施の形態１に示した錐体の解析モデルの一部が欠けた解析モデルを一例として説明したが、実施の形態２に示した錐体の解析モデルの一部が欠けた解析モデルとしてもよい。この場合も、上述した本実施の形態にかかるマンモグラフィ装置の解析モデルにおける画像再構成アルゴリズムと同様、（式）について、スペクトラム領域で変数変換を行う際の変数変換公式とその逆変換公式を変更することにより映像化関数を得ることができる。
【０３１６】
（変形例等）
以上、本発明に係る検査装置および検査方法として、マンモグラフィを例として、複数の実施の形態に基づいて説明したが、本発明は実施の形態に限定されるものではない。実施の形態に対して当業者が思いつく変形を施して得られる形態、および、複数の実施の形態における構成要素を任意に組み合わせて実現される別の形態も本発明に含まれる。
【０３１７】
例えば、上述した実施の形態では、マンモグラフィを例として説明しているため、対象物として生体、具体的には乳房を例に説明しているが、対象物は乳房に限らず、他の生体の構造であってもよい。また、生体に限らず、生体以外の物体構造であってもよい。例えば、円錐状のコンクリート支柱などの構造物であってもよい。この場合、対象物の形状に適合した解析モデルを適宜構築してもよい。
【０３１８】
また、本実施の形態では、波動としてマイクロ波を用いたが、マイクロ波に限らず、他の周波数帯域の電磁波または超音波であってもよい。また、本実施の形態では、マイクロ波を用いているため、所定の周波数を有する周期波を用いているが、波動として周期波に限らずパルス波を用いてもよい。
【０３１９】
また、上記の演算式、および、演算式の導出手順は、一例であって、別の演算式、および、別の導出手順が用いられてもよい。
【０３２０】
また、検査装置において、特定の処理部が実行する処理を別の処理部が実行してもよい。また、検査装置において処理を実行する順番が変更されてもよいし、複数の処理が並行して実行されてもよい。
【０３２１】
また、本発明の検査装置における対象物の内部情報の観測ステップは、コンピュータによって実行されることとしてもよい。そして、本発明は、散乱トモグラフィ方法に含まれるステップを、コンピュータに実行させるためのプログラムとして実現できる。さらに、本発明は、そのプログラムを記録したＣＤ−ＲＯＭ等の非一時的なコンピュータ読み取り可能な記録媒体として実現できる。
【０３２２】
また、検査装置に含まれる複数の構成要素は、集積回路であるＬＳＩとして実現されてもよい。これらの構成要素は、個別に１チップ化されてもよいし、一部または全てを含むように１チップ化されてもよい。ここでは、ＬＳＩとしたが、集積度の違いにより、ＩＣ（Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）、システムＬＳＩ、スーパーＬＳＩまたはウルトラＬＳＩと呼称されることもある。
【０３２３】
また、集積回路化の手法はＬＳＩに限るものではなく、専用回路または汎用プロセッサで実現してもよい。プログラムすることが可能なＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）、または、ＬＳＩ内部の回路セルの接続および設定を再構成可能なリコンフィギュラブル・プロセッサを利用してもよい。
【０３２４】
さらには、半導体技術の進歩または派生する別技術によりＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて、マンモグラフィ装置に含まれる構成要素の集積回路化を行ってもよい。
【０３２５】
以上、図面を参照して本発明の実施形態について述べたが、これらは本発明の例示であり、上記以外の様々な構成を採用することもできる。
【０３２６】
以下、参考形態の例を付記する。
１−１．  対象物に対して放射した波動の散乱波を解析して前記対象物の内部状態を検査する検査方法であって、
前記対象物に当接して前記対象物に波動を放射する送信アンテナ素子によって、前記波動を前記対象物に放射するステップと、
前記送信アンテナ素子と一体に設けられ、前記対象物に当接して、前記送信アンテナ素子から送信された前記波動が前記対象物において散乱した散乱波を受信する受信アンテナ素子によって、前記散乱波を受信するステップと、
前記受信アンテナ素子により受信した散乱波を示す散乱波データから、前記対象物の内部情報に関する画像を再構成するステップとを含み、
前記画像を再構成するステップにおいて、
前記対象物の内部情報に関する画像を再構成するための再構成関数を、前記対象物と同一の形状を有する３次元空間を対象としてあらかじめ設定し、
前記再構成関数の漸近式が満たす偏微分方程式を構築し、
前記散乱波データから、前記偏微分方程式を解くことにより得られる映像化関数を導出し、
前記映像化関数により、前記対象物の内部情報に関する画像を再構成する
検査方法。
１−２．  前記対象物は、円錐状の形状を有し、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物の円錐の母線に沿って移動される
１−１．に記載の検査方法。
１−３．  前記対象物は、曲線状の母線を有する略円錐状の形状を有し、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記曲線状の母線に沿って移動される
１−１．に記載の検査方法。
１−４．  前記対象物は、非対称の錐体の形状を有し、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記非対称の錐体の形状の母線に沿って移動する
１−１．に記載の検査方法。
１−５．  前記波動は、マイクロ波である
１−１．から１−４．のいずれか１つに記載の検査方法。
１−６．  前記波動は、パルス波または所定の周波数を有する周期波である
１−１．から１−４．のいずれか１つに記載の検査方法。
２−１．  対象物に対して放射した波動の散乱波を解析して前記対象物の内部状態を検査する検査装置であって、
前記対象物に当接し前記波動を前記対象物に放射する送信アンテナ素子と、
前記送信アンテナ素子と一体に設けられ、前記対象物に当接して前記送信アンテナ素子から放射された前記波動が前記対象物において散乱した散乱波を受信する受信アンテナ素子と、
前記受信アンテナ素子により受信した散乱波を示す散乱波データから、前記対象物の内部情報に関する画像を再構成する画像再構成部とを備え、
前記画像再構成部は、
前記対象物の内部情報に関する画像を再構成するための再構成関数を、前記対象物と同一の形状を有する３次元空間を対象としてあらかじめ設定し、
前記再構成関数の漸近式が満たす偏微分方程式を構築し、
前記散乱波データから、前記偏微分方程式を解くことにより得られる映像化関数を導出し、
前記映像化関数により、前記対象物の内部情報に関する画像を再構成する
検査装置。
２−２．  前記対象物は、円錐状の形状を有し、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物の円錐の母線に沿って移動される
２−１．に記載の検査装置。
２−３．  前記対象物は、曲線状の母線を有する略円錐状の形状を有し、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記曲線状の母線に沿って移動される
２−１．に記載の検査装置。
２−４．  前記対象物は、非対称の錐体の形状を有し、
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記非対称の錐体の形状の母線に沿って移動する
２−１．に記載の検査装置。
２−５．  前記波動は、マイクロ波である
２−１．から２−４．のいずれか１つに記載の検査装置。
２−６．  前記波動は、パルス波または所定の周波数を有する周期波である
２−１．から２−４．のいずれか１つに記載の検査装置。
【産業上の利用可能性】
【０３２７】
本発明に係る検査装置および検査方法は、フレキシブルな形状を有する対象物に対する検査に有用であり、例えば、マンモグラフィ装置等の医療機器に適用可能である。【０３２８】
この出願は、２０１５年９月２９日に出願された日本出願特願２０１５−１９２２１６号を基礎とする優先権を主張し、その開示の全てをここに取り込む。
【図１】
2017057524.tif 000104
【図２】
2017057524.tif 000105
【図３】
2017057524.tif 000106
【図４】
2017057524.tif 000107
【図５】
2017057524.tif 000108
【図６】
2017057524.tif 000109
【図７】
2017057524.tif 000110
【図８】
2017057524.tif 000111
【図９】
2017057524.tif 000112
【図１０】
2017057524.tif 000113
【図１１】
2017057524.tif 000114
【図１２】
2017057524.tif 000115
【手続補正書】
【提出日】平成29年2月23日(2017.2.23)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
波動を対象物に放射するステップと、
前記波動が前記対象物において散乱した散乱波を受信するステップと、
前記散乱波を示す散乱波データから、前記対象物の内部情報に関する画像を再構成するステップとを含み、
前記画像を再構成するステップにおいて、
前記散乱波データと形状を示す解析モデルとを用いて偏微分方程式を解くことにより、前記対象物の内部情報に関する画像を再構成するための再構成関数を導出し、
前記再構成関数の時間変数の極限値である映像化関数を導出し、前記映像化関数を用いて前記対象物の内部情報に関する画像を再構成し、
前記偏微分方程式は、
前記再構成関数が満たす方程式であり、
前記対象物内の互いに同一位置である送信点および受信点の位置を示す独立変数を用いて表され、
前記独立変数の数と同一の次元を持つ空間における各点において、前記散乱波が生じている場の関数である散乱場関数が解となる線形偏微分方程式であり、
前記放射するステップおよび前記受信するステップでは、前記対象物に前記波動を放射する送信アンテナ素子と前記散乱波を受信する受信アンテナ素子とが対をなして移動され、
モノスタティックな散乱トモグラフィを利用する画像化方法。
【請求項２】
（削除）
【請求項３】
前記解析モデルが示す形状は円錐であり、
前記映像化関数は、以下の（式Ａ）により表され、
前記画像を再構成するステップでは、前記映像化関数をθで積分することにより得られる以下の（式Ｂ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成する請求項１に記載の画像化方法。
【数１】
2017057524.tif 000116
【数２】
2017057524.tif 000117
2017057524.tif 000118
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記円錐の頂点を原点とし、前記原点から前記円錐の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、φＲはφ（ｘ，ｙ，０，ｋ）＝φＲ（ｙ，ｋ）δ（ｘ）が成り立つ関数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記円錐の軸を回転軸とする回転角度である。
【請求項４】
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち円錐の母線に相当する線に沿って移動される請求項３に記載の画像化方法。
【請求項５】
前記解析モデルが示す形状は曲線状の母線を有する略円錐であり、
前記映像化関数は、以下の（式Ｃ）により表され、
前記画像を再構成するステップでは、前記映像化関数をθで積分することにより得られる以下の（式Ｄ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成する請求項１に記載の画像化方法。
【数３】
2017057524.tif 000119
【数４】
2017057524.tif 000120
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記略円錐の頂点を原点とし、前記原点から前記略円錐の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、ａθは係数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記略円錐の軸を回転軸とする回転角度である。
【請求項６】
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち略円錐の母線に相当する線に沿って移動される請求項５に記載の画像化方法。
【請求項７】
前記解析モデルが示す形状は非対称の錐体であり、
前記映像化関数は、ξ、η、およびζがそれぞれ以下の（式Ｅ）を満たす変数であるとき以下の（式Ｆ）により表され、
前記画像を再構成するステップでは、前記映像化関数をθで積分することにより得られる以下の（式Ｇ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成する請求項１に記載の画像化方法。
【数５】
2017057524.tif 000121
【数６】
2017057524.tif 000122
【数７】
2017057524.tif 000123
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記錐体の頂点を原点とし、前記原点から前記錐体の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、φＲはφ（ｘ，ｙ，０，ｋ）＝φＲ（ｙ，ｋ）δ（ｘ）が成り立つ関数であり、ａθは係数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記錐体の軸を回転軸とする回転角度であり、αは前記Ｚ方向を基準とした傾き角度である。
【請求項８】
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち非対称の錐体の母線に相当する線に沿って移動される請求項７に記載の画像化方法。
【請求項９】
前記放射するステップでは、前記受信アンテナ素子と前記送信アンテナ素子とが一体に設けられ、前記受信アンテナ素子と前記送信アンテナ素子との間に電波吸収手段が位置しているプローブを用いて前記波動を前記対象物に放射し、
前記受信するステップでは、前記プローブを用いて前記散乱波を受信する請求項１および３〜８のいずれか一項に記載の画像化方法。
【請求項１０】
前記波動は、マイクロ波である
請求項１および３〜９のいずれか１項に記載の画像化方法。
【請求項１１】
前記波動は、パルス波または所定の周波数を有する周期波である
請求項１および３〜１０のいずれか１項に記載の画像化方法。
【請求項１２】
波動を対象物に放射する送信アンテナ素子と、
前記送信アンテナ素子から放射された前記波動が前記対象物において散乱した散乱波を受信する受信アンテナ素子と、
前記受信アンテナ素子により受信した前記散乱波を示す散乱波データから、前記対象物の内部情報に関する画像を再構成する画像再構成部とを備え、
前記画像再構成部は、
前記散乱波データと形状を示す解析モデルとを用いて偏微分方程式を解くことにより、前記対象物の内部情報に関する画像を再構成するための再構成関数を導出し、
前記再構成関数の時間変数の極限値である映像化関数を導出し、前記映像化関数を用いて前記対象物の内部情報に関する画像を再構成し、
前記偏微分方程式は、
前記再構成関数が満たす方程式であり、
前記対象物内の互いに同一位置である送信点および受信点の位置を示す独立変数を用いて表され、
前記独立変数の数と同一の次元を持つ空間における各点において、前記散乱波が生じている場の関数である散乱場関数が解となる線形偏微分方程式であり、
前記送信アンテナ素子と前記受信アンテナ素子とが対をなして移動され、
モノスタティックな散乱トモグラフィを利用する画像化装置。
【請求項１３】
（削除）
【請求項１４】
前記解析モデルが示す形状は円錐であり、
前記映像化関数は、以下の（式Ａ）により表され、
前記画像再構成部は、前記映像化関数をθで積分することにより得られる以下の（式Ｂ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成する請求項１２に記載の画像化装置。
【数８】
2017057524.tif 000124
【数９】
2017057524.tif 000125
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記円錐の頂点を原点とし、前記原点から前記円錐の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、φＲはφ（ｘ，ｙ，０，ｋ）＝φＲ（ｙ，ｋ）δ（ｘ）が成り立つ関数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記円錐の軸を回転軸とする回転角度である。
【請求項１５】
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち円錐の母線に相当する線に沿って移動される請求項１４に記載の画像化装置。
【請求項１６】
前記解析モデルが示す形状は曲線状の母線を有する略円錐であり、
前記映像化関数は、以下の（式Ｃ）により表され、
前記画像再構成部は、前記映像化関数をθで積分することにより得られる以下の（式Ｄ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成する請求項１２に記載の画像化装置。
【数１０】
2017057524.tif 000126
【数１１】
2017057524.tif 000127
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記略円錐の頂点を原点とし、前記原点から前記略円錐の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、ａθは係数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記略円錐の軸を回転軸とする回転角度である。
【請求項１７】
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち略円錐の曲線状の母線に相当する線に沿って移動される請求項１６に記載の画像化装置。
【請求項１８】
前記解析モデルが示す形状は非対称の錐体であり、
前記映像化関数は、ξ、η、およびζがそれぞれ以下の（式Ｅ）を満たす変数であるとき以下の（式Ｆ）により表され、
前記画像再構成部は、前記映像化関数をθで積分することにより得られる以下の（式Ｇ）により表される関数を用いて前記対象物の内部情報に関する画像を再構成する請求項１２に記載の画像化装置。
【数１２】
2017057524.tif 000128
【数１３】
2017057524.tif 000129
【数１４】
2017057524.tif 000130
ここで、ｘは前記送信点および前記受信点の位置のｘ座標であり、ｙは前記送信点および前記受信点の位置のｙ座標であり、ｚは前記送信点および前記受信点の位置のｚ座標であり、（Ｘ，Ｙ，Ｚ）は前記錐体の頂点を原点とし、前記原点から前記錐体の底面の中心に向かう方向をＺ方向とし、前記底面に平行で前記原点を通る方向をＸ方向とし、前記底面に平行で前記原点を通り前記Ｘ方向と直交する方向をＹ方向としたときの座標であり、ρは誘電率の関数であり、φは前記再構成関数であり、φＲはφ（ｘ，ｙ，０，ｋ）＝φＲ（ｙ，ｋ）δ（ｘ）が成り立つ関数であり、ａθは係数であり、ｋｘ、ｋｙ、ｋｚは波数のｘ、ｙ、ｚ成分であり、ｋは時間の波数であり、θは前記錐体の軸を回転軸とする回転角度であり、αは前記Ｚ方向を基準とした傾き角度である。
【請求項１９】
前記送信アンテナ素子および前記受信アンテナ素子は、一体となって前記対象物のうち非対称の錐体の母線に相当する線に沿って移動される請求項１８に記載の画像化装置。
【請求項２０】
前記受信アンテナ素子と前記送信アンテナ素子とが一体に設けられたプローブをさらに備え、
前記受信アンテナ素子と前記送信アンテナ素子との間には電波吸収手段が位置する請求項１２および１４〜１９のいずれか一項に記載の画像化装置。
【請求項２１】
前記波動は、マイクロ波である
請求項１２および１４〜２０のいずれか１項に記載の画像化装置。
【請求項２２】
前記波動は、パルス波または所定の周波数を有する周期波である
請求項１２および１４〜２１のいずれか１項に記載の画像化装置。
【手続補正書】
【提出日】平成30年4月12日(2018.4.12)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】０１９５
【補正方法】変更
【補正の内容】
【０１９５】
これに対し、本実施の形態にかかるマンモグラフィ装置１による観測では１〜１０ＧＨｚのＵＷＢ（ウルトラワイドバンド）マイクロ波を用いる。マイクロ波は、生体、特に、脂質等に対してマイクロ波の減衰が極めて小さいため、乳房中の欠陥組織の観測に有効である。また、装置構成がＭＲＩに比べると簡単であり、造影剤など他の物質も必要ないため、汎用的に用いることができる。さらに、本実施の形態にかかるマンモグラフィ装置１では、上記した逆散乱理論を用い、マイクロ波のモノスタティック散乱波データを基にして、乳房などの生体内部の３次元的構成要素を映像化するため、生体の内部情報（構成）を汎用的かつ高速に映像化することができる。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】０２１６
【補正方法】変更
【補正の内容】
【０２１６】
なお、本解析モデルでは、図１０に示すように、円錐体６０ｂの頂点を点Ｏ、円錐体６０ｂの底面の円の中心を点Ｏ'とする。点Ｏから点Ｏ'に向かう方向をＺ方向（Ｚ軸方向）とする。また、円錐体６０ｂの側面において、ある位置で円錐体６０ｂに接する平面を、接平面８０ｂとする。この接平面８０ｂがＺ軸と交差する点を含み、円錐体６０ｂの底面と平行な平面を基準面７０ｂとする。基準平面７０ｂにおける任意の一方向をＸ方向（Ｘ軸方向）、基準平面７０ｂにおけるＸ方向と垂直な方向をＹ方向（Ｙ軸方向）とする。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】０２２６
【補正方法】変更
【補正の内容】
【０２２６】
２．曲線状の母線を持つ円錐体を解析モデルとする、回転対称曲面上のモノスタティック逆散乱理論
以下、この微分作用素を求める方法について述べる。本実施の形態に係る解析モデルでは、図１０および図１１に示したように、実施の形態１の解析モデルにおける円錐体の母線が曲線状であり、この曲線状の母線に沿って、送信アンテナおよび受信アンテナが移動する。そして、アンテナアレイが、Ｚ軸を中心に回転移動もしてもよい。曲線状の母線の上ではｒ１、ｒ２のｘ、ｙ、ｚ座標は、ｒ１＝（ｘ,ｙ１,ｚ１）（＝（ｘ,ｙ,ｚ））、ｒ２＝（ｘ,ｙ２,ｚ２）（＝（ｘ,ｙ,ｚ））と表される。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】０２２８
【補正方法】変更
【補正の内容】
【０２２８】
図１１において、断面曲線Ｓ上の点ＰＩで放射され、点Ｐで反射した後、点ＰＩで受信される波動を考えると、関数
2017057524.tif 000131
は、次の（式３１）ように書くことができる。なお、ベクトルｒ１およびベクトルｒ２について、ｒ１＝ｒ２＝ｒ（ｘ，ｙ，ｚ）が成り立つ。
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】０２４４
【補正方法】変更
【補正の内容】
【０２４４】
ｚに関する偏微分をそれぞれＤｚと書くと以下の式を得る。
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】０２５４
【補正方法】変更
【補正の内容】
【０２５４】
ここで、積分方程式（式４３）の解を求める。曲面上の点ＰＩで測定された時系列データφ（ｘ,ｙＩ,ｚＩ,ｔ)をフーリエ変換した関数Φを次の（式４４）ように書くことができる。
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】０２５６
【補正方法】変更
【補正の内容】
【０２５６】
ここでｚＩは次の式を満たす。
【国際調査報告】
2017057524.tif 000132
2017057524.tif 000133