combined time series codes that I had in different projects:
-
crawtools
decimate
spectral
-
rptransient
: All except:decimate
(used one fromcrawtools
)read_mseed
(changed to independent function
-
wayne_obstools
fir_corr
=>timeseries
?Peterson_model
,plot_PPSD
,plot_response
,plot_sensitivity
Made sure that all of the equations correspond to equations in Bendat&Piersol (2010). Verified that we now get as good of results as the old Matlab code. Set up testing for the different sub-modules. Probably changed some method call parameters. Created a readthedocs page.
SpectralDensity
:- Added required argument
window_s
to creator
- Added required argument
- Renamed
TransferFunction
properties:values
tofrf
(frequency response function)uncert
touncertainty
channels
tochannel_names
TimeSpans
:- Now created using
spans
by default
- Now created using
DataCleaner
- always renames cleaned channels, putting information in the location code slot (second from last element in string separated by '.'s)
- Removed unused
fast_calc
parameter fromclean_sdf()
- Changed
SpectralDensity.from_stream()
z_threshold to apply to log10(spectra) rather than (spectra) (otherwise, "standard" zthreshold=3 rejects more than half of typical data) TransferFunction()
- now provides a "frequency response function" (
frf
) and acorrector
(for data cleaning): the two are the same only if all noise is on the output channel (see Bendat and Piersol chapter 6) - fixed (and simplified) calculation of uncertainty
- now provides a "frequency response function" (
DataCleaner
usesTransferFunction
'scorrector
property.
SpectralDensity
:- Added arguments
ts_starttime
andts_endtime
to creator - Added ``used()
,
unused()` methods and `window_seconds` property autospect()
and.crossspect()
now accept wildcards in channel names- Added
channel_name()
method - Added
**fig_kw
toplot_*()
methods
- Added arguments
TimeSpans
:- Added
spans
property - Added
invert()
method plot()
can include a trace or stream
- Added
- Updated documentation
DataCleaner
:- Created
CleanerString
class to help with channel names
- Created
Fix `TransferFunction` errorbar plotting bug
Change `MANIFEST.in` to recursively include subfiles/directories of decimate/
- Decimator.decimate() now returns same data.dtype by default
- TransferFunction now accepts wildcards for out_chan names
- Added 'outfile' to TransferFunction.plot()
- Added TransferFunction.put_response()
-
Renamed TransferFunction to ResponseFunction, changed internal and method names to correspond:
- .response -> instrument_response
- .frf -> .value
- .frf_wrt_counts -> .value_wrt_counts
- SpectralDensity.channel_response -> SpectralDensity.channel_instrument_response
- SpectralDensity.put_channel_response -> SpectralDensity.put_channel_instrument_response
- SpectralDensity.init(response=) -> SpectralDensity.init(instrument_response)
-
Internally renamed
tiskit
totiskitpy
Rewrote tracking of cleaning steps, a lot of internal work, including
new classes, but should be invisible when using the command-line codes
and mostly invisible when using the API.
- Added channel identification by
tiskitpy_id
, which includes cleaning information. Thetiskitpy_id
for uncleaned data is theseed_id
. - Added class CleanedStream and revised the guts of several classes, including
renaming
SpectralDensity.channel_names
toSpectralDensity.ids
. - CleanRotator class now has a property
variance_reduction
which gives the variance reduction obtained during init() - Added
SpectralDensity.plots()
and.plots_coherences()
, to compare multipleSpectralDensity
objects