使用python中@中從CSV文件生成TIFF文件

Question

我有這種格式許多CSV文件：

Latitude,Longitude,Concentration 
53.833399,-122.825257,0.021957 
53.837893,-122.825238,0.022642 
.... 

我的目標是產生基於這些文件中的信息（支持GeoTiff文件每CSV一個TIFF文件文件），最好使用python。這是幾年前在我正在進行的項目上完成的，但是他們之前的做法是如何失敗的。我所知道的是他們最有可能使用GDAL。

我試圖通過研究如何使用GDAL來做到這一點，但由於資源有限並且我不知道如何使用它，所以這並沒有給我帶來任何幫助。

有人可以幫助我嗎？

Answer 1

這裏有一個我適應你的情況的代碼。您需要將GDAL目錄中的所有* .exe添加到您的路徑中才能正常工作（在大多數情況下，它是C:\Program Files (x86)\GDAL）。

它使用gdal_grid.exe UTIL（見DOC這裏：http://www.gdal.org/gdal_grid.html）

如你所願的gdal_cmd變量適合您的需求，您可以修改。

import subprocess 
import os 

# your directory with all your csv files in it 
dir_with_csvs = r"C:\my_csv_files" 

# make it the active directory 
os.chdir(dir_with_csvs) 

# function to get the csv filenames in the directory 
def find_csv_filenames(path_to_dir, suffix=".csv"): 
    filenames = os.listdir(path_to_dir) 
    return [ filename for filename in filenames if filename.endswith(suffix) ] 

# get the filenames 
csvfiles = find_csv_filenames(dir_with_csvs) 

# loop through each CSV file 
# for each CSV file, make an associated VRT file to be used with gdal_grid command 
# and then run the gdal_grid util in a subprocess instance 
for fn in csvfiles: 
    vrt_fn = fn.replace(".csv", ".vrt") 
    lyr_name = fn.replace('.csv', '') 
    out_tif = fn.replace('.csv', '.tiff') 
    with open(vrt_fn, 'w') as fn_vrt: 
     fn_vrt.write('<OGRVRTDataSource>\n') 
     fn_vrt.write('\t<OGRVRTLayer name="%s">\n' % lyr_name) 
     fn_vrt.write('\t\t<SrcDataSource>%s</SrcDataSource>\n' % fn) 
     fn_vrt.write('\t\t<GeometryType>wkbPoint</GeometryType>\n') 
     fn_vrt.write('\t\t<GeometryField encoding="PointFromColumns" x="Longitude" y="Latitude" z="Concentration"/>\n') 
     fn_vrt.write('\t</OGRVRTLayer>\n') 
     fn_vrt.write('</OGRVRTDataSource>\n') 

    gdal_cmd = 'gdal_grid -a invdist:power=2.0:smoothing=1.0 -zfield "Concentration" -of GTiff -ot Float64 -l %s %s %s' % (lyr_name, vrt_fn, out_tif) 

    subprocess.call(gdal_cmd, shell=True)