How can I visualize 3D matrices with imagesc and cat functions ?

So I have this code in 2D which is working and I am visualizing the results using:

hhi = imagesc([1 N],[1 N], cat(3,Im_source/MMAX,Im_target/MMAX, zeros(size(Im_source))));

Now I am trying to move into 3D and I'd like to use the same syntaxis:

hhi = imagesc([1 N N],[1 N N],cat(4,Im_source/MMAX,Im_target/MMAX, zeros(size(Im_source))));

, but I am getting the following error:

Error using image
Indexed CData must be size [MxN], TrueColor CData must be size [MxNx3]

Error in imagesc (line 21)
    hh = image(varargin{:},'CDataMapping','scaled');

Error in prova_3D (line 46)
hhi = imagesc([1 N N],[1 N N],cat(4,Im_source/MMAX,Im_target/MMAX,
zeros(size(Im_source))));

In here, Im_source and Im_target are two 91-109-91 matrices, N is the maximum of the size of these matrices (109) and MMAX is the highest intensity pixel (~ 2e3). Later I'd like to update some parameters so the goal here is visualizing the update of those paramateres in these 3D matrices.

Thanks.