With the enrichment of GNSS data and the improvement in data processing accuracy, GNSS technology has been widely applied in fields such as crustal deformation. The Crustal Movement Observation Network of China (CMONOC) has provided decades of Global Navigation Satellite System (GNSS) data and related data products for crustal deformation research on the Chinese mainland. The coordinate time series of continuously observed reference stations contain abundant information on crustal movements. In contrast, the coordinate time series of periodically observed campaign stations have limited data, making it difficult to separate or remove instantaneous non-tectonic movements from the time series, as performed with reference stations, to obtain a stable and reliable crustal movement velocity field. To address this issue, this paper proposes a method to improve the estimation of crustal movement velocity at campaign stations using the information of neighboring reference stations. This method constructs a Delaunay triangulation of reference stations and fits the periodic movement of each campaign station using an inverse distance weighted interpolation algorithm based on the reference station information. The crustal movement velocity of the campaign stations is then estimated after removing the periodic movement. This method was verified by its application to the estimation of the vertical motion rate at some reference and campaign stations in Yunnan Province. The results show that the accuracy of vertical motion rate estimation for virtual and real campaign stations improved by an average of 24.4% and 9.6%, respectively, demonstrating the effectiveness of the improved method, which can be applied to estimate crustal movement velocity at campaign stations in other areas.