一种全极化SAR影像变化区域精确提取方法技术

本发明专利技术公开了一种全极化SAR影像变化区域精确提取方法，涉及遥感变化检测领域，基本步骤为：对变化前后的全极化SAR影像分别进行特征提取；利用特征影像计算得到差异影像；综合各差异影像将所有像素点逐一判定为变化像素点、未变化像素点或待确定像素点；利用变化像素点和未变化像素点，判定待确定像素点为变化像素点或未变化像素点；利用变化像素点提取全极化SAR影像上的变化区域。本发明专利技术综合考虑各特征的判定结果，逐像素进行全极化SAR影像变化检测，避免了用单一特征进行变化检测所带来的误检漏检，具有科学合理、易于实现、精度高等优点。

【技术实现步骤摘要】
一种全极化SAR影像变化区域精确提取方法
本专利技术涉及遥感变化检测领域，具体说是一种全极化SAR影像变化区域精确提取方法。
技术介绍
比值法和分类后比较法是利用全极化SAR影像进行变化检测的两种常用方法。比值法通过设置合适阈值进行区域划分，具有阈值选取困难，特征利用率低等问题；分类后比较法将变化前后的全极化SAR影像进行地物分类，通过比较变化前后地物分类图得到变化类型，其实验结果受地物分类误差累积影响，变化类型的识别精度受限。全极化SAR影像具有丰富的极化信息，可以计算得到多个特征，不同特征对同一地物表现不同，对地物是否发生变化的判定结果可能不一致，综合考虑各特征的判定可以避免特征利用率低和地物分类误差累积等问题。本专利技术提出了一种全极化SAR影像变化区域精确提取方法。
技术实现思路
(一)要解决的技术问题本专利技术提供一种全极化SAR影像变化区域精确提取方法，该方法综合考虑全极化SAR影像各特征在同一像素上的判定结果进行变化区域精确提取。(二)技术方案本专利技术包含以下步骤：(1)提取特征，对大小为M×N全极化SAR影像IA提取n个特征，得到特征影像Fi(i＝1,2,3....n)，对大小为M×N全极化SAR影像IB同样提取n个特征，得到特征影像Ei(i＝1,2,3....n)；其中IA为变化前的全极化SAR影像，IB为变化后的全极化SAR影像，i表示所提取的特征序号；n为13，13个特征分别为：极化度，极化总功率，一致性参数，熵，二阶矩，对比度，相关性，极化熵，平均散射角，反熵，表面散射，二次散射和体散射；(2)生成差异影像，根据步骤(1)中得到的特征影像，计算得到差异影像ri(i＝1,2,3,···,n)各像素点的像素值，即ri(s,t)＝Fi(s,t)/Ei(s,t)(s＝1,2,···,M；t＝1,2,···,N)其中ri表示第i个特征的差异影像，(s,t)表示差异影像上任意像素点，ri(s,t)为差异影像ri中像素点(s,t)的像素值，Fi(s,t)表示特征影像Fi中像素点(s,t)的像素值，Ei(s,t)表示特征影像Ei中像素点(s,t)的像素值；(3)根据步骤(2)中得到的差异影像ri，判定所有像素点为变化像素点、未变化像素点或待确定像素点；1)所有像素点均为判定对象，则任意像素点均可判定为三种可能事件：其中H1表示“像素点变化”，H2表示“像素点未变化”，H3表示“像素点待确定”；2)将每个特征的差异影像ri分别作为判别依据，求得所有像素点在每个判别依据下三种可能事件的信任程度，其计算公式如下：式中mi(H1)表示在第i个判别依据下，对“像素点变化”事件的信任程度；mi(H2)表示在第i个判别依据下，对“像素点未变化”事件的信任程度；mi(H3)表示在第i个判别依据下，对“像素点待确定”事件的信任程度；表示取其中的最大值；表示取其中的最大值，Ci为第i个差异影像中最大像素值，ci为第i个差异影像中最小像素值，α表示不确定度，其计算公式如下：3)将所有像素点分别在n个判别依据下进行判定之后，对n个判别依据进行合成，计算所有像素点在n个判别依据下三种可能事件的信任程度，其计算公式如下：式中m(H1)表示在n个判别依据下，对“像素点变化”事件的信任程度；m(H2)表示在n个判别依据下，对“像素点未变化”事件的信任程度；m(H3)表示在n个判别依据下，对“像素点待确定”事件的信任程度，K为冲突因子，计算公式如下：其中Ai(i＝1,2,···,n)表示在第i个判别依据下的可能事件，取值为H1,H2,H3；4)当m(H1)大于0.85时，判定像素点(s,t)发生变化，当m(H2)大于0.85时，判定像素点(s,t)未发生变化，当m(H1)小于0.85且m(H2)小于0.85时，判定像素点(s,t)为待确定像素点。(4)判定待确定像素点是否变化，根据步骤(2)中得到的差异影像ri和步骤(3)中所有像素点的判定结果，使用SVM模型对所有待确定像素点进行判定；(5)全极化SAR影像变化区域划分，结合步骤(3)和步骤(4)的判定结果，提取整幅全极化SAR影像的变化区域。(三)有益效果本专利技术的优点体现在：由于利用单一特征的差异影像进行变化检测，误检漏检现象严重，具有不确定性；且不同特征的差异影像判定变化区域具有冲突性。本专利技术综合考虑各特征的判定结果，逐像素进行全极化SAR影像变化检测，以提高检测精度。附图说明图1为本专利技术实施的步骤流程图。具体实施方式为使本专利技术的目的、内容、和优点更加清楚，下面结合附图和实施例，对本专利技术的具体实施方式作进一步详细描述：参照图1，本专利技术的具体实施步骤为：(1)提取特征，对大小为M×N全极化SAR影像IA提取n个特征，得到特征影像Fi(i＝1,2,3....n)，对大小为M×N全极化SAR影像IB同样提取n个特征，得到特征影像Ei(i＝1,2,3....n)；其中IA为变化前的全极化SAR影像，IB为变化后的全极化SAR影像，i表示所提取的特征序号；n为13，13个特征分别为：极化度，极化总功率，一致性参数，熵，二阶矩，对比度，相关性，极化熵，平均散射角，反熵，表面散射，二次散射和体散射；1)利用图像统计特性和协方差矩阵提取极化度，极化总功率，一致性参数三个特征；2)利用灰度共生矩阵提取熵、二阶矩、对比度、相关性四个特征；3)利用Cloude分解得到极化熵，平均散射角与反熵三个特征；4)利用Freeman分解得到表面散射、二次散射和体散射三个特征；上述选取的13个特征，对植被，建筑物，水体，裸地，道路这五类主要地物变化检测效果较好。(2)生成差异影像，根据步骤(1)中得到的特征影像，计算得到差异影像ri(i＝1,2,3,···,n)各像素点的像素值，即ri(s,t)＝Fi(s,t)/Ei(s,t)(s＝1,2,···,M；t＝1,2,···,N)其中ri表示第i个特征的差异影像，(s,t)表示差异影像上任意像素点，ri(s,t)为差异影像ri中像素点(s,t)的像素值，Fi(s,t)表示特征影像Fi中像素点(s,t)的像素值，Ei(s,t)表示特征影像Ei中像素点(s,t)的像素值。(3)根据步骤(2)中得到的差异影像ri，判定所有像素点为变化像素点、未变化像素点或待确定像素点；1)所有像素点均为判定对象，则任意像素点均可判定为三种可能事件：其中H1表示“像素点变化”，H2表示“像素点未变化”，H3表示“像素点待确定”；2)将每个特征的差异影像ri分别作为判别依据，求得所有像素点在每个判别依据下三种可能事件的信任程度，其计算公式如下：式中mi(H1)表示在第i个判别依据下，对“像素点变化”事件的信任程度；mi(H2)表示在第本文档来自技高网...

【技术保护点】
1.一种全极化SAR影像变化区域精确提取方法，其特征在于，包括以下步骤：/n(1)提取特征，对大小为M×N全极化SAR影像I

【技术特征摘要】
1.一种全极化SAR影像变化区域精确提取方法，其特征在于，包括以下步骤：
(1)提取特征，对大小为M×N全极化SAR影像IA提取n个特征，得到特征影像Fi(i＝1,2,3....n)，对大小为M×N全极化SAR影像IB同样提取n个特征，得到特征影像Ei(i＝1,2,3....n)；
其中IA为变化前的全极化SAR影像，IB为变化后的全极化SAR影像，i表示所提取的特征序号；
n为13，13个特征分别为：极化度，极化总功率，一致性参数，熵，二阶矩，对比度，相关性，极化熵，平均散射角，反熵，表面散射，二次散射和体散射；
(2)生成差异影像，根据步骤(1)中得到的特征影像，计算得到差异影像ri(i＝1,2,3,···,n)各像素点的像素值，即
ri(s,t)＝Fi(s,t)/Ei(s,t)(s＝1,2,···,M；t＝1,2,···,N)
其中ri表示第i个特征的差异影像，(s,t)表示差异影像上任意像素点，ri(s,t)为差异影像ri中像素点(s,t)的像素值，Fi(s,t)表示特征影像Fi中像素点(s,t)的像素值，Ei(s,t)表示特征影像Ei中像素点(s,t)的像素值；
(3)根据步骤(2)中得到的差异影像ri，判定所有像素点为变化像素点、未变化像素点或待确定像素点；
1)所有像素点均为判定对象，则任意像素点均可判定为三种可能事件：其中H1表示“像素点变化”，H2表示“像素点未变化”，H3表示“像素点待确定”；
2)将每个特征的差异影像ri分别作为判别依据，求得所有像素点在每个判别依据下三种可能事件的信任程度，其计算公式如下：

【专利技术属性】
技术研发人员：刘善伟，叶传龙，孙钦婷，赵吉祥，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东;37