Python：多个tiff转csv，并按照像元合并成一个csv

1）定义tiff转csv的函数

from osgeo import gdal import pandas as pd import numpy as np  gdal.AllRegister()  def tifftocsv(path):     filePath = path     dataset = gdal.Open(filePath)     adfGeoTransform = dataset.GetGeoTransform()     nXSize = dataset.RasterXSize  # 列数     nYSize = dataset.RasterYSize  # 行数     im_data = dataset.ReadAsArray(0,0,nXSize,nYSize)     index = []  # 纬度     columns = []  # 经度     for j in range(nYSize):         lat = adfGeoTransform[3] + j * adfGeoTransform[5]         index.append(lat)     for i in range(nXSize):         lon = adfGeoTransform[0] + i * adfGeoTransform[1]         columns.append(lon)     data = pd.DataFrame(im_data, index=index, columns=columns)     data_np=np.array(data)     data_np=data_np.reshape(5254,1)     data=pd.DataFrame(data_np)     return data    

2)计算得到一个DataFrame（作为后续合并的基础文件）

 aod_ab_num=6  for i in range(1,aod_ab_num+1):      AODpath='D:\\research\\5_AOD\\results\\AODresults\\nmf'+str(aod_ab_num)+'\\Abundance'+str(i)+'.tif'     AOD=tifftocsv(AODpath)     AOD.columns=['AOD'+str(i)]     if(i==1):         result=AOD     else:         frames=[result,AOD]         result = pd.concat(frames,axis=1)  

3）批量计算其余所有的变量

 VARlist=['NDVI','VOC','SO2','Precipitation','PM2.5','OC','NOX','LST','EVI','BC']  ab_num=6  for i in range(0,len(VARlist)):     var=str(VARlist[i])     print(var)     for j in range(1,ab_num+1):        # print(var)         #print(j)         VARpath='D:\\research\\5_AOD\\results\\'+var+'results\\nmf'+str(ab_num)+'\\Abundance'+str(j)+'.tif'        # print(VARpath)         VAR=tifftocsv(VARpath)         VAR.columns=[var+str(j)]         frame=[result,VAR]         result=pd.concat(frame,axis=1)         result=result.fillna(0)