The objective of this article is to propose a comprehensive solution for analyzing multidimensional non-life claim count data that exhibits time and cross-dependence, as well as zero inflation. To achieve this, we introduce a multivariate INAR(1) model, with the innovation term characterized by either a multivariate zero-inflated Poisson distribution or a multivariate zero-inflated hurdle Poisson distribution. Additionally, our modeling framework accounts for the impact of individual and coverage-specific covariates on the mean parameters of each model, thereby facilitating the computation of customized insurance premiums based on varying risk profiles. To estimate the model parameters, we employ a novel expectation-maximization (EM) algorithm. Our model demonstrates satisfactory performance in the analysis of European motor third-party liability claim count data.