专利名称：基于量子免疫克隆的sar图像变化检测方法
技术领域：
本发明属于图像处理领域，具体地说是一种划分聚类技术，可用于多时相SAR图 像变化区域的检测问题。
背景技术：
随着卫星遥感技术的发展，遥感数据飞速增加。其中的SAR数据由于不受大气环 境和和云层遮挡的影响成为一种重要的遥感数据。如何准确快速自动地从不同时相的SAR 图像中找出显著变化的区域则具有十分重要的意义。在SAR图像中寻找“非变化”和“变化”两分类的研究目前还处于初步阶段，大致 有两个路线，一条路线是分类后比较方法，也称后分类比较法，该方法先对两个时相的图像 进行独立分类，再对两幅分类图像进行逐象素的比较，最后得到变化检测图；另外一条路线 是差异图分类方法，该方法先对两个时相的图像进行逐象素的比较，比如逐象素点的差值、 比值、CVA等，再对比较得出的差异图像进行进一步地多种变换、概率分布等处理以达到两 分类，最后得到变化检测图。后分类比较法能够减少由于数据获取平台和环境的不同引起 的伪变化信息，不需要数据的辐射校正等复杂的预处理，目前较多的研究是沿着差异分类 模型这个路线进行的。差异图分类方法简单直观，得到的变化细节较为显著。构造差异影 像得到的结果与实际的变化和非变化趋势大体一致。按照策略现有的变化检测方法可以归结为七类算术运算法、变换法、分类法、高 级模型法、GIS方法、视觉分析法和其它方法。其中，在算术运算法中，应用比较广泛的是无 监督的分割方法，通常又叫做聚类方法。无监督分割方法一般可以分为两类层次聚类和划 分聚类，其中划分聚类通过最小化特定准则将数据集划分到不同的类属中，因此这类方法 可以看作是最优化问题，同时，图像的变化检测问题也可以视为组合优化问题，因此可以用 划分聚类方法来处理图像的变化检测问题。但是已有的优化方法对于处理优化问题的时候 往往耗时很长，并且在搜索过程中容易陷入局部最优，同时对于复杂图像的变化检测问题 往往会存在边缘定位不够准确的缺点，这样势必会影响到图像变化检测的区域一致性与边 缘保持的性能。
发明内容
本发明的目的在于克服上述已有技术的不足，提出一种基于量子免疫克隆聚类算 法的变化检测方法，以快速、有效地搜索到最优聚类中心，准确定位边缘，提高图像区域一 致性和边缘保持性能。实现本发明的技术方案是用聚类来实现变化检测，把变化区域的检测问题看作组 合优化问题，用量子免疫克隆算法计算搜索，使亲合度函数最大化的序列组合作为变化检 测结果，进而得到最结果。具体实现步骤如下(1)对输入的两时相图像分别采用中值滤波，得到滤波后的图像I1和I2 ；(2)对滤波后的图像I1和I2求对数比差异影像13，并将该I3的灰度值作为聚类数
4据集；(3)设置抗体规模N = 20、类别数k = 2和停机条件，随机产生初始量子抗
体2 (O =“2 “‘ ？ ’作为聚类数据集的初始聚类中心，Q(t)中的
P i P 2 … P m Jf，
^ = 2-m)都以等概率w々初始化，所述的停机条件包括最大亲合度值改变量的阈值 ε = 10e-5以及连续无法改进次数η = 5 ；(4)通过观察量子抗体Q(t)的状态，随机产生一个W，l]数，若随机产生的数大 于|<|2 ,则在相应的二进制位上取1，否则取0，以得到观测后的解~) = { …乂}，每个 《C/= 1,2,…，均表示长度为m的二进制串(X1XfXm)；(5)计算观测后的二进制抗体P (t)与聚类数据集的亲合度函数fk，保留当前种群 中的最优子抗体qbest;(6)将量子抗体Q(t)进行量子旋转门变异操作，得到量子种群Qm(t)；(7)将量子种群Qm(t)进行量子全干扰交叉重组操作，得到聚类中心Q。(t)；(8)将聚类中心Q。(t)观测成为新的二进制抗体P。(t)，计算ρ。(t)中每个子抗体 与聚类数据集的亲合度函数值f。；(9)对p。(t)进行选择操作，得到子代抗体ρ (t+1)；(10)判断子代抗体是否满足停机条件，如果满足该条件就将子代抗体中亲合度 最高的抗体对应的图像类属划分作为输出结果，否则返回步骤(4)，循环执行步骤(4) (10)，直到满足停机条件。本发明与现有技术相比具有如下优点1.在图像数据聚类过程中快速且有效地搜索到最优聚类中心，防止在进化过程中 陷入局部最优解。本发明由于采用了量子免疫克隆聚类算法来实现SAR图像变化检测，在操作中利 用量子编码的叠加性来构造抗体，使得作用在量子编码抗体上的操作具有高效的并行性， 为防止盲目的搜索，利用当前最优抗体的信息来控制变异，使种群以大概率向着优良模式 进化来加速收敛，并且有效地提高了搜索速度，使得该操作的时间复杂度降低。随着问题的 复杂求解能力不尽人意，在各个子群体间采用量子重组操作增强信息交流，在各个子群体 内部采用量子旋转门对抗体进行进化，并动态调整旋转角度，在全局搜索的同时兼顾局部， 有效防止了在进化过程中陷入局部最优解。2.图像变化检测效果好本发明由于采用了量子编码，量子重组，量子变异操作和精英选择策略操作，因而 具有比现有技术更好的图像变化检测效果。仿真实验表明，针对几幅变化检测SAR图像，其 变化检测结果的区域一致性，边缘准确性好，同时具有更低的变化检测错误率。
图1是本发明的SAR图像变化检测流程图；图2是现有量子旋转门的构造示意图；图3是用本发明方法对Bern城市水灾变化检测结果示意图4是用本发明方法对Ottawa地区水灾变化检测结果示意图。
具体实施例方式参照图1，本发明的变化检测过程如下步骤1，对待检测的两时相图像进行简化处理。这里选取形态学中最常用的工具之一中值滤波器，窗口大小为3x3，用该滤波器对 输入的两时相图像进行滤波处理，得到滤波后的图像I1和12，滤波的目的是去掉小的噪声 干扰以及对感知不重要的细节，对图像起到平滑作用，与经典的图像简化工具，如低通或中 通滤波器相比，中值滤波器的优势在于简化图像而不造成图像模糊或改变图像轮廓。步骤2，按如下步骤对滤波后的图像I1和I2求对数比差异影像I3，并将得到的I3 的灰度值作为聚类数据集。(2a)求待变化检测两时相图像的对数比差异影像Id = Ilog(I^l)-Iogd1+!)其中I' ！ = (I1 (i, j),l ^ i ^ 1,1 ^ j ^ J}, I2 = {I2(i, j),l ^ i ^ I, 1 ^ j ^ J};(2b)对差异影像进行归一化处理，得到对数比差异影像图I3 = 255* (I11-Imin) / (Imax-Imin)其中Imax = max (Id)，表示Id中最大灰度值，Imin = min (Id)，表示Id中最小灰度值。步骤3，设置初始化参数，生成初始量子抗体。(3a)按照量子免疫克隆算法设置抗体规模N = 20、类别数k = 2，同时给定停机 条件，该停机条件包括最大亲合度值改变量的阈值ε =10e_5以及连续无法改进次数η = 5 ；(3b)初始量子抗体Q(t)按如下定义随机产生，作为聚类数据集的初始聚类中心一个具有m个量子比特位的抗体Q(t)可以描述为Q (O =其中α，β表示一个量子比特位，是用一对复数归|Ψ>= α |0>+β | 1 >来定 义，这个复数表示一个量子比特位的状态可以取0或1，该量子比特位处于状态0和状态1 的概率分别是I α I2, I β I2，并且QiI2+! ^iI2 = Ki = 1,2,…，m)，这种表示方法能够表 征任意的线性叠加态。在本发明中，Q(t)中的= — 都以等概率▲初始化，意味着所有 可能的线性叠加态以相同的概率出现。步骤4，将初始量子抗体Q (t)观测成为二进制抗体ρ (t)。通过观察初始量子抗体Q (t)的状态，产生一组普通二进制解ρ (t)，其中在第t代 中P(t)可以描述为= K, '■··,<}其中<(/.= 1,2,…,《)是长度为m的二进制串(X1X2…xm)，它是由量子比特幅度|α,'|2或|A'|2 (^1,2,---,肌)得到的。产生二进制抗体(X1XfXm)的具体过程为随机产生一个
数，若随机产生的 数大于|<|2;则相应的二进制位上取1，否则取0。步骤5，计算每个观测后的二进制抗体ρ (t)与聚类数据集的亲合度函数fk，保留 当前群体中的最优抗体qbest。抗体亲合度函数定义为 其中，i表示类别，j表示样本点，μ 表示一个像素点样本j属于各个类别i的隶 属度，m表示模糊指数，d(Xj, Vi)表示第Vi个聚类中心到第Xj个像素样本之间的欧几里德距离。步骤6，将初始量子抗体Q(t)进行量子旋转门变异操作，得到新的量子种群 Qm⑴。在本发明中，各个状态间的转移是通过量子门变换矩阵实现的，在0、1编码的问 题中，设计了下面这种量子门变换矩阵来加速进化求优
「cos(0) -sin(0)"jU(0)=J,
sin(0) cos((9) _U( θ )表示量子旋转门变换矩阵，旋转变异的角度θ可由表量子旋转门变异角θ 变化表得到，用量子旋转门的旋转角度表征了量子抗体中的优化方向，进而方便的在变异 过程中加入最优抗体的信息，加快算法收敛。表1量子旋转门变异角θ变化表 表中Xi为当前抗体的第i位，best,为当前的最优抗体的第i位，均为二值编码， f(x)为亲合度函数，Δ Qi为旋转角度的大小，用以控制算法收敛的速度，S(CiJi)为旋转 角度的方向，用于保证算法的收敛。参照图2，本发明用一个二维坐标系来表示量子旋转门的结构，其中横坐标为α表示当前解取零的概率，纵坐标为βi，表示当前解取ι的概率，单位圆表示I ai|2+| β」2 = l(i = 1,2,…，m)的所包含的区域，当Xi = 0，besti = 1，f(x)彡f(best)时，为使当前 解收敛到一个具有更高亲合度的抗体，应增大当前解取0的概率，即要使I α」2变大，如果 (a,, ^j)在第一、三象限，θ应向顺时针方向旋转；如果(Cii, β,)在第二、四象限，θ应 向逆时针方向旋转。步骤7，将新的量子种群Qm (t)进行量子全干扰交叉重组操作，得到Q。(t)，并将其 作为新的聚类中心。本发明通过使用量子的相干特性构造了一种新的重组操作_ “全干扰交叉”，量子 抗体在实行这种重组操作之前的状态如表2量子重组操作前的抗体所示，实行量子重组后 各个抗体之间的状态如表3量子重组操作后的抗体所示。该量子重组操作采用对角线交叉的方式，让处于对角线上的抗体经过这种量子重 组方式之后处于同一行上面，这样保证了种群中的所有抗体均参与重组，这种量子重组可 以充分利用种群中的尽可能多的抗体的信息，从而改进普通交叉的局部性与片面性，在种 群进化出现早熟时，这种量子重组方法能够产生新的抗体，可以给进化过程注入新的动力。表2量子重组操作前的抗体 表3量子重组操作后的抗体 步骤8，将新的量子种群Q。(t)观测成为二进制抗体P。(t)，计算每个抗体与聚类数 据集的亲合度函数值f。；步骤9，对p。(t)进行选择操作，得到子代抗体ρ (t+1)。本发明采用的是精英选择策略，具体为，即在进化过程中，如果某一代中的最优解 的亲合度函数值优于当前代最优解的亲合度函数值，则当前代最优解就被该最优解所代 替，否则保留当前代最优解，从而可以保证在进化过程中每一代的最优解都不会丢失。选择操作的目的是选择出亲合度较高的抗体，这些抗体对应的差异图像的变化检 测的结果要好于亲合度低的抗体对应的差异图像的变化检测结果，采用精英选择策略可以保证某一代的最优解在整个进化过程中可以毫发无损地被保留下来。步骤10，输出图像类属划分结果的条件判断。按照输出差异影像图的变化检测结果时最佳的抗体亲合度改变量至少连续迭代 η次不变的原则，判断停机条件，如果第t代与t+Ι代最大亲合度值之差不大于阈值ε = lOe-5，则η = η+1，否则，η不变，如此反复迭代，直到满足η次不大于已设定的阈值ε = lOe-5，就将该抗体中亲合度最高的抗体对应的图像类属划分作为输出结果，否则返回步骤 (4)，循环执行步骤(4) (10)，直到满足输出类属划分结果的停机条件为止。通过以上十个步骤的操作，对于输入的待变化检测图像均可以输出一个最佳的变 化检测结果。本发明的效果可以通过以下仿真进一步说明1.仿真数据第一组的真实SAR数据是1999年4月和1999年5月由ERS 2所携带的SAR在瑞 士 Bern城区附近获得，该组数据变化部分是在Thim市、Bern市和Bern机场段Aare河的 洪水引起的，实验中所用图像大小为301X301像素，256灰度级，其中变化像元数为1155， 未变化像素元数为89446。第二组的真实SAR数据是1997年5月和1997年8月由Radarsat所携带的SAR 在加拿大Ottawa城区附近获得，实验中所用图像大小为290x350像素，256灰度级，其中变 化像元数为11952，未变化像素元数为82871。2.仿真内容及结果(1)仿真内容把以上仿真数据的第一时相原图和第二时相原图作为仿真对象，用本发明方法检 测第二时相原图相对于第一时相原图的变化情况。(2)仿真结果·对第一组真实SAR数据及实验结果如图3所示，其中图3 (a)为真实SAR第一时相 的原图，图3(b)为真实SAR第二时相的原图，图3(c)为两时相的变化检测参考图，图7(d) 为本发明的变化检测结果图，图3(e)为对比算法K均值算法变化检测结果图，图3(f)为对 比算法遗传算法变化检测结果图。对第二组真实SAR数据及实验结果如图4所示，其中图4 (a)为真实SAR第一时相 的原图，图4(b)为真实SAR第二时相的原图，图4(c)为两时相的变化检测参考图，图7(d) 为本发明的变化检测结果图，图4(e)为对比算法K均值算法变化检测结果图，图4(f)为对 比算法遗传算法变化检测结果图。表4给出了本发明实验结果数据。表4实验结果数据 3.仿真结果分析从图3和图4可以看出，本发明方法对于复杂图像的变化检测问题边缘定位准确， 图像变化检测的区域一致性与边缘保持的性能相对于K均值算法和遗传算法都有了很大 的提高。从表4可以得到，本发明对真实SAR数据，相对于K均值算法和遗传算法来说，误 检率和漏检率，错误检测率都处于较低的水平，总的来说，变化检测精度高。
权利要求
一种基于量子免疫克隆的SAR图像变化检测，包括以下步骤(1)对输入的两时相图像分别采用中值滤波，得到滤波后的图像I1和I2；(2)对滤波后的图像I1和I2求对数比差异影像I3，并将该I3的灰度值作为聚类数据集；(3)设置抗体规模N＝20、类别数k＝2和停机条件，随机产生初始量子抗体作为聚类数据集的初始聚类中心，Q(t)中的都以等概率初始化，所述的停机条件包括最大亲合度值改变量的阈值ε＝10e 5以及连续无法改进次数n＝5；(4)通过观察量子抗体Q(t)的状态，随机产生一个
数，若随机产生的数大于则在相应的二进制位上取1，否则取0，以得到观测后的解每个表示长度为m的二进制串(x1x2…xm)；(5)计算观测后的二进制抗体p(t)与聚类数据集的亲合度函数fk，保留当前种群中的最优子抗体qbest；(6)将量子抗体Q(t)进行量子旋转门变异操作，得到量子种群Qm(t)；(7)将量子种群Qm(t)进行量子全干扰交叉重组操作，得到聚类中心Qc(t)；(8)将聚类中心Qc(t)观测成为新的二进制抗体pc(t)，计算pc(t)中每个子抗体与聚类数据集的亲合度函数值fc；(9)对pc(t)进行选择操作，得到子代抗体p(t+1)；(10)判断子代抗体是否满足停机条件，如果满足该条件就将子代抗体中亲合度最高的抗体对应的图像类属划分作为输出结果，否则返回步骤(4)，循环执行步骤(4)～(10)，直到满足停机条件。FSA00000197006500011.tif,FSA00000197006500012.tif,FSA00000197006500013.tif,FSA00000197006500014.tif,FSA00000197006500015.tif,FSA00000197006500016.tif,FSA00000197006500017.tif
2.根据权利要求1所述的图像变化检测方法，其中步骤(5)所述的计算观测后的二进 制抗体P(t)与聚类数据集的亲合度函数fk，利用如下公式进行 其中，i表示类别，j表示样本点，μ ij表示一个像素点样本j属于各个类别i的隶属 度，m表示模糊指数，d(Xj, Vi)表示第Vi个聚类中心到第Xj个像素样本之间的欧几里德距罔。
3.根据权利要求1所述的图像变化检测方法，其中步骤(6)所述的将量子抗体Q(t)进 行量子旋转门变异操作，是在抗体Q(t)的每一个基本量子位按如下步骤进行(3a)查找量子旋转门变异角θ变化表，如果亲和度函数f (Xi)彡f(beSti)，查询得到 要旋转的角度Δ Qi，其中Xi为当前抗体种群Q(t)的第i个基本量子位，besti为当前的最 优抗体qbest的第i个基本量子位；(3b)根据得到的旋转角度Δ Qi，计算出新的量子位Xi—Mw:Xinew = XiW ( Δ θ,)"cos(A0,) -sin(Ae)1其中“仏幻二 .；‘； ；，为量子门变换矩阵。
4.根据权利要求1所述的图像变化检测方法，其中步骤(7)所述的将量子种群Qm(t) 进行量子全干扰交叉重组操作，是通过使用量子的相干特性构造全干扰交叉操作实现重组 算子，该操作采用对角线交叉的方式，使得种群中的所有抗体均参与交叉。
5.根据权利要求1所述的图像变化检测方法，其中步骤(9)所述的对p。(t)进行选择操 作，采用精英选择策略，即在进化过程中，如果某一代中的最优解的亲合度函数值优于当前 代最优解的亲合度函数值，则当前代最优解就被该最优解所代替，否则保留当前代最优解， 从而可以保证在进化过程中每一代的最优解都不会丢失。
全文摘要
本发明公开一种基于量子免疫克隆的SAR图像变化检测方法，主要解决已有优化方法耗时长，易陷入局部最优，对复杂图像边缘定位不准确的缺点。其步骤为(1)对两时相待变化检测图像进行滤波处理，求对数比差异影像图；(2)设置种群规模、类别数k及停机条件，随机产生量子抗体Q(t)作为初始聚类中心；(3)观测Q(t)成二进制抗体p(t)，计算每个抗体的亲合度fk，保留Q(t)最优抗体qbest；(4)对Q(t)进行变异操作得到Qm(t)；(5)重组Qm(t)得到Qc(t)；(6)观测Qc(t)成二进制抗体pc(t)，计算每个抗体的亲合度fc；(7)选择操作pc(t)，得到子代抗体；(8)若子代抗体满足停机条件，将子代抗体中亲合度最高的抗体对应图像类属划分作为输出结果。本发明具有变化检测精度高、边缘定位准确的优点，可用于对复杂图像的变化检测。
文档编号G01S13/90GK101908213SQ201010230980
公开日2010年12月8日 申请日期2010年7月16日 优先权日2010年7月16日
发明者刘芳, 吴娜娜, 吴建设, 尚荣华, 李阳阳, 焦李成, 缑水平, 马文萍 申请人:西安电子科技大学