silx-kit / pyFAI

Fast Azimuthal Integration in Python
Other
104 stars 94 forks source link

integrate2d with method=("no", "csr","python") generates a transpose matrix #2255

Closed EdgarGF93 closed 6 days ago

EdgarGF93 commented 2 weeks ago
import numpy as np
from pyFAI import load
from pyFAI.gui.jupyter import plot2d
data = np.load('/data/scisoft/edgar/pyfai_files/data_eiger2m.npy')
poni = '/data/scisoft/edgar/pyfai_files/poni_eiger.poni'
ai = load(poni)
res2d = ai.integrate2d(data=data, 
                       npt_rad=1000, 
                       npt_azim=360, 
                      )
res2d.intensity.shape
(360, 1000)

image

res2d_histogram = ai.integrate2d(data=data, 
                       npt_rad=1000, 
                       npt_azim=360, 
                       method=("no", "histogram", "python"),
                      )
res2d_histogram.intensity.shape
(360, 1000)

image

res2d_csr = ai.integrate2d(data=data, 
                       npt_rad=1000, 
                       npt_azim=360, 
                       method=("no", "csr", "python"),
                      )
res2d_csr.intensity.shape
(1000, 360)

image

kif commented 1 week ago

Affects all python-CSR integrators

kif commented 1 week ago

Affects also CSC integrator in python

kif commented 1 week ago

One of the issue is that test_csr checks integrate_legacy from cython withintegrate_ng from python ...