[发明专利]沙漠景观斑块分割方法及系统、电子设备及可读存储介质

权利要求书

1.一种沙漠景观斑块分割方法，其特征在于，包括：

S1：对输入的沙漠地区遥感影像分割得到光谱均质对象并计算出类相关格局特征，基于所述光谱均质对象与类相关格局特征，得到初始分割对象；

S2：基于所述初始分割对象的归一化植被指数与地形特征提取的种子点，分别得到绿洲斑块与缓起伏沙地斑块；将所述初始分割对象中的其他的未增长区域作为剩余初始对象；

S3：将所述剩余初始对象基于线格局特征拓展MRS，得到沙漠景观斑块。

2.根据权利要求1所述的沙漠景观斑块分割方法，其特征在于，所述S1具体包括：根据所述沙漠地区遥感影像的影像光谱特征与形状特征，利用MRS对所述沙漠地区遥感影像进行分割，得到光谱均质对象，利用ISODATA方法对沙漠地区遥感影像无监督聚类，计算出类相关格局特征，利用MRS对所述类相关格局特征与光谱均质对象进行合并，得到初始分割对象。

3.根据权利要求1所述的沙漠景观斑块分割方法，其特征在于，所述S2具体包括：

S21：选取初始分割对象中未被增长、且NDVI值大于第一设定阈值的作为绿洲种子点，若绿洲种子点与其相邻对象的NDVI绝对差值最小值小于第二给定阈值，则将绿洲种子点与具有最小绝对差值的第一相邻对象合并；

S22：判断所述初始分割对象中是否存在绿洲种子点或存在满足合并条件的第一相邻对象，若是，则执行S21，若否，则将得到的增长区域作为绿洲斑块并执行S23；

S23：选取初始分割对象中未被增长、且地形特征均值小于第三给定阈值的作为缓起伏沙地种子点，若缓起伏沙地种子点与其相邻对象的地形特征绝对差值最小值小于第四给定阈值，则将缓起伏沙地种子点与具有最小绝对差值的第二相邻对象合并，其中，所述地形特征包括坡度、起伏度和切割深度；

S24：判断所述初始分割对象中是否存在缓起伏沙地种子点或存在满足合并条件的第二相邻对象，若是，则执行B21，若否，则将得到的增长区域作为缓起伏沙地斑块并执行S25；

S25：将所述初始分割对象中的其他的未增长区域作为剩余初始对象。

4.根据权利要求1所述的沙漠景观斑块分割方法，其特征在于，所述S3具体包括：

S31：利用基于置信度的边缘检测法对数字高程模型进行沙丘走势线提取，矢量化边缘并进行后处理以去除边缘噪声，得到刻画沙丘基本形态的走势线；定义沙丘长度、沙丘最小间距、沙丘密度、沙丘方向、沙丘高度与沙丘弯曲度作为线格局特征，拉伸所有线格局特征至相同数值区间；

S32：定义与每个剩余初始对象相交的所有走势线为所述剩余初始对象的关联线集合，基于拓展MRS实现地理场景分割。

5.根据权利要求4所述的沙漠景观斑块分割方法，其特征在于，所述S32的具体实现方法包括：

S321：当O₁与O₂之间的异质度增量是O₁与其所有相邻对象的异质度增量的最小值、且小于第五给定阈值时，将O₁与O₂合并；

S322：判断所述剩余初始对象中是否存在满足合并条件的相邻对象，若是，则执行S321，若否，则得到的增长区域作为沙漠景观斑块。

其中，所述O₁与O₂之间的异质度增量的计算方法为：假设判断对象O₁与其相邻对象O₂是否合并为O₃，其中，O₁、O₂与O₃的关联线集合分别为L₁、L₂与L₃，则O₁与O₂合并前后异质度的变化即为异质度增量，可表示为：

f＝w_lsp*Δh_lsp+(1-w_lsp)*Δh_shape

式中，f为合并前后的异质度增量，由线格局异质度增量Δh_lsp与形状异质度增量Δh_shape加权构成，w_lsp为Δh_lsp所占权重。