Laguerre-Gaussian mode expansion for arbitrary optical fields using a subspace projection method