3D multiple linear regression along axis

Question

I would like to calculate a multiple linear regression of the form Y = aX1 + bX2 where Y.shape == X1.shape == X2.shape = (200,192,288). I would like to preform the mult