国家青藏高原科学数据中心
全球逐日0.05°时空连续地表温度数据集(2002-2022)
此代码仅用于该数据集处理
matlab编程批量处理hdf5格式转换为tif格式,并提取选定的一个子集
版本:arcgis10.2 matlab2020
参考:MATLAB hdf(h5)文件转成tif图片格式(+批量处理)
一、遇到问题
h5格式,说明文档里面有写空间位置。但是打开后没有设置投影信息,在arcgis里打开显示的范围特别大,影像是旋转过的还镜像的。定义投影后还是没有同矢量数据重叠,不知道要怎么处理
二、打开matlab2020
注意文件名不能和变量名相同,否则会报错文章来源:https://www.toymoban.com/news/detail-785569.html
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里打开 成功!
文章来源地址https://www.toymoban.com/news/detail-785569.html
到了这里,关于(Arcgis)matlab编程批量处理hdf5格式转换为tif格式的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!