Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP3987451B2 - Ground surface fluctuation measurement method using synthetic aperture radar - Google Patents
[go: Go Back, main page]

JP3987451B2 - Ground surface fluctuation measurement method using synthetic aperture radar - Google Patents

Ground surface fluctuation measurement method using synthetic aperture radar Download PDF

Info

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
Application number
JP2003099501A
Other languages
Japanese (ja)
Other versions
JP2004309178A (en
Inventor
敏実 水野
誠 山根
成樹 葛岡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Oyo Corp
ImageOne Co Ltd
Original Assignee
Oyo Corp
ImageOne Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Oyo Corp, ImageOne Co Ltd filed Critical Oyo Corp
Priority to JP2003099501A priority Critical patent/JP3987451B2/en
Publication of JP2004309178A publication Critical patent/JP2004309178A/en
Application granted granted Critical
Publication of JP3987451B2 publication Critical patent/JP3987451B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

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 wave reflecting structure 14 having a pair of radio wave reflectors 12a and 12b is installed on the ground surface 10 to be measured so as not to move. Then, a pulsed radio wave is emitted from the synthetic aperture radar 18 mounted on the artificial satellite 16, and the reflected waves from the radio wave reflectors 12a and 12b installed on the ground surface to be measured are received and the synthetic aperture is received. Interference processing. By using this result, the fluctuation amount of the ground surface is obtained. There is only one artificial satellite 16 used in this embodiment, and only one synthetic aperture radar 18 needs to be mounted thereon.
[0016]
In the method of the present invention, both the reflected wave data from the radio wave reflector 12a received by the artificial satellite 16 during the ascending orbit travel and the reflected wave data from the radio wave reflector 12b received by the same artificial satellite 16 during the descending orbit navigation. Is used. Therefore, as shown in the figure, a pair of radio wave reflectors 12a and 12b are attached to the radio wave reflecting structure 14 in opposite directions. In FIG. 1, A shows the state during the ascending orbit, and B shows the state during the descending orbit. Accordingly, one radio wave reflector 12a reflects the radio wave emitted by the synthetic aperture radar 18 on the ascending trajectory, and the other radio wave reflector 12b reflects the radio wave emitted by the synthetic aperture radar on the descending orbit. , Each orientation is adjusted.
[0017]
In the present embodiment, the radio wave reflecting structure 14 includes a tower-type base body 20 having a rigid structure that is firmly erected on the ground surface 10 to be measured, and a pair of the base-type base bodies 20 that are installed in opposite directions on the tower-type base body 20. It consists of radio wave reflectors 12a and 12b. The tower-type base body 20 has a rigid structure that does not bend even when wind or snow acts. For example, three pillars 22 are arranged so as to form a triangular prism, and shelves 24a, 24b is fixed, and one radio wave reflector 12a is attached to the middle shelf 24a, and the other radio wave reflector 12b is attached to the upper shelf 24b. When the tower base 20 is installed particularly in a snowy area, the length of the column 22 and the height of the shelves 24a and 24b are adjusted so that both radio wave reflectors 12a and 12b are not buried in the snow.
[0018]
In the method of the present invention, the artificial satellite 16 carrying the synthetic aperture radar 18 travels in a polar orbit. The synthetic aperture radar 18 emits a short pulse-shaped radio wave (for example, C band: frequency 5.3 GHz) toward the ground surface obliquely downward (for example, the lower left side) perpendicular to the traveling direction at regular intervals. Thus, for example, in an ascending orbit (orbit from the south side to the north side), radio waves are incident from the east side and reflected by one radio wave reflector 12a, and in a descending orbit (orbit from the north side to the south side) A radio wave enters and is reflected by the other radio wave reflector 12b. The emission direction of the radio wave from the synthetic aperture radar 18 mounted on the artificial satellite 16 is determined, and the orbit of the artificial satellite 16 is determined (orbital deviation can be corrected). Therefore, both trajectory data and image data are acquired, and necessary analysis is performed based on them.
[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 artificial satellite 16 on the ascending orbit to the measurement target point can be obtained with an accuracy of cm order or less. Similarly, by performing synthetic aperture interference processing on the reflected waves obtained in the descending orbit, the distance from the artificial satellite 16 on the descending orbit to the measurement target point can be obtained with an accuracy of cm order or less. Therefore, if there is a difference between the past measurement results, the difference is caused by the ground surface fluctuation during that time.
[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 spherical surface 30 having a radius La0 + da centered on the artificial satellite position Sa on the ascending orbit. Similarly, it is assumed that the linear distance between the satellite position Sd on the descending orbit before the change and the measurement target point P1 is Ld0 and changes to Ld0 + dd after the change. Then, the measurement target point after the change is located somewhere on the spherical surface 32 having the radius Ld0 + dd centered on the satellite position Sd on the descending orbit. Here, the circle 34 represents an intersection line between the one spherical surface 30 and the other spherical surface 32, and is a set of points that simultaneously satisfy the measurement value in the ascending orbit and the measurement value in the descending orbit. Therefore, the measurement target point after the change exists somewhere on the circle 34. However, it is not possible to specify the actual measurement target point by itself.
[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 inclination line direction 36 and the circle 34 is obtained from the measurement target point P1 before the change, and this point P2 indicates the actual measurement target point after the change. In this way, it is possible to measure the actual fluctuation amount along the ground surface of the measurement target point.
[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-type base body 20 having a rigid structure standing on the ground surface and a pair of tower-type base bodies 20 installed in opposite directions to each other. A structure composed of the radio wave reflectors 12a and 12b is desirable. The tower base 20 is preferably coated with a radio wave absorber or affixed to the pillars 22 and the shelves 24a and 24b, so that the tower type substrate 20 can be separated from parts other than the radio wave reflectors 12a and 12b. Unnecessary reflected waves can be suppressed. In addition, by attaching both radio wave reflectors 12a and 12b to the middle and upper stages, they can be installed at the same point, and by adjusting the height, observation can be continued throughout the four seasons without being buried in snow. It becomes.
[0024]
The radio wave reflector may be a radio wave reflector, but as shown in FIG. 3, three right-angled isosceles triangular metal plates 40 are combined in a triangular pyramid shape, and the right-angle portions of these metal plates 40 are joined by welding or the like. The corner reflector 42 is preferred. At the time of installation, the orientation and the like are adjusted so that the opening faces the position of the artificial satellite. As a result, the synthetic aperture radar radio wave can be returned toward the artificial satellite in a direction corresponding to the incident direction.
[0025]
Such a corner reflector 42 may cover the opening with a radio wave transmitting plate 44 as shown in FIG. 3A in order to prevent rain or snow from entering the opening. In that case, it is preferable that the plate material 44 is made of a water-repellent material or a water-repellent film is provided on the surface of the plate material 44 so that rain or snow does not adhere to the surface of the radio wave permeable plate material 44. . When the opening is not covered, as shown in FIG. 3B, a small-diameter drain hole 46 is formed in the bottom bottom of the corner reflector 42 so that rainwater and the like that have entered from the opening can be quickly discharged. preferable.
[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 SYMBOLS 10 Ground surface 12a, 12b used as measurement object Radio wave reflector 14 Radio wave reflection structure 16 Artificial satellite 18 Synthetic aperture radar 20 Tower type base | substrate

Claims (5)

飛翔体に搭載されている1台の合成開口レーダから発射されるパルス状の電波の地表面反射波を合成開口干渉処理することにより地表面の変動量を求める方法において、
計測対象となる地表面に剛構造の基体を設け、該基体に一対の電波反射体を、一方は昇交軌道上の合成開口レーダが発する電波を反射する向きに、他方は降交軌道上の合成開口レーダが発する電波を反射する向きに設置して、昇交軌道にて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.
地盤構造上特定される地表面の変動方向の条件は、地滑り変動が斜面の最大傾斜線に沿って発生するという条件であり、その条件を加えることで地滑り変動量を求める請求項1記載の合成開口レーダを用いる地表面変動量計測方法。The composition according to claim 1, wherein the condition of the direction of fluctuation of the ground surface specified by the ground structure is a condition that the landslide fluctuation occurs along the maximum slope line of the slope, and the amount of landslide fluctuation is obtained by adding the condition. Ground surface fluctuation measurement method using aperture radar. 請求項1又は2記載の合成開口レーダを用いた地表面変動量計測方法で用いる電波反射構造物であって、
地表面に立設される剛構造のタワー型基体と、該基体に互いに逆向きに設置される一対の電波反射体とを具備し、前記基体は電波吸収体で覆われ、電波反射体は積雪に埋没しない高さに調整された中段と上段にそれぞれ設置されている地表面変動量計測用の電波反射構造物。
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.
電波反射体は、3枚の直角二等辺三角形金属板を組み合わせ直角部分を互いに接合することで三角錐状にしたコーナリフレクタからなり、該コーナリフレクタの開口部が撥水性を呈する電波透過体で塞がれている請求項3記載の地表面変動量計測用の電波反射構造物。The radio wave reflector consists of a corner reflector made by combining three right-angled isosceles triangular metal plates and joining the right-angled parts to each other, and the opening of the corner reflector is closed with a radio wave-transmitting material exhibiting water repellency. The radio wave reflecting structure for measuring the amount of ground surface fluctuation according to claim 3, wherein 電波反射体は、3枚の直角二等辺三角形金属板を組み合わせ直角部分を互いに接合することで三角錐状にしたコーナリフレクタからなり、該コーナリフレクタの底部に排水穴が形成されている請求項3記載の地表面変動量計測用の電波反射構造物。The radio wave reflector is composed of a corner reflector formed by combining three right-angled isosceles triangular metal plates and joining the right-angle portions to each other, and a drain hole is formed at the bottom of the corner reflector. The radio wave reflecting structure for measuring the ground surface fluctuation amount described.
JP2003099501A 2003-04-02 2003-04-02 Ground surface fluctuation measurement method using synthetic aperture radar Expired - Lifetime JP3987451B2 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Cited By (2)

* Cited by examiner, † Cited by third party
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&#39;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