.was ending maybe it is something like that, not sure
I have been experimenting with the PPSD spectral_estimation.py code, adding some colorization to the various plot overlays such as the mean, mode, and percentile graphs. I want to add a little bit of optimization to the PPSD plot, as the overlays default to black. How would I go about becoming a first-time contributor to the project ?
nearest_samplethat sample might be left or right of the given time. If you need to have the same time and the samples dont align you could resample your data to have the exact same sampling points, but personally I would not do that and rather adjust the timing in the cross correlation function that you get as output
from obspy.clients.fdsn import Client client = Client('IRIS') inv = client.get_stations(network='BX', level='station') sta_codes_bx = [sta.code for net in inv for sta in net] for pick in event.picks: is_IR = pick.waveform_id.network_code == 'IR' should_be_BX = pick.waveform_id.station_code in sta_codes_bx if is_IR and should_be_BX: pick.waveform_id.network_code = 'BX' client.get_waveforms(....)
If I use the command scan inside a script.py it raises an error like this:
PermissionError: [WinError 5] Access denied: '/$Recycle.Bin\S-1-5-18'
on Ubuntu the error is almost the same.
What I do is the following:
from obspy.imaging.scripts.scan import scan
Trying to emulate what I do from the shell:
obspy-scan input_trace_path -v -o outputTracePath
From the shell I do not have any problem, from the script it doesn't work.
I need help with instrument response removal using poles and zeros. All the information I have about the seismometer is stated here (https://www.gfz-potsdam.de/en/section/geophysical-imaging/infrastructure/geophysical-instrument-pool-potsdam-gipp/pool-components/poles-and-zeros/3d-geophone) and below:
No Damping Resistor
Gs = 28.8 V/(m/s); Ds = 0.56; fs = 4.5 Hz
z1 = 0
z2 = 0
p1 = -15.88 23.43i
p2 = -15.88 -23.43i
a = 1.0
I can see the code example using PAZ dictionary here: https://docs.obspy.org/tutorial/code_snippets/seismometer_correction_simulation.html#using-a-paz-dictionary but not sure which of these parameters correspond to the ones used in the example.
There is also a response file available here: https://www.gfz-potsdam.de/fileadmin/gfz/sec22/pdf_doc/GIPP/RESP.XX.NS532..SLZ.HG6.375.100000.4_5B.28_8.txt and also on IRIS NRL: http://ds.iris.edu/NRL/sensors/hgs/RESP.XX.NS532..SLZ.HG6.375.9090.4_5B.27_7 but didn't work with ObsPy.
I kept getting "ValueError: No matching response information found" error when I use these files (after modifying the channel name, station name and start and end dates).
Kindly provide me some guidance on using the parameters or response file successfully.
if key == 'NOTE': try: value = [cleanup_and_decode_string(line) for line in value.split(self.line_terminator) if line] except: value = '' else: value = cleanup_and_decode_string(value)