appending to existing data
write 3d vs write 2d
document issues with writing compressed files (per #77)
- discuss and refer to topics
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. profile.update( dtype=rasterio.uint8, count=1, compress='lzw') 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.
GTiff is the only driver that supports writing directly to disk.
GeoTiffs use the
RasterUpdater and leverage the full capabilities
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
to produce the final output.
Some formats are known to produce invalid results using the
IndirectRasterUpdater. These formats will raise a
if you attempt to write to the. Currently this applies to the
driver but please let us know if you experience problems writing other formats.