JP3987451B2 - Ground surface fluctuation measurement method using synthetic aperture radar - Google Patents
Ground surface fluctuation measurement method using synthetic aperture radar Download PDFInfo
- Publication number
- JP3987451B2 JP3987451B2 JP2003099501A JP2003099501A JP3987451B2 JP 3987451 B2 JP3987451 B2 JP 3987451B2 JP 2003099501 A JP2003099501 A JP 2003099501A JP 2003099501 A JP2003099501 A JP 2003099501A JP 3987451 B2 JP3987451 B2 JP 3987451B2
- Authority
- JP
- Japan
- Prior art keywords
- radio wave
- ground surface
- synthetic aperture
- aperture radar
- fluctuation
- 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
- Radar Systems Or Details Thereof (AREA)
Description
【0001】
【発明の属する技術分野】
本発明は、人工衛星などの飛翔体に搭載されている合成開口レーダ(SAR)から発射されたパルス状電波の地表面反射波を合成開口干渉処理することにより計測対象地点の地表面に沿った変動量を計測する方法に関するものである。この技術は、特に限定されるものではないが、例えば地滑り予兆の調査などの微小な地表面変動量の遠隔計測に有用である。
【0002】
【従来の技術】
【特許文献1】
特開平7−72244号公報
【0003】
地表面を観測する技術の一つに、航空機や人工衛星等の飛翔体に搭載した合成開口レーダを用いる方法がある。これは、航行する飛翔体から一定間隔で進行方向に対して垂直斜め下方に短いパルス状の電波を発射し、地表面からの反射波を受信して、合成開口干渉処理と呼ばれる高度な信号処理を行うことにより、地表面の性状を観測する技術である。なお、この合成開口干渉処理とは、対象物からの反射波を2台の空間的に離れたアンテナで受信し、その2つの反射波信号の間の位相差を用いて、三角測量と類似の原理により対象物の方位や距離を求める方法である。
【0004】
合成開口レーダによる計測では、1基の飛翔体に複数の合成開口レーダを搭載する方式、あるいは複数の飛翔体にそれぞれ合成開口レーダを搭載する方式がある。その他、より一層簡便なシステムとして、1台の合成開口レーダを搭載した1基のみの飛翔体を用いる方式もある。この方式は、同一の合成開口レーダのみを使用することから、本質的に特性のばらつきがなく、計測精度の向上に有効である。
【0005】
この種の合成開口レーダは、小型のアンテナでも仮想的に大きなアンテナを用いた場合と同様の非常に高い分解能で観測が可能なこと、光学カメラと異なり昼夜間を問わず且つ全天候で観測が可能であること、などの優れた特徴がある。そこで、例えば特許文献1には、飛翔体に搭載した2台のアンテナを用いて同時に送信波を出力し同時に受信波を入力するようにして、地震予知、火山噴火予知、土石流予測等のように、数cm単位の精度で微妙な地形変動の抽出が要求される分野に適用可能な干渉型合成開口レーダ装置及び地形変動観測方式が開示されている。
【0006】
ところで、1台の合成開口レーダを搭載した1基のみの人工衛星を用いても、合成開口レーダから発射されたパルス状電波の地表面反射波を合成開口干渉処理することによって地表面上における特定の地点(計測対象地点)の変動量を計測することはできる。例えば、その計測対象地点に1個の電波反射体を設置し、昇交軌道あるいは降交軌道のいずれかで2回反射波を取得し、それらの反射波データで合成開口干渉処理を行う方法である。近年の技術開発の進展によって、計測精度はますます向上しており、cmオーダー若しくはそれ以下の計測精度が得られている。
【0007】
【発明が解決しようとする課題】
前記のように合成開口レーダによる地表面観測は、昼夜間を問わず全天候で計測が可能であり、高精度で行えるという優れた特徴がある。しかし、従来方法によって得られる地表面の変動量は、軌道上の人工衛星と電波反射体を結ぶ斜め方向での変動量でしかなく、地表面変動の有無やその程度は検出できるが、地表面に沿った実際の変動量を求めることはできない。そのため、地滑り予兆の観測など、地表面に沿った微小な変位を正確に計測することはできなかった。
【0008】
地表面上における特定の地点の変動を直接計測する方法として、GPS衛星を用いることも考えられる。GPS衛星を利用することで地表面上の特定地点の位置が直接求まることから、任意の時間経過時点でGPS衛星による位置計測を行い過去位置と比較すれば、原理的には計測対象となる地表面の変動量を直接求めることができる。しかし現時点では、GPS衛星による位置計測は精度が低く、地滑りの予兆観測に必要なcmオーダー若しくはそれ以下の要求精度には対応できない。また複数のGPS衛星を用いる計測となるため、衛星の組み合わせによる特性のばらつきなどの影響もあって、本質的に精度が上がらず、微小な地表面変動の計測には不向きである。
【0009】
本発明の目的は、飛翔体に搭載されている1台の合成開口レーダを用いて、地滑りの予兆など計測対象地点の地表面に沿った実際の変動量を高精度で求めることができる地表面変動量計測方法を提供することである。本発明の他の目的は、合成開口レーダを用いる地表面変動量計測方法の実施に適し、積雪時でも埋没することが無く、そのため季節を問わず継続して地表面変動を計測でき、また雪や雨水などの影響による反射特性の劣化が生じ難い電波反射構造物を提供することである。
【0010】
【課題を解決するための手段】
本発明は、飛翔体に搭載されている1台の合成開口レーダから発射されるパルス状の電波の地表面反射波を合成開口干渉処理することにより地表面の変動量を求める方法において、計測対象となる地表面に剛構造の基体を設け、該基体に一対の電波反射体を、一方は昇交軌道上の合成開口レーダが発する電波を反射する向きに、他方は降交軌道上の合成開口レーダが発する電波を反射する向きに設置して、昇交軌道にて2回取得した反射波同士を合成開口干渉処理することにより昇交軌道から計測対象地点までの距離を求めると共に、降交軌道にて2回取得した反射波同士を合成開口干渉処理することにより降交軌道から計測対象地点までの距離を求め、両方の距離データを同時に満たし且つ地盤構造上特定される地表面の変動方向の条件を加えることで計測対象地点の地表面に沿った過去位置からの変動量を求めることを特徴とする合成開口レーダを用いる地表面変動量計測方法である。
【0011】
ここで「飛翔体」とは人工衛星やスペースシャトルなど宇宙空間を極軌道で航行するプラットフォームを指す。本発明では、1基のみの飛翔体を用いて計測でき、その方が好ましい。しかし、昇交軌道と降交軌道とで別の飛翔体を用いて計測することもできる。
【0012】
地盤構造上特定される地表面の変動方向の条件としては、典型的には、地滑り変動が斜面の最大傾斜線に沿って発生するという条件があり、それによって地滑り変動量を求めることができる。その他、活断層の存在によって地表面の変動方向が限られている場合には、その方向を地表面の変動方向の条件として加えることになる。
【0013】
また本発明は、そのような合成開口レーダを用いた地表面変動量計測方法で用いる電波反射構造物である。この電波反射構造物は、地表面に立設される剛構造のタワー型基体と、該基体に互いに逆向きに設置される一対の電波反射体とを具備し、前記基体は電波吸収体で覆われ、電波反射体は積雪に埋没しない高さに調整された中段と上段にそれぞれ設置された構造をなしている。
【0014】
ここで電波反射体は、例えば3枚の直角二等辺三角形金属板を組み合わせ直角部分を互いに接合することで三角錐状にしたコーナリフレクタからなり、該コーナリフレクタの開口部が撥水性を呈する電波透過体で塞がれている構造とする。あるいは電波反射体は、3枚の直角二等辺三角形金属板を組み合わせ直角部分を互いに接合することで三角錐状にしたコーナリフレクタからなり、該コーナリフレクタの底部に排水穴を形成した構造でもよい。その他、平板状の電波反射鏡などの使用も可能である。
【0015】
【実施例】
図1は本発明に係る合成開口レーダを用いた地表面変動量計測方法の一実施例を示す説明図である。計測対象となる地表面10に一対の電波反射体12a,12bを備えた電波反射構造物14を動かないように設置する。そして、人工衛星16に搭載されている合成開口レーダ18からパルス状の電波を発射し、計測対象となる地表面に設置されている電波反射体12a,12bからの反射波を受信して合成開口干渉処理する。この結果を用いることによって地表面の変動量を求める。本実施例で使用する人工衛星16は1基のみであり、それには合成開口レーダ18が1台搭載されていればよい。
【0016】
本発明方法では、人工衛星16が昇交軌道航行時に受信した電波反射体12aからの反射波データと、同じ人工衛星16が降交軌道航行時に受信した電波反射体12bからの反射波データの両方を使用する。そのため、電波反射構造物14には、図示のように、一対の電波反射体12a,12bが互いに逆向きに取り付けられる。図1において、Aは昇交軌道時の状態を、Bは降交軌道時の状態を、それぞれ示している。従って、一方の電波反射体12aは昇交軌道上の合成開口レーダ18が発射する電波を反射し、他方の電波反射体12bは降交軌道上の合成開口レーダが発射する電波を反射するように、それぞれの向きが調整されている。
【0017】
本実施例では、電波反射構造物14は、計測対象となる地表面10に強固に立設した剛構造のタワー型基体20と、該タワー型基体20上に互いに逆向きに設置される一対の電波反射体12a,12bとからなる。タワー型基体20は、風雪などが作用しても撓まないような剛構造であり、例えば3本の支柱22を三角柱を形成するように配設して、その中段と上段とに棚24a,24bを固定し、一方の電波反射体12aを中段の棚24aに、他方の電波反射体12bを上段の棚24bに、それぞれ取り付ける。タワー型基体20は、特に積雪地帯に設置する場合には、両方の電波反射体12a,12bが積雪に埋没しないように支柱22の長さ、棚24a,24bの高さを調整する。
【0018】
本発明方法において、合成開口レーダ18を搭載している人工衛星16は極軌道を航行する。合成開口レーダ18は、一定間隔で進行方向に対して垂直斜め下方(例えば左側下方)に地表面に向けて短いパルス状の電波(例えばCバンド:周波数5.3GHz)を発射する。従って、例えば昇交軌道(衛星が南側から北側へ向かう軌道)では東側から電波が入射して一方の電波反射体12aで反射し、降交軌道(衛星が北側から南側へ向かう軌道)では西側から電波が入射して他方の電波反射体12bで反射する。人工衛星16に搭載されている合成開口レーダ18からの電波は発射方向が決まっており、人工衛星16は、その軌道も決まっている(軌道のずれは補正できる)。そこで、軌道データと画像データの両方を取得し、それらに基づいて必要な解析を行う。
【0019】
本発明では、昇交軌道で2回反射波データを取得し、降交軌道でも2回反射波データを取得する。そして、昇交軌道にて得られた反射波同士を合成開口干渉処理することにより、昇交軌道上の人工衛星16から計測対象地点までの距離がcmオーダーもしくはそれ以下の精度で求められる。同様に、降交軌道にて得られた反射波同士を合成開口干渉処理することにより、降交軌道上の人工衛星16から計測対象地点までの距離がcmオーダーもしくはそれ以下の精度で求められる。従って、過去の計測結果との間で差が生じていれば、その差分はその間の地表面変動に起因するものとなる。
【0020】
図2は、本発明方法により実際の地表面変動量の求め方を示す概念図である。変動前における計測対象地点(過去位置)をP1で表す。変動前における昇交軌道上の人工衛星位置Saと計測対象地点P1の間の直線距離をLa0とすると、ある時間が経過して変動した後の距離はLa0+daで表せる。このとき変動後の計測対象地点は、昇交軌道上の人工衛星位置Saを中心とする半径La0+daの球面30上のどこかに位置することになる。同様に、変動前における降交軌道上の人工衛星位置Sdと計測対象地点P1の間の直線距離をLd0とし、それが変動後にLd0+ddに変化したとする。すると、変動後の計測対象地点は、降交軌道上の人工衛星位置Sdを中心とする半径Ld0+ddの球面32上のどこかに位置することになる。ここで、円34は一方の球面30と他方の球面32との交線を表しており、昇交軌道での計測値と降交軌道での計測値を同時に満たしている点の集合である。従って、変動後の計測対象地点は、この円34上のどこかに存在することになる。しかし、それだけでは実際に変動した計測対象地点は特定できない。
【0021】
そこで本発明では、計測対象地点における地盤構造を勘案して、地盤構造上特定される地表面の変動方向を条件として加える。つまり、地表面が動く方向(あるいは動く可能性がある方向)を地盤工学的な見地から仮定する。ここでは、斜面における地滑り予兆の観測を目的としているために、計測対象地点は斜面上にあるものとする。すると、地表面の変動は、斜面の最大傾斜線に沿った方向(矢印36で示す方向)に生じると仮定できる。この条件を加えると、変動前の計測対象地点P1から最大傾斜線方向36と円34との交点P2が求まり、この点P2が実際の変動後の計測対象地点を示すことになる。このようにして、計測対象地点の地表面に沿った実際の変動量を計測することが可能となる。
【0022】
上記の説明では斜面の地滑り予兆の観測を目的としているために、計測対象地点が斜面上にあるものとし、地表面の変動は、斜面の最大傾斜線に沿って生じると仮定した。本発明方法を、その他、例えば活断層地帯での地表面変動量計測に適用する場合には、地表面変動が活断層に沿って生じると仮定することにより、変動後の計測対象地点を特定でき、その地表面に沿った実際の変動量を計測することができる。
【0023】
本発明方法で用いる電波吸収構造体としては、図1に示すように、地表面に立設される剛構造のタワー型基体20と、該タワー型基体20に互いに逆向きに設置される一対の電波反射体12a,12bとからなる構造が望ましい。前記タワー型基体20には、その支柱22や棚24a,24bなどに電波吸収体を塗布したり、あるいは電波吸収体を貼り付けるのが好ましく、それによって電波反射体12a,12b以外の部分からの不要反射波を抑制できる。また、両方の電波反射体12a,12bを中段と上段に取り付けることで同一地点に設置することができ、且つ高さ調整を行うことで積雪に埋没することなく四季を通して観測を継続することが可能となる。
【0024】
電波反射体としては、電波反射鏡でもよいが、図3に示すように、3枚の直角二等辺三角形金属板40を三角錐状に組み合わせ、それら各金属板40の直角部分を溶接などにより接合したコーナリフレクタ42が好ましい。設置の際には、その開口部が人工衛星位置を向くように方位などが調整される。これによって、合成開口レーダ電波を、その入射方向に一致する向きで人工衛星に向けて戻すことができる。
【0025】
このようなコーナリフレクタ42は、開口部からの雨や雪などの浸入を防止するために、図3のAに示すように開口部を電波透過性の板材44で覆うようにしてもよい。その場合、雨や雪などが電波透過性の板材44の表面に付着しないように、該板材44を撥水性の材料で作製したり、該板材44の表面に撥水性の膜を設けるのが好ましい。開口部を覆わない場合には、図3のBに示すように、コーナリフレクタ42の奥底部に小径の排水穴46を形成し、開口部から浸入した雨水などを速やかに排出可能とするのが好ましい。
【0026】
本発明では、電波反射体として、コーナリフレクタや電波反射鏡のような受動電波反射体のみならず、トランスポンダのような能動電波反射体も利用できる。その場合には、無指向性のアンテナを用いることで、1台のトランスポンダで昇交軌道及び降交軌道の両方の合成開口レーダ電波を反射するような構成も可能である。但し、トランスポンダは、受信した電波を増幅して送信する機能を有するものであるから、環境温度などによって電気的特性が変化する可能性もあり、位相差で計測する合成開口干渉処理において十分な精度を確保するためには、実際にはかなりの困難を伴う。そのような事情を勘案すると、実施例で説明したように、受動電波反射体、特にコーナリフレクタを利用する方法は安価であり、反射性能も良好で安定しており、設置も容易であるため好ましい。
【0027】
【発明の効果】
本発明は上記のように、合成開口レーダを用い、一対の電波反射体を用いて昇交軌道と降交軌道とで地表面反射波を合成開口干渉処理することで昇交軌道と降交軌道から測定対象地点までの距離をそれぞれ求め、更に地表面の変動方向の条件を加えて処理する方法であるから、地表面に沿った実際の変動量を正確に求めることができる。また本発明では、特に1台の人工衛星搭載合成開口レーダで計測を行うことができるために、計測値のばらつきが少なく、精度が向上し、しかも安価に実施できる利点がある。
【0028】
また本発明に係る電波反射構造体は、電波反射体が積雪時でも埋没することが無いように構成でき、そのため季節に無関係に継続的に地表面変動を計測することができる。また、コーナリフレクタの開口部を電波透過性の板体で覆ったり、コーナリフレクタの底部に排水穴を設けることで、雪や雨水などの影響を受け難くでき、反射特性の変化や劣化が生じ難くなり、良好な計測を長期間にわたり継続して実施することが可能となる。
【図面の簡単な説明】
【図1】本発明に係る地表面変動量計測方法の一実施例を示す説明図。
【図2】本発明方法による実際の地表面変動量の求め方を示す概念図。
【図3】電波反射体の構造例を示す説明図。
【符号の説明】
10 計測対象となる地表面
12a,12b 電波反射体
14 電波反射構造物
16 人工衛星
18 合成開口レーダ
20 タワー型基体[0001]
BACKGROUND OF THE INVENTION
The present invention is based on the ground surface of a measurement target point by performing synthetic aperture interference processing on the ground surface reflected wave of a pulsed radio wave emitted from a synthetic aperture radar (SAR) mounted on a flying object such as an artificial satellite. The present invention relates to a method for measuring a fluctuation amount. Although this technique is not specifically limited, For example, it is useful for the remote measurement of the minute ground surface fluctuation amount, such as investigation of a landslide sign.
[0002]
[Prior art]
[Patent Document 1]
Japanese Patent Laid-Open No. 7-72244
One technique for observing the ground surface is to use a synthetic aperture radar mounted on a flying object such as an aircraft or an artificial satellite. This is an advanced signal processing called synthetic aperture interference processing that emits a short pulse radio wave obliquely downward and perpendicular to the direction of travel from a navigating projectile and receives a reflected wave from the ground surface. This is a technique for observing the properties of the earth's surface. Note that this synthetic aperture interference processing is similar to triangulation using a phase difference between two reflected wave signals received from two spatially separated antennas. This is a method for obtaining the direction and distance of an object based on the principle.
[0004]
In the measurement by the synthetic aperture radar, there are a method in which a plurality of synthetic aperture radars are mounted on a single flying object, or a method in which a synthetic aperture radar is mounted on each of a plurality of flying objects. As another simpler system, there is a method using only one flying object equipped with one synthetic aperture radar. Since this method uses only the same synthetic aperture radar, there is essentially no variation in characteristics and is effective in improving measurement accuracy.
[0005]
This type of Synthetic Aperture Radar enables observation with very high resolution, even when using a small antenna, virtually the same as when a large antenna is used. Unlike optical cameras, it can be observed day and night and in all weather. It has excellent characteristics such as being. Therefore, for example, Patent Document 1 discloses that, for example, earthquake prediction, volcanic eruption prediction, debris flow prediction, etc., using two antennas mounted on a flying object, simultaneously outputting a transmission wave and simultaneously inputting a reception wave. In addition, an interferometric synthetic aperture radar apparatus and a terrain change observation method that are applicable to a field that requires extraction of subtle terrain fluctuations with an accuracy of several centimeters are disclosed.
[0006]
By the way, even if only one artificial satellite equipped with a single synthetic aperture radar is used, it is possible to identify on the ground surface by performing synthetic aperture interference processing on the ground surface reflected wave of the pulsed radio wave emitted from the synthetic aperture radar. It is possible to measure the fluctuation amount of the point (measurement target point). For example, a method is used in which one radio wave reflector is installed at the measurement target point, a reflected wave is acquired twice in either an ascending or descending orbit, and synthetic aperture interference processing is performed using the reflected wave data. is there. With the progress of technological development in recent years, measurement accuracy has been further improved, and measurement accuracy of cm order or less has been obtained.
[0007]
[Problems to be solved by the invention]
As described above, the surface observation by the synthetic aperture radar has an excellent feature that it can be measured in all weather regardless of day and night, and can be performed with high accuracy. However, the fluctuation amount of the ground surface obtained by the conventional method is only the fluctuation amount in the oblique direction connecting the satellite on the orbit and the radio wave reflector, and the presence or degree of the ground surface fluctuation can be detected. The actual amount of variation along the line cannot be determined. Therefore, it was not possible to accurately measure minute displacements along the ground surface, such as observing signs of landslides.
[0008]
As a method for directly measuring a change at a specific point on the ground surface, a GPS satellite may be used. Since the position of a specific point on the ground surface can be obtained directly by using a GPS satellite, if the position is measured by a GPS satellite at an arbitrary time point and compared with the past position, in principle, the position to be measured The amount of surface fluctuation can be obtained directly. However, at present, the position measurement by the GPS satellite is low in accuracy, and cannot meet the required accuracy of the order of cm or less necessary for predicting landslides. In addition, since measurement is performed using a plurality of GPS satellites, there is an influence such as variation in characteristics due to the combination of the satellites, and thus the accuracy is not essentially improved, and it is not suitable for measuring minute ground surface fluctuations.
[0009]
An object of the present invention is to use a single synthetic aperture radar mounted on a flying object, and to obtain an actual fluctuation amount along a ground surface of a measurement target point such as a sign of landslide with high accuracy. It is to provide a variation measurement method. Another object of the present invention is suitable for implementation of a method for measuring the amount of ground surface fluctuation using a synthetic aperture radar, and is not buried even during snowfall, so that ground surface fluctuation can be continuously measured regardless of the season. It is an object to provide a radio wave reflecting structure in which the reflection characteristics are hardly deteriorated due to the influence of rainwater or the like.
[0010]
[Means for Solving the Problems]
The present invention relates to a method for determining a ground surface fluctuation amount by performing synthetic aperture interference processing on a ground-surface reflected wave of a pulsed radio wave emitted from a single synthetic aperture radar mounted on a flying object. A ground structure is provided on the ground surface, a pair of radio wave reflectors are provided on the base, one is in the direction to reflect the radio wave emitted by the synthetic aperture radar on the ascending orbit, and the other is the synthetic aperture on the descending orbit Installed in the direction to reflect the radio waves emitted by the radar, and obtained the distance from the ascending orbit to the measurement target point by performing synthetic aperture interference processing of the reflected waves acquired twice in the ascending orbit, and descending orbit The distance from the descending trajectory to the measurement target point is obtained by performing synthetic aperture interference processing between the reflected waves acquired twice in step 1, and both the distance data are satisfied at the same time and the ground surface fluctuation direction specified on the ground structure is determined. The conditions It is ground surface variation measurement method using a synthetic aperture radar, characterized in that determining the amount of variation from the past position along the ground surface of the measurement object point by obtaining.
[0011]
Here, the “flying object” refers to a platform that navigates outer space in polar orbit, such as an artificial satellite or a space shuttle. In the present invention, measurement can be performed using only one flying object, which is preferable. However, the ascending and descending orbits can be measured using different flying objects.
[0012]
The condition of the fluctuation direction of the ground surface specified by the ground structure typically includes a condition that the landslide fluctuation occurs along the maximum slope line of the slope, and the amount of landslide fluctuation can be obtained thereby. In addition, when the fluctuation direction of the ground surface is limited due to the presence of active faults, this direction is added as a condition for the fluctuation direction of the ground surface.
[0013]
Moreover, this invention is a radio wave reflection structure used with the ground surface fluctuation amount measuring method using such a synthetic aperture radar. The radio wave reflecting structure includes a tower-type base body having a rigid structure standing on the ground surface and a pair of radio wave reflectors installed in opposite directions on the base body, and the base body is covered with a radio wave absorber. However, the radio wave reflector has a structure that is installed in the middle and upper stages, adjusted to a height not buried in the snow.
[0014]
Here, the radio wave reflector is composed of, for example, a corner reflector formed by combining three right-angled isosceles triangular metal plates and joining the right-angle portions to each other, and the opening of the corner reflector has a water repellency. The structure is closed by the body. Alternatively, the radio wave reflector may be composed of a corner reflector formed by combining three right-angled isosceles triangular metal plates and joining the right-angle portions to each other to form a triangular pyramid, and a drain hole is formed at the bottom of the corner reflector. In addition, it is possible to use a flat radio wave reflector.
[0015]
【Example】
FIG. 1 is an explanatory diagram showing an embodiment of a ground surface variation measuring method using a synthetic aperture radar according to the present invention. A radio
[0016]
In the method of the present invention, both the reflected wave data from the
[0017]
In the present embodiment, the radio
[0018]
In the method of the present invention, the
[0019]
In the present invention, the reflected wave data is acquired twice in the ascending orbit, and the reflected wave data is also acquired in the descending orbit. Then, by performing synthetic aperture interference processing on the reflected waves obtained in the ascending orbit, the distance from the
[0020]
FIG. 2 is a conceptual diagram showing how to determine the actual ground surface variation by the method of the present invention. A measurement target point (past position) before the change is represented by P1. Assuming that the linear distance between the satellite position Sa on the ascending orbit before the change and the measurement target point P1 is La0, the distance after the change after a certain time can be expressed by La0 + da. At this time, the changed measurement target point is located somewhere on the
[0021]
Therefore, in the present invention, the ground structure at the measurement target point is taken into consideration, and the fluctuation direction of the ground surface specified on the ground structure is added as a condition. That is, the direction in which the ground surface moves (or the direction in which it can move) is assumed from the viewpoint of geotechnical engineering. Here, since the purpose is to observe a sign of landslide on the slope, it is assumed that the measurement target point is on the slope. Then, it can be assumed that the fluctuation of the ground surface occurs in the direction along the maximum slope line of the slope (the direction indicated by the arrow 36). When this condition is added, the intersection P2 between the maximum
[0022]
In the above explanation, since the purpose is to observe the landslide sign of the slope, it is assumed that the measurement target point is on the slope, and that the ground surface changes along the maximum slope line of the slope. When the method of the present invention is applied to, for example, ground surface fluctuation amount measurement in an active fault zone, for example, it is possible to specify the measurement target point after the fluctuation by assuming that the ground surface fluctuation occurs along the active fault. The actual amount of fluctuation along the ground surface can be measured.
[0023]
As shown in FIG. 1, the radio wave absorption structure used in the method of the present invention includes a tower-
[0024]
The radio wave reflector may be a radio wave reflector, but as shown in FIG. 3, three right-angled isosceles
[0025]
Such a
[0026]
In the present invention, not only passive radio wave reflectors such as corner reflectors and radio wave reflectors but also active radio wave reflectors such as transponders can be used as radio wave reflectors. In that case, by using an omnidirectional antenna, it is also possible to use a single transponder to reflect synthetic aperture radar radio waves in both ascending and descending orbits. However, since the transponder has a function to amplify and transmit the received radio wave, the electrical characteristics may change depending on the environmental temperature, etc., and sufficient accuracy in the synthetic aperture interference processing that measures by the phase difference It is actually quite difficult to secure. In view of such circumstances, as described in the embodiment, a method using a passive radio wave reflector, particularly a corner reflector is preferable because it is inexpensive, has good and stable reflection performance, and is easy to install. .
[0027]
【The invention's effect】
As described above, the present invention uses a synthetic aperture radar and a pair of radio wave reflectors to perform ascending and descending orbits by performing synthetic aperture interference processing on the ground surface reflected waves in the ascending and descending orbits. In this method, the distance from the measurement target point to the measurement target point is obtained, and the condition of the fluctuation direction of the ground surface is further added and processed. Therefore, the actual fluctuation amount along the ground surface can be accurately obtained. In the present invention, in particular, since the measurement can be performed by one artificial satellite mounted synthetic aperture radar, there is an advantage that measurement values are less varied, the accuracy is improved, and the measurement can be performed at low cost.
[0028]
In addition, the radio wave reflecting structure according to the present invention can be configured so that the radio wave reflector is not buried even during snowfall, and therefore it is possible to continuously measure ground surface variations regardless of the season. In addition, by covering the opening of the corner reflector with a radio wave permeable plate or providing a drain hole at the bottom of the corner reflector, it is less susceptible to the effects of snow, rainwater, etc., and changes or deterioration of the reflection characteristics are less likely to occur. Therefore, it becomes possible to carry out good measurement continuously over a long period of time.
[Brief description of the drawings]
FIG. 1 is an explanatory diagram showing an embodiment of a ground surface variation measuring method according to the present invention.
FIG. 2 is a conceptual diagram showing how to determine the actual ground surface fluctuation amount by the method of the present invention.
FIG. 3 is an explanatory diagram showing a structure example of a radio wave reflector.
[Explanation of symbols]
DESCRIPTION OF
Claims (5)
計測対象となる地表面に剛構造の基体を設け、該基体に一対の電波反射体を、一方は昇交軌道上の合成開口レーダが発する電波を反射する向きに、他方は降交軌道上の合成開口レーダが発する電波を反射する向きに設置して、昇交軌道にて2回取得した反射波同士を合成開口干渉処理することにより昇交軌道から計測対象地点までの距離を求めると共に、降交軌道にて2回取得した反射波同士を合成開口干渉処理することにより降交軌道から計測対象地点までの距離を求め、両方の距離データを同時に満たし且つ地盤構造上特定される地表面の変動方向の条件を加えることで計測対象地点の地表面に沿った過去位置からの変動量を求めることを特徴とする合成開口レーダを用いる地表面変動量計測方法。In a method for determining the amount of fluctuation of the ground surface by performing synthetic aperture interference processing on the ground surface reflected wave of the pulsed radio wave emitted from one synthetic aperture radar mounted on the flying object,
A rigid base is provided on the ground surface to be measured, and a pair of radio wave reflectors are provided on the base, one in a direction to reflect radio waves emitted by the synthetic aperture radar on the ascending orbit, and the other on the descending orbit. Installed in a direction to reflect the radio wave emitted by the synthetic aperture radar, and by performing synthetic aperture interference processing of the reflected waves acquired twice in the ascending orbit, the distance from the ascending orbit to the measurement target point is obtained, Obtain the distance from the descending trajectory to the measurement target point by performing synthetic aperture interference processing between the reflected waves acquired twice in the alternating trajectory, and satisfy both distance data at the same time and change the ground surface specified in the ground structure A ground surface variation measuring method using a synthetic aperture radar, characterized in that a variation amount from a past position along a ground surface of a measurement target point is obtained by adding a direction condition.
地表面に立設される剛構造のタワー型基体と、該基体に互いに逆向きに設置される一対の電波反射体とを具備し、前記基体は電波吸収体で覆われ、電波反射体は積雪に埋没しない高さに調整された中段と上段にそれぞれ設置されている地表面変動量計測用の電波反射構造物。A radio wave reflecting structure used in a ground surface variation measuring method using the synthetic aperture radar according to claim 1 or 2,
A tower-type base body having a rigid structure standing on the ground surface and a pair of radio wave reflectors installed opposite to each other on the base body, the base body being covered with a radio wave absorber, and the radio wave reflector being snow covered Radio wave reflection structure for measuring ground surface fluctuations installed in the middle and upper tiers, adjusted to a height that is not buried in the ground.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003099501A JP3987451B2 (en) | 2003-04-02 | 2003-04-02 | Ground surface fluctuation measurement method using synthetic aperture radar |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003099501A JP3987451B2 (en) | 2003-04-02 | 2003-04-02 | Ground surface fluctuation measurement method using synthetic aperture radar |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004309178A JP2004309178A (en) | 2004-11-04 |
| JP3987451B2 true JP3987451B2 (en) | 2007-10-10 |
Family
ID=33463940
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2003099501A Expired - Lifetime JP3987451B2 (en) | 2003-04-02 | 2003-04-02 | Ground surface fluctuation measurement method using synthetic aperture radar |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3987451B2 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10006995B2 (en) | 2014-08-04 | 2018-06-26 | University Of Seoul Industry Cooperation Foundation | Method and apparatus for stacking multi-temporal MAI interferograms |
| KR102751028B1 (en) | 2024-09-30 | 2025-01-09 | 한국건설기술연구원 | Satellite SAR reflector for landfill monitoring |
Families Citing this family (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5207713B2 (en) * | 2007-11-29 | 2013-06-12 | 上田日本無線株式会社 | Reflector for millimeter wave radar |
| JP5401890B2 (en) * | 2008-09-25 | 2014-01-29 | 日本電気株式会社 | Interferometric synthetic aperture radar system, processing method, air vehicle control device, air vehicle control method, and program |
| KR101111689B1 (en) * | 2009-11-03 | 2012-02-14 | 연세대학교 산학협력단 | 3D surface displacement extraction method and device |
| IT1401469B1 (en) * | 2010-06-28 | 2013-07-26 | Tele Rilevamento Europa T R E S R L | SYSTEM FOR MEASURING THE MOVEMENT OF SOIL WITH IMAGES USING AT LEAST A REFLECTOR OF ELECTROMAGNETIC WAVES |
| CN101964135A (en) * | 2010-10-13 | 2011-02-02 | 南京泰通科技有限公司 | Device for monitoring landslide |
| FR2970818B1 (en) * | 2011-01-25 | 2014-01-24 | Commissariat Energie Atomique | PASSIVE ELECTROMAGNETIC WAVE REFLECTOR FOR LOCAL STRUCTURAL DEFORMATION MEASUREMENT AT EARTH SURFACE |
| CN102645650A (en) * | 2012-03-06 | 2012-08-22 | 北京北科安地科技发展有限公司 | Landslide dynamic identifying and monitoring technology based on synthetic aperture radar differential interferometry (D-InSAR) |
| KR101315243B1 (en) | 2013-03-07 | 2013-10-08 | 중앙항업(주) | Transformation methods of doppler frequency of radar images |
| KR101498131B1 (en) * | 2014-01-07 | 2015-03-04 | 이성 주식회사 | Underground Target Detection System and Method Based On Adaptive Signal Estimation |
| JP6314101B2 (en) * | 2015-02-26 | 2018-04-18 | 株式会社パスコ | Ground displacement measurement method corresponding to vegetation area using ground-mounted synthetic aperture radar |
| JP6555522B2 (en) * | 2015-09-01 | 2019-08-07 | 清水建設株式会社 | Information processing apparatus, information processing method, and program |
| JP6588322B2 (en) * | 2015-12-09 | 2019-10-09 | 東日本旅客鉄道株式会社 | Slope monitoring system and slope monitoring method |
| JP2017166949A (en) * | 2016-03-16 | 2017-09-21 | 国立研究開発法人防災科学技術研究所 | Terrain variation determination method and terrain variation determination system |
| JP6179911B1 (en) * | 2016-06-01 | 2017-08-16 | 国立研究開発法人宇宙航空研究開発機構 | Deformation degree determination method and modification degree determination system |
| KR102174582B1 (en) * | 2018-03-22 | 2020-11-05 | 주식회사 무림지앤아이 | Corner Reflector and Interferometric Synthetic Aperture Radar Using It |
| KR102142674B1 (en) * | 2019-08-01 | 2020-08-07 | 서울시립대학교 산학협력단 | Method and Apparatus for Synthetic Aperture Radar Phase Unwrapping based on SAR Offset Tracking Displacement Model |
| CN111781560B (en) * | 2020-08-07 | 2024-08-27 | 机械工业勘察设计研究院有限公司 | Adjustable corner reflector and mounting method thereof |
| CN113740844B (en) * | 2021-09-09 | 2024-04-02 | 雷添杰 | Dam three-dimensional deformation monitoring-oriented two-foundation radar combined observation method |
| CN119437025B (en) * | 2024-10-25 | 2025-10-28 | 中国长江电力股份有限公司 | Beidou+InSAR corner reflector deformation monitoring data fusion method and precision evaluation method |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6311885A (en) * | 1986-07-01 | 1988-01-19 | Natl Space Dev Agency Japan<Nasda> | Method for continuously observing specific area in earth observation |
| JPH03128876U (en) * | 1990-04-09 | 1991-12-25 | ||
| JP2590689B2 (en) * | 1993-06-14 | 1997-03-12 | 日本電気株式会社 | Interferometric synthetic aperture radar system and terrain change observation method |
| JPH0815426A (en) * | 1994-06-30 | 1996-01-19 | Nec Corp | Terrain alteration extracting device using interference type synthetic aperture radar |
| JP3301292B2 (en) * | 1995-12-12 | 2002-07-15 | 三菱電機株式会社 | Interferometric high-resolution radar device and terrain height measuring method using high-resolution radar device |
| JP3301297B2 (en) * | 1996-01-08 | 2002-07-15 | 三菱電機株式会社 | Surface deformation observation method using radar device, synthetic aperture radar device and transponder used in this method |
| JP2001091649A (en) * | 1999-09-20 | 2001-04-06 | Mitsubishi Electric Corp | Ground control device for precise geometric correction of synthetic aperture radar images |
-
2003
- 2003-04-02 JP JP2003099501A patent/JP3987451B2/en not_active Expired - Lifetime
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10006995B2 (en) | 2014-08-04 | 2018-06-26 | University Of Seoul Industry Cooperation Foundation | Method and apparatus for stacking multi-temporal MAI interferograms |
| KR102751028B1 (en) | 2024-09-30 | 2025-01-09 | 한국건설기술연구원 | Satellite SAR reflector for landfill monitoring |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2004309178A (en) | 2004-11-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3987451B2 (en) | Ground surface fluctuation measurement method using synthetic aperture radar | |
| Milan et al. | Interferometric evidence for the observation of ground backscatter originating behind the CUTLASS coherent HF radars | |
| Tomiyasu et al. | Synthetic aperture radar imaging from an inclined geosynchronous orbit | |
| CN105445730B (en) | A kind of Sea Current inverting Spaceborne SAR System and its method based on angle diversity | |
| CN102169173B (en) | Method for analyzing ambiguity of inclined geo-synchronization orbit synthetic aperture radar | |
| RU2446411C2 (en) | Method of monitoring displacements of earth's surface and deformation of structures on territory of mineral deposits | |
| US12607711B2 (en) | Synthetic aperture radar corner reflector | |
| CN114879197A (en) | Method for calibrating satellite DDM (distributed data management) in real time | |
| RU2411538C2 (en) | Method of determining error in measuring aircraft velocity with inertial navigation system and onboard navigation system for realising said method | |
| Aher | Synthetic aperture radar in Indian remote sensing | |
| KR102151362B1 (en) | Image decoding apparatus based on airborn using polar coordinates transformation and method of decoding image using the same | |
| Kolpuke et al. | Airborne multichannel UWB FMCW radar for snow depth measurements | |
| US8519886B2 (en) | System and method for using GPS interferometry to determine atmospheric conditions | |
| Reyhanigalangashi et al. | Ultra-wideband bistatic radar measurements of snow | |
| CA2480752A1 (en) | Method of reducing the doppler centroid in a coherent pulse radar system | |
| CN116165662A (en) | SAR antenna distance-to-double-pass directional diagram measuring method based on passive satellite | |
| Rouveure et al. | High resolution mapping of the environment with a ground-based radar imager | |
| Łabowski et al. | Inertial navigation system for radar terrain imaging | |
| Chicarella et al. | Improvement of GPR tracking by using inertial and GPS combined data | |
| Baskakov et al. | Problem of detecting space debris objects using multi-position radar system | |
| CN112455503A (en) | Train positioning method and device based on radar | |
| Rouveure et al. | Mapping of the environment with a high resolution ground-based radar imager | |
| JP2003302466A (en) | Antenna pattern correction method for synthetic aperture radar | |
| Bernhardt et al. | Bistatic observations of the ocean surface with HF radar, satellite and airborne receivers | |
| CN114280548A (en) | Satellite-borne downward-looking ice-detection synthetic aperture radar transmission path calculation method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20051116 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070629 |
|
| 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: 20070710 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070712 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 3987451 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100720 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100720 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130720 Year of fee payment: 6 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| 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 |
|
| 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 |