3D Cross Correlation

108 Ansichten (letzte 30 Tage)
Wayne
Wayne am 15 Feb. 2021
Bearbeitet: iman taghavi am 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 Kommentar
Quentin Mesnildrey
Quentin Mesnildrey am 10 Jan. 2022
Bearbeitet: Quentin Mesnildrey am 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

Melden Sie sich an, um zu kommentieren.

Akzeptierte Antwort

Pratyush Roy
Pratyush Roy am 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 Kommentar
iman taghavi
iman taghavi am 4 Jul. 2022
Bearbeitet: iman taghavi am 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

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Community Treasure Hunt

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

Start Hunting!

Translated by