sto cercando di reimplementare in pitone una funzione IDL:ridimensionare con media o Rebin una matrice 2D NumPy
http://star.pst.qub.ac.uk/idl/REBIN.html
che ridimensiona da un fattore intero array 2d facendo la media.
Ad esempio:
>>> a=np.arange(24).reshape((4,6))
>>> a
array([[ 0, 1, 2, 3, 4, 5],
[ 6, 7, 8, 9, 10, 11],
[12, 13, 14, 15, 16, 17],
[18, 19, 20, 21, 22, 23]])
desidero ridimensionare a (2,3) prendendo la media dei campioni pertinenti, i risultati attesi sarebbe:
>>> b = rebin(a, (2, 3))
>>> b
array([[ 3.5, 5.5, 7.5],
[ 15.5, 17.5, 19.5]])
cioè b[0,0] = np.mean(a[:2,:2]), b[0,1] = np.mean(a[:2,2:4])
e così via.
Credo che dovrei rimodellare una matrice a 4 dimensioni e quindi prendere la media sulla sezione corretta, ma non riuscire a capire l'algoritmo. Vuoi qualche suggerimento?
appena trovato ora questo è un duplicato di http://stackoverflow.com/questions/4624112/grouping-2d-numpy-array-in-average, tuttavia, non sono riuscito a trovare questo prima di utilizzare la funzione di ricerca in StackOverflow. –