CN109389571A - 一种遥感影像变化检测方法、装置及终端 - Google Patents
一种遥感影像变化检测方法、装置及终端 Download PDFInfo
- Publication number
- CN109389571A CN109389571A CN201710655794.9A CN201710655794A CN109389571A CN 109389571 A CN109389571 A CN 109389571A CN 201710655794 A CN201710655794 A CN 201710655794A CN 109389571 A CN109389571 A CN 109389571A
- Authority
- CN
- China
- Prior art keywords
- differential images
- group
- pixel
- remote sensing
- sensing 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 96
- 230000008859 change Effects 0.000 claims abstract description 91
- 238000012360 testing method Methods 0.000 claims abstract description 71
- 230000004927 fusion Effects 0.000 claims abstract description 41
- 239000011159 matrix material Substances 0.000 claims abstract description 35
- 238000005192 partition Methods 0.000 claims abstract description 30
- 238000005516 engineering process Methods 0.000 claims abstract description 17
- 230000000295 complement effect Effects 0.000 claims abstract description 16
- 238000000034 method Methods 0.000 claims description 40
- 238000004590 computer program Methods 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 14
- 238000003860 storage Methods 0.000 claims description 14
- 238000007689 inspection Methods 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000007499 fusion processing Methods 0.000 claims description 5
- 235000013399 edible fruits Nutrition 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 20
- 230000006870 function Effects 0.000 description 16
- 238000010586 diagram Methods 0.000 description 10
- 230000010354 integration Effects 0.000 description 5
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 4
- 239000010931 gold Substances 0.000 description 4
- 229910052737 gold Inorganic materials 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000017105 transposition Effects 0.000 description 3
- 239000004365 Protease Substances 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- POTUGHMKJGOKRI-UHFFFAOYSA-N ficin Chemical compound FI=CI=N POTUGHMKJGOKRI-UHFFFAOYSA-N 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 238000004080 punching Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- 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
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
-
- 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
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20224—Image subtraction
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Probability & Statistics with Applications (AREA)
- Image Analysis (AREA)
Abstract
本发明适用于遥感技术领域,提供了一种遥感影像变化检测方法、装置及终端,通过对两期遥感影像进行预处理后生成信息互补的差分影像组集合;对差分影像组集合进行模糊C均值聚类得到每组差分影像的模糊分割矩阵和变化检测图;根据得到的每组差分影像的模糊分割矩阵和变化检测图,通过相似度理论和模糊积分理论对差分影像组集合进行决策级融合,得到初始的变化检测结果;计算差分影像组集合在每个像元处的证据冲突程度,探测并定位初始变化检测结果中的证据冲突严重像元,并通过指示克里金理论对证据冲突严重像元重新进行分类。基于信息互补的差分影像组集合和指示克里金理论,实现了提高遥感变化检测的可靠性和变化检测结果的精度。
Description
技术领域
本发明属于遥感技术领域,尤其涉及一种遥感影像变化检测方法、装置及终端。
背景技术
遥感变化检测是通过遥感影像获取地表变化信息的重要手段,已广泛渗透到国民经济的各个领域。最常用的遥感变化检测方法是基于“差分影像”的方法,该方法首先通过对比两期遥感影像生成差分影像,然后通过分析差分影像获取地表变化信息。
由于遥感影像本身的复杂性,以及差异算子的局限性,使得生成的差分影像总是存在不确定性,影响了变化检测的效果。
现有基于融合技术对多差分影像进行融合以提高变化检测的精度,但对证据冲突严重像元的处理效果依然较差,检测精度不高。
发明内容
有鉴于此,本发明实施例提供了一种遥感影像变化检测方法、装置及终端,以降低现有差分影像变化检测的不确定性,提高遥感影像变化检测的可靠性和变化检测结果的精度。
本发明实施例的第一方面提供了一种遥感影像变化检测方法,包括:
获取两期遥感影像,并对所述遥感影像进行预处理;
根据预处理后的两期遥感影像生成信息互补的差分影像组集合;
对所述差分影像组集合进行模糊C均值聚类计算,得到每组差分影像的模糊分割矩阵和变化检测图;
根据所述的每组差分影像的模糊分割矩阵和变化检测图,通过相似度理论和模糊积分理论对所述差分影像组集合进行决策级融合,得到初始变化检测结果;
计算所述差分影像组集合在每个像元处的证据冲突程度;探测并定位所述初始变化检测结果中的证据冲突严重像元,通过指示克里金理论对所述证据冲突严重像元重新进行分类。
本发明实施例的第二方面提供了一种遥感影像变化检测装置,包括:
预处理单元,用于获取两期遥感影像,并对所述遥感影像进行预处理;
差分影像组集合生成单元,用于根据预处理后的两期遥感影像生成信息互补的差分影像组集合;
模糊C均值聚类计算单元,用于对所述差分影像组集合进行模糊C均值聚类计算,得到每组差分影像的模糊分割矩阵和变化检测图;
决策级融合单元,用于根据所述的每组差分影像的模糊分割矩阵和变化检测图,通过相似度理论和模糊积分理论对所述差分影像组集合进行决策级融合,得到初始变化检测结果;
证据冲突严重像元重分类处理单元,用于计算所述差分影像组集合在每个像元处的证据冲突程度;探测并定位所述初始变化检测结果中证据冲突严重像元,通过指示克里金理论对所述证据冲突严重像元重新进行分类。
本发明实施例的第三方面提供了一种终端,包括:
存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其中,所述处理器执行所述计算机程序时实现本发明实施例的第一方面提供的遥感影像变化检测方法的步骤。
其中,所述计算机程序包括:
预处理单元,用于获取两期遥感影像,并对所述遥感影像进行预处理;
差分影像组集合生成单元,用于根据预处理后的两期遥感影像生成信息互补的差分影像组集合;
模糊C均值聚类计算单元,用于对所述差分影像组集合进行模糊C均值聚类计算,得到每组差分影像的模糊分割矩阵和变化检测图;
决策级融合单元,用于根据所述的每组差分影像的模糊分割矩阵和变化检测图,通过相似度理论和模糊积分理论对所述差分影像组集合进行决策级融合,得到初始变化检测结果;
证据冲突严重像元重分类处理单元,用于计算所述差分影像组集合在每个像元处的证据冲突程度;探测并定位所述初始变化检测结果中证据冲突严重像元,通过指示克里金理论对所述证据冲突严重像元重新进行分类。
本发明实施例的第四方面提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其中,所述计算机程序被处理器执行时实现本发明实施例的第一方面提供的遥感影像变化检测方法的步骤。
其中,所述计算机程序包括:
预处理单元,用于获取两期遥感影像,并对所述遥感影像进行预处理;
差分影像组集合生成单元,用于根据预处理后的两期遥感影像生成信息互补的差分影像组集合;
模糊C均值聚类计算单元,用于对所述差分影像组集合进行模糊C均值聚类计算,得到每组差分影像的模糊分割矩阵和变化检测图;
决策级融合单元,用于根据所述的每组差分影像的模糊分割矩阵和变化检测图,通过相似度理论和模糊积分理论对所述差分影像组集合进行决策级融合,得到初始变化检测结果;
证据冲突严重像元重分类处理单元,用于计算所述差分影像组集合在每个像元处的证据冲突程度;探测并定位所述初始变化检测结果中的证据冲突严重像元,通过指示克里金理论对所述证据冲突严重像元重新进行分类。
本发明实施例与现有技术相比存在的有益效果是:首先从距离矩阵的角度,提出一种生成信息互补差分影像组集合的方案,充分提取两期遥感影像的变化信息;其次通过相似度理论和模糊积分理论实现对所述差分影像组集合的决策级融合;然后利用指示克里金技术对所述融合结果进行优化,重新分类融合结果中的证据冲突严重像元。通过以上三点实现了提高变化检测的可靠性和变化检测结果的精度。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种遥感影像变化检测方法的实现流程图;
图2是本发明实施例提供的遥感影像变化检测方法中步骤S104的具体实现流程图;
图3是本发明实施例提供的遥感影像变化检测方法中步骤S105的具体实现流程图;
图4是本发明实施例提供的一种利用指示克里金理论对证据冲突严重像元重新进行分类的方法具体实现流程图;
图5是一个7×7的克里金窗口示意图;
图6(a)是本发明实施例提供的T1时刻遥感影像第4波段的示意图;
图6(b)是本发明实施例提供的T2时刻遥感影像第4波段的示意图;
图6(c)是本发明实施例提供的T1时刻和T2时刻遥感影像的变化参考图;
图7(a)是模糊C均值根据差分影像DI1得到的变化检测图;
图7(b)是模糊C均值根据差分影像DI2得到的变化检测图;
图7(c)是模糊C均值根据差分影像DI3得到的变化检测图;
图7(d)是模糊C均值根据差分影像DI4得到的变化检测图;
图7(e)是模糊C均值根据差分影像DI5得到的变化检测图;
图7(f)是模糊C均值根据差分影像DI6得到的变化检测图;
图7(g)是通过融合技术投票法得到的变化检测结果示意图;
图7(h)是通过融合技术模糊积分法得到的变化检测结果示意图;
图7(i)是本发明实施例提供的遥感影像变化检测方法得到的变化检测结果示意图;
图8是本发明实施例提供的一种遥感影像变化检测装置的示意图;
图9是本发明实施例提供的一种终端的示意图。
具体实施方式
以下描述中,为了说明而不是为了限定,提出了诸如特定系统结构、技术之类的具体细节,以便透彻理解本发明实施例。然而,本领域的技术人员应当清楚,在没有这些具体细节的其他实施例中也可以实现本发明。在其他情况中,省略对众所周知的系统、装置、电路以及方法的详细说明,以免不必要的细节妨碍本发明的描述。
为了说明本发明所述的技术方案,下面通过具体实施例来进行说明。请参考图1,图1示出了本发明实施例提供的一种遥感影像变化检测方法的实现流程,详述如下:
在步骤S101中,获取两期遥感影像,并对所述遥感影像进行预处理。
在本发明实施例中,所述两期遥感影像分别为T1时刻和T2时刻的遥感影像,其中,T1≠T2。
所述预处理为遥感领域中的基础步骤,包括但不限于配准和辐射校正等。
在步骤S102中,根据预处理后的两期遥感影像生成信息互补的差分影像组集合。
在本发明实施例中,所述信息互补,是指所提取的信息具有互补性,能够充分反应两期遥感影像的变化。
在这里,上述步骤S102具体为,根据预处理后的两期遥感影像从距离矩阵的角度生成信息互补的差分影像组集合。所述差分影像组集合具体是指由六组信息互补的差分影像所组成的集合。
在这里,为便于描述,将进行预处理后的T1时刻的遥感影像记为X1,将进行预处理后的T2时刻的遥感影像记为X2。所述的进行预处理后的两期遥感影像X1和X2均包括B个波段,表示Xm的第b个波段,1<b≤B,m=1,2。
为充分提取两期遥感影像X2和X1的变化信息,并使所提取的信息具有互补性,本发明选取典型的距离度量生成六组差分影像:用DI={DI1,DI2,...,DI6}表示所生成的差分影像组集合,其中,DIn表示第n组差分影像,n=1,2,...,6。其中差分影像是逐像元计算的,用p(i,j)表示坐标值为(i,j)的像元。
其中,第一组差分影像DI1由使用最广泛的变化矢量分析技术生成,其本质上是根据欧氏距离公式以及两期遥感影像X2和X1所计算得到的欧氏距离:
这里DI1(i,j)表示第一组差分影像DI1在像元p(i,j)的取值。
第二组差分影像DI2为根据闵可夫斯基距离公式以及两期遥感影像X2和X1所计算得到的闵可夫斯基距离,即其中,p为变化参数,p≥1,在这里,优选p=10。
第三组差分影像DI3为根据马氏距离公式以及两期遥感影像X2和X1所计算得到的马氏距离,用符号ΣDX表示变化矢量的协方差矩阵,则这里(ΣDX)-1为ΣDX的逆矩阵,DX(i,j)T为DX(i,j)的转置。
第四组差分影像DI4为根据相关系数以及两期遥感影像X2和X1所计算得到的相关性距离:首先计算X2(i,j)和X1(i,j)的相关系数ρ(X2(i,j),X1(i,j)):
这里ρ(X2(i,j),X1(i,j))的取值范围为-1到1,为X2(i,j)各个波段的均值,为X1(i,j)各个波段的均值,1<b≤B,B为X1或X2中包含的波段个数;那么DI4(i,j)=1-ρ(X2(i,j),X1(i,j))。
第五组差分影像DI5通过将主成分分析理论应用到两期遥感影像X2和X1生成:首先通过比值算子求取两期遥感影像X2和X1的比值变化矢量然后计算该比值变化矢量RX的协方差矩阵ΣRX;假设ΣRX的特征根分别为β1,β2,...,βF,且满足β1≥β2≥...≥βF≥0,其中F表示矩阵RX的特征根个数,F>0。同时假设ΣRX的各个特征根对应的标准正交特征向量为γ1,γ2,...,γF,则比值变化矢量RX的第f个主成分为Yf=γfRXT,f=1,2,…F,其中,RXT为RX的转置;将各主成分的方差贡献率作为权值,对比值变化矢量RX的各主成分进行加权平均生成第五组差分影像DI5:其中,αf=βf/(β1+β2+…+βF),这里αf为第f个主成分的方差贡献率。
第六组差分影像DI6基于波谱形状信息,根据波谱曲线的波谱梯度差异生成:首先分别计算两期遥感影像X2和X1的波谱梯度矢量和其中m=1,2;b=1,2,…,B-1,这里λb为第b个波段的波长;然后结合欧氏距离公式以及波谱梯度矢量G1和G2生成第六组差分影像DI6,即
优选的,为避免各组差分影像间可能存在的不稳定性和不一致性,步骤S102还包括标准化过程:根据预设公式将每组差分影像的取值进行标准化使得其取值在预设取值区间范围之内。具体的,在得到差分影像组集合DI后,通过公式将每组差分影像的取值标准化到[0,255]区间,这里n=1,2,...,6,max(DIn)和min(DIn)分别表示第n组差分影像DIn的最大值和最小值,int表示取整运算。
在步骤S103中,对所述差分影像组集合进行模糊C均值聚类计算,得到每组差分影像的模糊分割矩阵和变化检测图。
在本发明实施例中,分别对所生成的六组差分影像DI={DI1,DI2,...,DI6}实施模糊C均值聚类计算,得到每组差分影像DIn的模糊分割矩阵Un和变化检测图CDn。
步骤S103是对差分影像组集合如所生成的六组差分影像进行决策级融合的准备阶段。在这里,模糊C均值(fuzzy c-means,FCM)聚类是遥感领域的基础技术,模糊分割矩阵变化检测图CDn=CDn(i,j)。其中,上标T表示转置,n=1,2,...,6,表示FCM根据每组差分影像DIn得到的像元p(i,j)属于未变化类wu的隶属度;表示FCM根据每组差分影像DIn得到的像元p(i,j)属于变化类wc的隶属度;变化检测图CDn在像元p(i,j)处的取值通过最大隶属度原则得到:其中CDn(i,j)=1表示根据第n组差分影像DIn像元p(i,j)被标记为变化类,CDn(i,j)=0表示根据第n组差分影像DIn像元p(i,j)被标记为未变化类。
在步骤S104中,根据所述的每组差分影像的模糊分割矩阵和变化检测图,通过相似度理论和模糊积分理论对所述差分影像组集合进行决策级融合,得到初始变化检测结果。
在本发明实施例中,根据所得到的每组差分影像的模糊分割矩阵Un和变化检测图CDn,利用相似度理论和模糊积分理论对所生成的六组差分影像进行决策级融合,得到初始变化检测结果。
在本发明实施例中,所述初始变化检测结果中,像元被分为变化类wc和未变化类wu这两种类别;每种类别中均包括证据冲突严重像元和证据几乎不冲突像元。
作为本发明的一个优选示例,采用Sugeno模糊积分对所生成的六组差分影像进行决策级融合:将所述六组差分影像视作为六个信源,为对其进行Sugeno模糊积分融合,需要首先计算每个信源DIn的Sugeno测度gλ(DIn)和信源DIn对每个像元p(i,j)的判断结果n=1,2,...,6。
图2示出了本发明实施例提供的遥感影像变化检测方法中步骤S104的具体实现流程。
参阅图2,所述步骤S104具体包括:
在步骤S201中,根据所述变化检测图CDn,通过相似度理论计算每组差分影像DIn的模糊测度以确定在融合过程中每组差分影像DIn的重要程度。
在本发明实施例中,根据变化检测图CDn,n=1,2,...,6,采用相似度理论计算每组差分影像DIn的模糊测度。
在这里,模糊测度特指Sugeno测度,每组差分影像DIn的Sugeno测度gλ(DIn),n=1,2,...,6,表示在融合过程中信源DIn的重要程度,即gλ(DIn)表示在融合过程中每组差分影像DIn的重要程度。
变化检测图CDn的相似度能够在很大程度上反映每组差分影像DIn的重要程度,本发明实施例中采用相似度来定义每组差分影像DIn的Sugeno测度。研究区内通常只有较少的像元发生变化,若使用检测图的所有像元来定义相似度,则会降低相似度对检测图差异的敏感性。为此,本发明实施例采用杰卡德相似度理论计算每组差分影像DIn的Sugeno测度gλ(DIn),n=1,2,...,6,即采用杰卡德距离定义不同检测图之间的相似度:
对变化检测图CDn和CDm,m,n=1,2,...,6且m≠n;用M00表示CDn和CDm的取值同时为0的像元数目,M01表示CDn的取值为0而CDm的取值为1的像元数目,M10表示CDn取值为1而CDm的取值为0的像元数目,M11表示CDn和CDm的取值同时为1的像元数目,那么根据杰卡德相似度理论计算得到CDn和CDm的相似度J(CDn,CDm)为:
然后通过计算相似度J(CDn,CDm)关于m的平均值可为变化检测图CDn定义一个总相似度J(CDn):
其中m,n=1,2,...,6且m≠n;变化检测图CDn的总相似度J(CDn)描述其与其他5组变化检测图的总一致程度,能够反映信源DIn的可靠性,本发明实施例中采用J(CDn)来定义差分影像DIn对融合结果的重要程度,令gλ(DIn)=J(CDn)。
在步骤S202中,根据所述模糊分割矩阵Un求解每组差分影像DIn对像元p(i,j)的判断结果。
在本发明实施例中,利用模糊分割矩阵来求解信源DIn对像元p(i,j)的判断结果,该判断结果包括像元被标记为变化类的概率和被标记为未变化类的概率的判断结果,即像元关于变化类或未变化类的判断结果。在这里,使用符号表示信源DIn对像元p(i,j)的判断结果,其中,表示根据信源DIn像元p(i,j)被标记为未变化类wu的概率,表示根据信源DIn像元p(i,j)被标记为变化类wc的概率。并令k∈{wu,wc},n=1,2,...,6,其中表示模糊C均值通过信源DIn得到的像元p(i,j)属于类别k的隶属度,这里k∈{wu,wc},wu表示未变化类,wc表示变化类。
在步骤S203中,根据所述重要程度gλ(DIn)和所述判断结果hn(i,j),通过模糊积分理论对所述差分影像组集合进行决策级融合,将融合后的结果作为初始变化检测结果。
在本发明实施例中,利用Sugeno积分对所生成的六组差分影像DI={DI1,DI2,...,DI6}进行决策级融合,得到初始变化检测结果。
在这里,将所述六组差分影像视作为六个信源,为便于描述,将差分影像组DIn对应描述为信源DIn,差分影像组集合对应描述为信源集合。
对任给像元p(i,j),其对应两种类型的数值:信源的重要程度和信源对像元的判断结果。Sugeno积分通过非线性策略对两种数值进行集成,并根据集成结果做出决策,具体的:
对未变化类wu,首先重新排序信源集合DI={DI1,DI2,...,DI6},使得其满足条件然后对排序后的信源集合DI={DI1,DI2,...,DI6}定义集合序列An={DI1,DI2,...,DIn},n=1,2,...,6,并通过下式计算集合An的Sugeno测度gλ(An):
其中,λ表示两个元素之间的相互作用程度,λ>-1且λ≠0,可通过求解方程得到λ。那么,Sugeno积分关于未变化类融合差分影像组集合的过程,即计算所述像元p(i,j)关于未变化类wu的模糊积分可通过下式计算得到:
其中,max表示取最大值运算,min表示取最小值运算。所述像元p(i,j)关于变化类wc的模糊积分可通过相同的计算步骤得到。
在得到像元p(i,j)关于未变化类wu和变化类wc的模糊积分和之后,通过下式对其分配类别,即将像元分为变化类或未变化类,从而得到初始变化检测结果:
其中,L(i,j)=wu表示将像元p(i,j)标记为未变化类;L(i,j)=wc表示将像元p(i,j)标记为变化类。
在这里,通过Sugeno积分能够有效地完成多源信息的融合,然而当不同信源提供的证据冲突严重时,通过Sugeno积分不能得到理想的变化检测结果。为此,需进一步处理初始变化检测结果中的证据冲突严重像元,对其重新分配类别,即对证据冲突严重像元重新进行判断以确定其是否发生了变化,将其重新划分为变化类像元或未变化类像元,以进一步提高变化检测的可靠性和变化检测结果的精度。
在步骤S105中,计算所述差分影像组集合在每个像元处的证据冲突程度;探测并定位所述初始变化检测结果中的证据冲突严重像元,通过指示克里金理论对所述证据冲突严重像元重新进行分类。
这里证据冲突程度指的是各个信源判断结果之间的不一致程度。例如,对任给像元p(i,j),可能三组信源认为其发生了变化,而另外三组信源认为其未发生变化,这是六组信源判断结果之间不一致最严重的情形,即证据冲突最严重的情形。
需要指出的是,步骤S105根据证据冲突程度将像元划分为证据冲突严重像元和证据几乎不冲突像元两部分;然后只对证据冲突严重像元重新分配类别,而对证据几乎不冲突像元,其初始变化检测结果被作为最终的变化检测结果。
图3示出了本发明实施例提供的遥感影像变化检测方法中步骤S105的具体实现流程。参阅图3,步骤S105包括:
在步骤S301中,计算任意两组差分影像在每个像元处的证据冲突程度,根据所述证据冲突程度得到所述差分影像组集合在每个像元处的总证据冲突程度。
在本发明实施例中,计算信源集合DI={DI1,DI2,...,DI6}在每个像元处的证据冲突程度。
任给像元p(i,j),首先计算两两信源之间的证据冲突程度,而后求取它们的平均值,并将所求取得到的平均值作为像元p(i,j)的总证据冲突程度。
用Kgh(i,j)表示信源DIg和信源DIh在像元p(i,j)处的证据冲突程度,g,h=1,2,...,6且g≠h;则Kgh(i,j)可通过下式计算:
在这里,表示FCM根据信源DIg得到的像元p(i,j)属于未变化类wu的隶属度;表示FCM通过信源DIg得到的像元p(i,j)属于变化类wc的隶属度;表示FCM通过信源DIh得到的像元p(i,j)属于未变化类wu的隶属度;表示FCM通过信源DIh得到的像元p(i,j)属于变化类wc的隶属度。
将计算得到的所有Kgh(i,j),g,h=1,2,...,6且g≠h,进行求平均值运算求解总证据冲突程度K(i,j):若K(i,j)的取值越大,则说明六个信源在像元p(i,j)处的证据冲突程度越高。本发明根据此性质探测并定位证据冲突严重像元。
在步骤S302中,根据所述总证据冲突程度将所述初始变化检测结果中的像元划分为两部分:证据冲突严重像元和证据几乎不冲突像元。
在本发明实施例中,通过自适应阈值技术探测并定位所述初始变化检测结果中的证据冲突严重像元,从而将初始变化检测结果中的像元划分为证据冲突严重像元和证据几乎不冲突像元两部分。
具体的,分别用FIu和FIc表示Sugeno积分融合结果,FIu表示未变化类像元的集合,FIc表示变化类像元的集合,即初始变化检测结果中的未变化类像元组成的集合FIu和变化类像元组成的集合FIc,并采用自适应阈值技术探测定位集合FIu和FIc中的证据冲突严重像元,其中集合FIu和FIc具体表示如下:
为叙述方便,用符号表示集合FIu中像元的总证据冲突程度K(i,j)组成的集合,并用和表示的均值和标准差;用符号表示集合FIc中像元的总证据冲突程度K(i,j)组成的集合,并用和表示集合的均值和标准差;将FIu中的证据冲突严重像元构成的集合记作和将FIc中的证据冲突严重像元构成的集合记作本发明实施例中采用均值和标准差来定义集合
这里Tu为正整数,且满足Tu取值越大,中的像元数目越少;Tu取值越小,中的像元数目越多。类似的,采用均值和标准差来定义集合
Tc为正整数,且满足Tc取值越大,中像元数目越少;Tc取值越小,中的像元数目越多。在这里,Tu和Tc是两个常数,优选的,在本发明实施例中Tu=2,Tc=6。
至此,初始变化检测结果被划分为三个区域:证据几乎不冲突的未变化类像元证据几乎不冲突的变化类像元和证据冲突严重像元
因证据几乎不冲突像元各个信源之间的一致程度较高,使得证据几乎不冲突像元的初始变化检测结果可信度较高,故将证据几乎不冲突像元的初始变化检测结果作为证据几乎不冲突像元的最终变化检测结果:即将集合中像元最终标记为未变化类像元,将集合中像元最终标记为变化类像元。而对集合中的像元,证据冲突比较严重,需重新判断其是否发生了变化,故此处将其类别重新标记为空,等待重新分类。
在步骤S303中,将所述的证据几乎不冲突像元的初始变化检测结果作为所述证据几乎不冲突像元的最终变化检测结果,通过指示克里金理论对所述证据冲突严重像元重新进行分类,得到所述证据冲突严重像元的最终变化检测结果。
在本发明实施例中,证据冲突严重像元包括初始变化检测结果中证据冲突严重的未变化类像元和证据冲突严重的变化类像元,其初始变化检测结果可能是不正确的,需要进一步判断其是否发生了变化,重新分配类别,从而得到证据冲突严重像元的最终变化检测结果。
具体的,利用指示克里金理论,对证据冲突严重像元重新进行分类。在这里,证据冲突严重像元指的是集合中的像元。对任给像元图4中示出了本发明实施例提供的一种利用指示克里金理论对证据冲突严重像元重新进行分类的方法的具体实现流程。
参阅图4,该流程包括:
在步骤S401中,通过指示克里金理论计算所述证据冲突严重像元属于未变化类wu的概率。
在本发明实施例中,利用指示克里金插值技术计算任给证据冲突严重像元p(i0,j0)属于未变化类wu的概率。
假设N(i0,j0)为像元p(i0,j0)的一个邻域,且N(i0,j0)满足:①N(i0,j0)以像元p(i0,j0)为中心;②则称N(i0,j0)为像元p(i0,j0)的一个克里金窗口。图5给出一个7×7的克里金窗口示意图。
为计算像元p(i0,j0)属于未变化类wu的概率,首先从未发生变化的角度根据步骤S302的结果定义指标变量Iu(i,j):通过步骤S302,一些像元被标记为未变化类,一些像元被标记为变化类,其余像元的类别为空;将已标记为未变化类像元的指标值设置为1,已标记为变化类像元的指标值设置为0;而对类别为空的像元,其有一半的可能发生了变化,一半的可能未发生变化,将其指标值设置为1/2,则Iu(i,j)具体表示如下:
然后根据所定义的指标变量Iu(i,j),利用指示克里金插值技术计算p(i0,j0)属于未变化类wu的概率
其中,N(i0,j0)表示像元p(i0,j0)的7×7的克里金窗口;为像元p(i,j)的克里金权重系数,通过求解下式给出的Iu(i,j)的克里金方程组得到:
其中,Cu(·,·)表示指标变量Iu(·,·)的协方差函数,Cu(p(i,j),p(g,h))表示像元p(i,j)和p(g,h)的协方差,Cu(p(i0,j0),p(g,h))表示像元p(i0,j0)和p(g,h)的协方差,μ为拉格朗日系数。
作为本发明的一个优选示例,计算指标变量Iu(·,·)的协方差函数时,采用简单有效的切可雪夫距离度量像元间的距离;同时假设Iu(·,·)各向同性,分别计算上、下、左、右、左上、左下、右上和右下八个方向的协方差函数,并将其平均值作为最终的协方差函数。
为避免发生概率为负的情况,通过克里金方程组得到克里金权重系数后对其进行如下处理:首先判断是否小于0,若小于0,则令否则保持其值不变;然后对大于0的进行归一化操作。
在步骤S402中,通过指示克里金理论计算所述证据冲突严重像元属于变化类wc的概率。
在本发明实施例中,利用指示克里金插值技术计算任给证据冲突严重像元p(i0,j0)属于变化类wc的概率。
与步骤S401类似,为计算像元p(i0,j0)属于变化类wc的概率,首先从发生变化的角度根据步骤S302的结果定义指标变量Ic(i,j):将已标记为变化类像元的指标值设置为1,已标记为未变化类像元的指标值设置为0;而将类别为空的像元的指标值设置为1/2。
那么,像元p(i0,j0)属于变化类wc的概率可利用指示克里金插值技术通过下式计算:
这里N(i0,j0)表示像元p(i0,j0)的7×7的克里金窗口,表示像元p(i,j)的克里金权重系数,通过求解下面的克里金方程组得到:
其中,Cc(·,·)表示指标变量Ic(·,·)的协方差函数,Cc(p(i,j),p(g,h))表示像元p(i,j)和p(g,h)的协方差,Cc(p(i0,j0),p(g,h))表示像元p(i0,j0)和p(g,h)的协方差,μ为拉格朗日系数。
在步骤S403中,通过最大概率原则对所述证据冲突严重像元进行重新分类。
在本发明实施例中,通过最大概率原则对证据冲突严重像元重新分配类别,确定该证据冲突严重像元属于变化类还是未变化类。具体的,对任给证据冲突严重像元p(i0,j0),得到其属于未变化类wu和变化类wc的概率和后,通过下式对其分配类别L(i0,j0):
其中L(i0,j0)=wu表示将像元p(i0,j0)标记为未变化类,L(i0,j0)=wc表示将像元p(i0,j0)标记为变化类。
通过步骤S104得到了初始的变化检测结果,初始变化检测结果包括证据几乎不冲突的像元和证据冲突严重像元:由于证据几乎不冲突像元的初始变化检测结果可信度较高,将证据几乎不冲突像元的初始变化检测结果作为最终的变化检测结果;步骤S105利用指示克里金理论对证据冲突严重像元重新分配类别,从而得到证据冲突严重像元的最终变化检测结果。至此,实现了变化检测任务。
为了说明本发明实施例提供的遥感影像变化检测方法的变化检测效果,图6(a)给出了T1时刻遥感影像的第4波段,图6(b)给出了T2时刻遥感影像的第4波段,图6(c)给出了T1时刻和T2时刻遥感影像的变化参考图,两期影像产生变化的主要原因是火灾。图7(a)给出了模糊C均值根据差分影像DI1得到的变化检测图,图7(b)给出了模糊C均值根据差分影像DI2得到的变化检测图,图7(c)给出了模糊C均值根据差分影像DI3得到的变化检测图,图7(d)给出了模糊C均值根据差分影像DI4得到的变化检测图,图7(e)给出了模糊C均值根据差分影像DI5得到的变化检测图,图7(f)给出了模糊C均值根据差分影像DI6得到的变化检测图,图7(g)给出了融合技术投票法的变化检测图,图7(h)给出了融合技术模糊积分法的变化检测图,图7(i)给出了本发明实施例提供的遥感影像变化检测方法的变化检测图。
通过比较图7的变化检测图和图6(c)的变化参考图,可以发现本发明实施例得到了最接近真实变化的变化检测结果。
表1进一步给出了上述九种变化检测技术的定量分析指标。其中A行表示模糊C均值根据差分影像DI1得到的变化检测结果;B行表示模糊C均值根据差分影像DI2得到的变化检测结果;C行表示模糊C均值根据差分影像DI3得到的变化检测结果;D行表示模糊C均值根据差分影像DI4得到的变化检测结果;E行表示模糊C均值根据差分影像DI5得到的变化检测结果;F行表示模糊C均值根据差分影像DI6得到的变化检测结果;G行表示投票法的变化检测结果;H行表示模糊积分法的变化检测结果;I行表示本发明实施例提供的遥感影像变化检测方法的变化检测结果。所述指标包括:1.虚检错误(将未变化类像元检测为变化类像元的像元个数);2.漏检错误(将变化类像元检测为未变化类像元的像元个数);3.总体错误(虚检错误与漏检错误之和);4.Kappa系数。
从表1的比对中可得,本发明实施例提供的变化检测结果明显优于其他八种变化检测技术:与前面的六种基于单差分影像的方法相比,本发明实施例的变化检测结果总体错误降低约1600—4500像元,Kappa系数提高约2%—6%;与后面两种经典融合算法相比,本发明实施例提供的变化检测结果总体错误降低约1700像元,Kappa系数提高约2%。因此,本发明实施例提供的遥感影像变化检测方法能够取得较优的变化检测结果。
表1
方法 | 虚检错误 | 漏检错误 | 总体错误 | Kappa系数 |
A | 4359 | 3217 | 7576 | 0.9114 |
B | 3865 | 3151 | 7016 | 0.9177 |
C | 2275 | 3248 | 5523 | 0.9332 |
D | 626 | 7771 | 8397 | 0.8937 |
E | 584 | 7348 | 7932 | 0.9000 |
F | 3297 | 2535 | 5832 | 0.9316 |
G | 2452 | 3178 | 5630 | 0.9330 |
H | 1458 | 4210 | 5668 | 0.9315 |
I | 1869 | 2060 | 3929 | 0.9535 |
在本发明实施例中,首先从距离矩阵的角度,提出一种生成信息互补差分影像组集合的方案,充分提取两期遥感影像的变化信息;其次通过相似度理论和模糊积分理论实现对所述差分影像组集合的决策级融合;然后利用指示克里金技术对所述融合结果进行优化,即重新分类融合结果中证据冲突严重像元,从而实现了提高变化检测的可靠性和变化检测结果的精度。
应理解,上述实施例中各步骤的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本发明实施例的实施过程构成任何限定。
对应于上文实施例所述的一种遥感影像变化检测方法,图8示出了本发明实施例提供的一种遥感影像变化检测装置的示意图,为了便于说明,仅示出了与本发明实施例相关的部分。
参照图8,该装置包括:
预处理单元81,用于获取两期遥感影像,并对所述遥感影像进行预处理;
差分影像组集合生成单元82,用于根据预处理后的两期遥感影像生成信息互补的差分影像组集合;
模糊C均值聚类计算单元83,用于对所述差分影像组集合进行模糊C均值聚类计算,得到每组差分影像的模糊分割矩阵和变化检测图;
决策级融合单元84,用于根据所述的每组差分影像的模糊分割矩阵和变化检测图,通过相似度理论和模糊积分理论对所述差分影像组集合进行决策级融合,得到初始变化检测结果;
证据冲突严重像元重分类处理单元85,用于计算所述差分影像组集合在每个像元处的证据冲突程度;探测并定位所述初始变化检测结果中的证据冲突严重像元,通过指示克里金理论对所述证据冲突严重像元重新进行分类。
具体的,所述差分影像组集合由六组信息互补的差分影像组成;
其中,所述差分影像组集合生成单元82具体包括:
第一差分影像组计算单元,用于根据预处理后的两期遥感影像,通过欧氏距离公式计算得到第一组差分影像;
第二差分影像组计算单元,用于根据预处理后的两期遥感影像,通过闵可夫斯基距离公式计算得到第二组差分影像;
第三差分影像组计算单元,用于根据预处理后的两期遥感影像,通过马氏距离公式计算得到第三组差分影像;
第四差分影像组计算单元,用于根据预处理后的两期遥感影像,通过相关性距离公式计算得到第四组差分影像;
第五差分影像组计算单元,用于根据预处理后的两期遥感影像,通过主成分分析理论生成第五组差分影像;
第六差分影像组计算单元,用于根据预处理后的两期遥感影像,通过基于波普曲线形状的波普梯度差异生成第六组差分影像。
优选的,所述遥感影像变化检测装置还包括:
标准化处理单元,用于根据预设公式将每组差分影像的取值进行标准化使得其取值在预设取值区间范围之内。
具体的,所述决策级融合单元84包括:
重要程度确认子单元,用于根据所述变化检测图,通过相似度理论计算每组差分影像的模糊测度以确定在融合过程中每组差分影像的重要程度;
判断结果求解子单元,用于根据所述模糊分割矩阵求解每组差分影像对像元的判断结果;
决策级融合子单元,用于根据所述重要程度和所述判断结果,通过模糊积分理论对所述差分影像组集合进行决策级融合,将融合后的结果作为初始变化检测结果。
具体的,所述证据冲突严重像元重分类处理单元85包括:
证据冲突程度计算子单元,用于计算任意两组差分影像在每个像元处的证据冲突程度,根据所述证据冲突程度得到所述差分影像组集合在每个像元处的总证据冲突程度;
证据冲突严重像元探测子单元,用于根据所述总证据冲突程度,通过自适应阈值技术将所述初始变化检测结果划分为两部分:证据几乎不冲突像元和证据冲突严重像元;
证据冲突严重像元重分类处理子单元,用于通过指示克里金技术对所述证据冲突严重像元重新进行分类。
其中,所述证据冲突严重像元重分类处理子单元具体用于:
通过指示克里金理论计算所述证据冲突严重像元属于未变化类的概率;
通过指示克里金理论计算所述证据冲突严重像元属于变化类的概率;
通过最大概率原则对所述证据冲突严重像元进行重新分类。
图9是本发明实施例提供的一种终端的示意图。如图9所示,该实施例的终端9包括:处理器90、存储器91以及存储在所述存储器91中并可在所述处理器90上运行的计算机程序92。所述处理器90执行所述计算机程序92时实现上述遥感影像变化检测方法实施例中的各个步骤,例如图1所示的步骤S101至S105。或者,所述处理器90执行所述计算机程序92时实现上述各装置实施例中各单元的功能,例如图8所示模块81至85的功能。
示例性的,所述计算机程序92可以被分割成一个或多个单元,所述一个或者多个单元被存储在所述存储器91中,并由所述处理器90执行,以完成本发明。所述一个或多个单元可以是能够完成特定功能的一系列计算机程序指令段,该指令段用于描述所述计算机程序92在所述终端9中的执行过程。例如,所述计算机程序92可以被分割成预处理单元81、差分影像组集合生成单元82、模糊C均值聚类计算单元83、决策级融合单元84、证据冲突严重像元重分类处理单元85,各单元具体功能如下:
预处理单元81,用于获取两期遥感影像,并对所述遥感影像进行预处理;
差分影像组集合生成单元82,用于根据预处理后的两期遥感影像生成信息互补的差分影像组集合;
模糊C均值聚类计算单元83,用于对所述差分影像组集合进行模糊C均值聚类计算,得到每组差分影像的模糊分割矩阵和变化检测图;
决策级融合单元84,用于根据所述的每组差分影像的模糊分割矩阵和变化检测图,通过相似度理论和模糊积分理论对所述差分影像组集合进行决策级融合,得到初始变化检测结果;
证据冲突严重像元重分类处理单元85,用于计算所述差分影像组集合在每个像元处的证据冲突程度;探测并定位所述初始变化检测结果中的证据冲突严重像元,通过指示克里金理论对所述证据冲突严重像元重新进行分类。
所述终端9可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所述终端9可包括,但不仅限于,处理器90、存储器91。本领域技术人员可以理解,图9仅仅是终端9的示例,并不构成对终端9的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件,例如所述终端还可以包括输入输出设备、网络接入设备、总线等。
所述处理器90可以是中央处理单元(Central Processing Unit,CPU),还可以是其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现成可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
所述存储器91可以是所述终端9的内部存储单元,例如终端9的硬盘或内存。所述存储器91也可以是所述终端9的外部存储设备,例如所述终端9上配备的插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,所述存储器91还可以既包括所述终端9的内部存储单元也包括外部存储设备。所述存储器91用于存储所述计算机程序以及所述终端所需的其他程序和数据。所述存储器91还可以用于暂时地存储已经输出或者将要输出的数据。
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各功能单元、模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能单元、模块完成,即将所述装置的内部结构划分成不同的功能单元或模块,以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。另外,各功能单元、模块的具体名称也只是为了便于相互区分,并不用于限制本申请的保护范围。上述系统中单元、模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述或记载的部分,可以参见其他实施例的相关描述。
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
在本发明所提供的实施例中,应该理解到,所揭露的装置/终端设备和方法,可以通过其他的方式实现。例如,以上所描述的装置/终端设备实施例仅仅是示意性的,例如,所述模块或单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口,装置或单元的间接耦合或通讯连接,可以是电性,机械或其他的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实现上述实施例方法中的全部或部分流程,也可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一计算机可读存储介质中,该计算机程序在被处理器执行时,可实现上述各个方法实施例的步骤。其中,所述计算机程序包括计算机程序代码,所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括:能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是,所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减,例如在某些司法管辖区,根据立法和专利实践,计算机可读介质不包括电载波信号和电信信号。
以上所述实施例仅用于说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围,均应包含在本发明的保护范围之内。
Claims (10)
1.一种遥感影像变化检测方法,其特征在于,包括:
获取两期遥感影像,并对所述遥感影像进行预处理;
根据预处理后的两期遥感影像生成信息互补的差分影像组集合;
对所述差分影像组集合进行模糊C均值聚类计算,得到每组差分影像的模糊分割矩阵和变化检测图;
根据所述的每组差分影像的模糊分割矩阵和变化检测图,通过相似度理论和模糊积分理论对所述差分影像组集合进行决策级融合,得到初始变化检测结果;
计算所述差分影像组集合在每个像元处的证据冲突程度;探测并定位所述初始变化检测结果中的证据冲突严重像元,通过指示克里金理论对所述证据冲突严重像元重新进行分类。
2.如权利要求1所述的遥感影像变化检测方法,其特征在于,所述差分影像组集合由六组信息互补的差分影像组成;
所述根据预处理后的两期遥感影像生成信息互补的差分影像组集合的步骤,具体为:
根据预处理后的两期遥感影像,通过欧氏距离公式计算得到第一组差分影像;
根据预处理后的两期遥感影像,通过闵可夫斯基距离公式计算得到第二组差分影像;
根据预处理后的两期遥感影像,通过马氏距离公式计算得到第三组差分影像;
根据预处理后的两期遥感影像,通过相关性距离公式计算得到第四组差分影像;
根据预处理后的两期遥感影像,通过主成分分析理论生成第五组差分影像;
根据预处理后的两期遥感影像,通过基于波普曲线形状的波普梯度差异生成第六组差分影像。
3.如权利要求1或2所述的遥感影像变化检测方法,其特征在于,所述根据所述的每组差分影像的模糊分割矩阵和变化检测图,通过相似度理论和模糊积分理论对所述差分影像组集合进行决策级融合,得到初始变化检测结果的步骤,包括:
根据所述变化检测图,通过相似度理论计算每组差分影像的模糊测度以确定在融合过程中每组差分影像的重要程度;
根据所述模糊分割矩阵求解每组差分影像对像元的判断结果;
根据所述重要程度和所述判断结果,通过模糊积分理论对所述差分影像组集合进行决策级融合,将融合后的结果作为初始变化检测结果。
4.如权利要求1所述的遥感影像变化检测方法,其特征在于,所述计算所述差分影像组集合在每个像元处的证据冲突程度;探测并定位所述初始变化检测结果中的证据冲突严重像元,通过指示克里金理论对所述的证据冲突严重像元重新进行分类的步骤,包括:
计算任意两组差分影像在每个像元处的证据冲突程度,根据所述证据冲突程度得到所述差分影像组集合在每个像元处的总证据冲突程度;
根据所述总证据冲突程度,通过自适应阈值技术将所述初始变化检测结果中的像元划分为两部分:证据几乎不冲突像元和证据冲突严重像元;
将所述证据几乎不冲突像元的初始变化检测结果作为所述证据几乎不冲突像元的最终变化检测结果,通过指示克里金理论对所述证据冲突严重像元重新进行分类,得到所述证据冲突严重像元的最终变化检测结果。
5.如权利要求1所述的遥感影像变化检测方法,其特征在于,所述的通过指示克里金理论对所述证据冲突严重像元重新进行分类的步骤,包括:
通过指示克里金理论计算所述证据冲突严重像元属于未变化类的概率;
通过指示克里金理论计算所述证据冲突严重像元属于变化类的概率;
通过最大概率原则对所述证据冲突严重像元进行重新分类。
6.一种遥感影像变化检测装置,其特征在于,包括:
预处理单元,用于获取两期遥感影像,并对所述遥感影像进行预处理;
差分影像组集合生成单元,用于根据预处理后的两期遥感影像生成信息互补的差分影像组集合;
模糊C均值聚类计算单元,用于对所述差分影像组集合进行模糊C均值聚类计算,得到每组差分影像的模糊分割矩阵和变化检测图;
决策级融合单元,用于根据所述的每组差分影像的模糊分割矩阵和变化检测图,通过相似度理论和模糊积分理论对所述差分影像组集合进行决策级融合,得到初始变化检测结果;
证据冲突严重像元重分类处理单元,用于计算所述差分影像组集合在每个像元处的证据冲突程度;探测并定位所述初始变化检测结果中的证据冲突严重像元,通过指示克里金理论对所述证据冲突严重像元重新进行分类。
7.如权利要求6所述的遥感影像变化检测装置,其特征在于,所述差分影像组集合由六组信息互补的差分影像组成;
所述差分影像组集合生成单元具体包括:
第一差分影像组计算单元,用于根据预处理后的两期遥感影像,通过欧氏距离公式计算得到第一组差分影像;
第二差分影像组计算单元,用于根据预处理后的两期遥感影像,通过闵可夫斯基距离公式计算得到第二组差分影像;
第三差分影像组计算单元,用于根据预处理后的两期遥感影像,通过马氏距离公式计算得到第三组差分影像;
第四差分影像组计算单元,用于根据预处理后的两期遥感影像,通过相关性距离公式计算得到第四组差分影像;
第五差分影像组计算单元,用于根据预处理后的两期遥感影像,通过主成分分析理论生成第五组差分影像;
第六差分影像组计算单元,用于根据预处理后的两期遥感影像,通过基于波普曲线形状的波普梯度差异生成第六组差分影像。
8.如权利要求6或7所述的遥感影像变化检测装置,其特征在于,所述决策级融合单元包括:
重要程度确认子单元,用于根据所述变化检测图,通过相似度理论计算每组差分影像的模糊测度以确定在融合过程中每组差分影像的重要程度;
判断结果求解子单元,用于根据所述模糊分割矩阵求解每组差分影像对像元的判断结果;
决策级融合子单元,用于根据所述重要程度和所述判断结果,通过模糊积分理论对所述差分影像组集合进行决策级融合,将融合后的结果作为初始变化检测结果。
9.一种终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至5任一项所述遥感影像变化检测方法的步骤。
10.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至5任一项所述遥感影像变化检测方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710655794.9A CN109389571B (zh) | 2017-08-03 | 2017-08-03 | 一种遥感影像变化检测方法、装置及终端 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710655794.9A CN109389571B (zh) | 2017-08-03 | 2017-08-03 | 一种遥感影像变化检测方法、装置及终端 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109389571A true CN109389571A (zh) | 2019-02-26 |
CN109389571B CN109389571B (zh) | 2020-03-17 |
Family
ID=65412228
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710655794.9A Active CN109389571B (zh) | 2017-08-03 | 2017-08-03 | 一种遥感影像变化检测方法、装置及终端 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109389571B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112016400A (zh) * | 2020-08-04 | 2020-12-01 | 香港理工大学深圳研究院 | 一种基于深度学习的单类目标检测方法、设备及存储介质 |
CN112767376A (zh) * | 2021-01-27 | 2021-05-07 | 三峡大学 | 引力模型优化条件随机场的多尺度融合图像变化检测方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050007466A1 (en) * | 2003-06-27 | 2005-01-13 | Katsuyuki Tsukui | Image sensing apparatus and control method thereof |
US20130336540A1 (en) * | 2012-06-14 | 2013-12-19 | Hitachi, Ltd. | Decomposition apparatus and method for refining composition of mixed pixels in remote sensing images |
CN103984946A (zh) * | 2014-05-23 | 2014-08-13 | 北京联合大学 | 一种基于K-means的高分辨率遥感地图道路提取方法 |
CN104867150A (zh) * | 2015-05-22 | 2015-08-26 | 武汉大学 | 遥感影像模糊聚类的波段修正变化检测方法及系统 |
CN105551031A (zh) * | 2015-12-10 | 2016-05-04 | 河海大学 | 基于fcm和证据理论的多时相遥感影像变化检测方法 |
CN106846246A (zh) * | 2016-12-16 | 2017-06-13 | 河海大学 | 一种基于对象的遥感影像超分辨率制图方法 |
CN106897679A (zh) * | 2017-02-13 | 2017-06-27 | 长江水利委员会长江科学院 | 一种基于改进模糊c均值聚类的语义变化检测方法及系统 |
-
2017
- 2017-08-03 CN CN201710655794.9A patent/CN109389571B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050007466A1 (en) * | 2003-06-27 | 2005-01-13 | Katsuyuki Tsukui | Image sensing apparatus and control method thereof |
US20130336540A1 (en) * | 2012-06-14 | 2013-12-19 | Hitachi, Ltd. | Decomposition apparatus and method for refining composition of mixed pixels in remote sensing images |
CN103984946A (zh) * | 2014-05-23 | 2014-08-13 | 北京联合大学 | 一种基于K-means的高分辨率遥感地图道路提取方法 |
CN104867150A (zh) * | 2015-05-22 | 2015-08-26 | 武汉大学 | 遥感影像模糊聚类的波段修正变化检测方法及系统 |
CN105551031A (zh) * | 2015-12-10 | 2016-05-04 | 河海大学 | 基于fcm和证据理论的多时相遥感影像变化检测方法 |
CN106846246A (zh) * | 2016-12-16 | 2017-06-13 | 河海大学 | 一种基于对象的遥感影像超分辨率制图方法 |
CN106897679A (zh) * | 2017-02-13 | 2017-06-27 | 长江水利委员会长江科学院 | 一种基于改进模糊c均值聚类的语义变化检测方法及系统 |
Non-Patent Citations (1)
Title |
---|
邓冰: "遥感影像信息度量方法研究", 《中国博士学位论文全文数据库(信息科技辑)》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112016400A (zh) * | 2020-08-04 | 2020-12-01 | 香港理工大学深圳研究院 | 一种基于深度学习的单类目标检测方法、设备及存储介质 |
CN112767376A (zh) * | 2021-01-27 | 2021-05-07 | 三峡大学 | 引力模型优化条件随机场的多尺度融合图像变化检测方法 |
CN112767376B (zh) * | 2021-01-27 | 2023-07-11 | 三峡大学 | 引力模型优化条件随机场的多尺度融合图像变化检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109389571B (zh) | 2020-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109871886B (zh) | 基于谱聚类的异常点比例优化方法、装置及计算机设备 | |
EP3077960B1 (en) | A method and system for computing distance measures on a quantum computer | |
De Cáceres et al. | The management of vegetation classifications with fuzzy clustering | |
Li et al. | An extension to rough c-means clustering based on decision-theoretic rough sets model | |
Mok et al. | A robust adaptive clustering analysis method for automatic identification of clusters | |
WO2021045818A1 (en) | Protecting machine learning models from privacy attacks | |
CN106203478A (zh) | 一种用于智能电表大数据的负荷曲线聚类方法 | |
BR102014023780B1 (pt) | Método para classificação multiclasse em cenários abertos e usos do mesmo | |
CN110211126A (zh) | 基于直觉模糊c均值聚类的图像分割方法 | |
CN108764726B (zh) | 根据规则对请求进行决策的方法及装置 | |
CN106326913A (zh) | 一种洗钱账户的确定方法及装置 | |
CN108206813B (zh) | 基于k均值聚类算法的安全审计方法、装置及服务器 | |
CN115130643A (zh) | 具有异质性的数据集的图形神经网络 | |
CN109672554A (zh) | 确定故障根因的方法及装置 | |
CN109389571A (zh) | 一种遥感影像变化检测方法、装置及终端 | |
CA3033201A1 (en) | Large scale social graph segmentation | |
CN109937421B (zh) | 用于预测特定项目所属类别的二类别分类方法以及利用该方法的计算设备 | |
Ren et al. | Integrated defense for resilient graph matching | |
CN106022359A (zh) | 基于有序信息熵的模糊熵空间聚类分析方法 | |
CN107480426A (zh) | 自迭代病历档案聚类分析系统 | |
CN110110628A (zh) | 一种频率综合器劣化的检测方法及检测设备 | |
CN113515519A (zh) | 图结构估计模型的训练方法、装置、设备及存储介质 | |
Xie et al. | Automatic neuron tracing in volumetric microscopy images with anisotropic path searching | |
Ramathilagam et al. | Extended fuzzy c-means: an analyzing data clustering problems | |
Nascimento et al. | Unsupervised initialization of archetypal analysis and proportional membership fuzzy clustering |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |