Least-squares collocation is a powerful method for combining interpolation, filtering and parameter determination in one single computational step. We show that the method is applicable to the computation of polar motion values from a very large set of basic observational data. In this study, we use the ILS observations from 1900 to 1978.