Simulate/Correct seismometer.
Parameters: |
|
---|---|
Returns: | The corrected data are returned as numpy.ndarray float64 array. float64 is chosen to avoid numerical instabilities. |
This function works in the frequency domain, where nfft is the next power of len(data) to avoid wrap around effects during convolution. The inverse of the frequency response of the seismometer (paz_remove) is convolved with the spectrum of the data and with the frequency response of the seismometer to simulate (paz_simulate). A 5% cosine taper is taken before simulation. The data must be detrended (e.g.) zero mean beforehand. If paz_simulate=None only the instrument correction is done. In the latter case, a broadband filter can be applied to the data trace using pre_filt. This restricts the signal to the valid frequency band and thereby avoids artefacts due to amplification of frequencies outside of the instrument’s passband (for a detailed discussion see Of Poles and Zeros, F. Scherbaum, Kluwer Academic Publishers).
Changed in version 0.5.1: The default for remove_sensitivity and simulate_sensitivity has been changed to True. Old deprecated keyword arguments paz, inst_sim, no_inverse_filtering have been removed.