Сопоставление 3D-шаблонов с использованием кросс-корреляции в Matlab

Вот минимальный пример того, что я пытаюсь сделать:

Создать 3D-матрицу

a(:,:,1)=[
    1 2 3 4 1;
    2 3 7 1 4;
    3 7 6 0 9;
    0 3 2 8 1;
    1 4 3 1 1]

a(:,:,2)=[
    1 7 3 4 2;
    2 9 2 3 1;
    1 4 7 7 0;
    1 2 3 4 1;
    0 9 3 3 9]

a(:,:,3)=[
    9 4 0 3 5;
    1 2 3 4 1;
    2 0 2 3 1;
    1 4 2 1 1;
    2 5 7 8 1]

a(:,:,4)=[
    2 3 5 2 0
    0 0 0 0 8
    5 2 7 9 8
    2 4 1 1 0
    6 3 8 7 9]

a(:,:,5)=[
    3 5 1 4 6;
    3 2 8 0 0;
    0 2 1 0 4;
    5 4 5 5 6;
    9 5 9 9 5]

Создать 3D-шаблон

b(:,:,1)=[
    9 4 0;
    1 2 3;
    2 0 2]

b(:,:,2)=[
    2 3 5;
    0 0 0;
    5 2 7]

b(:,:,3)=[
    3 5 1;
    3 2 8;
    0 2 1]

Вычислить взаимную корреляцию (3D-взаимную корреляцию в Matlab). Я думаю, что взаимная корреляция - это то же самое, что свертка с перевернутым шаблоном. Это правильно?

c=convn(a,b(end:-1:1,end:-1:1,end:-1:1));

Найти индексы наилучшего соответствия

[x y z] = ind2sub(size(c),find(c==max(c(:))));
x=x-(size(b, 1) - 1)/2
y=y-(size(b, 2) - 1)/2
z=z-(size(b, 3) - 1)/2

Я читал, что нужно вычесть половину размера шаблона из конечных координат, но у меня больше нет ссылки на страницу с этой информацией. Однако я думаю, что если бы кто-то этого не сделал, возвращаемые координаты были бы не в том месте, где был бы центр шаблона, а в углу шаблона.

В результате моего примера я ожидаю: x=2, y=2, z=4. Matlab говорит мне, что это x=4, y=4, z=4. Однако при изменении шаблона на

b(:,:,1)=[
    9 2 3;  
    4 7 7;
    2 3 4]

b(:,:,2)=[
    2 3 4;
    0 2 3;
    4 2 1]

b(:,:,3)=[
    0 0 0;
    2 7 9;
    4 1 1]

Я получаю правильный результат (x=3, y=3, z=3).

Что мне нужно изменить, чтобы всегда получать правильный результат?