CN105405133A - Remote sensing image alteration detection method - Google Patents
Remote sensing image alteration detection method Download PDFInfo
- Publication number
- CN105405133A CN105405133A CN201510742564.7A CN201510742564A CN105405133A CN 105405133 A CN105405133 A CN 105405133A CN 201510742564 A CN201510742564 A CN 201510742564A CN 105405133 A CN105405133 A CN 105405133A
- Authority
- CN
- China
- Prior art keywords
- sigma
- remote sensing
- omega
- sensing image
- image
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 32
- 230000004075 alteration Effects 0.000 title claims abstract description 11
- 238000000034 method Methods 0.000 claims abstract description 35
- 230000003287 optical effect Effects 0.000 claims abstract description 35
- 238000012549 training Methods 0.000 claims abstract description 23
- 230000011218 segmentation Effects 0.000 claims abstract description 15
- 238000012937 correction Methods 0.000 claims abstract description 14
- 238000010606 normalization Methods 0.000 claims abstract description 14
- 230000005855 radiation Effects 0.000 claims abstract description 14
- 230000003595 spectral effect Effects 0.000 claims abstract description 6
- 238000005457 optimization Methods 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005192 partition Methods 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 description 8
- 238000012360 testing method Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 238000012952 Resampling Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Image Processing (AREA)
- Measurement Of Radiation (AREA)
Abstract
The invention discloses a remote sensing image alteration detection method comprising the steps that two time phase high-resolution optical remote sensing images X1 and X2 are obtained; image registration is performed on X1 and X2; radiation normalization correction is performed on X1 and X2 by utilizing a multivariate alteration detection method; alteration vector amplitude XM and spectral angle information XSA are respectively calculated according to X1 and X2 after radiation normalization correction; an optimal segmentation threshold T is calculated by utilizing the Bayes theory and an expectation maximization algorithm according to the XM; a pseudo training sample area is selected according to T and XM; XM and XSA are combined to act as input of a core FCM, and optimal model parameter value selection is performed on a core FCM combining spatial neighborhood information model according to the pseudo training sample area; and the alteration area and the non-alteration area of the optical remote sensing images are determined by adopting the method of core FCM combining spatial neighborhood information according to the selected optimal model parameter values. The remote sensing image alteration detection method is more robust and higher in precision.
Description
Technical field
The present invention relates to remote sensing image change detection techniques field, particularly relate to a kind of remote sensing image variation detection method.
Background technology
Along with the continuous accumulation of multidate high-definition remote sensing data and the foundation in succession of spatial database, how to extract from these remotely-sensed datas and change detected information has become the important subject of remote sensing science and Geographical Information Sciences.According to the remote sensing image of the different phase of the same area, the information of the dynamic change such as city, environment can be extracted, for resource management and the department such as planning, environmental protection provide the foundation of science decision.China " 12 " implements strengthening to expand the high resolving power earth observation engineering that Eleventh Five-Year Plan has started enforcement, concern comprises high-definition remote sensing target and basic theory and the key technology research such as space environment signature analysis and high reliability automatic interpretation, is becoming the research focus solving national security and the great demand of socio-economic development.
It is exactly in remotely-sensed data never of the same period that the change of remote sensing image detects, and analyzes quantitatively and determines the feature that earth's surface changes and process.Scholars proposes many effective detection algorithms from different angles and applied research, as Change vector Analysis method (ChangeVectorAnalysis, CVA), clustering method etc. based on FuzzyC-means (FCM).Wherein, traditional change of the multidate remote optical sensing based on FCM detects, and how first to carry out CVA conversion, then carries out FCM cluster to the amplitude of diverse vector, and then obtains changing testing result.In such technology, the deficiency using FCM is only applicable to spherical or ellipsoid shape cluster, and to noise and outlier (Outlier) very sensitivity thereof.In addition, only use the amplitude of diverse vector, original multispectral information is not excavated fully, not steadily and surely, precision is not high.
For the problems referred to above, many scholars attempt by adding that in FCM objective function the constraint of different spatial neighborhoods solves, but the complicated and target prior imformation of high resolution image testing environment is deficient, and cause these algorithms all to there is certain limitation, precision is not high.For this reason, be necessary that the new High Resolution Visible Light Remote Sensing Imagery Change Detection technology of research effectively overcomes above-mentioned difficult point.
Summary of the invention
Technical matters to be solved by this invention is, a kind of remote sensing image variation detection method is provided, the method is a kind of multi-temporal remote sensing image change detecting method of combining the self-adaptive kernel FCM of CVA and SAM, and the present invention changes that testing result is more sane, precision is higher.
In order to solve the problems of the technologies described above, the invention provides a kind of remote sensing image variation detection method, comprising:
Obtain two phase high-resolution optical remote sensing image X
1and X
2;
To optical remote sensing image X
1and X
2carry out Image registration;
Utilize Multivariate alteration detection method to optical remote sensing image X
1and X
2carry out radiation normalization correction;
Optical remote sensing image X after correcting according to radiation normalization
1and X
2calculate diverse vector amplitude X respectively
mwith spectral modeling information X
sA;
According to diverse vector amplitude X
mbayes principle and EM algorithm is utilized to calculate optimum segmentation threshold value T;
According to optimum segmentation threshold value T and diverse vector amplitude X
mselect pseudo-training sample region;
By X
mand X
sAcombine the input as core FCM, according to described pseudo-training sample region, in conjunction with space neighborhood information model, the selection of optimization model parameter value is carried out to core FCM;
According to the optimization model parameter value selected, adopt core FCM in conjunction with the method for space neighborhood information, determine region of variation and the non-changing region of optical remote sensing image.
Implement the present invention, there is following beneficial effect: the present invention combines the diverse vector amplitude of multi-temporal remote sensing image and the spectral modeling mapping graph (SpectralAngleMapper of multidate, SAM) as the input of core FCM, again based on the method for core FCM in conjunction with space neighborhood information, obtain final change testing result.Wherein, the nuclear parameter etc. in core FCM objective function, is selected by the pseudo-training sample based on CVA technical limit spacing, and change testing result is more sane, precision is higher.
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 the schematic flow sheet of an embodiment of remote sensing image variation detection method provided by the invention;
Fig. 2 is original high resolution optical remote sensing image figure;
Fig. 3 is the experimental result comparison diagram of the inventive method and additive method
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.
Fig. 1 is the schematic flow sheet of an embodiment of remote sensing image variation detection method provided by the invention, and the present invention is a kind of multi-temporal remote sensing image change detecting method, is mainly applicable to high-resolution optical remote sensing image, as shown in Figure 1, the present invention includes step:
S101, obtain two phase high-resolution optical remote sensing image X
1and X
2.
Wherein, X
1, X
2two panel height resolution Optical remote sensing images of the different phase of the same area.
S102, to optical remote sensing image X
1and X
2carry out Image registration.
Concrete, step S102 specifically comprises step:
S1021, employing ENVI14.8 remote sensing software are to optical remote sensing image X
1and X
2carry out geometric approximate correction.
Geometric approximate correction concrete operation step is: (1) display reference images and image to be corrected; (2) ground control point GCPs is gathered; GCPs should be evenly distributed in entire image, and the number of GCPs is at least more than or equal to 9; (3) error of calculation; (4) multinomial model is selected; (5) bilinear interpolation is adopted to carry out resampling output.Bilinearity differential technique is wherein: if ask unknown function f in a value of P=(x, y), supposes that our known function f is at Q
11=(x
1, y
1), Q
12=(x
1, y
2), Q
21=(x
2, y
1), and Q
22=(x
2, y
2) value of four points.If select a coordinate system to make the coordinate of these four points be respectively (0,0), (0,1), (1,0) and (1,1), so bilinear interpolation formula just can be expressed as:
f(x,y)≈f(0,0)(1-x)(1-y)+f(1,0)x(1-y)+f(0,1)(1-x)y+f(1,1)xy。
S1022, utilize Auto-matching and Triangulation Method to the X after geometric approximate correction
1and X
2carry out geometric accurate correction.
Wherein, Triangulation Method is, incremental algorithm is adopted to build Delaunay triangulation network, to each triangle, utilize the geographic coordinate of the ranks on its three summits number reference images same place corresponding with it to determine the affine Transform Model parameter of this triangle interior, treat correcting image to correct, obtain the remote sensing image after correcting.
S103, utilize Multivariate alteration detection method (MultivariateAlterationDetection, MAD) to optical remote sensing image X
1and X
2carry out radiation normalization correction.
Concrete, step S103 specifically comprises step:
S1031, acquisition optical remote sensing image X
1and X
2the linear combination of each wave band brightness value, obtains the difference image that change information strengthens;
S1032, according to described difference image by threshold value determination region of variation and non-region of variation;
S1033, by the right mapping equation of two phase pixels corresponding to non-region of variation, complete relative detector calibration.
S104, according to radiation normalization correct after optical remote sensing image X
1and X
2calculate diverse vector amplitude X respectively
mwith spectral modeling information X
sA.
Concrete, step S104 comprises step:
S1041, according to radiation normalization correct after optical remote sensing image X
1and X
2calculate diverse vector amplitude X
m.
Wherein,
In formula, B represents the wave band number of each phase remote sensing image, and (i, j) is the coordinate of image, X
1brepresent X
1b wave band image, X
2brepresent X
2b wave band image;
S1042, according to radiation normalization correct after optical remote sensing image X
1and X
2calculate diverse vector amplitude X
m,
Wherein,
S105, according to diverse vector amplitude X
mbayes principle and EM algorithm (Expectation-Maximization, EM) is utilized to calculate optimum segmentation threshold value T.
Concrete, step S105 specifically comprises step:
S1051, employing EM algorithm estimate X
mimage does not change class ω
naverage m
nand variances sigma
n, change class ω
caverage m
cbe σ with variance
c, wherein,
In formula, t represents iterations, and t subscript represents value during the t time iteration of Current Content, such as,
represent m
nvalue during the t+1 time iteration, other representation classes seemingly,
represent
value during the t+1 time iteration,
I and J represents line number and the columns of image respectively,
represent X
mimage does not change class ω
nthe Gaussian distribution of obeying,
represent X
mimage changes class ω
cthe Gaussian distribution of obeying;
S1052, according to Bayes minimum error principle, solution formula
Obtain optimum segmentation threshold value T.
S106, according to optimum segmentation threshold value T and diverse vector amplitude X
mselect pseudo-training sample region.
Concrete, step S106 comprises step:
S1061, according to optimum segmentation threshold value T and diverse vector amplitude X
mselect not change the pseudo-training set sample of class and be
S1062, according to optimum segmentation threshold value T and diverse vector amplitude X
mselect to change the pseudo-training set sample of class and be
wherein, δ is X
m15% of dynamic range.
S107, by X
mand X
sAcombine the input as core FCM, according to described pseudo-training sample region, in conjunction with space neighborhood information model, the selection of optimization model parameter value is carried out to core FCM.
Concrete, step S107 specifically comprises step:
S1071, by X
mand X
sAcombine the input as core FCM, building core FCM in conjunction with space neighborhood information model is:
In formula, C is clusters number, and N is the sum of sample,
represent the fuzzy membership of kth sample for jth class cluster centre, m is the weighted index of degree of membership,
parameter alpha controls punishment effect,
for X
mlocal mean value image and X
sAthe combination of local mean value image,
S1072, setup parameter α and nuclear parameter g span, utilize pseudo-training sample set, search variability index C
indexfor the value of α and g time minimum is as optimization model parameter value.
Wherein, variability index C
index=D
index/ k
t,
k
trepresent the Kappa coefficient of model parameter on pseudo-training sample set, N
n(α, g) represents the non-changing number of pixels of the whole image utilizing the minimization of object function to obtain when given α and g; N
c(α, g) represents when given α and g, the change number of pixels of whole image; TN
n(α, g) represents when given α and g, the non-changing number of pixels that pseudo-training sample is concentrated; TN
c(α, g) represents when given α and g, the change number of pixels that pseudo-training sample is concentrated.
S108, according to the optimization model parameter value selected, adopt core FCM in conjunction with the method for space neighborhood information, determine region of variation and the non-changing region of optical remote sensing image.
Concrete, step S108 specifically comprises:
S1081, setting core FCM, in conjunction with the clusters number C=2 in space neighborhood information model, as the initial center not changing class and change class, select and diverse vector amplitude X
mthe vector that minimum value is corresponding with maximal value; If Weighting exponent m=2 of degree of membership, ε be greater than 0 constant, the value of parameter alpha and nuclear parameter g is selected described optimization model parameter value;
S1082, calculating X
m, X
sAlocal window average, window size is set to 3 × 3;
S1083, employing formula
Upgrade fuzzy partition matrix;
S1084, employing formula
Upgrade cluster centre;
S1085, repeat to upgrade fuzzy partition matrix and cluster centre until the cluster centre cluster of adjacent twice iteration is less than ε;
S1086, according to fuzzy partition matrix u
jkdetermine final change detection figure, obtain region of variation and the non-changing region of optical remote sensing image.
Effect of the present invention further illustrates by following experimental result and analysis:
Experimental data of the present invention is the multidate SPOT high-resolution image data in French Littoral area, and image size is 400 × 400, uses B1, B2 and B3 tri-wave bands.In order to verify validity of the present invention, change detecting method of the present invention and following change detecting method are compared:
(1) based on the EM method (CVA-EM) of CVA, [gondola BruzzoneL. etc. are at article " Automaticanalysisofdifferenceimageforunsupervisedchanged etection " (IEEETransactionsonGeoscienceandRemoteSensing, 2000,38 (3): 1171-1182.) detection method carried in].
(2) [Chensongchan etc. are at article " RobustImageSegmentationUsingFCMWithSpatialConstraintsBas edonNewKernel-InducedDistanceMeasure " (IEEETransactionsonSystems in conjunction with the sorting technique (FCM-S) of space neighborhood information for FCM, Man, andCybernetics-PartB:Cybernetics, 2004,34 (4): 1907-1916.) method carried in]
(3) the inventive method.
Detection perform false retrieval number FP, undetected several FN, total error number OE and Kappa coefficient four indexs are weighed.FP, FN and OE more close to 0, Kappa coefficient more close to 1, show that the performance of change detecting method is better.Testing result is as shown in table 1.From Fig. 2, Fig. 3 and table 1, the detection method performance that the present invention carries is better than other two kinds of detection methods, and this shows that the change detecting method that the present invention carries is effective.
The multidate SPOT5 remote sensing imagery change detection results contrast in table 1Littoral area
Method | FP | FN | OE | k |
CVA-EM | 7919 | 3882 | 11801 | 0.705 |
FCM-S | 1822 | 6928 | 8750 | 0.737 |
The inventive method | 2511 | 4689 | 7200 | 0.797 |
Desirable | 0 | 0 | 0 | 1 |
It should be noted that, in this article, term " comprises ", " comprising " or its any other variant are intended to contain comprising of nonexcludability, thus make to comprise the process of a series of key element, method, article or device and not only comprise those key elements, but also comprise other key elements clearly do not listed, or also comprise by the intrinsic key element of this process, method, article or device.When not more restrictions, the key element limited by statement " comprising ... ", and be not precluded within process, method, article or the device comprising this key element and also there is other identical element.
In several embodiments that the application provides, should be understood that, disclosed method can realize by another way.Such as, system embodiment described above is only schematic, professional can also recognize further, in conjunction with unit and the algorithm steps of each example of embodiment disclosed herein description, can realize with electronic hardware, computer software or the combination of the two, in order to the interchangeability of hardware and software is clearly described, generally describe composition and the step of each example in the above description according to function.These functions perform with hardware or software mode actually, depend on application-specific and the design constraint of technical scheme.Professional and technical personnel can use distinct methods to realize described function to each specifically should being used for, but this realization should not thought and exceeds scope of the present invention.Software module can be placed in the storage medium of other form any known in random access memory (RAM), internal memory, ROM (read-only memory) (ROM), electrically programmable ROM, electrically erasable ROM, register, hard disk, moveable magnetic disc, CD-ROM or technical field.
To the above-mentioned explanation of the disclosed embodiments, professional and technical personnel in the field are realized or uses the present invention.To be apparent for those skilled in the art to the multiple amendment of these embodiments, General Principle as defined herein can without departing from the spirit or scope of the present invention, realize in other embodiments.Therefore, the present invention can not be restricted to these embodiments shown in this article, but will meet the widest scope consistent with principle disclosed herein and features of novelty.
Claims (8)
1. a remote sensing image variation detection method, is characterized in that, comprising:
Obtain two phase high-resolution optical remote sensing image X
1and X
2;
To optical remote sensing image X
1and X
2carry out Image registration;
Utilize Multivariate alteration detection method to optical remote sensing image X
1and X
2carry out radiation normalization correction;
Optical remote sensing image X after correcting according to radiation normalization
1and X
2calculate diverse vector amplitude X respectively
mwith spectral modeling information X
sA;
According to diverse vector amplitude X
mbayes principle and EM algorithm is utilized to calculate optimum segmentation threshold value T;
According to optimum segmentation threshold value T and diverse vector amplitude X
mselect pseudo-training sample region;
By X
mand X
sAcombine the input as core FCM, according to described pseudo-training sample region, in conjunction with space neighborhood information model, the selection of optimization model parameter value is carried out to core FCM;
According to the optimization model parameter value selected, adopt core FCM in conjunction with the method for space neighborhood information, determine region of variation and the non-changing region of optical remote sensing image.
2. remote sensing image variation detection method as claimed in claim 1, is characterized in that, described to optical remote sensing image X
1and X
2carry out Image registration, specifically comprise:
Adopt ENVI14.8 remote sensing software to optical remote sensing image X
1and X
2carry out geometric approximate correction;
Utilize Auto-matching and Triangulation Method to the X after geometric approximate correction
1and X
2carry out geometric accurate correction.
3. remote sensing image variation detection method as claimed in claim 1, it is characterized in that, the described Multivariate alteration detection method that utilizes is to optical remote sensing image X
1and X
2carry out radiation normalization correction, specifically comprise:
Obtain optical remote sensing image X
1and X
2the linear combination of each wave band brightness value, obtains the difference image that change information strengthens;
According to described difference image by threshold value determination region of variation and non-region of variation;
By the mapping equation that the two phase pixels that non-region of variation is corresponding are right, complete relative detector calibration.
4. remote sensing image variation detection method as claimed in claim 1, is characterized in that, described according to the optical remote sensing image X after radiation normalization correction
1and X
2calculate diverse vector amplitude X respectively
mwith spectral modeling information X
sA, specifically comprise:
Optical remote sensing image X after correcting according to radiation normalization
1and X
2calculate diverse vector amplitude X
m, wherein,
In formula, B represents the wave band number of each phase remote sensing image, and (i, j) is the coordinate of image, X
1brepresent X
1b wave band image, X
2brepresent X
2b wave band image;
Optical remote sensing image X after correcting according to radiation normalization
1and X
2calculate diverse vector amplitude X
m, wherein,
5. remote sensing image variation detection method as claimed in claim 1, is characterized in that, described according to diverse vector amplitude X
mutilize Bayes principle and EM algorithm to calculate optimum segmentation threshold value T, specifically comprise:
EM algorithm is adopted to estimate X
mimage does not change class ω
naverage m
nand variances sigma
n, change class ω
caverage m
cbe σ with variance
c, wherein,
In formula, t represents iterations, and t subscript represents value during the t time iteration of Current Content,
I and J represents line number and the columns of image respectively,
represent X
mimage does not change class ω
nthe Gaussian distribution of obeying,
represent X
mimage changes class ω
cthe Gaussian distribution of obeying;
According to Bayes minimum error principle, solution formula
Obtain optimum segmentation threshold value T.
6. remote sensing image variation detection method as claimed in claim 1, is characterized in that, described according to optimum segmentation threshold value T and diverse vector amplitude X
mselect pseudo-training sample region, specifically comprise:
According to optimum segmentation threshold value T and diverse vector amplitude X
mselect not change the pseudo-training set sample of class and be
According to optimum segmentation threshold value T and diverse vector amplitude X
mselect to change the pseudo-training set sample of class and be
Wherein, δ is X
m15% of dynamic range.
7. remote sensing image variation detection method as claimed in claim 1, is characterized in that, described by X
mand X
sAcombine the input as core FCM, according to described pseudo-training sample region, in conjunction with space neighborhood information model, optimization model Selecting parameter carried out to core FCM, specifically comprise:
By X
mand X
sAcombine the input as core FCM, building core FCM in conjunction with space neighborhood information model is:
Setup parameter α and nuclear parameter g span, utilize pseudo-training sample set, search variability index C
indexfor the value of α and g time minimum is as optimization model parameter value; Wherein, variability index C
index=D
index/ k
t,
k
trepresent the Kappa coefficient of model parameter on pseudo-training sample set, N
n(α, g) represents the non-changing number of pixels of the whole image utilizing the minimization of object function to obtain when given α and g; N
c(α, g) represents when given α and g, the change number of pixels of whole image; TN
n(α, g) represents when given α and g, the non-changing number of pixels that pseudo-training sample is concentrated; TN
c(α, g) represents when given α and g, the change number of pixels that pseudo-training sample is concentrated.
8. remote sensing image variation detection method as claimed in claim 7, it is characterized in that the described optimization model parameter value according to selecting adopts core FCM in conjunction with the method for space neighborhood information, determine region of variation and the non-changing region of high-resolution optical remote sensing image, specifically comprise:
Setting core FCM, in conjunction with the clusters number C=2 in space neighborhood information model, as the initial center not changing class and change class, selects and diverse vector amplitude X
mthe vector that minimum value is corresponding with maximal value; If Weighting exponent m=2 of degree of membership, ε be greater than 0 constant, the value of parameter alpha and nuclear parameter g is selected described optimization model parameter value;
Calculate X
m, X
sAlocal window average, window size is set to 3 × 3;
Employing formula
Upgrade fuzzy partition matrix;
Employing formula
Upgrade cluster centre;
Repeat to upgrade fuzzy partition matrix and cluster centre until the cluster centre cluster of adjacent twice iteration is less than ε;
According to fuzzy partition matrix u
jkdetermine final change detection figure, obtain region of variation and the non-changing region of optical remote sensing image.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510742564.7A CN105405133B (en) | 2015-11-04 | 2015-11-04 | A kind of remote sensing image variation detection method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510742564.7A CN105405133B (en) | 2015-11-04 | 2015-11-04 | A kind of remote sensing image variation detection method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105405133A true CN105405133A (en) | 2016-03-16 |
CN105405133B CN105405133B (en) | 2018-01-19 |
Family
ID=55470600
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510742564.7A Expired - Fee Related CN105405133B (en) | 2015-11-04 | 2015-11-04 | A kind of remote sensing image variation detection method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105405133B (en) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106373120A (en) * | 2016-09-09 | 2017-02-01 | 河海大学 | Multi-temporal remote sensing image change detection method based on non-negative matrix decomposition and nucleus FCM |
CN106372612A (en) * | 2016-09-09 | 2017-02-01 | 河海大学 | Multi-temporal remote sensing image change detection method combining FCM with MRF model |
CN106384352A (en) * | 2016-09-09 | 2017-02-08 | 河海大学 | Multi-temporal remote sensing image change detection method based on fusion strategy and FCM |
CN106646469A (en) * | 2016-12-21 | 2017-05-10 | 中国科学院遥感与数字地球研究所 | SAR (Synthetic Aperture Radar) ship detection optimization method based on variation coefficient method |
CN107248172A (en) * | 2016-09-27 | 2017-10-13 | 中国交通通信信息中心 | A kind of remote sensing image variation detection method based on CVA and samples selection |
CN107346549A (en) * | 2017-06-09 | 2017-11-14 | 中国矿业大学 | A kind of multi-class change dynamic threshold detection method using remote sensing image multiple features |
CN107480634A (en) * | 2017-08-12 | 2017-12-15 | 天津市测绘院 | A kind of geographical national conditions ground mulching monitoring method based on multistage target classification |
CN107481235A (en) * | 2017-08-24 | 2017-12-15 | 河海大学 | The multi-temporal remote sensing image change detecting method that a kind of mathematical morphology filter converts with reference to card side |
CN107578040A (en) * | 2017-09-30 | 2018-01-12 | 中南大学 | A kind of house change detecting method based on Pulse Coupled Neural Network |
CN107688777A (en) * | 2017-07-21 | 2018-02-13 | 同济大学 | A kind of urban green space extracting method for cooperateing with multi-source Remote Sensing Images |
CN107992891A (en) * | 2017-12-01 | 2018-05-04 | 西安电子科技大学 | Based on spectrum vector analysis multi-spectral remote sensing image change detecting method |
CN109146933A (en) * | 2017-06-28 | 2019-01-04 | 中国石油化工股份有限公司 | Multiple dimensioned digital cores modeling method and computer readable storage medium |
CN109191503A (en) * | 2018-08-23 | 2019-01-11 | 河海大学 | Remote sensing image variation detection method and system based on condition random field |
CN109300115A (en) * | 2018-09-03 | 2019-02-01 | 河海大学 | A kind of multispectral high-resolution remote sensing image change detecting method of object-oriented |
CN109934799A (en) * | 2016-09-09 | 2019-06-25 | 河海大学 | Multidate difference image magnitude calculation and change detecting method |
CN110232302A (en) * | 2018-03-06 | 2019-09-13 | 香港理工大学深圳研究院 | A kind of change detecting method of integrated gray value, spatial information and classification knowledge |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1760888A (en) * | 2005-11-03 | 2006-04-19 | 复旦大学 | Method for recognizing change of earth's surface by using satellite SAR carried images at multiple time phases |
CN101950361A (en) * | 2010-09-06 | 2011-01-19 | 中国科学院遥感应用研究所 | Adaptive extraction method of remote sensing image thematic information based on spectrum matching degree |
CN103150580A (en) * | 2013-03-18 | 2013-06-12 | 武汉大学 | Method and device for Hyperspectral image semi-supervised classification |
CN104200458A (en) * | 2014-07-30 | 2014-12-10 | 浙江工业大学 | MeanShift based high-resolution remote sensing image segmentation distance measurement optimization method |
CN104751166A (en) * | 2013-12-30 | 2015-07-01 | 中国科学院深圳先进技术研究院 | Spectral angle and Euclidean distance based remote-sensing image classification method |
-
2015
- 2015-11-04 CN CN201510742564.7A patent/CN105405133B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1760888A (en) * | 2005-11-03 | 2006-04-19 | 复旦大学 | Method for recognizing change of earth's surface by using satellite SAR carried images at multiple time phases |
CN101950361A (en) * | 2010-09-06 | 2011-01-19 | 中国科学院遥感应用研究所 | Adaptive extraction method of remote sensing image thematic information based on spectrum matching degree |
CN103150580A (en) * | 2013-03-18 | 2013-06-12 | 武汉大学 | Method and device for Hyperspectral image semi-supervised classification |
CN104751166A (en) * | 2013-12-30 | 2015-07-01 | 中国科学院深圳先进技术研究院 | Spectral angle and Euclidean distance based remote-sensing image classification method |
CN104200458A (en) * | 2014-07-30 | 2014-12-10 | 浙江工业大学 | MeanShift based high-resolution remote sensing image segmentation distance measurement optimization method |
Non-Patent Citations (4)
Title |
---|
TAUQIR AHMED MOUGHAL AND FUSHENG YU: "An Automatic Unsupervised Method Based on Context-Sensitive Spectral Angle Mapper for Change Detection of Remote Sensing Images", 《ADMA2014》 * |
余先川 等: "一种基于光谱角和光谱距离自动加权融合分类方法", 《地质学刊》 * |
孙艳丽 等: "光谱角-欧氏距离的高光谱图像辐射归一化", 《遥感学报》 * |
赵磊 等: "基于模糊C均值聚类和了邻域分析的无监督多通道遥感图像变化检测", 《数据采集与处理》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106373120A (en) * | 2016-09-09 | 2017-02-01 | 河海大学 | Multi-temporal remote sensing image change detection method based on non-negative matrix decomposition and nucleus FCM |
CN106372612A (en) * | 2016-09-09 | 2017-02-01 | 河海大学 | Multi-temporal remote sensing image change detection method combining FCM with MRF model |
CN106384352A (en) * | 2016-09-09 | 2017-02-08 | 河海大学 | Multi-temporal remote sensing image change detection method based on fusion strategy and FCM |
CN109934799A (en) * | 2016-09-09 | 2019-06-25 | 河海大学 | Multidate difference image magnitude calculation and change detecting method |
CN109934799B (en) * | 2016-09-09 | 2021-06-25 | 南京工程学院 | Multi-time-phase difference image module value calculation and change detection method |
CN106373120B (en) * | 2016-09-09 | 2019-01-08 | 河海大学 | Multi-temporal remote sensing image change detecting method based on Non-negative Matrix Factorization and core FCM |
CN107248172A (en) * | 2016-09-27 | 2017-10-13 | 中国交通通信信息中心 | A kind of remote sensing image variation detection method based on CVA and samples selection |
CN106646469B (en) * | 2016-12-21 | 2019-01-29 | 中国科学院遥感与数字地球研究所 | SAR ship detection optimization method based on VC Method |
CN106646469A (en) * | 2016-12-21 | 2017-05-10 | 中国科学院遥感与数字地球研究所 | SAR (Synthetic Aperture Radar) ship detection optimization method based on variation coefficient method |
CN107346549A (en) * | 2017-06-09 | 2017-11-14 | 中国矿业大学 | A kind of multi-class change dynamic threshold detection method using remote sensing image multiple features |
CN107346549B (en) * | 2017-06-09 | 2020-04-14 | 中国矿业大学 | Multi-class change dynamic threshold detection method utilizing multiple features of remote sensing image |
CN109146933A (en) * | 2017-06-28 | 2019-01-04 | 中国石油化工股份有限公司 | Multiple dimensioned digital cores modeling method and computer readable storage medium |
CN109146933B (en) * | 2017-06-28 | 2020-12-01 | 中国石油化工股份有限公司 | Multi-scale digital core modeling method and computer-readable storage medium |
CN107688777A (en) * | 2017-07-21 | 2018-02-13 | 同济大学 | A kind of urban green space extracting method for cooperateing with multi-source Remote Sensing Images |
CN107480634A (en) * | 2017-08-12 | 2017-12-15 | 天津市测绘院 | A kind of geographical national conditions ground mulching monitoring method based on multistage target classification |
CN107481235A (en) * | 2017-08-24 | 2017-12-15 | 河海大学 | The multi-temporal remote sensing image change detecting method that a kind of mathematical morphology filter converts with reference to card side |
CN107578040A (en) * | 2017-09-30 | 2018-01-12 | 中南大学 | A kind of house change detecting method based on Pulse Coupled Neural Network |
CN107992891A (en) * | 2017-12-01 | 2018-05-04 | 西安电子科技大学 | Based on spectrum vector analysis multi-spectral remote sensing image change detecting method |
CN107992891B (en) * | 2017-12-01 | 2022-01-25 | 西安电子科技大学 | Multispectral remote sensing image change detection method based on spectral vector analysis |
CN110232302A (en) * | 2018-03-06 | 2019-09-13 | 香港理工大学深圳研究院 | A kind of change detecting method of integrated gray value, spatial information and classification knowledge |
CN109191503A (en) * | 2018-08-23 | 2019-01-11 | 河海大学 | Remote sensing image variation detection method and system based on condition random field |
CN109300115A (en) * | 2018-09-03 | 2019-02-01 | 河海大学 | A kind of multispectral high-resolution remote sensing image change detecting method of object-oriented |
CN109300115B (en) * | 2018-09-03 | 2021-11-30 | 河海大学 | Object-oriented multispectral high-resolution remote sensing image change detection method |
Also Published As
Publication number | Publication date |
---|---|
CN105405133B (en) | 2018-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105405133A (en) | Remote sensing image alteration detection method | |
Teegavarapu et al. | Geo-spatial grid-based transformations of precipitation estimates using spatial interpolation methods | |
Szcześniak et al. | A method for using street view imagery to auto-extract window-to-wall ratios and its relevance for urban-level daylighting and energy simulations | |
Yang et al. | Analyzing big time series data in solar engineering using features and PCA | |
CN105389817A (en) | Two-time phase remote sensing image change detection method | |
CN105551031A (en) | Multi-temporal remote sensing image change detection method based on FCM and evidence theory | |
CN103745459B (en) | Detection method of an unstructured point cloud feature point and extraction method thereof | |
CN102222313B (en) | Urban evolution simulation structure cell model processing method based on kernel principal component analysis (KPCA) | |
CN103456020B (en) | Based on the method for detecting change of remote sensing image of treelet Fusion Features | |
Chandra et al. | Small area prediction of counts under a non-stationary spatial model | |
CN104751185A (en) | SAR image change detection method based on mean shift genetic clustering | |
CN106372612A (en) | Multi-temporal remote sensing image change detection method combining FCM with MRF model | |
CN106447653B (en) | The multi-temporal remote sensing image change detecting method that card side based on space constraint converts | |
CN107481235A (en) | The multi-temporal remote sensing image change detecting method that a kind of mathematical morphology filter converts with reference to card side | |
Kim et al. | Automated classification of thermal defects in the building envelope using thermal and visible images | |
CN105095913A (en) | Remote sensing image classification method and system based on neighbor regular joint sparse representation | |
CN106650571B (en) | A kind of multi-temporal remote sensing image change detecting method converted based on adaptive card side | |
CN106373120B (en) | Multi-temporal remote sensing image change detecting method based on Non-negative Matrix Factorization and core FCM | |
CN104517286A (en) | SAR (synthetic aperture radar) image registration based on self-adaption threshold segmentation and combination optimization | |
CN112668461A (en) | Intelligent supervision system with wild animal identification function | |
Palacios-Rodríguez et al. | Generalized Pareto processes for simulating space-time extreme events: an application to precipitation reanalyses | |
CN106096622A (en) | Semi-supervised Classification of hyperspectral remote sensing image mask method | |
Fan | Research on deep learning energy consumption prediction based on generating confrontation network | |
CN105354845B (en) | A kind of semi-supervised change detecting method of remote sensing image | |
CN105160666B (en) | SAR image change detection based on Non-Stationary Analysis and condition random field |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180119 |