CN113064170A - Expansive soil area surface deformation monitoring method based on time sequence InSAR technology - Google Patents

Expansive soil area surface deformation monitoring method based on time sequence InSAR technology Download PDF

Info

Publication number
CN113064170A
CN113064170A CN202110335011.5A CN202110335011A CN113064170A CN 113064170 A CN113064170 A CN 113064170A CN 202110335011 A CN202110335011 A CN 202110335011A CN 113064170 A CN113064170 A CN 113064170A
Authority
CN
China
Prior art keywords
dem
time sequence
deformation
expansive soil
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
CN202110335011.5A
Other languages
Chinese (zh)
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.)
Changan University
Original Assignee
Changan University
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 Changan University filed Critical Changan University
Priority to CN202110335011.5A priority Critical patent/CN113064170A/en
Publication of CN113064170A publication Critical patent/CN113064170A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9023SAR image post-processing techniques combined with interferometric techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/885Radar or analogous systems specially adapted for specific applications for ground probing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a swelling soil area surface deformation monitoring method based on a time sequence I nSAR technology, which relates to the technical field of deformation monitoring. After a DEM (digital elevation model) product is updated iteratively, deformation monitoring is carried out on the Ankang airport by adopting an SBAS-I nSAR technology, the depth distribution of the filling expanded soil layer of the expanded soil area is obtained according to DEM error information and is used as a deformation auxiliary judgment means, and large-range and high-precision periodic time sequence deformation information of the expanded soil area is obtained. The used multi-source remote sensing data is simple, convenient and quick to obtain, the efficiency is high, and the cost of manpower and material resources is extremely low.

Description

Expansive soil area surface deformation monitoring method based on time sequence InSAR technology
Technical Field
The invention relates to the technical field of deformation monitoring, in particular to a method for monitoring surface deformation of an expansive soil area based on a time sequence InSAR technology.
Background
In recent years, the monitoring of surface deformation in expansive soil areas has been a hot topic due to the special hydraulic action characteristics of expansive soil bodies and the wide distribution of expansive soil. At present, GNSS measurement, leveling measurement and layered calibration technologies exist in earth surface deformation monitoring methods for expansive soil areas, the above technical schemes can obtain discrete point deformation with higher precision, but the method is limited by terrain factors due to higher labor and equipment cost, and large-area overall deformation distribution of the expansive soil areas is not easy to obtain, so that a high-precision, large-range and long-period deformation monitoring technology is needed.
Aiming at the problems in the prior art, the method for monitoring the deformation of the earth surface of the expansive soil area based on the time sequence InSAR technology is provided, and based on the remote sensing data acquired by multiple sources and multiple sensors, the external DEM data and the time sequence InSAR technology, the deep distribution and time sequence deformation information of the expansive soil layer filled in the expansive soil area are acquired, so that the deformation monitoring with low cost, high efficiency, high precision and large range is realized.
Disclosure of Invention
The invention aims to provide a method for monitoring surface deformation of an expansive soil area based on a time sequence InSAR technology, which is used for acquiring depth distribution and time sequence deformation information of a filling expansive soil layer of the expansive soil area based on remote sensing data acquired by a plurality of sources and sensors, external DEM data and the time sequence InSAR technology, and realizing deformation monitoring with low cost, high efficiency, high precision and large range.
The invention provides a method for monitoring surface deformation of an expansive soil area based on a time sequence InSAR technology, which comprises the following steps:
step S1: preprocessing SAR image data covering a research area by using precision rail data and a digital elevation model DEM;
step S2: selecting SAR image data with space-time baselines meeting a threshold value from a preprocessing result to generate an interference pair, resolving an error between DEM data and real-time terrain simulated by the SAR image data based on a least square criterion, and iteratively updating the obtained DEM error into original DEM data to form new DEM data;
step S3: and resolving the time sequence deformation distribution condition of the research area by using an SBAS-InSAR technology for the SAR image data and the new DEM data to obtain an earth surface deformation result of the expansive soil area.
Further, the preprocessing step in step S1 specifically includes:
step S101: selecting a reference image in the SAR image data;
step S102: removing the Doppler effect of the reference image based on the fine track data;
step S103: and registering and cutting the DEM and the reference image, and outputting a preprocessing result.
Further, the DEM error solution process in step S2 is as follows:
step S201: calculating the offset between interference pairs of the preprocessed data;
step S202: generating a differential interference pattern;
step S203: phase unwrapping the interference based on the MCF minimum cost stream;
step S204: calculating a refined base line and a DEM error;
step S205: and judging whether the DEM has errors, if so, returning to the step S201 to correct the errors of the DEM again, and if not, substituting the DEM data after iterative updating into the step S3 to perform time sequence deformation.
Further, the surface deformation calculation process based on SBAS-InSAR in step S3 includes:
step S301: calculating DEM offset of the interference pairs which accord with the space-time baseline threshold value to generate a differential interference graph;
step S302: calculating a coherence coefficient, and performing band-pass filtering according to the coherence coefficient;
step S303: performing phase unwrapping based on the MCF;
step S304: calculating a refined baseline, secondary differential interference and secondary filtering;
step S305: after filtering, phase unwrapping is carried out again;
step S306: removing the trend error;
step S307: calculating the time sequence deformation quantity;
step S308: and carrying out geocoding to obtain the time sequence deformation distribution condition.
Compared with the prior art, the invention has the following remarkable advantages:
the invention provides a method for monitoring surface deformation of an expansive soil area based on a time sequence InSAR technology. And after the DEM product is updated iteratively, deformation monitoring is carried out on the health airport by adopting a time sequence InSAR technology, and large-range and high-precision periodic time sequence deformation information of the expansive soil area is obtained. The used multi-source remote sensing data is simple, convenient and quick to obtain, the efficiency is high, and the cost of manpower and material resources is extremely low.
The invention provides a method for monitoring surface deformation of an expansive soil area based on a time sequence InSAR technology, which can provide a centimeter-level or even millimeter-level high-precision monitoring means for the expansive soil area. The method has the advantages that the high-precision iterative calculation result of the DEM deviation is overcome, the DEM error influence caused by construction in the expansive soil area is overcome, the result precision can be obviously improved, the deep distribution of the filling expansive soil layer in the expansive soil filling area can be obtained by calculating the DEM error, and the deep information of the expansive soil layer is used as an auxiliary judgment reference for the correctness of the deformation result. The large-range and long-period continuous monitoring of the expansive soil area is realized, deformation information is not omitted, and the analysis and interpretation of the whole deformation mechanism are facilitated.
Drawings
FIG. 1 is a flow chart of a monitoring technique provided by an embodiment of the present invention;
FIG. 2 is a DEM error distribution diagram provided by an embodiment of the invention;
FIG. 3 is a graph of cumulative sedimentation provided by an embodiment of the present invention;
fig. 4 is a diagram of a distribution of time-series deformation of feature points according to an embodiment of the present invention.
Detailed Description
The technical solutions of the embodiments of the present invention are clearly and completely described below with reference to the drawings in the present invention, and it is obvious that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, shall fall within the scope of protection of the present invention.
Referring to fig. 1-4, the invention provides a method for monitoring surface deformation of an expansive soil area based on a time sequence InSAR technology, which comprises the following steps:
step S1: utilizing the European Bureau to provide Sentinrl-1A precision rail data and a digital elevation model DEM to preprocess SAR image data provided by Sentinrl-1A covering a research area;
step S2: selecting SAR image data with space-time baselines meeting a threshold value from a preprocessing result to generate an interference pair, resolving an error between DEM data and real-time terrain simulated by the SAR image data based on a least square criterion, and iteratively updating the obtained DEM error into original DEM data to form new DEM data;
step S3: and (4) solving the time sequence deformation distribution condition of the research area by using an SBAS-InSAR technology for the SAR image data and the DEM data updated in the step S2 to obtain the earth surface deformation result of the expansive soil area.
Wherein the preprocessing step in step S1 specifically includes:
step S101: selecting a reference image in the SAR image data;
step S102: removing the Doppler effect of the reference image based on the fine track data;
step S103: and registering and cutting the DEM and the reference image, wherein the Insar technology can simulate the DEM, the external DEM generally refers to a DEM product provided by the existing mechanism, and a preprocessing result is output.
The DEM error resolving process in step S2 is as follows:
step S201: calculating the offset between interference pairs of the preprocessed data;
step S202: generating a differential interference pattern;
step S203: phase unwrapping the interference based on the MCF minimum cost stream;
step S204: calculating a refined base line and a DEM error;
step S205: and judging whether the DEM has errors, if so, returning to the step S201 to correct the errors of the DEM again, and if not, substituting the DEM data after iterative updating into the step S3 to perform time sequence deformation.
The surface deformation calculation process based on the SBAS-InSAR technology in step S3 includes:
step S301: calculating the offset of the interference pair according with the space-time baseline threshold value to generate a differential interference pattern;
step S302: calculating a coherence coefficient, and performing band-pass filtering according to the coherence coefficient;
step S303: performing phase unwrapping based on the MCF;
step S304: calculating a refined baseline, secondary differential interference and secondary filtering;
step S305: after filtering, phase unwrapping is carried out again;
step S306: removing the trend error;
step S307: calculating the time sequence deformation quantity;
step S308: and carrying out geocoding to obtain the time sequence deformation distribution condition.
Fig. 2 is a Digital Elevation Model (DEM) error distribution diagram of an example area (an expansive soil airport) calculated by the technical scheme, and multiband SAR (synthetic aperture radar) image data is adopted to obtain continuous deformation information with a larger range and high precision, calculate DEM error distribution, obtain expansive soil layer depth information of a filling area of the expansive soil filling area, use the expansive soil layer depth information as auxiliary judgment reference information for deformation result correctness, eliminate errors caused by the DEM, and have higher result reliability. The method is more suitable for monitoring the deformation of the expansive soil area, particularly the expansive soil area under the background of the excavation and filling engineering. As can be seen from FIG. 2, the DEM error distribution obtained by the small baseline set method is consistent with the excavation and filling data during construction.
FIG. 3 is a diagram of a distribution of accumulated deformation of the surface of an example area (an expansive soil airport) calculated by the technical scheme, wherein the distribution of the surface deformation of the example area can be clearly shown in FIG. 3, the maximum deformation of the area appears in the influence areas of expansive soil slopes on two sides of an airport runway, the maximum deformation reaches-11.6 mm, the distribution rule of the deformation of the expansive soil is matched with the distribution of fill described by DEM errors in FIG. 2, and the accuracy of the deformation result is verified.
Fig. 4 is a time sequence distribution diagram of deformation of feature points extracted from an example region (expansive soil airport) solved by the technical scheme, and the regularity of the deformation of expansive soil in the expansive soil filling region can be seen from fig. 4. As shown in fig. 4, deformation of each point in the fill area of the four types of areas has a certain rule: (1) the deformation trends are similar, and are the phenomenon of rising after deformation to the LOS direction; (2) the deformation and the rebound phenomenon of each point occur in similar time and all occur in rainy seasons with concentrated rainfall; (3) and (4) combining the DEM error estimation result to find that the deformation degree of the DEM error estimation result is in direct proportion to the filling amount.
The above disclosure is only for a few specific embodiments of the present invention, however, the present invention is not limited to the above embodiments, and any variations that can be made by those skilled in the art are intended to fall within the scope of the present invention.

Claims (4)

1. A method for monitoring surface deformation of an expansive soil area based on a time sequence InSAR technology is characterized by comprising the following steps:
step S1: preprocessing SAR image data covering a research area by using precision rail data and a digital elevation model DEM;
step S2: selecting SAR image data with space-time baselines meeting a threshold value from a preprocessing result to generate an interference pair, resolving an error between DEM data and real-time terrain simulated by the SAR image data based on a least square criterion, and iteratively updating the obtained DEM error into original DEM data to form new DEM data;
step S3: and (4) solving the time sequence deformation distribution condition of the research area by using an SBAS-InSAR technology for the SAR image data and the DEM data updated in the step S2 to obtain the earth surface deformation result of the expansive soil area.
2. The method for monitoring surface deformation of expansive soil area based on time sequence InSAR technology as claimed in claim 1, wherein the preprocessing step in step S1 specifically comprises:
step S101: selecting a reference image in the SAR image data;
step S102: removing the Doppler effect of the reference image based on the fine track data;
step S103: and registering and cutting the DEM and the reference image, and outputting a preprocessing result.
3. The method for monitoring surface deformation of expansive soil area based on time sequence InSAR technology as claimed in claim 1, wherein the DEM error solution process in step S2 is as follows:
step S201: calculating the offset between interference pairs of the preprocessed data;
step S202: generating a differential interference pattern;
step S203: phase unwrapping the interference based on the MCF minimum cost stream;
step S204: calculating a refined base line and a DEM error;
step S205: and judging whether the DEM has errors, if so, returning to the step S201 to correct the errors of the DEM again, and if not, substituting the DEM data after iterative updating into the step S3 to perform time sequence deformation.
4. The method for monitoring surface deformation of expansive soil area based on time sequence InSAR technology as claimed in claim 1, wherein the surface deformation solving process based on SBAS-InSAR technology in step S3 is:
step S301: calculating the offset of the interference pairs which accord with the space-time baseline threshold value to generate a differential interference pattern;
step S302: calculating a coherence coefficient, and performing band-pass filtering according to the coherence coefficient;
step S303: performing phase unwrapping based on the MCF;
step S304: calculating a refined baseline, secondary differential interference and secondary filtering;
step S305: after filtering, phase unwrapping is carried out again;
step S306: removing the trend error;
step S307: calculating the time sequence deformation quantity;
step S308: and carrying out geocoding to obtain the time sequence deformation distribution condition.
CN202110335011.5A 2021-03-29 2021-03-29 Expansive soil area surface deformation monitoring method based on time sequence InSAR technology Pending CN113064170A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110335011.5A CN113064170A (en) 2021-03-29 2021-03-29 Expansive soil area surface deformation monitoring method based on time sequence InSAR technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110335011.5A CN113064170A (en) 2021-03-29 2021-03-29 Expansive soil area surface deformation monitoring method based on time sequence InSAR technology

Publications (1)

Publication Number Publication Date
CN113064170A true CN113064170A (en) 2021-07-02

Family

ID=76564517

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110335011.5A Pending CN113064170A (en) 2021-03-29 2021-03-29 Expansive soil area surface deformation monitoring method based on time sequence InSAR technology

Country Status (1)

Country Link
CN (1) CN113064170A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113848551A (en) * 2021-09-24 2021-12-28 成都理工大学 Landslide depth inversion method using InSAR lifting rail deformation data
CN115205345A (en) * 2022-02-23 2022-10-18 西南交通大学 Land surface non-uniform deformation information extraction and management method for airport in sea reclamation area
CN115629384A (en) * 2022-12-08 2023-01-20 中南大学 Correction method of time sequence InSAR error and related equipment
CN116051620A (en) * 2023-04-03 2023-05-02 之江实验室 Method and system for estimating active layer thickness of frozen soil area based on InSAR technology

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110888130A (en) * 2019-10-30 2020-03-17 华东师范大学 Coal mine area ground surface deformation monitoring method based on lifting rail time sequence InSAR
CN111323776A (en) * 2020-02-27 2020-06-23 长沙理工大学 Method for monitoring deformation of mining area

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110888130A (en) * 2019-10-30 2020-03-17 华东师范大学 Coal mine area ground surface deformation monitoring method based on lifting rail time sequence InSAR
CN111323776A (en) * 2020-02-27 2020-06-23 长沙理工大学 Method for monitoring deformation of mining area

Non-Patent Citations (16)

* Cited by examiner, † Cited by third party
Title
刘广全: "基于SBAS-InSAR的丹巴县滑坡探测与监测" *
刘晓龙;张永红;李英会;孙广通;吴宏安;: "基于时间序列TerraSAR-X影像的线状地物沉降监测试验研究", 遥感信息, no. 04 *
卢亿雷: "合成孔径雷达卫星精密轨道对干涉测量的影响与分析" *
史绪国;张路;许强;赵宽耀;董杰;蒋厚军;廖明生;: "黄土台塬滑坡变形的时序InSAR监测分析", 武汉大学学报(信息科学版), no. 07 *
孟秀军: "PS-InSAR技术在西秦岭北缘断裂带地壳微小形变监测中的应用" *
张文婷: "基于InSAR技术的DEM重建理论与应用研究" *
朱仁义: "宽幅InSAR技术在地质灾害的综合形变监测应用研究" *
朱文峰;赵超英;张勤;康亚;: "利用InSAR识别与监测黑方台黄土滑坡", 测绘科学, no. 05 *
朱洪波: "基于时序InSAR技术的输电走廊地形形变提取方法研究" *
李伟轩: "时序InSAR技术的地面沉降监测应用研究" *
王霞迎: "地表形变时间序列InSAR监测法综述" *
蒋廷臣: "星载宽幅合成孔径雷达干涉测量形变监测理论与应用研究" *
袁悦: "基于SBAS-InSAR技术的海口地区地面沉降监测及机理分析" *
赵伟颖: "面向矿区沉降监测的InSAR技术应用研究" *
赵超英: "甘肃黑方台黄土滑坡InSAR识别、监测与失稳模式研究", pages 1 - 2 *
邵倩倩: "基于InSAR时序分析技术的任丘地区形变监测研究" *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113848551A (en) * 2021-09-24 2021-12-28 成都理工大学 Landslide depth inversion method using InSAR lifting rail deformation data
CN115205345A (en) * 2022-02-23 2022-10-18 西南交通大学 Land surface non-uniform deformation information extraction and management method for airport in sea reclamation area
CN115629384A (en) * 2022-12-08 2023-01-20 中南大学 Correction method of time sequence InSAR error and related equipment
CN116051620A (en) * 2023-04-03 2023-05-02 之江实验室 Method and system for estimating active layer thickness of frozen soil area based on InSAR technology

Similar Documents

Publication Publication Date Title
CN113064170A (en) Expansive soil area surface deformation monitoring method based on time sequence InSAR technology
Giordan et al. Morphological and kinematic evolution of a large earthflow: The Montaguto landslide, southern Italy
Yang et al. SRTM DEM and its application advances
Patro et al. Hydrodynamic modelling of a large flood‐prone river system in India with limited data
Tomás et al. Monitoring an earthfill dam using differential SAR interferometry: La Pedrera dam, Alicante, Spain
Shah et al. Modelling the effects of spatial variability in rainfall on catchment response. 1. Formulation and calibration of a stochastic rainfall field model
Sanabria et al. Subsidence activity maps derived from DInSAR data: Orihuela case study
CN104133996A (en) Ground settlement risk grade evaluation method based on cloud model and data field
Wang et al. Monitoring of large-scale landslides in Zongling, Guizhou, China, with improved distributed scatterer interferometric SAR time series methods
Yin et al. Long-term flood-hazard modeling for coastal areas using InSAR measurements and a hydrodynamic model: The case study of Lingang New City, Shanghai
Diao et al. Assessment of mining-induced damage to structures using InSAR time series analysis: a case study of Jiulong Mine, China
Liu et al. Urban subsidence monitoring by SBAS-InSAR technique with multi-platform SAR images: A case study of Beijing Plain, China
Uzielli et al. Risk analysis for the Ancona landslide—I: characterization of landslide kinematics
CN113251914A (en) Surface deformation prediction method combining InSAR technology and long-term memory neural network
Wang et al. Three-dimensional slope stability analysis using laser scanning and numerical simulation
CN114812496A (en) Regional ground settlement early warning method based on multi-source heterogeneous data
Nogherotto et al. An integrated hydrological and hydraulic modelling approach for the flood risk assessment over Po river basin
Engelbrecht et al. Detection and monitoring of surface subsidence associated with mining activities in the Witbank Coalfields, South Africa, using differential radar interferometry
Guo et al. Utilization of 3D laser scanning for stability evaluation and deformation monitoring of landslides
Li et al. Evolution of spatiotemporal ground deformation over 30 years in Xi’an, China, with multi-sensor SAR interferometry
Govil et al. Deformation monitoring of Surakachhar underground coal mines of Korba, India using SAR interferometry
Huang et al. Ancient landslide in Wanzhou District analysis from 2015 to 2018 based on ALOS-2 data by QPS-InSAR
An et al. Ground subsidence monitoring in based on UAV-LiDAR technology: a case study of a mine in the Ordos, China
Wu et al. Understanding the Slow Motion of the Wangjiashan Landslide in the Baihetan Reservoir Area (China) from Space‐Borne Radar Observations
CN114912075A (en) Soil humidity data missing filling method fusing spatial-temporal three-dimensional information of site and satellite remote sensing observation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20210702