Yes, it is compact: Consider a finite-dimensional complex Hilbert space and the space of all (bounded) linear operators on $H$, denoted by $\mathcal B(H)$. On this space, we can define an inner product  $\langle A,B\rangle:=\mathrm{Tr}A^\dagger B$ for $A,B\in \mathcal B(H)$.  The (sub-)set of density matrices is 
$$\mathcal S(H):=\left\{\rho\in \mathcal B(H)\,|\, \mathrm{Tr}\rho=1, \, \rho \geq 0\right\} \tag 1\quad .$$

Since $\mathcal S(H)$ is a subset of a finite-dimensional Hilbert space, we have to show that it is closed and bounded to prove compactness. 

To show that $\mathcal S(H)$ is closed, consider a convergent sequence of density operators $(\rho_n)_{n\in \mathbb N}$ with limit $\rho \in \mathcal B(H)$. We have to show that $\rho$ is a density operator. To this end, recall that inner products are continuous and thus $$0\leq \mathrm{Tr}\rho_nA \to \mathrm{Tr}\rho A\quad ,\tag 2$$

for all positive semi-definite $A\in \mathcal B(H)$, which implies that $\rho \geq 0$. Choosing $A=\mathbb I_H$, the identity operator on $H$, also shows that $\mathrm{Tr}\rho=1$, so $\rho\in\mathcal S(H)$.

Boundedness of $\mathcal S(H)$ in turn follows from the fact that for all $\rho \in S(H)$ it holds that

$$\mathrm{Tr}\rho^2 \leq 1 \tag 3\quad .$$