JP2000221266A - ３次元ボクセルデ―タ表示方法及び装置 - Google Patents
３次元ボクセルデ―タ表示方法及び装置Info
- Publication number
- JP2000221266A JP2000221266A JP11330966A JP33096699A JP2000221266A JP 2000221266 A JP2000221266 A JP 2000221266A JP 11330966 A JP11330966 A JP 11330966A JP 33096699 A JP33096699 A JP 33096699A JP 2000221266 A JP2000221266 A JP 2000221266A
- Authority
- JP
- Japan
- Prior art keywords
- absolute value
- maximum absolute
- voxel data
- dimensional voxel
- plane
- 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.)
- Granted
Links
Abstract
(57)【要約】
【課題】 媒質表面から媒質中へ放射した波動信号の
反射信号強度に基づいて媒質表面上の位置（ｘ，ｙ）と
反射時間ｔを座標（ｘ，ｙ，ｔ）として生成された３次
元ボクセルデータの表示方法及び装置に関し、ＳＮ比の
低い３次元ボクセルデータに対しても、簡単な２次元画
像表示で、媒質中に存在する物体等の把握を容易且つ的
確にできるようにする。 【解決手段】 ３次元ボクセルデータに対して、媒質表
面上の位置（ｘ，ｙ）毎に反射時間ｔ軸方向における振
幅値の最大絶対値を抽出する最大絶対値抽出工程と、最
大絶対値抽出工程で抽出された媒質表面上の位置（ｘ，
ｙ）毎の最大絶対値を所定平面上に配置する平面配置工
程と、所定平面上に配置された最大絶対値を表示する平
面表示工程とを実行する。
反射信号強度に基づいて媒質表面上の位置（ｘ，ｙ）と
反射時間ｔを座標（ｘ，ｙ，ｔ）として生成された３次
元ボクセルデータの表示方法及び装置に関し、ＳＮ比の
低い３次元ボクセルデータに対しても、簡単な２次元画
像表示で、媒質中に存在する物体等の把握を容易且つ的
確にできるようにする。 【解決手段】 ３次元ボクセルデータに対して、媒質表
面上の位置（ｘ，ｙ）毎に反射時間ｔ軸方向における振
幅値の最大絶対値を抽出する最大絶対値抽出工程と、最
大絶対値抽出工程で抽出された媒質表面上の位置（ｘ，
ｙ）毎の最大絶対値を所定平面上に配置する平面配置工
程と、所定平面上に配置された最大絶対値を表示する平
面表示工程とを実行する。
Description
【０００１】
【発明の属する技術分野】本発明は、媒質表面から媒質
中へ放射した波動信号の反射信号強度に基づいて媒質表
面上の位置（ｘ，ｙ）と反射時間ｔを座標（ｘ，ｙ，
ｔ）として生成された３次元ボクセルデータを表示する
３次元ボクセルデータ表示方法及び装置に関する。
中へ放射した波動信号の反射信号強度に基づいて媒質表
面上の位置（ｘ，ｙ）と反射時間ｔを座標（ｘ，ｙ，
ｔ）として生成された３次元ボクセルデータを表示する
３次元ボクセルデータ表示方法及び装置に関する。
【０００２】
【従来の技術】このような３次元ボクセルデータ表示
は、例えば、地中埋設物探査においては、地表面を移動
しながら、電磁波または音波による波動信号を地中へ放
射し、地中に存在する物体からの反射信号を受信し、そ
の受信した反射信号強度に対する地表面上の位置（ｘ，
ｙ）と反射時間ｔを座標（ｘ，ｙ，ｔ）とする３次元ボ
クセルデータを生成し、通常、かかる３次元ボクセルデ
ータを垂直あるいは水平断面に沿って得られる画像を多
数表示し、これらを比較検証するという手間並びに熟練
度を要する膨大な作業を行なっていた。例えば、図５に
示す埋設状況における３次元ボクセルデータを水平断面
に沿って表示した画像表示例を図１１に示す。この断面
表示の場合、一定の深さの断面しか表示できないため、
他の深さの埋設物の様子を把握するためには、深さを変
えて再表示しなければならない。
は、例えば、地中埋設物探査においては、地表面を移動
しながら、電磁波または音波による波動信号を地中へ放
射し、地中に存在する物体からの反射信号を受信し、そ
の受信した反射信号強度に対する地表面上の位置（ｘ，
ｙ）と反射時間ｔを座標（ｘ，ｙ，ｔ）とする３次元ボ
クセルデータを生成し、通常、かかる３次元ボクセルデ
ータを垂直あるいは水平断面に沿って得られる画像を多
数表示し、これらを比較検証するという手間並びに熟練
度を要する膨大な作業を行なっていた。例えば、図５に
示す埋設状況における３次元ボクセルデータを水平断面
に沿って表示した画像表示例を図１１に示す。この断面
表示の場合、一定の深さの断面しか表示できないため、
他の深さの埋設物の様子を把握するためには、深さを変
えて再表示しなければならない。
【０００３】
【発明が解決しようとする課題】上記した従来技術の問
題を解決すべく、第二の従来技術として、３次元データ
の全情報を如何に２次元画像として表示するかという手
法が、電子情報通信学会論文誌Ｄ（Vol.J71-D No.10 p
p.2002-2009）の「雪中レーダシステムのための３次元
データ表示方法」において提案されている。
題を解決すべく、第二の従来技術として、３次元データ
の全情報を如何に２次元画像として表示するかという手
法が、電子情報通信学会論文誌Ｄ（Vol.J71-D No.10 p
p.2002-2009）の「雪中レーダシステムのための３次元
データ表示方法」において提案されている。
【０００４】上記第二の従来技術の中の表面表示型とい
う手法を、図５に示す埋設状況における３次元ボクセル
データに適用した表示例を図１２〜図１４に示す。この
場合、適切な閾値を探る目的で閾値を、図１２、図１
３、図１４の順に高中低と変化させたものを３通り示し
ている。これより、閾値が高いときは埋設物からの反射
信号が途切れて断続的になり埋設物が判別しにくく、ま
た、閾値が低いときは、信号強度の弱い反射信号まで表
示できるものの、反射信号自体のＳＮ比が低いためノイ
ズ成分も同時に表示され視認性が悪化することが分か
る。
う手法を、図５に示す埋設状況における３次元ボクセル
データに適用した表示例を図１２〜図１４に示す。この
場合、適切な閾値を探る目的で閾値を、図１２、図１
３、図１４の順に高中低と変化させたものを３通り示し
ている。これより、閾値が高いときは埋設物からの反射
信号が途切れて断続的になり埋設物が判別しにくく、ま
た、閾値が低いときは、信号強度の弱い反射信号まで表
示できるものの、反射信号自体のＳＮ比が低いためノイ
ズ成分も同時に表示され視認性が悪化することが分か
る。
【０００５】次に、超音波診断等でも用いられている上
記第二の従来技術の中の積分型という手法を、同じく図
５に示す埋設状況における３次元ボクセルデータに対し
て反射時間ｔ軸方向に適用した表示例を図１５に示す。
この場合も、反射信号自体のＳＮ比が低く、埋設物の領
域が探査領域全体に対して比較的小さいため、積分効果
によりコントラストが低下してしまい、視認性が悪い。
他にも、差分和型、積和型等も提案されているが、これ
らの手法では隣接値との差分を用いるため、ＳＮ比の低
いデータの場合には細かいノイズが強調されてしまうと
いう問題がある。
記第二の従来技術の中の積分型という手法を、同じく図
５に示す埋設状況における３次元ボクセルデータに対し
て反射時間ｔ軸方向に適用した表示例を図１５に示す。
この場合も、反射信号自体のＳＮ比が低く、埋設物の領
域が探査領域全体に対して比較的小さいため、積分効果
によりコントラストが低下してしまい、視認性が悪い。
他にも、差分和型、積和型等も提案されているが、これ
らの手法では隣接値との差分を用いるため、ＳＮ比の低
いデータの場合には細かいノイズが強調されてしまうと
いう問題がある。
【０００６】本発明は、かかる実情に鑑みてなされたも
のであり、ＳＮ比の低い３次元ボクセルデータに対して
も、簡単な２次元画像表示で、媒質中に存在する物体等
の把握を容易且つ的確にできる３次元ボクセルデータ表
示方法及び装置を提供することを目的としている。
のであり、ＳＮ比の低い３次元ボクセルデータに対して
も、簡単な２次元画像表示で、媒質中に存在する物体等
の把握を容易且つ的確にできる３次元ボクセルデータ表
示方法及び装置を提供することを目的としている。
【０００７】
【課題を解決するための手段】この目的を達成するため
の本発明の第一の特徴構成は、特許請求の範囲の欄の請
求項１に記載の如く、本発明に係る３次元ボクセルデー
タ表示方法が、媒質表面から媒質中へ放射した波動信号
の反射信号強度に基づいて前記媒質表面上の位置（ｘ，
ｙ）と反射時間ｔを座標（ｘ，ｙ，ｔ）として生成され
た３次元ボクセルデータの表示方法であって、前記３次
元ボクセルデータに対して、前記媒質表面上の位置
（ｘ，ｙ）毎に反射時間ｔ軸方向における振幅値の最大
絶対値を抽出する最大絶対値抽出工程と、前記最大絶対
値抽出工程で抽出された前記媒質表面上の位置（ｘ，
ｙ）毎の前記最大絶対値を所定平面上に配置する平面配
置工程と、前記所定平面上に配置された前記最大絶対値
を表示する平面表示工程を実行する点にある。
の本発明の第一の特徴構成は、特許請求の範囲の欄の請
求項１に記載の如く、本発明に係る３次元ボクセルデー
タ表示方法が、媒質表面から媒質中へ放射した波動信号
の反射信号強度に基づいて前記媒質表面上の位置（ｘ，
ｙ）と反射時間ｔを座標（ｘ，ｙ，ｔ）として生成され
た３次元ボクセルデータの表示方法であって、前記３次
元ボクセルデータに対して、前記媒質表面上の位置
（ｘ，ｙ）毎に反射時間ｔ軸方向における振幅値の最大
絶対値を抽出する最大絶対値抽出工程と、前記最大絶対
値抽出工程で抽出された前記媒質表面上の位置（ｘ，
ｙ）毎の前記最大絶対値を所定平面上に配置する平面配
置工程と、前記所定平面上に配置された前記最大絶対値
を表示する平面表示工程を実行する点にある。
【０００８】ここで、ボクセルの振幅値とは、その座標
（ｘ，ｙ，ｔ）での反射信号強度或いは所定の信号処理
を行った場合はその処理後の信号強度を意味し、通常
は、装置構成等に依存して正負何れか一方またはその両
方の極性を取り得る。
（ｘ，ｙ，ｔ）での反射信号強度或いは所定の信号処理
を行った場合はその処理後の信号強度を意味し、通常
は、装置構成等に依存して正負何れか一方またはその両
方の極性を取り得る。
【０００９】同第二の特徴構成は、特許請求の範囲の欄
の請求項２に記載の如く、前記第一の特徴構成に加え
て、前記最大絶対値抽出工程で抽出された前記最大絶対
値の中から最大値と最小値を求め、その最大値と最小値
が夫々表示階調の上限及び下限となるように前記最大絶
対値を正規化する点にある。
の請求項２に記載の如く、前記第一の特徴構成に加え
て、前記最大絶対値抽出工程で抽出された前記最大絶対
値の中から最大値と最小値を求め、その最大値と最小値
が夫々表示階調の上限及び下限となるように前記最大絶
対値を正規化する点にある。
【００１０】同第三の特徴構成は、特許請求の範囲の欄
の請求項３に記載の如く、前記第一または第二の特徴構
成に加えて、前記最大絶対値抽出工程において、前記媒
質表面上の位置（ｘ，ｙ）毎に前記最大絶対値を抽出す
る際に、その最大絶対値をとる反射時間ｔＭＡＸ も併
せて抽出する点にある。
の請求項３に記載の如く、前記第一または第二の特徴構
成に加えて、前記最大絶対値抽出工程において、前記媒
質表面上の位置（ｘ，ｙ）毎に前記最大絶対値を抽出す
る際に、その最大絶対値をとる反射時間ｔＭＡＸ も併
せて抽出する点にある。
【００１１】同第四の特徴構成は、特許請求の範囲の欄
の請求項４に記載の如く、前記第一、第二または第三の
特徴構成に加えて、前記最大絶対値抽出工程において、
前記媒質表面上の位置（ｘ，ｙ）毎に前記最大絶対値を
抽出する際に、反射時間ｔの所定範囲内にある振幅値を
除外する点にある。
の請求項４に記載の如く、前記第一、第二または第三の
特徴構成に加えて、前記最大絶対値抽出工程において、
前記媒質表面上の位置（ｘ，ｙ）毎に前記最大絶対値を
抽出する際に、反射時間ｔの所定範囲内にある振幅値を
除外する点にある。
【００１２】同第五の特徴構成は、特許請求の範囲の欄
の請求項５に記載の如く、前記第一、第二、第三または
第四の特徴構成に加えて、前記最大絶対値抽出工程にお
いて、前記媒質表面上の位置（ｘ，ｙ）毎に前記最大絶
対値を抽出する際に、振幅値の極性が正または負の何れ
か一方のボクセルのみを対象とする点にある。
の請求項５に記載の如く、前記第一、第二、第三または
第四の特徴構成に加えて、前記最大絶対値抽出工程にお
いて、前記媒質表面上の位置（ｘ，ｙ）毎に前記最大絶
対値を抽出する際に、振幅値の極性が正または負の何れ
か一方のボクセルのみを対象とする点にある。
【００１３】同第六の特徴構成は、特許請求の範囲の欄
の請求項６に記載の如く、前記第一、第二、第三、第四
または第五の特徴構成に加えて、前記最大絶対値抽出工
程前に、前記反射信号強度からなる３次元ボクセルデー
タに対して合成開口処理若しくはマイグレーション処理
を施して、表示対象となる前記３次元ボクセルデータを
生成する点にある。
の請求項６に記載の如く、前記第一、第二、第三、第四
または第五の特徴構成に加えて、前記最大絶対値抽出工
程前に、前記反射信号強度からなる３次元ボクセルデー
タに対して合成開口処理若しくはマイグレーション処理
を施して、表示対象となる前記３次元ボクセルデータを
生成する点にある。
【００１４】同第七の特徴構成は、特許請求の範囲の欄
の請求項７に記載の如く、前記第六の特徴構成に加え
て、前記合成開口処理若しくはマイグレーション処理を
複数の媒質中の伝搬速度で実行して、前記伝搬速度毎に
前記３次元ボクセルデータを生成し、前記伝搬速度毎の
前記各３次元ボクセルデータに対して、前記最大絶対値
抽出工程と前記平面配置工程と前記平面表示工程を各別
に実行し、前記各伝搬速度毎の前記平面表示工程におけ
る前記３次元ボクセルデータの表示結果から、適切な伝
搬速度による処理結果を選択する点にある。
の請求項７に記載の如く、前記第六の特徴構成に加え
て、前記合成開口処理若しくはマイグレーション処理を
複数の媒質中の伝搬速度で実行して、前記伝搬速度毎に
前記３次元ボクセルデータを生成し、前記伝搬速度毎の
前記各３次元ボクセルデータに対して、前記最大絶対値
抽出工程と前記平面配置工程と前記平面表示工程を各別
に実行し、前記各伝搬速度毎の前記平面表示工程におけ
る前記３次元ボクセルデータの表示結果から、適切な伝
搬速度による処理結果を選択する点にある。
【００１５】同第八の特徴構成は、特許請求の範囲の欄
の請求項８に記載の如く、本発明に係る３次元ボクセル
データ表示装置が、媒質表面から媒質中へ放射した波動
信号の反射信号強度に基づいて前記媒質表面上の位置
（ｘ，ｙ）と反射時間ｔを座標（ｘ，ｙ，ｔ）として生
成された３次元ボクセルデータを表示する３次元ボクセ
ルデータ表示装置であって、前記３次元ボクセルデータ
に対して、前記媒質表面上の位置（ｘ，ｙ）毎に反射時
間ｔ軸方向における振幅値の最大絶対値を抽出する最大
絶対値抽出手段と、前記最大絶対値抽出手段で抽出され
た前記媒質表面上の位置（ｘ，ｙ）毎の前記最大絶対値
を所定平面上に配置する平面配置手段と、前記所定平面
上に配置された前記最大絶対値を表示する平面表示手段
とを備えてなる点にある。
の請求項８に記載の如く、本発明に係る３次元ボクセル
データ表示装置が、媒質表面から媒質中へ放射した波動
信号の反射信号強度に基づいて前記媒質表面上の位置
（ｘ，ｙ）と反射時間ｔを座標（ｘ，ｙ，ｔ）として生
成された３次元ボクセルデータを表示する３次元ボクセ
ルデータ表示装置であって、前記３次元ボクセルデータ
に対して、前記媒質表面上の位置（ｘ，ｙ）毎に反射時
間ｔ軸方向における振幅値の最大絶対値を抽出する最大
絶対値抽出手段と、前記最大絶対値抽出手段で抽出され
た前記媒質表面上の位置（ｘ，ｙ）毎の前記最大絶対値
を所定平面上に配置する平面配置手段と、前記所定平面
上に配置された前記最大絶対値を表示する平面表示手段
とを備えてなる点にある。
【００１６】この３次元ボクセルデータ表示装置は前記
第一から第七の特徴構成による夫々の３次元ボクセルデ
ータ表示方法に使用可能であり、その基本の作用効果は
上記第一の特徴構成による３次元ボクセルデータ表示方
法の作用効果と共通する。
第一から第七の特徴構成による夫々の３次元ボクセルデ
ータ表示方法に使用可能であり、その基本の作用効果は
上記第一の特徴構成による３次元ボクセルデータ表示方
法の作用効果と共通する。
【００１７】同第九の特徴構成は、特許請求の範囲の欄
の請求項９に記載の如く、本発明に係るコンピュータ読
み取り可能な記録媒体が、上記第一、第二、第三、第四
または第五の特徴構成の３次元ボクセルデータ表示方法
における前記最大絶対値抽出工程と前記平面配置工程を
１または２以上のコンピュータに実行させるためのプロ
グラムと、前記３次元ボクセルデータ表示方法における
前記平面表示工程或いは前記平面表示工程の内の前記最
大絶対値を被処理データとして前記コンピュータに設け
られた画像表示手段に出力する手順を前記コンピュータ
に実行させるためのプログラムとを少なくとも記録して
いる点にある。
の請求項９に記載の如く、本発明に係るコンピュータ読
み取り可能な記録媒体が、上記第一、第二、第三、第四
または第五の特徴構成の３次元ボクセルデータ表示方法
における前記最大絶対値抽出工程と前記平面配置工程を
１または２以上のコンピュータに実行させるためのプロ
グラムと、前記３次元ボクセルデータ表示方法における
前記平面表示工程或いは前記平面表示工程の内の前記最
大絶対値を被処理データとして前記コンピュータに設け
られた画像表示手段に出力する手順を前記コンピュータ
に実行させるためのプログラムとを少なくとも記録して
いる点にある。
【００１８】同第十の特徴構成は、特許請求の範囲の欄
の請求項１０に記載の如く、前記第九の特徴構成に加え
て、上記第六の特徴構成の３次元ボクセルデータ表示方
法における前記合成開口処理と前記マイグレーション処
理の少なくとも何れか一方を前記コンピュータに実行さ
せるためのプログラムを記録している点にある。
の請求項１０に記載の如く、前記第九の特徴構成に加え
て、上記第六の特徴構成の３次元ボクセルデータ表示方
法における前記合成開口処理と前記マイグレーション処
理の少なくとも何れか一方を前記コンピュータに実行さ
せるためのプログラムを記録している点にある。
【００１９】同第十一の特徴構成は、特許請求の範囲の
欄の請求項１１に記載の如く、前記第十の特徴構成に加
えて、上記第七の特徴構成の３次元ボクセルデータ表示
方法における前記合成開口処理若しくは前記マイグレー
ション処理を複数の媒質中の伝搬速度で実行し、前記伝
搬速度毎に前記３次元ボクセルデータを生成する処理を
前記コンピュータに実行させるためのプログラムと、前
記伝搬速度毎の前記各３次元ボクセルデータに対して、
前記最大絶対値抽出工程と前記平面配置工程と前記平面
表示工程或いは前記平面表示工程の内の前記最大絶対値
を被処理データとして前記コンピュータに設けられた画
像表示手段に出力する手順を各別に前記コンピュータに
実行させるためのプログラムとを記録している点にあ
る。
欄の請求項１１に記載の如く、前記第十の特徴構成に加
えて、上記第七の特徴構成の３次元ボクセルデータ表示
方法における前記合成開口処理若しくは前記マイグレー
ション処理を複数の媒質中の伝搬速度で実行し、前記伝
搬速度毎に前記３次元ボクセルデータを生成する処理を
前記コンピュータに実行させるためのプログラムと、前
記伝搬速度毎の前記各３次元ボクセルデータに対して、
前記最大絶対値抽出工程と前記平面配置工程と前記平面
表示工程或いは前記平面表示工程の内の前記最大絶対値
を被処理データとして前記コンピュータに設けられた画
像表示手段に出力する手順を各別に前記コンピュータに
実行させるためのプログラムとを記録している点にあ
る。
【００２０】以下に、これら特徴構成の作用並びに効果
について説明する。
について説明する。
【００２１】上記第一の特徴構成によれば、３次元ボク
セルデータに対して、各媒質表面上の位置（ｘ，ｙ）毎
に反射時間ｔ軸方向における振幅値の最大絶対値を抽出
し、各媒質表面上の位置（ｘ，ｙ）毎の前記最大絶対値
を所定平面上に配置して表示する。従って、３次元ボク
セルデータのＳＮ比が低く、媒質中の物体の領域が探査
領域全体に対して比較的小さくても、強い反射信号の領
域が局所的に存在すれば、高いコントラストで表示され
るため、物体の平面配置を高い視認性で容易に把握する
ことができるようになる。
セルデータに対して、各媒質表面上の位置（ｘ，ｙ）毎
に反射時間ｔ軸方向における振幅値の最大絶対値を抽出
し、各媒質表面上の位置（ｘ，ｙ）毎の前記最大絶対値
を所定平面上に配置して表示する。従って、３次元ボク
セルデータのＳＮ比が低く、媒質中の物体の領域が探査
領域全体に対して比較的小さくても、強い反射信号の領
域が局所的に存在すれば、高いコントラストで表示され
るため、物体の平面配置を高い視認性で容易に把握する
ことができるようになる。
【００２２】上記第二の特徴構成によれば、各媒質表面
上の位置（ｘ，ｙ）毎の前記最大絶対値を所定平面上に
配置して表示する際に、前記最大絶対値の中から最大値
と最小値を求め、その最大値と最小値が夫々表示階調の
上限及び下限となるように前記最大絶対値の正規化を行
なうため、更にコントラストを高め、視認性を向上する
ことができる。
上の位置（ｘ，ｙ）毎の前記最大絶対値を所定平面上に
配置して表示する際に、前記最大絶対値の中から最大値
と最小値を求め、その最大値と最小値が夫々表示階調の
上限及び下限となるように前記最大絶対値の正規化を行
なうため、更にコントラストを高め、視認性を向上する
ことができる。
【００２３】上記第三の特徴構成によれば、各媒質表面
上の位置（ｘ，ｙ）毎に反射時間ｔ軸方向における振幅
値の最大絶対値を抽出する際に、その最大絶対値をとる
反射時間ｔＭＡＸ も併せて抽出するため、強い信号を
反射する物体領域がどの反射時間ｔ（媒質表面からの距
離、地中埋設物探査の場合では埋設深さに相当。）に存
在するかを容易に把握することができるのである。
上の位置（ｘ，ｙ）毎に反射時間ｔ軸方向における振幅
値の最大絶対値を抽出する際に、その最大絶対値をとる
反射時間ｔＭＡＸ も併せて抽出するため、強い信号を
反射する物体領域がどの反射時間ｔ（媒質表面からの距
離、地中埋設物探査の場合では埋設深さに相当。）に存
在するかを容易に把握することができるのである。
【００２４】上記第四の特徴構成によれば、各媒質表面
上の位置（ｘ，ｙ）毎に反射時間ｔ軸方向における振幅
値の最大絶対値を抽出する際に、反射時間ｔの所定範囲
内にある振幅値を除外するため、地中探査の場合では地
表付近や地下水面等の媒質中の反射信号強度が極端に強
くなる領域を除外することができ、その結果、媒質中の
物体領域（地中埋設物探査の場合では埋設物領域）に対
する視認性を向上させることができる。
上の位置（ｘ，ｙ）毎に反射時間ｔ軸方向における振幅
値の最大絶対値を抽出する際に、反射時間ｔの所定範囲
内にある振幅値を除外するため、地中探査の場合では地
表付近や地下水面等の媒質中の反射信号強度が極端に強
くなる領域を除外することができ、その結果、媒質中の
物体領域（地中埋設物探査の場合では埋設物領域）に対
する視認性を向上させることができる。
【００２５】上記第五の特徴構成によれば、各媒質表面
上の位置（ｘ，ｙ）毎に反射時間ｔ軸方向における振幅
値の最大絶対値を抽出する際に、振幅値の極性が正また
は負の何れか一方のボクセルのみを対象とすることによ
って、媒質中の物体からの反射信号のＳＮ比の高い方の
極性の振幅値についてのみ抽出するため、媒質中の物体
領域に対する視認性を向上させることができる。これ
は、地中埋設物探査の場合、放射信号の埋設物での反射
係数が埋設物の材質によってその極性が異なることによ
るもので、例えば、金属管の場合は反射係数が負で、樹
脂管や空洞の場合は正の反射係数となる。
上の位置（ｘ，ｙ）毎に反射時間ｔ軸方向における振幅
値の最大絶対値を抽出する際に、振幅値の極性が正また
は負の何れか一方のボクセルのみを対象とすることによ
って、媒質中の物体からの反射信号のＳＮ比の高い方の
極性の振幅値についてのみ抽出するため、媒質中の物体
領域に対する視認性を向上させることができる。これ
は、地中埋設物探査の場合、放射信号の埋設物での反射
係数が埋設物の材質によってその極性が異なることによ
るもので、例えば、金属管の場合は反射係数が負で、樹
脂管や空洞の場合は正の反射係数となる。
【００２６】上記第六の特徴構成によれば、合成開口処
理若しくはマイグレーション処理を施すことによって、
媒質表面に平行な（ｘ，ｙ）面内の分解能を向上するこ
とができる。また、反射信号強度の生データから構成さ
れる原３次元ボクセルデータに、合成開口処理若しくは
マイグレーション処理を施して得られた３次元ボクセル
データは、物体の埋設状況等に近似できるデータ（深度
スケールに変換された情報）となるため、視認性が一層
向上するのである。
理若しくはマイグレーション処理を施すことによって、
媒質表面に平行な（ｘ，ｙ）面内の分解能を向上するこ
とができる。また、反射信号強度の生データから構成さ
れる原３次元ボクセルデータに、合成開口処理若しくは
マイグレーション処理を施して得られた３次元ボクセル
データは、物体の埋設状況等に近似できるデータ（深度
スケールに変換された情報）となるため、視認性が一層
向上するのである。
【００２７】上記第七の特徴構成によれば、複数の媒質
中の伝搬速度を用いて前記合成開口処理若しくはマイグ
レーション処理を実行して、伝搬速度毎に生成した３次
元ボクセルデータを平面表示工程において表示し、その
表示結果のフォーカスの良さを判断することにより、媒
質中の伝搬速度が未知な場合であっても、容易に適切な
伝搬速度による処理結果を選択することが可能となるの
で、媒質表面に平行な（ｘ，ｙ）面内における高い分解
能を確保することができるのである。
中の伝搬速度を用いて前記合成開口処理若しくはマイグ
レーション処理を実行して、伝搬速度毎に生成した３次
元ボクセルデータを平面表示工程において表示し、その
表示結果のフォーカスの良さを判断することにより、媒
質中の伝搬速度が未知な場合であっても、容易に適切な
伝搬速度による処理結果を選択することが可能となるの
で、媒質表面に平行な（ｘ，ｙ）面内における高い分解
能を確保することができるのである。
【００２８】上記第九の特徴構成によれば、前記記録媒
体に記録されたプログラムが一旦前記コンピュータの主
記憶領域にインストールされると、その主記憶領域上の
プログラムを実行することにより、前記コンピュータに
上記第一乃至第五の特徴構成の３次元ボクセルデータ表
示方法を実行させることができる。結果として、上記第
一乃至第五の特徴構成の３次元ボクセルデータ表示方法
と同じ作用効果が発揮される。
体に記録されたプログラムが一旦前記コンピュータの主
記憶領域にインストールされると、その主記憶領域上の
プログラムを実行することにより、前記コンピュータに
上記第一乃至第五の特徴構成の３次元ボクセルデータ表
示方法を実行させることができる。結果として、上記第
一乃至第五の特徴構成の３次元ボクセルデータ表示方法
と同じ作用効果が発揮される。
【００２９】また、上記第十または第十一の特徴構成に
よれば、前記記録媒体に記録されたプログラムが一旦前
記コンピュータの主記憶領域にインストールされると、
その主記憶領域上のプログラムを実行することにより、
前記コンピュータに上記第六または第七の特徴構成の３
次元ボクセルデータ表示方法を実行させることができ
る。結果として、上記第六または第七の特徴構成の３次
元ボクセルデータ表示方法と同じ作用効果が発揮され
る。
よれば、前記記録媒体に記録されたプログラムが一旦前
記コンピュータの主記憶領域にインストールされると、
その主記憶領域上のプログラムを実行することにより、
前記コンピュータに上記第六または第七の特徴構成の３
次元ボクセルデータ表示方法を実行させることができ
る。結果として、上記第六または第七の特徴構成の３次
元ボクセルデータ表示方法と同じ作用効果が発揮され
る。
【００３０】
【発明の実施の形態】以下に、本発明を地中埋設物の３
次元探査に適用した場合における実施の形態を、図面に
基づいて説明する。図１に示すように、本発明に係る３
次元ボクセルデータ表示装置を具備した３次元探査装置
は、送受信手段である送受信機１０と、送受信機１０で
得られた信号を処理するデータ解析装置２０とを、主な
機器として備えて構成されている。そして、本願にあっ
ては、データ解析装置２０における解析処理にその特徴
がある。
次元探査に適用した場合における実施の形態を、図面に
基づいて説明する。図１に示すように、本発明に係る３
次元ボクセルデータ表示装置を具備した３次元探査装置
は、送受信手段である送受信機１０と、送受信機１０で
得られた信号を処理するデータ解析装置２０とを、主な
機器として備えて構成されている。そして、本願にあっ
ては、データ解析装置２０における解析処理にその特徴
がある。
【００３１】図１に示すように、媒質である土壌１にガ
ス等の流体を配送する鋼管などの物体２が埋設されてお
り、送受信機１０とデータ解析装置２０を備えた探査装
置３が地表面を移動しながら、物体２の埋設位置を探査
する。その移動方向は、ｘ方向である。そして、本願の
ように３次元ボクセルデータを得る場合は、ｘ方向のデ
ータの収集を終了した後、ｙ方向（図１の表裏方向）に
所定量の移動を繰り返しながら、ｘ方向データを逐次、
収集する。尚、図１に示す物体２は、探査対象である物
体を模式的に例示したものであり、例えば、図５に示す
ような埋設状況における埋設管３０の一部を示すもので
ある。
ス等の流体を配送する鋼管などの物体２が埋設されてお
り、送受信機１０とデータ解析装置２０を備えた探査装
置３が地表面を移動しながら、物体２の埋設位置を探査
する。その移動方向は、ｘ方向である。そして、本願の
ように３次元ボクセルデータを得る場合は、ｘ方向のデ
ータの収集を終了した後、ｙ方向（図１の表裏方向）に
所定量の移動を繰り返しながら、ｘ方向データを逐次、
収集する。尚、図１に示す物体２は、探査対象である物
体を模式的に例示したものであり、例えば、図５に示す
ような埋設状況における埋設管３０の一部を示すもので
ある。
【００３２】送受信機１０は、例えば１００ＭＨｚ〜１
ＧＨｚの図２（ロ）（１）に例示する単発のパルス信号
を送信回路１３で発生し、送信アンテナ１１より電磁波
として土壌１に放射する。例えば、図２（イ）に例示す
るように物体２表面上を移動した場合、送信アンテナ１
１より放射された電磁波の中の土壌に入射した入射波４
は物体２表面で反射散乱し、その中の反射波５が受信ア
ンテナ１２で受信された後、受信回路１４において、図
２（ロ）（２）に例示するような受信信号として復調増
幅される（この図において単一の線が一定位置で時間差
をおいて受信される複数の受信信号群に対応する）。送
信アンテナ１１より放射され、受信アンテナ１２で受信
されるまでの時間差（これが実質上の反射時間）ｔは土
壌１の表面から物体２までの距離と土壌１の比誘電率ε
または電磁波の伝搬速度より一義的に決定される。図１
に示す場合にあっては、送信アンテナ１１と受信アンテ
ナ１２は一定間隔で地表面に対向して配置される。ｘ方
向移動は、物体２を横切るように行われることとなる。
ＧＨｚの図２（ロ）（１）に例示する単発のパルス信号
を送信回路１３で発生し、送信アンテナ１１より電磁波
として土壌１に放射する。例えば、図２（イ）に例示す
るように物体２表面上を移動した場合、送信アンテナ１
１より放射された電磁波の中の土壌に入射した入射波４
は物体２表面で反射散乱し、その中の反射波５が受信ア
ンテナ１２で受信された後、受信回路１４において、図
２（ロ）（２）に例示するような受信信号として復調増
幅される（この図において単一の線が一定位置で時間差
をおいて受信される複数の受信信号群に対応する）。送
信アンテナ１１より放射され、受信アンテナ１２で受信
されるまでの時間差（これが実質上の反射時間）ｔは土
壌１の表面から物体２までの距離と土壌１の比誘電率ε
または電磁波の伝搬速度より一義的に決定される。図１
に示す場合にあっては、送信アンテナ１１と受信アンテ
ナ１２は一定間隔で地表面に対向して配置される。ｘ方
向移動は、物体２を横切るように行われることとなる。
【００３３】図１に示すように、送受信機１０には、受
信回路１４の増幅部の利得を時間差ｔに応じて変調する
信号強度変調手段１５が設けられており、時間差ｔが長
くなるにつれて土壌１を伝搬するパルス信号の損失が大
きくなり、受信信号強度（反射信号強度と同じ）が減衰
するのを振幅補正し、時間差ｔ、つまりは反射時間ｔの
増加に対して急激に減衰しない受信信号強度分布を得る
構成とされている。この構成により、後の信号処理に必
要な信号強度を確保できる。
信回路１４の増幅部の利得を時間差ｔに応じて変調する
信号強度変調手段１５が設けられており、時間差ｔが長
くなるにつれて土壌１を伝搬するパルス信号の損失が大
きくなり、受信信号強度（反射信号強度と同じ）が減衰
するのを振幅補正し、時間差ｔ、つまりは反射時間ｔの
増加に対して急激に減衰しない受信信号強度分布を得る
構成とされている。この構成により、後の信号処理に必
要な信号強度を確保できる。
【００３４】次に、受信信号が送られるデータ解析装置
２０について、図１及び図３に基づいて説明する。デー
タ解析装置２０は、マイクロコンピュータや半導体メモ
リ等によって構成されるデータ処理部２１と、外部から
の操作指示を入力するマウスやキーボード等の入力部２
２と、各処理段階での画像データや出力結果を表示する
ＣＲＴモニタや液晶ディスプレイ等の表示部２３を備え
て構成されている。更に、各処理段階でのデータや出力
結果等を保管格納する磁気ディスク等の外部補助記憶部
２４を備えている。
２０について、図１及び図３に基づいて説明する。デー
タ解析装置２０は、マイクロコンピュータや半導体メモ
リ等によって構成されるデータ処理部２１と、外部から
の操作指示を入力するマウスやキーボード等の入力部２
２と、各処理段階での画像データや出力結果を表示する
ＣＲＴモニタや液晶ディスプレイ等の表示部２３を備え
て構成されている。更に、各処理段階でのデータや出力
結果等を保管格納する磁気ディスク等の外部補助記憶部
２４を備えている。
【００３５】図３に示すように、このデータ処理部２１
は、受信回路１４から入力してくる受信信号を、媒質表
面上の位置（ｘ，ｙ）と時間ｔとの関係において整理、
処理する３次元ボクセルデータ生成手段２５を備えてい
る。この３次元ボクセルデータ生成手段２５は、以降の
処理で使用される３次元ボクセルデータを生成するため
のものであり、受信信号強度ｓをそのまま媒質表面上の
位置（ｘ，ｙ）と時間ｔとの関数とする原３次元ボクセ
ルデータｓ（ｘ，ｙ，ｔ）を生成する。更に、必要な場
合、この原３次元ボクセルデータｓ（ｘ，ｙ，ｔ）は、
マイグレーション処理されて、マイグレーション処理済
の新たな３次元ボクセルデータＳ（ｘ，ｙ，ｔ）が生成
される。これらのデータは３次元構造を取るため、ｓ
（ｘ，ｙ，ｔ）、Ｓ（ｘ，ｙ，ｔ）は共に、３次元ボク
セルデータに属し、３次元ボクセルデータ生成手段２５
で生成されるデータ群である。ここで、マイグレーショ
ン処理とは、媒質の表面において得られる移動方向情報
（空間、深度＝０、時間の情報）を、波の伝播を代表す
る波動方程式に基づいて、フーリエ、逆フーリエ変換手
法を利用して媒質の深度方向の情報（空間、深度、時間
＝０の情報）に変換する公知の手法である。この手法
は、空間座標としてｘ方向のみを対象とする場合、ｘが
媒質表面上の観測ライン、ｚが媒質中に向けて正の方向
をとる深度、ｔが伝搬時間である場合、波動の場をｕ
（ｘ，ｚ，ｔ）で表し、ｕ（ｘ，ｚ，ｔ）を３次元フー
リエ変換したものをＵ（ξ，η，ω）とすると、マイグ
レーション法はｔ＞０に対して得られた観測データｕ
（ｘ，０，ｔ）（レーダ画像）から時刻ｔ＝０における
深さ方向の場ｕ（ｘ，ｚ，０）（深度断面）を求めるも
のである。即ち、マイグレーション法の一種であるフェ
ーズ・シフト法では、以下のような処理を行う。 １．観測されたデータをｘとｔに関して２次元フーリエ
変換しＵ（ξ，０，ω）を求める。 ２．求めたい深度の１ラインを、Ｕ（ξ，０，ω）から
求める。 ３．２．の計算を深度を更新しながら繰り返し、深度断
面全体を求める。 一方、マイグレーション法の一種であるＦ−Ｋマイグレ
ーション法では、以下のような処理を行う。 １．観測されたデータをｘとｔに関して２次元フーリエ
変換しＵ（ξ，０，ω）を求める。 ２．周波数領域上で深度断面のフーリエ変換の値を求め
る。 ３．この値を、ξ，ηに関して２次元逆フーリエ変換
し、ｕ（ｘ，ｚ，０）を得る。 このようにして、ｔ＝０における断面構造を得ることが
できる。ここでは、空間座標としてｘ方向のみに関する
説明をしたが、本願のようにｘ、ｙ方向を共に対象とす
る場合も同様に取扱うことができる。
は、受信回路１４から入力してくる受信信号を、媒質表
面上の位置（ｘ，ｙ）と時間ｔとの関係において整理、
処理する３次元ボクセルデータ生成手段２５を備えてい
る。この３次元ボクセルデータ生成手段２５は、以降の
処理で使用される３次元ボクセルデータを生成するため
のものであり、受信信号強度ｓをそのまま媒質表面上の
位置（ｘ，ｙ）と時間ｔとの関数とする原３次元ボクセ
ルデータｓ（ｘ，ｙ，ｔ）を生成する。更に、必要な場
合、この原３次元ボクセルデータｓ（ｘ，ｙ，ｔ）は、
マイグレーション処理されて、マイグレーション処理済
の新たな３次元ボクセルデータＳ（ｘ，ｙ，ｔ）が生成
される。これらのデータは３次元構造を取るため、ｓ
（ｘ，ｙ，ｔ）、Ｓ（ｘ，ｙ，ｔ）は共に、３次元ボク
セルデータに属し、３次元ボクセルデータ生成手段２５
で生成されるデータ群である。ここで、マイグレーショ
ン処理とは、媒質の表面において得られる移動方向情報
（空間、深度＝０、時間の情報）を、波の伝播を代表す
る波動方程式に基づいて、フーリエ、逆フーリエ変換手
法を利用して媒質の深度方向の情報（空間、深度、時間
＝０の情報）に変換する公知の手法である。この手法
は、空間座標としてｘ方向のみを対象とする場合、ｘが
媒質表面上の観測ライン、ｚが媒質中に向けて正の方向
をとる深度、ｔが伝搬時間である場合、波動の場をｕ
（ｘ，ｚ，ｔ）で表し、ｕ（ｘ，ｚ，ｔ）を３次元フー
リエ変換したものをＵ（ξ，η，ω）とすると、マイグ
レーション法はｔ＞０に対して得られた観測データｕ
（ｘ，０，ｔ）（レーダ画像）から時刻ｔ＝０における
深さ方向の場ｕ（ｘ，ｚ，０）（深度断面）を求めるも
のである。即ち、マイグレーション法の一種であるフェ
ーズ・シフト法では、以下のような処理を行う。 １．観測されたデータをｘとｔに関して２次元フーリエ
変換しＵ（ξ，０，ω）を求める。 ２．求めたい深度の１ラインを、Ｕ（ξ，０，ω）から
求める。 ３．２．の計算を深度を更新しながら繰り返し、深度断
面全体を求める。 一方、マイグレーション法の一種であるＦ−Ｋマイグレ
ーション法では、以下のような処理を行う。 １．観測されたデータをｘとｔに関して２次元フーリエ
変換しＵ（ξ，０，ω）を求める。 ２．周波数領域上で深度断面のフーリエ変換の値を求め
る。 ３．この値を、ξ，ηに関して２次元逆フーリエ変換
し、ｕ（ｘ，ｚ，０）を得る。 このようにして、ｔ＝０における断面構造を得ることが
できる。ここでは、空間座標としてｘ方向のみに関する
説明をしたが、本願のようにｘ、ｙ方向を共に対象とす
る場合も同様に取扱うことができる。
【００３６】この目的から、データ処理部２１は、３次
元ボクセルの状態で得られているデータを、マイグレー
ション処理できるマイグレーション処理手段２６を備え
ている。更に、マイグレーション処理の代わりに、同様
に公知の方法である合成開口処理を施すものとしてもよ
い。この場合、マイグレーション処理手段２６の代わり
に合成開口処理手段が備えられることとなる。この状況
を図に括弧書きで示した。
元ボクセルの状態で得られているデータを、マイグレー
ション処理できるマイグレーション処理手段２６を備え
ている。更に、マイグレーション処理の代わりに、同様
に公知の方法である合成開口処理を施すものとしてもよ
い。この場合、マイグレーション処理手段２６の代わり
に合成開口処理手段が備えられることとなる。この状況
を図に括弧書きで示した。
【００３７】更に、このデータ処理部２１は、前記３次
元ボクセルデータ生成手段２５により生成された３次元
ボクセルデータｓ（ｘ，ｙ，ｔ）または前記マイグレー
ション処理手段２６により生成された３次元ボクセルデ
ータＳ（ｘ，ｙ，ｔ）に対して、前記媒質表面上の位置
（ｘ，ｙ）毎に反射時間ｔ軸方向における振幅値の最大
絶対値を抽出する最大絶対値抽出手段２７と、前記最大
絶対値抽出手段２７で抽出された前記媒質表面上の位置
（ｘ，ｙ）毎の前記最大絶対値を所定平面上に配置する
平面配置手段２８と、その平面配置された前記最大絶対
値を出力処理して前記表示部２３に表示する平面表示手
段２９とを備えている。
元ボクセルデータ生成手段２５により生成された３次元
ボクセルデータｓ（ｘ，ｙ，ｔ）または前記マイグレー
ション処理手段２６により生成された３次元ボクセルデ
ータＳ（ｘ，ｙ，ｔ）に対して、前記媒質表面上の位置
（ｘ，ｙ）毎に反射時間ｔ軸方向における振幅値の最大
絶対値を抽出する最大絶対値抽出手段２７と、前記最大
絶対値抽出手段２７で抽出された前記媒質表面上の位置
（ｘ，ｙ）毎の前記最大絶対値を所定平面上に配置する
平面配置手段２８と、その平面配置された前記最大絶対
値を出力処理して前記表示部２３に表示する平面表示手
段２９とを備えている。
【００３８】ここで、前記データ処理部２１を構成する
前記３次元ボクセルデータ生成手段２５、前記マイグレ
ーション処理手段２６或いは合成開口処理手段、前記最
大絶対値抽出手段２７、前記平面配置手段２８、前記平
面表示手段２９は、上述したマイクロコンピュータや半
導体メモリ等によって構成されるコンピュータを前記各
手段として機能させるためのプログラムが、前記コンピ
ュータの主記憶領域にインストールされ実行可能な状態
におかれることにより実現される。尚、主記憶領域と
は、前記プログラムを前記コンピュータの中央演算処理
装置がそのプログラムを逐次実行可能な状態に記録して
いる所定の記憶装置内の領域をいい、通常の記憶装置で
は、中央演算処理装置が高速にそのプログラムを読み出
せる半導体メモリとその実行プログラムを保存する不揮
発性の磁気ディスクメモリとの階層構造となっている
が、特に主記憶領域を構成する記憶装置が単体で高速ア
クセス性と不揮発性を兼ね備えている場合は、当該階層
構造は必ずしも必要ない。
前記３次元ボクセルデータ生成手段２５、前記マイグレ
ーション処理手段２６或いは合成開口処理手段、前記最
大絶対値抽出手段２７、前記平面配置手段２８、前記平
面表示手段２９は、上述したマイクロコンピュータや半
導体メモリ等によって構成されるコンピュータを前記各
手段として機能させるためのプログラムが、前記コンピ
ュータの主記憶領域にインストールされ実行可能な状態
におかれることにより実現される。尚、主記憶領域と
は、前記プログラムを前記コンピュータの中央演算処理
装置がそのプログラムを逐次実行可能な状態に記録して
いる所定の記憶装置内の領域をいい、通常の記憶装置で
は、中央演算処理装置が高速にそのプログラムを読み出
せる半導体メモリとその実行プログラムを保存する不揮
発性の磁気ディスクメモリとの階層構造となっている
が、特に主記憶領域を構成する記憶装置が単体で高速ア
クセス性と不揮発性を兼ね備えている場合は、当該階層
構造は必ずしも必要ない。
【００３９】次に、本発明に係る３次元ボクセルデータ
表示方法の一実施形態を、図４に示すデータ処理部２１
での典型的なデータ処理手順のフローチャートに基づい
て説明する。
表示方法の一実施形態を、図４に示すデータ処理部２１
での典型的なデータ処理手順のフローチャートに基づい
て説明する。
【００４０】１．３次元ボクセルデータ生成工程（ＳＴ
１） この工程は、３次元探査装置３を移動しながら、データ
を収集し、これを処理して、以降の処理に使用される３
次元ボクセルデータｓ（ｘ，ｙ，ｔ）若しくはＳ（ｘ，
ｙ，ｔ）を生成する工程である。この工程は、ディジタ
ル化された受信信号強度より、物体を含む土壌１の断面
画像を、アンテナ１１及び１２の媒質表面上の位置
（ｘ，ｙ）と反射波５の物体２からの反射時間ｔ（実際
は所定の入射信号を発振してから反射信号が受信アンテ
ナにいたるまでの時間）を座標（ｘ，ｙ，ｔ）とする原
３次元ボクセルデータｓ（ｘ，ｙ，ｔ）として取り込む
工程であり、前記３次元ボクセルデータ生成手段２５に
より処理される。ここで、受信信号はその強度に応じて
複数階調で輝度表示され、信号強度の正値を白（輝度
大）、信号強度の負値を黒（輝度小）、信号強度０を中
間階調として取り込まれる。この階調は、具体的には、
８ビット（２５６）階調で表現され、階調１２８が反射
信号強度の振幅値０で、１２９以上の階調で振幅値が正
値で、１２７以下の階調で振幅値が負値となっている。
更に具体的には、ディジタル化された受信信号は、Ａ／
Ｄ変換処理されたときの量子化ビット幅で、媒質表面上
の位置（ｘ，ｙ）と反射波５の物体２からの反射時間ｔ
で決定される座標（ｘ，ｙ，ｔ）がアドレス信号として
エンコードされ、複数階調の原３次元ボクセルデータｓ
（ｘ，ｙ，ｔ）としてデータ処理部内のメモリ２１ａの
所定の領域に格納される。
１） この工程は、３次元探査装置３を移動しながら、データ
を収集し、これを処理して、以降の処理に使用される３
次元ボクセルデータｓ（ｘ，ｙ，ｔ）若しくはＳ（ｘ，
ｙ，ｔ）を生成する工程である。この工程は、ディジタ
ル化された受信信号強度より、物体を含む土壌１の断面
画像を、アンテナ１１及び１２の媒質表面上の位置
（ｘ，ｙ）と反射波５の物体２からの反射時間ｔ（実際
は所定の入射信号を発振してから反射信号が受信アンテ
ナにいたるまでの時間）を座標（ｘ，ｙ，ｔ）とする原
３次元ボクセルデータｓ（ｘ，ｙ，ｔ）として取り込む
工程であり、前記３次元ボクセルデータ生成手段２５に
より処理される。ここで、受信信号はその強度に応じて
複数階調で輝度表示され、信号強度の正値を白（輝度
大）、信号強度の負値を黒（輝度小）、信号強度０を中
間階調として取り込まれる。この階調は、具体的には、
８ビット（２５６）階調で表現され、階調１２８が反射
信号強度の振幅値０で、１２９以上の階調で振幅値が正
値で、１２７以下の階調で振幅値が負値となっている。
更に具体的には、ディジタル化された受信信号は、Ａ／
Ｄ変換処理されたときの量子化ビット幅で、媒質表面上
の位置（ｘ，ｙ）と反射波５の物体２からの反射時間ｔ
で決定される座標（ｘ，ｙ，ｔ）がアドレス信号として
エンコードされ、複数階調の原３次元ボクセルデータｓ
（ｘ，ｙ，ｔ）としてデータ処理部内のメモリ２１ａの
所定の領域に格納される。
【００４１】次に、このようにして得られた原３次元ボ
クセルデータｓ（ｘ，ｙ，ｔ）が前記マイグレーション
処理手段２６によってマイグレーション処理され、以降
の処理の対象となる３次元ボクセルデータＳ（ｘ，ｙ，
ｔ）が生成される。このようにして得られたマイグレー
ション処理後の３次元ボクセルデータＳ（ｘ，ｙ，ｔ）
が原データｓ（ｘ，ｙ，ｔ）と置換され、３次元ボクセ
ルデータ生成手段２５の出力として、以降の工程に引き
渡される。
クセルデータｓ（ｘ，ｙ，ｔ）が前記マイグレーション
処理手段２６によってマイグレーション処理され、以降
の処理の対象となる３次元ボクセルデータＳ（ｘ，ｙ，
ｔ）が生成される。このようにして得られたマイグレー
ション処理後の３次元ボクセルデータＳ（ｘ，ｙ，ｔ）
が原データｓ（ｘ，ｙ，ｔ）と置換され、３次元ボクセ
ルデータ生成手段２５の出力として、以降の工程に引き
渡される。
【００４２】以上の原３次元ボクセルデータｓ（ｘ，
ｙ，ｔ）及び３次元ボクセルデータＳ（ｘ，ｙ，ｔ）を
生成する処理は、前記３次元ボクセルデータ生成手段２
５と前記マイグレーション処理手段２６を実現するプロ
グラムを前記コンピュータの主記憶領域にインストール
し、前記コンピュータ上で実行することにより実行され
る。
ｙ，ｔ）及び３次元ボクセルデータＳ（ｘ，ｙ，ｔ）を
生成する処理は、前記３次元ボクセルデータ生成手段２
５と前記マイグレーション処理手段２６を実現するプロ
グラムを前記コンピュータの主記憶領域にインストール
し、前記コンピュータ上で実行することにより実行され
る。
【００４３】２．最大絶対値抽出工程（ＳＴ２） ３次元ボクセルデータ生成工程において生成された３次
元ボクセルデータＳ（ｘ，ｙ，ｔ）に対して、前記最大
絶対値抽出手段２７が各媒質表面上の位置（ｘ，ｙ）毎
に反射時間ｔ軸方向における振幅値を前記メモリ２１ａ
から逐次読み出し、その絶対値の最大値である最大絶対
値を抽出する。各媒質表面上の位置（ｘ，ｙ）毎に最大
絶対値を抽出したときに、その最大絶対値をとる反射時
間ｔＭＡＸ も併せて抽出する。これにより、処理後に
おいても、強い信号を反射する埋設物領域がどの反射時
間ｔ（埋設深さに相当）に存在するかを容易に把握する
ことが可能となる。
元ボクセルデータＳ（ｘ，ｙ，ｔ）に対して、前記最大
絶対値抽出手段２７が各媒質表面上の位置（ｘ，ｙ）毎
に反射時間ｔ軸方向における振幅値を前記メモリ２１ａ
から逐次読み出し、その絶対値の最大値である最大絶対
値を抽出する。各媒質表面上の位置（ｘ，ｙ）毎に最大
絶対値を抽出したときに、その最大絶対値をとる反射時
間ｔＭＡＸ も併せて抽出する。これにより、処理後に
おいても、強い信号を反射する埋設物領域がどの反射時
間ｔ（埋設深さに相当）に存在するかを容易に把握する
ことが可能となる。
【００４４】３．平面配置工程（ＳＴ３） 前記最大絶対値抽出工程において抽出された各媒質表面
上の位置（ｘ，ｙ）毎の最大絶対値を前記平面配置手段
２８によって所定平面上に配置する。
上の位置（ｘ，ｙ）毎の最大絶対値を前記平面配置手段
２８によって所定平面上に配置する。
【００４５】実際には、各媒質表面上の位置（ｘ，ｙ）
毎に、最大絶対値抽出工程と平面配置工程を連続して実
行することにより効率的に処理される。具体的には、各
媒質表面上の位置（ｘ，ｙ）毎に抽出した最大絶対値と
反射時間ｔＭＡＸ を前記メモリ２１ａの所定領域に格
納する。全ての媒質表面上の位置（ｘ，ｙ）に対して両
工程を実行すると、最終的に、各媒質表面上の位置
（ｘ，ｙ）毎の最大絶対値が所定平面上に配置されるこ
とになる。このとき、最大絶対値の格納領域をｔ＝０と
割り当てることにより、当該平面が３次元ボクセルデー
タＳ（ｘ，ｙ，ｔ）の最上面に配置される結果となる。
毎に、最大絶対値抽出工程と平面配置工程を連続して実
行することにより効率的に処理される。具体的には、各
媒質表面上の位置（ｘ，ｙ）毎に抽出した最大絶対値と
反射時間ｔＭＡＸ を前記メモリ２１ａの所定領域に格
納する。全ての媒質表面上の位置（ｘ，ｙ）に対して両
工程を実行すると、最終的に、各媒質表面上の位置
（ｘ，ｙ）毎の最大絶対値が所定平面上に配置されるこ
とになる。このとき、最大絶対値の格納領域をｔ＝０と
割り当てることにより、当該平面が３次元ボクセルデー
タＳ（ｘ，ｙ，ｔ）の最上面に配置される結果となる。
【００４６】従って、前記最大絶対値抽出手段２７と前
記平面配置手段２８は、上記した前記最大絶対値抽出工
程と前記平面配置工程を各媒質表面上の位置（ｘ，ｙ）
毎に連続して実行するためのプログラムを前記コンピュ
ータの主記憶領域にインストールし、前記コンピュータ
上で実行することにより実現される。
記平面配置手段２８は、上記した前記最大絶対値抽出工
程と前記平面配置工程を各媒質表面上の位置（ｘ，ｙ）
毎に連続して実行するためのプログラムを前記コンピュ
ータの主記憶領域にインストールし、前記コンピュータ
上で実行することにより実現される。
【００４７】４．平面表示工程（ＳＴ４） 前記平面表示手段２９が、平面配置工程において平面配
置された最大絶対値を出力処理して前記表示部２３に表
示する。具体的には、前記平面表示手段２９は、前記コ
ンピュータとして汎用のパーソナルコンピュータ等を使
用する場合、前記パーソナルコンピュータ等が基本的に
備えている所定の画像データを前記表示部２３に表示す
るための画像表示手段と、その画像表示手段に、前記平
面配置された最大絶対値を前記画像データとして出力す
る前処理手段とで構成される。前記画像表示手段と前記
前処理手段は、前記コンピュータをそれらの手段として
機能させるプログラムを前記コンピュータの主記憶領域
にインストールし、前記コンピュータ上で実行すること
により実現される。
置された最大絶対値を出力処理して前記表示部２３に表
示する。具体的には、前記平面表示手段２９は、前記コ
ンピュータとして汎用のパーソナルコンピュータ等を使
用する場合、前記パーソナルコンピュータ等が基本的に
備えている所定の画像データを前記表示部２３に表示す
るための画像表示手段と、その画像表示手段に、前記平
面配置された最大絶対値を前記画像データとして出力す
る前処理手段とで構成される。前記画像表示手段と前記
前処理手段は、前記コンピュータをそれらの手段として
機能させるプログラムを前記コンピュータの主記憶領域
にインストールし、前記コンピュータ上で実行すること
により実現される。
【００４８】図５に示す埋設管を探査した３次元ボクセ
ルデータを処理した結果について、本工程において画像
表示した例を図６に示す。図６に示した表示画像は、従
来技術に比べて視認性が改善され、埋設状況の平面配置
が容易に把握できることが分かる。
ルデータを処理した結果について、本工程において画像
表示した例を図６に示す。図６に示した表示画像は、従
来技術に比べて視認性が改善され、埋設状況の平面配置
が容易に把握できることが分かる。
【００４９】〔別実施形態〕以下に別実施形態を説明す
る。
る。
【００５０】〈１〉前記最大絶対値抽出工程において、
各媒質表面上の位置（ｘ，ｙ）毎に反射時間ｔ軸方向に
おける最大絶対値を抽出する際に、反射時間ｔの所定範
囲内にある振幅値を除外するのも好ましい実施の形態で
ある。このようにして抽出した最大絶対値を、上記実施
例と同様に平面表示工程で画像表示した例を図７に示
す。尚、図７に示す実施例では、３次元ボクセルデータ
Ｓ（ｘ，ｙ，ｔ）の反射時間ｔの範囲を８ｎｓ〜７８ｎ
ｓ（埋設深さに換算して０．４ｍ〜３．９ｍ）の振幅値
のみを対象として前記最大絶対値の抽出を行なってい
る。図７より、深い埋設管の弱い反射信号がより明解に
把握できるようになったことが分かる。
各媒質表面上の位置（ｘ，ｙ）毎に反射時間ｔ軸方向に
おける最大絶対値を抽出する際に、反射時間ｔの所定範
囲内にある振幅値を除外するのも好ましい実施の形態で
ある。このようにして抽出した最大絶対値を、上記実施
例と同様に平面表示工程で画像表示した例を図７に示
す。尚、図７に示す実施例では、３次元ボクセルデータ
Ｓ（ｘ，ｙ，ｔ）の反射時間ｔの範囲を８ｎｓ〜７８ｎ
ｓ（埋設深さに換算して０．４ｍ〜３．９ｍ）の振幅値
のみを対象として前記最大絶対値の抽出を行なってい
る。図７より、深い埋設管の弱い反射信号がより明解に
把握できるようになったことが分かる。
【００５１】また、各媒質表面上の位置（ｘ，ｙ）毎に
反射時間ｔ軸方向における最大絶対値を抽出する際に、
振幅値の極性が正または負の何れか一方のボクセルのみ
を対象とすることもＳＮ比の改善という観点から好まし
い。
反射時間ｔ軸方向における最大絶対値を抽出する際に、
振幅値の極性が正または負の何れか一方のボクセルのみ
を対象とすることもＳＮ比の改善という観点から好まし
い。
【００５２】〈２〉前記３次元ボクセルデータ生成工程
において、前記原３次元ボクセルデータｓ（ｘ，ｙ，
ｔ）に対して、複数の地中伝搬速度（媒質中の伝搬速
度）を用いて複数の合成開口処理若しくはマイグレーシ
ョン処理を施して前記複数の伝搬速度毎の３次元ボクセ
ルデータＳ（ｘ，ｙ，ｔ）を生成し、夫々に対して前記
最大絶対値抽出工程と前記平面配置工程と前記平面表示
工程を各別に実行するようにしても構わない。上記処理
は、具体的には、予め所定の記憶領域に記憶してある前
記複数の地中伝搬速度を逐次読み出し、或いは、前記複
数の地中伝搬速度が相互に一定の関係をもって離間して
いる場合は、その関係を表す関係式に基づいて一つの地
中伝搬速度を算出して、その地中伝搬速度を入力変数と
して、前記合成開口処理若しくは前記マイグレーション
処理用のプログラムに与え、当該プログラムを実行させ
前記３次元ボクセルデータＳ（ｘ，ｙ，ｔ）を生成し、
前記最大絶対値抽出工程と前記平面配置工程と前記平面
表示工程を上記した要領で夫々のプログラムを実行し
て、その結果を画面表示する。そして、この一連の処理
は、前記複数の地中伝搬速度に対して順番に実行するプ
ログラムを、前記マイグレーション処理手段２６（また
は、合成開口処理手段）、前記最大絶対値抽出手段２
７、前記平面配置手段２８、前記平面表示手段２９を実
現する各プログラムとともに前記コンピュータの主記憶
領域にインストールし、前記コンピュータ上で実行する
ことにより実現される。
において、前記原３次元ボクセルデータｓ（ｘ，ｙ，
ｔ）に対して、複数の地中伝搬速度（媒質中の伝搬速
度）を用いて複数の合成開口処理若しくはマイグレーシ
ョン処理を施して前記複数の伝搬速度毎の３次元ボクセ
ルデータＳ（ｘ，ｙ，ｔ）を生成し、夫々に対して前記
最大絶対値抽出工程と前記平面配置工程と前記平面表示
工程を各別に実行するようにしても構わない。上記処理
は、具体的には、予め所定の記憶領域に記憶してある前
記複数の地中伝搬速度を逐次読み出し、或いは、前記複
数の地中伝搬速度が相互に一定の関係をもって離間して
いる場合は、その関係を表す関係式に基づいて一つの地
中伝搬速度を算出して、その地中伝搬速度を入力変数と
して、前記合成開口処理若しくは前記マイグレーション
処理用のプログラムに与え、当該プログラムを実行させ
前記３次元ボクセルデータＳ（ｘ，ｙ，ｔ）を生成し、
前記最大絶対値抽出工程と前記平面配置工程と前記平面
表示工程を上記した要領で夫々のプログラムを実行し
て、その結果を画面表示する。そして、この一連の処理
は、前記複数の地中伝搬速度に対して順番に実行するプ
ログラムを、前記マイグレーション処理手段２６（また
は、合成開口処理手段）、前記最大絶対値抽出手段２
７、前記平面配置手段２８、前記平面表示手段２９を実
現する各プログラムとともに前記コンピュータの主記憶
領域にインストールし、前記コンピュータ上で実行する
ことにより実現される。
【００５３】この結果、前記各伝搬速度毎の前記平面表
示工程における前記３次元ボクセルデータＳ（ｘ，ｙ，
ｔ）の前記最大絶対値の表示結果から、適切な伝搬速度
による処理結果を選択することができるのである。地中
伝搬速度を０．５０Ｃ0 、０．３３Ｃ0 、０．２３Ｃ0
の３通りに設定した実施結果を図８、図９、図１０の順
に示す。ここで、Ｃ0 は真空中での電磁波の伝搬速度で
あり、各伝搬速度に対応する地中の比誘電率は４、９、
１９である。図８、図９、図１０より、地中伝搬速度が
０．３３Ｃ0 （比誘電率９）の時のフォーカス（図９）
が最も良いので、この処理結果を適切な地中伝搬速度に
よる処理結果と判断することができる。
示工程における前記３次元ボクセルデータＳ（ｘ，ｙ，
ｔ）の前記最大絶対値の表示結果から、適切な伝搬速度
による処理結果を選択することができるのである。地中
伝搬速度を０．５０Ｃ0 、０．３３Ｃ0 、０．２３Ｃ0
の３通りに設定した実施結果を図８、図９、図１０の順
に示す。ここで、Ｃ0 は真空中での電磁波の伝搬速度で
あり、各伝搬速度に対応する地中の比誘電率は４、９、
１９である。図８、図９、図１０より、地中伝搬速度が
０．３３Ｃ0 （比誘電率９）の時のフォーカス（図９）
が最も良いので、この処理結果を適切な地中伝搬速度に
よる処理結果と判断することができる。
【００５４】〈３〉上記実施形態では、前記３次元ボク
セルデータ生成工程、前記最大絶対値抽出工程、前記平
面配置工程、及び、前記平面表示工程の各工程を処理す
る前記３次元ボクセルデータ生成手段２５、前記マイグ
レーション処理手段２６（または、合成開口処理手
段）、前記最大絶対値抽出手段２７、前記平面配置手段
２８、前記平面表示手段２９（前記画像表示手段と前記
前処理手段）は、前記コンピュータとその上で実行可能
な各手段に対応したプログラム、つまり汎用ハードウェ
アと専用ソフトウェアとの組み合わせで構成されていた
が、前記各手段或いは一部の手段を夫々専用のハードウ
ェアで構成しても、または、前記手段の中の一部の機能
を専用のハードウェアで構成しても構わない。また、前
記プログラムを実行可能なコンピュータは、必ずしも一
台のコンピュータである必要はない。つまり、一部のプ
ログラムをその処理に適切なアーキテクチャの特別なコ
ンピュータを使用して当該プログラムを実行させるよう
にしても構わない。
セルデータ生成工程、前記最大絶対値抽出工程、前記平
面配置工程、及び、前記平面表示工程の各工程を処理す
る前記３次元ボクセルデータ生成手段２５、前記マイグ
レーション処理手段２６（または、合成開口処理手
段）、前記最大絶対値抽出手段２７、前記平面配置手段
２８、前記平面表示手段２９（前記画像表示手段と前記
前処理手段）は、前記コンピュータとその上で実行可能
な各手段に対応したプログラム、つまり汎用ハードウェ
アと専用ソフトウェアとの組み合わせで構成されていた
が、前記各手段或いは一部の手段を夫々専用のハードウ
ェアで構成しても、または、前記手段の中の一部の機能
を専用のハードウェアで構成しても構わない。また、前
記プログラムを実行可能なコンピュータは、必ずしも一
台のコンピュータである必要はない。つまり、一部のプ
ログラムをその処理に適切なアーキテクチャの特別なコ
ンピュータを使用して当該プログラムを実行させるよう
にしても構わない。
【００５５】〈４〉次に、本発明に係るコンピュータ読
み取り可能な記録媒体について説明する。上記したよう
に、本発明に係る３次元ボクセルデータ表示方法は、前
記コンピュータを、前記３次元ボクセルデータ生成手段
２５、前記マイグレーション処理手段２６（または、合
成開口処理手段）、前記最大絶対値抽出手段２７、前記
平面配置手段２８、前記平面表示手段２９（前記画像表
示手段と前記前処理手段）として機能させるためのプロ
グラムによって実行するには、前記プログラムを先ず前
記コンピュータの主記憶領域にインストールする必要が
あることから、これらのインストール前のプログラムを
前記コンピュータが読み取り可能な記録媒体、例えば、
ＣＤ−ＲＯＭ（ｃｏｍｐａｃｔ ｄｉｓｃ ｒｅａｄ
ｏｎｌｙ ｍｅｍｏｒｙ）等の可搬型の記録媒体に保存
しておき、３次元ボクセルデータ表示装置として機能さ
せる各コンピュータに前記記録媒体を介して前記プログ
ラムを前記主記憶領域にインストールする。更には、当
該インストールをコンピュータ・ネットワークを介して
行う場合は、前記インストール前の各プログラムを所定
のサーバ上の記憶媒体（ハードディスク装置等）に前記
３次元ボクセルデータ表示装置として機能させるコンピ
ュータから読み出し可能に記憶しておいても構わない。
み取り可能な記録媒体について説明する。上記したよう
に、本発明に係る３次元ボクセルデータ表示方法は、前
記コンピュータを、前記３次元ボクセルデータ生成手段
２５、前記マイグレーション処理手段２６（または、合
成開口処理手段）、前記最大絶対値抽出手段２７、前記
平面配置手段２８、前記平面表示手段２９（前記画像表
示手段と前記前処理手段）として機能させるためのプロ
グラムによって実行するには、前記プログラムを先ず前
記コンピュータの主記憶領域にインストールする必要が
あることから、これらのインストール前のプログラムを
前記コンピュータが読み取り可能な記録媒体、例えば、
ＣＤ−ＲＯＭ（ｃｏｍｐａｃｔ ｄｉｓｃ ｒｅａｄ
ｏｎｌｙ ｍｅｍｏｒｙ）等の可搬型の記録媒体に保存
しておき、３次元ボクセルデータ表示装置として機能さ
せる各コンピュータに前記記録媒体を介して前記プログ
ラムを前記主記憶領域にインストールする。更には、当
該インストールをコンピュータ・ネットワークを介して
行う場合は、前記インストール前の各プログラムを所定
のサーバ上の記憶媒体（ハードディスク装置等）に前記
３次元ボクセルデータ表示装置として機能させるコンピ
ュータから読み出し可能に記憶しておいても構わない。
【００５６】ここで、前記記録媒体には、前記コンピュ
ータを前記最大絶対値抽出手段２７と前記平面配置手段
２８、及び、前記平面表示手段２９の内の前記前処理手
段として機能させるためのプログラムが少なくとも記憶
されている必要がある。前記３次元ボクセルデータ生成
手段２５、前記マイグレーション処理手段２６（また
は、合成開口処理手段）、前記平面表示手段２９の内の
前記画像表示手段は、既存のプログラムや前記コンピュ
ータの基本機能として具備されているものが使用可能で
ある。但し、汎用コンピュータにおいて必ずしもこれら
のプログラムがインストール済とは限らないので、同じ
記憶媒体にこれらのプログラムの一部または全部を記憶
していても構わない。
ータを前記最大絶対値抽出手段２７と前記平面配置手段
２８、及び、前記平面表示手段２９の内の前記前処理手
段として機能させるためのプログラムが少なくとも記憶
されている必要がある。前記３次元ボクセルデータ生成
手段２５、前記マイグレーション処理手段２６（また
は、合成開口処理手段）、前記平面表示手段２９の内の
前記画像表示手段は、既存のプログラムや前記コンピュ
ータの基本機能として具備されているものが使用可能で
ある。但し、汎用コンピュータにおいて必ずしもこれら
のプログラムがインストール済とは限らないので、同じ
記憶媒体にこれらのプログラムの一部または全部を記憶
していても構わない。
【００５７】更に、上記別実施形態〈２〉で説明した、
複数の地中伝搬速度に対する前記一連の処理を行うプロ
グラムを前記記録媒体に記録しても構わない。
複数の地中伝搬速度に対する前記一連の処理を行うプロ
グラムを前記記録媒体に記録しても構わない。
【００５８】尚、記憶媒体は必ずしも物理的に単体の記
憶装置である必要はない。例えば、複数のＣＤ−ＲＯＭ
や、コンピュータ・ネットワーク上に分散された複数の
サーバ上に前記各プログラムが分散して配置されていて
も構わない。結局、前記３次元ボクセルデータ表示装置
として機能させるコンピュータの主記憶領域に必要なプ
ログラムがインストールされ得るからである。
憶装置である必要はない。例えば、複数のＣＤ−ＲＯＭ
や、コンピュータ・ネットワーク上に分散された複数の
サーバ上に前記各プログラムが分散して配置されていて
も構わない。結局、前記３次元ボクセルデータ表示装置
として機能させるコンピュータの主記憶領域に必要なプ
ログラムがインストールされ得るからである。
【図１】３次元ボクセルデータ表示装置を備えた探査装
置のブロック構成図
置のブロック構成図
【図２】送信信号と受信信号の波形説明図
【図３】データ解析装置の機能ブロック図
【図４】本発明の３次元ボクセルデータ表示方法を示す
フローチャート
フローチャート
【図５】本発明の３次元ボクセルデータ表示方法に使用
した探査データを採取した場所の埋設状況を示す説明図
した探査データを採取した場所の埋設状況を示す説明図
【図６】本発明の３次元ボクセルデータ表示方法により
得られた処理結果を示す表示画像例
得られた処理結果を示す表示画像例
【図７】本発明の別実施形態による３次元ボクセルデー
タ表示方法により得られた処理結果を示す表示画像例
タ表示方法により得られた処理結果を示す表示画像例
【図８】本発明の他の別実施形態による３次元ボクセル
データ表示方法により得られた処理結果を示す表示画像
例（地中の比誘電率は４）
データ表示方法により得られた処理結果を示す表示画像
例（地中の比誘電率は４）
【図９】本発明の他の別実施形態による３次元ボクセル
データ表示方法により得られた処理結果を示す表示画像
例（地中の比誘電率は９）
データ表示方法により得られた処理結果を示す表示画像
例（地中の比誘電率は９）
【図１０】本発明の他の別実施形態による３次元ボクセ
ルデータ表示方法により得られた処理結果を示す表示画
像例（地中の比誘電率は１９）
ルデータ表示方法により得られた処理結果を示す表示画
像例（地中の比誘電率は１９）
【図１１】従来技術の一水平断面における処理結果を示
す表示画像例
す表示画像例
【図１２】従来技術の表面表示型で表示した場合の処理
結果を示す表示画像例
結果を示す表示画像例
【図１３】従来技術の表面表示型で表示した場合の処理
結果を示す表示画像例
結果を示す表示画像例
【図１４】従来技術の表面表示型で表示した場合の処理
結果を示す表示画像例
結果を示す表示画像例
【図１５】従来技術の積分型で表示した場合の処理結果
を示す表示画像例
を示す表示画像例
１ 土壌（媒質） ２ 物体 ３ ３次元探査装置 ４ 入射波 ５ 反射波 １０ 送受信機 ２０ データ解析装置 ２１ データ処理部 ２１ａ メモリ ２２ 入力部 ２３ 表示部 ２４ 外部補助記憶部 ２５ ３次元ボクセルデータ生成手段 ２６ マイグレーション処理手段 ２７ 最大絶対値抽出手段 ２８ 平面配置手段 ２９ 平面表示手段 ３０ 埋設管
─────────────────────────────────────────────────────
【手続補正書】
【提出日】平成１１年１１月２２日（１９９９．１１．
２２）
２２）
【手続補正１】
【補正対象書類名】図面
【補正対象項目名】図６
【補正方法】変更
【補正内容】
【図６】
【手続補正２】
【補正対象書類名】図面
【補正対象項目名】図７
【補正方法】変更
【補正内容】
【図７】
【手続補正３】
【補正対象書類名】図面
【補正対象項目名】図８
【補正方法】変更
【補正内容】
【図８】
【手続補正４】
【補正対象書類名】図面
【補正対象項目名】図９
【補正方法】変更
【補正内容】
【図９】
【手続補正５】
【補正対象書類名】図面
【補正対象項目名】図１０
【補正方法】変更
【補正内容】
【図１０】
【手続補正６】
【補正対象書類名】図面
【補正対象項目名】図１１
【補正方法】変更
【補正内容】
【図１１】
【手続補正７】
【補正対象書類名】図面
【補正対象項目名】図１２
【補正方法】変更
【補正内容】
【図１２】
【手続補正８】
【補正対象書類名】図面
【補正対象項目名】図１３
【補正方法】変更
【補正内容】
【図１３】
【手続補正９】
【補正対象書類名】図面
【補正対象項目名】図１４
【補正方法】変更
【補正内容】
【図１４】
【手続補正１０】
【補正対象書類名】図面
【補正対象項目名】図１５
【補正方法】変更
【補正内容】
【図１５】
Claims (11)
- 【請求項１】 媒質表面から媒質中へ放射した波動信号
の反射信号強度に基づいて前記媒質表面上の位置（ｘ，
ｙ）と反射時間ｔを座標（ｘ，ｙ，ｔ）として生成され
た３次元ボクセルデータの表示方法であって、 前記３次元ボクセルデータに対して、前記媒質表面上の
位置（ｘ，ｙ）毎に反射時間ｔ軸方向における振幅値の
最大絶対値を抽出する最大絶対値抽出工程と、前記最大
絶対値抽出工程で抽出された前記媒質表面上の位置
（ｘ，ｙ）毎の前記最大絶対値を所定平面上に配置する
平面配置工程と、前記所定平面上に配置された前記最大
絶対値を表示する平面表示工程を実行する３次元ボクセ
ルデータ表示方法。 - 【請求項２】 前記最大絶対値抽出工程で抽出された前
記最大絶対値の中から最大値と最小値を求め、その最大
値と最小値が夫々表示階調の上限及び下限となるように
前記最大絶対値を正規化することを特徴とする請求項１
記載の３次元ボクセルデータ表示方法。 - 【請求項３】 前記最大絶対値抽出工程において、前記
媒質表面上の位置（ｘ，ｙ）毎に前記最大絶対値を抽出
する際に、その最大絶対値をとる反射時間ｔＭＡＸ も
併せて抽出することを特徴とする請求項１または２記載
の３次元ボクセルデータ表示方法。 - 【請求項４】 前記最大絶対値抽出工程において、前記
媒質表面上の位置（ｘ，ｙ）毎に前記最大絶対値を抽出
する際に、反射時間ｔの所定範囲内にある振幅値を除外
することを特徴とする請求項１、２または３記載の３次
元ボクセルデータ表示方法。 - 【請求項５】 前記最大絶対値抽出工程において、前記
媒質表面上の位置（ｘ，ｙ）毎に前記最大絶対値を抽出
する際に、振幅値の極性が正または負の何れか一方のボ
クセルのみを対象とすることを特徴とする請求項１、
２、３または４記載の３次元ボクセルデータ表示方法。 - 【請求項６】 前記最大絶対値抽出工程前に、前記反射
信号強度からなる３次元ボクセルデータに対して合成開
口処理若しくはマイグレーション処理を施して、表示対
象となる前記３次元ボクセルデータを生成することを特
徴とする請求項１、２、３、４または５記載の３次元ボ
クセルデータ表示方法。 - 【請求項７】 前記合成開口処理若しくはマイグレーシ
ョン処理を複数の媒質中の伝搬速度で実行して、前記伝
搬速度毎に前記３次元ボクセルデータを生成し、前記伝
搬速度毎の前記各３次元ボクセルデータに対して、前記
最大絶対値抽出工程と前記平面配置工程と前記平面表示
工程を各別に実行し、前記各伝搬速度毎の前記平面表示
工程における前記３次元ボクセルデータの表示結果か
ら、適切な伝搬速度による処理結果を選択することを特
徴とする請求項６記載の３次元ボクセルデータ表示方
法。 - 【請求項８】 媒質表面から媒質中へ放射した波動信号
の反射信号強度に基づいて前記媒質表面上の位置（ｘ，
ｙ）と反射時間ｔを座標（ｘ，ｙ，ｔ）として生成され
た３次元ボクセルデータを表示する３次元ボクセルデー
タ表示装置であって、 前記３次元ボクセルデータに対して、前記媒質表面上の
位置（ｘ，ｙ）毎に反射時間ｔ軸方向における振幅値の
最大絶対値を抽出する最大絶対値抽出手段と、前記最大
絶対値抽出手段で抽出された前記媒質表面上の位置
（ｘ，ｙ）毎の前記最大絶対値を所定平面上に配置する
平面配置手段と、前記所定平面上に配置された前記最大
絶対値を表示する平面表示手段とを備えてなる３次元ボ
クセルデータ表示装置。 - 【請求項９】 請求項１、２、３、４または５記載の３
次元ボクセルデータ表示方法における前記最大絶対値抽
出工程と前記平面配置工程を１または２以上のコンピュ
ータに実行させるためのプログラムと、前記３次元ボク
セルデータ表示方法における前記平面表示工程或いは前
記平面表示工程の内の前記最大絶対値を被処理データと
して前記コンピュータに設けられた画像表示手段に出力
する手順を前記コンピュータに実行させるためのプログ
ラムとを少なくとも記録したコンピュータ読み取り可能
な記録媒体。 - 【請求項１０】 請求項６記載の３次元ボクセルデータ
表示方法における前記合成開口処理と前記マイグレーシ
ョン処理の少なくとも何れか一方を前記コンピュータに
実行させるためのプログラムを記録した請求項９記載の
コンピュータ読み取り可能な記録媒体。 - 【請求項１１】 請求項７記載の３次元ボクセルデータ
表示方法における前記合成開口処理若しくは前記マイグ
レーション処理を複数の媒質中の伝搬速度で実行し、前
記伝搬速度毎に前記３次元ボクセルデータを生成する処
理を前記コンピュータに実行させるためのプログラム
と、前記伝搬速度毎の前記各３次元ボクセルデータに対
して、前記最大絶対値抽出工程と前記平面配置工程と前
記平面表示工程或いは前記平面表示工程の内の前記最大
絶対値を被処理データとして前記コンピュータに設けら
れた画像表示手段に出力する手順を各別に前記コンピュ
ータに実行させるためのプログラムとを記録した請求項
１０記載のコンピュータ読み取り可能な記録媒体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP33096699A JP3409001B2 (ja) | 1998-11-24 | 1999-11-22 | ３次元ボクセルデータ表示方法及び装置 |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP33245998 | 1998-11-24 | ||
| JP10-332459 | 1998-11-24 | ||
| JP33096699A JP3409001B2 (ja) | 1998-11-24 | 1999-11-22 | ３次元ボクセルデータ表示方法及び装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2000221266A true JP2000221266A (ja) | 2000-08-11 |
| JP3409001B2 JP3409001B2 (ja) | 2003-05-19 |
Family
ID=26573688
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP33096699A Expired - Lifetime JP3409001B2 (ja) | 1998-11-24 | 1999-11-22 | ３次元ボクセルデータ表示方法及び装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3409001B2 (ja) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002031537A1 (fr) * | 2000-10-12 | 2002-04-18 | Cos Co., Ltd. | Procede pour etalonner la permittivite relative d'un detecteur electromagnetique et detecteur electromagnetique |
| JP2004301610A (ja) * | 2003-03-31 | 2004-10-28 | Secom Co Ltd | 地中空洞検出装置 |
| JP2006071416A (ja) * | 2004-09-01 | 2006-03-16 | Anritsu Corp | 印刷はんだ検査装置 |
| WO2006046483A1 (ja) * | 2004-10-26 | 2006-05-04 | Matsushita Electric Industrial Co., Ltd. | 車載用レーダ装置 |
| JP2008128974A (ja) * | 2006-11-24 | 2008-06-05 | Matsushita Electric Works Ltd | 物体探知装置 |
| JP2019100882A (ja) * | 2017-12-04 | 2019-06-24 | 三菱電機株式会社 | 地中空洞検出装置および地中空洞検出方法 |
| JP2019168464A (ja) * | 2017-08-03 | 2019-10-03 | 本多電子株式会社 | 海図画像表示装置 |
-
1999
- 1999-11-22 JP JP33096699A patent/JP3409001B2/ja not_active Expired - Lifetime
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002031537A1 (fr) * | 2000-10-12 | 2002-04-18 | Cos Co., Ltd. | Procede pour etalonner la permittivite relative d'un detecteur electromagnetique et detecteur electromagnetique |
| JP2004301610A (ja) * | 2003-03-31 | 2004-10-28 | Secom Co Ltd | 地中空洞検出装置 |
| JP2006071416A (ja) * | 2004-09-01 | 2006-03-16 | Anritsu Corp | 印刷はんだ検査装置 |
| WO2006046483A1 (ja) * | 2004-10-26 | 2006-05-04 | Matsushita Electric Industrial Co., Ltd. | 車載用レーダ装置 |
| JP2008128974A (ja) * | 2006-11-24 | 2008-06-05 | Matsushita Electric Works Ltd | 物体探知装置 |
| JP2019168464A (ja) * | 2017-08-03 | 2019-10-03 | 本多電子株式会社 | 海図画像表示装置 |
| JP2019100882A (ja) * | 2017-12-04 | 2019-06-24 | 三菱電機株式会社 | 地中空洞検出装置および地中空洞検出方法 |
| JP7030496B2 (ja) | 2017-12-04 | 2022-03-07 | 三菱電機株式会社 | 地中空洞検出装置および地中空洞検出方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3409001B2 (ja) | 2003-05-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2000013037A1 (fr) | Procede de recherche tridimensionnel, procede d'affichage de donnees de voxels tridimensionnelles, et dispositif de realisation de ces procedes | |
| CN106446919A (zh) | 一种探地雷达双曲线目标快速检测方法 | |
| JP3409001B2 (ja) | ３次元ボクセルデータ表示方法及び装置 | |
| CN103439707A (zh) | 一种探地雷达加窗距离偏移成像方法 | |
| CN108761449B (zh) | 一种钢筋回波干扰下的灾害目标成像方法 | |
| Zhao et al. | A Comprehensive Horizon‐Picking Method on Subbottom Profiles by Combining Envelope, Phase Attributes, and Texture Analysis | |
| JP3409002B2 (ja) | ３次元探査方法及び装置 | |
| JPH116879A (ja) | ３次元探査方法及び装置 | |
| JP2961336B2 (ja) | 不可視物体探査方法 | |
| JP4059609B2 (ja) | ３次元ボクセルデータ表示方法 | |
| CN115201816A (zh) | 一种基于改进后向投影的三维探地雷达量化成像方法 | |
| JP3263752B2 (ja) | 埋設物探知に用いる受信信号の処理方法 | |
| JP3005353B2 (ja) | 隠蔽場所の断面検出装置 | |
| JPH09211121A (ja) | 探査方法及び装置 | |
| Wang et al. | Characteristic parameters extraction method of hidden Karst Cave from Borehole radar signal | |
| US10859695B2 (en) | Acoustic system and method for characterizing granular media | |
| Hu et al. | Traveltime tomography of crosshole ground-penetrating radar based on an arctangent functional with compactness constraints | |
| JP2000075025A (ja) | ３次元探査方法及び装置 | |
| Alshamy et al. | Evaluation of GPR Detection for buried objects material with different depths and scanning angles | |
| JPH11271440A (ja) | ３次元探査方法及び装置 | |
| JP3223897B2 (ja) | 地中レーダ信号処理装置 | |
| Srinivas et al. | Buried pipes detection through ground penetrating radar images | |
| Cooper et al. | The use of textural analysis to locate features in geophysical data | |
| JP2893010B1 (ja) | 探査方法及び装置 | |
| JPH10141907A (ja) | 探査解析方法及び探査解析装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| TRDD | Decision of grant or rejection written | ||
| R150 | Certificate of patent or registration of utility model |
Ref document number: 3409001 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120314 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150314 Year of fee payment: 12 |
|
| EXPY | Cancellation because of completion of term |