[發明專利]一種遙感影像變化檢測方法

有效

• 申請號: 202010147402.X
• 申請日: 2020-03-05
• 公開（公告）號: CN111340792B
• 公開（公告）日: 2022-04-12
• 發明（設計）人: 聶倩;趙艷福;陳立波;廖佳;張偉;七珂珂
• 申請（專利權）人: 寧波市測繪和遙感技術研究院
• 主分類號: G06T7/00
• 分類號: G06T7/00;G06T7/33;G06T7/13;G06T7/11;G06K9/62;G06V10/762
• 代理公司: 寧波誠源專利事務所有限公司 33102
• 代理人: 陳洪娜;方寧

1.一種遙感影像變化檢測方法，其特征在于：包括以下步驟：

步驟1、輸入同一區域、不同時相的兩幅多光譜遙感影像，分別記為遙感影像X和遙感影像Y；

步驟2、對遙感影像X和遙感影像Y進行影像配準，得到影像配準后的遙感影像X和遙感影像Y；

步驟3、對影像配準后的遙感影像X和遙感影像Y進行相對輻射歸一化校正，得到遙感影像X’和遙感影像Y’；

步驟4、將遙感影像X’和遙感影像Y’中的前時相遙感影像作為參考影像，后時相遙感影像作為待檢測影像，使用迭代加權多元變化檢測算法IR-MAD確定出待檢測影像相對于參考影像中的變化像元，并將所有變化像元組成變化差異圖；

步驟5、利用均值漂移算法對變化差異圖進行聚類處理，得到變化檢測圖像；

具體步驟為：

步驟5-1、在變化差異圖中任意選取一個像元x，選取以該像元x為中心的一個窗口；

步驟5-2、計算該像元x的均值漂移向量m_h(x)；計算公式為：

其中，x_i是以像元x為中心的窗口中的第i個像元點，i＝1、2、3…n，n為以像元x為中心的窗口中所有像元點的總數；‖‖為求歐式距離；為計算k(y)的梯度，k(y)為高斯核函數，h為高斯核函數的帶寬；

步驟5-3、判斷|m_h(x)-x|ε是否成立，ε為預設的誤差閾值，如成立，則該像元x為收斂點，并轉至步驟5-4；如不成立，則使x＝m_h(x)，并轉至步驟5-2，將更新后的x作為新的像元重新計算出收斂點；

步驟5-4、按照步驟5-1～步驟5-3中相同的方法依次計算出變化差異圖中每個像元的收斂點；

步驟5-5、將具有相同收斂點的像元點歸為同一類，即將歸為同一類的像元點劃為同一個分割子區域，得到分割后的圖像，該圖像為變化檢測圖像；

步驟6、采用形態學腐蝕和膨脹算法對變化檢測圖像進行圖像處理，得到最終的變化檢測圖像。

2.根據權利要求1所述的遙感影像變化檢測方法，其特征在于：所述步驟3中相對輻射歸一化校正使用的方法為偽不變特征法。

3.根據權利要求1所述的遙感影像變化檢測方法，其特征在于：所述步驟6中對變化檢測圖像進行圖像處理的步驟為：先對變化檢測圖像進行腐蝕運算，再對腐蝕之后得到的圖像進行膨脹運算，該膨脹運算之后得到的圖像為最終的變化檢測圖像。