petermakus / seismic Goto Github PK

Problem

Currently, CorrDB does not differ between data from various locations.

Proposed solution

The probably easiest way to implement this into the current framework would be to create an additional subgroup in the hdf5 files, in the sense that the current hdf5 tree would be updated to net0-net1/stat0-stat1/loc0-loc1/cha0-cha1.

Important is to allow for backward compatibility with the old hdf5 structure. probably only possible by allowing for both file structures in the params.yaml.

When attempting to load any correlation from julday 360 (which year doesn't matter), the following error occurs:
ValueError: unconverted data remains: 1

This is caused by an issue with obspy's UTCDateTime and should be fixed with their next release (1.23). However, until then it should be possible to implement a dirty fix.

Hello Peter, I have a question about the reference trace for the dv/v calculations. The documentation on the reference trace says: We can opt for a single reference trace for the whole period or multiple reference traces. Is it also possible to have a single reference trace that doesn't span over the whole time window that the dv/v is calculated over? Thank you!

When remove_response=True, Correlator.pxcorr sometimes acts up, when processing using mpi. This seems to be caused by a process that opens an hdf5 file at CorrelationDataBase._write() and freezes while doing so. Then, the next process wanting to write to that file crashes as files cannot be opened by several processes.

Avoid duplicates

• I searched existing issues

Bug Summary

In the examples, if I run the tutorial "Plot a time-dependent spectrogram"

from seismic.plot.plot_spectrum import plot_spct_series
from matplotlib import pyplot as plt

f, t, S = sc.compute_spectrogram('X9', 'IR1', 'HHE', starttime, endtime, 7200, freq_max=12.5)
fig = plt.figure(figsize=(9, 7))
plot_spct_series(S, f, t, log_scale=True)

the module sc is not found

NameError Traceback (most recent call last)
Cell In[1], line 4
1 from seismic.plot.plot_spectrum import plot_spct_series
2 from matplotlib import pyplot as plt
----> 4 f, t, S = sc.compute_spectrogram('X9', 'IR1', 'HHE', starttime, endtime, 7200, freq_max=12.5)
5 fig = plt.figure(figsize=(9, 7))
6 plot_spct_series(S, f, t, log_scale=True)

NameError: name 'sc' is not defined

Code to Reproduce

from seismic.plot.plot_spectrum import plot_spct_series
from matplotlib import pyplot as plt

f, t, S = sc.compute_spectrogram('X9', 'IR1', 'HHE', starttime, endtime, 7200, freq_max=12.5)
fig = plt.figure(figsize=(9, 7))
plot_spct_series(S, f, t, log_scale=True)

Error Traceback

NameError                                 Traceback (most recent call last)
Cell In[1], line 4
      1 from seismic.plot.plot_spectrum import plot_spct_series
      2 from matplotlib import pyplot as plt
----> 4 f, t, S = sc.compute_spectrogram('X9', 'IR1', 'HHE', starttime, endtime, 7200, freq_max=12.5)
      5 fig = plt.figure(figsize=(9, 7))
      6 plot_spct_series(S, f, t, log_scale=True)

NameError: name 'sc' is not defined

SeisMIC Version?

last

Operating System?

OS/X Ventura

Python Version?

Python 3.11.5

Installation Method?

conda

Problem

Hello, Peter:
I am very interested in your work and appreciate the program you have released. But I have noticed that in certain frequency bands, especially in the high-frequency range, the calculated cross-correlation function shows some issues at Lag Time = -MAX, as shown in the figure below.
I am wondering if this is related to the use of the frequency domain method for calculating cross-correlation in the "correlate.py" file's "_pxcorr_matrix" function. Therefore, I am seeking your advice to understand whether this issue might affect the dvv results for high-frequency signals.

Proposed solution

No response

Problem

Lately, in particular, scipy and numpy have the habit of changing their syntax.

Proposed solution

Fix version codes in the various dependency files.

Hello, I collected some data by myself with a three-component seismograph, the original format of the data is DLD, the data can be converted into miniseed format using sololite, but there is no network information in the data, may I ask how to edit the network when I edit and sort the files? Can i define one myself?