(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文件fori=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']);% 构建完整的输出文件路径% 输出为GeoTIFFgeotiffwrite(output_file,data,R);% 将数据和地理空间参考对象写入输出文件(GeoTIFF格式)end