JPH02126176A - 埋設物探査装置の深度方向分解能改善方法 - Google Patents
埋設物探査装置の深度方向分解能改善方法Info
- Publication number
- JPH02126176A JPH02126176A JP63278407A JP27840788A JPH02126176A JP H02126176 A JPH02126176 A JP H02126176A JP 63278407 A JP63278407 A JP 63278407A JP 27840788 A JP27840788 A JP 27840788A JP H02126176 A JPH02126176 A JP H02126176A
- Authority
- JP
- Japan
- Prior art keywords
- waveform
- reflected
- moving average
- buried object
- received
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims description 15
- 238000001514 detection method Methods 0.000 claims 1
- 238000002592 echocardiography Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 6
- 238000005070 sampling Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000003014 reinforcing effect Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
(産業上の利用分野)
本発明は、パルス状電磁波を使用して地中埋設物や、コ
ンクリート構造物内の鉄筋等を探査する技術に関し、特
に深度方向の分解能を向上する方法に係る。
ンクリート構造物内の鉄筋等を探査する技術に関し、特
に深度方向の分解能を向上する方法に係る。
(従来技術)
パルス状電磁波を使用したこの種の装置は、−般に以下
のように構成されている。
のように構成されている。
第1図に本体装置の構成ブロック図を示す。ここで、1
は送信アンテナ、2は受信アンテナ、3はパルス状電磁
波を作成するパルス送信器、4は受信器、5はトリガ回
路、6はサンプリング回路、7は370回路、8はA/
D変換器、9はCPU、10はメモリ、11はビデオメ
モリ、12はCRT等の表示装置、13は探査物体表面
、14は探査物体、15は探査対象物標である。
は送信アンテナ、2は受信アンテナ、3はパルス状電磁
波を作成するパルス送信器、4は受信器、5はトリガ回
路、6はサンプリング回路、7は370回路、8はA/
D変換器、9はCPU、10はメモリ、11はビデオメ
モリ、12はCRT等の表示装置、13は探査物体表面
、14は探査物体、15は探査対象物標である。
装置を動作させるには、トリガ回路5からの信号に同期
したパルス状電磁波をパルス送信器3で発生させ、探査
物体表面13上の送信アンテナlより探査物体14へ発
射する。探査物体14中に探査対象物標15が存在すれ
ば、それによって反射されたエコーが受信アンテナ2に
より受信機4に受信される。受信された反射波形はサン
プリング回路6により低周波へ変換され、STC回路7
で深度による電磁波減衰の補正をされた後、A/D変換
器8によりデジタル信号に変換される。デジタル信号と
なった受信波は、CPU9によりメモリ10に記憶され
、測定者の必要に応じた表示形式で、ビデオメモリ11
に書き込まれ、表示装力波形を、(b)に受信された反
射信号の振幅へ一時間軸−にへの一表示例を示す、16
は探査物体表面13からの反射波、17は探査対象物標
15からの反射波である。一般にこの種の装置は、至近
距離の物標を測定対象としているため、送信されるパル
ス状電磁波は(a)に示すように半値幅数ナノ秒程度で
ある。このため送信波の持つ周波数帯域は、数十メガヘ
ルツから数百メガヘルツ程度と広く分布している。また
同時に、深度方向に接近して存在する二つの探査対象物
標を分離して判断できる最少分解能力も、送信されるパ
ルス状電磁波の幅により決定される。
したパルス状電磁波をパルス送信器3で発生させ、探査
物体表面13上の送信アンテナlより探査物体14へ発
射する。探査物体14中に探査対象物標15が存在すれ
ば、それによって反射されたエコーが受信アンテナ2に
より受信機4に受信される。受信された反射波形はサン
プリング回路6により低周波へ変換され、STC回路7
で深度による電磁波減衰の補正をされた後、A/D変換
器8によりデジタル信号に変換される。デジタル信号と
なった受信波は、CPU9によりメモリ10に記憶され
、測定者の必要に応じた表示形式で、ビデオメモリ11
に書き込まれ、表示装力波形を、(b)に受信された反
射信号の振幅へ一時間軸−にへの一表示例を示す、16
は探査物体表面13からの反射波、17は探査対象物標
15からの反射波である。一般にこの種の装置は、至近
距離の物標を測定対象としているため、送信されるパル
ス状電磁波は(a)に示すように半値幅数ナノ秒程度で
ある。このため送信波の持つ周波数帯域は、数十メガヘ
ルツから数百メガヘルツ程度と広く分布している。また
同時に、深度方向に接近して存在する二つの探査対象物
標を分離して判断できる最少分解能力も、送信されるパ
ルス状電磁波の幅により決定される。
しかし、実際に受信される波形は、送信アンテナ1や受
信アンテナ2の特性や、探査物体14内部での減衰によ
り帯域制限されたものであり、fblのようにリンギン
グを伴って受信される。このリンギングによって送信波
形の幅が充分に狭いものであっても、実際の時間軸方向
最小分解能力は低下してしまう、つまり探査物体14内
部の深度方向に接近した2物標を分離、判断できない欠
点が生じてしまうのである。
信アンテナ2の特性や、探査物体14内部での減衰によ
り帯域制限されたものであり、fblのようにリンギン
グを伴って受信される。このリンギングによって送信波
形の幅が充分に狭いものであっても、実際の時間軸方向
最小分解能力は低下してしまう、つまり探査物体14内
部の深度方向に接近した2物標を分離、判断できない欠
点が生じてしまうのである。
従来よりこの欠点を解決するために、多くのリンギング
除去方法が考案されてきた。例えば、受信される反射波
形は、送信アンテナlからの送信波形と探査物体14内
部の反射係数を時間軸上で積和演算した結果であり、こ
れは周波数軸上で扱うと互いの積算となることを利用し
て、周波数帯域で除算を行い反射係数のみを取り出す方
法である。
除去方法が考案されてきた。例えば、受信される反射波
形は、送信アンテナlからの送信波形と探査物体14内
部の反射係数を時間軸上で積和演算した結果であり、こ
れは周波数軸上で扱うと互いの積算となることを利用し
て、周波数帯域で除算を行い反射係数のみを取り出す方
法である。
また一方では、時間軸」二の理想波形に対するフィルタ
特性を最小2乗法により算出する方法なども考案されて
いる。
特性を最小2乗法により算出する方法なども考案されて
いる。
(発明が解決しようとする課題)
しかし、いずれの方法においても、時間軸上から周波数
軸上への変換や、求めたフィルタと受信波形の積和演算
には膨大な演算時間を必要とするため、強力な演算能力
を持たないこの種の装置では実現が困難である。
軸上への変換や、求めたフィルタと受信波形の積和演算
には膨大な演算時間を必要とするため、強力な演算能力
を持たないこの種の装置では実現が困難である。
(課題を解決するための手段)
本発明はこれらの欠点を解決するため、受信波形とその
差分波形及び移動平均波形を使用して受信波形の処理を
行う。即ちサンプリングされた受信波形とそのサンプリ
ング信号間の差分波形、及び最少10点程度から最大で
も反射波形長程度の範囲での移動平均に対して、各々受
信波形とピーク値を同一にするために必要な値の175
倍から5倍程度の波形倍率と、最大でも反射波形長以内
の遅れ、または進み位相差を持たせて加算することによ
り、反射波形に付随するリンギングを除去するようにし
たもので以下実施例につき図面により詳細に説明する。
差分波形及び移動平均波形を使用して受信波形の処理を
行う。即ちサンプリングされた受信波形とそのサンプリ
ング信号間の差分波形、及び最少10点程度から最大で
も反射波形長程度の範囲での移動平均に対して、各々受
信波形とピーク値を同一にするために必要な値の175
倍から5倍程度の波形倍率と、最大でも反射波形長以内
の遅れ、または進み位相差を持たせて加算することによ
り、反射波形に付随するリンギングを除去するようにし
たもので以下実施例につき図面により詳細に説明する。
(実施例)
第2図山)に示したサンプリングされた受信波形より、
対象物標からの反射エコーのみを取り出した波形f (
t)を第3図talに、その各サンプリング点の変化量
を示すため下記の式で求めた差分波形r ’ (t)を
(blに、 f ’(t)=f(t)l(t+1) ・・・・・
・・・・(1)下記の式で求めた、N区間内での移動平
均波形F(t)をtelに示す。
対象物標からの反射エコーのみを取り出した波形f (
t)を第3図talに、その各サンプリング点の変化量
を示すため下記の式で求めた差分波形r ’ (t)を
(blに、 f ’(t)=f(t)l(t+1) ・・・・・
・・・・(1)下記の式で求めた、N区間内での移動平
均波形F(t)をtelに示す。
ここで、第3図に示した波形は、各々最大値で正規化し
て示している。
て示している。
つぎに、上記で求めた3種の波形を各々適度な係数倍し
、適度な位相差を持たせて加算すると(4)に示すよう
な波形G (t)を得ることができる。
、適度な位相差を持たせて加算すると(4)に示すよう
な波形G (t)を得ることができる。
G(t)=f(t)+α−f ’ (t+ r) +β
・FM (t+φ)・・・・・・(3) ここで、tはサンプリング時間、α、βは波形倍率、τ
、φは位相差、Nは移動平均範囲である。
・FM (t+φ)・・・・・・(3) ここで、tはサンプリング時間、α、βは波形倍率、τ
、φは位相差、Nは移動平均範囲である。
希望する波形Q (t)を算出する各係数の探索は、全
ての組合せについて一定の評価関数のもとで良否の度合
を判定しても良いし、例えば最急傾斜法等の最適化アル
ゴリズムを使用して高速に算出することも可能である。
ての組合せについて一定の評価関数のもとで良否の度合
を判定しても良いし、例えば最急傾斜法等の最適化アル
ゴリズムを使用して高速に算出することも可能である。
その結果求められた波形G (t)をfd+に示すが、
明らかに元受信波形(a)に比較して時間軸つまり探査
深度方向の最小分解能が向上している。この実施例の場
合、反射波形の極性判定性を重視した評価関数によって
各係数を探索したため(diの波形となっているが、他
の目的に応じた評価関数を使用して最適な他の変形倍率
、位相差を探索することも可能である。
明らかに元受信波形(a)に比較して時間軸つまり探査
深度方向の最小分解能が向上している。この実施例の場
合、反射波形の極性判定性を重視した評価関数によって
各係数を探索したため(diの波形となっているが、他
の目的に応じた評価関数を使用して最適な他の変形倍率
、位相差を探索することも可能である。
ここで、移動平均波形の計算範囲Nは最少10点程度か
ら最大でも反射波形製程度の範囲、差分波形、移動平均
波形の倍率係数α、βは、受信波形とピーク値を同一に
するために必要な値の175倍から5倍程度、同様に位
相差T、φは、最大遅れ、進みとも反射波形製程度の値
である。
ら最大でも反射波形製程度の範囲、差分波形、移動平均
波形の倍率係数α、βは、受信波形とピーク値を同一に
するために必要な値の175倍から5倍程度、同様に位
相差T、φは、最大遅れ、進みとも反射波形製程度の値
である。
本発明の他の実施例として、上記実施例のように差分波
形、移動平均波形のみに倍率係数を持たせるだけでなく
、 十 G(t)=α−f(t)+β−f’(t+r)+y−F
(tXφ) −(41上式のように、反射波形自体に
も倍率係数を設けることも当然可能である。これは、前
記の(3)式の演算と等価であるが、係数を三波形に設
けることでそれらの係数α、β、γを整数値としても圧
縮の度合があまり悪化しないことがわかっている。
形、移動平均波形のみに倍率係数を持たせるだけでなく
、 十 G(t)=α−f(t)+β−f’(t+r)+y−F
(tXφ) −(41上式のように、反射波形自体に
も倍率係数を設けることも当然可能である。これは、前
記の(3)式の演算と等価であるが、係数を三波形に設
けることでそれらの係数α、β、γを整数値としても圧
縮の度合があまり悪化しないことがわかっている。
そこで、受信波形の入力源となるA/D変換器8の出力
が一般に整数値であることから、全乗算が整数計算のみ
で実行できる特徴を持つ。このことは、大きな演算能力
を必要とする小数乗算を使用せず、より小さな演算能力
の装置上での実現性を高める事となる。
が一般に整数値であることから、全乗算が整数計算のみ
で実行できる特徴を持つ。このことは、大きな演算能力
を必要とする小数乗算を使用せず、より小さな演算能力
の装置上での実現性を高める事となる。
本発明の効果を示すために、第4図(alに反射係数値
、fblにその場合の受信波形、fc)に第3図(d)
を算出したものと同一の係数N、α、β、γ、τ。
、fblにその場合の受信波形、fc)に第3図(d)
を算出したものと同一の係数N、α、β、γ、τ。
φを使用して処理を行った結果を示す。ここで、第3図
(atは3個の反射物体の位置及び反射係数を示してい
るが、各反射物体の距離間隔18が第4図ta+で示し
た単独の受信波形の継続時間と比較して短いために、単
独の受信波形第3図(alと反射係数(alの積和演算
結果である受信波形中)では、三個の物体からの反射波
形19が一塊となってしまい、分離した三個の反射波形
の合成であるとは判断できない。しかし、本処理の実行
後telでは、三個の物体からの反射波形20が、反射
極性の正逆をも含めて正確に判断できる。
(atは3個の反射物体の位置及び反射係数を示してい
るが、各反射物体の距離間隔18が第4図ta+で示し
た単独の受信波形の継続時間と比較して短いために、単
独の受信波形第3図(alと反射係数(alの積和演算
結果である受信波形中)では、三個の物体からの反射波
形19が一塊となってしまい、分離した三個の反射波形
の合成であるとは判断できない。しかし、本処理の実行
後telでは、三個の物体からの反射波形20が、反射
極性の正逆をも含めて正確に判断できる。
さらに、リンギング除去効果は低下するが、受信波形に
よっては、対雑音性を高めるために、雑音の影響が比較
的大きい差分波形を用いないで実行することも可能であ
る。
よっては、対雑音性を高めるために、雑音の影響が比較
的大きい差分波形を用いないで実行することも可能であ
る。
また、以上の実施例記載は、ディジタルサンプリングさ
れた波形より、差分波形と移動平均波形を求めることで
実施されているが、アナログ信号に対して微分回路と積
分回路、遅延回路、増幅器、アナログ加算器等により実
現しても本発明の請求範囲を越えるものではない。
れた波形より、差分波形と移動平均波形を求めることで
実施されているが、アナログ信号に対して微分回路と積
分回路、遅延回路、増幅器、アナログ加算器等により実
現しても本発明の請求範囲を越えるものではない。
(発明の効果)
以上説明したように、本発明は、主として加算と常数の
積算を中心とした演算で受信波形のリンギング除去が行
えるため、従来同様の効果を得るために必要であった膨
大な演算時間を大幅に短縮し、より筒便小型な装置への
搭載も可能となり、この種の探査装置の深度方向の最小
分解能を向上させ、近接した物標に対しての誤った判定
を防止し、瞬時に適切な判定を下せる利点を持つ。
積算を中心とした演算で受信波形のリンギング除去が行
えるため、従来同様の効果を得るために必要であった膨
大な演算時間を大幅に短縮し、より筒便小型な装置への
搭載も可能となり、この種の探査装置の深度方向の最小
分解能を向上させ、近接した物標に対しての誤った判定
を防止し、瞬時に適切な判定を下せる利点を持つ。
第1図は−・船釣なパルス状電磁波を使用した本第3図
(alは第2図中)に示したサンプリングされた受信波
形より、対象物標からの反射エコーのみを取り出したも
の、山)はその各サンプリング点の変化量を示した差分
波形、(C1は適度な区間内での移動平均を示した移動
平均波形、(d)は(a)の反射波形に対し、■)の差
分波形、(C)の移動平均波形を各々適度な係数倍して
適度な位相差を持たせて加算した結果、第4図(δ)は
反射係数値、(blはその場合の受信波形、tc+は処
理を行った結果を示す。 l・・・送信アンテナ、2・・・受信アンテナ、3・・
・パルス状電磁波を作成するパルス送信器、4・・・受
信器、5・・・トリガ回路、6・・・サンプリング回路
、7−5 T C回路、8−A / D変換器、9−C
P U 。 10・・・メモリ、11・・・ビデオメモリ、12・・
・CRT等の表示装置、13・・・探査物体表面、14
・・・探査物体、15・・・探査対象物標、16・・・
探査物体表面13からの反射波、17・・・探査対象物
標15からの反射波形、18・・・各反射物体の距離間
隔、19・・・受信時の三個の物体からの反射波形、2
0・・・本処理の実行後の三個の物体からの反射波形。 第4図 (a) (c) 3図 (b) (d)
(alは第2図中)に示したサンプリングされた受信波
形より、対象物標からの反射エコーのみを取り出したも
の、山)はその各サンプリング点の変化量を示した差分
波形、(C1は適度な区間内での移動平均を示した移動
平均波形、(d)は(a)の反射波形に対し、■)の差
分波形、(C)の移動平均波形を各々適度な係数倍して
適度な位相差を持たせて加算した結果、第4図(δ)は
反射係数値、(blはその場合の受信波形、tc+は処
理を行った結果を示す。 l・・・送信アンテナ、2・・・受信アンテナ、3・・
・パルス状電磁波を作成するパルス送信器、4・・・受
信器、5・・・トリガ回路、6・・・サンプリング回路
、7−5 T C回路、8−A / D変換器、9−C
P U 。 10・・・メモリ、11・・・ビデオメモリ、12・・
・CRT等の表示装置、13・・・探査物体表面、14
・・・探査物体、15・・・探査対象物標、16・・・
探査物体表面13からの反射波、17・・・探査対象物
標15からの反射波形、18・・・各反射物体の距離間
隔、19・・・受信時の三個の物体からの反射波形、2
0・・・本処理の実行後の三個の物体からの反射波形。 第4図 (a) (c) 3図 (b) (d)
Claims (4)
- (1)パルス状電磁波を使用した埋設物探査装置におい
て、サンプリングされた受信波形と、そのサンプリング
信号間の差分波形、及び最小10点程度から最大でも反
射波形長程度の範囲での移動平均波形に対して、各々受
信波形とピーク値を同一にするために必要な値の1/5
倍から5倍程度の波形倍率と、最大でも反射波形長以内
の遅れ、または進み位相差を持たせて加算することによ
り、反射波形に付随するリンギングを除去するようにし
た埋設物探査装置の深度方向分解能改善方法。 - (2)前記波形倍率を受信波形に対しても決定し、各波
形倍率を整数値とした請求項(1)記載の埋設物探査装
置の深度方向分解能改善方法。 - (3)前記差分波形と移動平均波形を共にあるいは一方
を、差分波形を微分波形により、移動平均波形を積分波
形によってリンギングを除去した請求項(1)、(2)
記載の埋設物探査装置の深度方向分解能改善方法。 - (4)前記差分波形を使用せずにリンギングを除去した
請求項(1)、(2)記載の埋設物探査装置の深度方向
分解能改善方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP63278407A JPH02126176A (ja) | 1988-11-02 | 1988-11-02 | 埋設物探査装置の深度方向分解能改善方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP63278407A JPH02126176A (ja) | 1988-11-02 | 1988-11-02 | 埋設物探査装置の深度方向分解能改善方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JPH02126176A true JPH02126176A (ja) | 1990-05-15 |
Family
ID=17596918
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP63278407A Pending JPH02126176A (ja) | 1988-11-02 | 1988-11-02 | 埋設物探査装置の深度方向分解能改善方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH02126176A (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5363106A (en) * | 1992-03-27 | 1994-11-08 | Asahi Glass Company, Ltd. | Method for evaluating internal quality of fused cast refractories |
JP2021012168A (ja) * | 2019-07-09 | 2021-02-04 | オムロン株式会社 | 埋設物検出装置および埋設物検出方法 |
JP2021012154A (ja) * | 2019-07-09 | 2021-02-04 | オムロン株式会社 | 埋設物検出装置および埋設物検出方法 |
-
1988
- 1988-11-02 JP JP63278407A patent/JPH02126176A/ja active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5363106A (en) * | 1992-03-27 | 1994-11-08 | Asahi Glass Company, Ltd. | Method for evaluating internal quality of fused cast refractories |
JP2021012168A (ja) * | 2019-07-09 | 2021-02-04 | オムロン株式会社 | 埋設物検出装置および埋設物検出方法 |
JP2021012154A (ja) * | 2019-07-09 | 2021-02-04 | オムロン株式会社 | 埋設物検出装置および埋設物検出方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0211219B1 (en) | Apparatus for measuring ultrasonic charateristics of a medium | |
US4896116A (en) | Pulse radar method and apparatus for detecting an object | |
EP0557660B1 (en) | Signal processor | |
US5644314A (en) | Portable geophysical system using an inverse collocation-type metehodology | |
US4403314A (en) | Active detection system using simultaneous multiple transmissions | |
JP2005530164A (ja) | 物体検出用システムにおける妨害を抑制するための方法 | |
JPH07191125A (ja) | 送信器の位置決定方法及び装置 | |
US4222050A (en) | Moving target indication radar | |
EP0240201B1 (en) | Method and apparatus for combining waveforms | |
US5949739A (en) | Sonar bearing estimation of extended targets | |
JPH06294864A (ja) | レーダ装置 | |
JP2005128011A (ja) | パルス圧縮処理装置 | |
JPH02126176A (ja) | 埋設物探査装置の深度方向分解能改善方法 | |
Babur et al. | Design and implementation of cross-channel interference suppression for polarimetric LFM-CW radar | |
CA2063132A1 (en) | Echo ranging system | |
JP2002303645A (ja) | 周波数計測装置、周波数計測方法およびレーダ装置 | |
KR101524550B1 (ko) | 표적 속도에 따른 도플러 효과를 보상하는 고속 lfm 표적 검출 방법 및 장치 | |
WO2011058527A1 (en) | Method and apparatus for processing sonar signals | |
JP3528580B2 (ja) | 物体計測装置 | |
Shutko et al. | The Doppler Radar Signals Processing | |
JP2586958B2 (ja) | パルス波を用いた媒質内探査方法 | |
JP2697627B2 (ja) | 相関処理装置および合成開口レーダの画像再生処理装置 | |
JP2007033326A (ja) | 測角装置 | |
GB2110439A (en) | Digital data processing | |
JPS62165174A (ja) | レ−ダ装置 |