GDAL RasterizeLayer不会将所有多边形刻录到Raster吗？
我正在尝试使用GDAL的RasterizeLayer将shapefile刻录到栅格。给定特定的像素大小，我从另一个shapefile中预先创建了一个感兴趣的区域栅格。然后，此AOI将用作所有后续栅格化（相同数量的列和行，相同的投影和地理转换）的基础。 但是，当我基于相同的像素大小和投影将形状刻录到自己的栅格时，会发生问题。下面的链接（没有足够的代表来发布图像），以棕褐色显示原始shapefile，以及RasterizeLayer刻录数据的深粉红色。浅粉红色是深粉红色栅格数据的nodata值。灰色是完成形状文件刻录的AOI。 给定shapefile多边形的范围，我希望看到底部两个角处的刻录值，以及显示的数据下方的两个像素。但是，显然不是这种情况。 以下是我用来生成这些代码的代码。所有形状都是使用QGIS创建的，并且都在同一投影中创建。（应该注意的是，所示图片中的网格只是用来说明我所使用的像素大小。） from osgeo import ogr from osgeo import gdal aoi_uri = 'AOI_Raster.tif' aoi_raster = gdal.Open(aoi_uri) def new_raster_from_base(base, outputURI, format, nodata, datatype): cols = base.RasterXSize rows = base.RasterYSize projection = base.GetProjection() geotransform = base.GetGeoTransform() bands = base.RasterCount driver = gdal.GetDriverByName(format) new_raster = driver.Create(str(outputURI), cols, rows, bands, datatype) new_raster.SetProjection(projection) …