Shift the data from the time domain to the frequency domain, possibly with an FFT:
https://en.wikipedia.org/wiki/Fast_Fourier_transform
Some micros (Renasys?) have instructions that make it easier ton perform.
Another possible way to look at the problem is through a "multi-channe (MCA)l" analyzer. Where you might have an acceleration(x) and a magnitude. Not really sure how to handle that. An MCA is used to detect X-Ray events. Your essentially looking at the number and the energy of the event:
https://en.wikipedia.org/wiki/Energy-dispersive_X-ray_spectroscopy. Analysis software basically strips off the bottom "junk" and then tries to find a peak using the FWHM (Full Width Half Maximum) method.
Or even spectral analysys:
https://www.electro-tech-online.com/custompdfs/2012/04/Spectrum-Analysis.pdf
The thing that might be wierd with your situation is that the the accelerometer data seems n-dimensional, so maybe collect it as a phasor.
Hey, you just asked for "ideas". I'm not sure what might work. In my head, I'm seeing a 3-D surface. # if events, direction, magnatude.
So, I'm seeing binning the levels. e.g discrete levels. (window dector) so to speak.
Dividing the direction into binnable levels 1 degree apart?)
And counting the events,
So you have a "surface". The surface can be displayed in "intensity modulated", "color modulated" or "y-modulated" form. The "Y-modulated" data could also be shown in color.
Hey, you asked for ideas. I did what I could.