[發(fā)明專利]一種面向光場光譜數(shù)據(jù)亞像元定位的豐度修正方法及裝置

有效

• 申請?zhí)枺?/td>: 202010752242.1
• 申請日: 2020-07-30
• 公開（公告）號: CN111899300B
• 公開（公告）日: 2022-05-31
• 發(fā)明（設計）人: 蘇麗娟;袁艷;徐岳
• 申請（專利權）人: 北京航空航天大學
• 主分類號: G06T7/70
• 分類號: G06T7/70;G06V10/74;G06K9/62
• 代理公司: 北京清亦華知識產(chǎn)權代理事務所(普通合伙) 11201
• 代理人: 戴冬瑾

1.一種面向光場光譜數(shù)據(jù)亞像元定位的豐度修正方法，其特征在于，包括：

對待修正的混合像元，確定第k類地物類別對應的亞像元生長點的坐標值；

根據(jù)所述亞像元生長點的坐標值，計算所述混合像元內(nèi)每個亞像元到所述亞像元生長點的距離；

按照從小到大的順序?qū)λ龌旌舷裨獌?nèi)每個亞像元到所述亞像元生長點的距離進行排序，并根據(jù)排序結果依次將每個距離對應的亞像元在光場光譜成像條件下混合像元內(nèi)的豐度系數(shù)進行累加，當累加結果大于所述第k類地物類別在所述混合像元內(nèi)的豐度值時，停止累加，獲取參與累加的亞像元數(shù)目作為所述第k類地物類別對應的亞像元數(shù)量；

根據(jù)所述亞像元數(shù)量，利用亞像元定位算法進行亞像元定位。

2.如權利要求1所述的方法，其特征在于，所述對待修正的混合像元，確定第k類地物類別對應的亞像元生長點的坐標值，包括：

獲取所述混合像元的鄰域像元；

以所述混合像元的中心作為坐標軸原點，建立橫縱坐標系；

根據(jù)所述第k類地物類別在所述鄰域像元中的豐度值，計算得到所述第k類地物類別在所述橫縱坐標系中對應的亞像元生長點的橫坐標值x^k和縱坐標值y^k；

其中，采用如下公式計算得到所述亞像元生長點的橫坐標值x^k和縱坐標值y^k：

其中，S是重構尺度，表示第k類地物類別在鄰域像元P_a-1,b-1內(nèi)的豐度值。

3.如權利要求2所述的方法，其特征在于，所述根據(jù)所述亞像元生長點的坐標值，計算所述混合像元內(nèi)每個亞像元到所述亞像元生長點的距離，包括：

采用如下公式計算所述混合像元內(nèi)每個亞像元到所述亞像元生長點的距離：

其中，表示混合像元內(nèi)亞像元p_i,j到所述亞像元生長點的距離，x_i,j表示所述亞像元p_i,j的橫坐標值，y_i,j表示所述亞像元p_i,j的縱坐標值。

4.如權利要求3所述的方法，其特征在于，所述按照從小到大的順序?qū)λ龌旌舷裨獌?nèi)每個亞像元到所述亞像元生長點的距離進行排序，并根據(jù)排序結果依次將每個距離對應的亞像元在光場光譜成像條件下混合像元內(nèi)的豐度系數(shù)進行累加，包括：

采用如下公式將每個距離對應的亞像元在光場光譜成像條件下混合像元內(nèi)的豐度系數(shù)進行累加：

其中，表示第k個地物類別的豐度系數(shù)當前累加結果，表示豐度系數(shù)的上一累加結果，C_i,j表示(i,j)位置處的亞像元在光場光譜豐度數(shù)據(jù)中的貢獻系數(shù)。

5.一種面向光場光譜數(shù)據(jù)亞像元定位的豐度修正裝置，其特征在于，包括：

確定模塊，用于對待修正的混合像元，確定第k類地物類別對應的亞像元生長點的坐標值；

計算模塊，用于根據(jù)所述亞像元生長點的坐標值，計算所述混合像元內(nèi)每個亞像元到所述亞像元生長點的距離；

處理模塊，用于按照從小到大的順序?qū)λ龌旌舷裨獌?nèi)每個亞像元到所述亞像元生長點的距離進行排序，并根據(jù)排序結果依次將每個距離對應的亞像元在光場光譜成像條件下混合像元內(nèi)的豐度系數(shù)進行累加，當累加結果大于所述第k類地物類別在所述混合像元內(nèi)的豐度值時，停止累加，獲取參與累加的亞像元數(shù)目作為所述第k類地物類別對應的亞像元數(shù)量；

定位模塊，用于根據(jù)所述亞像元數(shù)量，利用亞像元定位算法進行亞像元定位。