The numerical values of the components of the structure constants can be not only "written down". They may be explicitly calculated from the explicit form of the generators. The values that you referred to correspond to a particular basis that is the most widespread convention, the Gell-Mann matrices:
They're the basis of all 3x3 Hermitian traceless matrices. Two of them, the 3rd and 8th, are taken to be diagonal. The remaining 6 are off-diagonal matrices that only have $+1,+1$; or $+i,-i$ at the $ij$- and $ji$-location of the matrix. The order is so that we do $ij=12,13,23$ in this order, and for each value of $ij$, we first write the matrix with $1,1$ and then with $+i,-i$, with the $-i$ above the diagonal.
The first three Gell-Mann matrices are just the three Pauli matrices for $SU(2)$. Then we're adding the remaining five matrices – out of them, the last one is the diagonal one, up to the overall sign, it's the unique diagonal matrix that is traceless and orthogonal to all previous matrices (including and especially the diagonal, third Gell-Mann or Pauli matrix).
All of the Gell-Mann matrices are normalized so that the trace of their square equals two.
If you just substitute $T^a = \lambda^a/2$ for your generators, you may simply calculate all the commutators and identify them as another Gell-Mann matrix.
Any basis of the 8-dimensional space of $SU(3)$ generators will contain either "random numbers" or entries of different kinds. For example, in a comment, marmot mentioned that you may choose the basis from the Cartan algebra and then the eigenstates corresponding to some roots etc. In this basis, the structure constants will still look somewhat random and you will have to distinguish two types of the indices $a,b,c=1,2,\dots 8$. But there will be a nicer logic in them than in the Gell-Mann basis.