Объедините два массива numpy в 4-м измерении

У меня есть два массива numpy с тремя измерениями (3 x 4 x 5), и я хочу объединить их, чтобы результат имеет четыре измерения (3 х 4 х 5 х 2). В Matlab это можно сделать с помощью cat (4, a, b) , но не с помощью Numpy.

Например:

a = ones((3,4,5))
b = ones((3,4,5))
c = concatenate((a,b), axis=3) # error!

Чтобы уточнить, я хочу, чтобы c [:,:,:, 0] и c [:,:,:, 1] соответствовали исходные два массива.