The angle you want is $$\tan^{-1}\frac{v_{SN}}{v_{WE}}$$
As for the quadrants, many programming languages have a two-argument arctangent function precisely so as to get you into the right quadrant without any effort and also handle the pathological cases of “due north” and “due south”. It is quite likely that Python does too.