After a decade of design and construction, South Africa's SKA-MID precursor MeerKAT has begun its science operations. To make full use of the widefield capability of the array, it is imperative that we have an accurate model of the primary beam of its antennas. We have used an available L-band full-polarization astro-holographic observation and electromagnetic simulation to create sparse representations of the beam using principal components and Zernike polynomials. The spectral behaviour of the spatial coefficients has been modelled using discrete cosine transform. We have provided the Zernike-based model over a diameter of 10 degrees in an associated software tool that can be useful for direction dependent calibration and imaging. The model is more accurate for the diagonal elements of the beam Jones matrix and at lower frequencies. As we get more accurate beam measurements and simulations in the future, especially for the cross-polarization patterns, our pipeline can be used to create more accurate sparse representations of MeerKAT beam.

Modern interferometric imaging relies on advanced calibration that incorporates direction-dependent effects. Their increasing number of antennas (e.g. in LOFAR, VLA, MeerKAT/SKA) and sensitivity are often tempered with the accuracy of their calibration. Beam accuracy drives particularly the capability for high dynamic range imaging (HDR - contrast > 1:$10^6$). The Radio Interferometric Measurement Equation (RIME) proposes a refined calibration framework for wide field of views (i.e. beyond the primary lobe and first null) using beam models. We have used holography data taken on 12 antennas of the Very Large Array (VLA) with two different approaches: a `data-driven' representation derived from Principal Component Analysis (PCA) and a projection on the Zernike polynomials. We determined sparse representations of the beam to encode its spatial and spectral variations. For each approach, we compressed the spatial and spectral distribution of coefficients using low-rank approximations. The spectral behaviour was encoded with a Discrete Cosine Transform (DCT). We compared our modelling to that of the Cassbeam software which provides a parametric model of the antenna and its radiated field. We present comparisons of the beam reconstruction fidelity vs. `compressibility'. We found that the PCA method provides the most accurate model. In the case of VLA antennas, we discuss the frequency ripple over L-band which is associated with a standing wave between antenna reflectors. The results are a series of coefficients that can easily be used `on-the-fly' in calibration pipelines to generate accurate beams at low computing costs.