Applying cross power spectral density in Python to array with invalid values

Question

I want to apply cross power spectral density (CPSD) in Python using the scipy.signal.csd function:

https://docs.scipy.org/doc/scipy/reference/generated/scipy