manoch-iris
commented
1 year ago Dear Elisabetta,
Thank you for your contacting us. That is definitely a package conflict issue. Unfortunately I will be out until July 5th and I will check it then. Would you please make sure all your Python packages are up to date? Also, please confirm that you have not changed any of the parameter values in the param file. After these checks, if you still have the same issue, please reconfirm the versions of the Python modules and we will test and get back to you after July 4th. Thanks
On Wed, Jun 28, 2023 at 8:34 PM elidana @.***> wrote:
Dear EarthScope data management center,
I am trying to use one of your great tools to perform some noise analysis of our data. I am interested in using the PDF/PSD package of the noise-toolkit, and I was successfully using that last year.
I was trying again today, and the main change since last time I was using is some upgrades on my local host.
When I run the "default" command suggested in the wiki quickinstall, I get the following error, which I think is related to the version of matplotlilb.
This is the command and error
@. noise-toolkit]$ python bin/ntk_computePSD.py net=TA sta=O18A loc=DASH chan=BHZ start=2008-08-14T12:00:00 end=2008-08-14T13:00:00 xtype=period plot=1 [INFO] script: ntk_computePSD.py v.2.1.0 8 args: ['bin/ntk_computePSD.py', 'net=TA', 'sta=O18A', 'loc=DASH', 'chan=BHZ', 'start=2008-08-14T12:00:00', 'end=2008-08-14T13:00:00', 'xtype=period', 'plot=1'] [INFO] Use Internet to get metadata: True [INFO] accessing logon client as @.
**** START ****
[INFO] Requesting TA.--.O18A.BHZ from 2008-08-14T12:00:00 to 2008-08-14T13:00:00 [INFO] sending request to fedcatalog: http://service.iris.edu/irisws/fedcatalog/1/query?net=TA&sta=O18A&loc=--&cha=BHZ&targetservice=dataselect&level=channel&format=request&start=2008-08-14T12:00:00&end=2008-08-14T13:00:00&includeoverlaps=false&nodata=404 [INFO] from the IRISDMC,http://ds.iris.edu data center [INFO] Initiating fedcatalog request for IRISDMC [INFO] dataselect service is http://service.iris.edu/fdsnws/dataselect/1/ [INFO] Requesting data from IRISDMC_0 via http://service.iris.edu
[INFO] Reading TA.O18A.--.BHZ from 2008-08-14 12:00:00.0 to 2008-08-14 13:00:00.0 from FDSN stream [INFO] ntk_computePSD.py 1 Trace(s) in Stream: TA.O18A..BHZ | 2008-08-14T12:00:00.000000Z - 2008-08-14T13:00:00.000000Z | 40.0 Hz, 144001 samples 16384 Traceback (most recent call last): File "/home/elidana/tmp/aws-workshop/noise-toolkit/bin/ntk_computePSD.py", line 610, in
power, freq = csd(tr.data, tr.data, NFFT=nfft, noverlap=nfft * windlap, Fs=1. / delta, scale_byfreq=True) ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/home/elidana/.local/lib/python3.11/site-packages/matplotlib/mlab.py", line 640, in csd Pxy, freqs, = _spectral_helper(x=x, y=y, NFFT=NFFT, Fs=Fs, ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/home/elidana/.local/lib/python3.11/site-packages/matplotlib/mlab.py", line 382, in _spectral_helper result = _stride_windows(x, NFFT, noverlap) ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "/home/elidana/.local/lib/python3.11/site-packages/matplotlib/mlab.py", line 258, in _stride_windows return np.lib.stride_tricks.sliding_window_view( ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ TypeError: slice indices must be integers or None or have an index method I am using matplotlib 3.7.1 and Python 3.11.3 version.
Happy to provide further details if this is of any help.
Thanks a lot for any advice you are happy to share!
Elisabetta
— Reply to this email directly, view it on GitHub https://github.com/iris-edu/noise-toolkit/issues/4, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABL5L2UC5NH77OTIMGCGUYDXNTEQJANCNFSM6AAAAAAZX2F2I4 . You are receiving this because you are subscribed to this thread.Message ID: @.***>
-- Manochehr Bahavar Data Scientist @.***
elidana
Dear EarthScope data management center,
I am trying to use one of your great tools to perform some noise analysis of seismic data. I am interested in using the PDF/PSD package of the noise-toolkit, and I was successfully using that last year.
I was trying again today, and the main change since last time I was using is some upgrades on my local host.
When I run the "default" command suggested in the wiki quickinstall, I get the following error, which I think is related my version of matplotlilb.
This is the command and error
I am using
matplotlib 3.7.1andPython 3.11.3version.Happy to provide further details if this is of any help.
Thanks a lot for any advice you are happy to share!
Elisabetta