0

我需要从文件中提取栅格（存储为 numpy 数组）。在非常流行的OGR Cookbook之后，我正在阅读 OGR 层（geojson），然后对向量进行光栅化。我使用 GDAL 的ReadAsArray()函数读取了该数组。这一切都很好，我可以对它做各种麻木的事情。但是，GDAL 会自动写出我创建的 GDAL 数据集，因为它会在程序结束后自动取消引用。我不需要/不想输出这个文件，因为它在磁盘上没有用，我只需要内存中的数据。你怎么能防止这种情况发生？

我试过不调用FlushCache()函数，但文件最终仍然得到输出。代码：

...

    # Create the destination data source
    target = gdal.GetDriverByName('GTiff').Create(output_raster_path, source_raster.RasterXSize, source_raster.RasterYSize, 1, gdal.GDT_UInt16)
    target.SetGeoTransform(source_raster.GetGeoTransform())
    target.SetProjection(source_raster.GetProjection())
    band = target.GetRasterBand(1)
    band.SetNoDataValue(no_data_value)
    gdal.RasterizeLayer(target, [1], source_layer, options=["ATTRIBUTE=BuildingID"])
    raster = band.ReadAsArray()
    return raster

之后，一旦程序完成，就会将 geotiff 写入output_raster_path，我刚刚将其设置为“temp.tif”。

4

1 回答 1

1

你可以使用 In-Memory Driver 来做类似的事情。

mem_drv = gdal.GetDriverByName('MEM')
target = mem_drv.Create('', source_raster.RasterXSize, source_raster.RasterYSize, 1, gdal.GDT_UInt16)
于 2019-03-14T18:40:48.050 回答