Many sulphur-bearing species have been detected in different astronomical environments and have allowed to derive important information about the chemical and physical composition of interstellar regions. In particular, these species have also been showed to trace and probe hot-core environment time evolution. Among the most prominent sulphur-bearing molecules, SO, sulphur monoxide radical, is one of the more ubiquitous and abundant, observed also in its isotopic substituted species such as $^{34}$SO and S$^{18}$O. Due to the importance of this simple diatomic system and to face the challenge of modern radioastronomical facilities, an extension to THz range of the rare isotopologues of sulphur monoxide has been performed. High-resolution rotational molecular spectroscopy has been employed to extend the available dataset of four isotopic species, SO, $^{34}$SO, S$^{17}$O, and S$^{18}$O up to the 1.5 THz region. The frequency coverage and the spectral resolution of our measurements allowed a better constraint of the molecular constants of the four species considered, focusing especially for the two oxygen substituted isotopologues. Our measurements were also employed in an isotopically invariant fit including all available pure rotational and ro-vibrational transitions for all SO isotopologues, thus enabling accurate predictions for rotational transitions at higher frequencies. Comparison with recent works performed on the same system are also provided, showing the quality of our experiment and the improvement of the datasets for all the species here considered. Transition frequencies for this system can now be used with confidence by the astronomical community well into the THz spectral region.