iris-edu / noise-toolkit Goto Github PK

View Code? Open in Web Editor NEW

30.0 6.0 10.0 18.36 MB

A code repository for the IRIS DMC Noise Toolkit data product.

Home Page: https://ds.iris.edu/ds/products/noise-toolkit/

Python 100.00%

noise-toolkit iris-dmc waveform-data microseism-energy psd

noise-toolkit's People

Stargazers

Watchers

Forkers

gomexnez awacero dgcornwell winteryao7 lukman-abdn19 yaoshangheng kzt206 elidana catmou dp01101100

noise-toolkit's Issues

PSD computation

HI,

I am using noise toolkit to compute PSDs, it was possible for me to produce hourly and daily PSD plots using window length = 3600s and 86400s respectively in the parameter file. But when I tried the same method to generate monthly plot, it didn't work. I would be very happy and thankful to get workable suggestions on how to compute monthly PSDs and a summarised general plot, thank you

"fileTag" command line argument for NTK

My data is split into multiple (daily) directories in the noise_toolkit/data directory, and when I want to run the scripts for a new day I need to manually change the param file.
Could a command line argument be added to change the filetag directory and filename to point to the preferred file location? This way I may be able to automate the process as well.

ntk_computePSD.py list slicing error

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

[elidana@fedora 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 [email protected]

******************** 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 <module>
    power, freq = csd(tr.data, tr.data, NFFT=nfft, noverlap=nfft * windlap, Fs=1. / delta, scale_by_freq=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

Sir, how to use this toolkit to plot this picture that has lots of PSDs?

Sir, When I plot the picture,just like this

But I want to plot picture like this:

Sir ,could you tell me how to use the toolkit to get it?

Best Wishes!

Questions

Dear All,
I started to use the noise-toolkit to compute polarization, while working downloading from IRIS ok, however we have some offline data and we would like to apply your tool in our off-line data.

1 hour and 24 hour long data are located at:
D:\Inves\ML\noise-toolkit-master\data\MSEED\

My first question is related to param configuration file on computePolarization.py:

requestClient = "FILES"
fromFileOnly = True
fileTag = os.path.join(r'D:\Inves\ML\noise-toolkit-master\data\MSEED\' , "mseed", "*.mseed")

Is it ok?

Another question is related to my XML response files, I have separated files for each channel, however I saw over the examples that in just one XML you managed the three channels, to fulfill this requirement I am using the [http://service.iris.edu/fdsnws/station/1/query?sta=SOEV&startafter=2020-06-01&endbefore=2020-06-02&minlatitude=-20&maxlatitude=-9&minlongitude=-70&maxlongitude=-63&level=channel] page to 'merge' all channel in one XML, is it ok?

Thanks in advance.

Tonino

Recommend Projects

React

A declarative, efficient, and flexible JavaScript library for building user interfaces.
Vue.js

🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
Typescript

TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
TensorFlow

An Open Source Machine Learning Framework for Everyone
Django

The Web framework for perfectionists with deadlines.
Laravel

A PHP framework for web artisans
D3

Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

javascript

JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
web

Some thing interesting about web. New door for the world.
server

A server is a program made to process requests and deliver data to clients.
Machine learning

Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Visualization

Some thing interesting about visualization, use data art
Game

Some thing interesting about game, make everyone happy.

Recommend Org

Facebook

We are working to build community through open source technology. NB: members must have two-factor auth.
Microsoft

Open source projects and samples from Microsoft.
Google

Google ❤️ Open Source for everyone.
Alibaba

Alibaba Open Source for everyone
D3

Data-Driven Documents codes.
Tencent

China tencent open source team.