gdal的矢量栅格化接口GDALRasterizeLayers使用(一)
作者：朱金灿
来源：http://blog.csdn.net/clever101
gdal库提供了一个矢量栅格化接口函数GDALRasterizeLayers。今天让我们学习一下这个接口函数的使用。该函数的原型如下：
CPLErr GDALRasterizeLayers ( GDALDatasetH hDS,
int nBandCount,
int * panBandList,
int nLayerCount,
OGRLayerH * pahLayers,
GDALTransformerFunc pfnTransformer,
void * pTransformArg,
double * padfLayerBurnValues,
char ** papszOptions,
GDALProgressFunc pfnProgress,
void * pProgressArg
)
下面解释下该接口的参数：
hDS —— 输出的栅格数据，注意该数据必须以update模式打开
nBandCount —— 栅格数据要更新的波段数
panBandList —— 要更新的波段列表
nLayerCount —— 矢量图层数，对应参数矢量图层数组pahLayers的元素个数
pahLayers —— 矢量图层数组
pfnTransformer —— 将几何图形转换为栅格图像的像素/行列的转换器，如果该参数为空，则在内部创建一个。
pTransformArg ——将几何图形转换为栅格图像的像素/行列的转换器所用到的参数
padfLayerBurnValues—— 指定各个输出波段的输出像素值，因此一定是nBandCount个
papszOptions —— 控制栅格化的一系列选项值，是最复杂的一个参数，包括如下选项：
"ATTRIBUTE":
指定属性字段中的字段值作为栅格值写入栅格文件中，该值将输出到所有输出波段中。假如该值指定了，padfLayerBurnValues参数的值将失效并且padfLayerBurnValues参数可以设置为NULL。
"CHUNKYSIZE":
指定该运行操作的块的高度。该值越大所需的计算时间越小。如果该值没有设置或者设置为0则由GDAL的缓存大小根据公式：缓存所占的字节数/扫描函数的字节数得到。所以该值不会超出缓存的大小。
"ALL_TOUCHED":
设置为TRUE表示所有的像素接触到矢量中线或多边形（全部参与矢量化？），否则只是多边形中心或被brezenhams line algorithm算法（注：brezenhams line algorithm算法为一有名的矢量栅格化算法）选中的部分。默认值为FALSE。该选项暂时搞不明白，知道的同学请留言告诉我。
"BURN_VALUE_FROM":
gdal的函数说明说该参数是用于设置几何图形的Z值，就是高程值。我搞不明白设置该值对栅格化有何影响。
"MERGE_ALG":
设置重写或增加新值到栅格数据中。选择REPLACE为 重写，选择ADD为添加一个新值到已存在的栅格数据中。默认值为REPLACE。
pfnProgress —— gdal的进度函数
pProgressArg —— 传递给进度函数的参数值
返回值:
CE_None 表示处理成功，CE_Failure表示有错误发生。
值得注意的是，gdal的命令行程序gdal_rasterize也可以实现矢量栅格化，但该程序并不是直接调用GDALRasterizeLayers函数的，有兴趣的朋友可以看看gdal_rasterize的实现。计划以后提供一个GDALRasterizeLayers的使用例子程序。
参考文献：
1. gdal_alg.hFile Reference
2. GDAL源码剖析（四）之命令行程序说明二，作者李民录。