JP2017166949A - Terrain variation determination method and terrain variation determination system - Google Patents

Terrain variation determination method and terrain variation determination system Download PDF

Info

Publication number
JP2017166949A
JP2017166949A JP2016051994A JP2016051994A JP2017166949A JP 2017166949 A JP2017166949 A JP 2017166949A JP 2016051994 A JP2016051994 A JP 2016051994A JP 2016051994 A JP2016051994 A JP 2016051994A JP 2017166949 A JP2017166949 A JP 2017166949A
Authority
JP
Japan
Prior art keywords
terrain
sar image
microwave reflection
determination method
image data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2016051994A
Other languages
Japanese (ja)
Inventor
庄一郎 内山
Shoichiro Uchiyama
庄一郎 内山
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.)
National Research Institute for Earth Science and Disaster Prevention (NIED)
Original Assignee
National Research Institute for Earth Science and Disaster Prevention (NIED)
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 National Research Institute for Earth Science and Disaster Prevention (NIED) filed Critical National Research Institute for Earth Science and Disaster Prevention (NIED)
Priority to JP2016051994A priority Critical patent/JP2017166949A/en
Publication of JP2017166949A publication Critical patent/JP2017166949A/en
Pending legal-status Critical Current

Links

Images

Abstract

PROBLEM TO BE SOLVED: To provide a terrain variation determination method capable of simply determining a variation state of a terrain without requiring a labor or cost.SOLUTION: The terrain variation determination method includes: a database preparation step of storing an observation object point in which a microwave reflection device array is installed, in association with coordinates; an observation SAR image database acquisition step (S101); an attention point selection step (S103); a coordinate acquisition step (S104) of referring to the database and acquiring coordinates of a selected attention point; a pixel specification step (S105) of specifying pixels in the observation SAR image data with respect to the acquired coordinates; a reflection intensity calculation step (S106) of calculating reflection intensities of the specified pixels; and a determination step (S107) of determining presence/absence of terrain variation at the attention point on the basis of the calculated reflection intensities.SELECTED DRAWING: Figure 6

Description

本発明は、合成開口レーダーを搭載した飛翔体から発射されるマイクロ波を反射するマイクロ波反射デバイスアレイを用いた地形変動判定方法、及びこのような地形変動判定方法を実行する地形変動判定システムに関する。   The present invention relates to a landform change determination method using a microwave reflection device array that reflects microwaves emitted from a flying object equipped with a synthetic aperture radar, and a landform change determination system that executes such a landform change determination method. .

ALOS2や航空機に搭載された合成開口レーダー(SAR;Synthetic Aperture Radar)は、人工衛星や航空機などの飛翔体に搭載したデバイスから発射され、地表付近で散乱したマイクロ波を観測することにより、地表面の状態を調査するものである。合成開口レーダーで取得された観測SAR画像データは、画像ピクセルの反射強度を用いた土地分類(森林、農耕地、水域、等)や二時期のデータを用いた干渉SARによる変動量測定に用いられる。   Synthetic Aperture Radar (SAR) mounted on ALOS2 and aircraft observes microwaves emitted from devices mounted on flying objects such as artificial satellites and aircraft and scattered near the surface of the earth. Is to investigate the state of Observed SAR image data acquired by the synthetic aperture radar is used for land classification (forest, cultivated land, water area, etc.) using the reflection intensity of the image pixel, and measurement of fluctuation amount by interference SAR using two-time data. .

例えば、特許文献1(特開2004−309178号公報)には、飛翔体16に搭載されている1台の合成開口レーダー18から発射されるパルス状の電波の地表面反射波を合成開口干渉処理することにより地表面の変動量を求める方法であって、計測対象となる地表面10に剛構造のタワー型基体20を設け、一方の電波反射体12aは昇交軌道上の、他方の電波反射体12bは降交軌道上の合成開口レーダー電波を反射するように設置して、昇交軌道及び降交軌道にて2回ずつ取得した反射波同士を合成開口干渉処理することで昇交軌道と降交軌道から測定対象地点までの距離をそれぞれ求め、更に地盤構造上特定される地表面の変動方向の条件を加えることで計測対象地点の過去位置からの実際の変動量を求めることが開示されている。
特開2004−309178号公報
For example, in Patent Document 1 (Japanese Patent Laid-Open No. 2004-309178), a synthetic aperture interference process is performed on a ground-surface reflected wave of a pulsed radio wave emitted from one synthetic aperture radar 18 mounted on a flying object 16. In this method, the fluctuation amount of the ground surface is obtained by providing a tower-type base body 20 having a rigid structure on the ground surface 10 to be measured, and one of the radio wave reflectors 12a is reflected on the other side of the ascending orbit. The body 12b is installed so as to reflect the synthetic aperture radar radio wave on the descending orbit, and the reflected waves obtained twice each in the ascending orbit and descending orbit are subjected to the synthetic aperture interference process to obtain the ascending orbit. It is disclosed that the actual fluctuation amount from the past position of the measurement target point is calculated by obtaining the distance from the descending trajectory to the measurement target point and adding the conditions of the ground surface fluctuation direction specified by the ground structure. ing
JP 2004-309178 A

特許文献1記載の技術によって地形変動を求めようとすると、干渉SAR画像データの取得のために複数回のSAR観測を行う必要があり、手間や時間がかかる、という問題があった。さらに、特許文献1記載の技術では、電波反射体を装備した、剛構造を有するタワー型の基体を建設する必要があり、大きなコストを要する。このため、地表の移動量の精密な測定を行うには適しているが、広域の山地斜面に広く設置することは難しい、という問題があった。   If the terrain variation is to be obtained by the technique described in Patent Document 1, it is necessary to perform SAR observation a plurality of times in order to acquire interferometric SAR image data, and there is a problem that it takes time and effort. Furthermore, in the technique described in Patent Document 1, it is necessary to construct a tower-type base body having a rigid structure equipped with a radio wave reflector, which requires a large cost. Therefore, although it is suitable for precise measurement of the amount of movement of the ground surface, there is a problem that it is difficult to install widely on a wide mountain slope.

この発明は、上記のような課題を解決するものであって、本発明に係る地形変動判定方法は、マイクロ波反射デバイスアレイを設置した観測対象ポイントと、当該ポイントの座標とを関連づけて記憶するデータベースを準備するデータベース準備ステップと、観測されたSAR画像データを取得する観測SAR画像データ取得ステップと、着目するポイントを選択する着目ポイント選択ステップと、前記データベース準備ステップで準備されたデータベースを参照し、前記着目ポイント選択ステップで選択された着目ポイントの座標を取得する座標取得ステップと、前記座標取得ステップで取得された座標に対する、観測SAR画像データ中のピクセルを特定するピクセル特定ステップと、前記ピクセル特定ステップで特定されたピクセルの反射強度を算出する反射強度算出ステップと、前記反射強度算出ステップで算出された反射強度に基づいて、着目ポイントにおける地形変動の有無を判定する判定ステップと、を有することを特徴とする。   The present invention solves the above-described problems, and the landform change determination method according to the present invention stores the observation target point on which the microwave reflection device array is installed and the coordinates of the point in association with each other. A database preparation step for preparing a database, an observation SAR image data acquisition step for acquiring observed SAR image data, a point of interest selection step for selecting a point of interest, and the database prepared in the database preparation step A coordinate acquisition step for acquiring the coordinates of the target point selected in the target point selection step, a pixel specifying step for specifying a pixel in the observed SAR image data with respect to the coordinates acquired in the coordinate acquisition step, and the pixel Of the pixel identified in the particular step And reflection intensity calculation step of calculating the elevation strength, based on the reflection intensity calculated by the reflected intensity calculating step, characterized by having a, a determining step of determining the presence or absence of terrain variations in the target point.

また、本発明に係る地形変動判定方法は、前記マイクロ波反射デバイスアレイは、複数の単位マイクロ波反射デバイスからなることを特徴とする。   In the terrain variation determination method according to the present invention, the microwave reflection device array includes a plurality of unit microwave reflection devices.

また、本発明に係る地形変動判定方法は、前記マイクロ波反射デバイスアレイは、観測SAR画像データにおいて所定の反射強度を有するピクセルとして現れることを特徴とする。   The landform change determination method according to the present invention is characterized in that the microwave reflection device array appears as a pixel having a predetermined reflection intensity in the observed SAR image data.

また、本発明に係る地形変動判定システムは、前記いずれかに記載の地形変動判定方法を実行する地形変動判定システムである。   The terrain variation determination system according to the present invention is a terrain variation determination system that executes any one of the terrain variation determination methods described above.

本発明に係る地形変動判定方法及び地形変動判定システムでは、広域の山地斜面などに広く設置することが可能なマイクロ波反射デバイスアレイを用いると共に、干渉SAR画像データを取得する必要もなく、観測SAR画像データを取得した直後に地形変動を判定することができるので、このような本発明に係る地形変動判定方法及び地形変動判定システムによれば、台風や豪雨、雪崩などの大規模な自然現象によって地滑りなどの地形変動が発生したか否か、さらに、地滑りなどの地形変動が発生した地点がどこであるかを、手間や時間やコストをかけることなく、迅速かつ簡便に地形の変動状態を判定することが可能となる。   In the terrain variation determination method and the terrain variation determination system according to the present invention, a microwave reflection device array that can be widely installed on a wide mountain slope or the like is used, and it is not necessary to acquire interference SAR image data. Since the terrain variation can be determined immediately after acquiring the image data, according to the terrain variation determination method and the terrain variation determination system according to the present invention, due to a large-scale natural phenomenon such as a typhoon, heavy rain, and avalanche. Determine whether the landform change such as landslide has occurred, and where the landform change occurred such as landslide is quickly and easily without any effort, time and cost. It becomes possible.

本発明の実施形態に係る地形変動判定方法を実行するシステムを構成するコンピューターの一例を示す図である。It is a figure which shows an example of the computer which comprises the system which performs the landform change determination method which concerns on embodiment of this invention. 本発明の実施形態に係る地形変動判定方法で用いる単位マイクロ波反射デバイス1の一例を示す斜視図である。It is a perspective view which shows an example of the unit microwave reflective device 1 used with the landform change determination method which concerns on embodiment of this invention. 本発明に係る地形変動判定方法の準備フローを示す図である。It is a figure which shows the preparation flow of the landform fluctuation | variation determination method which concerns on this invention. 観測対象ポイントに複数の単位マイクロ波反射デバイス1を配置しマイクロ波反射デバイスアレイ7とした様子を示す図である。It is a figure which shows a mode that the some unit microwave reflection device 1 is arrange | positioned in the observation object point, and it was set as the microwave reflection device array 7. FIG. 本発明の実施形態に係る地形変動判定方法で用いられるデータベース40のデータ構成の一例を示す図である。It is a figure which shows an example of the data structure of the database 40 used with the landform change determination method which concerns on embodiment of this invention. 地形変動判定処理のフローチャートを示す図である。It is a figure which shows the flowchart of a terrain variation determination process. 地形変動判定処理のフローチャートにおけるデータ処理イメージを示す図である。It is a figure which shows the data processing image in the flowchart of a terrain variation determination process. 観測対象ポイントAに設置された複数の単位マイクロ波反射デバイス1が地滑りにより散乱した状態を示す図である。It is a figure which shows the state which the several unit microwave reflective device 1 installed in the observation object point A scattered by the landslide.

以下、本発明の実施の形態を図面を参照しつつ説明する。図1は本発明の実施形態に係る地形変動判定方法を実行するシステムを構成するコンピューターの一例を示す図である。   Hereinafter, embodiments of the present invention will be described with reference to the drawings. FIG. 1 is a diagram illustrating an example of a computer constituting a system for executing a landform change determination method according to an embodiment of the present invention.

図1において、10はシステムバス、11はCPU(Central Processing Unit)、12はRAM(Random Access Memory)、13はROM(Read Only Memory)、14は外部情報機器との通信を司る通信制御部、15はキーボードコントローラなどの入力制御部、16は出力制御部、17は外部記憶装置制御部、18はキーボード、ポインティングデバイス、マウスなどの入力機器からなる入力部、19は印刷装置などの出力部、20はHDD(Hard Disk Drive)等の外部記憶装置、21はグラフィック制御部、22はディスプレイ装置をそれぞれ示している。   In FIG. 1, 10 is a system bus, 11 is a CPU (Central Processing Unit), 12 is a RAM (Random Access Memory), 13 is a ROM (Read Only Memory), 14 is a communication control unit that controls communication with an external information device, 15 is an input control unit such as a keyboard controller, 16 is an output control unit, 17 is an external storage device control unit, 18 is an input unit composed of input devices such as a keyboard, pointing device, and mouse, 19 is an output unit such as a printing device, Reference numeral 20 denotes an external storage device such as an HDD (Hard Disk Drive), 21 denotes a graphic control unit, and 22 denotes a display device.

図1において、CPU11は、ROM13内のプログラム用ROM、或いは、大容量の外部記憶装置20に記憶されたプログラム等に応じて、外部機器と通信することでデータを検索・取得したり、また、図形、イメージ、文字、表等が混在した出力データの処理を実行したり、更に、外部記憶装置20に格納されているデータベースの管理を実行したり、などといった演算処理を行うものである。   In FIG. 1, the CPU 11 retrieves and acquires data by communicating with an external device in accordance with a program ROM stored in the ROM 13 or a program stored in the large-capacity external storage device 20. Processing of output data in which graphics, images, characters, tables, etc. are mixed is executed, and management of a database stored in the external storage device 20 is further executed.

また、CPU11は、システムバス10に接続される各デバイスを統括的に制御する。ROM13内のプログラム用ROMあるいは外部記憶装置20には、CPU11の制御用の基本プログラムであるオペレーティングシステムプログラム(以下OS)等が記憶されている。また、ROM13あるいは外部記憶装置20には出力データ処理等を行う際に使用される各種データが記憶されている。メインメモリーであるRAM12は、CPU11の主メモリ、ワークエリア等として機能する。   Further, the CPU 11 comprehensively controls each device connected to the system bus 10. The program ROM in the ROM 13 or the external storage device 20 stores an operating system program (hereinafter referred to as OS) that is a basic program for controlling the CPU 11. The ROM 13 or the external storage device 20 stores various data used when performing output data processing or the like. The RAM 12 as a main memory functions as a main memory, work area, and the like for the CPU 11.

入力制御部15は、キーボードや不図示のポインティングデバイスからの入力部18を制御する。また、出力制御部16は、プリンタなどの出力部19の出力制御を行う。   The input control unit 15 controls the input unit 18 from a keyboard or a pointing device (not shown). The output control unit 16 performs output control of the output unit 19 such as a printer.

外部記憶装置制御部17は、ブートプログラム、各種のアプリケーション、フォントデータ、ユーザーファイル、編集ファイル、プリンタドライバ等を記憶するHDD(Hard Disk Drive)や、或いはフロッピーディスク(FD)等の外部記憶装置20へのアクセスを制御する。本発明の地形変動判定方法を実現するシステムプログラムは、上記のような外部記憶装置20に記憶されている。また、グラフィック制御部21は、ディスプレイ装置22に表示する情報を描画処理するための構成である。   The external storage device control unit 17 is an external storage device 20 such as an HDD (Hard Disk Drive) or a floppy disk (FD) that stores a boot program, various applications, font data, user files, editing files, printer drivers, and the like. Control access to. A system program for realizing the landform change determination method of the present invention is stored in the external storage device 20 as described above. The graphic control unit 21 is configured to perform drawing processing on information to be displayed on the display device 22.

また、通信制御部14は、ネットワークを介して、外部機器と通信を制御するものであり、これによりシステムが必要とするデータを、インターネットやイントラネット上の外部機器が保有するデータベースから取得したり、外部機器に情報を送信したりすることができるように構成される。   Further, the communication control unit 14 controls communication with an external device via a network, thereby acquiring data required by the system from a database held by an external device on the Internet or an intranet, It is configured to be able to send information to an external device.

外部記憶装置20には、CPU11の制御プログラムであるオペレーティングシステムプログラム(以下OS)以外に、本発明の地形変動判定システムをCPU11上で動作させるシステムプログラム、及びこのシステムプログラムで用いるデータなどがインストールされ保存・記憶されている。   In addition to an operating system program (hereinafter referred to as OS) that is a control program for the CPU 11, the external storage device 20 is installed with a system program for operating the terrain variation determination system of the present invention on the CPU 11, data used in the system program, and the like. Saved and remembered.

本発明の地形変動判定方法を実現するシステムプログラムで利用されるデータとしては、基本的には外部記憶装置20に保存されていることが想定されているが、場合によっては、これらのデータを、通信制御部14を介してインターネットやイントラネット上の外部機器から取得するように構成することも可能である。また、本発明の地形変動判定方法を実現するシステムプログラムで利用されるデータを、USBメモリやCD、DVDなどの各種メディアから取得するように構成することもできる。   The data used in the system program for realizing the landform change determination method of the present invention is basically assumed to be stored in the external storage device 20, but in some cases, these data may be It is also possible to configure to acquire from an external device on the Internet or an intranet via the communication control unit 14. The data used in the system program for realizing the landform change determination method of the present invention can also be obtained from various media such as a USB memory, a CD, and a DVD.

次に、以上のように構成されるコンピューターで実行することが可能な、本発明に係る地形変動判定方法をより詳細に説明する。   Next, the landform change determination method according to the present invention that can be executed by the computer configured as described above will be described in more detail.

まず、本発明に係る地形変動判定方法において、観測対象となる箇所に実際に設置するマイクロ波の反射デバイスについて説明する。マイクロ波反射デバイスは、人工衛星や航空機などの飛翔体に搭載された合成開口レーダーから出射されるマイクロ波を反射するものである。   First, a microwave reflection device actually installed at a location to be observed in the landform change determination method according to the present invention will be described. The microwave reflection device reflects microwaves emitted from a synthetic aperture radar mounted on a flying object such as an artificial satellite or an aircraft.

マイクロ波反射デバイスで反射されたマイクロ波は、合成開口レーダーによって得られる観測SAR画像データの画像の一部に自然地形では通常生じない強さを持った輝点として写り込む。   The microwave reflected by the microwave reflection device is reflected as a bright spot having a strength that does not normally occur in natural terrain in a part of the image of the observed SAR image data obtained by the synthetic aperture radar.

本発明に係る地形変動判定方法においては、マイクロ波反射デバイスとして、複数の単位マイクロ波反射デバイス1をアレイ状に配置したマイクロ波反射デバイスアレイ7を用いることを基本としている。複数の単位マイクロ波反射デバイス1をアレイとして使用することにより、マイクロ波反射デバイスアレイ7は、合成開口レーダーで取得される観測SAR画像データにおいて、自然地形では通常は生じない、所定の反射強度を有するピクセルとして現れるように調整される。マイクロ波反射デバイスアレイ7を構成する単位マイクロ波反射デバイス1の数は、地形変動を監視しようとする観測ポイントにおける地形の斜度(傾斜角度)、開空率(全天を100とした時の空が占める割合)、植生の有無、植生の高さ等によって調整する。   In the landform change determination method according to the present invention, a microwave reflection device array 7 in which a plurality of unit microwave reflection devices 1 are arranged in an array is used as a microwave reflection device. By using a plurality of unit microwave reflection devices 1 as an array, the microwave reflection device array 7 has a predetermined reflection intensity that does not normally occur in natural terrain in observed SAR image data acquired by a synthetic aperture radar. Adjusted to appear as having pixels. The number of unit microwave reflection devices 1 constituting the microwave reflection device array 7 is as follows. The slope of the terrain (inclination angle) and the open space ratio (when the sky is 100) Adjust according to the proportion of the sky), the presence or absence of vegetation, the height of vegetation, etc.

ここで、観測SAR画像データ上により強い反射強度を持った輝点を得るために、マイクロ波反射デバイスアレイ7を設置する敷地の大きさは、観測SAR画像データの地上分解能と同等かそれよりも小さいほうが好ましい。   Here, in order to obtain a bright spot having stronger reflection intensity on the observed SAR image data, the size of the site where the microwave reflection device array 7 is installed is equal to or more than the ground resolution of the observed SAR image data. Smaller is preferable.

ここで、複数の単位マイクロ波反射デバイス1をアレイとして使用することにより、一つまたは複数の単位マイクロ波反射デバイス1が、生物による着巣や、落石による破壊等の、不測の理由によるマイクロ波反射能力の喪失が生じた場合でも、観測SAR画像データ中の観測対象ポイントにおいて自然地形では通常は生じない、所定の反射強度を有するピクセルとして現れる限り、地形変動の有無を判定することができる。   Here, by using a plurality of unit microwave reflection devices 1 as an array, one or a plurality of unit microwave reflection devices 1 may become microwaves due to unforeseen reasons such as nests by living organisms or destruction by falling rocks Even if the loss of reflection capability occurs, it is possible to determine the presence or absence of terrain variation as long as it appears as a pixel having a predetermined reflection intensity that does not normally occur in natural terrain at the observation target point in the observed SAR image data.

ここで、地滑りなどの地形変動によってアレイ中の一個または数個の単位マイクロ波反射デバイス1が失われたり、土砂に埋没したり、破壊されたりしても、一つ以上の単位マイクロ波反射デバイス1が残存していれば、観測SAR画像データにおいて自然地形では通常生じない十分に強い反射強度を有するピクセルが得られる可能性が高くなる。これにより、地形変動後も自然地形では通常生じない十分に強い反射強度を有するピクセルが得られる場合は、地形変動の向きや水平方向の変動量を推定することも可能である。   Here, even if one or several unit microwave reflection devices 1 in the array are lost, buried in the earth, or destroyed due to terrain fluctuation such as landslide, one or more unit microwave reflection devices If 1 remains, there is a high possibility that a pixel having sufficiently strong reflection intensity that does not normally occur in natural terrain in the observed SAR image data can be obtained. As a result, when a pixel having sufficiently strong reflection intensity that does not normally occur in natural terrain after terrain change is obtained, it is also possible to estimate the direction of terrain change and the amount of change in the horizontal direction.

ここで、本発明に係る地形変動判定方法においては、マイクロ波反射デバイスとして、複数の単位マイクロ波反射デバイス1を用いることを一応基本とはしているが、1つの単位マイクロ波反射デバイス1を用いるようにすることもできる。図2に、あらゆる方向から照射されたマイクロ波を反射することを企図したマイクロ波反射デバイスの一例を挙げる。   Here, in the landform change determination method according to the present invention, it is basically assumed that a plurality of unit microwave reflection devices 1 are used as the microwave reflection devices, but one unit microwave reflection device 1 is used. It can also be used. FIG. 2 shows an example of a microwave reflection device intended to reflect microwaves irradiated from all directions.

図2は本発明の実施形態に係る地形変動判定方法で用いる単位マイクロ波反射デバイス1の一例を示す斜視図である。本発明に係る地形変動判定方法においては、図2に示す単位マイクロ波反射デバイス1以外に他の形態のものも用いることができる。   FIG. 2 is a perspective view showing an example of the unit microwave reflection device 1 used in the landform change determination method according to the embodiment of the present invention. In the terrain variation determination method according to the present invention, other forms besides the unit microwave reflection device 1 shown in FIG. 2 can be used.

単位マイクロ波反射デバイス1は、アルミニウムなどの金属板などからなる、複数の反射板3によって構成されている。単位マイクロ波反射デバイス1は、3面の反射面を有するコーナーリフレクターが、8個組み合わされたような形状を有している。すなわち、単位マイクロ波反射デバイス1の上方においては4方に向いて開口している3面コーナーリフレクターが4個配され、下方においては4方に向いて開口している3面コーナーリフレクターが4個配されたような構造を有している。   The unit microwave reflection device 1 includes a plurality of reflection plates 3 made of a metal plate such as aluminum. The unit microwave reflecting device 1 has a shape in which eight corner reflectors having three reflecting surfaces are combined. That is, four three-surface corner reflectors that open in four directions are arranged above the unit microwave reflection device 1, and four three-surface corner reflectors that open in four directions below. It has a structure as arranged.

上記のような単位マイクロ波反射デバイス1は、不図示の基台上に載置されるようにして、観測対象となる場所に複数設置される。また、単位マイクロ波反射デバイス1は、不図示の補強手段などによって、補強されることも好ましい形態である。このような補強手段としては、反射面と、当該反射面と隣接する反射面との間に筋交いを設けることが有効である。また、不図示のマイクロ波を透過するレドーム(覆い)によって雨や土、落ち葉、生物の着巣などによる単位マイクロ波反射デバイス1の汚損を防止することも好ましい設置形態である。   A plurality of the unit microwave reflection devices 1 as described above are installed at a place to be observed so as to be placed on a base (not shown). The unit microwave reflection device 1 is also preferably reinforced by a reinforcing means (not shown). As such a reinforcing means, it is effective to provide a brace between the reflecting surface and the reflecting surface adjacent to the reflecting surface. It is also a preferable installation form to prevent the unit microwave reflection device 1 from being soiled by rain, soil, fallen leaves, living organisms' nest, etc., by a radome (cover) that transmits microwaves (not shown).

次に、上記のような単位マイクロ波反射デバイス1を観測対象となるポイントに実際に設置するなどして、本発明に係る地形変動判定方法を準備する工程について説明する。   Next, a process of preparing the landform change determination method according to the present invention by actually installing the unit microwave reflection device 1 as described above at a point to be observed will be described.

図3は本発明に係る地形変動判定方法の準備フローを示す図である。このような準備フローは人手によって実行されるものである。   FIG. 3 is a diagram showing a preparation flow of the landform change determination method according to the present invention. Such a preparation flow is executed manually.

図3のステップS10で、本発明に係る地形変動判定方法の準備フローが開始されると続いてステップS20に進む。   When the preparation flow of the landform change determination method according to the present invention is started in step S10 of FIG. 3, the process proceeds to step S20.

ステップS20では、観測SAR画像データの少なくとも1つのピクセルにおいて、所定以上の反射強度が得られる個数の単位マイクロ波反射デバイス1を観測対象ポイントに設置する工程を実施する。本発明に係る地形変動判定方法においては、観測対象ポイントに、複数の単位マイクロ波反射デバイス1を設置して、これらからなるマイクロ波反射デバイスアレイ7を構成することが一つの特徴となっている。   In step S20, a step of installing the number of unit microwave reflection devices 1 that can obtain a reflection intensity equal to or higher than a predetermined value at the observation target point in at least one pixel of the observation SAR image data is performed. One feature of the landform change determination method according to the present invention is that a plurality of unit microwave reflection devices 1 are installed at an observation target point, and a microwave reflection device array 7 including these units is configured. .

図4は観測対象ポイントに複数の単位マイクロ波反射デバイス1を配置しマイクロ波反射デバイスアレイ7とした様子を示す図である。図4に示す観測対象ポイントの例では、山の斜面に4つの単位マイクロ波反射デバイス1からなるマイクロ波反射デバイスアレイ7を配置している。ここで、マイクロ波反射デバイスアレイ7を構成するための単位マイクロ波反射デバイス1の数は任意である。ただし、マイクロ波反射デバイスアレイ7は、合成開口レーダーで取得される観測SAR画像データ中の1つのピクセルにおいて、マイクロ波の反射強度値が所定以上となるような数の単位マイクロ波反射デバイス1で構成する。   FIG. 4 is a diagram showing a state in which a plurality of unit microwave reflection devices 1 are arranged at the observation target point to form a microwave reflection device array 7. In the example of the observation target point illustrated in FIG. 4, a microwave reflection device array 7 including four unit microwave reflection devices 1 is disposed on a slope of a mountain. Here, the number of unit microwave reflection devices 1 for constituting the microwave reflection device array 7 is arbitrary. However, the microwave reflection device array 7 is the number of unit microwave reflection devices 1 such that the reflection intensity value of the microwaves is greater than or equal to a predetermined value in one pixel in the observed SAR image data acquired by the synthetic aperture radar. Configure.

ステップS30では、全ての観測ポイントについて、マイクロ波反射デバイスアレイ7を設置したか否かを判定する。この判定がNOであればステップS50に進み、次の観測対象ポイントに着目する。続いてステップS10に戻り、マイクロ波反射デバイスアレイ7の設置を行う。一方、ステップS30の判定がYESであればステップS40に進む。   In step S30, it is determined whether or not the microwave reflection device array 7 has been installed for all observation points. If this judgment is NO, it will progress to Step S50 and will pay attention to the next observation object point. Then, it returns to step S10 and installation of the microwave reflection device array 7 is performed. On the other hand, if the determination in step S30 is YES, the process proceeds to step S40.

ステップS40では、観測対象ポイントと、その緯度、経度を対応づけて記憶したデータベース40を構築する。経度、緯度はGNSS(Global Navigation Satellite System)やトータルステーションなどによって計測する。この他、地形変動の推定に資する以下の諸情報を追加してもよい。例えば、マイクロ波反射デバイスアレイ7を構成するために用いた単位マイクロ波反射デバイス1の個数データと、観測対象ポイントの諸条件として設置した地形の斜度(傾斜角)と、開空率(全天を100とした時の空が占める割合)と、植生の有無と、空を覆う樹高の高い植生の種類とその高さ、が挙げられる。図5はこのようなデータベース40のデータ構成の一例を示す図である。ここに挙げる他に、単位マイクロ波反射デバイスの製造年月日、製造者、コーナーリフレクター開口面積、前回の観測SAR画像データの観測結果など、マイクロ波反射デバイスアレイ7および観測対象ポイントを維持・管理するために有益な情報をデータベースに追加してもよい。   In step S40, the database 40 that stores the observation target point and the latitude and longitude in association with each other is constructed. The longitude and latitude are measured by a GNSS (Global Navigation Satellite System), a total station, or the like. In addition, the following information that contributes to estimation of topographical changes may be added. For example, the number data of the unit microwave reflection devices 1 used to configure the microwave reflection device array 7, the inclination (tilt angle) of the terrain installed as various conditions of the observation target point, and the open space ratio (total The ratio of the sky when the sky is 100), the presence or absence of vegetation, and the type and height of vegetation with a high tree height covering the sky. FIG. 5 is a diagram showing an example of the data structure of such a database 40. In addition to the above, the microwave reflection device array 7 and observation target points are maintained and managed, such as the date of manufacture of the unit microwave reflection device, the manufacturer, the opening area of the corner reflector, and the observation results of the previous observation SAR image data. Useful information may be added to the database.

図5に示す例では、観測対象ポイントであるA、B、C、D、・・・にマイクロ波反射デバイスアレイ7を設置する例を示している。それぞれのポイントにおいて、座標データとなるGNSSデータと、マイクロ波反射デバイスアレイ7を構成するために用いた単位マイクロ波反射デバイス1の数のデータとが関連づけられデータ化されている。   In the example shown in FIG. 5, the example which installs the microwave reflective device array 7 in A, B, C, D, ... which are observation object points is shown. At each point, the GNSS data serving as coordinate data and the data of the number of unit microwave reflection devices 1 used to configure the microwave reflection device array 7 are associated and converted into data.

なお、本実施形態では、GNSSデータとして緯度、経度からなる座標データを用いるようにしたが、他の基準の座標データを用いるようにしてもよい。また、データベース40には、その他のデータ−例えば、ポイントの地名などのデータ−を関連づけて、格納するようにしてもよい。この他、ステップS40の説明で挙げた関連項目の情報も含めることができる。   In the present embodiment, coordinate data including latitude and longitude is used as GNSS data, but other reference coordinate data may be used. Further, the database 40 may store other data in association with each other, for example, data such as place names of points. In addition, the information on the related items mentioned in the description of step S40 can also be included.

ステップS40で準備されたデータベース40は、外部記憶装置20に記憶され、本発明に係る地形変動判定処理に用いる。   The database 40 prepared in step S40 is stored in the external storage device 20, and is used for the landform change determination process according to the present invention.

ステップS60で、本発明に係る地形変動判定方法の準備フローが終了する。   In step S60, the preparation flow of the landform change determination method according to the present invention ends.

次に、以上のような準備フローを経て実行される、本発明に係る地形変動判定処理について説明する。図6は地形変動判定処理のフローチャートを示す図である。このようなフローチャートは、例えば図1に示されるコンピューターによって実行されるものである。   Next, the terrain variation determination process according to the present invention, which is executed through the preparation flow as described above, will be described. FIG. 6 is a flowchart of the terrain variation determination process. Such a flowchart is executed by, for example, the computer shown in FIG.

また、図7は地形変動判定処理のフローチャートにおけるデータ処理イメージを示す図である。図7では着目ポイントとしてポイントAが選定された場合における、ステップS103、S104、S105、S106でのデータ処理イメージを例示している。   FIG. 7 is a diagram showing a data processing image in the flowchart of the landform change determination process. FIG. 7 illustrates data processing images in steps S103, S104, S105, and S106 when point A is selected as the point of interest.

図6において、ステップS100で変動判定処理が開始されると、ステップS101に進み、合成開口レーダーで取得される観測SAR画像データの読み込みが実行される。続くステップS102では、データベース40の読み込みが実行される。   In FIG. 6, when the variation determination process is started in step S100, the process proceeds to step S101, and reading of observation SAR image data acquired by the synthetic aperture radar is executed. In the subsequent step S102, the database 40 is read.

ステップS103では、地形変動判定処理を行おうとする個別の着目ポイントが選定される。(図7参照。)
次にステップS104では、データベース40を参照して、選定された着目ポイントの座標データを取得し、続くステップS105では、当該座標データに対応する観測SAR画像データ中のピクセルを特定し、次のステップS106で、当該ピクセルの反射強度を算出する。ピクセルの反射強度の算出にあたっては、当該ピクセルが持つ反射強度データを参照することにより求めることができる。(いずれのステップも図7参照。)
ここで、ステップS105では、観測SAR画像データ自身のもつ誤差(ピクセルの位置のズレ)を考慮するため、ステップS104で取得した着目ポイントの座標データを中心として半径5ピクセル程度の範囲内に存在するすべてのピクセルの反射強度を算出する。反射強度を算出する範囲の設定は、使用する観測SAR画像の地上分解能や性能によって調整する。観測ポイントの座標と、豪雨や地滑り、雪崩等の地形変化の発生していないことが明らかなSAR観測日の観測SAR画像データから得られる観測ポイントの座標とを比較することで、反射強度を算出する範囲の参考値が得られる。
In step S103, an individual point of interest for which the terrain variation determination process is to be performed is selected. (See Figure 7.)
In step S104, the coordinate data of the selected point of interest is acquired with reference to the database 40. In the subsequent step S105, the pixel in the observed SAR image data corresponding to the coordinate data is specified, and the next step In S106, the reflection intensity of the pixel is calculated. In calculating the reflection intensity of a pixel, it can be obtained by referring to the reflection intensity data of the pixel. (See Figure 7 for all steps.)
Here, in step S105, in order to consider the error (pixel position deviation) of the observed SAR image data itself, it exists within a range of about 5 pixels in radius with the coordinate data of the point of interest acquired in step S104 as the center. Calculate the reflection intensity of all pixels. The setting of the range for calculating the reflection intensity is adjusted according to the ground resolution and performance of the observation SAR image to be used. The reflection intensity is calculated by comparing the coordinates of the observation point with the coordinates of the observation point obtained from the observed SAR image data on the SAR observation day where it is clear that topographic changes such as heavy rain, landslides, and avalanches have not occurred. The reference value of the range to be obtained is obtained.

ステップS107では、算出された反射強度が予め定められた所定値以上であるか否かが判定される。   In step S107, it is determined whether or not the calculated reflection intensity is greater than or equal to a predetermined value.

ステップS107における判定結果がYESであるときには、マイクロ波反射デバイスアレイ7は設置初期と変わらないものと考えられるので、ステップS110に進み、着目ポイントの地形変動は無いものと判定する。   When the determination result in step S107 is YES, it is considered that the microwave reflection device array 7 is not changed from the initial setting, so the process proceeds to step S110, and it is determined that there is no landform change at the point of interest.

一方、ステップS107における判定結果がNOであるときには、マイクロ波反射デバイスアレイ7は設置初期から何かしらの変動があったものと考えられるので、ステップS108に進み、着目ポイントの地形変動が有るものと判定する。着目ポイントの地形変動が有るものと判定された場合は、該当ポイントの観測SAR画像データをグラフィックとして表示したり、出力装置(出力部19)によって印刷したり、通信装置(通信制御部14)によって通知したりすることができる。   On the other hand, when the determination result in step S107 is NO, it is considered that the microwave reflecting device array 7 has undergone some change from the initial installation stage, and therefore, the process proceeds to step S108, where it is determined that there is a landform change at the point of interest. To do. When it is determined that there is a terrain variation at the point of interest, the observed SAR image data of the point is displayed as a graphic, printed by the output device (output unit 19), or by the communication device (communication control unit 14). Can be notified.

図8は観測対象ポイントAに設置された複数の単位マイクロ波反射デバイス1が地滑り等の地形変化により散乱した状態を示す図である。マイクロ波反射デバイスアレイ7を構成する個々の単位マイクロ波反射デバイス1が図8に示すような状態となると、観測SAR画像データ中の該当ピクセルで、所定の反射強度を得ることができなくなる。これにより、ステップS108のような判定を行うことが可能となる。不図示だが、単位マイクロ波反射デバイス1が地滑り等の地形変化によって埋没したり、破壊されたり、失われたりしても図8の状態と同様に、観測SAR画像データ中の該当ピクセルで、所定の反射強度を得ることができなくなる。   FIG. 8 is a diagram illustrating a state in which a plurality of unit microwave reflection devices 1 installed at the observation target point A are scattered due to a terrain change such as a landslide. When the individual unit microwave reflection devices 1 constituting the microwave reflection device array 7 are in a state as shown in FIG. 8, it becomes impossible to obtain a predetermined reflection intensity at the corresponding pixel in the observed SAR image data. Thereby, it is possible to perform the determination as in step S108. Although not shown, even if the unit microwave reflection device 1 is buried, destroyed, or lost due to a terrain change such as a landslide, it is determined at a predetermined pixel in the observed SAR image data as in the state of FIG. The reflection intensity cannot be obtained.

ステップS109では、全ての着目ポイントについて処理が実行されたか否かが判定される。ステップS109での判定がNOであるときには、ステップS111に進み、次の着目ポイントを選定し、ステップS104に戻る。一方、ステップS109での判定がYESであるときには、ステップS112に進み、変動判定処理を終了する。   In step S109, it is determined whether or not processing has been executed for all the points of interest. When the determination in step S109 is NO, the process proceeds to step S111, the next point of interest is selected, and the process returns to step S104. On the other hand, when the determination in step S109 is YES, the process proceeds to step S112 to end the variation determination process.

以上、本発明に係る地形変動判定方法及び地形変動判定システムでは、広域の山地斜面などに広く設置することが可能なマイクロ波反射デバイスアレイ7を用いると共に、干渉SAR画像データを取得する必要もなく、観測SAR画像データを取得した直後に地形変動を判定することができるので、このような本発明に係る地形変動判定方法及び地形変動判定システムによれば、台風や豪雨、雪崩などの大規模な自然現象によって地滑りなどの地形変動が発生したか否か、さらに、地滑りなどの地形変動が発生した地点がどこであるかを、手間や時間やコストをかけることなく、迅速かつ簡便に地形の変動状態を判定することが可能となる。   As described above, in the terrain variation determination method and the terrain variation determination system according to the present invention, the microwave reflection device array 7 that can be widely installed on a wide mountain slope is used, and there is no need to acquire interference SAR image data. Since the terrain variation can be determined immediately after the observation SAR image data is acquired, according to the terrain variation determination method and the terrain variation determination system according to the present invention, a large-scale such as typhoon, heavy rain, avalanche, etc. Whether or not landslides or other terrain changes have occurred due to natural phenomena, and where the landslide or other terrain changes have occurred, the state of terrain changes quickly and easily without any effort, time, or cost. Can be determined.

1・・・単位マイクロ波反射デバイス
3・・・反射板
7・・・マイクロ波反射デバイスアレイ
10・・・システムバス
11・・・CPU(Central Processing Unit)
12・・・RAM(Random Access Memory)
13・・・ROM(Read Only Memory)
14・・・通信制御部
15・・・入力制御部
16・・・出力制御部
17・・・外部記憶装置制御部
18・・・入力部
19・・・出力部
20・・・外部記憶装置
21・・・グラフィック制御部
22・・・ディスプレイ装置
40・・・データベース
DESCRIPTION OF SYMBOLS 1 ... Unit Microwave reflection device 3 ... Reflector 7 ... Microwave reflection device array 10 ... System bus 11 ... CPU (Central Processing Unit)
12 ... RAM (Random Access Memory)
13 ... ROM (Read Only Memory)
DESCRIPTION OF SYMBOLS 14 ... Communication control part 15 ... Input control part 16 ... Output control part 17 ... External storage device control part 18 ... Input part 19 ... Output part 20 ... External storage device 21 ... Graphic control unit 22 ... Display device 40 ... Database

Claims (4)

マイクロ波反射デバイスアレイを設置した観測対象ポイントと、当該ポイントの座標とを関連づけて記憶するデータベースを準備するデータベース準備ステップと、
観測されたSAR画像データを取得する観測SAR画像データ取得ステップと、
着目するポイントを選択する着目ポイント選択ステップと、
前記データベース準備ステップで準備されたデータベースを参照し、前記着目ポイント選択ステップで選択された着目ポイントの座標を取得する座標取得ステップと、
前記座標取得ステップで取得された座標に対する、観測SAR画像データ中のピクセルを特定するピクセル特定ステップと、
前記ピクセル特定ステップで特定されたピクセルの反射強度を算出する反射強度算出ステップと、
前記反射強度算出ステップで算出された反射強度に基づいて、着目ポイントにおける地形変動の有無を判定する判定ステップと、
を有することを特徴とする地形変動判定方法。
A database preparation step of preparing a database for storing the observation target point where the microwave reflection device array is installed and the coordinates of the point in association with each other;
An observation SAR image data acquisition step of acquiring observed SAR image data;
A point of interest selection step for selecting a point of interest;
A coordinate acquisition step of referring to the database prepared in the database preparation step and acquiring the coordinates of the target point selected in the target point selection step;
A pixel specifying step of specifying a pixel in the observed SAR image data with respect to the coordinates acquired in the coordinate acquiring step;
A reflection intensity calculating step of calculating the reflection intensity of the pixel specified in the pixel specifying step;
Based on the reflection intensity calculated in the reflection intensity calculation step, a determination step for determining the presence or absence of terrain fluctuation at the point of interest;
A method for determining terrain variation, comprising:
前記マイクロ波反射デバイスアレイは、複数の単位マイクロ波反射デバイスからなることを特徴とする請求項1に記載の地形変動判定方法。 The terrain fluctuation determination method according to claim 1, wherein the microwave reflection device array includes a plurality of unit microwave reflection devices. 前記マイクロ波反射デバイスアレイは、観測SAR画像データにおいて所定の反射強度を有するピクセルとして現れることを特徴とする請求項1又は請求項2に記載の地形変動判定方法。 The landform change determination method according to claim 1 or 2, wherein the microwave reflection device array appears as a pixel having a predetermined reflection intensity in the observed SAR image data. 前記請求項1乃至請求項3のいずれか1項に記載の地形変動判定方法を実行する地形変動判定システム。 A terrain variation determination system for executing the terrain variation determination method according to any one of claims 1 to 3.
JP2016051994A 2016-03-16 2016-03-16 Terrain variation determination method and terrain variation determination system Pending JP2017166949A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016051994A JP2017166949A (en) 2016-03-16 2016-03-16 Terrain variation determination method and terrain variation determination system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016051994A JP2017166949A (en) 2016-03-16 2016-03-16 Terrain variation determination method and terrain variation determination system

Publications (1)

Publication Number Publication Date
JP2017166949A true JP2017166949A (en) 2017-09-21

Family

ID=59913361

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016051994A Pending JP2017166949A (en) 2016-03-16 2016-03-16 Terrain variation determination method and terrain variation determination system

Country Status (1)

Country Link
JP (1) JP2017166949A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2019215252A (en) * 2018-06-13 2019-12-19 株式会社安藤・間 Dynamic monitoring system, and dynamic monitoring method
CN114677401A (en) * 2022-03-28 2022-06-28 珠江水利委员会珠江水利科学研究院 Water body extraction method and system based on polarization radar self-image features

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02170081A (en) * 1988-10-28 1990-06-29 H Bain John Jr Expandable oval radar reflector
JPH09189762A (en) * 1996-01-08 1997-07-22 Mitsubishi Electric Corp Method for observing variation of ground surface using radar and synthetic aperture radar and transponder therefor
US6466156B1 (en) * 1999-02-26 2002-10-15 Totalforsvarets Forskningsinstitut Method of detecting objects that change with time by means of a SAR radar
JP2003139845A (en) * 2001-11-02 2003-05-14 Omron Corp Invader detector and radio wave reflector
JP2004309178A (en) * 2003-04-02 2004-11-04 Oyo Corp Ground surface variation measuring method using synthetic aperture radar
JP2006349476A (en) * 2005-06-15 2006-12-28 Toyota Motor Corp Radar reflector
JP2011007373A (en) * 2009-06-24 2011-01-13 Ihi Aerospace Co Ltd Folding corner reflector
JP2011214975A (en) * 2010-03-31 2011-10-27 Tokyo Electric Power Co Inc:The Object displacement monitoring device and object displacement monitoring method

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02170081A (en) * 1988-10-28 1990-06-29 H Bain John Jr Expandable oval radar reflector
JPH09189762A (en) * 1996-01-08 1997-07-22 Mitsubishi Electric Corp Method for observing variation of ground surface using radar and synthetic aperture radar and transponder therefor
US6466156B1 (en) * 1999-02-26 2002-10-15 Totalforsvarets Forskningsinstitut Method of detecting objects that change with time by means of a SAR radar
JP2003139845A (en) * 2001-11-02 2003-05-14 Omron Corp Invader detector and radio wave reflector
JP2004309178A (en) * 2003-04-02 2004-11-04 Oyo Corp Ground surface variation measuring method using synthetic aperture radar
JP2006349476A (en) * 2005-06-15 2006-12-28 Toyota Motor Corp Radar reflector
JP2011007373A (en) * 2009-06-24 2011-01-13 Ihi Aerospace Co Ltd Folding corner reflector
JP2011214975A (en) * 2010-03-31 2011-10-27 Tokyo Electric Power Co Inc:The Object displacement monitoring device and object displacement monitoring method

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2019215252A (en) * 2018-06-13 2019-12-19 株式会社安藤・間 Dynamic monitoring system, and dynamic monitoring method
JP7113673B2 (en) 2018-06-13 2022-08-05 株式会社安藤・間 Movement monitoring system and movement monitoring method
CN114677401A (en) * 2022-03-28 2022-06-28 珠江水利委员会珠江水利科学研究院 Water body extraction method and system based on polarization radar self-image features

Similar Documents

Publication Publication Date Title
Wulder et al. Lidar sampling for large-area forest characterization: A review
Santillan et al. Vertical accuracy assessment of 30-M resolution ALOS, ASTER, and SRTM global DEMS over Northeastern Mindanao, Philippines
Miliaresis et al. Vertical accuracy of the SRTM DTED level 1 of Crete
JP6130101B2 (en) Method and system for generating weather and ground reflection information
Goodwin et al. Assessing the repeatability of terrestrial laser scanning for monitoring gully topography: A case study from Aratula, Queensland, Australia
Nyström et al. Prediction of tree biomass in the forest–tundra ecotone using airborne laser scanning
Forsmoo et al. Structure from motion photogrammetry in ecology: Does the choice of software matter?
WO2016136775A1 (en) Laser measuring system, reflection target body, and laser measuring method
JP6745169B2 (en) Laser measuring system and laser measuring method
Altunel Evaluation of TanDEM-X 90 m digital elevation model
JP2017166949A (en) Terrain variation determination method and terrain variation determination system
Singh et al. Moving from plot-based to hillslope-scale assessments of savanna vegetation structure with long-range terrestrial laser scanning (LR-TLS)
Miller et al. A robust surface matching technique for coastal geohazard assessment and management
KR101036823B1 (en) Radar observational density analysis system and analysis method of the same
González et al. Relative height accuracy estimation method for InSAR-based DEMs
Fraser et al. Towards precise drone-based measurement of elevation change in permafrost terrain experiencing thaw and thermokarst
Spore et al. Collection, processing, and accuracy of mobile terrestrial lidar survey data in the coastal environment
Gupta et al. Digital elevation model
Ebaid Accuracy enhancement of SRTM and ASTER DEMs using weight estimation regression model
Díaz Vilariño et al. Scan planning optimization for outdoor archaeological sites
Rusli et al. Accuracy assessment of DEM from UAV and TanDEM-X imagery
Leslar et al. Error analysis of a mobile terrestrial LiDAR system
Soycan et al. Three dimensional modeling of a forested area using an airborne light detection and ranging method
Lubczonek Location determination of radar sensors by using LIDAR data
JP2020165921A (en) Information processing device and control program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180525

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190313

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190320

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190425

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190925

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20200701