Mastering Geospatial Analysis with Python

Reading and writing raster data with Rasterio

After covering how to read and write various vector data formats in Python, we'll now do the same for raster data. We'll start with the Rasterio library and have a look at how we can read and write raster data. Open up a new Jupyter Notebook where you have access to the Rasterio library and type the following code:

In: import rasterio        dataset = rasterio.open(r"C:\data\gdal\NE\50m_raster\NE1_50M_SR_W    \NE1_50M_SR_W.tif")

This imports the rasterio library and opens a GeoTIFF file. We can now perform some simple data description commands, such as printing the number of image bands.

Raster images contain either a single or multiple bands. All bands are contained in a single file, with each band covering ...