import os import glob in_pathfiles = "C:\Data\H5_Data_LST\*" out_pathfiles = "C:\Data\Out_Data_LST\\" myfiles = glob.glob(in_pathfiles) #list of all files for f in myfiles: print(f) filename = f.split("\\")[-1] f_out = out_pathfiles + filename + ".tiff" f_rep = out_pathfiles + filename + "_rep.tiff" os.system('gdal_translate -of GTiff -a_srs \ "+proj=geos +h=35785831 +a=6378169 +b=6356583.8 +no_defs" \ -a_ullr -5568748.27576 5568748.27576 5568748.27576 -5568748.27576 \ "HDF5:'+ f + '://LST" '+ f_out) os.system('gdalwarp -ot Float32 -s_srs "+proj=geos +h=35785831 +a=6378169 +b=6356583.8 +no_defs" -t_srs EPSG:4326 -r near -of GTiff ' + f_out + ' ' + f_rep)