物候提取模型优化
1. 前提
之前我写过一篇使用Savitzky-Golay处理遥感数据,获取地面物候信息的MATLAB代码:Python地理数据处理 十七:植被物候提取和分析(Savitzky-Golay)
但是发现存在一些问题:对物候开始期和物候结束期的阈值设定都为相同值,但是有些研究中,需要根据不同的环境条件,设置不同的物候阈值,如下:
动态阈值法模型如下所示:
2. MATLAB代码
SOS阈值设置为0.2,EOS阈值设置为0.5:
% 读取文件夹下的子文件夹列表
folderPath = 'D:\dataset\phenology';
subFolders