國(guó)家青藏高原科學(xué)數(shù)據(jù)中心
中國(guó)區(qū)域1km無(wú)縫地表溫度數(shù)據(jù)集(2002-2020)
此代碼僅用于該數(shù)據(jù)集處理
matlab編程批量處理hdf4格式轉(zhuǎn)換為tif格式,并提取選定的一個(gè)子集
版本:arcgis10.2 matlab2020
參考:MATLAB hdf(h5)文件轉(zhuǎn)成tif圖片格式(+批量處理)
一、遇到問(wèn)題
hdf格式,說(shuō)明文檔里面有寫空間位置。但是打開(kāi)后沒(méi)有設(shè)置投影信息,在arcgis里打開(kāi)顯示的范圍特別大,影像是旋轉(zhuǎn)過(guò)的還鏡像的。定義投影后還是沒(méi)有同矢量數(shù)據(jù)重疊,不知道要怎么處理
二、打開(kāi)matlab2020
注意文件名不能和變量名相同,否則會(huì)報(bào)錯(cuò)文章來(lái)源:http://www.zghlxwxcb.cn/news/detail-785567.html
clc % 清除命令窗口內(nèi)容
clear % 清除工作空間變量
close all % 關(guān)閉所有圖形窗口
% 輸入文件夾和輸出文件夾路徑
input_folder = 'C:\Users\53446\Desktop\1\hdf'; % 輸入文件夾路徑
output_folder = 'C:\Users\53446\Desktop\1\tif'; % 輸出文件夾路徑
% 獲取輸入文件夾下所有文件的信息
files = dir(fullfile(input_folder, '*.hdf')); % 獲取輸入文件夾下所有.hdf文件的信息
% 遍歷每個(gè)文件并執(zhí)行操作
for i = 1:length(files)
% 讀取數(shù)據(jù)
filename = fullfile(input_folder, files(i).name); % 構(gòu)建完整的輸入文件路徑
data = hdfread(filename, '/Lst'); % 從.hdf文件中讀取/Lst數(shù)據(jù)集
data = permute(data,[2,1]); % 轉(zhuǎn)置數(shù)據(jù),此數(shù)據(jù)集需要
data = flipud(data); % 上下翻轉(zhuǎn)數(shù)據(jù),此數(shù)據(jù)集需要
% 創(chuàng)建空間參考對(duì)象
R = georasterref('RasterSize', size(data), 'Latlim', [double(18.33) double(53.49)], 'Lonlim', [double(73.55) double(134.99)]); % 創(chuàng)建一個(gè)空間參考對(duì)象,指定了柵格大小和經(jīng)緯度范圍
% 寫入輸出文件
output_filename = fullfile(output_folder, [files(i).name, '.tif']); % 構(gòu)建完整的輸出文件路徑
geotiffwrite(output_filename, data, R); % 將數(shù)據(jù)和空間參考對(duì)象寫入輸出文件
end
三、arcgis打開(kāi) 成功!
文章來(lái)源地址http://www.zghlxwxcb.cn/news/detail-785567.html
到了這里,關(guān)于(Arcgis)matlab編程批量處理hdf4格式轉(zhuǎn)換為tif格式的文章就介紹完了。如果您還想了解更多內(nèi)容,請(qǐng)?jiān)谟疑辖撬阉鱐OY模板網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關(guān)文章,希望大家以后多多支持TOY模板網(wǎng)!