Polynomial Chaos Expansion for Operator Learning