Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

annotations on waveform fits plots #221

Closed
carltape opened this issue Aug 23, 2023 · 6 comments
Closed

annotations on waveform fits plots #221

carltape opened this issue Aug 23, 2023 · 6 comments

Comments

@carltape
Copy link
Member

Please consider some additional annotations for the MTUQ waveform fits plots. This excerpt from Figure S2 of Alvizuri and Tape (2018), generated with the older code, includes some items that we found useful.
Screen Shot 2023-08-22 at 4 01 21 PM
It shows the length of the P window and surface wave windows (1s -- not used here; 500 s), the number of stations used (N = 34), the number of P window measurements (Np = 0), and number of surface wave measurements (Ns = 102). Ideally it would list the number of first-motion polarities used in the header, but this information is listed as the 4th text row on the left. So "1 (0.74)" represents the measured polarity (1 = up) and the theoretical amplitude factor for this particular mechanisms (0.74). A theoretical value near 0 would represent a near-nodal station. Thanks.

@rmodrak
Copy link
Member

rmodrak commented Aug 23, 2023

Here's one idea-- I think we could make the following two functions "swappable", so to speak.

add_trace_labels()
add_station_labels()

Then users could customize annotations however they like by passing these functions as arguments to the higher level waveform plotting functions. Perhaps this could be a workable solution?

@carltape
Copy link
Member Author

That flexibility would be nice, too. Does it address the overall header line -- the one in the figure with "Pwin 1 Swin 500 N 34 Np 0 Ns 102"? The one item that should not be optional is the length of the time windows, since (without a scale bar) the viewer cannot tell how much time is displayed in the waveforms. (They can guess from the bandpass listed.) On a related note, given the recent tests with source duration, it might be good to list some text related to what kind of stf (or hdur) is used. And: One cannot tell the difference between a constrained inversion (e.g., deviatoric or double couple) and an unconstrained one. (One could use gamma delta 7 0* for a deviatoric inversion, to distinguish it from gamma delta = 7 0, which would indicate a FMT run that happens to be deviatoric.) Sorry, lots of different thoughts here. The overarching item is that by looking at the waveform fits plot, we should be able to tell most of the key choices for the inversion run.

@rmodrak
Copy link
Member

rmodrak commented Aug 23, 2023

  • I believe the length of the time windows is now included in the third line of the MTUQ header. Implementing scale bars is also welcome, if anyone wants to.

  • Perhaps the STF information could be added to the second line of the MTUQ header? A challenge is that sometimes the line gets cut off (e.g. if the Earth model has a long name)

  • Agreed, an asterisk to denote a fixed parameter is a great idea. We would just have to think about how to implement this, since the current plotting function arguments may not provide the necessary information.

  • For me, certain closely spaced text can be difficult to read (an issue since childhood), so I like the idea of allowing customization-- allowing both easy-to-read headers and more detailed headers, depending on preference.

@carltape
Copy link
Member Author

Ah, you're right about the window lengths! I see them here.
These are "when time permits" items, but good to have the bullets. Thanks.

@rmodrak
Copy link
Member

rmodrak commented Aug 27, 2023

Taking a stab at customizable annotations here:
7aa099a

@rmodrak
Copy link
Member

rmodrak commented Aug 29, 2023

Customizable trace and station annotations added in #224

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants