CN105242279A - Landslide body change detection method based on laser radar technology - Google Patents
Landslide body change detection method based on laser radar technology Download PDFInfo
- Publication number
- CN105242279A CN105242279A CN201510454755.3A CN201510454755A CN105242279A CN 105242279 A CN105242279 A CN 105242279A CN 201510454755 A CN201510454755 A CN 201510454755A CN 105242279 A CN105242279 A CN 105242279A
- Authority
- CN
- China
- Prior art keywords
- sliding mass
- model
- slope
- laser radar
- change
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
- Alarm Systems (AREA)
Abstract
The present invention discloses a landslide body change detection method based on a laser radar. According to the method, in a landslide disaster, especially for the change detection of a landslide body at a complex terrain, the change of position and attitude of the landslide body at a certain time period is detected automatically, once the change of a potential landslide body is detected and a change range is outlined, the field staff is reminded to pay attention and take safety measures or leave a danger zone. The method has the advantages of low cost, easy operation, fast computing, high accuracy and convenient detection, a scientific basis is provided for landslide disaster monitoring and evaluation, the security is provided for natural disaster rescue, and economic and social benefits are produced.
Description
Technical field
The present invention relates to landslide monitoring field, a kind of sliding mass change detecting method based on laser radar technique in room.
Background technology
Along with the acceleration of socio-economic development, the aggravation of economic activity, to the excessive exploitation of natural resources, the sharp increase of the size of population, the deterioration of ecologic environment, make China's disaster present the development trend be on the rise, the work of fighting disasters and ensuring adequate disaster relief faces a severe challenge.
Landslide disaster is the sliding phenomenon that slope Rock And Soil occurs along used logical shear breakage, and its mechanism is that on a certain glide plane, shear stress has exceeded caused by the shearing strength in this face.In Southwest China, particularly southwestern Hills, the most basic topography and landform character is exactly that massif is numerous, mountain shape is steep, cheuch river is dispersed throughout among massif, mutually cuts with it, thus forms numerous Inclination body with enough sliding spaces and cut surface, the extensive pacing items that there is landslide and occur, landslide disaster is quite frequent.In addition, earthquake is very large on the impact on landslide, is first that the strong effect of earthquake makes the inner structure of slope soil stone occur to destroy and change, original structural plane open split, lax, it is very disadvantageous for adding that underground water also has the unexpected rising of larger change, particularly underground water table or reduces slope stability.In addition, the generation of a violent earthquake is often along with many aftershocks, and under the vibratory impulse repeatedly of seismic force, slope soil stone body just more easily deforms, and finally will develop into landslide.Such as, 5.12 Wenchuan violent earthquakes cause huge life and property loss to the country and people, to the landslide disaster that local geologic structure causes considerable damage to cause, still threaten the productive life of local people so far, and this threat also will continue within the following long duration.In addition, the landslide of causing due to tectonic structure, rainfall also happens occasionally throughout the country, the Chongqing WuLong, Yunnan mountain of papers, prestige, landslide, Fengqing etc. that within such as 2009, occur.The feature of landslide disaster occurs suddenly, and destructive power is huge, and particularly in the rescue operations of landslide, tumbledown Rolling Stone, rock slope, upper side slope and crag moment constitute a threat to masses suffering from disaster and rescue personnel.These landslides produce harm to buildings, farmland, road and underground utilities etc. or destroy, and cause a series of ecological environment problem.
For landslide, done a lot of research both at home and abroad, the method on traditional detection landslide has: the earth precision measurement, GPS layout of the monitoring points, InSAR technology etc., but there is certain shortcoming while having some superiority.Such as, traditional geodetic surveying mode and GPS layout of the monitoring points are point type detection means, are difficult to the information obtaining space.And InSAR technology is as the technology of high precision monitor Ground Deformation, easily produces phase misalignment when surface relief is larger and be correlated with, thus affect computational accuracy.
As the laser radar technique of one of Spatial Information Technology, be not subject to the impact of weather and surrounding physical environment, and image data precision and efficiency aspect are all obviously better than alternate manner, at present in disaster detection of dynamic such as landslides, and forecast aspect embodies very important effect and advantage.
Summary of the invention
The object of this invention is to provide a kind of sliding mass change detecting method based on laser radar, the method is in earthquake and landslide disaster, especially to complicated landform dive sliding mass carry out change detect, automatically the position of certain time period sliding mass and the change of attitude generation is detected, once detect the change that potential sliding mass occurs, and draw a circle to approve variation range and the contingent change of next step sliding mass of prediction, remind field personnel to note, take safety measures or leave hazardous location.
For achieving the above object, technical solution of the present invention is:
Based on a sliding mass change detecting method for laser radar, the change that potential sliding mass carries out is detected; It is characterized in that, comprise concrete steps as follows:
A) two phases or the two phases above laser radar point cloud data of sliding mass study area different time is chosen;
B) first to steps A) in choose study area laser radar point cloud data carry out pre-service, the complexity according to landform carries out resampling to laser radar point cloud data, improves high efficiency and accuracy that the later stage calculates;
C) to step B) cloud data after resampling carries out coordinate conversion, the cloud data coordinate system that laser radar obtains is the world geodetic system [XYZ] under WGS84, sliding mass cloud data after resampling is converted to local local coordinate system [xyz], so that unified;
D) after cinematic data being transformed into local coordinate system, iterative linear least square interpolation filter is carried out to it, ground point and non-ground points is sorted out after filtering, after ground point is and rejects vegetation, buildings, real terrain surface point, filtered ground point can simulate sliding mass model more truly;
E) step D) ground point that obtains after filtering, to the DEM of ground point cloud data construct sliding mass study area;
F) two phases dem datas of sliding mass study area are calculated obtain gradient model, the slope aspect model of its each issue DEM respectively, from whole potential landslide areas according to the dangerous landslide areas of Analyzing on Size of the gradient and slope aspect;
G) to step e) the sliding mass study area DEM that obtains after process, formed if this sliding mass is physiographic relief, DEM is utilized to obtain the level line of study area, if this sliding mass is formed after manual intervention (side slope of such as surface mine, the side slope of dam), characteristic curve at the bottom of the Po Ding slope of sliding mass, the characteristic curve of step and level line can be obtained according to the feature of landform;
H) to step e), F) and G) in DEM, gradient model, slope aspect model, the contour line model of sliding mass study area, characteristic curve model at the bottom of Po Ding slope carries out doing difference and calculates, and obtains the characteristic curve difference variation model of the difference in height variation model of sliding mass in certain period, algebraic difference between adjacent gradients variation model, slope aspect difference variation model, level line difference variation model, characteristic curve difference variation model at the bottom of Po Ding slope, step;
I) to step H) in difference in height variation model estimation sliding mass landslide rock earth volume, for the later stage landslide cleaning transport correct rock earth volume data are provided;
J) according to step H) in the various variation models of sliding mass that obtain, change is carried out to the change of sliding mass and detects.
Described sliding mass change detecting method, is characterized in that, described potential sliding mass, is rock slope, upper side slope, crag and Rolling Stone.
The described sliding mass change detecting method based on laser radar, it is characterized in that, described A) step laser radar point cloud data can be any one of the radar data on airborne radar data and ground, quantitatively be at least two phase scan-datas, present patent application have employed the radar data of two phase territorial laser scannings.
Described sliding mass change detecting method, it is characterized in that, described B) step, the data volume of Laser Radar Scanning is huge, pre-service resampling determines sliding mass cloud data dot density according to the complexity of actual landform, dot density is too high, can affect the speed of process, dot density is too low, well can not describe terrain model, therefore need to select suitable dot density, the resampling of cloud data is carried out by the method for curvature analysis, should data volume unsuitable excessive, also farthest to retain the unique point of landform.
Described sliding mass change detecting method, it is characterized in that, described C) point cloud model of sliding mass study area is local coordinate system from WGS84 coordinate conversion in step, it no matter is spaceborne, airborne and ground laser radar, its coordinate system is all WGS84 coordinate system, and being converted to local coordinate system can be fused in unified coordinate system for other system.
Described sliding mass change detecting method, it is characterized in that, described D) and E) in step, the earth's surface of sliding mass can cover vegetation or other atural object, by obtaining the model that sliding mass presses close to earth's surface most after iterative linear least square interpolation filter, the height value of culture point is higher than height value topocentric around it, after interpolation, the regression criterion in laser footpoint elevation and matching face disobeys normal distribution, the regression criterion of culture point be all on the occasion of and deviation is larger, topocentric regression criterion is less and may for negative value. and the method is iterative computation, first apply thicker graticule mesh interpolation and generate initial fit surface, in fact this surface is a face between true ground (DTM) and atural object surface (DSM), carry out interpolation process again, go down to obtain final DEM with this iteration.
Described sliding mass change detecting method, it is characterized in that, described E) DEM that generates in step is high-resolution DEM, at present, the precision of the DEM obtained by laser radar is relatively higher with other modes, F) and G) extract in step gradient model, slope aspect model, contour line model, characteristic curve model at the bottom of Po Ding slope, step characteristic curve model be all high-precision sliding mass model.
The high precision sliding mass aspect of model comprises:
A) sliding mass gradient model can be potential landslide areas from sliding mass study area whole observation to which, according to the size of the gradient, sliding mass survey region is classified: the gradient is less than 30 degree of regions belonging to more stable, the gradient belongs to the region that landslide easily occurs between 30 degree and 50 degree, and the gradient is greater than 50 degree and belongs to the region very easily occurring to come down;
B) the landslide trend of sliding mass after coming down can be observed according to sliding mass slope aspect model, and the region that landslide ground is easily piled up;
C) generation of characteristic curve model, step-feature line model at the bottom of the contour line model of sliding mass, Po Ding slope is mainly for the sliding mass formed under different mechanism, there is self-assembling formation, have manually to fill out and dig formation, take wherein one or several sliding mass analytical models.
Sliding mass change detecting method feature comprises:
A) sliding mass gradient difference variation model does difference by two phase sliding mass gradient models to obtain, the slope change of interior sliding mass study area during this period of time, can analyze the change of sliding mass in this period according to sliding mass gradient difference variation model, which block region there occurs landslide;
B) sliding mass slope aspect difference variation model does difference by two phase sliding mass slope aspect models to obtain, the slope aspect change of interior sliding mass study area during this period of time, the change trend of sliding mass in this period can be analyzed according to sliding mass slope aspect difference variation model, and next step Landslide deposit district;
C) the difference in height variation model of sliding mass does difference by two phase DEM to obtain, and analyze from the H angle of block mold, which block region height increases, and which block region height reduces;
D) the difference variation model of characteristic curve at the bottom of the level line difference variation model of sliding mass, Po Ding slope, the difference variation model of step-feature line are the contour line model of sliding mass model two phase respectively, characteristic curve model, step-feature line model at the bottom of Po Ding slope compare, the change of sliding mass is observed from trickle linear angles, can carry out cutting that is arbitrarily angled and position to model where necessary, the change of sliding mass is observed in segmentation.
Described sliding mass change detecting method, is characterized in that, I) calculating of rock earth volume in step involves the landslip treatmant in later stage, and transport talcum in time, rational scientific arrangement is carried out to the progress of whole duration.
Described sliding mass change detecting method, it is characterized in that, H), G), H) and multiple sliding mass parameter model I) in step comprehensively obtain the method that sliding mass change detects, the change detecting method from entirety to local has been carried out in the change proposed for sliding mass.
Accompanying drawing explanation
Fig. 1 is a kind of sliding mass change detecting method process flow diagram based on laser radar of the present invention;
Fig. 2 is the two phases satellite image superposition panorama sketch (blue travelling belt is the position after sliding, and grey travelling belt is the position before sliding, and is landslide areas in red wire frame) in sliding mass Disaster Study district of the present invention;
Fig. 3 is landslide disaster rescue site;
Fig. 4 a is the first phase in 2013 sliding mass cloud data after pre-service resampling;
Fig. 4 b is the second phase in 2014 localized landslip body cloud data after pre-service resampling;
Fig. 5 is design sketch before and after sliding mass part filter;
Fig. 6 is the DEM obtained after sliding mass filtering;
Fig. 7 a is the contour map that dem data extracts;
Fig. 7 b is the feature line chart at the bottom of Po Ding slope extracted;
Fig. 7 c is the feature line chart of the step extracted;
Fig. 8 a is the gradient illustraton of model (gradient size 0-87.6 degree) of sliding mass laser point cloud in 2014;
Fig. 8 b is the slope aspect illustraton of model (slope aspect size 0-180 degree) of sliding mass laser point cloud in 2014;
Fig. 9 a is according to the location diagram (light grey model is laser point cloud model in 2013, and Dark grey model is laser point cloud model in 2014) that is stacked together with two issues of coordinate grid;
Fig. 9 b is difference in height variation model figure (upper figure is overall condition, and figure below is southern nation situation);
Fig. 9 c is algebraic difference between adjacent gradients variation model figure (slope change size-65.9-87.3 degree);
Fig. 9 d is level line variation model figure;
Fig. 9 e is step-feature line variation model figure.
Embodiment
The present invention is a kind of sliding mass change detecting method based on laser radar, and the enforcement of this research method is formed primarily of the laser scanning data of the different phase of the same area being no less than for two phases and computer system.On the whole, the power consumption of this change detecting method is very low, computational accuracy is high, counting yield is high, has following several large feature:
Feature 1: detection method of the present invention needs supporting two nucleus modules, i.e. laser radar point cloud data processing module and change detection module.Laser radar point cloud data processing module can realize:
A) because laser radar point cloud data obtains the data of the bulk redundancy comprising study area and non-study district usually, need to carry out pre-service and resampling to data, reject the redundant points beyond study area, and the dot density of general laser radar point cloud data is higher, actual conditions according to landform carry out resampling to data according to curvature analysis, have both farthest remained terrain feature point, take into account the accuracy requirement of data, have compressed data volume again, improve the efficiency that later stage change detects;
B) conversion of study area cloud data local coordinate system, laser radar device self is with the GPS locating module at WGS84 coordinate system, the data of all collections are WGS84 coordinate, for realizing and the unification of local coordinate system, boolean's sand seven-parameter transformation model is adopted to carry out space coordinate conversion to cloud data;
C) putting cloud filtering is to realize being separated of ground point and non-ground points, by iterative linear least square interpolation filter, this filtering method is the improvement based on Linear Iterative Method, relative to other filtering methods, this filtering method not only precision is obviously better than other modes, and the efficiency of filtering is very high, and the data volume of some cloud is with millions of, operand is quite large, can obtain the truest ground DEM after filtering.
Change detection module is made up of following functions:
A) obtaining characteristic curve model, step-feature line model at the bottom of the gradient model of sliding mass study area, slope aspect model, contour line model, Po Ding slope by processing later DEM, obtaining some basic models of analysis of landslide body;
B) by A more than two phases or two phases) in gradient model, slope aspect model, contour line model, characteristic curve model, step-feature line model at the bottom of Po Ding slope obtain characteristic curve variation model, step-feature line variation model at the bottom of algebraic difference between adjacent gradients variation model, slope aspect difference variation model, equation of equal altitude variation model, Po Ding slope;
C) the difference in height variation model of sliding mass study area, obtained by two phases or DEM of many phases mathematic interpolation, can obtain the rock earth volume of overall landslide areas or part landslide areas according to difference in height variation model, the landslip treatmant for the later stage provides science reference data accurately.
Feature 2: present invention employs " method based on discrete curvature is analyzed " to characteristic curve at the bottom of the resampling of original laser radar cloud data and Po Ding slope, the mentioning of step-feature line.Laser radar point cloud data be all discrete, do not have well-regulated point, in the place that a Curvature varying is less, more less points are retained during resampling, and large at Curvature varying, the place that topography variation is large, will retain more some during resampling, scientific and technological high-precision analog topographical surface model can be ensured, both have compressed the size of data volume, and remained again the feature of landform to the full extent, data will improve the efficiency that sliding mass change detects upon compression.Must be the place that terrain feature is undergone mutation in the place that these discrete point curvatures are undergone mutation, connect with line the unique point that these curvature undergo mutation and just generate characteristic curve at the bottom of Po Ding slope, step-feature line.
Feature 3: present invention employs " method of iterative linear least square interpolation filter " and carry out filtering, this filtering is the improvement carried out on the basis of linear iteration, the successive ignition matching iteration face that the method is not only considered, first apply thicker graticule mesh interpolation and generate initial fit surface, in fact this surface is a face between true ground (DTM) and atural object surface (DSM), carry out interpolation process again, go down to obtain final DEM with this iteration, and treatment effeciency is very high, no matter large again data volume, substantially the processing time of filtering be all with point/second is for chronomere.
Feature 4: core technology of the present invention is the sliding mass change detecting method of " based on sliding mass running parameter collection ", changes single in the past or several parameter, or the simple change relying on the movement of indivedual a certain amount of calibration point to carry out sliding mass detects.The sliding mass change detecting method of " based on sliding mass running parameter collection " is from tentatively obtaining gradient model parameter, slope aspect model parameter, contour line model parameter, characteristic curve model parameter at the bottom of Po Ding slope, step-feature line model parameter carries out the detection analysis of sliding mass, and obtain algebraic difference between adjacent gradients variation model parameter further, slope aspect difference variation model parameter, level line difference variation model parameter, characteristic curve variation model parameter at the bottom of Po Ding slope, step-feature line variation model parameter, difference in height variation model parameter forms sliding mass change detection parameter collection, detection of dynamic is carried out to the change of sliding mass, the change of analysis of landslide body.
Feature 5: after the most sliding mass of method that sliding mass change detects comes down, rock can change a lot on position and attitude, pile up in a certain regional extent, estimate by difference in height variation model the amount that the rock earthwork is piled up, the timely process for next step provides science data accurately.
Embodiment:
For 2013-2014 years domestic certain glory-hole landslides, utilize this research method to carry out change to this sliding mass (comprising rock slope, upper side slope, crag and Rolling Stone) and detect.
A) pre-service of laser radar point cloud data: the two phase sliding mass cloud datas obtained for territorial laser scanning carry out pre-service, some flying spots mainly in processing scan data of pre-service and assorted point, and the resampling of a cloud is carried out according to the method for curvature analysis, territorial laser scanning is while acquisition sliding mass cloud data, also other culture points beyond landslide area can be scanned, culture point beyond these sliding masses can disturb the analysis of sliding mass, the object of resampling ensures that the size of sliding mass cloud data amount is enough little, the speed calculated can be ensured like this, but the accuracy requirement of data must be ensured equally, according to the method for resampling of curvature analysis, in the mode of maximum reservation terrain feature, resampling is carried out to cloud data, pretreated two issues certificates are as Fig. 4 a, 4b,
B) coordinate conversion of cloud data: the point cloud model of sliding mass study area is local coordinate system from WGS84 coordinate conversion, ground laser radar self has positioning system, receive gps signal to position self, the data coordinate system of its scanning is all the coordinate under GWS84 coordinate system, change according to seven parameter coordinate conversion, being converted to local coordinate system can be fused in unified coordinate system for other system, and the precision after conversion is assessed, the conversion accuracy error of conversion is in ± 0.001m;
C) filtering of cloud data: the model (before and after Fig. 5 sliding mass part filter design sketch) pressing close to earth's surface by obtaining sliding mass after iterative linear least square interpolation filter most, the height value of culture point is higher than height value topocentric around it, after interpolation, the regression criterion in laser footpoint elevation and matching face disobeys normal distribution, the regression criterion of culture point be all on the occasion of and deviation is larger, topocentric regression criterion is less and may for negative value. and the method is iterative computation, first apply thicker graticule mesh interpolation and generate initial fit surface, in fact this surface is a face between true ground (DTM) and atural object surface (DSM), carry out interpolation process again, go down to obtain final DEM (Fig. 6) with this iteration,
D) the characteristic curve model of the level line of sliding mass cloud data, characteristic curve at the bottom of Po Ding slope, step: the contour map 7a extracted according to sliding mass DEM, the characteristic curve of method extraction at the bottom of Po Ding slope, the characteristic curve of step that adopt discrete curvature to analyze, there is curvature in each point of sliding mass, usual Curvature varying difference is between points all not too large, when the place that large change place must be sliding mass side slope existing characteristics line occurs curvature, according to characteristic curve Fig. 7 b, 7c of the characteristic curve at the bottom of the thinking extraction Po Ding slope of curvature mutation matching, step;
E) gradient, the slope aspect model of sliding mass: to sliding mass DEM vector to raster conversion, pixel size after rasterizing is 0.1m (can according to the size of actual conditions oneself setting pixel), then calculate the gradient of sliding mass, gradient model that slope aspect numerical value obtains sliding mass (Fig. 8 a), slope aspect model (Fig. 8 b);
F) the difference in height variation model of sliding mass, slope change model, slope aspect variation model, characteristic curve variation model at the bottom of Po Ding slope, step-feature line variation model: the corresponding sliding mass discrepancy in elevation position after the change of sliding mass landslide has corresponding change (two issues are according to the location diagram 9a be stacked together) naturally, two issues obtain difference in height variation model (Fig. 9 b) according to carrying out mathematic interpolation, the gradient, slope aspect difference variation model (Fig. 9 c) is by the two phase gradients, slope aspect model carries out mathematic interpolation, characteristic curve variation model at the bottom of Po Ding slope, step-feature line variation model carries out by the feature obtained of two phases the characteristic curve variation characteristic (Fig. 9 d) that topological relation analysis obtains opposite position,
G) come down rock Earthwork Calculation: after the change of sliding mass landslide, according to the change calculations local of the discrepancy in elevation or the talcum earth volume of entirety;
H) by characteristic curve, step-feature line at the bottom of the gradient, slope aspect, level line, Po Ding slope, and the change of discrepancy in elevation running parameter to sliding mass detects, and changes in the past that the change to sliding mass detects by single or several parameter.
Claims (10)
1., based on a sliding mass change detecting method for laser radar technique, the change that potential sliding mass carries out is detected; It is characterized in that, comprise concrete steps as follows:
A) two phases or the two phases above laser radar point cloud data of sliding mass study area different time is chosen;
B) first to steps A) in choose study area laser radar point cloud data carry out pre-service, the complexity according to landform carries out resampling to laser radar point cloud data, considers the high efficiency into the later stage calculates and accuracy;
C) to step B) cloud data after resampling carries out coordinate conversion, the cloud data coordinate system that laser radar obtains is the world geodetic system [XYZ] under WGS84, sliding mass cloud data after resampling is converted to local local coordinate system [xyz], so that unified;
D) convert after local coordinate system, iterative linear least square interpolation filter is carried out to it, after filtering, sorts out ground point and non-ground points, after ground point is and rejects vegetation, buildings, real topographical surface point, filtered ground point can simulate sliding mass model more truly;
E) step D) ground point that obtains after filtering, to the DEM of ground point cloud data construct sliding mass study area;
F) two phases dem datas of sliding mass study area are calculated obtain gradient model, the slope aspect model of its each issue DEM respectively, from whole potential landslide areas according to the dangerous landslide areas of Analyzing on Size of the gradient and slope aspect;
G) to step e) the sliding mass study area DEM that obtains after process, formed if this sliding mass is physiographic relief, DEM is utilized to obtain the level line of study area, if this sliding mass is formed after manual intervention (side slope of such as surface mine, the side slope of dam), characteristic curve at the bottom of the Po Ding slope of sliding mass, the characteristic curve of step and level line can be obtained according to the feature of landform;
H) to step e), F) and G) in DEM, gradient model, slope aspect model, the contour line model of sliding mass study area, characteristic curve model at the bottom of Po Ding slope carries out doing difference and calculates, and obtains the characteristic curve difference variation model of the difference in height variation model of sliding mass in certain period, algebraic difference between adjacent gradients variation model, slope aspect difference variation model, level line difference variation model, characteristic curve difference variation model at the bottom of Po Ding slope, step;
I) to step H) in difference in height variation model estimation sliding mass landslide rock earth volume, for the later stage landslide cleaning transport correct rock earth volume data are provided;
J) according to step H) in the various variation models of sliding mass that obtain, change is carried out to the change of sliding mass and detects.
2. a kind of sliding mass change detecting method based on laser radar technique as claimed in claim 1, it is characterized in that, described sliding mass, is rock slope, upper side slope, crag and Rolling Stone.
3. a kind of sliding mass change detecting method based on laser radar technique as claimed in claim 1, it is characterized in that, described A) step laser radar point cloud data can be any one of radar data on spaceborne radar data, airborne radar data and ground, is quantitatively at least the radar data of two phases laser scanning.
4. a kind of sliding mass change detecting method based on laser radar technique as claimed in claim 1, it is characterized in that, described B) step, the data volume of Laser Radar Scanning is huge, pre-service resampling determines sliding mass cloud data dot density according to the complexity of actual landform, dot density is too high, can affect the speed of process, dot density is too low, well can not describe terrain model, need to select suitable dot density, the resampling of cloud data is carried out by the method for curvature analysis, both ensured that data volume was unsuitable excessive, also farthest retain the unique point of landform.
5. a kind of sliding mass change detecting method based on laser radar technique as claimed in claim 3, it is characterized in that, described C) point cloud model of sliding mass study area is local coordinate system from WGS84 coordinate conversion in step, no matter airborne and ground laser radar, its coordinate system is all GWS84 coordinate system, and being converted to local coordinate system can be fused in unified coordinate system with other measuring systems.
6. a kind of sliding mass change detecting method based on laser radar technique as claimed in claim 1, it is characterized in that, described D) and E) in step, the earth's surface of sliding mass can cover vegetation or other atural object, by obtaining the model that sliding mass presses close to earth's surface most after iterative linear least square interpolation filter, the height value of culture point is higher than height value topocentric around it, after interpolation, the regression criterion in laser footpoint elevation and matching face disobeys normal distribution, the regression criterion of culture point be all on the occasion of and deviation is larger, topocentric regression criterion is less and may for negative value. and the method is iterative computation, first apply thicker graticule mesh interpolation and generate initial fit surface, in fact this surface is a face between true ground (DTM) and atural object surface (DSM), carry out interpolation process again, go down to obtain final DEM with this iteration.
7. a kind of sliding mass change detecting method based on laser radar technique as claimed in claim 1, it is characterized in that, described E) DEM that generates in step is high-resolution DEM, at present, the precision of the DEM obtained by laser radar is relatively higher with other modes, the characteristic curve model of the gradient model F) and G) extracted in step, slope aspect model, contour line model, characteristic curve model at the bottom of Po Ding slope, step is all high-precision sliding mass model, comprising:
A) sliding mass gradient model can be potential landslide areas from sliding mass study area whole observation to which, according to the size of the gradient, sliding mass survey region is classified: the gradient is less than 30 degree of regions belonging to more stable, the gradient belongs to the region that landslide easily occurs between 30 degree and 50 degree, and the gradient is greater than 50 degree and belongs to the region very easily occurring to come down;
B) the landslide trend of sliding mass after coming down can be observed according to sliding mass slope aspect model, and the region that landslide ground is easily piled up;
C) generation of characteristic curve model, step-feature line model at the bottom of the contour line model of sliding mass, Po Ding slope is mainly for the sliding mass formed under different mechanism, there is self-assembling formation, have manually to fill out and dig formation, take wherein one or several sliding mass analytical models.
8. a kind of sliding mass change detecting method based on laser radar technique as claimed in claim 1, is characterized in that, described H) and various sliding masses change detection models I) in step, comprising:
A) sliding mass gradient difference variation model does difference by two phase sliding mass gradient models to obtain, the slope change of interior sliding mass study area during this period of time, can analyze the change of sliding mass in this period according to sliding mass gradient difference variation model, which region there occurs landslide;
B) sliding mass slope aspect difference variation model does difference by two phase sliding mass slope aspect models to obtain, the slope aspect change of interior sliding mass study area during this period of time, the change trend of sliding mass in this period can be analyzed according to sliding mass slope aspect difference variation model, and next step Landslide deposit district;
C) the difference in height variation model of sliding mass does difference by two phase DEM to obtain, and analyze from the H angle of block mold, which block region height increases, and which block region height reduces;
D) the difference variation model of characteristic curve at the bottom of the level line difference variation model of sliding mass, Po Ding slope, the difference variation model of step-feature line are the contour line model of sliding mass model two phase respectively, characteristic curve model, step-feature line model at the bottom of Po Ding slope compare, the change of sliding mass is observed from trickle linear angles, can carry out cutting that is arbitrarily angled and position to model where necessary, the change of sliding mass is observed in piecewise.
9. a kind of sliding mass change detecting method based on laser radar technique as claimed in claim 1, it is characterized in that, the calculating of the rock earth volume I) in step, effectively can support that the landslide cubic meter of stone in later stage clears planning, the timely transhipment landslide cubic meter of stone, carries out rational scientific arrangement to the progress of whole duration.
10. a kind of sliding mass change detecting method based on laser radar technique as claimed in claim 1, it is characterized in that, H), G), H) and multiple sliding mass parameter model I) in step comprehensively obtain the method that sliding mass change detects, the change detecting method from entirety to local has been carried out in the change proposed for sliding mass.
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510454755.3A CN105242279B (en) | 2015-07-30 | 2015-07-30 | A kind of slip mass change detecting method based on laser radar technique |
| CN201810308179.5A CN108594261A (en) | 2015-07-30 | 2015-07-30 | A kind of slip mass change detecting method based on laser radar technique |
| CN201810308178.0A CN108594260A (en) | 2015-07-30 | 2015-07-30 | A kind of slip mass change detecting method based on laser radar technique |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510454755.3A CN105242279B (en) | 2015-07-30 | 2015-07-30 | A kind of slip mass change detecting method based on laser radar technique |
Related Child Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810308179.5A Division CN108594261A (en) | 2015-07-30 | 2015-07-30 | A kind of slip mass change detecting method based on laser radar technique |
| CN201810308178.0A Division CN108594260A (en) | 2015-07-30 | 2015-07-30 | A kind of slip mass change detecting method based on laser radar technique |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105242279A true CN105242279A (en) | 2016-01-13 |
| CN105242279B CN105242279B (en) | 2018-07-31 |
Family
ID=55039982
Family Applications (3)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810308178.0A Withdrawn CN108594260A (en) | 2015-07-30 | 2015-07-30 | A kind of slip mass change detecting method based on laser radar technique |
| CN201510454755.3A Active CN105242279B (en) | 2015-07-30 | 2015-07-30 | A kind of slip mass change detecting method based on laser radar technique |
| CN201810308179.5A Withdrawn CN108594261A (en) | 2015-07-30 | 2015-07-30 | A kind of slip mass change detecting method based on laser radar technique |
Family Applications Before (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810308178.0A Withdrawn CN108594260A (en) | 2015-07-30 | 2015-07-30 | A kind of slip mass change detecting method based on laser radar technique |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810308179.5A Withdrawn CN108594261A (en) | 2015-07-30 | 2015-07-30 | A kind of slip mass change detecting method based on laser radar technique |
Country Status (1)
| Country | Link |
|---|---|
| CN (3) | CN108594260A (en) |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106370225A (en) * | 2016-08-18 | 2017-02-01 | 中国科学院、水利部成都山地灾害与环境研究所 | Rapid surveying and imaging method for accumulated layer landslide |
| CN106777734A (en) * | 2016-12-27 | 2017-05-31 | 湖南师范大学 | A kind of mountain region Prediction of Landslide |
| CN107831484A (en) * | 2017-09-29 | 2018-03-23 | 深圳市行者机器人技术有限公司 | A kind of ladder data capture method and device based on laser radar |
| CN107860375A (en) * | 2017-11-03 | 2018-03-30 | 广州地理研究所 | A kind of landslide disaster volume rapid extracting method based on three-dimensional laser scanning technique |
| CN108090284A (en) * | 2017-12-19 | 2018-05-29 | 建基工程咨询有限公司 | Application based on laser scanning modeling reverse Engineering Technology in construction Supervision |
| CN108168526A (en) * | 2017-12-15 | 2018-06-15 | 长安大学 | A kind of Loess Landslide EARLY RECOGNITION method based on surface subsidence information |
| CN108363050A (en) * | 2018-02-09 | 2018-08-03 | 长江水利委员会长江科学院 | A kind of dyke destructive test method based on on-board LiDAR data |
| CN108876799A (en) * | 2018-06-12 | 2018-11-23 | 杭州视氪科技有限公司 | A kind of real-time step detection method based on binocular camera |
| CN109556673A (en) * | 2018-11-22 | 2019-04-02 | 中铁四局集团有限公司 | A kind of earthwork calculation amount method and system based on unmanned plane |
| CN109584515A (en) * | 2018-12-21 | 2019-04-05 | 武汉万集信息技术有限公司 | Method for early warning, device and the readable storage medium storing program for executing of massif disaster |
| CN110415331A (en) * | 2019-06-14 | 2019-11-05 | 中国电建集团昆明勘测设计研究院有限公司 | A kind of profile or hole based on point cloud data quickly detects differentiation and hole repair method |
| CN112327698A (en) * | 2020-11-05 | 2021-02-05 | 叶远 | Flood disaster early warning system and method based on Internet of things |
| CN112462361A (en) * | 2020-02-26 | 2021-03-09 | 苏州锐思突破电子科技有限公司 | Slope radar early warning algorithm |
| CN114333247A (en) * | 2022-01-07 | 2022-04-12 | 威海兴海源网具有限公司 | Disaster detection early warning system |
| CN114882366A (en) * | 2022-05-26 | 2022-08-09 | 广州市城市规划勘测设计研究院 | Three-dimensional scene catastrophe monitoring and early warning method |
| CN115471144A (en) * | 2022-11-14 | 2022-12-13 | 四川省公路规划勘察设计研究院有限公司 | Debris flow monitoring and early warning method, device and medium based on multi-source data fusion |
| CN115615344A (en) * | 2022-12-06 | 2023-01-17 | 中国地质环境监测院(自然资源部地质灾害技术指导中心) | Surface deformation monitoring method and device, electronic equipment and storage medium |
| CN117952823A (en) * | 2024-03-26 | 2024-04-30 | 青岛慧拓智能机器有限公司 | Automatic extraction, classification and updating method for slope top and bottom lines of open-pit mining area |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110412617B (en) * | 2019-09-06 | 2023-06-30 | 李娜 | Unmanned aerial vehicle rescue method based on self-feedback laser radar scanning and application |
| CN111445057B (en) * | 2020-02-27 | 2023-06-23 | 重庆大学 | Displacement gradient method for predicting landslide slip direction and moment |
| CN111983678B (en) * | 2020-07-22 | 2024-02-23 | 中海石油(中国)有限公司深圳分公司 | Method for rapidly evaluating development potential of deep water sand body |
| CN113935202B (en) * | 2021-12-17 | 2022-03-01 | 四川省交通勘察设计研究院有限公司 | River channel evolution early warning method based on laser radar |
| CN114791273B (en) * | 2022-06-24 | 2022-09-13 | 中国铁道科学研究院集团有限公司铁道建筑研究所 | InSAR deformation monitoring result interpretation method for landslide |
| CN117782226B (en) * | 2024-02-23 | 2024-05-14 | 四川省能源地质调查研究所 | Mine safety early warning system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102005104A (en) * | 2009-09-02 | 2011-04-06 | 吴立新 | Remote and rapid monitoring and alarming device and method for displacement and gesture of sliding mass |
| CN102680971A (en) * | 2012-05-22 | 2012-09-19 | 南京泰通科技有限公司 | Radar device for monitoring mountain landslide |
-
2015
- 2015-07-30 CN CN201810308178.0A patent/CN108594260A/en not_active Withdrawn
- 2015-07-30 CN CN201510454755.3A patent/CN105242279B/en active Active
- 2015-07-30 CN CN201810308179.5A patent/CN108594261A/en not_active Withdrawn
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102005104A (en) * | 2009-09-02 | 2011-04-06 | 吴立新 | Remote and rapid monitoring and alarming device and method for displacement and gesture of sliding mass |
| CN102680971A (en) * | 2012-05-22 | 2012-09-19 | 南京泰通科技有限公司 | Radar device for monitoring mountain landslide |
Non-Patent Citations (2)
| Title |
|---|
| 刘锦程: ""三维激光扫描技术在滑坡监测中的应用研究"", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 辛麒: ""基于机载激光雷达数据构建DEM的精度分析"", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106370225A (en) * | 2016-08-18 | 2017-02-01 | 中国科学院、水利部成都山地灾害与环境研究所 | Rapid surveying and imaging method for accumulated layer landslide |
| CN106777734A (en) * | 2016-12-27 | 2017-05-31 | 湖南师范大学 | A kind of mountain region Prediction of Landslide |
| CN107831484A (en) * | 2017-09-29 | 2018-03-23 | 深圳市行者机器人技术有限公司 | A kind of ladder data capture method and device based on laser radar |
| CN107860375A (en) * | 2017-11-03 | 2018-03-30 | 广州地理研究所 | A kind of landslide disaster volume rapid extracting method based on three-dimensional laser scanning technique |
| CN108168526B (en) * | 2017-12-15 | 2020-09-01 | 长安大学 | Loess landslide early-stage identification method based on ground settlement information |
| CN108168526A (en) * | 2017-12-15 | 2018-06-15 | 长安大学 | A kind of Loess Landslide EARLY RECOGNITION method based on surface subsidence information |
| CN108090284A (en) * | 2017-12-19 | 2018-05-29 | 建基工程咨询有限公司 | Application based on laser scanning modeling reverse Engineering Technology in construction Supervision |
| CN108363050A (en) * | 2018-02-09 | 2018-08-03 | 长江水利委员会长江科学院 | A kind of dyke destructive test method based on on-board LiDAR data |
| CN108363050B (en) * | 2018-02-09 | 2020-05-19 | 长江水利委员会长江科学院 | Embankment damage detection method based on airborne LiDAR data |
| CN108876799A (en) * | 2018-06-12 | 2018-11-23 | 杭州视氪科技有限公司 | A kind of real-time step detection method based on binocular camera |
| CN108876799B (en) * | 2018-06-12 | 2021-07-02 | 杭州视氪科技有限公司 | Real-time step detection method based on binocular camera |
| CN109556673A (en) * | 2018-11-22 | 2019-04-02 | 中铁四局集团有限公司 | A kind of earthwork calculation amount method and system based on unmanned plane |
| CN109584515A (en) * | 2018-12-21 | 2019-04-05 | 武汉万集信息技术有限公司 | Method for early warning, device and the readable storage medium storing program for executing of massif disaster |
| CN110415331B (en) * | 2019-06-14 | 2022-08-19 | 中国电建集团昆明勘测设计研究院有限公司 | Method for rapidly detecting, judging and repairing contour or hole based on point cloud data |
| CN110415331A (en) * | 2019-06-14 | 2019-11-05 | 中国电建集团昆明勘测设计研究院有限公司 | A kind of profile or hole based on point cloud data quickly detects differentiation and hole repair method |
| CN112462361A (en) * | 2020-02-26 | 2021-03-09 | 苏州锐思突破电子科技有限公司 | Slope radar early warning algorithm |
| CN112462361B (en) * | 2020-02-26 | 2023-11-21 | 苏州锐思突破电子科技有限公司 | Slope radar early warning algorithm |
| CN112327698A (en) * | 2020-11-05 | 2021-02-05 | 叶远 | Flood disaster early warning system and method based on Internet of things |
| CN114333247A (en) * | 2022-01-07 | 2022-04-12 | 威海兴海源网具有限公司 | Disaster detection early warning system |
| CN114882366A (en) * | 2022-05-26 | 2022-08-09 | 广州市城市规划勘测设计研究院 | Three-dimensional scene catastrophe monitoring and early warning method |
| CN115471144A (en) * | 2022-11-14 | 2022-12-13 | 四川省公路规划勘察设计研究院有限公司 | Debris flow monitoring and early warning method, device and medium based on multi-source data fusion |
| CN115471144B (en) * | 2022-11-14 | 2023-04-07 | 四川省公路规划勘察设计研究院有限公司 | Debris flow monitoring and early warning method, device and medium based on multi-source data fusion |
| CN115615344A (en) * | 2022-12-06 | 2023-01-17 | 中国地质环境监测院(自然资源部地质灾害技术指导中心) | Surface deformation monitoring method and device, electronic equipment and storage medium |
| CN117952823A (en) * | 2024-03-26 | 2024-04-30 | 青岛慧拓智能机器有限公司 | Automatic extraction, classification and updating method for slope top and bottom lines of open-pit mining area |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105242279B (en) | 2018-07-31 |
| CN108594260A (en) | 2018-09-28 |
| CN108594261A (en) | 2018-09-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105242279A (en) | Landslide body change detection method based on laser radar technology | |
| CN101221246B (en) | Remote sensing and quantizing reconnaissance method of snowslide | |
| CN113283802B (en) | Landslide risk assessment method for complicated and difficult mountain areas | |
| Barreca et al. | A pilot GIS database of active faults of Mt. Etna (Sicily): A tool for integrated hazard evaluation | |
| CN102635059B (en) | Bridge investigation method | |
| CN111582597B (en) | Method and equipment for predicting landslide hazard of power transmission line | |
| Guo et al. | Geohazard recognition and inventory mapping using airborne LiDAR data in complex mountainous areas | |
| CN102938095A (en) | Mining subsidence monitoring and early warning method based on multi-source data | |
| CN102306233A (en) | River basin landslide space-time predicting method under rainfall effect | |
| Somma et al. | Application of laser scanning for monitoring coastal cliff instability in the Pozzuoli Bay, Coroglio site, Posillipo Hill, Naples | |
| Cappadonia et al. | Preliminary results of a geomorphological and DInSAR characterization of a recently identified Deep-Seated Gravitational Slope Deformation in Sicily (Southern Italy) | |
| CN115526544A (en) | Green mine construction planning method and system based on multi-type data | |
| Zhang et al. | Spatial patterns and controlling factors of the evolution process of karst depressions in Guizhou province, China | |
| Liu et al. | Processing outcrop point clouds to 3D rock structure using open source software | |
| Narciso et al. | Developing a site-conditions map for seismic hazard assessment in Portugal | |
| Magliulo | Soil erosion susceptibility maps of the Janare Torrent Basin (southern Italy) | |
| Sajikumar et al. | Integrated remote sensing and GIS approach for groundwater exploration using analytic hierarchy process (AHP) technique | |
| Ventura et al. | Monitoring and structural significance of ground deformations at Campi Flegrei supervolcano (Italy) from the combined 2D and 3D analysis of PS-InSAR, geophysical, geological and structural data | |
| Gesch | An inventory of topographic surface changes: the value of multi-temporal elevation data for change analysis and monitoring | |
| Mala et al. | Earthquake Microzonation Study on Batubesi Dam of Nuha, East Luwu, South Sulawesi, Indonesia | |
| Stead et al. | Rock slope characterization and geomechanical modelling | |
| Ma et al. | A comprehensive approach to evaluate coastal dune evolution in Haiyang, China | |
| Kodge et al. | Elevation contour analysis and water body extraction for finding water scarcity locations using DEM | |
| She et al. | Preliminary Results and Analyses of Post-Earthquake Geological Hazards in Jiuzhaigou Based on Airborne Lidar and Imagery | |
| Huang | Multidisciplinary investigation of surface deformation above salt domes in Houston, Texas |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information |
Inventor after: Gan Xiaohong Inventor after: Xie Linger Inventor after: Zhu Jianchao Inventor after: Gu Jungang Inventor before: Wang Zhi |
|
| CB03 | Change of inventor or designer information | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20180510 Address after: 315324 Mei fang Road, Dongqian Lake holiday resort, Ningbo, Zhejiang, 12 Applicant after: Ningbo Traffic Construction Engineering Test Center Co., Ltd. Address before: 110004 No. 265 Lane 3, Lane 3, Cultural District, Heping District, Shenyang, Liaoning. Applicant before: Wang Zhi |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |