Add GDALDataset::[Get|Compute]InterBandCovarianceMatrix() and corresponding C and Python API

Open rouault opened this issue 3 weeks ago • 1 comments

This generates and reads STATISTICS_COVARIANCES band metadata items as done by ArcGIS:

https://pro.arcgis.com/en/pro-app/3.4/tool-reference/spatial-analyst/how-band-collection-statistics-works.htm

which is consistant with the results of numpy.cov(): https://numpy.org/devdocs//reference/generated/numpy.cov.html

Dec 04 '25 22:12 rouault

coverage: 71.407% (+0.006%) from 71.401% when pulling d390084a1499a46898b22773aa8751151af9f8b9 on rouault:covariance_matrix into 379bec35a297acc2c984e647f62d9e303a3013a0 on OSGeo:master.

Dec 08 '25 01:12 coveralls