发明名称 基于人工免疫多目标聚类的SAR图像变化检测方法
摘要 本发明公开了一种基于人工免疫多目标聚类的SAR图像变化检测方法,主要解决SAR图像变化检测结果准确率及效率低的问题。其实现步骤是:(1)读入两时相SAR图像;(2)对两时相SAR图像构造差异图;(3)对差异图进行基于灰度值的自适应免疫多目标聚类,将其分为变化类、非变化类与待识别类;(4)对待识别类进行基于非下采样小波变换的免疫克隆多目标聚类,得到待识别类的一组聚类中心;(5)根据一组聚类中心对待识别类进行最小距离分类,得到一组变化检测结果图;(6)计算变化检测结果图的目标函数值;(7)根据目标函数值选出最小的目标函数值;(8)将最小的目标函数值所对应的变化检测结果图作为最终的检测结果。本发明方法具有检测效率高及检测精度高的优点。
申请公布号 CN103366365B 申请公布日期 2016.05.25
申请号 CN201310241979.7 申请日期 2013.06.18
申请人 西安电子科技大学 发明人 尚荣华;焦李成;齐丽萍;王爽;吴建设;公茂果;李阳阳;马晶晶;马文萍
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于人工免疫多目标聚类的SAR图像变化检测方法,包括如下步骤:(1)读入已配准和已校正的两幅不同时间同一地点的SAR图像<img file="FDA0000932426790000011.GIF" wi="55" he="79" />和<img file="FDA0000932426790000012.GIF" wi="86" he="78" />(2)根据这两幅SAR图像<img file="FDA0000932426790000013.GIF" wi="53" he="77" />和<img file="FDA0000932426790000014.GIF" wi="55" he="78" />构造差异图D:<maths num="0001"><math><![CDATA[<mrow><mi>D</mi><mo>=</mo><mfrac><mrow><munder><mo>&Sigma;</mo><mrow><mi>x</mi><mo>&Element;</mo><mi>N</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></munder><mi>m</mi><mi>i</mi><mi>n</mi><mo>(</mo><msub><mi>I</mi><msub><mi>t</mi><mn>1</mn></msub></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>I</mi><msub><mi>t</mi><mn>2</mn></msub></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>)</mo></mrow><mrow><munder><mo>&Sigma;</mo><mrow><mi>x</mi><mo>&Element;</mo><mi>N</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></munder><mi>m</mi><mi>a</mi><mi>x</mi><mo>(</mo><msub><mi>I</mi><msub><mi>t</mi><mn>1</mn></msub></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>I</mi><msub><mi>t</mi><mn>2</mn></msub></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>)</mo></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000932426790000015.GIF" wi="541" he="213" /></maths>其中,N(x)表示像素点x周围的一个邻域像素集,<img file="FDA0000932426790000016.GIF" wi="113" he="71" />是SAR图像<img file="FDA0000932426790000017.GIF" wi="52" he="69" />中像素点x的灰度值,<img file="FDA0000932426790000018.GIF" wi="118" he="70" />是SAR图像<img file="FDA0000932426790000019.GIF" wi="55" he="70" />中像素点x的灰度值;(3)对差异图D进行基于灰度值的自适应免疫多目标聚类,自适应地将差异图D分为变化类B、非变化类F与待识别类S;(4)对待识别类S进行基于非下采样小波变换的免疫克隆多目标聚类,得到待识别类S的一组聚类中心E;(5)根据一组聚类中心E对待识别类S进行最小距离分类,得到一组变化检测结果图G,G={g′<sub>l</sub>,l=1,...,a},a是变化检测结果图G的总个数,g′<sub>l</sub>是变化检测结果图G中的第l幅变化检测结果图;(6)根据变化检测结果图G={g′<sub>l</sub>,l=1,...,a},计算目标函数值F,F={f′<sub>l</sub>,l=1,...,a},f′<sub>l</sub>是第l幅变化检测结果图g′<sub>l</sub>的目标函数值,表示为:<maths num="0002"><math><![CDATA[<mrow><msup><msub><mi>f</mi><mi>l</mi></msub><mo>&prime;</mo></msup><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>r</mi><mo>=</mo><mn>0</mn></mrow><mn>1</mn></munderover><mfrac><msub><mi>N</mi><mi>r</mi></msub><mrow><mi>M</mi><mo>&times;</mo><mi>N</mi></mrow></mfrac><munder><mo>&Sigma;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>&Element;</mo><msub><mi>R</mi><mi>r</mi></msub></mrow></munder><msup><mrow><mo>(</mo><mi>D</mi><mo>(</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mo>)</mo><mo>-</mo><msub><mi>&mu;</mi><mi>r</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow>]]></math><img file="FDA00009324267900000110.GIF" wi="668" he="135" /></maths>其中,M×N是第l幅变化检测结果图g'<sub>l</sub>的图像大小,D(i,j)是差异图D中第i行第j列的像素点的灰度值,r=0,1,当r=0时,N<sub>r</sub>是第l幅变化检测结果图g'<sub>l</sub>中非变化类像素点的总个数,R<sub>r</sub>是第l幅变化检测结果图g′<sub>l</sub>中非变化类像素点的集合,μ<sub>r</sub>是第l幅变化检测结果图g′<sub>l</sub>中非变化类像素点的灰度均值,当r=1时,N<sub>r</sub>是第l幅变化检测结果图g'<sub>l</sub>中变化类像素点的总个数,R<sub>r</sub>是第l幅变化检测结果图g'<sub>l</sub>中变化类像素点的集合,μ<sub>r</sub>是第l幅变化检测结果图g′<sub>l</sub>中变化类像素点的灰度均值;(7)根据目标函数值F={f′<sub>l</sub>,l=1,...,a},选出最小的目标函数值f′<sub>z</sub>,将最小的目标函数值f′<sub>z</sub>所对应的第z幅变化检测结果图g'<sub>z</sub>作为最终的检测结果。
地址 710071 陕西省西安市太白南路2号