Factorization of singular integral operators with a Carleman shift via factorization of matrix functions