An integrative approach based on probabilistic modelling and statistical inference for morpho-statistical characterization of astronomical data