Writing Datasets


  • appending to existing data

  • context manager

  • write 3d vs write 2d

  • document issues with writing compressed files (per #77)

  • discuss and refer to topics
    • creation options

    • transforms

    • dtypes

    • block windows

Opening a file in writing mode is a little more complicated than opening a text file in Python. The dimensions of the raster dataset, the data types, and the specific format must be specified.

Here’s an example of basic rasterio functionality. An array is written to a new single band TIFF.

# Register GDAL format drivers and configuration options with a
# context manager.
with rasterio.Env():

    # Write an array as a raster band to a new 8-bit file. For
    # the new file's profile, we start with the profile of the source
    profile = src.profile

    # And then change the band count to 1, set the
    # dtype to uint8, and specify LZW compression.

    with rasterio.open('example.tif', 'w', **profile) as dst:
        dst.write(array.astype(rasterio.uint8), 1)

# At the end of the ``with rasterio.Env()`` block, context
# manager exits and all drivers are de-registered.

Writing data mostly works as with a Python file. There are a few format- specific differences.

Supported Drivers

GTiff is the only driver that supports writing directly to disk. GeoTiffs use the RasterUpdater and leverage the full capabilities of the GDALCreate function. We highly recommend using GeoTiff driver for writing as it is the best-tested and best-supported format.

Some other formats that are writable by GDAL can also be written by Rasterio. These use an IndirectRasterUpdater which does not create directly but uses a temporary in-memory dataset and GDALCreateCopy to produce the final output.

Some formats are known to produce invalid results using the IndirectRasterUpdater. These formats will raise a RasterioIOError if you attempt to write to the. Currently this applies to the netCDF driver but please let us know if you experience problems writing other formats.