CN103047970B - Method for determining damage boundary of subsided farmland - Google Patents

Method for determining damage boundary of subsided farmland Download PDF

Info

Publication number
CN103047970B
CN103047970B CN201210552842.9A CN201210552842A CN103047970B CN 103047970 B CN103047970 B CN 103047970B CN 201210552842 A CN201210552842 A CN 201210552842A CN 103047970 B CN103047970 B CN 103047970B
Authority
CN
China
Prior art keywords
information
boundary
subsidence
area
border
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.)
Active
Application number
CN201210552842.9A
Other languages
Chinese (zh)
Other versions
CN103047970A (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.)
China University of Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology Beijing CUMTB
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 China University of Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN201210552842.9A priority Critical patent/CN103047970B/en
Publication of CN103047970A publication Critical patent/CN103047970A/en
Application granted granted Critical
Publication of CN103047970B publication Critical patent/CN103047970B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention relates to a technical method for determining the damage boundary of a subsided farmland and belongs to the technical field of surveying and mapping as well as environmental monitoring. The method comprises the following steps: selecting corresponding remote sensing data and non-remote sensing data to determine a preliminary area of subsidence; extracting subsided water area boundary information in the subsidence area by using multi-scale wavelet transform edge detection, establishing an index model of subsidence land boundary buffering, buffering the boundary of a subsided water area by using the model and obtaining subsidence boundary information of the area to be surveyed and ground subsidence depth information; obtaining DEM (Digital Elevation Model) elevation information through a method for analyzing photogrammetric survey; determining and supplementing the subsidence boundary by using elevation variation information data of the area to be surveyed, wherein the elevation variation information data of the area to be surveyed are obtained through difference analysis of DEMs at different time phases; and then eliminating the subsidence land boundary information caused by coal mining activities; and finally obtaining subsidence land boundary and subsidence depth information of the area to be surveyed. According to the invention, the boundary of the subsided farmland can be rapidly determined and basic support can be provided to compensation and treatment of the subsidence land.

Description

A kind of method of determining subsiding infringement border, arable land
Technical field
The invention belongs to mapping and environmental monitoring technology field, particularly determine the technical method on the infringement border, arable land that subsides.
Background technology
Coal is the Chinese topmost energy, and along with developing rapidly of China's mining cause, it develops when the national economic development is played a significant role, and the problem of thereupon bringing also displays.A large amount of extraction of underground coal resource cause the balance of goaf superincumbent stratum to be destroyed, produce distortion, cause surface subsidence damage in various degree, the form of surface subsidence damage has: subsidence basin, collapse doline, step, surface cracks etc., in deifferent regions.China, its degree of subsiding is also different, such as in Eastern Plain district being the Chinese compound main producing region of coal-grain, the destruction enormous amount that underground exploitation is ploughed to ground, there is very serious surface collapse in many mining areas, and cause serious mining area ponding problem, the proportion that wherein infringement is ploughed is larger, Cultivated Land Area Decrease, quality degradation, badly influence the agricultural production of this area, the people's daily life.
According to existing mining subsidence and three unders in coal mining (refer under buildings, under railway and coal mining under water-bodies method) standard, subsidence basin refers to the region that submergence depth is greater than 10mm.And for ploughing, subsiding is greater than 10mm and agricultural production is not affected, Wei Fei disturbed area, this region, the region that crop yield is exerted an influence that only subsides is only the infringement region of ploughing.
The boundary demarcation that mining area destruction is ploughed at present lacks unified standard, cause in ground dispute is subsided in processing and cannot follow on the one hand, cause workers and peasants' aggravation of contradictions, also affected on the other hand the enforcement of land reclamation work, even coal industry is produced and also caused immeasurable loss.
Summary of the invention
The object of the invention is for solving an above-mentioned difficult problem, a kind of method of determining subsiding infringement border, arable land is proposed, be intended to make up subside and damage the deficiency that arable land boundary demarcation lacks unified standard, this method combines subside buffer model and DEM differential technique based on remote sensing ponding edge, can determine fast the border that subsides and plough, for compensation and the improvement on the ground that subsides provide shoring of foundation.
A kind of technical method of determining subsiding infringement border, arable land that the present invention proposes, is characterized in that, comprises the following steps:
1) choosing of data source and determining of the infringement regional extent of subsiding to be measured: the start time point occurring according to district to be measured infringement and carry out the time point of determining work in border and set section detection time obtains remotely-sensed data (comprise satellite, aerial remote sensing images and satellite, AEROSAT image cubic phase to) and the topographic map data (non-remotely-sensed data) of location, mining area phase when more than at least 2 different in this time period; Distribute and coal seam information and mining information (as mining engineering plan) according to the mining area coal obtaining, carry out prediction of mining subsidence, draw depression isogram, determine the preliminary scope in the infringement region that subsides to be measured;
2) subsided water area boundary information extracts: described remotely-sensed data is carried out to geometry correction, image co-registration, image enhancement processing, colored synthetic, image co-registration processing, and according to step 1) in determine preliminary range boundary, from remote sensing images, delete the image-region outside this preliminary range boundary; Utilize the method for multi-scale wavelet transformation rim detection to detect the border, waters in remote sensing images, and combination Typical Land survey information on the spot, determine and supplement extracting subsided water area boundary information; Described terrain data is carried out to geometry correction processing, and according to step 1) in the infringement preliminary scope in region of subsiding separation and Extraction in topographic map data of determining go out subsided water area boundary information;
3) the subside foundation of ground border buffer model and obtaining of surface collapse depth information: (mining area is all provided with subsiding observation station first to obtain the Ground Subsidence Monitoring data in district to be measured, regularly the observation station in research station is observed and obtained data with total powerstation), and according to settlement monitoring data, set up the exponential model of border, the ground buffering of subsiding by numerical simulation:
d i = d 0 ( 1 - ln F i ln f i )
Wherein, F irepresent the degree of subsiding that distance is subsided in central point (the region geometric center point of subsiding) different distance, f 0the size of synthesis index of representation space entity generally must participate in computing, d after maximal value standardization 0represent maximum effect distance that this subsides, d ibe illustrated in the actual range at this Mou Dianjugai center in maximum effect distance of center that subsides; According to described border buffer model, to step 2) in the subsided water area boundary information that extracts carry out border buffered, thereby obtain preliminary scope subside boundary information and surface collapse depth information;
4) foundation and the differential analysis of subsidence area ground elevation model DEM to be measured: if analytic application photogrammetric survey method obtains preliminary scope to satellite, AEROSAT image cubic phase ground elevation model DEM(to carrying out Treatment Analysis is non-remotely-sensed data, need from topomap, isolate contour map layer, then in Geographic Information System program program, generate corresponding DEM); While utilizing obtain more than at least 2 different, the DEM of the phase difference processing (differential analysis) that superposes obtains preliminary scope elevation change information data;
5) border, subsidence area is revised: in Geographic Information System program, utilize step 4) preliminary scope elevation change information data and step 3) in preliminary scope subside boundary information and surface collapse depth information data superpose, comprehensively, the border that subsides is determined and is supplemented; Then contrast depression isogram, weed out the ground boundary information that subsides that non-coal mining activity causes, thereby realized step 3) in the object revised of result, finally obtain subsidence area range boundary and surface collapse depth information.
The present invention mainly has following technological merit:
1) adopt multi-source data to realize rim detection and the information extraction on the ponding border that subsides;
2) obtaining on the Ground Subsidence Monitoring data basis at sedimentation station, district to be measured, set up the exponential model of border, the ground buffering of subsiding by numerical simulation, this model is applied to the testing result on subsided water area border, boundary information and the surface collapse depth information on the ground that can obtain fast subsiding;
3) applied satellite and aviation image picture pair, adopt the method for analytical photogrammetry to generate DEM, and the district to be measured elevation change information data of using the DEM differential analysis of obtained at least 2 phases when different to obtain are determined and supplement the border that subsides, then contrast depression isogram, weed out the ground boundary information that subsides that non-coal mining activity causes, realize the object that border buffering result is revised, guaranteed to cushion correctness and the accuracy of the result of subsiding.
Accompanying drawing explanation
Fig. 1 is the technology of the present invention FB(flow block);
Fig. 2 is the waters boundary graph that the method by multi-scale wavelet transformation rim detection of the present embodiment obtains;
Fig. 3 be the present embodiment with after the buffer model buffered of border subside ground areal map;
The analytical Photogrammetric Evaluation process surface chart of Fig. 4 the present embodiment;
Fig. 5 is the DEM figure in district to be measured of the present embodiment;
Revised ground bounds and the surface collapse depth information figure of subsiding of Fig. 6 the present embodiment;
Specific implementation method
The present invention propose a kind of determine subsiding plough infringement border technical method, the inventive method overall procedure block diagram as shown in Figure 1, by reference to the accompanying drawings and embodiment be described in detail as follows (take Yanzhou Mining District Yanzhou coalfield as example):
1) data source choose and to be measured subside infringement regional extent determine: Yanzhou Mining District is established in 1958, also has so far the history of more than 50 years, and the cycle of subsiding is longer, subside as obtained Spatio-temporal Evolution situation, so just need to select multi-source data research.The present embodiment is set Duan Weicong nineteen ninety-five to 2007 detection time year topographic map data, at this moment the remotely-sensed data of choosing this location, mining area in section between has November 28 nineteen ninety-five, on August 21st, 2,000 2 phase LandsatTM remote sensing images (30m resolution), the aviation image figure (0.5m resolution) of 2005, the IRSP5 image (2.5m resolution) in October, 2007 and the while bar landsat data in mutually; Topographic map data (non-remotely-sensed data) comprises calendar year 2001 each mining area 1:2000 topomap, the surface-underground contrast plan in mining area, present landuse map, clay detailed Investigation data topographic map data etc.; Distribute and coal seam information and mining information (mining engineering plan) carry out prediction of mining subsidence according to the mining area coal obtaining, draw depression isogram, determine that to be measured subsiding damage the preliminary scope in region;
2) subsided water area boundary information extracts: the remotely-sensed data that step 1) is obtained is carried out geometry correction, image co-registration, image enhancement processing, colored synthetic, image co-registration processing, and according to step 1) in the preliminary range boundary determined from remote sensing images, delete the image-region outside this preliminary range boundary; Utilize the method for multi-scale wavelet transformation rim detection to detect the border, waters in remote sensing images, and combination Typical Land survey information on the spot, determine and supplement and (be illustrated in figure 2 waters Boundary Detection result images extracting subsided water area boundary information, dotted line is to subside to damage the preliminary scope in region, and solid line is subsided water area border); Non-remotely-sensed data is carried out to geometry correction processing, and according to step 1) in the infringement preliminary scope in region of subsiding separation and Extraction in topographic map data of determining go out subsided water area boundary information;
3) the subside foundation of ground border buffer model and obtaining of surface collapse depth information: (mining area is all provided with subsiding observation station first to obtain the Ground Subsidence Monitoring data of preliminary scope, regularly the observation station in research station is observed and obtained data with total powerstation), and according to settlement monitoring data, set up the exponential model of border, the ground buffering of subsiding by numerical simulation:
d i = d 0 ( 1 - ln F i ln f i )
Wherein, F irepresent the degree of subsiding that distance is subsided in central point (the region geometric center point of subsiding) different distance, f 0the size of synthesis index of representation space entity generally must participate in computing, d after maximal value standardization 0represent maximum effect distance that this subsides, d ibe illustrated in the actual range at this Mou Dianjugai center in maximum effect distance of center that subsides.According to above-mentioned border buffer model, to step 2) in extract subsided water area boundary information carry out border buffered, thereby obtain preliminary range boundary information and surface collapse depth information and (be illustrated in figure 3 the ground areal map that subsides after the buffer model buffered of border, grey segment is the buffering results area of subsiding with subsiding, dotted line frame is to subside to damage the preliminary scope in region, and base map is remote sensing images figure);
4) foundation and differential analysis of subsidence area ground elevation model DEM to be measured: analytic application photogrammetric survey method (if Fig. 4 is image processing program surface chart) if to satellite, AEROSAT image cubic phase to carry out Treatment Analysis obtain 4, the ground of preliminary scope when different the elevation model DEM(of phase be non-remotely-sensed data, need from topomap, isolate contour map layer, then in Geographic Information System program program, generate corresponding DEM).Figure 5 shows that the preliminary scope DEM figure in district to be measured, the depth of its color represents the variation of elevation size.Utilize the DEM of obtain 4 phases when the different difference processing (differential analysis) that superposes to obtain subsidence area to be measured preliminary scope elevation change information data;
5) border, subsidence area is revised: in Geographic Information System program, utilize step 4) district to be measured preliminary scope elevation change information data and step 3) in the preliminary scope in district to be measured subside boundary information and surface collapse depth information data superpose, comprehensively, the border that subsides is determined and is supplemented; Then contrast depression isogram, weed out the ground boundary information that subsides that non-coal mining activity causes, thereby realized step 3) in the object revised of result, geographic coverage border and surface collapse depth information (be illustrated in figure 6 revised ground bounds and the surface collapse depth information figure of subsiding, different colors represents the different depth rangees that subsides) finally obtain subsiding.

Claims (1)

1. a method of determining the infringement border, arable land that subsides, is characterized in that, comprises the following steps:
1) choosing of data source and determining of the infringement regional extent of subsiding to be measured: damage the start time point occurring and carry out the time point of determining work in border and set section detection time according to district to be measured, the remotely-sensed data of phase obtain more than at least 2 difference in location, mining area in this time period time, comprise satellite, aerial remote sensing images and satellite, AEROSAT image cubic phase pair, and topographic map data; Distribute and coal seam information and mining information according to the mining area coal obtaining, carry out prediction of mining subsidence, draw depression isogram, determine the preliminary scope in the infringement region that subsides to be measured;
2) subsided water area boundary information extracts: described remotely-sensed data is carried out to geometry correction, image co-registration, image enhancement processing, colored synthetic, image co-registration processing, and according to step 1) in determine preliminary range boundary, from remote sensing images, delete the image-region outside this preliminary range boundary; Utilize the method for multi-scale wavelet transformation rim detection to detect the border, waters in remote sensing images, and combination Typical Land survey information on the spot, determine and supplement extracting subsided water area boundary information; Described terrain data is carried out to geometry correction processing, and according to step 1) in the infringement preliminary scope in region of subsiding separation and Extraction in topographic map data of determining go out subsided water area boundary information;
3) the subside foundation of ground border buffer model and obtaining of surface collapse depth information: first obtain the Ground Subsidence Monitoring data in district to be measured, and according to settlement monitoring data, by numerical simulation sets up subside the exponential model that cushions of border:
d i = d 0 ( 1 - ln F i ln f i )
Wherein, F irepresent the degree of subsiding that distance is subsided in central point different distance, f ithe size of synthesis index of representation space entity participates in computing, d after maximal value standardization 0represent maximum effect distance that this subsides, d ibe illustrated in the actual range at this Mou Dianjugai center in maximum effect distance of center that subsides; According to described border buffer model, to step 2) in the subsided water area boundary information that extracts carry out border buffered, thereby obtain preliminary scope subside boundary information and surface collapse depth information;
4) foundation and the differential analysis of subsidence area ground elevation model DEM to be measured: analytic application photogrammetric survey method obtains the ground elevation model DEM of preliminary scope to carrying out Treatment Analysis to satellite, AEROSAT image cubic phase, to topographic map data, first from topomap, isolate contour map layer, then in Geographic Information System program program, generate corresponding DEM; While utilizing obtain more than at least 2 different, the DEM of the phase difference processing that superposes obtains preliminary scope elevation change information data;
5) border, subsidence area is revised: in Geographic Information System program, utilize step 4) preliminary scope elevation change information data and step 3) in preliminary scope subside boundary information and surface collapse depth information data superpose, comprehensively, the border that subsides is determined and is supplemented; Then contrast depression isogram, weed out the ground boundary information that subsides that non-coal mining activity causes, thereby realized step 3) in the object revised of result, finally obtain subsidence area range boundary and surface collapse depth information.
CN201210552842.9A 2012-12-18 2012-12-18 Method for determining damage boundary of subsided farmland Active CN103047970B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210552842.9A CN103047970B (en) 2012-12-18 2012-12-18 Method for determining damage boundary of subsided farmland

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210552842.9A CN103047970B (en) 2012-12-18 2012-12-18 Method for determining damage boundary of subsided farmland

Publications (2)

Publication Number Publication Date
CN103047970A CN103047970A (en) 2013-04-17
CN103047970B true CN103047970B (en) 2014-05-21

Family

ID=48060702

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210552842.9A Active CN103047970B (en) 2012-12-18 2012-12-18 Method for determining damage boundary of subsided farmland

Country Status (1)

Country Link
CN (1) CN103047970B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107449404A (en) * 2017-09-12 2017-12-08 中煤航测遥感集团有限公司 DLG collecting methods and device

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103472037B (en) * 2013-09-10 2015-12-09 淮南矿业(集团)有限责任公司 The monitoring water environment method of depression pools zone and device
CN103940407B (en) * 2014-02-13 2016-09-07 鲁东大学 One extracts coombe method based on landform and Remote Sensing Image Fusion technology
CN106097130A (en) * 2016-06-22 2016-11-09 中国矿业大学(北京) Coal-mining subsidence based on depression parameter is ploughed and is damaged the determination method on border
CN108469255B (en) * 2018-02-25 2021-03-05 长江空间信息技术工程有限公司(武汉) River and lake historical water level calculation method
CN109558859B (en) * 2018-12-13 2020-11-24 中国地质大学(武汉) Mining area distribution information extraction method and system based on DInSAR and DCNN
CN111101941B (en) * 2020-01-07 2021-04-30 安徽省交通航务工程有限公司 Boundary optimization-based shallow advanced treatment method for excavation depth pad of high-diving-level coal mining subsidence area
CN115063546B (en) * 2022-06-08 2024-04-02 中南大学 Method for delineating range of open-air-to-underground security pillar
CN115713531B (en) * 2023-01-05 2023-05-09 山东环宇地理信息工程有限公司 InSAR-based earth surface image data processing system
CN118071843A (en) * 2023-11-20 2024-05-24 南京佳格耕耘科技有限公司 Farmland mapping system based on remote sensing data

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007188117A (en) * 2006-01-11 2007-07-26 Asahi Koyo Kk Cave-in area extraction method, device and program
CN101334278A (en) * 2008-08-05 2008-12-31 中国水电顾问集团华东勘测设计研究院 Digital remote sense geological mapping process and device
CN102494667A (en) * 2011-11-30 2012-06-13 同济大学 Characterizing method of land subsidence
CN102609782A (en) * 2012-01-10 2012-07-25 中国矿业大学(北京) Optimization method of simultaneous harvest and reclamation opportunity based on mining subsidence simulation and second ploughing rate
CN102680972A (en) * 2012-06-04 2012-09-19 中国神华能源股份有限公司 Method and device for monitoring surface deformation and data processing equipment
CN102706328A (en) * 2012-04-28 2012-10-03 中国神华能源股份有限公司 Method and device for detecting subsidence area of ground and data processing equipment
CN102758622A (en) * 2012-04-28 2012-10-31 中国神华能源股份有限公司 Monitoring method for earth surface environment damage in underground mining area

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007188117A (en) * 2006-01-11 2007-07-26 Asahi Koyo Kk Cave-in area extraction method, device and program
CN101334278A (en) * 2008-08-05 2008-12-31 中国水电顾问集团华东勘测设计研究院 Digital remote sense geological mapping process and device
CN102494667A (en) * 2011-11-30 2012-06-13 同济大学 Characterizing method of land subsidence
CN102609782A (en) * 2012-01-10 2012-07-25 中国矿业大学(北京) Optimization method of simultaneous harvest and reclamation opportunity based on mining subsidence simulation and second ploughing rate
CN102706328A (en) * 2012-04-28 2012-10-03 中国神华能源股份有限公司 Method and device for detecting subsidence area of ground and data processing equipment
CN102758622A (en) * 2012-04-28 2012-10-31 中国神华能源股份有限公司 Monitoring method for earth surface environment damage in underground mining area
CN102680972A (en) * 2012-06-04 2012-09-19 中国神华能源股份有限公司 Method and device for monitoring surface deformation and data processing equipment

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
中巴资源卫星影像HR数据在煤矿矿区地面塌陷调查中的应用;陈文平,等;《测绘与空间地理信息》;20120228;第35卷(第2期) *
基于小波变换的GPS快速精密定位;郭秋英,等;《煤炭学报》;20071130;第32卷(第11期) *
王钦军,等.矿山地面塌陷的高分辨率遥感识别与边界提取.《国土资源遥感》.2011,(第3期,总第90期), *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107449404A (en) * 2017-09-12 2017-12-08 中煤航测遥感集团有限公司 DLG collecting methods and device
CN107449404B (en) * 2017-09-12 2019-12-06 中煤航测遥感集团有限公司 DLG data acquisition method and device

Also Published As

Publication number Publication date
CN103047970A (en) 2013-04-17

Similar Documents

Publication Publication Date Title
CN103047970B (en) Method for determining damage boundary of subsided farmland
Han et al. Extraction and analysis of geological lineaments combining a DEM and remote sensing images from the northern Baoji loess area
Zhang et al. Investigating slow-moving landslides in the Zhouqu region of China using InSAR time series
CN103049655B (en) Mining area land ecological damage information acquisition method based on multiple data sources
DeLong et al. Multitemporal ALSM change detection, sediment delivery, and process mapping at an active earthflow
CN101221246B (en) Remote sensing and quantizing reconnaissance method of snowslide
CN102426395B (en) Slip mass volume parameter measuring method based on three-dimensional geologic remote sensing interpretation technology
CN102938095A (en) Mining subsidence monitoring and early warning method based on multi-source data
CN105444730A (en) Time-space characteristic and cross-border mining identification method for multi-source data monitoring mining area deformation
Ghadiry et al. A new GIS-based model for automated extraction of Sand Dune encroachment case study: Dakhla Oases, western desert of Egypt
CN104123470A (en) Method for optimizing land subsidence monitoring net
Deffontaines et al. Active interseismic shallow deformation of the Pingting terraces (Longitudinal Valley–Eastern Taiwan) from UAV high-resolution topographic data combined with InSAR time series
CN111932591B (en) Method and system for remote sensing intelligent extraction of typical geological disasters
CN103700063B (en) Topography integration quick mapping method based on high definition satellite image
Infante et al. Differential SAR interferometry technique for control of linear infrastructures affected by ground instability phenomena
Pandey et al. Analysing topographical changes in open cast coal-mining region of Patratu, Jharkhand using CARTOSAT-I Stereopair satellite images
CN101788664A (en) Land degradation response unit-oriented high-spectrum land degradation information extracting method
CN117475314B (en) Geological disaster hidden danger three-dimensional identification method, system and medium
Pigazzi et al. Geology of the area of the Piuro 1618 event (Val Bregaglia, Italian Central Alps): the setting of a catastrophic historical landslide
CN114674277A (en) Method for monitoring surface subsidence of deep mining and evaluating land reclamation suitability by combining whole-field survey lines
Cigna et al. Exploiting satellite SAR for archaeological prospection and heritage site protection
Palamara et al. An evaluation of airborne laser scan data for coalmine subsidence mapping
Hsu et al. Seasonal, long-term, and short-term deformation in the Central Range of Taiwan induced by landslides
CN116863029A (en) Earthquake disaster emergency drawing method and device based on multi-source remote sensing data
CN112166688B (en) Method for monitoring desert and desertification land based on minisatellite

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant