I have been working on that with some porphyrins. I found a couple of papers where they simply do:
Spin density = alpha population - beta population.
(DOI: 10.1007/s00894-011-1279-x)
In Gaussian, request an NBO analysis with the keyword 'pop=nbo'. Then in your log file check for 'Summary of Natural Population Analysis:' Which looks something like this: