The coseismic surface three-dimensional deformation fields were significant for quantified the geometric and kinematic characteristics of earthquake rupture faults. Traditional geodetic techniques were constrained by their intrinsic limitations, such as InSAR, which could only extract far-field deformation fields due to its incoherence, and GNSS, which could only acquire deformations at discrete points. The recently developing optical pixel correlation technique based on high resolution remote sensing images could acquire the near-field coseismic horizontal deformation. In this study, the vertical displacement field of the 2022 Mw 6.6 Menyuan earthquake was obtained by differencing the pre- and post-earthquake DEMs from GaoFen-7 stereo satellite images. Further, the line-of-sight (LOS) and azimuthal (AZI) directions far-field deformations by InSAR, the horizontal near-field deformation by optical pixel correlation based on pre- and post-earthquake GF-2/7 images, and the vertical deformation by differencing the pre- and post-earthquake DEMs, were incorporated to comprehensively solve the complete three-dimensional deformation fields of the 2022 Mw 6.6 Menyuan earthquake. Firstly, the vertical displacement field by differencing DEMs indicated there were significant vertical displacements about 1 m at the bend region, which was induced by the local compressive stress. In addition, at the epicenter on the middle segment of ruptured Lenglongling fault (LLLF), the maximum lifting occurred on the southern sides of the main and secondary faults exceeded 2 m. Secondly, the three-dimensional deformation fields solved by multiple deformation data demonstrated that the near-field deformation field calculated by optical pixel correlation method could quantify the displacements distributed over the rupture fault zone, which were not available on the InSAR deformation maps. Finally, the surface 2D strain derived from the displacement maps calculated by optical pixel correlation revealed high strain concentration on the rupture fault zone. Our study focused on the complete surface three-dimensional deformation estimation from multiple far- and near-field deformation data, and provided a new perspective for a deeper understanding of the characteristics of coseismic surface deformation and the rupture pattern of the fault.