This paper proposes a family of weighted- \(\theta \) high-order ADI schemes for sine-Gordon equations in high dimensions in combination with the discretization of a three-level linearized scheme in temporal direction and a classical compact difference scheme in spatial direction. The unique solvability, convergence and stability of the difference scheme in two dimensions are strictly proved with the convergence order four in space and order two in time under the \(L^\infty \) -norm. More remarkably, the present numerical scheme can be extended to cope with the three-dimensional case feasibly. Several benchmark problems in two and three dimensions verify theoretical results and the capacity in simulating the dynamic behaviour of ring solitons. Specially, several novel numerical isosurfaces in three dimensions are presented for the first time.
