CN104239419A - Discontinuous subnet connecting method and device oriented to time sequence InSAR - Google Patents
Discontinuous subnet connecting method and device oriented to time sequence InSAR Download PDFInfo
- Publication number
- CN104239419A CN104239419A CN201410410540.7A CN201410410540A CN104239419A CN 104239419 A CN104239419 A CN 104239419A CN 201410410540 A CN201410410540 A CN 201410410540A CN 104239419 A CN104239419 A CN 104239419A
- Authority
- CN
- China
- Prior art keywords
- subnet
- subnetwork
- discontiguous
- subnets
- discontinuous
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Data Exchanges In Wide-Area Networks (AREA)
Abstract
The invention provides a discontinuous subnet connecting method and device oriented to a time sequence InSAR. The method comprises the following steps: numbering all discontinuous Delaunay triangulation network subnets and apexes thereof in a regional radar image; finding all point targets and subnets to which the point targets belong in a preset initial step length L specific to four boundary points of each subnet; calculating the distances between all point targets and the boundary points of a current subnet, sequencing the distances from small to large, calculating the maximum model correlation coefficient of each edge in sequence, judging that the edge is a valid edge when the correlation coefficient is greater than a given threshold value, otherwise, continually calculating a next edge; traversing all the subnets, and then recounting the number of subnets, wherein if only one subnet exists, expansion link calculation is not required; if more than 1 subnet is counted, repeating the steps to perform extension link by taking 2L, 3L, 4L, 5L, ..., nL as step lengths till circulation stops. By adopting the method and the device, the total quantity of the subnets can be reduced gradually, rapid connection of discontinuous subnets is realized, and the surface deformation monitoring in a large area is facilitated.
Description
Technical field
The present invention relates to InSAR (Synthetic Aperture Radar Interferometry, synthetic aperture radar interferometry) technical field, particularly relate to a kind of discontiguous subnetwork method of attachment towards time series InSAR and device.
Background technology
Synthetic-aperture radar (SAR, Synthetic Aperture Radar) is the most important earth observation technology grown up the 1950's.But SAR technology can only obtain the two-dimensional signal of earth's surface object, lack the ability obtaining impact point elevation information and monitoring objective miniature deformation.Synthetic aperture radar interferometry technology (the InSAR that interfere measurement technique and SAR combine with technique are formed, Synthetic Aperture Radar Interferometry) provide the completely new approach obtaining Three Dimensional Ground information, it is observed by two slave antennas or by common antenna twice parallel observation simultaneously, obtain the complex pattern pair of the same view in ground, according to the phase differential of ground each point in two width complex patterns, draw the path length difference of each point microwave in twice imaging, thus obtain the elevation information of terrain object.Based on the development of InSAR technology, Differential Interferometric Synthetic Aperture Radar (Differential Synthetic Aperture Radar Interferometry, DInSAR) technology be to two width more than interferogram or a width ground digital elevation model (DEM is added to a width interferogram, Digital Elevation Model) carry out a kind of technology of reprocessing, the interferometric phase that it causes by removing landform, obtain the information about Ground Deformation, in fields such as volcano monitoring, earthquake displacement measurement, land subsidence, there is great application potential.
DInSAR technology based on repeat track is easily done by space dephasing, time dephasing is done and the impact of the factor such as atmospheric interference, is difficult to the practical application carrying out normalization.The reason causing space dephasing to do comprises large vertical parallax and large deformation gradient.Time dephasing is done, and causes because twice imaging moment environment inconsistent causes the scattering properties of scatterer in the pixel of same ground to change.Variable atmospheric conditions may cause interfering phase delay inconsistent on image at two width, thus cause distortion measurement error.
In order to overcome these restrictions of traditional DInSAR technology, from the end of the nineties in last century, some new InSAR treatment technologies are suggested.The common feature of these technology is: process based on time series SAR image, the object of process is not whole pixels of view picture image, but wherein there is stable scattering properties thus high relevant pixel subset can be kept in long period interval, namely high coherent point.These technology may be summarized to be two classes generally: interfere (Permanent Scatterer or Persistent Scatterer Interferometry with Permanent scatterers, or PS-InSAR) be the single main image time series InSAR technology of representative and the how main image time series InSAR technology that is representative with Small Baseline Subset technology (Small baseline subset interferometry, or SBAS InSAR).For sake of convenience, these two kinds of technology are referred to as time series InSAR technology.Time series InSAR technology all has good immunity to above-mentioned three kinds of restraining factors, replaced at present traditional DInSAR technology in volcano, earthquake, landslide, the field such as land subsidence widely applied.Because urban area has intensive natural point target, as street lamp, roof etc., therefore time series InSAR technology is most widely used in the measurement of urban ground subsidence.
In time series InSAR technology, after having extracted all point targets, usually utilize Delaunay triangulation network to connect all point targets, set up linear phase model solve deformation parameter by consecutive point phase difference.But affect by the restriction of air correlation distance and model coefficient of coherence threshold value, Delaunay triangulation network usually can not intactly connect all point targets.When point target quality is good and density is larger, as city, each point target can be connected preferably by Delaunay networking, realize resolving of deformation parameter; But when point target is second-rate or relatively far apart, as non-city, easily there is a large amount of isolated subnet, at this moment just cannot resolve deformation data.
For the problems referred to above, current researchist proposes two kinds of solutions, and one is point target free connecting method, the institute in certain distance is connected a little between two, thus greatly adds the limit number of the triangulation network; One is complex network method, first it utilize Delaunay triangulation network to the quick networking of all point targets, extracts the deformation parameter of city point of density target, and then connecting within the scope of air correlation distance does not allly have estimative point target, set up complex network, recalculate deformation parameter.Above-mentioned two kinds of methods all can improve the annexation between point target, but it is evident that, the calculated amount that they bring is very huge, is not suitable for large regions mass data processing.How fast and steadily tie point targeted subnet, to realize the monitoring of large regions Ground Deformation, this is those skilled in the art's key issues urgently to be resolved hurrily.
Summary of the invention
The embodiment of the present invention provides a kind of discontiguous subnetwork method of attachment towards time series InSAR and device, to solve the sane connectivity problem of discontinuous Delaunay triangulation network, thus better realizes the monitoring of large regions Ground Deformation.
On the one hand, embodiments provide a kind of discontiguous subnetwork method of attachment towards time series InSAR, the described discontiguous subnetwork method of attachment towards time series InSAR comprises:
The all discontinuous Delaunay triangulation network subnet of one band radar image and described all discontinuous Delaunay triangulation network subnet summits are numbered;
To four frontier points of each subnet in described all discontinuous Delaunay triangulation network subnets, search all point targets within the scope of preset initial step length L and the subnet belonging to described all point targets;
Calculate the distance of described all point targets and current sub network frontier point, and by described apart from ascending sequence, calculate the maximum model coefficient of coherence on every bar limit successively, think the effective edge of this limit for being connected with current sub network by the subnet belonging to this point target when this maximum model coefficient of coherence is greater than given threshold value, otherwise continue to calculate next limit;
After having traveled through all subnets, again add up subnet number, if this subnet number only has 1, whole subnet successful connection all has been described, then connected without the need to expansion and calculate; If this subnet number is more than 1, the isolated subnet be not connected in addition is described, more successively with 2L, 3L, 4L, 5L ... nL is step-length, repeats above-mentioned steps and carries out expansion connection, until circulation stops, now all-ones subnet has all completed to expand and connected; If also there is the subnet of some more than nL of being separated by with major networks, the subnet of described more than nL of being separated by with major networks is abandoned, wherein n be greater than 1 natural number.
The embodiment of the present invention, by this multi-layer extended mode, can gradually reduce subnet sum, thus realize the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, in an embodiment of the present invention, described all discontinuous Delaunay triangulation network subnet of a band radar image and described all discontinuous Delaunay triangulation network subnet summits to be numbered, comprise: by described all discontinuous Delaunay triangulation network subnet number according to point target order from small to large, be numbered successively, respectively with subnet 1, subnet 2, subnet 3 ... mark, until complete the traversal of all point targets; Record the point target and this point target numbering that contain in each subnet simultaneously.The embodiment of the present invention, by above-mentioned numbering, furthermore achieved that the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, in an embodiment of the present invention, described preset initial step length L adjusts according to the density of the resolution of SAR image and all point targets, but maximum step-length nL is less than or equal to air correlation distance, i.e. 3000m; Described preset initial step length L can be set to 50 times of resolution: as when SAR pixel resolution is 10m, and this preset initial step length L can be set to 500m; When SAR pixel resolution is 20m, this preset initial step length L can be set to 1000m; When resolution is 3m, this preset initial step length L can be set to 150m.The embodiment of the present invention, by the set-up mode of above-mentioned preset initial step length L, by this multi-layer extended mode, can gradually reduce subnet sum further, thus realize the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, in an embodiment of the present invention, described maximum model coefficient of coherence, reflect relative deformation speed and relative altitude error to the fitting degree of consecutive point target, its maximal value is 1, and minimum value is 0; The given threshold value of described maximum model coefficient of coherence is set to 0.7, when the maximum model coefficient of coherence on certain limit is greater than 0.7, thinks that this limit is reliable annexation, two subnets are connected; Otherwise this limit as the fillet of two subnets, can not need the maximum model coefficient of coherence continuing to calculate next limit, until complete the circulation of all limits.The embodiment of the present invention, by the set-up mode of above-mentioned maximum model coefficient of coherence, realizes the quick connection of discontiguous subnetwork further, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, in an embodiment of the present invention, when preset initial step length is set to L, described after having traveled through all subnets, again add up subnet number, if this subnet number only has 1, whole subnet successful connection all is described, then connect without the need to expansion and calculate; If this subnet number is more than 1, the isolated subnet be not connected in addition is described, more successively with 2L, 3L, 4L, 5L ... nL is step-length, repeats above-mentioned steps and carries out expansion connection, until circulation stops, and now all-ones subnet successful connection all; If also there is the subnet of some more than nL of being separated by with major networks, the subnet of described more than nL of being separated by with major networks is abandoned, comprise: after completing the search of the all-ones subnet in L, again add up discontiguous subnetwork number, and described discontiguous subnetwork summit is numbered; If subnet number only has 1, whole subnet successful connection all is described, then connects without the need to expansion and calculate; If subnet number is more than 1, then be step-length with 2L, repeat above-mentioned steps, adjoining subnetworks search is carried out to each subnet; By that analogy, when new subnet number is still more than 1, then expand step-size in search distance further, gradually with 3L, 4L, 5L ... 3000m is step-length cycle calculations, carries out subnet link, until circulation stops, now all-ones subnet has completed expansion connection all; If also there is the subnet of some more than 3000m of being separated by with major networks, the subnet of described more than 3000m of being separated by with major networks is abandoned.The embodiment of the present invention connects the mode calculated by above-mentioned expansion, by this multi-layer extended mode, can gradually reduce subnet sum, thus realize the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
On the other hand, embodiments provide a kind of discontiguous subnetwork coupling arrangement towards time series InSAR, the described discontiguous subnetwork coupling arrangement towards time series InSAR comprises:
Numbered cell, for being numbered all discontinuous Delaunay triangulation network subnet of a band radar image and described all discontinuous Delaunay triangulation network subnet summits;
Search unit, for four frontier points to each subnet in described all discontinuous Delaunay triangulation network subnets, search all point targets within the scope of preset initial step length L and the subnet belonging to described all point targets;
Computing unit, for calculating the distance of described all point targets and current sub network frontier point, and by described apart from ascending sequence, calculate the maximum model coefficient of coherence on every bar limit successively, think the effective edge of this limit for being connected with current sub network by the subnet belonging to this point target when this maximum model coefficient of coherence is greater than given threshold value, otherwise continue to calculate next limit;
Traversal Unit, for after having traveled through all subnets, has again added up subnet number, if this subnet number only has 1, whole subnet successful connection all has been described, then connect without the need to expansion and calculate; If this subnet number is more than 1, the isolated subnet be not connected in addition is described, more successively with 2L, 3L, 4L, 5L ... nL is step-length, repeats above-mentioned steps and carries out expansion connection, until circulation stops, now all-ones subnet has all completed to expand and connected; If also there is the subnet of some more than nL of being separated by with major networks, the subnet of described more than nL of being separated by with major networks is abandoned, wherein n be greater than 1 natural number.
The embodiment of the present invention, by this multi-layer extended mode, can gradually reduce subnet sum, thus realize the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, in an embodiment of the present invention, described numbered cell, be further used for by described all discontinuous Delaunay triangulation network subnet number according to point target order from small to large, be numbered successively, respectively with subnet 1, subnet 2, subnet 3 ... mark, until complete the traversal of all point targets; Record the point target and this point target numbering that contain in each subnet simultaneously.The embodiment of the present invention, by above-mentioned numbering, furthermore achieved that the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, in an embodiment of the present invention, described preset initial step length L adjusts according to the density of the resolution of SAR image and all point targets, but maximum step-length nL is less than or equal to air correlation distance, i.e. 3000m; Described preset initial step length L can be set to 50 times of resolution: as when SAR pixel resolution is 10m, and this preset initial step length L can be set to 500m; When SAR pixel resolution is 20m, this preset initial step length L can be set to 1000m; When resolution is 3m, this preset initial step length L can be set to 150m.The embodiment of the present invention, by the set-up mode of above-mentioned preset initial step length L, by this multi-layer extended mode, can gradually reduce subnet sum, thus realize the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, in an embodiment of the present invention, described maximum model coefficient of coherence, reflect relative deformation speed and relative altitude error to the fitting degree of consecutive point target, its maximal value is 1, and minimum value is 0; The given threshold value of described maximum model coefficient of coherence is set to 0.7, when the maximum model coefficient of coherence on certain limit is greater than 0.7, thinks that this limit is reliable annexation, two subnets are connected; Otherwise this limit as the fillet of two subnets, can not need the maximum model coefficient of coherence continuing to calculate next limit, until complete the circulation of all limits.The embodiment of the present invention, by the set-up mode of above-mentioned maximum model coefficient of coherence, realizes the quick connection of discontiguous subnetwork further, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, in an embodiment of the present invention, when preset initial step length is set to L, described Traversal Unit, is further used for, after completing the search of the all-ones subnet in L, again adding up discontiguous subnetwork number, and being numbered described discontiguous subnetwork summit; If subnet number only has 1, whole subnet successful connection all is described, then connects without the need to expansion and calculate; If subnet number is more than 1, then be step-length with 2L, repeat above-mentioned steps, adjoining subnetworks search is carried out to each subnet; By that analogy, when new subnet number is still more than 1, then expand step-size in search distance further, gradually with 3L, 4L, 5L ... 3000m is step-length cycle calculations, carries out subnet link, until circulation stops, and now all-ones subnet successful connection all; If also there is the subnet of some more than 3000m of being separated by with major networks, the subnet of described more than 3000m of being separated by with major networks is abandoned.The embodiment of the present invention connects the mode calculated by above-mentioned expansion, by this multi-layer extended mode, can gradually reduce subnet sum, thus realize the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Technique scheme has following beneficial effect: the discontiguous subnetwork method of attachment towards time series InSAR described in the embodiment of the present invention comprises: all discontinuous Delaunay triangulation network subnet of a band radar image and described all discontinuous Delaunay triangulation network subnet summits are numbered; To four frontier points of each subnet in described all discontinuous Delaunay triangulation network subnets, search all point targets within the scope of preset initial step length L and the subnet belonging to described all point targets; Calculate the distance of described all point targets and current sub network frontier point, and by described apart from ascending sequence, calculate the maximum model coefficient of coherence on every bar limit successively, think the effective edge of this limit for being connected with current sub network by the subnet belonging to this point target when this maximum model coefficient of coherence is greater than given threshold value, otherwise continue to calculate next limit; After having traveled through all subnets, again add up subnet number, if this subnet number only has 1, whole subnet successful connection all has been described, then connected without the need to expansion and calculate; If this subnet number is more than 1, the isolated subnet be not connected in addition is described, more successively with 2L, 3L, 4L, 5L ... nL is step-length, repeats above-mentioned steps and carries out expansion connection, until circulation stops, and now all-ones subnet successful connection all; If also there is the subnet of some more than nL of being separated by with major networks, the subnet of described more than nL of being separated by with major networks is abandoned, wherein n be greater than 1 natural number.The embodiment of the present invention reaches following technique effect: by this multi-layer extended mode, can gradually reduce subnet sum, thus realizes the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is a kind of discontiguous subnetwork method of attachment process flow diagram towards time series InSAR of the embodiment of the present invention;
Fig. 2 is a kind of discontiguous subnetwork coupling arrangement structural representation towards time series InSAR of the embodiment of the present invention;
Fig. 3 is the average amplitude figure of test site, the Hebei province 24 scape European Space Agency ENVISAT ASAR image that application example of the present invention uses:
Fig. 4 is a kind of discontiguous subnetwork method of attachment process flow diagram towards time series InSAR of application example of the present invention;
Fig. 5 is the test site discontinuous Delaunay triangulation network subnet schematic diagram that application example of the present invention uses;
Fig. 6 be application example of the present invention with 500 meters for step-length, maximum search distance is 3000 meters, utilizes the network diagram that multi-layer quick connecting method obtains;
Fig. 7 is the linear earth's surface rate of deformation schematic diagram (unit: millimeter/year) of 5179 point targets that application example of the present invention obtains based on multi-layer quick connecting method;
Fig. 8 is the network diagram utilizing the method for attachment of prior art complex network to obtain;
Fig. 9 is the linear earth's surface rate of deformation schematic diagram (unit: millimeter/year) of 5141 point targets that prior art obtains based on complex network method of attachment.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.
As shown in Figure 1, be a kind of discontiguous subnetwork method of attachment process flow diagram towards time series InSAR of the embodiment of the present invention, the described discontiguous subnetwork method of attachment towards time series InSAR comprises:
101, all discontinuous Delaunay triangulation network subnet of a band radar image and described all discontinuous Delaunay triangulation network subnet summits are numbered;
102, to four frontier points of each subnet in described all discontinuous Delaunay triangulation network subnets, all point targets within the scope of preset initial step length L and the subnet belonging to described all point targets is searched;
103, the distance of described all point targets and current sub network frontier point is calculated, and by described apart from ascending sequence, calculate the maximum model coefficient of coherence on every bar limit successively, think the effective edge of this limit for being connected with current sub network by the subnet belonging to this point target when this maximum model coefficient of coherence is greater than given threshold value, otherwise continue to calculate next limit;
104, after having traveled through all subnets, again add up subnet number, if this subnet number only has 1, whole subnet successful connection all is described, then connect without the need to expansion and calculate; If this subnet number is more than 1, the isolated subnet be not connected in addition is described, more successively with 2L, 3L, 4L, 5L ... nL is step-length, repeats above-mentioned steps and carries out expansion connection, until circulation stops, and now all-ones subnet successful connection all; If also there is the subnet of some more than nL of being separated by with major networks, the subnet of described more than nL of being separated by with major networks is abandoned, wherein n be greater than 1 natural number.
The embodiment of the present invention, by this multi-layer extended mode, can gradually reduce subnet sum, thus realize the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, described all discontinuous Delaunay triangulation network subnet of a band radar image and described all discontinuous Delaunay triangulation network subnet summits to be numbered, comprise: by described all discontinuous Delaunay triangulation network subnet number according to point target order from small to large, be numbered successively, respectively with subnet 1, subnet 2, subnet 3 ... mark, until complete the traversal of all point targets; Record the point target and this point target numbering that contain in each subnet simultaneously.The embodiment of the present invention, by above-mentioned numbering, furthermore achieved that the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, described preset initial step length L adjusts according to the density of the resolution of SAR image and all point targets, but maximum step-length nL is less than or equal to air correlation distance, i.e. 3000m; Described preset initial step length L can be set to 50 times of resolution: as when SAR pixel resolution is 10m, and this preset initial step length L can be set to 500m; When SAR pixel resolution is 20m, this preset initial step length L can be set to 1000m; When resolution is 3m, this preset initial step length L can be set to 150m.The embodiment of the present invention, by the set-up mode of above-mentioned preset initial step length L, by this multi-layer extended mode, can gradually reduce subnet sum, thus realize the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, described maximum model coefficient of coherence, reflect relative deformation speed and relative altitude error to the fitting degree of consecutive point target, its maximal value is 1, and minimum value is 0; The given threshold value of described maximum model coefficient of coherence is set to 0.7, when the maximum model coefficient of coherence on certain limit is greater than 0.7, thinks that this limit is reliable annexation, two subnets are connected; Otherwise this limit as the fillet of two subnets, can not need the maximum model coefficient of coherence continuing to calculate next limit, until complete the circulation of all limits.The embodiment of the present invention, by the set-up mode of above-mentioned maximum model coefficient of coherence, realizes the quick connection of discontiguous subnetwork further, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, when preset initial step length is set to L, described after having traveled through all subnets, again add up subnet number, if this subnet number only has 1, whole subnet successful connection all is described, then connect without the need to expansion and calculate; If this subnet number is more than 1, the isolated subnet be not connected in addition is described, more successively with 2L, 3L, 4L, 5L ... nL is step-length, repeats above-mentioned steps and carries out expansion connection, until circulation stops, and now all-ones subnet successful connection all; If also there is the subnet of some more than nL of being separated by with major networks, the subnet of described more than nL of being separated by with major networks is abandoned, comprise: after completing the search of the all-ones subnet in L, again add up discontiguous subnetwork number, and described discontiguous subnetwork summit is numbered; If subnet number only has 1, whole subnet successful connection all is described, then connects without the need to expansion and calculate; If subnet number is more than 1, then be step-length with 2L, repeat above-mentioned steps, adjoining subnetworks search is carried out to each subnet; By that analogy, when new subnet number is still more than 1, then expand step-size in search distance further, gradually with 3L, 4L, 5L ... 3000m is step-length cycle calculations, carries out subnet link, until circulation stops, and now all-ones subnet successful connection all; If also there is the subnet of some more than 3000m of being separated by with major networks, the subnet of described more than 3000m of being separated by with major networks is abandoned.The embodiment of the present invention connects the mode calculated by above-mentioned expansion, by this multi-layer extended mode, can gradually reduce subnet sum, thus realize the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
As shown in Figure 2, be a kind of discontiguous subnetwork coupling arrangement structural representation towards time series InSAR of the embodiment of the present invention, the described discontiguous subnetwork coupling arrangement towards time series InSAR comprises:
Numbered cell 21, for being numbered all discontinuous Delaunay triangulation network subnet of a band radar image and described all discontinuous Delaunay triangulation network subnet summits;
Search unit 22, for four frontier points to each subnet in described all discontinuous Delaunay triangulation network subnets, search all point targets within the scope of preset initial step length L and the subnet belonging to described all point targets;
Computing unit 23, for calculating the distance of described all point targets and current sub network frontier point, and by described apart from ascending sequence, calculate the maximum model coefficient of coherence on every bar limit successively, think the effective edge of this limit for being connected with current sub network by the subnet belonging to this point target when this maximum model coefficient of coherence is greater than given threshold value, otherwise continue to calculate next limit;
Traversal Unit 24, for after having traveled through all subnets, has again added up subnet number, if this subnet number only has 1, whole subnet successful connection all has been described, then connect without the need to expansion and calculate; If this subnet number is more than 1, the isolated subnet be not connected in addition is described, more successively with 2L, 3L, 4L, 5L ... nL is step-length, repeats above-mentioned steps and carries out expansion connection, until circulation stops, and now all-ones subnet successful connection all; If also there is the subnet of some more than nL of being separated by with major networks, the subnet of described more than nL of being separated by with major networks is abandoned, wherein n be greater than 1 natural number.
The embodiment of the present invention, by this multi-layer extended mode, can gradually reduce subnet sum, thus realize the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, described numbered cell 21, be further used for by described all discontinuous Delaunay triangulation network subnet number according to point target order from small to large, be numbered successively, respectively with subnet 1, subnet 2, subnet 3 ... mark, until complete the traversal of all point targets; Record the point target and this point target numbering that contain in each subnet simultaneously.The embodiment of the present invention, by above-mentioned numbering, furthermore achieved that the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, described preset initial step length L adjusts according to the density of the resolution of SAR image and all point targets, but maximum step-length nL is less than or equal to air correlation distance, i.e. 3000m; Described preset initial step length L can be set to 50 times of resolution: as when SAR pixel resolution is 10m, and this preset initial step length L can be set to 500m; When SAR pixel resolution is 20m, this preset initial step length L can be set to 1000m; When resolution is 3m, this preset initial step length L can be set to 150m.Inventive embodiments, by the set-up mode of above-mentioned preset initial step length L, by this multi-layer extended mode, can gradually reduce subnet sum, thus realize the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, described maximum model coefficient of coherence, reflect relative deformation speed and relative altitude error to the fitting degree of consecutive point target, its maximal value is 1, and minimum value is 0; The given threshold value of described maximum model coefficient of coherence is set to 0.7, when the maximum model coefficient of coherence on certain limit is greater than 0.7, thinks that this limit is reliable annexation, two subnets are connected; Otherwise this limit as the fillet of two subnets, can not need the maximum model coefficient of coherence continuing to calculate next limit, until complete the circulation of all limits.The embodiment of the present invention, by the set-up mode of above-mentioned maximum model coefficient of coherence, realizes the quick connection of discontiguous subnetwork further, is particularly conducive to the Ground Deformation monitoring of large regions.
Preferably, when preset initial step length is set to L, described Traversal Unit 24, is further used for, after completing the search of the all-ones subnet in L, again adding up discontiguous subnetwork number, and being numbered described discontiguous subnetwork summit; If subnet number only has 1, whole subnet successful connection all is described, then connects without the need to expansion and calculate; If subnet number is more than 1, then be step-length with 2L, repeat above-mentioned steps, adjoining subnetworks search is carried out to each subnet; By that analogy, when new subnet number is still more than 1, then expand step-size in search distance further, gradually with 3L, 4L, 5L ... 3000m is step-length cycle calculations, carries out subnet link, until circulation stops, and now all-ones subnet successful connection all; If also there is the subnet of some more than 3000m of being separated by with major networks, the subnet of described more than 3000m of being separated by with major networks is abandoned.The embodiment of the present invention connects the mode calculated by above-mentioned expansion, by this multi-layer extended mode, can gradually reduce subnet sum, thus realize the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Technique scheme has following beneficial effect: the discontiguous subnetwork method of attachment towards time series InSAR described in the embodiment of the present invention comprises: all discontinuous Delaunay triangulation network subnet of a band radar image and described all discontinuous Delaunay triangulation network subnet summits are numbered; To four frontier points of each subnet in described all discontinuous Delaunay triangulation network subnets, search all point targets within the scope of preset initial step length L and the subnet belonging to described all point targets; Calculate the distance of described all point targets and current sub network frontier point, and by described apart from ascending sequence, calculate the maximum model coefficient of coherence on every bar limit successively, think the effective edge of this limit for being connected with current sub network by the subnet belonging to this point target when this maximum model coefficient of coherence is greater than given threshold value, otherwise continue to calculate next limit; After having traveled through all subnets, again add up subnet number, if this subnet number only has 1, whole subnet successful connection all has been described, then connected without the need to expansion and calculate; If this subnet number is more than 1, the isolated subnet be not connected in addition is described, more successively with 2L, 3L, 4L, 5L ... nL is step-length, repeats above-mentioned steps and carries out expansion connection, until circulation stops, and now all-ones subnet successful connection all; If also there is the subnet of some more than nL of being separated by with major networks, the subnet of described more than nL of being separated by with major networks is abandoned, wherein n be greater than 1 natural number.The embodiment of the present invention reaches following technique effect: by this multi-layer extended mode, can gradually reduce subnet sum, thus realizes the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Below lift application example to be described:
As shown in Figure 3, be the average amplitude figure of test site, the Hebei province 24 scape European Space Agency ENVISAT ASAR image that application example of the present invention uses, resolution is about 20 meters, and line number is 740 pixels, and columns is 600 pixels, is below described in detail:
As shown in Figure 4, be a kind of discontiguous subnetwork method of attachment process flow diagram towards time series InSAR of application example of the present invention, described method comprises:
401, all discontinuous Delaunay triangulation network subnet of a band radar image and described all discontinuous Delaunay triangulation network subnet summits are numbered.
First the point target of all Delaunay triangulation network Subnetwork connections is numbered, arranges according to the order of Column Row, can obtain a little 1, point 2, point 3 ..., put N number of point target;
Secondly according to point target number order from small to large, search for the subnet be connected with this point target successively, to its numbering, identify, until complete the traversal of all point targets with subnet 1, subnet 2, subnet 3 ..., subnet M respectively;
In the process to subnet search numbering, record point target and numbering thereof that this subnet contains simultaneously.
402, to four frontier points of each subnet in described all discontinuous Delaunay triangulation network subnets, search at the subnet with 500m belonging to all point targets in radius and described all point targets.
Cycle calculations from subnet 1, for current sub network, finds the frontier point of this subnet, i.e. the intersection point of maximum boundary rectangle and subnet with maximum boundary rectangle, under normal circumstances, can obtain 4 frontier points; But when this subnet only comprises two or three summits, then the frontier point number of this subnet will be less than 4.
Then, centered by each frontier point, take 500m as radius, in the space search window of 1000m*1000m, search the point target of non-present subnet, and record the subnet number belonging to each point target.
Optionally, search radius suitably can adjust according to different resolution SAR image here, increases or reduces.
403, the distance of described all point targets and current sub network frontier point is calculated, and by described apart from ascending sequence, calculate the maximum model coefficient of coherence on every bar limit successively, think the effective edge of this limit for being connected with current sub network by the subnet belonging to this point target when this maximum model coefficient of coherence is greater than given threshold value, otherwise continue to calculate next limit.
By the coordinate of each point target, calculate the point target and current border point distance that search, and by ascending sequence, calculate the maximum model coefficient of coherence on every bar limit successively.Model coefficient of coherence circular is as follows:
In formula, γ is model coefficient of coherence; M is that interference image is to number; δ φ
diffor the second order difference phase place of consecutive point target; δ φ
modelfor the linear model phase place of adjacent high coherent point; (x
m, y
m), (x
n, y
n) be the position coordinates of two adjacent high coherent point; T
ibe the time basis of the i-th width differential phase figure, T
i=t
m-t
s.
In formula (1), consecutive point target (x
m, y
m), (x
n, y
n) between second order difference phase place can be expressed as:
Wherein, λ is radar carrier wavelength, and it is long that r is that target arrives radar sensor oblique distance, and b is vertical parallax, and θ is radar incident angle, and ξ is vertical error, v
1be respectively Linear deformation rate, t
m, t
srepresent the time of the main and auxiliary image capturing that interference image is right.(x
m, y
m), (x
n, y
n) be the position coordinates on this summit, limit two, T
ibe the time basis of the i-th width interferogram, T
i=t
m-t
s, β is permanent set phase place, and α is atmospheric effect phase place, and n is noise phase.
In formula (1), the linear model phase place between consecutive point target can be expressed as:
Wherein, Δ v
1with Δ ξ is opposite linear rate of deformation between two high coherent point and relative altitude error respectively.
Formula (1) obtains maximum model coefficient of coherence by optimum solving method, and it reflects relative deformation speed and relative altitude error to the fitting degree of consecutive point target, its maximal value is 1, and minimum value is 0.Usually the threshold value of this coefficient is set to 0.7, and when the maximum model coefficient of coherence on certain limit is greater than 0.7, think that this limit is reliable annexation, two subnets can be connected; Otherwise this limit as the fillet of two subnets, can not need the maximum model coefficient of coherence continuing to calculate next limit, until complete the circulation of all limits.
404, after having traveled through all subnets, again add up subnet number, if this subnet number only has 1, whole subnet successful connection all is described, then connect without the need to expansion and calculate; If this subnet number is more than 1, the isolated subnet be not connected in addition is described, more successively with 1000m, 1500m, 2000m, 2500m, 3000m for step-length, repeat above-mentioned steps carry out expansion connect, until circulation stopping, now all-ones subnet successful connection all; If also there is the subnet of some more than 3000m of being separated by with major networks, the subnet of described more than 3000m of being separated by with major networks is abandoned.
After completing the search of the all-ones subnet in 500m (or other distances), again add up discontiguous subnetwork number, and its summit is numbered.When new subnet number is greater than 1, then be step-length with 1000m, repeat step 103 and step 104, adjoining subnetworks search is carried out to each subnet; By that analogy, then expand step-size in search distance further, gradually with 1500m, 2000m, 2500m, 3000m for step-length, carry out subnet expansion and connect, until circulation stops, in cyclic process, if when new subnet number is 1, then stop expansion connecting.
If when final subnet number only has 1, all-ones subnet successful connection all is described; Otherwise the subnet that also there are some more than 3000m of being separated by with major networks is described, these subnets will no longer be considered.
Optionally, step-size in search can adjust according to the density of the resolution of SAR image and all point targets, but maximum step-length should not exceed air correlation distance, i.e. 3000m.
This application example said method has following beneficial effect: by multi-layer extended mode, can gradually reduce subnet number, thus realizes the quick connection of discontiguous subnetwork.
In order to validity and the superiority of technical solution of the present invention are described better, now following comparative analysis is carried out to application example of the present invention and existing complex network method of attachment: see Fig. 3, the average amplitude figure of test site, the Hebei province 24 scape European Space Agency ENVISAT ASAR image that application example of the present invention uses, resolution is about 20 meters, line number is 740 pixels, and columns is 600 pixels.As shown in Figure 5, be the discontinuous Delaunay triangulation network subnet in test site that application example of the present invention uses, have 349 subnets, 5219 point targets.As shown in Figure 6, be application example of the present invention with 500 meters for step-length, maximum search distance is 3000 meters, utilizes the network that multi-layer quick connecting method obtains.Connected by subnet expansion, 14 minutes consuming time, finally obtain 9 subnets, 10041 limits.As shown in Figure 7, the linear earth's surface rate of deformation (unit: millimeter/year) of 5179 point targets obtained based on multi-layer quick connecting method for application example of the present invention.Fig. 8 is the network utilizing the method for attachment of prior art complex network to obtain.Within the scope of 3000 meters, connected by complex network, 43 minutes consuming time, finally obtain 9 subnets, 13000 limits.Fig. 9 is the linear earth's surface rate of deformation (unit: millimeter/year) of 5141 point targets that prior art obtains based on complex network method of attachment.Can see, the subnet number that two kinds of methods finally connect is the same, be all 9 subnets, but the efficiency of multi-layer quick connecting method is higher, is consuming timely only 14 minutes, and the point target number integrated is more, is 5179; And complex network method of attachment is consuming time much larger than multi-layer quick connecting method, be 43 minutes, the point target number integrated is also less, is 5141.This shows, for discontinuous Delaunay triangulation network subnet, multi-layer quick connecting method of the present invention is better than the method for attachment of prior art complex network.
As the interference baseline parameter that following table 1 is the 24 scape European Space Agency ENVISAT ASAR images that application example of the present invention adopts.
Table 1
If following table 2 is the multi-layer quick connecting method of application example of the present invention employing and existing complex network method of attachment performance comparison.
| Relatively content | Multi-layer Subnetwork connection | Complex network connects |
| The limit number connected | 10041 | 13000 |
| Final subnet number | 9 | 9 |
| Integrated counts | 5179 | 5141 |
| (point) consuming time | 14 | 43 |
Table 2
Can see that complex network is due within the scope of 3000 meters by above-mentioned table 1 and table 2, whole point targets of search non-present subnet, the limit number that it obtains far exceedes the limit number that multi-layer Subnetwork connection obtains, and corresponding calculated amount is larger, obviously consuming time, take 43 minutes altogether; And multi-layer Subnetwork connection method, by progressively outside expanded search subnet, subnet number can be reduced fast, calculate consuming time less, take 14 minutes, be only complex network consuming time 32.56%.
As fully visible, application example multi-layer Subnetwork connection method of the present invention not only will more than complex network at finally integrated point target number, the performance of its aspect consuming time is also far superior to complex network, therefore the method is more suitable for the quick connection of discontiguous subnetwork, is particularly conducive to the Ground Deformation monitoring of large regions.
Those skilled in the art can also recognize the various illustrative components, blocks (illustrative logical block) that the embodiment of the present invention is listed, unit, and step can pass through electronic hardware, computer software, or both combinations realize.For the replaceability (interchangeability) of clear displaying hardware and software, above-mentioned various illustrative components (illustrative components), unit and step have universally described their function.Such function is the designing requirement realizing depending on specific application and whole system by hardware or software.Those skilled in the art for often kind of specifically application, can use the function described in the realization of various method, but this realization can should not be understood to the scope exceeding embodiment of the present invention protection.
Various illustrative logical block described in the embodiment of the present invention, or unit can pass through general processor, digital signal processor, special IC (ASIC), field programmable gate array or other programmable logic device, discrete gate or transistor logic, discrete hardware components, or the design of above-mentioned any combination realizes or operates described function.General processor can be microprocessor, and alternatively, this general processor also can be any traditional processor, controller, microcontroller or state machine.Processor also can be realized by the combination of calculation element, such as digital signal processor and microprocessor, multi-microprocessor, and a Digital Signal Processor Core combined by one or more microprocessor, or other similar configuration any realizes.
The software module that method described in the embodiment of the present invention or the step of algorithm directly can embed hardware, processor performs or the combination of both.Software module can be stored in the storage medium of other arbitrary form in RAM storer, flash memory, ROM storer, eprom memory, eeprom memory, register, hard disk, moveable magnetic disc, CD-ROM or this area.Exemplarily, storage medium can be connected with processor, with make processor can from storage medium reading information, and write information can be deposited to storage medium.Alternatively, storage medium can also be integrated in processor.Processor and storage medium can be arranged in ASIC, and ASIC can be arranged in user terminal.Alternatively, processor and storage medium also can be arranged in the different parts in user terminal.
In one or more exemplary design, the above-mentioned functions described by the embodiment of the present invention can realize in the combination in any of hardware, software, firmware or this three.If realized in software, these functions can store on the medium with computer-readable, or are transmitted on the medium of computer-readable with one or more instruction or code form.Computer readable medium comprises computer storage medium and is convenient to make to allow computer program transfer to the telecommunication media in other place from a place.Storage medium can be that any general or special computer can the useable medium of access.Such as, such computer readable media can include but not limited to RAM, ROM, EEPROM, CD-ROM or other optical disc storage, disk storage or other magnetic storage device, or other anyly may be used for carrying or store the medium that can be read the program code of form with instruction or data structure and other by general or special computer or general or special processor.In addition, any connection can be properly termed computer readable medium, such as, if software is by a concentric cable, fiber optic cables, twisted-pair feeder, Digital Subscriber Line (DSL) or being also comprised in defined computer readable medium with wireless way for transmittings such as such as infrared, wireless and microwaves from a web-site, server or other remote resource.Described video disc (disk) and disk (disc) comprise Zip disk, radium-shine dish, CD, DVD, floppy disk and Blu-ray Disc, and disk is usually with magnetic duplication data, and video disc carries out optical reproduction data with laser usually.Above-mentioned combination also can be included in computer readable medium.
Above-described embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only the specific embodiment of the present invention; the protection domain be not intended to limit the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (10)
1. towards a discontiguous subnetwork method of attachment of time series InSAR, it is characterized in that, the described discontiguous subnetwork method of attachment towards time series InSAR comprises:
The all discontinuous Delaunay triangulation network subnet of one band radar image and described all discontinuous Delaunay triangulation network subnet summits are numbered;
To four frontier points of each subnet in described all discontinuous Delaunay triangulation network subnets, search all point targets within the scope of preset initial step length L and the subnet belonging to described all point targets;
Calculate the distance of described all point targets and current sub network frontier point, and by described apart from ascending sequence, calculate the maximum model coefficient of coherence on every bar limit successively, think the effective edge of this limit for being connected with current sub network by the subnet belonging to this point target when this maximum model coefficient of coherence is greater than given threshold value, otherwise continue to calculate next limit;
After having traveled through all subnets, again add up subnet number, if this subnet number only has 1, whole subnet successful connection all has been described, then connected without the need to expansion and calculate; If this subnet number is more than 1, the isolated subnet be not connected in addition is described, more successively with 2L, 3L, 4L, 5L ... nL is step-length, repeats above-mentioned steps and carries out expansion connection, until circulation stops, and now all-ones subnet successful connection all; If also there is the subnet of some more than nL of being separated by with major networks, the subnet of described more than nL of being separated by with major networks is abandoned, wherein n be greater than 1 natural number.
2. as claimed in claim 1 towards the discontiguous subnetwork method of attachment of time series InSAR, it is characterized in that, described all discontinuous Delaunay triangulation network subnet of a band radar image and described all discontinuous Delaunay triangulation network subnet summits to be numbered, to comprise:
By described all discontinuous Delaunay triangulation network subnet number according to point target order from small to large, be numbered successively, identify, until complete the traversal of all point targets with subnet 1, subnet 2, subnet 3 ... subnet M respectively; Record the point target that contains in each subnet and this point target numbering simultaneously, wherein M be greater than 1 natural number.
3. as claimed in claim 1 towards the discontiguous subnetwork method of attachment of time series InSAR, it is characterized in that, described preset initial step length L adjusts according to the density of the resolution of SAR image and all point targets, but maximum step-length nL is less than or equal to air correlation distance, i.e. 3000m; Described preset initial step length L is set to 50 times of image resolution: when SAR pixel resolution is 10m, and this preset initial step length L is set to 500m; When SAR pixel resolution is 20m, this preset initial step length L is set to 1000m; When resolution is 3m, this preset initial step length L is set to 150m.
4. as claimed in claim 1 towards the discontiguous subnetwork method of attachment of time series InSAR, it is characterized in that, described maximum model coefficient of coherence, reflect relative deformation speed and relative altitude error to the fitting degree of consecutive point target, its maximal value is 1, and minimum value is 0; The given threshold value of described maximum model coefficient of coherence is set to 0.7, when the maximum model coefficient of coherence on certain limit is greater than 0.7, thinks that this limit is reliable annexation, two subnets are connected; Otherwise this limit as the fillet of two subnets, can not need the maximum model coefficient of coherence continuing to calculate next limit, until complete the circulation of all limits.
5. as claimed in claim 1 towards the discontiguous subnetwork method of attachment of time series InSAR, it is characterized in that, when preset initial step length is set to L, described after having traveled through all subnets, again subnet number is added up, if this subnet number only has 1, whole subnet successful connection all is described, then connects without the need to expansion and calculate; If this subnet number is more than 1, the isolated subnet be not connected in addition is described, more successively with 2L, 3L, 4L, 5L ... nL is step-length, repeats above-mentioned steps and carries out expansion connection, until circulation stops, now all-ones subnet has all completed to expand and connected; If also there is the subnet of some more than nL of being separated by with major networks, the subnet of described more than nL of being separated by with major networks is abandoned, comprising:
After completing the search of the all-ones subnet in L, again add up discontiguous subnetwork number, and described discontiguous subnetwork summit is numbered; If subnet number only has 1, whole subnet successful connection all is described, then connects without the need to expansion and calculate; If subnet number is more than 1, then be step-length with 2L, repeat above-mentioned steps, adjoining subnetworks search is carried out to each subnet; By that analogy, when new subnet number is still more than 1, then expand step-size in search distance further, gradually with 3L, 4L, 5L ... 3000m is step-length cycle calculations, carries out subnet link, until circulation stops, now all-ones subnet has completed expansion connection all; If also there is the subnet of some more than 3000m of being separated by with major networks, the subnet of described more than 3000m of being separated by with major networks is abandoned.
6. towards a discontiguous subnetwork coupling arrangement of time series InSAR, it is characterized in that, the described discontiguous subnetwork coupling arrangement towards time series InSAR comprises:
Numbered cell, for being numbered all discontinuous Delaunay triangulation network subnet of a band radar image and described all discontinuous Delaunay triangulation network subnet summits;
Search unit, for four frontier points to each subnet in described all discontinuous Delaunay triangulation network subnets, search all point targets within the scope of preset initial step length L and the subnet belonging to described all point targets;
Computing unit, for calculating the distance of described all point targets and current sub network frontier point, and by described apart from ascending sequence, calculate the maximum model coefficient of coherence on every bar limit successively, think the effective edge of this limit for being connected with current sub network by the subnet belonging to this point target when this maximum model coefficient of coherence is greater than given threshold value, otherwise continue to calculate next limit;
Traversal Unit, for after having traveled through all subnets, has again added up subnet number, if this subnet number only has 1, whole subnet successful connection all has been described, then connect without the need to expansion and calculate; If this subnet number is more than 1, the isolated subnet be not connected in addition is described, more successively with 2L, 3L, 4L, 5L ... nL is step-length, repeats above-mentioned steps and carries out expansion connection, until circulation stops, now all-ones subnet has all completed to expand and connected; If also there is the subnet of some more than nL of being separated by with major networks, the subnet of described more than nL of being separated by with major networks is abandoned, wherein n be greater than 1 natural number.
7. as claimed in claim 6 towards the discontiguous subnetwork coupling arrangement of time series InSAR, it is characterized in that, described numbered cell, be further used for by described all discontinuous Delaunay triangulation network subnet number according to point target order from small to large, be numbered successively, identify, until complete the traversal of all point targets with subnet 1, subnet 2, subnet 3 ... subnet M respectively; Record the point target that contains in each subnet and this point target numbering simultaneously, wherein M be greater than 1 natural number.
8. as claimed in claim 6 towards the discontiguous subnetwork coupling arrangement of time series InSAR, it is characterized in that, described preset initial step length L adjusts according to the density of the resolution of SAR image and all point targets, but maximum step-length nL is less than or equal to air correlation distance, i.e. 3000m; Described preset initial step length L is set to 50 times of image resolution: when SAR pixel resolution is 10m, and this preset initial step length L is set to 500m; When SAR pixel resolution is 20m, this preset initial step length L is set to 1000m; When resolution is 3m, this preset initial step length L is set to 150m.
9. as claimed in claim 6 towards the discontiguous subnetwork coupling arrangement of time series InSAR, it is characterized in that, described maximum model coefficient of coherence, reflect relative deformation speed and relative altitude error to the fitting degree of consecutive point target, its maximal value is 1, and minimum value is 0; The given threshold value of described maximum model coefficient of coherence is set to 0.7, when the maximum model coefficient of coherence on certain limit is greater than 0.7, thinks that this limit is reliable annexation, two subnets are connected; Otherwise this limit as the fillet of two subnets, can not need the maximum model coefficient of coherence continuing to calculate next limit, until complete the circulation of all limits.
10. as claimed in claim 6 towards the discontiguous subnetwork coupling arrangement of time series InSAR, it is characterized in that, when preset initial step length is set to L, described Traversal Unit, be further used for after completing the search of the all-ones subnet in L, again add up discontiguous subnetwork number, and described discontiguous subnetwork summit is numbered; If subnet number only has 1, whole subnet successful connection all is described, then connects without the need to expansion and calculate; If subnet number is more than 1, then be step-length with 2L, repeat above-mentioned steps, adjoining subnetworks search is carried out to each subnet; By that analogy, when new subnet number is still more than 1, then expand step-size in search distance further, gradually with 3L, 4L, 5L ... 3000m is step-length cycle calculations, carries out subnet link, until circulation stops, now all-ones subnet has completed expansion connection all; If also there is the subnet of some more than 3000m of being separated by with major networks, the subnet of described more than 3000m of being separated by with major networks is abandoned.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410410540.7A CN104239419B (en) | 2014-08-20 | 2014-08-20 | A kind of discontiguous subnetwork connection method and device towards time series InSAR |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410410540.7A CN104239419B (en) | 2014-08-20 | 2014-08-20 | A kind of discontiguous subnetwork connection method and device towards time series InSAR |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104239419A true CN104239419A (en) | 2014-12-24 |
| CN104239419B CN104239419B (en) | 2017-11-28 |
Family
ID=52227479
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410410540.7A Expired - Fee Related CN104239419B (en) | 2014-08-20 | 2014-08-20 | A kind of discontiguous subnetwork connection method and device towards time series InSAR |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104239419B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105842691A (en) * | 2016-03-17 | 2016-08-10 | 中国测绘科学研究院 | Automatic adjustment and inlay method of large area InSAR deformation result and device thereof |
| CN110673145A (en) * | 2019-10-24 | 2020-01-10 | 中国地质大学（北京） | InSAR (interferometric synthetic Aperture Radar) surface deformation monitoring method and system based on discontinuous coherence |
| CN111308468A (en) * | 2019-11-27 | 2020-06-19 | 北京东方至远科技股份有限公司 | Method for automatically identifying deformation risk area based on In SAR technology |
| CN111859786A (en) * | 2020-07-03 | 2020-10-30 | 安徽理工大学 | Full-scale gradient mining subsidence D-InSAR three-dimensional prediction method for improving dynamic prediction model constraint |
| CN113008202A (en) * | 2021-03-31 | 2021-06-22 | 中铁上海设计院集团有限公司 | Ground settlement monitoring method integrating different synthetic aperture radar interferometry |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102629390A (en) * | 2012-02-23 | 2012-08-08 | 中国测绘科学研究院 | Mass airborne LiDAR point cloud Delaunay triangulation network parallel construction method and apparatus thereof |
-
2014
- 2014-08-20 CN CN201410410540.7A patent/CN104239419B/en not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102629390A (en) * | 2012-02-23 | 2012-08-08 | 中国测绘科学研究院 | Mass airborne LiDAR point cloud Delaunay triangulation network parallel construction method and apparatus thereof |
Non-Patent Citations (2)
| Title |
|---|
| ，O MORA，JJ MALLORQUI，A BROQUETAS: "Linear and nonlinear terrain deformation maps from a reduced set of interferometric SAR images", 《IEEE TRANSACTIONS ON GEOSCIENCE & REMOTE SENSING》 * |
| W XU，I CUMMING: "A region-growing algorithm for InSAR phase unwrapping", 《IEEE TRANSACTIONS ON GEOSCIENCE & REMOTE SENSING》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105842691A (en) * | 2016-03-17 | 2016-08-10 | 中国测绘科学研究院 | Automatic adjustment and inlay method of large area InSAR deformation result and device thereof |
| CN105842691B (en) * | 2016-03-17 | 2018-01-12 | 中国测绘科学研究院 | The automation adjustment and method for embedding and device of big region InSAR deformation results |
| CN110673145A (en) * | 2019-10-24 | 2020-01-10 | 中国地质大学（北京） | InSAR (interferometric synthetic Aperture Radar) surface deformation monitoring method and system based on discontinuous coherence |
| CN111308468A (en) * | 2019-11-27 | 2020-06-19 | 北京东方至远科技股份有限公司 | Method for automatically identifying deformation risk area based on In SAR technology |
| CN111859786A (en) * | 2020-07-03 | 2020-10-30 | 安徽理工大学 | Full-scale gradient mining subsidence D-InSAR three-dimensional prediction method for improving dynamic prediction model constraint |
| CN113008202A (en) * | 2021-03-31 | 2021-06-22 | 中铁上海设计院集团有限公司 | Ground settlement monitoring method integrating different synthetic aperture radar interferometry |
| CN113008202B (en) * | 2021-03-31 | 2023-02-21 | 中铁上海设计院集团有限公司 | Ground settlement monitoring method integrating different synthetic aperture radar interferometry |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104239419B (en) | 2017-11-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102608584B (en) | Time sequence InSAR (Interferometric Synthetic Aperture Radar) deformation monitoring method and device based on polynomial inversion model | |
| CN104239419A (en) | Discontinuous subnet connecting method and device oriented to time sequence InSAR | |
| CN101976468B (en) | Method and system for visualizing multiresolution dynamic landform | |
| CN106600661B (en) | The method for accurately generating segmental arc geologic section | |
| Liu et al. | Urban subsidence monitoring by SBAS-InSAR technique with multi-platform SAR images: A case study of Beijing Plain, China | |
| CN102636809B (en) | Method for generating spreading angle domain common image point gathers | |
| CN105355042A (en) | Road network extraction method based on taxi GPS | |
| CN105759311A (en) | Near-real time earthquake source position positioning method | |
| CN108053477A (en) | The Numerical Methods of deformation in a kind of pipeline | |
| CN107884772A (en) | Method and device of the optimal interference images of time series InSAR to selection | |
| CN103713287A (en) | Elevation reestablishing method and device based on coprime of multiple base lines | |
| WO2023098441A1 (en) | Method and apparatus for inverting crustal structure of passive continental margin on basis of subsidence in stratigraphic record | |
| CN103824510A (en) | Voronoi diagram-based electronic map point element screening method | |
| CN107484139A (en) | A kind of car networking Cooperative Localization Method and device based on geographical location information | |
| CN108845358B (en) | Tomography and the recognition methods of structural anomaly body and device | |
| CN104375172A (en) | Volcanic underlayer structural configuration correct imaging method | |
| CN105425286A (en) | Earthquake time-travelling acquisition method and crosshole earthquake time-travelling tomography method based on the earthquake time-travelling acquisition method | |
| CN103698810A (en) | Hybrid network minimum travel time ray tracing tomography method | |
| CN112099082B (en) | Seismic folding wave travel time inversion method for coplanar element common azimuth gather | |
| CN103592627B (en) | Determine method and the device of vibroseis combination center | |
| Sam et al. | Estimation of dune celerity and sand flux in part of West Rajasthan, Gadra area of the Thar Desert using temporal remote sensing data | |
| CN104849751B (en) | The method of Prestack seismic data imaging | |
| CN103778657B (en) | Space partition based sound beam tracing method | |
| CN103365911B (en) | Map space indexing method and system based on two dimension partitioned structure | |
| CN105842691A (en) | Automatic adjustment and inlay method of large area InSAR deformation result and device thereof |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171128 Termination date: 20180820 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |