We study the problem of compressed sensing for quaternionic Fourier matrices as arising in color representation of images. We will show that such matrices are allowing a sparse reconstruction by means of an \(l_1\)-minimization with high probability. Examples of sparse sampling of color images are provided.