Optimal characterization of Gaussian channels using photon-number-resolving detectors