代码分享-NC多时相数据利用Matlab代码提取

NC多时相数据利用Matlab代码提取

ncFileBasePath = ‘’; % nc 文件的基础路径
tiffOutFileBasePath = ‘’; % tif 文件的输出路径
kenyaLatlim = [-4 4];
kenyaLonlim = [33 42];
% 获取 nc 文件列表>> % 获取 nc 文件列表
ncFiles = dir(fullfile(ncFileBasePath, ‘*.nc’));
% 循环处理每个 nc 文件
for fileIndex = 1:length(ncFiles)
% 完整文件名
filename = ncFiles(fileIndex).name;
% 不带后缀的文件名
filenameWithoutSuffix = filename(1:find(filename == ‘.’) - 1);
% nc 文件的完整路径
fileFullPath = fullfile(ncFileBasePath, filename);
% 读取 SPEI 数据
ETaSets = ncread(fileFullPath, ‘spei’);
% 获取数据集的大小信息
sizeInfo = size(ETaSets);% 循环处理每个子集
for subsetIndex = 1:sizeInfo(3)% 获取单个子集数据
MonthData = ETaSets(:, :, subsetIndex);% 将 NaN 值替换为特定值
MonthData(isnan(MonthData)) = -8888;% 创建 tif 文件的输出路径和文件名
tifOutputFullPath = fullfile(tiffOutFileBasePath, …
[filenameWithoutSuffix, sprintf(‘%02d’, subsetIndex), ‘.tif’]);% 对数据进行处理和转换
MonthData = rot90(MonthData); % 旋转数据
Refference = georasterref(‘RasterSize’, size(MonthData), …
‘Latlim’, kenyaLatlim, ‘Lonlim’, kenyaLonlim); % 创建参考信息
Refference.ColumnsStartFrom = ‘north’; % 指定列的起始位置
geotiffwrite(tifOutputFullPath, MonthData, Refference); % 写入为 geotiff 文件
end
end