# 读取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)
好好生活,珍惜爱你的和你爱的人!
版权归作者 小白是哪个小白_ 所有,转载、引用请注明链接出处。