（Arcgis）matlab编程批量处理hdf5格式转换为tif格式

国家青藏高原科学数据中心
全球逐日0.05°时空连续地表温度数据集（2002-2022）

此代码仅用于该数据集处理
matlab编程批量处理hdf5格式转换为tif格式，并提取选定的一个子集
版本：arcgis10.2 matlab2020
参考：MATLAB hdf(h5)文件转成tif图片格式（+批量处理）

一、遇到问题

h5格式，说明文档里面有写空间位置。但是打开后没有设置投影信息，在arcgis里打开显示的范围特别大，影像是旋转过的还镜像的。定义投影后还是没有同矢量数据重叠，不知道要怎么处理

二、打开matlab2020

注意文件名不能和变量名相同，否则会报错

clc % 清除命令窗口内容
clear % 清除工作空间变量
close all % 关闭所有图形窗口

% HDF5文件所在的文件夹路径
input_folder = 'C:\Users\53446\Desktop\1\h5';

% GeoTIFF文件输出文件夹路径
output_folder = 'C:\Users\53446\Desktop\1\tif';

% 创建输出文件夹
mkdir(output_folder); % 在指定路径下创建文件夹

% 获取输入文件夹中所有的HDF5文件
file_list = dir(fullfile(input_folder, '*.h5')); % 获取输入文件夹路径下所有.h5文件的信息

% 循环处理每个HDF5文件
for i = 1:length(file_list)

    input_file = fullfile(input_folder, file_list(i).name); % 构建完整的输入文件路径
    
    h5info(input_file); % 显示HDF5文件的结构信息
    
    h5disp(input_file, '/'); % 显示HDF5文件中所有group的信息
    
    data = h5read(input_file, '/LST_Day_CMG'); % 从HDF5文件中读取/LST_Day_CMG数据集
    
    % 对矩阵进行处理
    data = permute(data, [1, 2]); % 转置数据，此数据集需要
    data = flipud(data); % 上下翻转数据，此数据集需要
    
    R = georasterref('RasterSize', size(data), 'Latlim', [double(-90) double(90)], 'Lonlim', [double(-180) double(180)]); % 创建一个地理空间参考对象，指定了栅格大小和经纬度范围
    
    % 构建输出文件路径
    [~, filename, ~] = fileparts(file_list(i).name); % 获取文件名（不包括扩展名）
    output_file = fullfile(output_folder, [filename, '.tif']); % 构建完整的输出文件路径
    
    % 输出为GeoTIFF
    geotiffwrite(output_file, data, R); % 将数据和地理空间参考对象写入输出文件（GeoTIFF格式）
end

三、arcgis里打开 成功！