3D Cross Correlation

84 次查看（过去 30 天）
Wayne
Wayne 2021-2-15
编辑： iman taghavi 2022-7-4
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 个评论
Quentin Mesnildrey
Quentin Mesnildrey 2022-1-10
编辑：Quentin Mesnildrey 2022-1-10
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

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Pratyush Roy
Pratyush Roy 2021-2-19
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 个评论
iman taghavi
iman taghavi 2022-7-4
编辑：iman taghavi 2022-7-4
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

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Correlation and Convolution 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by