Inclined threading dislocations (TDs) piercing the oriented free surface of a crystal are currently observed after growth of oriented thin films on substrates. Up to date the unique way to treat their anisotropic elastic properties nearby the free surface region is to use the integral formalism, which assumes no dislocation core size and needs numerical double integrations. In a first stage of the work, a new and alternative approach to the integral formalism is developed using double Fourier series and the concept of a finite core size, which is often observed in high-resolution transmission electron microscopy. In a second stage, the integral formalism and the Fourier series approaches are applied to the important case of a TD piercing the basal free surface of a hexagonal crystal. For this particular geometry, easy-to-use expressions are derived and compared to a third approach previously known for a plate-like crystal. Finally, the numerical interest and the convergence of these approaches are tested using the basal free surface of the GaN compound, in particular for TDs with Burgers vectors c and (a + c).