发明名称 基于多源目标特征支持的高光谱遥感图像变化检测方法
摘要 本发明公开一种基于多源目标特征和机器学的高光谱遥感图像变化检测方法,包括:可见光遥感图像变化检测;可见光与高光谱遥感图像配准,将可见光变化检测得到的不变化目标变换至高光谱遥感图像;基于变换至高光谱遥感图像中的不变化目标,选择训练样本;高光谱子空间学;在学得到的子空间内进行高光谱遥感图像变化检测。本发明的方法基于多源目标特征,通过可见光遥感图像确定不变化的高光谱样本用以学高光谱最佳判别子空间,在该子空间中进行高光谱遥感图像变化检测,性能稳定,适用范围广。
申请公布号 CN102819838B 申请公布日期 2015.03.25
申请号 CN201210247785.3 申请日期 2012.07.17
申请人 北京市遥感信息研究所;中国科学院自动化研究所 发明人 江碧涛;樊彬;张秀玲;余晓刚;潘春洪;李京龙;胡世仓
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 中科专利商标代理有限责任公司 11021 代理人 宋焰琴
主权项 一种基于多源目标特征和机器学习的高光谱遥感图像变化检测方法,其中输入数据为同一时间同一地点的可见光和高光谱遥感图像以及另一时间同一地点的可见光和高光谱遥感图像,该方法包括以下步骤:步骤S1:对可见光遥感图像进行变化检测,得到初始的检测结果,即区分出变化目标和不变化目标;步骤S2:通过可见光与高光谱遥感图像配准,将步骤S1中可见光遥感图像变化检测得到的不变化目标变换至高光谱遥感图像;步骤S3:基于可见光遥感图像变化检测的结果,选择置信度高,即变化概率低的不变化目标在高光谱遥感图像上对应的高光谱特征作为训练样本;步骤S4:根据步骤S3得到的训练样本集合,进行高光谱子空间学习;步骤S5:在高光谱子空间中对高光谱遥感图像进行变化检测;其中,在步骤S2,采用可见光与高光谱遥感图像配准的方法,将可见光遥感图像变化检测得到的不变化目标变换至高光谱遥感图像,得到高光谱遥感图像中对应的不变像素;在步骤S4,基于可见光遥感图像变化检测以及可见光与高光谱遥感图像配准的结果,利用可见光遥感图像上变化概率最低的20%不变化目标在对应高光谱遥感图像上的高光谱特征,进行子空间学习;在步骤S4,子空间学习的目的是使得进行子空间投影后,训练样本之间的距离尽可能小,而同时使得所有训练样本的方差尽可能大;步骤S4中,进行子空间学习所定义的目标函数为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>min</mi><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msup><mrow><mo>(</mo><msup><mi>W</mi><mi>T</mi></msup><mrow><mo>(</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><mrow><mo>(</mo><msup><mrow><mo>(</mo><msup><mi>W</mi><mi>T</mi></msup><msub><mi>X</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msup><mi>W</mi><mi>T</mi></msup><msub><mi>Y</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mi>min</mi><mfrac><mrow><msup><mi>W</mi><mi>T</mi></msup><mi>BW</mi></mrow><mrow><msup><mi>W</mi><mi>T</mi></msup><mi>AW</mi></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000585401600000021.GIF" wi="941" he="182" /></maths>其中<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>B</mi><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><mrow><mo>(</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>i</mi></msub><mo>)</mo></mrow><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000585401600000022.GIF" wi="585" he="115" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><mrow><mo>(</mo><msub><mi>X</mi><mi>i</mi></msub><msubsup><mi>X</mi><mi>i</mi><mi>T</mi></msubsup></mrow><mrow><mo>+</mo><mrow><msub><mi>Y</mi><mi>i</mi></msub><msubsup><mi>Y</mi><mi>i</mi><mi>T</mi></msubsup><mo>)</mo></mrow></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000585401600000023.GIF" wi="516" he="118" /></maths>N为训练样本的个数,X<sub>i</sub>是第i个样本在变化前高光谱遥感图像上的高光谱特征,Y<sub>i</sub>是第i个样本在变化后高光谱遥感图像上的高光谱特征,W为待求解的投影矩阵,通过拉格朗日乘数法,求解该目标函数得到:AW=λBW,这是一个广义特征值问题,λ是特征值,W是对应的特征向量,满足该目标函数的最优嵌入子空间由前K个最大的广义特征值所对应的特征向量组成,其中K是子空间的维数。
地址 100192 北京市海淀区清河小营东路2号院