MSNoise January 2014

msnoise@lists.observatory.be

5 participants
7 discussions

6 stations added to station list, but only 2 stations cross-correlated.

by Abbott, Robert E

Thomas, My previous thread about the AK.WAT data is only of academic interest to me as it was just a test data set downloaded form the IRIS DMC. I will be happy to work with you to find out why it is happening, however. The good news is that I am able to go through the process all the way through the stacking portion on my own data. The bad news is that it only uses 2 of the 6 stations. When I populate the station table, all six stations are located and are verified to be in the table. When I scan the archive, however, only the first two stations are present in the resultant 'data_availability' table. All the mseed files were produced identically and are in BUD format. Thanks, Rob +----------------------------+ Robert E. Abbott, Ph.D. Sandia National Laboratories Geophysics Department MS 0750 P.O. Box 5800 Albuquerque, NM 87185-0750 (505) 845-0266 +----------------------------+

11 years, 12 months

Empty Filters Table

by Abbott, Robert E

Thomas, I get the same error referenced by Xiao in this thread, http://mailman-as.oma.be/pipermail/msnoise/2013-October/000033.html. Namely, in the Configurator UI, I cannot seem to edit the filters table and this leads to the errors in s02new_jobs.py (I believe). Traceback (most recent call last): File "s02new_jobs.py", line 89, in <module> if day != old_day and day_pairs != []: NameError: name 'day' is not defined I did not see a definitive fix in the following emails. If I have to edit the filters table manually (in MySQL I assume), can you supply me with an example to guide me? A screenshot is attached. I am on OS X 10.8.5. Thanks, Rob +----------------------------+ Robert E. Abbott, Ph.D. Sandia National Laboratories Geophysics Department MS 0750 P.O. Box 5800 Albuquerque, NM 87185-0750 (505) 845-0266 +----------------------------+

12 years, 1 month

Re: [MSNoise] MSNoise Digest, Vol 3, Issue 11

by Esteban Chaves

Hi Rob, The way to add the filters in your input parameters is using a sqlite manager for firefox: https://addons.mozilla.org/en-US/firefox/addon/sqlite-manager/ Then you can open the msnoise.sqlite file, which is in your main directory, and edit your filters. I think, this error is not relative to the filters. You have to check that your stations and network are well established on this input parameters and then your waveforms must be identify as "T" instead "D". Traceback (most recent call last): File "s02new_jobs.py", line 89, in <module> if day != old_day and day_pairs != []: NameError: name 'day' is not defined Esteban. Esteban J. Chaves PhD Student in Seismology Earth and Planetary Sciences Earth and Marine Science Building office C317B University of California, Santa Cruz On Mon, Jan 27, 2014 at 12:16 PM, <msnoise-request(a)mailman-as.oma.be> wrote: > Send MSNoise mailing list submissions to > msnoise(a)mailman-as.oma.be > > To subscribe or unsubscribe via the World Wide Web, visit > http://mailman-as.oma.be/mailman/listinfo/msnoise > or, via email, send a message with subject or body 'help' to > msnoise-request(a)mailman-as.oma.be > > You can reach the person managing the list at > msnoise-owner(a)mailman-as.oma.be > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of MSNoise digest..." > > > Today's Topics: > > 1. Empty Filters Table (Abbott, Robert E) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 27 Jan 2014 20:15:29 +0000 > From: "Abbott, Robert E" <reabbot(a)sandia.gov> > To: "Python Package for Monitoring Seismic Velocity Changes using > Ambient Seismic Noise" <msnoise(a)mailman-as.oma.be> > Subject: [MSNoise] Empty Filters Table > Message-ID: <CF0C0C6F.4DAF%reabbot(a)sandia.gov> > Content-Type: text/plain; charset="us-ascii" > > Thomas, > > I get the same error referenced by Xiao in this thread, > http://mailman-as.oma.be/pipermail/msnoise/2013-October/000033.html. > > Namely, in the Configurator UI, I cannot seem to edit the filters table > and this leads to the errors in s02new_jobs.py (I believe). > > Traceback (most recent call last): > File "s02new_jobs.py", line 89, in <module> > if day != old_day and day_pairs != []: > NameError: name 'day' is not defined > > I did not see a definitive fix in the following emails. If I have to edit > the filters table manually (in MySQL I assume), can you supply me with an > example to guide me? > > A screenshot is attached. I am on OS X 10.8.5. > > Thanks, > Rob > > +----------------------------+ > Robert E. Abbott, Ph.D. > Sandia National Laboratories > Geophysics Department MS 0750 > P.O. Box 5800 > Albuquerque, NM 87185-0750 > (505) 845-0266 > +----------------------------+ > >

12 years, 1 month

Re: [MSNoise] MSNoise Digest, Vol 3, Issue 7

by Esteban Chaves

Dear Thomas, This problem is relative to a DAS configuration. Here is a program called msmod to change the miniseed headers: https://seiscode.iris.washington.edu/search?utf8=%E2%9C%93&q=msmod&all_word… Esteban, Esteban J. Chaves PhD Student in Seismology Earth and Planetary Sciences Earth and Marine Science Building office C317B University of California, Santa Cruz On Wed, Jan 22, 2014 at 7:54 AM, <msnoise-request(a)mailman-as.oma.be> wrote: > Send MSNoise mailing list submissions to > msnoise(a)mailman-as.oma.be > > To subscribe or unsubscribe via the World Wide Web, visit > http://mailman-as.oma.be/mailman/listinfo/msnoise > or, via email, send a message with subject or body 'help' to > msnoise-request(a)mailman-as.oma.be > > You can reach the person managing the list at > msnoise-owner(a)mailman-as.oma.be > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of MSNoise digest..." > > > Today's Topics: > > 1. Re: [EXTERNAL] Re: Possible to use data from ref2mseed? > (Abbott, Robert E) > 2. Re: [EXTERNAL] Re: Possible to use data from ref2mseed? > (Thomas Lecocq) > 3. Re: [EXTERNAL] Re: Possible to use data from ref2mseed? > (Abbott, Robert E) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 22 Jan 2014 15:20:20 +0000 > From: "Abbott, Robert E" <reabbot(a)sandia.gov> > To: "Python Package for Monitoring Seismic Velocity Changes using > Ambient Seismic Noise" <msnoise(a)mailman-as.oma.be> > Subject: Re: [MSNoise] [EXTERNAL] Re: Possible to use data from > ref2mseed? > Message-ID: <CF052D8D.4BDC%reabbot(a)sandia.gov> > Content-Type: text/plain; charset="iso-8859-1" > > The filenames take the serial number of the DAS automatically. I do not > think this behavior can be changed during conversion The header values > are populated by a parameter file input by the user. Below is a few lines > of the parameter file: > > #das;refchan;refstrm;netcode;station;channel;loccode;encoding;samplerate;ga > in > BACB; 1; 1; PA; CE1; HHZ; 00; INT32; 125; x1 > BACB; 2; 1; PA; CE1; HHN; 00; INT32; 125; x1 > BACB; 3; 1; PA; CE1; HHE; 00; INT32; 125; x1 > B9F0; 1; 1; PA; R1A; HHZ; 00; INT32; 125; x1 > > > I can change the 'station' parameter to match the 'das' parameter, or > write a script to batch rename the files to match the 'station'. I think I > would prefer to do the latter, as the DAS serial number has no geographic > significance, unlike the station value. For instance R1A, R1B are on Ring > 1 of the array and R2A, R2B are Ring 2, etc. CE1 is the center element of > the array. > > Do you want me to send you renamed data? > > Best, > Rob > > > > +----------------------------+ > > Robert E. Abbott, Ph.D. > Sandia National Laboratories > Geophysics Department MS 0750 > P.O. Box 5800 > Albuquerque, NM 87185-0750 > (505) 845-0266 > +----------------------------+ > > > > > > On 1/22/14 1:43 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> wrote: > > >Hi Rob, > > > >I can read those files. I'm a little concerned by the fact the station > >name in the filename is not the station name in the actual data > >contained in the file. That is strange, and this could be a problem. > >Indeed, the station table is populated from the filenames, and the > >archive is scanned searching for selected stations in the configuration. > >But, the scan reads data within the file and stores metadata in the > >database, and these metadata have different station names. > > > >Do you have an idea why the names don't match ? > > > >Thomas > > > >Le 20/01/2014 23:30, Abbott, Robert E a ?crit : > >> Thomas, > >> > >> New files are on the way to you. I changed conversion algorithm from > >> 'ref2mseed' to 'rt2ms' (also from PASSCAL, but more modern). rt2ms > >> produces much the same naming scheme except it has a YYYY format instead > >> of a YY, among other changes. > >> > >> Hope this helps, > >> Rob > >> > >> +----------------------------+ > >> Robert E. Abbott, Ph.D. > >> Sandia National Laboratories > >> Geophysics Department MS 0750 > >> P.O. Box 5800 > >> Albuquerque, NM 87185-0750 > >> (505) 845-0266 > >> +----------------------------+ > >> > >> > >> > >> > >> > >> On 1/14/14 2:35 PM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> wrote: > >> > >>> Rob, > >>> > >>> Le 14.01.2014 21:40, Abbott, Robert E a ?crit : > >>>> Thanks, Thomas. I am on travel right now, so I will have to get you > >>>> the parameters Monday. There is a finite chance I screwed up and gave > >>>> you SEG-Y files, they are named similarly. If that is the case, I > >>>> apologize for the confusion. > >>> I tried reading SEG-Y but it doesn't work either. Let's wait till > >>>Monday > >>> to check ! > >>> > >>> Cheers > >>> > >>> Thomas > >>> > >>>> Sent from my iPad > >>>> > >>>>> On Jan 14, 2014, at 10:27 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> > >>>>> wrote: > >>>>> > >>>>> Hi Rob, > >>>>> > >>>>> OK, got your files. I can modify the code to automatically find the > >>>>> stations, that's OK. > >>>>> > >>>>> But, I can't read them in Python. Although you state it's MiniSeed > >>>>> (and I believe you :-) ), I can't find the proper parameters to read > >>>>> them. Did you pass any argument to the converter ? Are you able to > >>>>> open them using another tool ? > >>>>> > >>>>> Cheers, > >>>>> > >>>>> Thomas > >>>>> > >>>>> Le 10.01.2014 18:18, Abbott, Robert E a ?crit : > >>>>>> Thomas, > >>>>>> Thanks for very much for your help. I will send you some data, but 3 > >>>>>> days, > >>>>>> 3 stations is quite a large data volume as my sample rate is 200 Hz. > >>>>>> Do > >>>>>> you have a dropbox.com account? That has worked in the past for me. > >>>>>> Alternatively, I could send just an hour or two of data per day... > >>>>>> Thanks, > >>>>>> Rob > >>>>>> +----------------------------+ > >>>>>> Robert E. Abbott, Ph.D. > >>>>>> Sandia National Laboratories > >>>>>> Geophysics Department MS 0750 > >>>>>> P.O. Box 5800 > >>>>>> Albuquerque, NM 87185-0750 > >>>>>> (505) 845-0266 > >>>>>> +----------------------------+ > >>>>>>> On 1/10/14 9:06 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> wrote: > >>>>>>> Hi Rob, > >>>>>>> Yes, it is definitively possible. Highly dependant on your python > >>>>>>> skills, I'm afraid... > >>>>>>> The data structure will be something like : > >>>>>>> data_structure['ref2mseed'] = "RDAY.LOC/YEAR.DAY.HOUR.MIN.STA.CHAN" > >>>>>>> but: > >>>>>>> I assume yy is a 2digit year format ? Then YEAR will not work, we > >>>>>>> need > >>>>>>> to put YY and define it in the reader (in populate station table) > >>>>>>> I set STA to your instrument number, which I suppose are located at > >>>>>>> a > >>>>>>> unique position within this very archive. > >>>>>>> I can help writing the reader part (populate station table). I'm > >>>>>>>not > >>>>>>> quite sure how the scan_archive process will (dis)like such a bad > >>>>>>> structure (no net/ or sta/ folders)... > >>>>>>> If you provide me some data, say 3 days, 3 stations, I can test the > >>>>>>> code > >>>>>>> on it... > >>>>>>> Best regards from Brussels, > >>>>>>> Thomas > >>>>>>> Le 09/01/2014 22:43, Abbott, Robert E a ?crit : > >>>>>>>> Greetings, > >>>>>>>> Is it possible to use MSNoise on data generated from the PASSCAL > >>>>>>>> Instrument Center utility "ref2mseed" (the data is from Reftek 130 > >>>>>>>> digitizers). The data is named like this: > >>>>>>>> Rddd.0S/yy.ddd.hh.mm.ss.iiii.c > >>>>>>>> where > >>>>>>>> yy -- year of the record > >>>>>>>> ddd -- julian day of the record > >>>>>>>> S -- data stream of the instrument (usually 1) > >>>>>>>> hh -- hour the record was recorded > >>>>>>>> mm -- minute the record was recorded > >>>>>>>> iiii -- instrument number > >>>>>>>> c -- channel number > >>>>>>>> How would I format the data_structure parameter> > >>>>>>>> Thanks, > >>>>>>>> Rob > >>>>>>>> +----------------------------+ > >>>>>>>> Robert E. Abbott, Ph.D. > >>>>>>>> Sandia National Laboratories > >>>>>>>> Geophysics Department MS 0750 > >>>>>>>> P.O. Box 5800 > >>>>>>>> Albuquerque, NM 87185-0750 > >>>>>>>> (505) 845-0266 > >>>>>>>> +----------------------------+ > >>>>>>>> _______________________________________________ > >>>>>>>> MSNoise mailing list > >>>>>>>> MSNoise(a)mailman-as.oma.be > >>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>>>>> _______________________________________________ > >>>>>>> MSNoise mailing list > >>>>>>> MSNoise(a)mailman-as.oma.be > >>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>>>> _______________________________________________ > >>>>>> MSNoise mailing list > >>>>>> MSNoise(a)mailman-as.oma.be > >>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>>> _______________________________________________ > >>>>> MSNoise mailing list > >>>>> MSNoise(a)mailman-as.oma.be > >>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>> _______________________________________________ > >>>> MSNoise mailing list > >>>> MSNoise(a)mailman-as.oma.be > >>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>> _______________________________________________ > >>> MSNoise mailing list > >>> MSNoise(a)mailman-as.oma.be > >>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >> _______________________________________________ > >> MSNoise mailing list > >> MSNoise(a)mailman-as.oma.be > >> http://mailman-as.oma.be/mailman/listinfo/msnoise > > > >_______________________________________________ > >MSNoise mailing list > >MSNoise(a)mailman-as.oma.be > >http://mailman-as.oma.be/mailman/listinfo/msnoise > > > > ------------------------------ > > Message: 2 > Date: Wed, 22 Jan 2014 16:27:08 +0100 > From: Thomas Lecocq <thomas.lecocq(a)oma.be> > To: msnoise(a)mailman-as.oma.be > Subject: Re: [MSNoise] [EXTERNAL] Re: Possible to use data from > ref2mseed? > Message-ID: <52DFE34C.8020206(a)oma.be> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Rob, > Le 22/01/2014 16:20, Abbott, Robert E a ?crit : > > The filenames take the serial number of the DAS automatically. I do not > > think this behavior can be changed during conversion The header values > > are populated by a parameter file input by the user. Below is a few lines > > of the parameter file: > > > > > #das;refchan;refstrm;netcode;station;channel;loccode;encoding;samplerate;ga > > in > > BACB; 1; 1; PA; CE1; HHZ; 00; INT32; 125; x1 > > BACB; 2; 1; PA; CE1; HHN; 00; INT32; 125; x1 > > BACB; 3; 1; PA; CE1; HHE; 00; INT32; 125; x1 > > B9F0; 1; 1; PA; R1A; HHZ; 00; INT32; 125; x1 > > > > > > I can change the 'station' parameter to match the 'das' parameter, or > > write a script to batch rename the files to match the 'station'. I think > I > > would prefer to do the latter, as the DAS serial number has no geographic > > significance, unlike the station value. For instance R1A, R1B are on Ring > > 1 of the array and R2A, R2B are Ring 2, etc. CE1 is the center element of > > the array. > Yes, that would make sense, even for pure archive maintenance > considerations > > > > Do you want me to send you renamed data? > Yes, go ahead ! > > Thom > > > > > Best, > > Rob > > > > > > > > +----------------------------+ > > > > Robert E. Abbott, Ph.D. > > Sandia National Laboratories > > Geophysics Department MS 0750 > > P.O. Box 5800 > > Albuquerque, NM 87185-0750 > > (505) 845-0266 > > +----------------------------+ > > > > > > > > > > > > On 1/22/14 1:43 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> wrote: > > > >> Hi Rob, > >> > >> I can read those files. I'm a little concerned by the fact the station > >> name in the filename is not the station name in the actual data > >> contained in the file. That is strange, and this could be a problem. > >> Indeed, the station table is populated from the filenames, and the > >> archive is scanned searching for selected stations in the configuration. > >> But, the scan reads data within the file and stores metadata in the > >> database, and these metadata have different station names. > >> > >> Do you have an idea why the names don't match ? > >> > >> Thomas > >> > >> Le 20/01/2014 23:30, Abbott, Robert E a ?crit : > >>> Thomas, > >>> > >>> New files are on the way to you. I changed conversion algorithm from > >>> 'ref2mseed' to 'rt2ms' (also from PASSCAL, but more modern). rt2ms > >>> produces much the same naming scheme except it has a YYYY format > instead > >>> of a YY, among other changes. > >>> > >>> Hope this helps, > >>> Rob > >>> > >>> +----------------------------+ > >>> Robert E. Abbott, Ph.D. > >>> Sandia National Laboratories > >>> Geophysics Department MS 0750 > >>> P.O. Box 5800 > >>> Albuquerque, NM 87185-0750 > >>> (505) 845-0266 > >>> +----------------------------+ > >>> > >>> > >>> > >>> > >>> > >>> On 1/14/14 2:35 PM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> wrote: > >>> > >>>> Rob, > >>>> > >>>> Le 14.01.2014 21:40, Abbott, Robert E a ?crit : > >>>>> Thanks, Thomas. I am on travel right now, so I will have to get you > >>>>> the parameters Monday. There is a finite chance I screwed up and gave > >>>>> you SEG-Y files, they are named similarly. If that is the case, I > >>>>> apologize for the confusion. > >>>> I tried reading SEG-Y but it doesn't work either. Let's wait till > >>>> Monday > >>>> to check ! > >>>> > >>>> Cheers > >>>> > >>>> Thomas > >>>> > >>>>> Sent from my iPad > >>>>> > >>>>>> On Jan 14, 2014, at 10:27 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be > > > >>>>>> wrote: > >>>>>> > >>>>>> Hi Rob, > >>>>>> > >>>>>> OK, got your files. I can modify the code to automatically find the > >>>>>> stations, that's OK. > >>>>>> > >>>>>> But, I can't read them in Python. Although you state it's MiniSeed > >>>>>> (and I believe you :-) ), I can't find the proper parameters to read > >>>>>> them. Did you pass any argument to the converter ? Are you able to > >>>>>> open them using another tool ? > >>>>>> > >>>>>> Cheers, > >>>>>> > >>>>>> Thomas > >>>>>> > >>>>>> Le 10.01.2014 18:18, Abbott, Robert E a ?crit : > >>>>>>> Thomas, > >>>>>>> Thanks for very much for your help. I will send you some data, but > 3 > >>>>>>> days, > >>>>>>> 3 stations is quite a large data volume as my sample rate is 200 > Hz. > >>>>>>> Do > >>>>>>> you have a dropbox.com account? That has worked in the past for > me. > >>>>>>> Alternatively, I could send just an hour or two of data per day... > >>>>>>> Thanks, > >>>>>>> Rob > >>>>>>> +----------------------------+ > >>>>>>> Robert E. Abbott, Ph.D. > >>>>>>> Sandia National Laboratories > >>>>>>> Geophysics Department MS 0750 > >>>>>>> P.O. Box 5800 > >>>>>>> Albuquerque, NM 87185-0750 > >>>>>>> (505) 845-0266 > >>>>>>> +----------------------------+ > >>>>>>>> On 1/10/14 9:06 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> wrote: > >>>>>>>> Hi Rob, > >>>>>>>> Yes, it is definitively possible. Highly dependant on your python > >>>>>>>> skills, I'm afraid... > >>>>>>>> The data structure will be something like : > >>>>>>>> data_structure['ref2mseed'] = > "RDAY.LOC/YEAR.DAY.HOUR.MIN.STA.CHAN" > >>>>>>>> but: > >>>>>>>> I assume yy is a 2digit year format ? Then YEAR will not work, we > >>>>>>>> need > >>>>>>>> to put YY and define it in the reader (in populate station table) > >>>>>>>> I set STA to your instrument number, which I suppose are located > at > >>>>>>>> a > >>>>>>>> unique position within this very archive. > >>>>>>>> I can help writing the reader part (populate station table). I'm > >>>>>>>> not > >>>>>>>> quite sure how the scan_archive process will (dis)like such a bad > >>>>>>>> structure (no net/ or sta/ folders)... > >>>>>>>> If you provide me some data, say 3 days, 3 stations, I can test > the > >>>>>>>> code > >>>>>>>> on it... > >>>>>>>> Best regards from Brussels, > >>>>>>>> Thomas > >>>>>>>> Le 09/01/2014 22:43, Abbott, Robert E a ?crit : > >>>>>>>>> Greetings, > >>>>>>>>> Is it possible to use MSNoise on data generated from the PASSCAL > >>>>>>>>> Instrument Center utility "ref2mseed" (the data is from Reftek > 130 > >>>>>>>>> digitizers). The data is named like this: > >>>>>>>>> Rddd.0S/yy.ddd.hh.mm.ss.iiii.c > >>>>>>>>> where > >>>>>>>>> yy -- year of the record > >>>>>>>>> ddd -- julian day of the record > >>>>>>>>> S -- data stream of the instrument (usually 1) > >>>>>>>>> hh -- hour the record was recorded > >>>>>>>>> mm -- minute the record was recorded > >>>>>>>>> iiii -- instrument number > >>>>>>>>> c -- channel number > >>>>>>>>> How would I format the data_structure parameter> > >>>>>>>>> Thanks, > >>>>>>>>> Rob > >>>>>>>>> +----------------------------+ > >>>>>>>>> Robert E. Abbott, Ph.D. > >>>>>>>>> Sandia National Laboratories > >>>>>>>>> Geophysics Department MS 0750 > >>>>>>>>> P.O. Box 5800 > >>>>>>>>> Albuquerque, NM 87185-0750 > >>>>>>>>> (505) 845-0266 > >>>>>>>>> +----------------------------+ > >>>>>>>>> _______________________________________________ > >>>>>>>>> MSNoise mailing list > >>>>>>>>> MSNoise(a)mailman-as.oma.be > >>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>>>>>> _______________________________________________ > >>>>>>>> MSNoise mailing list > >>>>>>>> MSNoise(a)mailman-as.oma.be > >>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>>>>> _______________________________________________ > >>>>>>> MSNoise mailing list > >>>>>>> MSNoise(a)mailman-as.oma.be > >>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>>>> _______________________________________________ > >>>>>> MSNoise mailing list > >>>>>> MSNoise(a)mailman-as.oma.be > >>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>>> _______________________________________________ > >>>>> MSNoise mailing list > >>>>> MSNoise(a)mailman-as.oma.be > >>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>> _______________________________________________ > >>>> MSNoise mailing list > >>>> MSNoise(a)mailman-as.oma.be > >>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>> _______________________________________________ > >>> MSNoise mailing list > >>> MSNoise(a)mailman-as.oma.be > >>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >> _______________________________________________ > >> MSNoise mailing list > >> MSNoise(a)mailman-as.oma.be > >> http://mailman-as.oma.be/mailman/listinfo/msnoise > > _______________________________________________ > > MSNoise mailing list > > MSNoise(a)mailman-as.oma.be > > http://mailman-as.oma.be/mailman/listinfo/msnoise > > > > ------------------------------ > > Message: 3 > Date: Wed, 22 Jan 2014 15:53:58 +0000 > From: "Abbott, Robert E" <reabbot(a)sandia.gov> > To: "Python Package for Monitoring Seismic Velocity Changes using > Ambient Seismic Noise" <msnoise(a)mailman-as.oma.be> > Subject: Re: [MSNoise] [EXTERNAL] Re: Possible to use data from > ref2mseed? > Message-ID: <CF053777.4BF8%reabbot(a)sandia.gov> > Content-Type: text/plain; charset="Windows-1252" > > Data is on its way to you? > > +----------------------------+ > Robert E. Abbott, Ph.D. > Sandia National Laboratories > Geophysics Department MS 0750 > P.O. Box 5800 > Albuquerque, NM 87185-0750 > (505) 845-0266 > +----------------------------+ > > > > > > On 1/22/14 8:27 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> wrote: > > >Rob, > >Le 22/01/2014 16:20, Abbott, Robert E a ?crit : > >> The filenames take the serial number of the DAS automatically. I do not > >> think this behavior can be changed during conversion The header values > >> are populated by a parameter file input by the user. Below is a few > >>lines > >> of the parameter file: > >> > >> > >>#das;refchan;refstrm;netcode;station;channel;loccode;encoding;samplerate; > >>ga > >> in > >> BACB; 1; 1; PA; CE1; HHZ; 00; INT32; 125; x1 > >> BACB; 2; 1; PA; CE1; HHN; 00; INT32; 125; x1 > >> BACB; 3; 1; PA; CE1; HHE; 00; INT32; 125; x1 > >> B9F0; 1; 1; PA; R1A; HHZ; 00; INT32; 125; x1 > >> > >> > >> I can change the 'station' parameter to match the 'das' parameter, or > >> write a script to batch rename the files to match the 'station'. I > >>think I > >> would prefer to do the latter, as the DAS serial number has no > >>geographic > >> significance, unlike the station value. For instance R1A, R1B are on > >>Ring > >> 1 of the array and R2A, R2B are Ring 2, etc. CE1 is the center element > >>of > >> the array. > >Yes, that would make sense, even for pure archive maintenance > >considerations > >> > >> Do you want me to send you renamed data? > >Yes, go ahead ! > > > >Thom > > > >> > >> Best, > >> Rob > >> > >> > >> > >> +----------------------------+ > >> > >> Robert E. Abbott, Ph.D. > >> Sandia National Laboratories > >> Geophysics Department MS 0750 > >> P.O. Box 5800 > >> Albuquerque, NM 87185-0750 > >> (505) 845-0266 > >> +----------------------------+ > >> > >> > >> > >> > >> > >> On 1/22/14 1:43 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> wrote: > >> > >>> Hi Rob, > >>> > >>> I can read those files. I'm a little concerned by the fact the station > >>> name in the filename is not the station name in the actual data > >>> contained in the file. That is strange, and this could be a problem. > >>> Indeed, the station table is populated from the filenames, and the > >>> archive is scanned searching for selected stations in the > >>>configuration. > >>> But, the scan reads data within the file and stores metadata in the > >>> database, and these metadata have different station names. > >>> > >>> Do you have an idea why the names don't match ? > >>> > >>> Thomas > >>> > >>> Le 20/01/2014 23:30, Abbott, Robert E a ?crit : > >>>> Thomas, > >>>> > >>>> New files are on the way to you. I changed conversion algorithm from > >>>> 'ref2mseed' to 'rt2ms' (also from PASSCAL, but more modern). rt2ms > >>>> produces much the same naming scheme except it has a YYYY format > >>>>instead > >>>> of a YY, among other changes. > >>>> > >>>> Hope this helps, > >>>> Rob > >>>> > >>>> +----------------------------+ > >>>> Robert E. Abbott, Ph.D. > >>>> Sandia National Laboratories > >>>> Geophysics Department MS 0750 > >>>> P.O. Box 5800 > >>>> Albuquerque, NM 87185-0750 > >>>> (505) 845-0266 > >>>> +----------------------------+ > >>>> > >>>> > >>>> > >>>> > >>>> > >>>> On 1/14/14 2:35 PM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> wrote: > >>>> > >>>>> Rob, > >>>>> > >>>>> Le 14.01.2014 21:40, Abbott, Robert E a ?crit : > >>>>>> Thanks, Thomas. I am on travel right now, so I will have to get you > >>>>>> the parameters Monday. There is a finite chance I screwed up and > >>>>>>gave > >>>>>> you SEG-Y files, they are named similarly. If that is the case, I > >>>>>> apologize for the confusion. > >>>>> I tried reading SEG-Y but it doesn't work either. Let's wait till > >>>>> Monday > >>>>> to check ! > >>>>> > >>>>> Cheers > >>>>> > >>>>> Thomas > >>>>> > >>>>>> Sent from my iPad > >>>>>> > >>>>>>> On Jan 14, 2014, at 10:27 AM, "Thomas Lecocq" > >>>>>>><thomas.lecocq(a)oma.be> > >>>>>>> wrote: > >>>>>>> > >>>>>>> Hi Rob, > >>>>>>> > >>>>>>> OK, got your files. I can modify the code to automatically find the > >>>>>>> stations, that's OK. > >>>>>>> > >>>>>>> But, I can't read them in Python. Although you state it's MiniSeed > >>>>>>> (and I believe you :-) ), I can't find the proper parameters to > >>>>>>>read > >>>>>>> them. Did you pass any argument to the converter ? Are you able to > >>>>>>> open them using another tool ? > >>>>>>> > >>>>>>> Cheers, > >>>>>>> > >>>>>>> Thomas > >>>>>>> > >>>>>>> Le 10.01.2014 18:18, Abbott, Robert E a ?crit : > >>>>>>>> Thomas, > >>>>>>>> Thanks for very much for your help. I will send you some data, > >>>>>>>>but 3 > >>>>>>>> days, > >>>>>>>> 3 stations is quite a large data volume as my sample rate is 200 > >>>>>>>>Hz. > >>>>>>>> Do > >>>>>>>> you have a dropbox.com account? That has worked in the past for > >>>>>>>>me. > >>>>>>>> Alternatively, I could send just an hour or two of data per day... > >>>>>>>> Thanks, > >>>>>>>> Rob > >>>>>>>> +----------------------------+ > >>>>>>>> Robert E. Abbott, Ph.D. > >>>>>>>> Sandia National Laboratories > >>>>>>>> Geophysics Department MS 0750 > >>>>>>>> P.O. Box 5800 > >>>>>>>> Albuquerque, NM 87185-0750 > >>>>>>>> (505) 845-0266 > >>>>>>>> +----------------------------+ > >>>>>>>>> On 1/10/14 9:06 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> > wrote: > >>>>>>>>> Hi Rob, > >>>>>>>>> Yes, it is definitively possible. Highly dependant on your python > >>>>>>>>> skills, I'm afraid... > >>>>>>>>> The data structure will be something like : > >>>>>>>>> data_structure['ref2mseed'] = > >>>>>>>>>"RDAY.LOC/YEAR.DAY.HOUR.MIN.STA.CHAN" > >>>>>>>>> but: > >>>>>>>>> I assume yy is a 2digit year format ? Then YEAR will not work, we > >>>>>>>>> need > >>>>>>>>> to put YY and define it in the reader (in populate station table) > >>>>>>>>> I set STA to your instrument number, which I suppose are located > >>>>>>>>>at > >>>>>>>>> a > >>>>>>>>> unique position within this very archive. > >>>>>>>>> I can help writing the reader part (populate station table). I'm > >>>>>>>>> not > >>>>>>>>> quite sure how the scan_archive process will (dis)like such a bad > >>>>>>>>> structure (no net/ or sta/ folders)... > >>>>>>>>> If you provide me some data, say 3 days, 3 stations, I can test > >>>>>>>>>the > >>>>>>>>> code > >>>>>>>>> on it... > >>>>>>>>> Best regards from Brussels, > >>>>>>>>> Thomas > >>>>>>>>> Le 09/01/2014 22:43, Abbott, Robert E a ?crit : > >>>>>>>>>> Greetings, > >>>>>>>>>> Is it possible to use MSNoise on data generated from the PASSCAL > >>>>>>>>>> Instrument Center utility "ref2mseed" (the data is from Reftek > >>>>>>>>>>130 > >>>>>>>>>> digitizers). The data is named like this: > >>>>>>>>>> Rddd.0S/yy.ddd.hh.mm.ss.iiii.c > >>>>>>>>>> where > >>>>>>>>>> yy -- year of the record > >>>>>>>>>> ddd -- julian day of the record > >>>>>>>>>> S -- data stream of the instrument (usually > >>>>>>>>>>1) > >>>>>>>>>> hh -- hour the record was recorded > >>>>>>>>>> mm -- minute the record was recorded > >>>>>>>>>> iiii -- instrument number > >>>>>>>>>> c -- channel number > >>>>>>>>>> How would I format the data_structure parameter> > >>>>>>>>>> Thanks, > >>>>>>>>>> Rob > >>>>>>>>>> +----------------------------+ > >>>>>>>>>> Robert E. Abbott, Ph.D. > >>>>>>>>>> Sandia National Laboratories > >>>>>>>>>> Geophysics Department MS 0750 > >>>>>>>>>> P.O. Box 5800 > >>>>>>>>>> Albuquerque, NM 87185-0750 > >>>>>>>>>> (505) 845-0266 > >>>>>>>>>> +----------------------------+ > >>>>>>>>>> _______________________________________________ > >>>>>>>>>> MSNoise mailing list > >>>>>>>>>> MSNoise(a)mailman-as.oma.be > >>>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>>>>>>> _______________________________________________ > >>>>>>>>> MSNoise mailing list > >>>>>>>>> MSNoise(a)mailman-as.oma.be > >>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>>>>>> _______________________________________________ > >>>>>>>> MSNoise mailing list > >>>>>>>> MSNoise(a)mailman-as.oma.be > >>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>>>>> _______________________________________________ > >>>>>>> MSNoise mailing list > >>>>>>> MSNoise(a)mailman-as.oma.be > >>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>>>> _______________________________________________ > >>>>>> MSNoise mailing list > >>>>>> MSNoise(a)mailman-as.oma.be > >>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>>> _______________________________________________ > >>>>> MSNoise mailing list > >>>>> MSNoise(a)mailman-as.oma.be > >>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>> _______________________________________________ > >>>> MSNoise mailing list > >>>> MSNoise(a)mailman-as.oma.be > >>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>> _______________________________________________ > >>> MSNoise mailing list > >>> MSNoise(a)mailman-as.oma.be > >>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >> _______________________________________________ > >> MSNoise mailing list > >> MSNoise(a)mailman-as.oma.be > >> http://mailman-as.oma.be/mailman/listinfo/msnoise > > > >_______________________________________________ > >MSNoise mailing list > >MSNoise(a)mailman-as.oma.be > >http://mailman-as.oma.be/mailman/listinfo/msnoise > > > > ------------------------------ > > _______________________________________________ > MSNoise mailing list > MSNoise(a)mailman-as.oma.be > http://mailman-as.oma.be/mailman/listinfo/msnoise > > > End of MSNoise Digest, Vol 3, Issue 7 > ************************************* >

12 years, 1 month

Possible to use data from ref2mseed?

by Abbott, Robert E

Greetings, Is it possible to use MSNoise on data generated from the PASSCAL Instrument Center utility "ref2mseed" (the data is from Reftek 130 digitizers). The data is named like this: Rddd.0S/yy.ddd.hh.mm.ss.iiii.c where yy -- year of the record ddd -- julian day of the record S -- data stream of the instrument (usually 1) hh -- hour the record was recorded mm -- minute the record was recorded iiii -- instrument number c -- channel number How would I format the data_structure parameter> Thanks, Rob +----------------------------+ Robert E. Abbott, Ph.D. Sandia National Laboratories Geophysics Department MS 0750 P.O. Box 5800 Albuquerque, NM 87185-0750 (505) 845-0266 +----------------------------+

12 years, 1 month

Re: [MSNoise] MSNoise Digest, Vol 3, Issue 3

by Esteban Chaves

Hello Rob and Thomas I'm using Reftek's RT-130, and what I do is convert the data from reftek to miniseed using rt_mseed from the Reftek utilities. In case you need to correct the header of your miniseed data you can use a program called msmod. Then, you can write an script to change the external name of your data, for example: your miniseed: Rddd.0S/yy.ddd.hh.mm.ss.iiii.c ; and asuming a station called INDI and XX as network; Script: for i in yy.ddd.hh*; do j=`echo $i | sed 's/yy.ddd.hh*/XX.INDI.00.HHZ.yy.ddd.hh/'` mv $i $j done Esteban. UCSC. Esteban J. Chaves PhD Student in Seismology Earth and Planetary Sciences Earth and Marine Science Building office C317B University of California, Santa Cruz On Wed, Jan 15, 2014 at 4:00 AM, <msnoise-request(a)mailman-as.oma.be> wrote: > Send MSNoise mailing list submissions to > msnoise(a)mailman-as.oma.be > > To subscribe or unsubscribe via the World Wide Web, visit > http://mailman-as.oma.be/mailman/listinfo/msnoise > or, via email, send a message with subject or body 'help' to > msnoise-request(a)mailman-as.oma.be > > You can reach the person managing the list at > msnoise-owner(a)mailman-as.oma.be > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of MSNoise digest..." > > > Today's Topics: > > 1. Re: [EXTERNAL] Re: Possible to use data from ref2mseed? > (Thomas Lecocq) > 2. Re: [EXTERNAL] Re: Possible to use data from ref2mseed? > (Abbott, Robert E) > 3. Re: [EXTERNAL] Re: Possible to use data from ref2mseed? > (Thomas Lecocq) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 14 Jan 2014 18:26:43 +0100 > From: Thomas Lecocq <thomas.lecocq(a)oma.be> > To: Python Package for Monitoring Seismic Velocity Changes using > Ambient Seismic Noise <msnoise(a)mailman-as.oma.be> > Subject: Re: [MSNoise] [EXTERNAL] Re: Possible to use data from > ref2mseed? > Message-ID: <b79674c9011399c24e9a6ec42f9d195c(a)oma.be> > Content-Type: text/plain; charset=UTF-8; format=flowed > > Hi Rob, > > OK, got your files. I can modify the code to automatically find the > stations, that's OK. > > But, I can't read them in Python. Although you state it's MiniSeed (and > I believe you :-) ), I can't find the proper parameters to read them. > Did you pass any argument to the converter ? Are you able to open them > using another tool ? > > Cheers, > > Thomas > > Le 10.01.2014 18:18, Abbott, Robert E a ?crit?: > > Thomas, > > > > Thanks for very much for your help. I will send you some data, but 3 > > days, > > 3 stations is quite a large data volume as my sample rate is 200 Hz. Do > > you have a dropbox.com account? That has worked in the past for me. > > Alternatively, I could send just an hour or two of data per day... > > > > Thanks, > > Rob > > > > > > > > +----------------------------+ > > Robert E. Abbott, Ph.D. > > Sandia National Laboratories > > Geophysics Department MS 0750 > > P.O. Box 5800 > > Albuquerque, NM 87185-0750 > > (505) 845-0266 > > +----------------------------+ > > > > > > > > > > > > On 1/10/14 9:06 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> wrote: > > > >> Hi Rob, > >> > >> Yes, it is definitively possible. Highly dependant on your python > >> skills, I'm afraid... > >> > >> The data structure will be something like : > >> > >> data_structure['ref2mseed'] = "RDAY.LOC/YEAR.DAY.HOUR.MIN.STA.CHAN" > >> > >> but: > >> > >> I assume yy is a 2digit year format ? Then YEAR will not work, we need > >> to put YY and define it in the reader (in populate station table) > >> I set STA to your instrument number, which I suppose are located at a > >> unique position within this very archive. > >> > >> I can help writing the reader part (populate station table). I'm not > >> quite sure how the scan_archive process will (dis)like such a bad > >> structure (no net/ or sta/ folders)... > >> > >> If you provide me some data, say 3 days, 3 stations, I can test the > >> code > >> on it... > >> > >> Best regards from Brussels, > >> > >> Thomas > >> > >> Le 09/01/2014 22:43, Abbott, Robert E a ?crit : > >>> Greetings, > >>> > >>> Is it possible to use MSNoise on data generated from the PASSCAL > >>> Instrument Center utility "ref2mseed" (the data is from Reftek 130 > >>> digitizers). The data is named like this: > >>> > >>> Rddd.0S/yy.ddd.hh.mm.ss.iiii.c > >>> where > >>> yy -- year of the record > >>> ddd -- julian day of the record > >>> S -- data stream of the instrument (usually 1) > >>> hh -- hour the record was recorded > >>> mm -- minute the record was recorded > >>> iiii -- instrument number > >>> c -- channel number > >>> > >>> How would I format the data_structure parameter> > >>> > >>> Thanks, > >>> Rob > >>> > >>> +----------------------------+ > >>> Robert E. Abbott, Ph.D. > >>> Sandia National Laboratories > >>> Geophysics Department MS 0750 > >>> P.O. Box 5800 > >>> Albuquerque, NM 87185-0750 > >>> (505) 845-0266 > >>> +----------------------------+ > >>> > >>> _______________________________________________ > >>> MSNoise mailing list > >>> MSNoise(a)mailman-as.oma.be > >>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >> > >> _______________________________________________ > >> MSNoise mailing list > >> MSNoise(a)mailman-as.oma.be > >> http://mailman-as.oma.be/mailman/listinfo/msnoise > > > > _______________________________________________ > > MSNoise mailing list > > MSNoise(a)mailman-as.oma.be > > http://mailman-as.oma.be/mailman/listinfo/msnoise > > > ------------------------------ > > Message: 2 > Date: Tue, 14 Jan 2014 20:40:20 +0000 > From: "Abbott, Robert E" <reabbot(a)sandia.gov> > To: "Python Package for Monitoring Seismic Velocity Changes using > Ambient Seismic Noise" <msnoise(a)mailman-as.oma.be> > Subject: Re: [MSNoise] [EXTERNAL] Re: Possible to use data from > ref2mseed? > Message-ID: <B30BCD1F-6E70-44B7-8C85-DD6954C361EC(a)sandia.gov> > Content-Type: text/plain; charset="iso-8859-1" > > Thanks, Thomas. I am on travel right now, so I will have to get you the > parameters Monday. There is a finite chance I screwed up and gave you SEG-Y > files, they are named similarly. If that is the case, I apologize for the > confusion. > > Sent from my iPad > > > On Jan 14, 2014, at 10:27 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> > wrote: > > > > Hi Rob, > > > > OK, got your files. I can modify the code to automatically find the > stations, that's OK. > > > > But, I can't read them in Python. Although you state it's MiniSeed (and > I believe you :-) ), I can't find the proper parameters to read them. Did > you pass any argument to the converter ? Are you able to open them using > another tool ? > > > > Cheers, > > > > Thomas > > > > Le 10.01.2014 18:18, Abbott, Robert E a ?crit : > >> Thomas, > >> Thanks for very much for your help. I will send you some data, but 3 > days, > >> 3 stations is quite a large data volume as my sample rate is 200 Hz. Do > >> you have a dropbox.com account? That has worked in the past for me. > >> Alternatively, I could send just an hour or two of data per day... > >> Thanks, > >> Rob > >> +----------------------------+ > >> Robert E. Abbott, Ph.D. > >> Sandia National Laboratories > >> Geophysics Department MS 0750 > >> P.O. Box 5800 > >> Albuquerque, NM 87185-0750 > >> (505) 845-0266 > >> +----------------------------+ > >>> On 1/10/14 9:06 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> wrote: > >>> Hi Rob, > >>> Yes, it is definitively possible. Highly dependant on your python > >>> skills, I'm afraid... > >>> The data structure will be something like : > >>> data_structure['ref2mseed'] = "RDAY.LOC/YEAR.DAY.HOUR.MIN.STA.CHAN" > >>> but: > >>> I assume yy is a 2digit year format ? Then YEAR will not work, we need > >>> to put YY and define it in the reader (in populate station table) > >>> I set STA to your instrument number, which I suppose are located at a > >>> unique position within this very archive. > >>> I can help writing the reader part (populate station table). I'm not > >>> quite sure how the scan_archive process will (dis)like such a bad > >>> structure (no net/ or sta/ folders)... > >>> If you provide me some data, say 3 days, 3 stations, I can test the > code > >>> on it... > >>> Best regards from Brussels, > >>> Thomas > >>> Le 09/01/2014 22:43, Abbott, Robert E a ?crit : > >>>> Greetings, > >>>> Is it possible to use MSNoise on data generated from the PASSCAL > >>>> Instrument Center utility "ref2mseed" (the data is from Reftek 130 > >>>> digitizers). The data is named like this: > >>>> Rddd.0S/yy.ddd.hh.mm.ss.iiii.c > >>>> where > >>>> yy -- year of the record > >>>> ddd -- julian day of the record > >>>> S -- data stream of the instrument (usually 1) > >>>> hh -- hour the record was recorded > >>>> mm -- minute the record was recorded > >>>> iiii -- instrument number > >>>> c -- channel number > >>>> How would I format the data_structure parameter> > >>>> Thanks, > >>>> Rob > >>>> +----------------------------+ > >>>> Robert E. Abbott, Ph.D. > >>>> Sandia National Laboratories > >>>> Geophysics Department MS 0750 > >>>> P.O. Box 5800 > >>>> Albuquerque, NM 87185-0750 > >>>> (505) 845-0266 > >>>> +----------------------------+ > >>>> _______________________________________________ > >>>> MSNoise mailing list > >>>> MSNoise(a)mailman-as.oma.be > >>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>> _______________________________________________ > >>> MSNoise mailing list > >>> MSNoise(a)mailman-as.oma.be > >>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >> _______________________________________________ > >> MSNoise mailing list > >> MSNoise(a)mailman-as.oma.be > >> http://mailman-as.oma.be/mailman/listinfo/msnoise > > _______________________________________________ > > MSNoise mailing list > > MSNoise(a)mailman-as.oma.be > > http://mailman-as.oma.be/mailman/listinfo/msnoise > > > ------------------------------ > > Message: 3 > Date: Tue, 14 Jan 2014 22:35:50 +0100 > From: Thomas Lecocq <thomas.lecocq(a)oma.be> > To: Python Package for Monitoring Seismic Velocity Changes using > Ambient Seismic Noise <msnoise(a)mailman-as.oma.be> > Subject: Re: [MSNoise] [EXTERNAL] Re: Possible to use data from > ref2mseed? > Message-ID: <a911a3a8a074e094f158064841ac1133(a)oma.be> > Content-Type: text/plain; charset=UTF-8; format=flowed > > Rob, > > Le 14.01.2014 21:40, Abbott, Robert E a ?crit?: > > Thanks, Thomas. I am on travel right now, so I will have to get you > > the parameters Monday. There is a finite chance I screwed up and gave > > you SEG-Y files, they are named similarly. If that is the case, I > > apologize for the confusion. > > I tried reading SEG-Y but it doesn't work either. Let's wait till Monday > to check ! > > Cheers > > Thomas > > > > > Sent from my iPad > > > >> On Jan 14, 2014, at 10:27 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> > >> wrote: > >> > >> Hi Rob, > >> > >> OK, got your files. I can modify the code to automatically find the > >> stations, that's OK. > >> > >> But, I can't read them in Python. Although you state it's MiniSeed > >> (and I believe you :-) ), I can't find the proper parameters to read > >> them. Did you pass any argument to the converter ? Are you able to > >> open them using another tool ? > >> > >> Cheers, > >> > >> Thomas > >> > >> Le 10.01.2014 18:18, Abbott, Robert E a ?crit : > >>> Thomas, > >>> Thanks for very much for your help. I will send you some data, but 3 > >>> days, > >>> 3 stations is quite a large data volume as my sample rate is 200 Hz. > >>> Do > >>> you have a dropbox.com account? That has worked in the past for me. > >>> Alternatively, I could send just an hour or two of data per day... > >>> Thanks, > >>> Rob > >>> +----------------------------+ > >>> Robert E. Abbott, Ph.D. > >>> Sandia National Laboratories > >>> Geophysics Department MS 0750 > >>> P.O. Box 5800 > >>> Albuquerque, NM 87185-0750 > >>> (505) 845-0266 > >>> +----------------------------+ > >>>> On 1/10/14 9:06 AM, "Thomas Lecocq" <thomas.lecocq(a)oma.be> wrote: > >>>> Hi Rob, > >>>> Yes, it is definitively possible. Highly dependant on your python > >>>> skills, I'm afraid... > >>>> The data structure will be something like : > >>>> data_structure['ref2mseed'] = "RDAY.LOC/YEAR.DAY.HOUR.MIN.STA.CHAN" > >>>> but: > >>>> I assume yy is a 2digit year format ? Then YEAR will not work, we > >>>> need > >>>> to put YY and define it in the reader (in populate station table) > >>>> I set STA to your instrument number, which I suppose are located at > >>>> a > >>>> unique position within this very archive. > >>>> I can help writing the reader part (populate station table). I'm not > >>>> quite sure how the scan_archive process will (dis)like such a bad > >>>> structure (no net/ or sta/ folders)... > >>>> If you provide me some data, say 3 days, 3 stations, I can test the > >>>> code > >>>> on it... > >>>> Best regards from Brussels, > >>>> Thomas > >>>> Le 09/01/2014 22:43, Abbott, Robert E a ?crit : > >>>>> Greetings, > >>>>> Is it possible to use MSNoise on data generated from the PASSCAL > >>>>> Instrument Center utility "ref2mseed" (the data is from Reftek 130 > >>>>> digitizers). The data is named like this: > >>>>> Rddd.0S/yy.ddd.hh.mm.ss.iiii.c > >>>>> where > >>>>> yy -- year of the record > >>>>> ddd -- julian day of the record > >>>>> S -- data stream of the instrument (usually 1) > >>>>> hh -- hour the record was recorded > >>>>> mm -- minute the record was recorded > >>>>> iiii -- instrument number > >>>>> c -- channel number > >>>>> How would I format the data_structure parameter> > >>>>> Thanks, > >>>>> Rob > >>>>> +----------------------------+ > >>>>> Robert E. Abbott, Ph.D. > >>>>> Sandia National Laboratories > >>>>> Geophysics Department MS 0750 > >>>>> P.O. Box 5800 > >>>>> Albuquerque, NM 87185-0750 > >>>>> (505) 845-0266 > >>>>> +----------------------------+ > >>>>> _______________________________________________ > >>>>> MSNoise mailing list > >>>>> MSNoise(a)mailman-as.oma.be > >>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>>> _______________________________________________ > >>>> MSNoise mailing list > >>>> MSNoise(a)mailman-as.oma.be > >>>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >>> _______________________________________________ > >>> MSNoise mailing list > >>> MSNoise(a)mailman-as.oma.be > >>> http://mailman-as.oma.be/mailman/listinfo/msnoise > >> _______________________________________________ > >> MSNoise mailing list > >> MSNoise(a)mailman-as.oma.be > >> http://mailman-as.oma.be/mailman/listinfo/msnoise > > _______________________________________________ > > MSNoise mailing list > > MSNoise(a)mailman-as.oma.be > > http://mailman-as.oma.be/mailman/listinfo/msnoise > > > ------------------------------ > > _______________________________________________ > MSNoise mailing list > MSNoise(a)mailman-as.oma.be > http://mailman-as.oma.be/mailman/listinfo/msnoise > > > End of MSNoise Digest, Vol 3, Issue 3 > ************************************* >

12 years, 1 month

Re: [MSNoise] Consult about 07.plot_dtt.py

by Thomas Lecocq

Hi Esteban, Sorry for not replying earlier, your email got fw to my AGU-related folder automatically, dunno why... Anyway, I cc the mailing-list so your problem could help others: Le 10/01/2014 08:24, Esteban Chaves a écrit : > Hi Thomas, > > How it’s going? Happy new year. > > After run successfully the msnoise using the example data you gave us > at the workshop, > I was trying to plot the results using the script 07.plot_dtt.py but I > got the next error: > > File "07.plot_dtt.py", line 14, in <module> > import matplotlib.pyplot as plt > File > "/Users/ech/anaconda/lib/python2.7/site-packages/matplotlib/pyplot.py", line > 98, in <module> > _backend_mod, new_figure_manager, draw_if_interactive, _show = > pylab_setup() > File > "/Users/ech/anaconda/lib/python2.7/site-packages/matplotlib/backends/__init__.py", > line 28, in pylab_setup > globals(),locals(),[backend_name],0) > File > "/Users/ech/anaconda/lib/python2.7/site-packages/matplotlib/backends/backend_wxagg.py", > line 6, in <module> > import backend_wx # already uses wxversion.ensureMinimal('2.8') > File > "/Users/ech/anaconda/lib/python2.7/site-packages/matplotlib/backends/backend_wx.py", > line 53, in <module> > raise ImportError(missingwxversion) > ImportError: Matplotlib backend_wx and backend_wxagg require > wxversion, which was not found. > Well, it looks like the WX backend is missing on your machine (backends are what render your graphic interfaces). So, try to delete line 3 (https://github.com/ROBelgium/MSNoise/blob/master/07.plot_dtt.py#L3) of the file and leave matplotlib to its default config, which, I guess, will use the PySide or Qt4 backends. > > I’ll be very grateful if you can help me to try to figure it out. Have > you created other scripts to plot the results? Yes, by the end of next week, they will be on github too ! Cheers, Thomas > > Thanks in advance, > > Esteban > >

12 years, 1 month

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

MSNoise January 2014