[gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

[Alex HighViz]

________________________________
From: Alex Highviz <alexhighviz at hotmail.com>
Sent: 22 June 2018 18:27:14
To: even.rouault at spatialys.com
Subject: Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark



On 14/06/2018 21:22, Even Rouault wrote:
>
> [ ]
> - 1000x1000 is somewhat small. Perhaps benchmark on larger rasters
>
> [ ]
> - it would be interested if you could also benchmark a VRT with Python numpy-
> based pixel function, and enable Numba. See
> http://www.gdal.org/gdal_vrttut.html#gdal_vrttut_derived_python (or just pure
> numpy based computation + Numba)
>
> [ ]
> Even
>
Hi Even,

I tried a few of your suggestions and am pleasantly surprised. These
results are for 10000  x 10000 matrices

VRT with C++ function : 10.1 seconds

VRT with NUMPY function (no Numba): 9.2 seconds

GDAL Direct using GetLockedBlockRef: 8.7 seconds

GDAL Direct using RasterIO: 30.7 seconds

Pronto Raster: 10.8 seconds

Just FYI
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20180622/7371aaae/attachment.html>