Image change pixel level detection method and system
附图说明 图1为本发明的整体模型框架图； 图2为本发明中残差空谱协同注意力模块的结构示意图； 图3为本发明中图像语义分割网络的结构图。 技术领域 本发明涉及图像变化检测技术领域，具体涉及一种图像变化像素级检测方法与系统。 具体实施方式 下面结合附图对本发明的一种优选实施方式作详细的说明。 本实施例公开了一种图像变化像素级检测方法，具体包括如下步骤： S1、选取图像并进行预处理，获得具有成对图像的训练数据集和测试数据集，具体包括： S11、准备训练数据集和测试数据集，数据集中的图像包括但不限于遥感图像、地震图像； S12、先将历史图像和当前图像位置对齐，对于每张图像，进行像素归一化处理，得到图像X；其中历史图像即当前图像对应的变化前图像，也可称为参考图像。 S2、搭建图像变化检测模型，包括特征提取子网络、融合模块、残差空谱协同注意力模块ResS2At，特征提取子网络通过图像语义分割网络得到，图像语义分割网络包括下采样路径和上采样路径，所述特征提取子网络为图像语义分割网络的下采样路径；特征提取子网络与融合模块组成多尺度特征差分融合子网络。具体如下： S21、定义残差空谱协同注意力模块： 残差空谱协同注意力模块ResS2At（Residual Spatio-Spectral Attention）的结构如图2所示。对于给定的输入特征图组Ma0，依次经过两个N1模块，得到特征图组Ma1，其中N1模块由卷积层、批归一化和ReLU激活函数组成，卷积核尺寸为3×3，步长为1，填充为1。 特征组图Ma1输入到谱注意力模块CAM，CAM输出结果与特征图组Ma1按通道相乘后得到特征图组Ma2；其中CAM针对特征图组Ma1的各个通道计算注意力强度，对于特征图组，表示实数域，w、h、c分别表示特征图宽、高以及特征图组通道数，将特征图组Ma1中的每张特征图的最大值组成通道描述向量，将Ma1中每张特征图的平均值组成通道描述向量，和分别通过全连接层，得到向量，和相加后并输入到Sigmoid激活函数，得到的结果即为CAM的输出结果。CAM输出结果为的向量，向量的每个元素对应Ma1中每张特征图的通道注意力强度，上述的按通道相乘指的是每个通道的注意力强度与对应特征图所有像素值相乘。 特征图组Ma2输入到空间注意力模块SAM，SAM输出结果与Ma2按空间相乘后得到特征图组Ma3；其中SAM针对特征图组Ma2的各个空间坐标计算注意力强度；对特征图组Ma2的同一空间坐标上的像素点求取最大值，得到空间描述矩阵，对特征图组Ma2的同一空间坐标上的像素点求取平均值，得到空间描述矩阵，将和进行特征拼接后依次输入卷积层和Sigmoid激活函数，得到的结果即为SAM的输出结果。SAM输出结果为的矩阵，矩阵的每个元素对应Ma2中的每个空间坐标的空间注意力强度，上述按空间相乘指的是每个空间坐标的注意力强度与对应空间坐标所有像素值相乘。 最后，Ma3与Ma1相加得到特征图组Ma4，Ma4与Ma0相加后并经过ReLU激活函数，得到残差空谱协同注意力模块的最终输出。 S22、定义特征提取子网络： 特征提取子网络由预训练的图像语义分割网络得到，是图像语义分割网络的一部分，图像语义分割网络的结构如图3所示。图像语义分割网络是由基于结合小波变换的上采样路径（解码部分）、下采样路径（编码部分）和残差空谱协同注意力模块构建的U型深度网络分割模型。 图像语义分割网络的下采样路径整体如下：原始图像X依次经过N2模块、N1模块输出特征图组Md1，Md1依次经过N5模块、N1模块、N1模块，输出特征图组Md2，Md2依次经过N5模块、N1模块、N1模块输出特征图组Md3，Md3依次经过N5模块、N1模块、N1模块输出特征图组Md4，Md4依次经过N5模块、N1模块输出特征图组Md5。 图像语义分割网络的上采样路径整体如下：Md5经过N6模块后与Md4按通道合并后，依次经过N1模块、N1模块输出特征图组Mu4，Mu4经过N6模块后与Md3按通道合并，再依次经过N1模块、N1模块输出特征图组Mu3，Mu3经过N6模块后与Md2按通道合并，再依次经过N1模块、N1模块输出特征图组Mu2，Mu2经过N6模块后与Md1按通道合并，再依次经过N1模块、N1模块、N4模块、N3模块，输出分割结果图R。 特征提取子网络即图像语义分割网络的下采样路径。 N1模块由卷积层、批归一化和ReLU激活函数组成；N2模块由卷积层、批归一化组成，N3模块由卷积层、SoftMax激活函数组成，N4模块即为残差空谱协同注意力模块ResS2At；N5模块为离散小波变换DWT，用于下采样；N6模块为小波逆变换IWT，用于上采样；上述卷积核可设置为尺寸为3×3，步长为1，填充为1，DWT和IWT可设置为使用Haar小波的2DDWT。具体的，使用Haar小波的2D DWT对于给定的特征图Md进行变换，使用四个步长为2的卷积滤波器，将特征图Md分解成四个子带图像，即、、和，其中表示平均信息，表示水平信息，表示垂直信息，表示对角线信息，为低通滤波器，、、为高通滤波器，表示卷积运算符，表示因子为2的标准下采样运算符；使用Haar小波的2D IWT是2D DWT的逆变化，可将输入的四个子带图重构为原始图像。 本实施例中，下采样路径由四个下采样模块组成，四个下采样模块从浅至深形式依次如下： Conv(1,64,3)+BN+Conv(64,64,3)+BN+ReLU+DWT； Conv(256,128,3)+BN+ReLU+Conv(128,128,3)+BN+ReLU+DWT； Conv(512,256,3)+BN+ReLU+Conv(256,256,3)+BN+ReLU+DWT； Conv(1024,512,3)+BN+ReLU+Conv(512,512,3)+BN+ReLU+DWT。 其中，Conv()表示卷积操作，其内三个参数分别为输入通道数、输出通道数、方形卷积核边长，Conv()的卷积步长为1，填充为1，以保证卷积后特征图大小不变；BN表示批归一化；ReLU为ReLU激活函数；特征图组经过DWT后，特征图大小缩小为原来的一半，通道数增加到原来的4倍。 上采样路径由一个卷积层N1模块和四个上采样模块组成，从浅至深形式依次如下： Conv(2048,1024,3)+BN+ReLU； IWT+Conv(768,512,3)+BN+ReLU+Conv(512,512,3)+BN+ReLU； IWT+Conv(384,256,3)+BN+ReLU+Conv(256,256,3)+BN+ReLU； IWT+Conv(192,128,3)+BN+ReLU+Conv(128,128,3)+BN+ReLU； IWT+Conv(96,64,3)+BN+ReLU+Conv(64,64,3)+BN+ReLU+ResS2At+Conv(64,1,3)+SoftMax。 其中，Conv的卷积步长为1，填充为1，特征图组经过IWT后，特征图大小增大为原来的2倍，通道数缩减到原来的1/4倍；SoftMax为SoftMax激活函数。 S23、定义多尺度特征差分融合子网络： 多尺度特征差分融合子网络由特征提取子网络和融合模块组成，由于特征提取子网络从预训练的图像语义分割网络中提取，因此特征提取子网络在多尺度特征差分融合子网络中是不可训练的。 如图1所示，将成对的历史图像Xr和当前图像Xc，输入到固定的特征提取子网络，分别得到与历史图像Xr对应的特征图组Mdr1、Mdr2、Mdr3、Mdr4、Mdr5，与当前图像对应的特征图组Mdc1、Mdc2、Mdc3、Mdc4、Mdc5，这里的Md对应到S22下采样路径中特征图组。(Mdr1,Mdc1)，(Mdr2,Mdc2)，(Mdr3,Mdc3)，(Mdr4,Mdc4)，(Mdr5,Mdc5)成组输入到N7模块，得到不同尺度的特征差分图组D1至D5，D1至D5经过融合模块按通道进行合并，得到多通道的特征差分图组，其中，N7模块为差分后上采样模块，Conc是Concatenate的缩写，表示特征连接合并；多通道的特征差分图组依次经过N4模块、N4模块和N8模块，即可得到历史图像Xr和当前图像Xc的变化图。该变化图为灰度图，发生变化像素点接近1，即亮度较大，无变化像素点接近0，即亮度较低。本实施例中，多通道的特征差分图组需要重复经过两次残差空谱协同注意力模块ResS2At，N8模块由1x1卷积层和Sigmoid激活函数组成。 S3、定义损失函数并进行网络训练： 首先训练图像语义分割网络，以遥感图像树木变化检测为例，训练图像语义分割网络时采用的标签为：树木、非树木，图像语义分割网络训练完毕后提取其下采样路径得到特征提取子网络，图像语义分割网络训练过程中采用交叉熵损失函数来表达像素分类误差。 将预训练的特征提取子网络并入多尺度特征差分融合子网络并固定特征提取子网络所有参数不参与训练，仅训练多尺度特征差分融合子网络中除去特征提取子网络的其余部分，标签为人工给定的变化图，其中变化图、历史图像、当前图像尺寸一致，变化图为灰度图，像素亮度最大表示发生变化，像素亮度最小表示未发生变化；多尺度特征差分融合子网络训练中采用交叉熵损失函数来表达像素级变化检测误差。 S4、模型应用： 把测试数据集中成对的历史图像和当前图像输入到训练好的图像变化检测模型中，得到像素级的变化图。 本发明还公开了一种图像变化像素级检测系统，本发明中的系统与方法对应，方法中的具体方案同样适用于系统。 对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内，不应将权利要求中的任何附图标记视为限制所涉及的权利要求。 此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立技术方案，说明书的这种叙述方式仅仅是为了清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。 背景技术 图像变化检测作为计算机视觉中重要的一个方向，利用计算机技术对不同时刻拍摄到的同一地理区域的图像信息进行识别，并从中提取最为关键的视觉特征信息作为具体特征，最终确定变化信息，其研究的目的是找到感兴趣的变化信息，滤除干扰因素出现的不相干变化信息。该技术在资源和环境检测、自然灾害的监测评估、测绘以及军事方面极为重要，如对森林或者植被的动态变化检测，对农作物的生长情况、城镇变化(街道、建筑物等)的实时监测，对军事战略目标进行动态监视等方面，也可用于时延地震变化检测，以发现油藏变化情况。 近年来，随着人工智能技术的不断发展，基于深度学习的图像变化检测技术也不断涌现，大量研究表明，基于深度学习方法构建的模型在图像差异检测问题中能够取得良好的性能，但由于不确定因素的影响，如光照、季节、气候、拍摄角度等因素产生的其他变化，这种伪变化差异会对最终结果产生负面影响，正常去除伪变化需人工剔除，耗时耗力，效率低，局限性较大。同时，许多网络使用深层的特征信息，对细节纹理（如小物体）的检测存在不确定性，即小物体只有较少的像素信息，在下采样的过程当中容易丢失本身的主要特征信息。此外，在小样本情况下的图像差异检测很难达到满意的效果，会存在偏差以及缺乏泛化能力的问题。 发明内容 为解决上述技术问题，本发明提供一种图像变化像素级检测方法与系统。 为解决上述技术问题，本发明采用如下技术方案： 一种图像变化像素级检测方法，将采集到的历史图像和当前图像输入到完成训练的图像变化检测模型中，得到像素级的变化图； 图像变化检测模型的构建和训练方法，包括以下步骤： 步骤一：对多个成对的历史图像Xr和当前图像Xc进行预处理后，得到训练数据集； 步骤二、构建基于图像语义分割网络的图像变化检测模型： 图像变化检测模型包括特征提取子网络、融合模块、残差空谱协同注意力模块； 特征提取子网络通过图像语义分割网络得到，图像语义分割网络包括下采样路径和上采样路径，所述特征提取子网络为图像语义分割网络的下采样路径；特征提取子网络与融合模块组成多尺度特征差分融合子网络； 将训练数据集中成对的历史图像Xr和当前图像Xc输入到图像变化检测模型中时，历史图像Xr和当前图像Xc先输入到特征提取子网络，分别得到特征图组和特征图组,1≤i≤n，n为特征提取子网络输出的特征图组总数；将特征图组对(Mdri,Mdci)进行差分后上采样，将得到的不同尺度的特征差分图组Di输入至融合模块并按通道进行合并，得到多通道的特征差分图组，多通道的特征差分图组依次经过残差空谱协同注意力模块、卷积层和Sigmoid激活函数，得到历史图像Xr与当前图像Xc之间的变化图； 步骤三、定义损失函数并进行网络训练： 对图像语义分割网络进行训练，训练完成后提取图像语义分割网络的下采样路径得到预训练的特征提取子网络，图像语义分割网络训练过程中采用交叉熵损失函数表达像素分类误差； 将预训练的特征提取子网络并入多尺度特征差分融合子网络，在多尺度特征差分融合子网络训练过程中，特征提取子网络所有参数保持固定，标签为人工给定的变化图；其中变化图、历史图像和当前图像尺寸一致；变化图为灰度图，像素亮度最大表示发生变化，像素亮度最小表示未发生变化；多尺度特征差分融合子网络训练时，采用交叉熵损失函数表达像素级变化检测误差。 进一步地，图像包括遥感图像和地震图像。 进一步地，步骤一中，对多个成对的历史图像Xr和当前图像Xc进行预处理时，将历史图像和当前图像的位置对齐，并进行像素归一化处理。 进一步地，所述图像语义分割网络为结合小波变换的上采样路径、下采样路径和残差空谱协同注意力模块组成的U型网络。 进一步地，所述残差空谱协同注意力模块依次包括两个N1模块、谱注意力模块CAM、空间注意力模块SAM； 步骤二中，输入到残差空谱协同注意力模块中的多通道的特征差分图组记为特征图组Ma0； 特征图组Ma0依次经过两个N1模块，得到特征图组Ma1，其中N1模块包括卷积层、批归一化和ReLU激活函数； 特征组图Ma1输入到谱注意力模块CAM，CAM输出结果与特征图组Ma1按通道相乘后得到特征图组Ma2；其中CAM针对特征图组Ma1的各个通道计算注意力强度，对于特征图组，表示实数域，w、h、c分别表示特征图组Ma1中特征图的宽、高以及特征图组Ma1通道数，将特征图组Ma1中的每张特征图的最大值组成通道描述向量，将Ma1中每张特征图的平均值组成通道描述向量，和分别通过全连接层，得到向量，和相加后并输入到Sigmoid激活函数，得到的结果即为CAM的输出结果； 特征图组Ma2输入到空间注意力模块SAM，SAM输出结果与Ma2按空间相乘后得到特征图组Ma3；其中SAM针对特征图组Ma2的各个空间坐标计算注意力强度；对特征图组Ma2的同一空间坐标上的像素点求取最大值，得到空间描述矩阵，对特征图组Ma2的同一空间坐标上的像素点求取平均值，得到空间描述矩阵，将和进行特征拼接后依次输入卷积层和Sigmoid激活函数，得到的结果即为SAM的输出结果； 最后，Ma3与Ma1相加得到特征图组Ma4，Ma4与Ma0相加后并经过ReLU激活函数，得到残差空谱协同注意力模块的最终输出。 进一步地，历史图像Xr和当前图像Xc统一记为图像X； 图像语义分割网络的下采样路径如下：图像X依次经过N2模块、N1模块输出特征图组，依次经过N5模块、N1模块、N1模块，输出特征图组Md2，Md2依次经过N5模块、N1模块、N1模块输出特征图组Md3，Md3依次经过N5模块、N1模块、N1模块输出特征图组Md4，Md4依次经过N5模块、N1模块输出特征图组Md5； 图像语义分割网络的上采样路径如下：Md5经过N6模块后与Md4按通道合并后，再依次经过N1模块、N1模块输出特征图组Mu4，Mu4经过N6模块后与Md3按通道合并，再依次经过N1模块、N1模块输出特征图组Mu3，Mu3经过N6模块后与Md2按通道合并，再依次经过N1模块、N1模块输出特征图组Mu2，Mu2经过N6模块后与Md1按通道合并，再依次经过N1模块、N1模块、N4模块、N3模块； N1模块由卷积层、批归一化和ReLU激活函数组成；N2模块由卷积层、批归一化组成；N3模块由卷积层、SoftMax激活函数组成；N4模块即为残差空谱协同注意力模块ResS2At；N5模块为离散小波变换，用于下采样；N6模块为小波逆变换，用于上采样。 进一步地，下采样路径包括四个下采样模块，四个下采样模块从浅至深形式依次如下： Conv(1,64,3)+BN+Conv(64,64,3)+BN+ReLU+DWT； Conv(256,128,3)+BN+ReLU+Conv(128,128,3)+BN+ReLU+DWT； Conv(512,256,3)+BN+ReLU+Conv(256,256,3)+BN+ReLU+DWT； Conv(1024,512,3)+BN+ReLU+Conv(512,512,3)+BN+ReLU+DWT； 其中，Conv()表示卷积操作，Conv()内的三个参数分别为输入通道数、输出通道数、方形卷积核边长；BN表示批归一化；ReLU为ReLU激活函数；DWT为离散小波变换；特征图组经过DWT后，特征图组内的各特征图大小缩小为原来的一半，通道数增加到原来的四倍。 进一步地，上采样路径由一个N1模块和四个上采样模块组成，从浅至深形式依次如下： Conv(2048,1024,3)+BN+ReLU； IWT+Conv(768,512,3)+BN+ReLU+Conv(512,512,3)+BN+ReLU； IWT+Conv(384,256,3)+BN+ReLU+Conv(256,256,3)+BN+ReLU； IWT+Conv(192,128,3)+BN+ReLU+Conv(128,128,3)+BN+ReLU； IWT+Conv(96,64,3)+BN+ReLU+Conv(64,64,3)+BN+ReLU+ResS2At+Conv(64,1,3)+SoftMax； 其中，Conv()表示卷积操作，Conv()内的三个参数分别为输入通道数、输出通道数、方形卷积核边长；BN表示批归一化；ReLU为ReLU激活函数；IWT为小波逆变换；特征图组经过IWT后，特征图组内的各特征图大小增大为原来的2倍，通道数缩减到原来的1/4倍；SoftMax为SoftMax激活函数。 一种图像变化像素级检测系统，包括： 数据获取模块：对多个成对的历史图像Xr和当前图像Xc进行预处理后，得到训练数据集； 网络构建模块：构建基于图像语义分割网络的图像变化检测模型；图像变化检测模型包括特征提取子网络、融合模块、残差空谱协同注意力模块；特征提取子网络通过图像语义分割网络得到，图像语义分割网络包括下采样路径和上采样路径，所述特征提取子网络为图像语义分割网络的下采样路径；特征提取子网络与融合模块组成多尺度特征差分融合子网络；将训练数据集中成对的历史图像Xr和当前图像Xc输入到图像变化检测模型中时，历史图像Xr和当前图像Xc输入到特征提取子网络，分别得到特征图组Mdri和特征图组Mdci,1≤i≤n，n为特征提取子网络输出的特征图组总数；将特征图组对(Mdri,Mdci)进行差分后上采样，将得到的不同尺度的特征差分图组Di输入至融合模块并按通道进行合并，得到多通道的特征差分图组，多通道的特征差分图组依次经过残差空谱协同注意力模块、卷积层和Sigmoid激活函数，得到历史图像Xr与当前图像Xc之间的变化图； 网络训练模块：对图像语义分割网络进行训练，训练完成后提取图像语义分割网络的下采样路径得到预训练的特征提取子网络，图像语义分割网络训练过程中采用交叉熵损失函数来表达像素分类误差；将预训练的特征提取子网络并入多尺度特征差分融合子网络，在多尺度特征差分融合子网络训练过程中，特征提取子网络所有参数保持固定，标签为人工给定的变化图；其中变化图、历史图像、当前图像尺寸一致；变化图为灰度图，像素亮度最大表示发生变化，像素亮度最小表示未发生变化；多尺度特征差分融合子网络训练时，采用交叉熵损失函数来表达像素级变化检测误差； 推理模块：将实际采集到的历史图像和当前图像输入到完成训练的图像变化检测模型中，得到像素级的变化图。 本发明中的系统与方法对应，方法中的优选方案同样适用于系统。 与现有技术相比，本发明的有益技术效果是： 1、训练图像语义分割网络以获取特征提取子网络，得到具有目标导向的特征图，能够显著消除非目标的影响；2、小波变换及其逆变换的引入使网络的上下采样信息损失减小；3、多种注意力机制的引入能够有效提升特征图融合性能；4、实现目标变化的像素级检测。 The invention relates to the technical field of image change detection, and discloses an image change pixel-level detection method and system, and the method comprises the following steps: carrying out the preprocessing of a plurality of paired historical images Xr and a current image Xc, and forming a training data set; constructing an image change detection model based on an image semantic segmentation network; defining a loss function and carrying out network training; inputting the collected historical image and the current image into the trained image change detection model to obtain a pixel-level change graph, and realizing pixel-level detection of target change; according to the method, the feature extraction sub-network is obtained through the pre-trained image semantic segmentation network, the target-oriented feature map is obtained, and the influence of non-targets can be remarkably eliminated; the loss of upper and lower sampling information of the network is reduced by introducing wavelet transform and inverse transform; the introduction of various attention mechanisms can effectively improve the feature map fusion performance. 1.一种图像变化像素级检测方法，其特征在于，将采集到的历史图像和当前图像输入到完成训练的图像变化检测模型中，得到像素级的变化图； 图像变化检测模型的构建和训练方法，包括以下步骤： 步骤一：对多个成对的历史图像Xr和当前图像Xc进行预处理后，得到训练数据集； 步骤二、构建基于图像语义分割网络的图像变化检测模型： 图像变化检测模型包括特征提取子网络、融合模块、残差空谱协同注意力模块； 特征提取子网络通过图像语义分割网络得到，图像语义分割网络包括下采样路径和上采样路径，所述特征提取子网络为图像语义分割网络的下采样路径；特征提取子网络与融合模块组成多尺度特征差分融合子网络； 将训练数据集中成对的历史图像Xr和当前图像Xc输入到图像变化检测模型中时，历史图像Xr和当前图像Xc先输入到特征提取子网络，分别得到特征图组和特征图组,1≤i≤n，n为特征提取子网络输出的特征图组总数；将特征图组对(Mdri,Mdci)进行差分后上采样，将得到的不同尺度的特征差分图组Di输入至融合模块并按通道进行合并，得到多通道的特征差分图组，多通道的特征差分图组依次经过残差空谱协同注意力模块、卷积层和Sigmoid激活函数，得到历史图像Xr与当前图像Xc之间的变化图； 步骤三、定义损失函数并进行网络训练： 对图像语义分割网络进行训练，训练完成后提取图像语义分割网络的下采样路径得到预训练的特征提取子网络，图像语义分割网络训练过程中采用交叉熵损失函数表达像素分类误差； 将预训练的特征提取子网络并入多尺度特征差分融合子网络，在多尺度特征差分融合子网络训练过程中，特征提取子网络所有参数保持固定，标签为人工给定的变化图；其中变化图、历史图像和当前图像尺寸一致；变化图为灰度图，像素亮度最大表示发生变化，像素亮度最小表示未发生变化；多尺度特征差分融合子网络训练时，采用交叉熵损失函数表达像素级变化检测误差。 2.根据权利要求1所述的图像变化像素级检测方法，其特征在于：图像包括遥感图像和地震图像。 3.根据权利要求1所述的图像变化像素级检测方法，其特征在于：步骤一中，对多个成对的历史图像Xr和当前图像Xc进行预处理时，将历史图像和当前图像的位置对齐，并进行像素归一化处理。 4.根据权利要求1所述的图像变化像素级检测方法，其特征在于：所述图像语义分割网络为结合小波变换的上采样路径、下采样路径和残差空谱协同注意力模块组成的U型网络。 5.根据权利要求1所述的图像变化像素级检测方法，其特征在于：所述残差空谱协同注意力模块依次包括两个N1模块、谱注意力模块CAM、空间注意力模块SAM； 步骤二中，输入到残差空谱协同注意力模块中的多通道的特征差分图组记为特征图组Ma0； 特征图组Ma0依次经过两个N1模块，得到特征图组Ma1，其中N1模块包括卷积层、批归一化和ReLU激活函数； 特征组图Ma1输入到谱注意力模块CAM，CAM输出结果与特征图组Ma1按通道相乘后得到特征图组Ma2；其中CAM针对特征图组Ma1的各个通道计算注意力强度，对于特征图组，表示实数域，w、h、c分别表示特征图组Ma1中特征图的宽、高以及特征图组Ma1通道数，将特征图组Ma1中的每张特征图的最大值组成通道描述向量，将Ma1中每张特征图的平均值组成通道描述向量，和分别通过全连接层，得到向量，和相加后并输入到Sigmoid激活函数，得到的结果即为CAM的输出结果； 特征图组Ma2输入到空间注意力模块SAM，SAM输出结果与Ma2按空间相乘后得到特征图组Ma3；其中SAM针对特征图组Ma2的各个空间坐标计算注意力强度；对特征图组Ma2的同一空间坐标上的像素点求取最大值，得到空间描述矩阵，对特征图组Ma2的同一空间坐标上的像素点求取平均值，得到空间描述矩阵，将和进行特征拼接后依次输入卷积层和Sigmoid激活函数，得到的结果即为SAM的输出结果； 最后，Ma3与Ma1相加得到特征图组Ma4，Ma4与Ma0相加后并经过ReLU激活函数，得到残差空谱协同注意力模块的最终输出。 6.根据权利要求1所述的图像变化像素级检测方法，其特征在于：历史图像Xr和当前图像Xc统一记为图像X； 图像语义分割网络的下采样路径如下：图像X依次经过N2模块、N1模块输出特征图组，依次经过N5模块、N1模块、N1模块，输出特征图组Md2，Md2依次经过N5模块、N1模块、N1模块输出特征图组Md3，Md3依次经过N5模块、N1模块、N1模块输出特征图组Md4，Md4依次经过N5模块、N1模块输出特征图组Md5； 图像语义分割网络的上采样路径如下：Md5经过N6模块后与Md4按通道合并后，再依次经过N1模块、N1模块输出特征图组Mu4，Mu4经过N6模块后与Md3按通道合并，再依次经过N1模块、N1模块输出特征图组Mu3，Mu3经过N6模块后与Md2按通道合并，再依次经过N1模块、N1模块输出特征图组Mu2，Mu2经过N6模块后与Md1按通道合并，再依次经过N1模块、N1模块、N4模块、N3模块； N1模块由卷积层、批归一化和ReLU激活函数组成；N2模块由卷积层、批归一化组成；N3模块由卷积层、SoftMax激活函数组成；N4模块即为残差空谱协同注意力模块ResS2At；N5模块为离散小波变换，用于下采样；N6模块为小波逆变换，用于上采样。 7.根据权利要求6所述的图像变化像素级检测方法，其特征在于，下采样路径包括四个下采样模块，四个下采样模块从浅至深形式依次如下： Conv(1,64,3)+BN+Conv(64,64,3)+BN+ReLU+DWT； Conv(256,128,3)+BN+ReLU+Conv(128,128,3)+BN+ReLU+DWT； Conv(512,256,3)+BN+ReLU+Conv(256,256,3)+BN+ReLU+DWT； Conv(1024,512,3)+BN+ReLU+Conv(512,512,3)+BN+ReLU+DWT； 其中，Conv()表示卷积操作，Conv()内的三个参数分别为输入通道数、输出通道数、方形卷积核边长；BN表示批归一化；ReLU为ReLU激活函数；DWT为离散小波变换；特征图组经过DWT后，特征图组内的各特征图大小缩小为原来的一半，通道数增加到原来的四倍。 8.根据权利要求6所述的图像变化像素级检测方法，其特征在于：上采样路径由一个N1模块和四个上采样模块组成，从浅至深形式依次如下： Conv(2048,1024,3)+BN+ReLU； IWT+Conv(768,512,3)+BN+ReLU+Conv(512,512,3)+BN+ReLU； IWT+Conv(384,256,3)+BN+ReLU+Conv(256,256,3)+BN+ReLU； IWT+Conv(192,128,3)+BN+ReLU+Conv(128,128,3)+BN+ReLU； IWT+Conv(96,64,3)+BN+ReLU+Conv(64,64,3)+BN+ReLU+ResS2At+Conv(64,1,3)+SoftMax； 其中，Conv()表示卷积操作，Conv()内的三个参数分别为输入通道数、输出通道数、方形卷积核边长；BN表示批归一化；ReLU为ReLU激活函数；IWT为小波逆变换；特征图组经过IWT后，特征图组内的各特征图大小增大为原来的2倍，通道数缩减到原来的1/4倍；SoftMax为SoftMax激活函数。 9.一种图像变化像素级检测系统，其特征在于，包括： 数据获取模块：对多个成对的历史图像Xr和当前图像Xc进行预处理后，得到训练数据集； 网络构建模块：构建基于图像语义分割网络的图像变化检测模型；图像变化检测模型包括特征提取子网络、融合模块、残差空谱协同注意力模块；特征提取子网络通过图像语义分割网络得到，图像语义分割网络包括下采样路径和上采样路径，所述特征提取子网络为图像语义分割网络的下采样路径；特征提取子网络与融合模块组成多尺度特征差分融合子网络；将训练数据集中成对的历史图像Xr和当前图像Xc输入到图像变化检测模型中时，历史图像Xr和当前图像Xc输入到特征提取子网络，分别得到特征图组Mdri和特征图组Mdci,1≤i≤n，n为特征提取子网络输出的特征图组总数；将特征图组对(Mdri,Mdci)进行差分后上采样，将得到的不同尺度的特征差分图组Di输入至融合模块并按通道进行合并，得到多通道的特征差分图组，多通道的特征差分图组依次经过残差空谱协同注意力模块、卷积层和Sigmoid激活函数，得到历史图像Xr与当前图像Xc之间的变化图； 网络训练模块：对图像语义分割网络进行训练，训练完成后提取图像语义分割网络的下采样路径得到预训练的特征提取子网络，图像语义分割网络训练过程中采用交叉熵损失函数来表达像素分类误差；将预训练的特征提取子网络并入多尺度特征差分融合子网络，在多尺度特征差分融合子网络训练过程中，特征提取子网络所有参数保持固定，标签为人工给定的变化图；其中变化图、历史图像、当前图像尺寸一致；变化图为灰度图，像素亮度最大表示发生变化，像素亮度最小表示未发生变化；多尺度特征差分融合子网络训练时，采用交叉熵损失函数来表达像素级变化检测误差； 推理模块：将实际采集到的历史图像和当前图像输入到完成训练的图像变化检测模型中，得到像素级的变化图。