The covariance matrix function is characterized in this paper for a Gaussian or elliptically contoured vector random field that is stationary, isotropic, and mean square continuous on the compact two-point homogeneous space. Necessary and sufficient conditions are derived for a symmetric and continuous matrix function to be an isotropic covariance matrix function on all compact two-point homogeneous spaces. It is also shown that, for a symmetric and continuous matrix function with compact support, if it makes an isotropic covariance matrix function in the Euclidean space, then it makes an isotropic covariance matrix function on the sphere or the real projective space.