Looking back on my math education, I noticed that even though the trace of an endomorphism came up a lot, I'd be hard pressed to give a good description of what the trace really means. I'm looking for some good intuition about its meaning.
To elaborate on what I'm looking for: if I forgot the rigorous definition of the determinant, I could rebuild it from scratch without reference to a basis because I know that it's supposed to measure the change in volume and orientation of a parallelepiped under a linear transformation. For the same reason, I can quickly tell that it is independent of basis, multiplicative, and determines wether the endomorphism is injective or not, all without doing any calculations. I want something similar for the trace. It doesn't need to be geometric, but I want to know what the trace tells us about how the endomorphism acts.

One way to think about it is that the trace is the unique linear map with the property
$$\operatorname{tr}(|u\rangle\langle v|) = \langle u | v\rangle $$
In particular in the case of a rank 1 matrix $A=\lambda|u\rangle\langle v|$ we see that the trace in some sense measures "how much $\ker(A)^\perp$ is aligned with $\operatorname{Im}(A)$". Does this viewpoint carry over to matrices that are not rank-1?
Consider a higher rank matrix $A=\sum_{i=1}^k \sigma_i|u_i\rangle\langle v_i|$. Again we may assume $u_i$ and $v_i$ are normalized (you may notice that when $k=n$ and the $u_i$ and $v_j$ are orthogonal this is the SVD of $A$) then
$$\operatorname{tr}(A)= \operatorname{tr}(\sum_i \sigma_i|u_i\rangle\langle v_i|) = \sum \sigma_i\langle v_i|u_i\rangle $$
Which is a weighted sum of how much the orthogonal complements of the kernels of the rank-1 components align with their images. Due to linearity, it does not matter how we represent $A$ as a sum of rank-1 matrices.
Notable special cases:
if $A$ has an orthogonal eigenbasis, then $A=\sum_i \lambda_i |v_i\rangle\langle v_i|$ and so $\operatorname{tr}(A)=\sum_i\lambda_i\langle v_i|v_i\rangle = \sum_i \lambda_i$. Here the orthogonal complements and the images of the rank 1 components are perfectly aligned.
For a projection matrix $P$ (i.e. $P^2=P$) we have $\operatorname{tr}(P)=\dim \operatorname{Im}(P)$. Since $P$ acts like an identity on the subspace is projects onto, again the orthogonal complements of the kernel and the image are perfectly aligned.
For a nilpotent matrix $N$, the trace is zero. In fact, we can write $N$ as a sum of rank-1 matrices where the orthogonal complement of the kernel is orthogonal to the image for each of them. (can be proven for example via Schur decomposition)