________________________________ From: Alex Highviz <alexhigh...@hotmail.com> Sent: 22 June 2018 18:27:14 To: even.roua...@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
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev
