recon_shape = (tomo1.shape[0], tomo1.shape[2], tomo1.shape[2])