[FEATURE] add --asarray option to gdal_calc.py, which will enable it to handle more than 26 input files and simplify its usage

[xg1990]

When the option --asarray is toggled, value of the input files will be stored in a list named A like this:

[

[ value_of_Dataset1_Band1, value_of_Dataset1_Band2, ..., ]

[ value_of_Dataset2_Band1, value_of_Dataset2_Band2, ..., ]

....

[ value_of_Datasetn_Band1, value_of_Datasetn_Band2, ..., ]

]

Besides, all the trailing arguments will be treated as input files. And the alpha list is not necessary any more.

For example, user can calculate the maximum of the 1st band of all the TIFF files in some folder by:

gdal_calc.py --asarray --calc="numpy.max([a[0] for a in A])" ./some_folder/*.tif

or calculate the average of all the bands by:

gdal_calc.py --asarray  --calc="sum([sum(a) for a in A]) / sum([len(a) for a in A])" ./some_folder/*.tif

Simple operation such as addition is also possible:

gdal_calc.py --asarray --calc="A[0][0] + A[1][0]"  F1.tif F2.tif

Most importantly, many numpy functions can be useful now.

[xg1990]

This is the modified file

[Even Rouault]

It would also be necessary to enhance swig/python/scripts/gdal_calc.dox and autotest/pyscripts/test_gdal_calc.py

[etourigny]

sorry Even I don't have time to look into this, in fact I will be unavailable for a pretty long time.

[xg1990]

Replying to rouault:

It would also be necessary to enhance swig/python/scripts/gdal_calc.dox and autotest/pyscripts/test_gdal_calc.py

Maybe I can help. But I'm not sure what to do with them.

[Jukka Rahkonen]

Related discussion in the gdal-dev list: ​http://thread.gmane.org/gmane.comp.gis.gdal.devel/40781

[Even Rouault]

This ticket has been automatically closed because Trac is no longer used for GDAL bug tracking, since the project has migrated to GitHub. If you believe this ticket is still valid, you may file it to ​https://github.com/OSGeo/gdal/issues if it is not already reported there.