# 讀取TRMM影像
data = []
g = os.walk("C:/Users/KAI/Desktop/PYTHON作業/衛星降雨數據TRMM3B43_分辨率0.25度")
for path,d,filelist in g:
for filename in filelist:
if filename.endswith('TIF'):
str = (os.path.join(path, filename)).replace("\\", "/")
dataset = gdal.Open(str)
im_width = dataset.RasterXSize #柵格矩陣的列數
im_height = dataset.RasterYSize #柵格矩陣的行數
im_geotrans = dataset.GetGeoTransform() #仿射矩陣,左上角像素的大地座標和像素分辨率
im_proj = dataset.GetProjection() #地圖投影信息,字符串表示
im_data = dataset.ReadAsArray(0,0,im_width,im_height)
data.append(int(im_data))
print(np.array(data).shape)
或者
#---------------------------獲取目錄下的所有作業文件----------------------------
file_list = os.listdir('E:\\課程PPT\\earth system science\\assignment2')
ass_file = []
for i in file_list:
a,b = os.path.splitext(i)
if b == '.xlsx' or b == '.py' or b == '' :
continue
else:
ass_file.append(i)
好好生活,珍惜愛你的和你愛的人!
版權歸作者 小白是哪個小白_ 所有,轉載、引用請註明鏈接出處。