geotransform,projection,noDataValue = Readxy('F://hourly_maps_raster//liang//2018-01-01_00.tif')
#获取源文件的行...filepath = os.path.join(dirpath,filename)
purename = filename.replace('.tif','') #获得除去扩展名的文件名...= np.array(filedata)
filedata = Read(filepath) #将2010年的13幅图像数据存入filedata中..._ == "__main__":
rows,cols,geotransform,projection,noDataValue = Readxy('D:/minxinan/o3idw/0.tif')
#获取源文件的行...(filepath) #将2010年的13幅图像数据存入filedata中
count+=1
np.add(filesum,filedata,filesum) #求13幅图像相应栅格值的和