如何使用GDAL栅格化（python API）按属性指定刻录值？

Question

我正在使用GeoTiff模板将形状文件转换为GeoTiff，同时通过刻录输出值。我想输出的是一个.tif，其中刻录值对应于给定属性的值。我发现这是烧焦到与我的属性字段中的值不对应的奇怪值。以下是我目前拥有的内容：gdal.RasterizeLayer()ATTRIBUTEgdal.RasterizeLayer()    gdalformat = 'GTiff'    datatype = gdal.GDT_Byte    # Open Shapefile    shapefile = ogr.Open(self.filename)    shapefile_layer = shapefile.GetLayer()    # Get projection info from reference image    image = gdal.Open(ref_image, gdal.GA_ReadOnly)    output = gdal.GetDriverByName(gdalformat).Create(output_tif, image.RasterXSize, image.RasterYSize, 1, datatype,                                                     options=['COMPRESS=DEFLATE'])    output.SetProjection(image.GetProjectionRef())    output.SetGeoTransform(image.GetGeoTransform())    # Write data to band 1    band = output.GetRasterBand(1)    band.SetNoDataValue(0)    gdal.RasterizeLayer(output, [1], shapefile_layer, options=['ATTRIBUTE=FCode'])    # Close datasets    band = None    output = None    image = None    shapefile = None    # Build image overviews    subprocess.call("gdaladdo --config COMPRESS_OVERVIEW DEFLATE " + output_tif + " 2 4 8 16 32 64", shell=True)发生的情况是，输出.tif为每个属性正确分配了不同的刻录值，但该值与属性值不对应。例如，输入属性值 FCode=46006 将变为刻录值 182（不清楚为什么！我尝试添加和删除该选项，并添加和删除的“ 3D”选项。不影响输出刻录值。'COMPRESS=DEFLATE'gdal.RasterizeLayer()您可以在此处查看输入形状文件和属性值：输入 .shp以及具有不正确值的输出，如下所示：输出栅格

Answer 1

我自己修复了这个问题，将类型更改为.gdal.GDT_Int32