It seems that tf.nn.convolution should be able to do 4D convolutions, but I haven't been able to succeed in creating a Keras layer to use this function.
I have tried using Keras Lambda layer to wrap the tf.nn.convolution function but maybe somebody else has a better idea?
I want to exploit the higher-dimensional structure of my data, so reshaping might not be able to capture the nature of my dataset.