T的协方 差；5：ft和Xtf为两个随机变量之间的互协方差；
[0035] 所述特征值计算单元，被配置为计算所述总体方差矩阵的广义特征方 程，并在单位方差约束条件下计算广义特征方程
和
的特征向量a和b，其中，P为对应的广义特征方程的特征根；
[0036] 所述变化参数计算单元，被配置为根据公式M=U-V计算所选择的所述分割图斑 的变化参数M;其中，U=aT*F，V=bVr;
[0037] 所述判定参数计算单元，被配置为通过卡方分布概率密度函数对参数Z进行赋权 计算，
其中，K为每组所述随机变量F或T中的变量个数，为第i个变化 参数M的变量标准差；
[0038] 所述图斑变化判定单元，被配置为判断所述参数Z是否小于等于预设数值，若是， 则判定当前所选择的分割图斑为变化图斑。
[0039] 本发明的有益效果包括：本发明提供的一种遥感影像变化检测的方法及系统，其 通过使用已确定的时相遥感影像及其包含的图斑对要检测的时相遥感影像进行分割，从而 将所处理的时相遥感影像划分出多个规整，且与变化地物的边界能够吻合的图斑。该方法 能够避免传统的多时相高分辨率遥感影像的变化检测中由类内光谱差异引起的细碎图斑 造成的椒盐噪声。其能够更加合理地集中并识别变化信息，提高图斑变化检测的效率及数 据处理的速度，同时也就提高了地物变化检测的效率。明显减低误检和漏检，使检测精度得 到显著提高。
【附图说明】
[0040]图1为本发明一种遥感影像变化检测的方法的一具体实施例的流程图；
[0041] 图2为本发明一种遥感影像变化检测的方法的另一具体实施例的流程图；
[0042]图3为使用基于像元的典型相关分析多元变化检测方法对实验区域的检测结果 图像；
[0043]图4为使用本方法的遥感影像变化检测的方法进行检测结果图像；
[0044] 图5为实验区域地表的真实变化图；
[0045] 图6为实验区域变化检测混淆矩阵对比图；
[0046] 图7 (al)为基于像元的典型相关分析多元变化检测方法的前一时相遥感影像局 部分析示意图；
[0047] 图7化1)为使用本方法的前一时相遥感影像局部示意图；
[0048] 图7 (a2)为与图7 (al)对应的基于像元的典型相关分析多元变化检测方法的后一 时相遥感影像局部分析示意图；
[0049] 图7化2)为与图7化1)对应的使用本方法的后一时相遥感影像局部分析示意图；
[0050] 图8 (al)为基于像元的典型相关分析多元变化检测方法的前一时相遥感影像混 淆矩阵局部分析示意图；
[0051] 图8化1)为使用本方法的前一时相遥感影像混淆矩阵局部示意图；
[0052] 图8 (a2)为与图7 (al)对应的基于像元的典型相关分析多元变化检测方法的后一 时相遥感影像混淆矩阵局部分析示意图；
[0053] 图8化2)为与图7化1)对应的使用本方法的后一时相遥感影像混淆矩阵局部分析 不意图；
[0054] 图9为本发明一种遥感影像变化检测的系统的一具体实施例的结构示意图；
[0055]图10为本发明一种遥感影像变化检测的系统的另一具体实施例的结构示意图；
[0056] 图11为本发明一种遥感影像变化检测的系统的一具体实施例的变化判断模块构 成示意图。
【具体实施方式】
[0057]为了使本发明的目的、技术方案及优点更加清楚明白，W下结合附图对本发明的 遥感影像变化检测的方法及系统的【具体实施方式】进行说明。应当理解，此处所描述的具体 实施例仅仅用W解释本发明，并不用于限定本发明。
[005引本发明一实施例的遥感影像变化检测的方法，如图1所示，包括W下步骤：
[0059]S100，W前一时相遥感影像相对应的矢量图斑为基准对后一时相遥感影像进行分 害d，获取所述后一时相遥感影像的分割图斑。
[0060] 本步骤中，使用前一时相遥感影像的相对应的基期矢量数据为基准，将其与要进 行变化判断的后一时相遥感影像进行套合。可通过设定分割尺度、紧致度W及光谱指数等 参数，在所述前一时相遥感影像相对应的矢量图斑边界限制下，对所述后一时相遥感影像 进行分割，获取同质像斑。采用运种基于已知数据对后期获取的影像进行分割，可使分割结 果兼顾前后两个时相的地物特征，获取到与已知矢量图斑相对应的较为完整、规则的分割 图斑。能够避免椒盐现象。且能够保证分割后获取的图斑"类间光谱相异"和"类内光谱同 质"。
[0061] 本领域技术人员可W理解，前一时相遥感影像和后一时相遥感影像是地理位置相 对应的两帖影像，但是两帖影像在时间上是不同的，前一时相遥感影像是前期获取的，影像 图斑已经确定的遥感影像。而后一时相遥感影像是后期在对地表是否发生变化进行判断 时，所新获取的时相遥感影像。且前一时相遥感影像中的矢量图斑的数量一般为多个，而基 于其进行分割的后一时相遥感影像分割完成后也能分割出多个分割图斑。
[0062] S200,选定一个所述分割图斑，提取所选定的所述分割图斑对应的前一时相遥感 影像区域的全部或部分遥感影像特征构成前一空间向量特征，并根据所述前一空间向量特 征构建前一空间向量。其中，所述遥感影像特征包括光谱特征和纹理特征等。
[0063] S300,提取步骤S200中所选定前述分割图斑对应的后一时相遥感影像区域的与 所述前一空间向量特征相对应的特征后一空间向量特征，并根据所述后一空间特征向量构 建后一至间向量。
[0064] 运两步骤中，W获得的分割结果为基准，将所选择的分割图斑与对应的前一时相 和后一时相遥感影像进行叠置分析，再分别对两个影像图斑的特征进行提取；且可提取两 个影像图斑的全部光谱特征和纹理特征，再通过特征优化选取合适的特征分别构建后一空 间向量和前一空间向量。当然，此处是为了根据提取的特征对后一时相遥感影像的分割图 斑是否发生了变化进行判断，因此，分割图斑中所最后选择的后一空间特征向量特征和相 对应的分割图斑中最后选择的前一空间特征向量特征是相同的，即两者所选择的为相同的 光谱特征和纹理特征中的几个。
[0065] 另外，在提取了选定的分割图斑各自的光谱特征和纹理特征后，可通过特征优化 选取合适的特征分别构建前一空间特征向量和后一空间特征向量。特征优化的原则就是在 最邻近分类器中w距离测度为依据，通过比较已提取的特征，计算不同特征之间的最小距 离，从而得到最佳特征组合，优化之后即可得到最优特征。通过确定合理的波段及最优特征 构建前一空间特征向量和后一空间特征向量，避免数据冗余，提高计算速度和效率。
[0066]S400,通过比较所述后一空间向量和所述前一空间向量判断所选择的分割图斑相 对所述前一时相遥感影像对应区域是否发生变化，若发生了变化，则判定所述后一时相遥 感影像对应区域发生了变化。
[0067] 本发明实施例的遥感影像变化检测的方法，其通过使用已确定的时相遥感影像及 其包含的图斑对要检测的时相遥感影像进行分割，从而将所处理的时相遥感影像划分出多 个规整，且与变化地物的边界能够吻合的图斑。该方法能够避免传统的多时相高分辨率遥 感影像的变化检测中由类内光谱差异引起的细碎图斑造成的椒盐噪声。其能够更加合理地 集中并识别变化信息，提高图斑变化检测的效率及数据处理的速度，同时也就提高了地物 变化检测的效率。明显减低误检和漏检，使检测精度得到显著提高。
[0068] 其中，步骤S400,通过比较所述后一空间向量和所述前一空间向量判断所选择的 分割图斑是否发生变化中，当所选择的分割图斑发生的变化超过预设阔值时，判定所述后 一时相遥感影像中所选择的所述分割图斑发生了变化。
[0069] 本实施例中，设定变化的预设值，并将数据处理完成后的变化大小与预设的预设 变化阔值进行比较，在发生的变化超过设定的预设变化阔值时，判定所处理的分割图斑发 生了变化。如此设置，可忽略发生微小变化的图斑，根据需要的地物变化情况设定相应的阔 值则可检测到相应发生变化的图斑。如在只需要确定地物是否发生巨大改变时，则可设置 较大的预设变化阔值，当需要将地物较小的改变也检测到时，则需要设定较小的所述预设 变化阔值，从而可W检测到发生微小变化的分割图斑。
[0070] 由前面的描述可知，步骤S100对后一时相遥感影像进行分割后会得到多个分割 图斑，而后续的步骤S200~S400是对其中一个分割图斑是否为变化图斑进行的判断。本 领域技术人员可W理解，在根据时相遥感影像对地物变化进行检测时，一般需要判断所处 理