[gdal-dev] Raster algebra with Gdal and Python

[Greg Coats]

Yes, GDAL can successfully process images with 15,000 columns by 15,000 rows. I have used GDAL on a computer with only 3 GB RAM to process an image with Image Width: 260,000 Image Length: 195,000. What version of GAL are you using? Are you using GDAL 1.7.0? What operating system are you using? How much RAM is available to GDAL? Greg

On Jan 31, 2010, at 1:40 PM, Alexander Bruy wrote:

> Second problem: processing large rasters. I get error
> 
> File "C:/Documents and Settings/alex/.qgis//python/plugins\rastercalc\rastercalcutils.py", line 162, in layerAsArray
>    array = gdalData.ReadAsArray()
> File "C:\OSGeo4W\apps\gdal-16\pymod\osgeo\gdal.py", line 603, in ReadAsArray
>    return gdalnumeric.DatasetReadAsArray( self, xoff, yoff, xsize, ysize )
> File "C:\OSGeo4W\apps\gdal-16\pymod\osgeo\gdal_array.py", line 108, in DatasetReadAsArray
>    return numpy.concatenate( array_list )
> MemoryError
> 
> Large rasters is landsats with 6 bands, 30m/pixel, scene size is aprox. 185x185 km, aprox. 9000x9000 pixels. On disk one file occupied ~500-700 Mb. Some files
> can be larger, up to 15000x15000 pix. and more.
> Is it possible to work with large files?