I am working with a large (30GB) 3 band geotiff. I am using the following code to attempt to import the bands into python
import gdal import numpy as np raster_path = "C:/Path/To/Image.tif" raster_dataset = gdal.OpenEx(raster_path, gdal.GA_ReadOnly) geo_transform = raster_dataset.GetGeoTransform() proj = raster_dataset.GetProjectionRef() bands_data =  for b in range(1, raster_dataset.RasterCount+1): band = raster_dataset.GetRasterBand(b) bands_data.append(band.ReadAsArray())
However, when I try to iterate over the bands, all my memory gets used and it won't load for hours.
Any suggestions regarding other ways to handle large spatial data in python? I have used ArcGIS and R in the past, but am relatively new to python. I would prefer to work from command line, if possible.