Rasterio

Rasterio is a Python library for reading and writing geospatial raster data using a clean, idiomatic API based on Numpy arrays.

Rasterio bridges the gap between the powerful GDAL C++ library and the Python ecosystem. It allows developers to handle GeoTIFFs, satellite imagery, and terrain models as standard N-dimensional arrays (ndarrays) while maintaining critical spatial metadata like Coordinate Reference Systems (CRS) and affine transforms. The library supports Python 3.12+ and GDAL 3.8+, offering a robust command line interface (rio) for quick inspections and a Pythonic context manager for safe file handling. Whether you are calculating NDVI from Landsat 8 bands or reprojecting datasets to EPSG:4326, Rasterio provides the efficiency of C with the readability of modern Python.

https://rasterio.readthedocs.io/