Skip to content

Including S preevent noise in the denominator during wiener deconvolution #22

Open
@jaredbryan881

Description

@jaredbryan881

I'm just curious why the denominator in the Wiener deconvolution doesn't include preevent noise spectra from the S components as in Audet (2010).

RfPy/rfpy/rfdata.py

Lines 737 to 744 in 5b4b26b

# Auto and cross spectra
Spp = np.sum(np.real(Fp*np.conjugate(Fp)), axis=0)
Sd1p = np.sum(Fd1*np.conjugate(Fp), axis=0)
Sd2p = np.sum(Fd2*np.conjugate(Fp), axis=0)
Snn = np.sum(np.real(Fn*np.conjugate(Fn)), axis=0)
# Denominator
Sdenom = Spp + Snn

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions