We propose a novel estimation procedure for certain spectral distributions associated with a class of high dimensional linear time series. The processes under consideration are of the form $X_t = \sum_{\ell=0}^\infty \mathbf{A}_\ell Z_{t-\ell}$ with iid innovations $(Z_t)$. The key structural assumption is that the coefficient matrices and the variance of the innovations are simultaneously diagonalizable in a common orthonormal basis. We develop a strategy for estimating the joint spectral distribution of the coefficient matrices and the innovation variance by making use of the asymptotic behavior of the eigenvalues of appropriately weighted integrals of the sample periodogram. Throughout we work under the asymptotic regime $p,n \to \infty$, such that $p/n\to c \in (0,\infty)$, where $p$ is the dimension and $n$ is the sample size. Under this setting, we first establish a weak limit for the empirical distribution of eigenvalues of the aforementioned integrated sample periodograms. This result is proved using techniques from random matrix theory, in particular the characterization of weak convergence by means of the Stieltjes transform of relevant distributions. We utilize this result to develop an estimator of the joint spectral distribution of the coefficient matrices, by minimizing an $L^κ$ discrepancy measure, for $κ\geq 1$, between the empirical and limiting Stieltjes transforms of the integrated sample periodograms. This is accomplished by assuming that the joint spectral distribution is a discrete mixture of point masses. We also prove consistency of the estimator corresponding to the $L^2$ discrepancy measure. We illustrate the methodology through simulations and an application to stock price data from the S\&P 500 series.