In this paper, the velocities of 154 stations obtained in 2001 and 2003 GPS survey campaigns are applied to formulate a continuous velocity field by the least-squares collocation method. The strain rate field obtained by the least-sqares collocation method shows more clear deformation patterns than that of the conventional discrete triangle method. The significant deformation zones obtained re mainly located in three places, to the north of Tangshan, between Tianjing and Shijiazhuang, and to the north of Datong, which agree with the places of the Holocene active deformation zones obtained by geological investigations. The maximum shear strain rate is located at latitude 38.6◦N and longitude 116.8◦E, with a magnitude of 0.13ppm/a. The strain rate field obtained can be used for earthquake prediction research in the North China Basin.