Implementation of the CTBTO/IDC polarization analysis for three-component stations.
The feature extraction as implemented at CTBTO IDC for single three-component stations determines back azimuth (station to origin) and slowness, including the uncertainties for both of these. In the IDC source code and database, the back azimuth is referred to as only azimuth.
The algorithm computes polarization attributes for a three-component station using a modification to the Jurkevics  algorithm. Some of these attributes are then used to determine detection azimuth (seazp = P-type azimuth in degrees), detection slowness and azimuth/slowness uncertainties (inang1 = emergence (incidence) angle and rect = rectilinearity).
A fixed noise window of 9.5 seconds ([-30s;-20.5s] with respect to trigger time) and a signal window of 5.5 seconds ([-4s;1.5s] with respect to trigger time) is used. The signal window is subdivided into intervals of 1.5s length which overlap by 50%.
De-mean data according to mean of noise window.
Apply cosine ramp to noise data and filter the entire data window.
Rotate three components into ZNE space.
Compute 3x3 covariance matrix for each interval.
Extract eigenvalues and compute parameters including rectilinearity.
Choose the result set with the largest rectilinearity.
In addition to the extracted back azimuth and slowness values the rectilinearity
is added as a comment to the resulting pick. The comment ID is
DFX:rectilinearity and the comment is the value in string representation.