I'm attempting to make use of wavelet scattering to analyse my signal samples.
The example IQ Wav file I am using is as follows: https://www.dropbox.com/s/dd6fr4va4alpazj/move_x_movedist_1_speed_25k_sample_6.wav?dl=0
Currently, following a guide from the Kymatio documentation, I have the following code:
import scipy.io.wavfile
import numpy as np
import matplotlib.pyplot as plt
from kymatio.numpy import Scattering1D
path = r"move_x_movedist_1_speed_25k_sample_6.wav"
# Read in the sample WAV file
fs, x = scipy.io.wavfile.read(path)
x = x.T
print(fs)
# Once the recording is in memory, we normalise it to +1/-1
x = x / np.max(np.abs(x))
print(x)
# Set up parameters for the scattering transform
## number of samples, T
N = x.shape[-1]
print(N)
## Averaging scale as power of 2, 2**J, for averaging
## scattering scale of 2**6 = 64 samples
J = 6
## No. of wavelets per octave (resolve frequencies at
## a resolution of 1/16 octaves)
Q = 16
# Create object to compute scattering transform
scattering = Scattering1D(J, N, Q)
# Compute scattering transform of our signal sample
Sx = scattering(x)
# Extract meta information to identify scattering coefficients
meta = scattering.meta()
# Zeroth-order
order0 = np.where(meta['order'] == 0)
# First-order
order1 = np.where(meta['order'] == 1)
# Second-order
order2 = np.where(meta['order'] == 2)
#%%
# Plot original signal
plt.figure(figsize=(8, 2))
plt.plot(x)
plt.title('Original Signal')
plt.show()
# Plot zeroth-order scattering coefficient (average of
# original signal at scale 2**J)
plt.figure(figsize=(8,8))
plt.subplot(3, 1, 1)
plt.plot(Sx[order0][0])
plt.title('Zeroth-Order Scattering')
# Plot first-order scattering coefficient (arrange
# along time and log-frequency)
plt.subplot(3, 1, 2)
plt.imshow(Sx[order1], aspect='auto')
plt.title('First-order scattering')
# Plot second-order scattering coefficient (arranged
# along time but has two log-frequency indicies -- one
# first- and one second-order frequency. Both are mixed
# along the vertical axis)
plt.subplot(3, 1, 3)
plt.imshow(Sx[order2], aspect='auto')
plt.title('Second-order scattering')
plt.show()
The print statements for the sample rate, the normalised X (data) and the number of samples are as follows:
2000000
[[-0.65671316 0.40170009]
[-0.67349608 0.50152572]
[-0.62685266 0.54555362]
...
[-0.59829991 -0.11006975]
[-0.64930253 -0.03116827]
[-0.6619442 0.05557977]]
2
There seems to be 2 samples as the main signal was split up into chunks of 2s using pydub's make_chunks and exporting them as wav format. The sample rate is 2MHz as these recordings are of RF signals stored in IQ WAV files.
However, the problem here is that I am getting the following error, from which I read seems to potentially be an issue with my J value (averaging scale):
Traceback (most recent call last): File "analyse.py", line 38, in scattering = Scattering1D(J, T, Q) File "C:\Users\xwb18152\AppData\Roaming\Python\Python38\site-packages\kymatio\scattering1d\frontend\numpy_frontend.py", line 19, in init ScatteringBase1D.build(self) File "C:\Users\xwb18152\AppData\Roaming\Python\Python38\site-packages\kymatio\scattering1d\frontend\base_frontend.py", line 55, in build min_to_pad = compute_minimum_support_to_pad( File "C:\Users\xwb18152\AppData\Roaming\Python\Python38\site-packages\kymatio\scattering1d\utils.py", line 125, in compute_minimum_support_to_pad _, _, _, t_max_phi = scattering_filter_factory( File "C:\Users\xwb18152\AppData\Roaming\Python\Python38\site-packages\kymatio\scattering1d\filter_bank.py", line 676, in scattering_filter_factory psi_f[0] = morlet_1d( File "C:\Users\xwb18152\AppData\Roaming\Python\Python38\site-packages\kymatio\scattering1d\filter_bank.py", line 135, in morlet_1d morlet_f *= get_normalizing_factor(morlet_f, normalize=normalize) File "C:\Users\xwb18152\AppData\Roaming\Python\Python38\site-packages\kymatio\scattering1d\filter_bank.py", line 157, in get_normalizing_factor raise ValueError('Zero division error is very likely to occur, ' + ValueError: Zero division error is very likely to occur, aborting computations now.
Any help would be appreciated to understand what's going on!
EDIT: I can get to "the next step" if I use Q=1,J=1, however I get the following warning:
warnings.warn('Signal support is too small to avoid border effects'
And when trying to plot I get the error:
File "analyse.py", line 100, in <module>
plt.imshow(Sx[order1], aspect='auto')
[....]
raise TypeError("Invalid shape {} for image data"
TypeError: Invalid shape (2, 3, 1) for image data
