3D Cross Correlation

107 views (last 30 days)
Wayne
Wayne on 15 Feb 2021
Edited: iman taghavi on 4 Jul 2022
Hi there,
I have two 3D datasets with the first two (X,Y) indexes representing spatial positions, and the last index (Z) representing time. Since the two datasets were collected separately and there could be some positional shifts/errors between the two measurements, I want cross-correlate the two matrices such that I can shift one of the dataset accordingly in all three indexes to match the other one.
I found this https://uk.mathworks.com/matlabcentral/fileexchange/61468-xcorr3 which does 3D xcorr, but it does the correlation between a 3D matrix and a 1D signal. In my case, I would like to cross-correlate between two 3D matrices.
Does anyone have any suggestions? Thank you.
  1 Comment
Quentin Mesnildrey
Quentin Mesnildrey on 10 Jan 2022
Edited: Quentin Mesnildrey on 10 Jan 2022
Dear Wayne,
I hope this message reaches you well. I am facing the same situation and would be very interested in getting your feedback.
My objective is to compare 2 trajectories (defined by 2D coordinates (x,y) over time), which seems pretty close to what you mentionned.
Best regards,
Quentin

Sign in to comment.

Accepted Answer

Pratyush Roy
Pratyush Roy on 19 Feb 2021
Hi Wayne,
Since cross-correlation between two arrays can be assumed to be convolution between the first array and the flipped version of the second array, we can perform N-dimensional convolution between the first array and the flipped second array as a workaround. The convn function might be used to obtain the 3-D convolution. The following code snippet demonstrates the use of convn :
Z = convn(X,Y(end:-1:1,end:-1:1,end:-1:1)) %Here the indices vary from end to 1 as a result of flipping
Hope this helps!
Regards,
Pratyush.
  1 Comment
iman taghavi
iman taghavi on 4 Jul 2022
Edited: iman taghavi on 4 Jul 2022
Just a minor comment! If the signals are complex then the cross correlation would be:
Z = convn(X,conj(Y(end:-1:1,end:-1:1,end:-1:1))); %Here the indices vary from end to 1 as a result of flipping

Sign in to comment.

More Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!