CN106093938A - A kind of mining area based on manual corner reflector side-play amount deformation monitoring method - Google Patents
A kind of mining area based on manual corner reflector side-play amount deformation monitoring method Download PDFInfo
- Publication number
- CN106093938A CN106093938A CN201610326486.7A CN201610326486A CN106093938A CN 106093938 A CN106093938 A CN 106093938A CN 201610326486 A CN201610326486 A CN 201610326486A CN 106093938 A CN106093938 A CN 106093938A
- Authority
- CN
- China
- Prior art keywords
- image
- corner reflector
- point
- play amount
- initial position
- 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
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
- G01S13/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B15/00—Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons
- G01B15/06—Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons for measuring the deformation in a solid
-
- 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
- G01S13/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention provides a kind of mining area based on manual corner reflector side-play amount deformation monitoring method, belong to Mine Monitoring field.Utilize corner reflector strong scattering feature on SAR intensity map, and impulse response characteristics, by there being limit for length unit impulse response wave filter (FIR, Finite Impulse Response) corner reflector strength signal carried out the above over-sampling of Radix Achyranthis Bidentatae, ask for the intensity peak location variation identified on time series image, as corner reflector point deformation data.This invention solves the difficult point of InSAR technical monitoring mining area big magnitude Deformation Field, and monitoring accuracy is higher compared with SAR strength offsets amount track algorithm.
Description
Technical field
The invention belongs to Mine Monitoring field, supervise particularly to the deformation of a kind of mining area based on manual corner reflector side-play amount
Survey method.
Background technology
Synthetic aperture radar interferometry (Synthetic Aperture Radar Interferometry, InSAR) this
One emerging earth observation technology, by synthetic aperture radar (Synthetic Aperture Radar, SAR) round-the-clock, whole day
Time, the advantage such as high accuracy and covering on a large scale, the nearly more than ten years are developed rapidly, are widely used to all kinds of geological disaster prison
Survey: urban ground subsidence, earthquake, volcano, landslide etc., also widely applied in terms of Coal Mining Subsidence monitoring simultaneously.
Most domestic mining causes Ground Deformation to have, and deformation magnitude is big, non-linear on the time, do not connect
Continuous, and the feature such as coverage is little, this makes all kinds of InSAR technology based on phase information face technical bottleneck, main
Be cannot to solve actual largest deformation amount and contradiction between the largest deformation gradient that InSAR can monitor.It is additionally based on SAR strong
Degree side-play amount tracking technique (Intensity offset-tracking) proposes, although solving InSAR technology can monitor deformation
The restricted problem of gradient, but it is affected by other factors such as image resolutions, and calculation accuracy is relatively low.
During realizing the present invention, inventor finds that prior art at least there is problems in that
Existing InSAR technology cannot break through its restriction that can monitor deformation gradient, accurately obtains mining area major part depression
The deformation quantity in district and maximum subsidence amount；It is additionally based on the side-play amount tracking technique of SAR intensity map by image resolution, earth's surface phase
It is relatively low that the factors such as Guan Xing affect its precision resolved.The problems referred to above can affect the accurate estimation to mining subsidence parameter further,
And then impact is to the grasp of Coal Mining Subsidence rule in specific region.
Summary of the invention
In order to solve problem of the prior art, the invention provides a kind of mining area based on manual corner reflector side-play amount shape
Change monitoring method, described mining area based on manual corner reflector side-play amount deformation monitoring method, including:
From based on the whole image of seasonal effect in time series being chosen with reference to image, determine in whole image in addition to reference image
System deviation amount between residue image and reference image；
Longitude and latitude according to known corner reflector point and height, determine described angle based on RD geo-location model
Reflector points position in whole images, is set to the initial position of described corner reflector point, based on initial bit by described position
Put, open centered by initial position in whole images and take a certain size window and carry out over-sampling process, acquisition process after strong
Degree peak；
Based on the intensity peak got, in conjunction with before the system deviation amount that determines carry out difference processing, determine described angle
The pixel side-play amount of reflector points, according to pixel side-play amount, determines the real offset of described corner reflector point.
Optionally, from based on the whole image of seasonal effect in time series being chosen with reference to image, determine in whole image except with reference to shadow
The system deviation amount between residue image and reference image outside Xiang, including:
Based on time reference line, Space Baseline and the optimum main image choosing method of doppler centroid difference, from whole shadows
Select in Xiang with reference to image；Based on formula 1, determine with reference to image with surplus in conjunction with satellite orbit data or big window intensity method of correlation
System deviation amount between remaining image,
Wherein (Ri,Ai) represent target in the i-th scape image middle-range descriscent and orientation to position, Δ (R, A)iRepresent same item
It is marked on the i-th scape image and with reference to alternate position spike between image, a0,…,a5；b0,…,b5For coefficient to be asked.
Optionally, choosing corresponding function model with reference to image is
Wherein, K represents the quantity of image,Represent interfere to vertical parallax, Tk,mRepresent interfere to time base
Line,Represent that doppler centroid is poor, Bc、TcAnd fcBeing respectively corresponding critical condition, α, β and θ are respectively corresponding finger
The number factor, γmRepresenting overall correlation coefficient, minimum is optimum.
Optionally, according to longitude and latitude and the height of known corner reflector point, determine based on RD geo-location model
Described corner reflector point position in whole images, is set to the initial position of described corner reflector point by described position, based on
Initial position, opens centered by initial position in whole images and takes a certain size window and carry out over-sampling process, and acquisition processes
After intensity peak position, including:
Determine longitude and latitude and the height of corner reflector point, based on location model, determine that described corner reflector point is at whole shadows
Initial position in Xiang, described initial position includes initial line number and initial column number；
Centered by initial position, in reference to image, choose the window of default size, use lowpass digital filter pair
All image carries out over-sampling process in orientation the most respectively to distance, obtains after sampling with reference to the intensity peak position in image
Put；
Intensity peak position in residue image after same acquisition sampling.
Optionally, the function model determining described corner reflector point initial position in whole images is
Wherein R is the satellite oblique distance to ground,It is respectively ground point and the coordinate vector of satellite point,For satellite
Velocity,For satellite to ground distance vector, λ is wavelength, fdopFor the Doppler frequency of P point signal, XP、YP、ZPFor P
The three-dimensional coordinate of point, h is the geodetic height of P point, a and b is respectively major semiaxis and the semi-minor axis of ellipsoid.
Optionally, based on the intensity peak got, in conjunction with before the system deviation amount that determines carry out difference processing, determine
The pixel side-play amount of described corner reflector point, according to pixel side-play amount, determines the real offset of described corner reflector point, bag
Include:
Based on intensity peak, in conjunction with before the residue image that determines relative to the system deviation amount with reference to image, will be initial
Position and system deviation amount carry out difference processing, obtain the renewal position of described corner reflector point；
From whole images, choose the first scape image, the described corner reflector o'clock position in the first scape image has been set to
Beginning position, the intensity peak position extracted in whole image in the intensity peak position of corner reflector point, with the first scape image is entered
Row difference processing, determine described corner reflector point distance to orientation pixel side-play amount upwards；
Based on described corner reflector point distance to orientation Pixel size upwards, determine described corner reflector point away from
Descriscent and orientation actual displacement amount upwards.
The technical scheme that the present invention provides has the benefit that
By low pass filter, corner reflector strength signal is carried out the above over-sampling of Radix Achyranthis Bidentatae, ask on time series image
The intensity peak location variation identified, as corner reflector point deformation data.This invention solves InSAR technical monitoring
The difficult point of mining area big magnitude Deformation Field, and monitoring accuracy is higher compared with SAR strength offsets amount track algorithm.
Accompanying drawing explanation
In order to be illustrated more clearly that technical scheme, the required accompanying drawing used in embodiment being described below
It is briefly described, it should be apparent that, the accompanying drawing in describing below is only some embodiments of the present invention, general for this area
From the point of view of logical technical staff, on the premise of not paying creative work, it is also possible to obtain other accompanying drawing according to these accompanying drawings.
Fig. 1 is the techniqueflow schematic diagram that the present invention provides；
The acquisition image coverage that Fig. 2 (a) provides for the present invention；
The corner reflector that Fig. 2 (b) provides for present invention display on SAR intensity map；
The TerraSAR data list that the experiment that Fig. 2 (c) provides for the present invention uses；
The identification of the 18207-1 corner reflector that Fig. 3 (a) provides for the present invention；
The identification of the 18207-2 corner reflector that Fig. 3 (b) provides for the present invention；
The 18207-1 corner reflector point contrast that Fig. 4 (a) provides for the present invention；
The 18207-2 corner reflector point contrast that Fig. 4 (b) provides for the present invention；
Four kinds of SAR methods that Fig. 5 (a) provides for the present invention are with GPS difference average ratio relatively；
Four kinds of SAR methods that Fig. 5 (b) provides for the present invention compare with GPS difference root-mean-square error；
The precision comparison of two kinds of offset methods that Fig. 5 (c) provides for the present invention；
Two kinds of concrete numerical value of offset method precision that Fig. 5 (d) provides for the present invention.
Detailed description of the invention
Structure and advantage for making the present invention are clearer, make the structure of the present invention further below in conjunction with accompanying drawing
Describe.
Embodiment one
The invention provides a kind of mining area based on manual corner reflector side-play amount deformation monitoring method, as it is shown in figure 1, this
Monitoring method, including:
11, from based on the whole image of seasonal effect in time series being chosen with reference to image, determine in whole image except with reference to image it
System deviation amount between outer residue image and reference image.
12, according to longitude and latitude and the height of known corner reflector point, determine described based on RD geo-location model
Corner reflector point position in whole images, is set to the initial position of described corner reflector point, based on initially by described position
Position, opens centered by initial position in whole images and takes a certain size window and carry out over-sampling process, after acquisition processes
Intensity peak position.
13, based on the intensity peak got, in conjunction with before the system deviation amount that determines carry out difference processing, determine described
The pixel side-play amount of corner reflector point, according to pixel side-play amount, determines the real offset of described corner reflector point.
In step 11, from based on the whole image of seasonal effect in time series being chosen with reference to image, determine in whole image except reference
The system deviation amount between residue image and reference image outside image, including:
101, based on time reference line, Space Baseline and the optimum main image choosing method of doppler centroid difference, from entirely
Portion's image is selected with reference to image.
102, based on formula 1, determine with reference to image and residue image in conjunction with satellite orbit data or image intensity method of correlation
Between system deviation amount,
Wherein (Ri,Ai) represent target in the i-th scape image middle-range descriscent and orientation to position, Δ (R, A)iRepresent same item
It is marked on the i-th scape image and with reference to alternate position spike between image, a0,…,a5；b0,…,b5For coefficient to be asked.
Optionally, in order to realize accurately choosing with reference to image, the Selection of Function model of reference is
Wherein, K represents the quantity of image,Represent interfere to vertical parallax, Tk,mRepresent interfere to time base
Line,Represent that doppler centroid is poor, Bc、TcAnd fcBeing respectively corresponding critical condition, α, β and θ are respectively corresponding finger
The number factor, γmRepresenting overall correlation coefficient, minimum is optimum.
Here Selection of Function model use Kampes et al. propose consider time reference line, Space Baseline and
The optimum main image choosing method of doppler centroid difference, chooses a scape image as ginseng from time series n scape SAR image
Examine image.
In step 12, according to longitude and latitude and the height of known corner reflector point, true based on RD geo-location model
Fixed described corner reflector point position in whole images, is set to the initial position of described corner reflector point, base by described position
In initial position, open centered by initial position in whole images and take a certain size window and carry out over-sampling process, at acquisition
Intensity peak position after reason, including:
201, determine longitude and latitude and the height of corner reflector point, based on location model, determine that described corner reflector point is entirely
Initial position in portion's image, described initial position includes initial line number and initial column number.According to known corner reflector point longitude and latitude
Degree and geodetic height, based on tight RD location model, calculate corner reflector point position (ranks number) in all images, anti-as angle
Emitter point initial position.Longitude and latitude and elevation information in this step can pass through global positioning system (Global
Positioning System, GPS) obtain, it is also possible to according to other alignment systems, such as GPS
(GLOBAL NAVIGATION SATELLITE SYSTEM, GLONASS), Beidou satellite navigation system (BeiDou
Navigation Satellite System, BDS) etc. alignment system obtain.Location model is typical distance-Doppler
(RD:Range-Doppler) model.
202, centered by initial position, in reference to image, choose the window of default size, use low-pass digital filter
Device the most respectively carries out over-sampling process in orientation to distance to whole images, obtains the intensity in reference to image after sampling
Peak.Centered by corner reflector point initial position, SAR intensity image takes a certain size window and (typically goes up bottom left
The right side respectively takes 3,5 or 7, is as the criterion with given corner reflector positional accuracy, along with the raising window of precision can suitably reduce), use low
Logical FIR (Finite Impulse Response) digital filter to intensity map in window in orientation to distance to entering respectively
Row over-sampling (in view of the precision of the location recognition of corner reflector intensity peak, over-sampling multiple >=100), after asking for over-sampling
Intensity peak position also stores.All images are carried out successively aforesaid operations, finally tries to achieve corner reflector in all images
The position (R (i), Azi (i)) at peak value place0, i=1,2 ... n.
203, the intensity peak position in residue image after same acquisition sampling.By step 202 is asked for (R (i),
Azi(i))0Bring (formula 2) into, try to achieve corner reflector point position respectively from image relative to the system deviation amount of main image, then
By the position of corner reflector identified for the first time and system deviation amount difference, thus obtain registration to all images after main image
In new corner reflector position (R (i), Azi (i))1, i=1,2 ... n.
Wherein it is determined that the function model of the initial position that described corner reflector point is in whole images is
Wherein R is the satellite oblique distance to ground,It is respectively ground point and the coordinate vector of satellite point,For satellite
Velocity,For satellite to ground distance vector, λ is wavelength, fdopFor the Doppler frequency of P point signal, XP、YP、ZPFor P
The three-dimensional coordinate of point, h is the geodetic height of P point, a and b is respectively major semiaxis and the semi-minor axis of ellipsoid.
In step 13, based on the intensity peak got, in conjunction with before the system deviation amount that determines carry out difference processing, really
The pixel side-play amount of fixed described corner reflector point, according to pixel side-play amount, determines the real offset of described corner reflector point, bag
Include:
301, based on intensity peak, in conjunction with before the residue image that determines relative to reference to the system deviation amount of image, will
Initial position and system deviation amount carry out difference processing, obtain the renewal position of described corner reflector point.
302, from whole images, the first scape image is chosen, the described corner reflector o'clock position in the first scape image is fixed
For original position, extract the intensity peak position in the intensity peak position of corner reflector point, with the first scape image in whole image
Put and carry out difference processing, determine described corner reflector point distance to orientation pixel side-play amount upwards.
303, based on described corner reflector point distance to orientation Pixel size upwards, determine described corner reflector point
Distance to orientation actual displacement amount upwards.
In force, in the first scape image, corner reflector point position is as the original position of this corner reflector, by other
Difference is upwards done with the first scape image in distance to orientation in the corner reflector intensity peak position of all image identifications, asks with this
Obtain this corner reflector point distance to orientation to pixel side-play amount.Computing formula is as follows:
(Δ R (i), Δ Azi (i))=(R (i), Azi (i))1-(R(1),Azi(1))1, i=2 ... n.
According to corner reflector position distance to orientation to Pixel size, ask for corner reflector distance to and side
Position to actual displacement amount.Computing formula is as follows:
(DR(i),DAzi(i))=(Δ R (i), Δ Azi (i)) diag (r, azi), i=2 ... n.
Wherein (r, azi) be image distance to orientation to Pixel size, (DR(i),DAzi(i)) represent the i-th time
Sequence corner reflector respectively distance to orientation to deformation quantity.
It should be noted that the whole images including reference image and residue image are that synthetic aperture radar gets
Image.Needing at mining area established angle reflector before this, corner reflector is arranged on around weak scatterers as far as possible.
The present invention use the mode of " point-to-point " eliminate system deviation amount, it is to avoid Image registration resampling interpolation causes
The impact of phase error, improves the accuracy of identification of corner reflector point peak value, when accurately asking for corner reflector deformation to greatest extent
Between sequence.It should be noted that when asking for system deviation amount multinomial, deformation region mask will be fallen, be not involved in system inclined
Shifting amount fitting of a polynomial.It addition, system deviation amount is polynomial asks for may be used without other method: interference coefficient method, maximum
Interference spectrum method, Phase contrast averagely fluctuate function method etc..Additionally, this invention make use of corner reflector geographical position more accurately
Confidence ceases, if known corner reflector approximate range, then can be by the visual interpretation of SAR intensity image be obtained at the beginning of corner reflector
The beginning position feasibility of (corner reflector strong scattering feature determine the method).
The invention provides a kind of side-play amount algorithm based on manual corner reflector identification and carry out mining area big magnitude deformation prison
The method surveyed.Utilize corner reflector strong scattering feature on SAR intensity map, and the impulse response meeting 2D SINC function is special
Levy, by low-pass FIR filter, corner reflector strength signal is carried out the above over-sampling of Radix Achyranthis Bidentatae, ask for knowing on time series image
Other intensity peak location variation, as corner reflector point deformation data.This invention solves InSAR technical monitoring ore deposit
The difficult point of district's big magnitude Deformation Field, and monitoring accuracy is higher compared with SAR strength offsets amount track algorithm.
In order to improve the readability of the present invention, specific embodiments is to be provided with Seam, Gujiao Mining Area, the Taiyuan work of corner reflector
Describe in detail as a example by making face monitoring.
In survey region with a varied topography, totally in medium mountain area, locally landform is precipitous, rock partial denudation, gentle slope and low
Having waste soil at Wa, much the most " V "-shaped cheuch is.Experimental Area is Lan Tun the ore deposit 18207 and 18a203 work surface of Seam, Gujiao Mining Area,
Two work surface mining of adjacent, recovery time has certain overlap.18207 working face extraction times were-2012 years 10 in May, 2011
Month, mining depth is about 354～355m.Coal seam thickness 2.3～3.65, moves towards long 1254m, is inclined to wide 201 meters, and seam inclination is about
It it is 5 °；The 18a203 working face extraction time is in June, 2012 (in June, 2014), mining depth 313～493m so far, averagely adopts
Thickness about 1.46 meters, moves towards long 817 meters, is inclined to wide 198 meters, and seam inclination is about 4 °.
Experimental data uses the 7 scape TerraSAR from April, 2012 in July, 2012 that DLR (space office of Germany) provides
Data, distance to and azimuth resolution about 2 meters, image coverage is shown in Fig. 2 (a), in Experimental Area install corner reflector
Display on SAR intensity map is shown in that Fig. 2 (b), detailed data list are shown in Fig. 2 (c).
It is as follows that the present invention is embodied as step:
(1) the choosing of main image: combine formula (1), finally chooses image 2012-04-04 as with reference to image.
(2) side-play amount multinomial is calculated: use described in step 2, based on DIFF module in GAMMA software, finally try to achieve
Remaining 6 scape image and the system deviation amount multinomial with reference to image 2012-04-04.The step for carry out based on cutting image, ginseng
Examine the cutting scope unsuitable excessive (when region is excessive, being readily incorporated the system deviation amount error that regional area calculates) of image, this
Shown in secondary test main image term of reference such as Fig. 2 (b), in figure, red dotted line frame is research position, mining area, and intensity is ginseng after cutting
Examine image capturing range.
(3) corner reflector point location in image: according to GPS measure 18207-1 corner reflector dot position information (longitude,
Latitude and elevation) and each image parameters file (comprising the essential informations such as image track), primary Calculation corner reflector point shadow
Position (pixel position rounds) in Xiang.
(4) corner reflector position resolves: in (3) step centered by corner reflector initial position, upper and lower at SAR intensity image
Left and right respectively takes 7 pixels, uses low pass FIR (Finite Impulse Response) digital filter to intensity map in window
In orientation to distance to carrying out 300 times of over-samplings respectively, intensity peak position storing, 18207-after asking for over-sampling
Shown in 1 corner reflector point identification such as Fig. 3 (a).
(5) (R (i), the Azi (i)) that will ask in step 40Bring (formula 2) into, try to achieve corner reflector point position respectively from
Image is relative to the system deviation amount of main image, then by the position of corner reflector identified for the first time and system deviation amount difference,
Corner reflector position new in all images after thus obtaining registration extremely main image, the 18207-1 position redefined.
(6) using first scape image definition corner reflector position as the original position of this corner reflector, other is owned
Difference is upwards done with the first scape image in distance to orientation in the corner reflector intensity peak position of image identification, tries to achieve this with this
Corner reflector point distance to orientation to pixel side-play amount.Final 18207-1 point pixel side-play amount oblique distance depends on to difference value
Secondary it is: 0.000,0.043,0.087,0.103,0.120,0.100,0.107,0.167；Orientation is followed successively by difference value:
0.00 ,-0.027 ,-0.037 ,-0.043 ,-0.063 ,-0.170 ,-0.147 ,-0.127.
(7) according to corner reflector position oblique distance to orientation to unit pixel represent distance, ask for corner reflector
Distance to orientation to the real offset of movement.The TerraSAR image oblique distance that this test is selected to resolution is
0.909m, azimuth resolution is 1.965m, and the corner reflector thus calculated is followed successively by deformation quantity in oblique distance: 0m, 0.085m,
0.171m, 0.202m, 0.236m, 0.197m, 0.210m；Azimuth resolution is: 0m ,-0.053m ,-0.073m ,-
(, on the occasion of representing away from SAR satellite, negative value is near SAR satellite for oblique distance for 0.085m ,-0.124m ,-0.334m ,-0.289m；Orientation
To on the occasion of represent along satellite flight direction, negative value represent deviate from satellite flight direction).
(8) repeating step 3-7, again try to achieve 18207-2 point corner reflector deformation data, wherein 18207-2 angle is anti-
Shown in emitter point location recognition such as Fig. 3 (b).
For verifying the application effect of the present invention, contrasted with the additive method of InSAR/SAR technology, including sequential
SBAS technology, D-InSAR mono-interferogram superimposing technique, IPTA technology based on coherent point and strength offsets amount tracking technique.
Five kinds of methods corner reflector point position deformation data contrast as shown in Figure 4 with GPS sequential (figure compares for vertical shape
Become, be to project to contrast with GPS vertical strain in vertical direction to deformation by the oblique distance that above-mentioned 5 kinds of technology calculate), Fig. 4
For SBAS, D-InSAR phase stacking technique (Chain), side-play amount tracking technique (Offset-tracking) and based on corner reflection
The timing offset technology of device identification technology goes out deformation data and GPS monitoring knot in study area at two corner reflector points
Really comparison diagram.The accuracy assessment result of Fig. 4 is as shown in Figure 5.Fig. 4 and Fig. 5 can be seen that when big magnitude deformation, based on partially
Move metering method better than InSAR monitoring method based on phase place generally, with GPS monitor value result closer to, and anti-based on angle
The offset method of emitter identification is better than strength offsets amount tracking on the whole, its average closest to GPS observed quantity, mean square
Root error is minimum, preferably supplements the deficiency of InSAR and SAR strength offsets amount tracking technique, illustrates that the present invention is feasible
Reliably.
The foregoing is only embodiments of the invention, not in order to limit the present invention, all in the spirit and principles in the present invention
Within, any modification, equivalent substitution and improvement etc. made, should be included within the scope of the present invention.
Claims (6)
1. mining area based on a manual corner reflector side-play amount deformation monitoring method, it is characterised in that described based on artificial angle
The mining area deformation monitoring method of reflectors offset, including:
From based on the whole image of seasonal effect in time series being chosen with reference to image, determine the residue in addition to reference to image in whole image
System deviation amount between image and reference image；
Longitude and latitude according to known corner reflector point and height, determine described corner reflection based on RD geo-location model
Device point position in whole images, is set to the initial position of described corner reflector point by described position, based on initial position,
All images are opened centered by initial position and take a certain size window and carry out over-sampling process, acquisition process after intensity peak
Position；
Based on the intensity peak got, in conjunction with before the system deviation amount that determines carry out difference processing, determine described corner reflection
The pixel side-play amount of device point, according to pixel side-play amount, determines the real offset of described corner reflector point.
Mining area based on manual corner reflector side-play amount the most according to claim 1 deformation monitoring method, it is characterised in that
From based on the whole image of seasonal effect in time series being chosen with reference to image, determine the residue image in addition to reference to image in whole image
And the system deviation amount between reference image, including:
Based on time reference line, Space Baseline and the optimum main image choosing method of doppler centroid difference, from whole images
Select with reference to image；
Based on formula 1, determine with reference between image and residue image being in conjunction with satellite orbit data or image intensity method of correlation
System side-play amount,
Wherein (Ri,Ai) represent target in the i-th scape image middle-range descriscent and orientation to position, Δ (R, A)iRepresent that target of the same name exists
Alternate position spike between i-th scape image and reference image, a0,…,a5；b0,…,b5For coefficient to be asked.
Mining area based on manual corner reflector side-play amount the most according to claim 2 deformation monitoring method, it is characterised in that:
The function model choosing correspondence with reference to image is
Wherein, K represents the quantity of image,Represent interfere to vertical parallax, Tk,mRepresent interfere to time reference line,Represent that doppler centroid is poor, Bc、TcAnd fcBeing respectively corresponding critical condition, α, β and θ are respectively corresponding index
The factor, γmRepresenting overall correlation coefficient, minimum is optimum.
Mining area based on manual corner reflector side-play amount the most according to claim 1 deformation monitoring method, it is characterised in that
Longitude and latitude according to known corner reflector point and height, determine described corner reflector point based on RD geo-location model
Position in whole images, is set to the initial position of described corner reflector point by described position, based on initial position, all
Image is opened centered by initial position and takes a certain size window and carry out over-sampling process, acquisition process after intensity peak position
Put, including:
Determine longitude and latitude and the height of corner reflector point, based on location model, determine that described corner reflector point is in whole images
Initial position, described initial position includes initial line number and initial column number；
Centered by initial position, in reference to image, choose the window of default size, use lowpass digital filter to all
Image carries out over-sampling process in orientation the most respectively to distance, obtains the intensity peak position in reference to image after sampling
Put；
Intensity peak position in residue image after same acquisition sampling.
Mining area based on manual corner reflector side-play amount the most according to claim 4 deformation monitoring method, it is characterised in that:
The function model determining described corner reflector point initial position in whole images is
Wherein R is the satellite oblique distance to ground,It is respectively ground point and the coordinate vector of satellite point,For satellite velocities
Vector,For satellite to ground distance vector, λ is wavelength, fdopFor the Doppler frequency of P point signal, XP、YP、ZPFor P point
Three-dimensional coordinate, h is the geodetic height of P point, a and b is respectively major semiaxis and the semi-minor axis of ellipsoid.
Mining area based on manual corner reflector side-play amount the most according to claim 1 deformation monitoring method, it is characterised in that
Based on the intensity peak got, in conjunction with before the system deviation amount that determines carry out difference processing, determine described corner reflector point
Pixel side-play amount, according to pixel side-play amount, determine the real offset of described corner reflector point, including:
Based on intensity peak, in conjunction with before the residue image that determines relative to the system deviation amount with reference to image, by initial position
Carry out difference processing with system deviation amount, obtain the renewal position of described corner reflector point；
From whole images, choose the first scape image, the described corner reflector o'clock position in the first scape image is set to start bit
Putting, the intensity peak position extracted in whole image in the intensity peak position of corner reflector point, with the first scape image is poor
Divisional processing, determine described corner reflector point distance to orientation pixel side-play amount upwards；
Based on described corner reflector point distance to orientation Pixel size upwards, determine described corner reflector point distance to
With orientation actual displacement amount upwards.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610326486.7A CN106093938A (en) | 2016-05-17 | 2016-05-17 | A kind of mining area based on manual corner reflector side-play amount deformation monitoring method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610326486.7A CN106093938A (en) | 2016-05-17 | 2016-05-17 | A kind of mining area based on manual corner reflector side-play amount deformation monitoring method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106093938A true CN106093938A (en) | 2016-11-09 |
Family
ID=57230039
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610326486.7A Pending CN106093938A (en) | 2016-05-17 | 2016-05-17 | A kind of mining area based on manual corner reflector side-play amount deformation monitoring method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106093938A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106526593A (en) * | 2016-12-19 | 2017-03-22 | 国家测绘地理信息局卫星测绘应用中心 | Sub-pixel-level corner reflector automatic positioning method based on SAR rigorous imaging model |
| CN108446516A (en) * | 2018-04-04 | 2018-08-24 | 长安大学 | A kind of construction fracture and surface subsidence deformation decomposition method based on inversion of the force |
| CN111142119A (en) * | 2020-01-10 | 2020-05-12 | 中国地质大学（北京） | Mine geological disaster dynamic identification and monitoring method based on multi-source remote sensing data |
| CN111366925A (en) * | 2020-03-27 | 2020-07-03 | 长安大学 | SAR offset two-dimensional deformation time sequence calculation method and system |
| CN111522005A (en) * | 2020-06-05 | 2020-08-11 | 中国有色金属长沙勘察设计研究院有限公司 | Deformation monitoring and terrain reconstruction method |
| CN112986949A (en) * | 2021-04-21 | 2021-06-18 | 北京东方至远科技股份有限公司 | SAR high-precision time sequence deformation monitoring method and device for diagonal reflector |
| CN113296136A (en) * | 2020-02-24 | 2021-08-24 | 香港理工大学深圳研究院 | Multi-target deformation monitoring method and device and receiver |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001221856A (en) * | 2000-02-09 | 2001-08-17 | Univ Tokyo | Doppler earth speedometer |
| CN103217665A (en) * | 2013-04-11 | 2013-07-24 | 河海大学 | Measurable and orientable foundation radar angle reflector |
| CN103576149A (en) * | 2013-06-05 | 2014-02-12 | 河海大学 | Foundation interference radar three-dimensional deformation extraction method based on amplitude information |
| CN104122553A (en) * | 2014-07-23 | 2014-10-29 | 中国国土资源航空物探遥感中心 | Regional ground settlement monitoring method based on multiple track and long strip CTInSAR (coherent target synthetic aperture radar interferometry) |
-
2016
- 2016-05-17 CN CN201610326486.7A patent/CN106093938A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001221856A (en) * | 2000-02-09 | 2001-08-17 | Univ Tokyo | Doppler earth speedometer |
| CN103217665A (en) * | 2013-04-11 | 2013-07-24 | 河海大学 | Measurable and orientable foundation radar angle reflector |
| CN103576149A (en) * | 2013-06-05 | 2014-02-12 | 河海大学 | Foundation interference radar three-dimensional deformation extraction method based on amplitude information |
| CN104122553A (en) * | 2014-07-23 | 2014-10-29 | 中国国土资源航空物探遥感中心 | Regional ground settlement monitoring method based on multiple track and long strip CTInSAR (coherent target synthetic aperture radar interferometry) |
Non-Patent Citations (3)
| Title |
|---|
| 李平湘 等: "《雷达干涉测量原理与应用》", 31 December 2006 * |
| 牛玉芬: "SAR/InSAR技术用于矿区探测与形变监测研究", 《中国优秀硕士学位论文全文数据库工程科技I辑》 * |
| 陈强: "基于永久散射体雷达差分干涉探测区域地表形变的研究", 《万方数据库》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106526593A (en) * | 2016-12-19 | 2017-03-22 | 国家测绘地理信息局卫星测绘应用中心 | Sub-pixel-level corner reflector automatic positioning method based on SAR rigorous imaging model |
| CN106526593B (en) * | 2016-12-19 | 2019-01-11 | 国家测绘地理信息局卫星测绘应用中心 | Sub-pixel-level corner reflector automatic positioning method based on the tight imaging model of SAR |
| CN108446516A (en) * | 2018-04-04 | 2018-08-24 | 长安大学 | A kind of construction fracture and surface subsidence deformation decomposition method based on inversion of the force |
| CN108446516B (en) * | 2018-04-04 | 2021-09-07 | 长安大学 | Structural fracture and ground settlement deformation decomposition method based on force source inversion |
| CN111142119A (en) * | 2020-01-10 | 2020-05-12 | 中国地质大学（北京） | Mine geological disaster dynamic identification and monitoring method based on multi-source remote sensing data |
| CN111142119B (en) * | 2020-01-10 | 2021-08-17 | 中国地质大学（北京） | Mine geological disaster dynamic identification and monitoring method based on multi-source remote sensing data |
| CN113296136A (en) * | 2020-02-24 | 2021-08-24 | 香港理工大学深圳研究院 | Multi-target deformation monitoring method and device and receiver |
| CN113296136B (en) * | 2020-02-24 | 2024-03-19 | 香港理工大学深圳研究院 | Multi-target deformation monitoring method, device and receiver |
| CN111366925A (en) * | 2020-03-27 | 2020-07-03 | 长安大学 | SAR offset two-dimensional deformation time sequence calculation method and system |
| CN111366925B (en) * | 2020-03-27 | 2022-11-22 | 长安大学 | SAR offset two-dimensional deformation time sequence calculation method and system |
| CN111522005A (en) * | 2020-06-05 | 2020-08-11 | 中国有色金属长沙勘察设计研究院有限公司 | Deformation monitoring and terrain reconstruction method |
| CN112986949A (en) * | 2021-04-21 | 2021-06-18 | 北京东方至远科技股份有限公司 | SAR high-precision time sequence deformation monitoring method and device for diagonal reflector |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106093938A (en) | A kind of mining area based on manual corner reflector side-play amount deformation monitoring method | |
| Telling et al. | Review of Earth science research using terrestrial laser scanning | |
| Barbarella et al. | Monitoring of large landslides by Terrestrial Laser Scanning techniques: field data collection and processing | |
| Baltsavias et al. | Digital surface modelling by airborne laser scanning and digital photogrammetry for glacier monitoring | |
| Liu et al. | A complete high-resolution coastline of Antarctica extracted from orthorectified Radarsat SAR imagery | |
| CN102323582B (en) | Autonomous orbit determination method for satellite based on synthetic aperture radar | |
| CN114444158B (en) | Underground roadway deformation early warning method and system based on three-dimensional reconstruction | |
| Fokker et al. | Subsidence in the Dutch Wadden Sea | |
| CN106556834A (en) | It is a kind of to concentrate accurate extraction Ground Vertical Deformation method from two parallel orbit SAR datas | |
| CN110044327A (en) | A kind of infrastructure Monitoring method of the subsidence and system based on SAR data and GNSS data | |
| CN109471104A (en) | A method of obtaining earth's surface three-dimensional amount of movement from two parallel orbit SAR datas | |
| Battaglia et al. | Monitoring volcanic deformation | |
| Liu et al. | Correction of positional errors and geometric distortions in topographic maps and DEMs using a rigorous SAR simulation technique | |
| Baker | Investigating the dynamics of basaltic volcano magmatic systems with space geodesy | |
| Barratt et al. | Survey at Wroxeter using carrier phase, differential GPS surveying techniques | |
| Whitehead | An integrated approach to determining short-term and long-term patterns of surface change and flow characteristics for a polythermal arctic glacier | |
| Catalão et al. | Photogrammetric analysis of the coastal erosion in the Algarve (Portugal)/Analyse photogrammétrique de l'érosion côtiere en Algarve (Portugal) | |
| Nolan et al. | Which are the highest peaks in the US Arctic? Fodar settles the debate | |
| Xu | Three-dimensional virtual geology: photorealistic outcrops, and their acquisition, visualization and analysis | |
| Hlaváčová et al. | Sentinel-1 INSAR processing of corner reflector information in the northern-bohemian coal basin | |
| CN114199189B (en) | Mining subsidence monitoring method combining unmanned plane and DINSAR technology | |
| White et al. | Evolution of the USACE National Coastal Mapping Program and Implementation of Small Unmanned Aerial Systems | |
| Luppichini et al. | A METHOD BASED ON BEACH PROFILE ANALYSIS FOR SHORELINE IDENTIFICATION | |
| Rowe | An examination of structural constraints on rockfall behaviour using LiDAR data | |
| CN114384585A (en) | Mountain land underground explosion absolute positioning method based on relative position and minimum burial depth |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161109 |
|
| WD01 | Invention patent application deemed withdrawn after publication |