Rational functions play a fundamental role in systems engineering for modelling, identification, and control applications. In this paper we extend the framework by Lindquist and Picci for obtaining such models from the circulant trigonometric moment problems, from the one-dimensional to the multidimensional setting in the sense that the spectrum domain is multidimensional. We consider solutions to weighted entropy functionals, and show that all rational solutions of certain bounded degree can be characterized by these. We also consider identification of spectra based on simultaneous covariance and cepstral matching, and apply this theory for image compression. This provides an approximation procedure for moment problems where the moment integral is over a multidimensional domain, and is also a step towards a realization theory for random fields.