This paper discusses the use of robust geostatistical methods on a data set of rainfall measurements in Switzerland. The variables are detrended via non-parametric estimation penalized with a smoothing parameter. The optimal trend is computed with a smoothing parameter based on cross-validation. Then, the variogram is estimated by a highly robust estimator of scale. The parametric variogram model is fitted by generalized least squares, thus taking account of the variance-covariance structure of the variogram estimates. Comparison of kriging with the initial measurements is completed and yields interesting results. All these computations are done with the software S+SpatialStats, extended with new functions in S+ that are made available.