[发明专利]一种火后森林砍伐分布识别的制图方法及装置

权利要求书

1.一种火后森林砍伐分布识别的制图方法，其特征在于，包括：

获取并基于目标区域的遥感影像数据，计算所述遥感影像数据对应的燃烧指数、穗帽变换亮度指数和穗帽变换湿度指数；

设置影像合成的起始日期范围，获取多个年份中所述起始日期范围内的所有目标区域遥感影像，并对所述所有目标区域遥感影像进行影像合成，得到长时间序列年度无云最优影像；

基于预设的燃烧指数阈值和差分燃烧指数阈值，从所述长时间序列年度无云最优影像中提取出未过火森林区域和过火森林区域，对所述未过火森林区域进行像元筛选，得到所述长时间序列年度无云最优影像中的成熟森林像元；

逐年获取并基于所述成熟森林像元的穗帽变换亮度指数和穗帽变换湿度指数对应的均值及标准差，计算不同年份中所述过火森林区域中每个像元对应的标准穗帽变换亮度指数和标准穗帽变换湿度指数，并根据所述标准穗帽变换亮度指数和所述标准穗帽变换湿度指数，得到所述过火森林区域中每个像元对应的森林扰动指数；

获取火后未砍伐样区和火后砍伐样区，对所述火后未砍伐样区和所述火后砍伐样区进行分层随机抽样处理，得到火后未砍伐样区样本像元和火后砍伐样区样本像元，计算并基于每个火后未砍伐样区样本像元和每个火后砍伐样区样本像元对应的样本森林扰动指数，设置发生砍伐森林扰动指数阈值；

根据所述森林扰动指数和所述发生砍伐森林扰动指数阈值，逐年对所述过火森林区域中每个像元进行分类，得到不同年份对应的初始火后砍伐森林区域，并生成初始火后森林砍伐时空分布图像。

2.如权利要求1所述的一种火后森林砍伐分布识别的制图方法，其特征在于，生成初始火后森林砍伐时空分布图像后，还包括：

基于图像差分法计算相邻年份所有火后未砍伐样区样本像元和所述火后砍伐样区样本像元的差分森林扰动指数，根据所述差分森林扰动指数，设置差分扰动指数阈值；

基于所述过火森林区域中每个像元对应的森林扰动指数，计算相邻年份中所述过火森林区域中每个像元对应的差分森林扰动指数；

根据所述发生砍伐森林扰动指数阈值和所述差分森林扰动指数，逐年对所述过火森林区域中每个像元进行分类，得到火后森林初次砍伐区域及火后森林二次砍伐区域，并对应生成火后森林初次砍伐时空分布图像和火后森林二次砍伐时空分布图像。

3.如权利要求1所述的一种火后森林砍伐分布识别的制图方法，其特征在于，得到不同年份对应的初始火后砍伐森林区域后，还包括：

将所述初始火后砍伐森林区域中的每个像元设置为火后砍伐森林像元，获取当前年份中所有火后砍伐森林像元对应的第一森林扰动指数，同时获取预设年份中所有火后砍伐森林像元对应的第二森林扰动指数；

判断所述第一森林扰动指数和所述第二森林扰动指数的扰动差值是否在预设差值范围内，若否，则认为所述火后砍伐森林像元分类错误，并将所述火后砍伐森林像元修正为火后未砍伐森林像元，若是，则认为所述火后砍伐森林像元分类正确。

4.如权利要求1所述的一种火后森林砍伐分布识别的制图方法，其特征在于，基于预设的燃烧指数阈值和差分燃烧指数阈值，从所述长时间序列年度无云最优影像中提取出未过火森林区域和过火森林区域，具体包括：

获取所述长时间序列年度无云最优影像中每个像元对应的燃烧指数，将所述燃烧指数与所述燃烧指数阈值进行对比；

若所述燃烧指数大于所述燃烧指数阈值，则将所述像元设置为森林像元，整合所有森林像元，得到森林区域，若所述燃烧指数小于或等于所述燃烧指数阈值，则将所述像元设置为非森林像元，整合所有非森林像元，得到非森林区域；

计算相邻年份中所述森林区域的每个像元对应的差分燃烧指数，将所述差分燃烧指数与所述差分燃烧指数阈值进行对比；

若所述差分燃烧指数小于或等于所述差分燃烧指数阈值，则将所述森林像元设置为未过火森林像元，整合所有未过火森林像元，得到未过火森林区域，若所述燃烧指数大于或等于所述燃烧指数阈值，则将所述森林像元设置为过火森林像元，整合所有过火森林像元，得到过火森林区域。