[发明专利]一种多指数决策的浒苔混合像元分解方法

权利要求书

1.一种多指数决策的浒苔混合像元分解方法，其特征在于，包括以下步骤：

(1)获取浒苔研究区遥感影像Im，并对其进行预处理，得到处理后的影像I，大小为row×col；

其中，预处理步骤为：辐射定标、几何校正、图像裁剪、陆地掩膜和云掩膜；

(2)计算影像I的归一化植被指数NDVI、比值植被指数RVI、差值植被指数DVI和增强植被指数EVI，得到各植被指数的灰度图g_index，分别为g_NDVI、g_RVI、g_DVI和g_EVI，并采用最大类间方差法对所得灰度图进行阈值分割，得到各二值图G_index，分别为G_NDVI、G_RVI、G_DVI和G_EVI；

其中，G_index中各像元的取值为0或1，0代表该处为海水，1代表该处为浒苔；

(3)利用步骤(2)所得的各二值图G_index对I中各像元pixel(i,j)进行分类，将其归类至像元集合P，包括浒苔像元集合P_algae、海水像元集合P_seawater和待定像元集合P_candidate；

其中，(i,j)代表第i行的第j列，i＝1,2,…,row，j＝1,2,…,col；

(4)利用步骤(3)所得的浒苔像元集合P_algae和海水像元集合P_seawater分别获取浒苔端元集合E_algae和海水端元集合E_seawater，再分别获取用于浒苔丰度计算的浒苔端元DVI阈值DVI_algae和海水端元DVI阈值DVI_seawater；

(5)利用步骤(4)所得的浒苔端元DVI阈值DVI_algae和海水端元DVI阈值DVI_seawater计算P_candidate中各待定像元的浒苔丰度根据丰度值将其归类至P_algae或P_seawater以对P_algae和P_seawater中像元进行更新；

其中，各待定像元的浒苔丰度的计算公式如下：

式中，为像元pixel(i,j)的DVI值；

(6)利用步骤(5)所得的更新后的浒苔像元集合P_algae和海水像元集合P_seawater，根据步骤(4)分别对浒苔端元集合E_algae和海水端元集合E_seawater中端元进行更新，并分别利用更新后的E_algae和E_seawater对浒苔端元DVI阈值DVI_algae和海水端元DVI阈值DVI_seawater进行更新；

(7)利用步骤(6)所得的更新后的浒苔端元DVI阈值DVI_algae和海水端元DVI阈值DVI_seawater计算得到I中各像元pixel(i,j)的浒苔丰度

2.根据权利要求1所述的一种多指数决策的浒苔混合像元分解方法，其特征在于：所述步骤(3)中，各像元pixel(i,j)的分类过程如下：

式中，G_NDVI(i,j)、G_DVI(i,j)、G_EVI(i,j)和G_RVI(i,j)分别为G_NDVI、G_DVI、G_EVI和G_RVI(i,j)处像元值。