<p>The Kumamoto earthquake is analyzed mainly with InSAR data combined with strong earthquake and GPS data, using a variety of joint InSAR technology methods and multisource data solution methods and comprehensively considering the normalization and weighting of multisource data. The three-dimensional (3D) deformation field is determined. The results show that the joint solution of multisource data can improve the accuracy of the 3D solution deformation results to a certain extent. From the 3D solution results, the maximum east-west deformation caused by the 2016 Kumamoto earthquake is approximately 2 m; the north-south direction mainly manifests expansion and stretching; the northwestern side subsides vertically, with a maximum subsidence of 2 m; and the southeastern side is uplifted. The horizontal deformation characteristics show that the earthquake is dominated by right-lateral strike-slip; the strike is NE-SW, the dip of the seismogenic fault is nearly vertical, and the Futagawa fault has a few normal fault properties. By analyzing the coseismic 3D deformation field, the seismogenic fault can be better understood, which provides a foundation for studying seismic mechanisms.</p>