Projection methods for high numerical aperture phase retrieval