JP4073987B2 - 地中レーダ装置の出力表示方法 - Google Patents
地中レーダ装置の出力表示方法 Download PDFInfo
- Publication number
- JP4073987B2 JP4073987B2 JP31088497A JP31088497A JP4073987B2 JP 4073987 B2 JP4073987 B2 JP 4073987B2 JP 31088497 A JP31088497 A JP 31088497A JP 31088497 A JP31088497 A JP 31088497A JP 4073987 B2 JP4073987 B2 JP 4073987B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- frequency
- desired depth
- penetrating radar
- ground penetrating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Radar Systems Or Details Thereof (AREA)
Description
【発明の属する技術分野】
本発明は、地中レーダ装置の出力表示方法に関し、特に地中に電磁波を送信しその反射波を受信して埋設管などを探査する地中レーダ装置の出力表示方法に関するものである。
【0002】
【従来の技術】
従来、この種の地中レーダ装置の出力を表示する表示方法として、図11に示すように、地中レーダ装置の側線(走査線)を含む断面の映像、すなわちBモード映像が中心であった。
図12は地中レーダ装置の原理を示す説明図である。
電磁波を所定の周期で所定周波数帯域幅だけ周波数を掃引することにより周波数変調(FM)して地中に送信するFMCWレーダ装置の場合、ある時点における送信波と、そのとき受信した反射波の振幅の比と位相の差を抽出する。
【0003】
そして、得られた信号をフーリエ変換により時間領域の信号に変換し、送信から受信までの遅延時間とその振幅に基づいて、埋設物までの距離を求めるものとなっている。
実際には、図13に示すように、地中レーダ装置の測線に沿ったトラバース方向の各地点で、電磁波の送受信が行われ、横軸が送信から受信までの遅延時間を示し、縦軸が振幅を示す信号波形、すなわちAモード波形が得れる。
【0004】
このようにして得られた各地点のAモード波形を並べて、横軸をトラバース方向とし、縦軸に遅延時間すなわち地表面からの深さをとって示すことにより、地中レーダ装置の側線(走査線)を含む断面の映像、すなわちBモード映像が得られるものとなっていた。
しかし、平面的広がりのあるもの、例えば、埋設管のように地表面に対して平行に埋設されているものを探査対象物とする場合には、その埋設状況が把握し難く、図12に示すように、地表面に沿った映像、すなわちCモード映像のほうが埋設状況を把握しやすい。
【0005】
【発明が解決しようとする課題】
しかしながら、このような従来の地中レーダ装置の出力表示方法では、地中レーダ装置の側線を含む断面を示すBモード映像であるため、多数の測線から得られた双曲線状のBモード映像を水平に切り出した場合には、1つの埋設管が分割された不自然な映像となり、分解能が低下するという問題点があった。
また、得られた信号の振幅および位相情報をベクトル合成処理する合成開口法などの信号処理技術を用いた場合であっても、地中の伝搬速度などの非一様性のために、埋設物のフォーカス点が上下し、切り出した面上にフォーカスした点があるとは限らず、同一深さの連続したCモード映像が安定して得られないという問題点があった。
本発明はこのような課題を解決するためのものであり、地表面に沿った所定深さの鮮明な地中断面映像が得られる地中レーダ装置の出力表示方法を提供することを目的としている。
【0006】
【課題を解決するための手段】
このような目的を達成するために、本発明による地中レーダ装置の出力表示方法は、周波数領域の出力データをフーリエ変換処理することにより時間領域の時系列データに変換し、得られた時間領域の時系列データのうち所望深さより浅い深さに対応する時間領域のデータを除去することにより修正時系列データを生成し、得られた修正時系列データを逆フーリエ変換処理することにより複素スペクトラムデータに変換し、得られた複素スペクトラムデータのうち所望深さに対応する特定周波数(利用者が指定する周波数であり、通常、地中レーダ装置の周波数帯域の中心周波数を用いる)の複素スペクトラムデータを抽出し、抽出した特定周波数の複素スペクトラムデータを用いて合成開口処理を行うことにより所望深さのフォーカス像を生成し、地中レーダ装置の所定トラバース方向に沿った複数の測線ごとに得られた所望深さのフォーカス像を合成することにより、所望深さでの埋設状況を示す平面映像を表示するようにしたものである。
したがって、浅い部分のデータを含まない修正時系列データから得られた複素スペクトラムデータのうち特定周波数(利用者が指定する周波数であり、通常、地中レーダ装置の周波数帯域の中心周波数を用いる)の複素スペクトラムデータを用いて合成開口処理が行われて所望深さのフォーカス像が生成され、地中レーダ装置の所定トラバース方向に沿った複数の測線ごとに得られた所望深さのフォーカス像が合成されて、所望深さでの埋設状況を示す平面映像が表示される。
【0008】
【発明の実施の形態】
次に、本発明について図面を参照して説明する。
図1は本発明の一実施の形態である地中レーダ装置の出力表示方法による信号処理を示すフローチャート、図2は一般的なFMCW方式の地中レーダ装置を示すブロック図である。
まず、図2を参照して、本発明が適用されるFMWC方式の地中レーダ装置について説明する。
【0009】
図2において、操作部1からの所定操作入力に応じて、制御部8は送信信号発生部3を制御して、所定周期で例えばゼロから数GHzまで掃引することにより周波数変調(FM)した送信信号を発生させる。
この送信信号は、方向性結合器からなるカプラ4に入力されて、後述する参照信号と分岐され、送信アンテナ5から地中に送信される。
【0010】
地中に埋設管が存在した場合、送信された電磁波の一部が、その埋設管により地表面側に反射する。
この反射波は、受信アンテナ6で受信され、受信信号として混合器からなる検波回路7に入力される。
ここで、受信信号は、カプラ4から分岐された参照信号と混合され、両信号の振幅の積と位相の差のみが抽出される。
なお、振幅の積が抽出されるが、参照信号と送信信号の振幅が既知の場合には、それらの比すなわち反射率に換算できる。
【0011】
この場合、埋設管からの反射波は、送信アンテナから受信アンテナまでの伝搬経路が長く、その分だけ位相が遅れるため、掃引周波数に対しては交流成分として得られ、この交流成分の周波数が埋設管の深さに比例するものとなる。
検波回路7からの出力は、制御部8のIFアンプにより増幅され、信号処理部9に入力される。
【0012】
信号処理部9では、制御部8からの出力信号に対して、図1に示すような手順で各種信号処理が行われる。
これにより、地表から所望の深さにおけるCモード映像が生成され、画像表示部10の画面に表示出力される。
【0013】
次に、図1を参照して、本発明の動作として、信号処理部における信号処理について説明する。
まず、制御部8のIFアンプから出力された信号が取得され(ステップ11)、フーリエ変換処理により時系列データへ変換される(ステップ12)。
図3は制御部からの出力信号の例を示す波形図であり、横軸が埋設管の深さに対応する周波数を示し、縦軸が信号の振幅を示している。
【0014】
この信号がフーリエ変換処理されて、図4に示すように、横軸が送信から受信までの遅延時間に対応する時間軸に変換された時系列データ、すなわちAモード波形が生成される。
次に、この時系列データに対して、利用者により設定入力された所望の深さ(z方向座標値)に応じた時間軸ウインドウ処理が実施され、図5に示すような、修正時系列データが生成される(ステップ13)。
【0015】
ここでいう時間軸ウインドウ処理とは、各時系列データのうち、所望の深さ(z方向)より浅い部分のデータを削除する処理である。
後述の合成開口処理において、固定的な特定周波数を用いた場合に、焦点深度が深くなって鮮明な像が得られるが、その欠点として、地表面付近の物体からの強い反射信号が非フォーカスの状態で重なってしまい、結果として得られる像がぼやける。
【0016】
したがって、時間軸ウィンドウ処理により、地表面付近の物体からの強い反射信号が存在する部分、例えば所望の深さ(z方向)より浅い部分のデータを削除することにより、鮮明な像が得られる。
FMCW方式において、それぞれ深さL1 ,L2 なる地点に2つの物体が埋設されており、送信信号の反射率をそれぞれr1 ,r2 とすると、周波数ωにおける受信信号は、検波された後、数1のように表される。
【0017】
【数1】
【0018】
ここで、周波数ωを一定(すなわち固定周波数)とすると、これら1つの信号を分離することは不可能となる。
一方、周波数ωを掃引した場合は、時間の関数ω(t)=atで表されることから、この場合の受信信号は、数2のように、それぞれの深さL1 ,L2 に対応した2の異なった周波数の信号の和で表され、高域通過フィルタなどにより、分離可能となる。
【0019】
【数2】
【0020】
続いて、信号処理部9では、この修正時系列データを逆フーリエ変換処理することにより、図6に示すような、修正複素スペクトラムデータを生成する(ステップ14)。
そして、この修正複素スペクトラムデータのうち、所望の深さに対応する特定周波数における修正複素スペクトラムデータ(実部および虚部)を抽出する(ステップ15)。
【0021】
次に、この固定的な特定周波数成分を用いて合成開口処理を実施する(ステップ16)。
通常、波動関数φは、x,y,z,kの4つの変数を持つ関数で与えられるが、本発明において、前述の時間軸ウインドウ処理を行い、特定周波数成分のみを処理対象としているため、変数kは定数と考えることができる。
【0022】
したがって、波動関数φ(x,y,z,k)はφ(x,y,z)となる。
この波動関数が満たす方程式は、時間に依存する因子をe2ωt とすると、数3のようになる。
ここで、φ(x,y,z)に対してフーリエ展開を行うと、数4となる。
【0023】
【数3】
【0024】
【数4】
【0025】
この数4を方程式数3に代入して、数5が得られ、Qに対する一般解は、数6のように表される。
【0026】
【数5】
【0027】
【数6】
【0028】
ここで、数6からは2つの解が得られるが、1つはz軸の負方向に埋設管が存在するときの解で、もう1つは正方向に埋設管が存在するときの解となる。
したがって、2つの解のうち、実際に物理的な意味を有するのは、第1項であることから、波動関数φ(x,y,z)の一般解は、数7,8のようになる。
【0029】
【数7】
【0030】
【数8】
【0031】
さらに、z=0の測定値でaを置き換え、zの値を埋設管の深さz=z0 に合わせることにより、z=0で測定したデータφ(x,y,0)は、z=z0 にある埋設管のすぐ上で測定した場合と等価な値、すなわちφ(x,y,z0 )に変換される。
【0032】
このように、固定的な特定周波数成分を用いて合成開口処理を行うことにより、深い焦点深度が得られ、任意の測線(y=固定)では、図7に示すように、比較的大きな幅の深さで、焦点のあった埋設管のフォーカス像が得られ、所望の深さと埋設管の深さに多少差があっても鮮明な像が得られる。
【0033】
なお、Cモード映像を得たい所望の深さの特定周波数は、利用者により指定されるものであり、通常、地中レーダ装置の周波数帯域の中心周波数が用いられる。
また、深い位置を探査する場合には、比較的低めの周波数を用い、浅い位置を探査する場合には、比較的高めの周波数を用いることにより、地中による電磁波の減衰の影響を抑制できる。
【0034】
以上の処理を全測線で繰り返し実施し(ステップ17)、yが異なる各測線ごとに得られたデータφ(x,y,z0 )を合成して平面図を生成することにより、図8に示すような、深さz=z0 におけるx,y平面上の強度分布、すなわちCモード画像が得られる。
【0035】
したがって、Cモードのための計測を行う場合には、図9に示すように、トラバース方向(x軸方向)に並行した複数本の測線に沿って、各地点で計測を行い、得られたデータに対して図1に示す信号処理を行うことにより、図10に示すような、埋設状況が把握し易い、水平平面に沿ったCモード画像が表示出力される。
【0036】
なお、以上の説明では、探査対象物として、水道管やガス管などの埋設管を例にしたが、これに限定されるものではなく、他の探査対象物であっても、前述と同様の作用効果を奏する。
【0037】
【発明の効果】
以上説明したように、本発明は、時間領域の時系列データのうち所望深さより浅い深さに対応する時間領域のデータを除去することにより修正時系列データを生成し、得られた修正時系列データを逆フーリエ変換処理して得られた複素スペクトルデータのうち、所望深さに対応する特定周波数(利用者が指定する周波数であり、通常、地中レーダ装置の周波数帯域の中心周波数を用いる)の複素スペクトラムデータを用いて合成開口処理を行うことにより所望深さのフォーカス像を生成し、地中レーダ装置の所定トラバース方向に沿った複数の測線ごとに得られた所望深さのフォーカス像を合成して表示するようにしたので、従来のように、測線を含む断面を示すBモード映像と比較して、所望深さでの埋設状況を示すCモード映像を得ることができ、探査した地中の埋設状況を的確に把握できるとともに、焦点深度が深くなり、地中の伝搬速度などの非一様性のために埋設物のフォーカス点が上下しても、鮮明な画像が得られる。これに加え、地表面からの強い反射および所望深さより浅い地点にある物体からの強い反射信号による非フォーカス映像の重複を回避でき、さらに鮮明なCモード映像が得られる。
【図面の簡単な説明】
【図1】 本発明の一実施の形態による地中レーダ装置の出力表示方法を示すフローチャートである。
【図2】 FMWC方式の地中レーダ装置を示すブロック図である。
【図3】 レーダ出力データ例を示す説明図である。
【図4】 時系列データ例を示す説明図である。
【図5】 修正時系列データ例を示す説明図である。
【図6】 修正複素スペクトラムデータ例を示す説明図である。
【図7】 任意測線でのフォーカス像の例を示す説明図である。
【図8】 所定深さにおける埋設管像の例を示す説明図である。
【図9】 Cモードのための計測例を示す説明図である。
【図10】 本発明により得られるCモード映像例を示す説明図である。
【図11】 地中レーダ装置の出力表示方法を示す説明図である。
【図12】 地中レーダ装置の原理を示す説明図である。
【図13】 従来の出力表示例を示す説明図である。
【符号の説明】
1…操作部、2…エンコーダ、3…送信信号発生部、4…カプラ、5…送信アンテナ、6…受信アンテナ、7…検波回路、8…制御部、9…信号処理部、10…画像表示部。
Claims (1)
- 電磁波を所定の周期で所定周波数帯域幅だけ周波数を掃引することにより周波数変調して地中に送信し、その送信波と受信した反射波の振幅の比と位相の差を示す周波数領域の出力データを取得し、地表面に沿って走査することにより各地点で得られた複数の出力データに基づいて、地中埋設物の埋設状況を示す映像を出力する地中レーダ装置において、
周波数領域の出力データをフーリエ変換処理することにより時間領域の時系列データに変換し、
得られた時間領域の時系列データのうち所望深さより浅い深さに対応する時間領域のデータを除去することにより修正時系列データを生成し、
得られた修正時系列データを逆フーリエ変換処理することにより複素スペクトラムデータに変換し、
得られた複素スペクトラムデータのうち所望深さに対応する特定周波数の複素スペクトラムデータを抽出し、
抽出した特定周波数の複素スペクトラムデータを用いて合成開口処理を行うことにより所望深さのフォーカス像を生成し、
地中レーダ装置の所定トラバース方向に沿った複数の測線ごとに得られた所望深さのフォーカス像を合成することにより、所望深さでの埋設状況を示す平面映像を表示することを特徴とする地中レーダ装置の出力表示方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP31088497A JP4073987B2 (ja) | 1997-11-12 | 1997-11-12 | 地中レーダ装置の出力表示方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP31088497A JP4073987B2 (ja) | 1997-11-12 | 1997-11-12 | 地中レーダ装置の出力表示方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH11142513A JPH11142513A (ja) | 1999-05-28 |
JP4073987B2 true JP4073987B2 (ja) | 2008-04-09 |
Family
ID=18010551
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP31088497A Expired - Lifetime JP4073987B2 (ja) | 1997-11-12 | 1997-11-12 | 地中レーダ装置の出力表示方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4073987B2 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0122357D0 (en) * | 2001-09-15 | 2001-11-07 | Secr Defence | Sub-surface radar imaging |
JP5579360B2 (ja) * | 2007-08-24 | 2014-08-27 | 三菱電機株式会社 | レーダ装置 |
-
1997
- 1997-11-12 JP JP31088497A patent/JP4073987B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JPH11142513A (ja) | 1999-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6573855B1 (en) | Three-dimensional questing method, three-dimensional voxel data displaying method, and device therefor | |
US4896116A (en) | Pulse radar method and apparatus for detecting an object | |
JP4344050B2 (ja) | 媒体中の物体の探査方法 | |
EP1342079B1 (en) | Frequency domain processing of acoustic micro imaging signals | |
Olver et al. | FMCW radar for hidden object detection | |
KR101865016B1 (ko) | 지표투과레이다 자료 기반의 공동 검출 방법 및 공동 검출 장치 | |
CN104244838A (zh) | 用于超声波成像的方法和装置 | |
CA1258536A (en) | Method and apparatus for combining waveforms | |
JP5814556B2 (ja) | 信号処理装置 | |
KR20180092793A (ko) | 3d 지표투과 레이더 신호 처리 장치 및 방법 | |
JP4073987B2 (ja) | 地中レーダ装置の出力表示方法 | |
Al-Nuaimy et al. | Automatic detection of hyperbolic signatures in ground-penetrating radar data | |
JP2528148B2 (ja) | 地下埋設物の探知方法および装置 | |
CN112924964A (zh) | 用于步进频率雷达信号的数据融合方法、装置及存储介质 | |
JP6829656B2 (ja) | 超音波診断装置 | |
JP3407988B2 (ja) | 地中埋設物探査装置 | |
JPH0390856A (ja) | 超音波顕微鏡 | |
JP2643318B2 (ja) | 超音波による3次元物体形状認識方法 | |
JP2551952B2 (ja) | 不可視物体の探知方法 | |
JP2003107164A (ja) | 開口合成探査装置 | |
Kachanov et al. | Frequency separation of signals and interferences in noise-immune ultrasonic testing of articles manufactured from complexly structured materials | |
JPH04204080A (ja) | 探査レーダの信号処理装置 | |
JPH09113630A (ja) | 埋設物探査装置 | |
JPH07101235B2 (ja) | 物体探査装置 | |
JPH03165289A (ja) | 地中レーダ装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20041008 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20060831 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060905 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070522 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070723 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20080122 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080124 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110201 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120201 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140201 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140201 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150201 Year of fee payment: 7 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
EXPY | Cancellation because of completion of term |