森林结构参数提取
于 2022-02-13 发布
文件大小:630.10 kB
0 140
下载积分: 2
下载次数: 1
代码说明：
提取森林结构参数 #读取冠层高度模型CHM raster = gdal.Open("./chm.tif") banddataraster = raster.GetRasterBand(1) dataraster = banddataraster.ReadAsArray() #对CHM进行高斯滤波，平滑数据 dataraster_gau = gaussian_filter(dataraster,sigma=1) #寻找CHM中的局部最大值作为分水岭的标记点 local_maxi = peak_local_max(dataraster_gau, indices=False) markers = ndi.label(local_maxi)[0] #利用分水岭算法进行分割，labels变量中存储了每个分割结果 labels = watershed(-dataraster_gau, markers, mask= dataraster_gau[:]>5) #利用matplotlib出图查看初步结果
下载说明：请别用迅雷下载，失败请重下，重下不扣分！
发表评论
