森林结构参数提取

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

watershed.py
chm.tif