Predicting the ground state spin of transition metal complexes is a challenging task. Previous attempts have been focused on specific regions of chemical space, whereas a more general automated approach is required to process crystallographic structures for high-throughput quantum chemistry computations. In this work, we developed a method to predict ground state spins of transition metal complexes. We started by constructing a dataset which contains 2,032 first row transition metal complexes taken from experimental crystal structures and their computed ground state spins. This dataset showed large chemical diversity in terms of metals, metal oxidation states, coordination geometries, and ligands. Then, we analyzed the trends between structural and electronic features of the complexes and their ground state spins, and put forward an empirical spin state assignment model. We also used simple descriptors to build a statistical model with 97% predictive accuracy across the board. With this, we provide a practical and automated way to determine the ground state spin of transition metal complex from structure, enabling the high-throughput exploration of crystallographic repositories.