This paper proposes a morphostatistical characterization of the galaxy distribution through spatial statistical modelling based on inhomogeneous Gibbs point processes. The galaxy distribution is supposed to exhibit two components. The first one is related to the major geometrical features exhibited by the observed galaxy field, here, its corresponding filamentary pattern. The second one is related to the interactions exhibited by the galaxies. Gibbs point processes are statistical models able to integrate these two aspects in a probability density, controlled by some parameters. Several such models are fitted to real observational data via the ABC shadow algorithm. This algorithm provides simultaneous parameter estimation and posterior-based inference, hence allowing the derivation of the statistical significance of the obtained results.