python矢量转栅格

矢量转栅格

• 使用gdal和ogr
• 参考https://gis.stackexchange.com/questions/212795/rasterizing-shapefiles-with-gdal-and-python。这个答案得到的结果是上下翻转的，不能直接使用
• 本人修改后的代码
• gdal.RasterizeLayer的options一些注意的地方：

1. ALL_TOUCHED=TRUE 表示所有与矢量相交的像元都赋值
2. "ATTRIBUTE=%s"%field 表示栅格的值为field字段的值，如果不加这条表示矢量转换为一个值

import gdal
import ogr
import gdalconst
def shp2Raster(shp,templatePic,output,field