资 源 简 介
提取森林结构参数
#读取冠层高度模型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出图查看初步结果
