MATLAB 地理データ処理 25: 植生季節学の抽出と分析モデルの最適化 (Savitzky-Golay)

1. 前提条件

以前、Savitzky-Golay を使用してリモート センシング データを処理し、地上生物季節学情報の MATLAB コードを取得する記事「Python 地理データ処理 17: 植生季節学の抽出と分析 (Savitzky-Golay)」を書きましたが、
いくつかの問題があることがわかりました。生物季節学の期間と生物季節学の終わりの閾値設定は同じ値ですが、一部の研究では、次のように環境条件の違いに応じて異なる季節学的閾値を設定する必要があります。

動的しきい値メソッドのモデルは次のとおりです。

2. MATLAB コード

SOS しきい値は 0.2 に設定され、EOS しきい値は 0.5 に設定されます。

% 读取文件夹下的子文件夹列表
folderPath = 'D:\dataset\phenology';
subFolders = dir(folderPath);
subFolders = subFolders([subFolders.isdir]);
subFolders = subFolders(3:end); % 排除"."和".."文件夹

% 循环处理每个子文件夹
for i = 1:length(subFolders)
	fprintf('d=%d\n', i)
    subFolder = subFolders(i).name;
    subsubFolder = subFolder(1:4);
    subFolderPath = fullfile(folderPath, subFolder);
    
    % 读取第一个tif影像获取基本信息
    tifFile = dir(fullfile(subFolderPath, '*.tif'));
    tifFilePath = fullfile(subFolderPath, tifFile(1).name);
    [a, R] = geotiffread(tifFilePath);
    info = geotiffinfo(tifFilePath);