[MSNoise] [EXTERNAL] Re: MSNoise Digest, Vol 3, Issue 7
Thomas Lecocq
thomas.lecocq at oma.be
Thu Jan 23 09:37:15 UTC 2014
Hi Esteban & Rob,
Yeah ! Glad to see users helping other users, that's great !
Rob, let us know how it goes !
Cheers,
Thomas
Le 22/01/2014 18:27, Abbott, Robert E a écrit :
> Thanks Esteban!
>
> I was able to write BUD-style files with the following command no problem
> following 'rt2ms':
>
> msmod *.m -A ./%n/%s/%s.%n.%l.%c.%Y.%j
>
> I guess that solves my problem. So no new 'rt2ms' style need to be written
> as this is trivial to implement.
>
>
> +----------------------------+
> 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 9:58 AM, "Esteban Chaves" <echfisica at gmail.com> wrote:
>
>> 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_wor
>> ds=&all_words=1&titles_only=&projects=1&commit=Submit
>>
>> 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 at mailman-as.oma.be>
>> wrote:
>>
>>> Send MSNoise mailing list submissions to
>>> msnoise at 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 at mailman-as.oma.be
>>>
>>> You can reach the person managing the list at
>>> msnoise-owner at 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 at sandia.gov>
>>> To: "Python Package for Monitoring Seismic Velocity Changes using
>>> Ambient Seismic Noise" <msnoise at mailman-as.oma.be>
>>> Subject: Re: [MSNoise] [EXTERNAL] Re: Possible to use data from
>>> ref2mseed?
>>> Message-ID: <CF052D8D.4BDC%reabbot at 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 at 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 at 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 at 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 at 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 at mailman-as.oma.be
>>>>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>>>>> _______________________________________________
>>>>>>>>>> MSNoise mailing list
>>>>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>>>> _______________________________________________
>>>>>>>>> MSNoise mailing list
>>>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>>> _______________________________________________
>>>>>>>> MSNoise mailing list
>>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>> _______________________________________________
>>>>>>> MSNoise mailing list
>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>> _______________________________________________
>>>>>> MSNoise mailing list
>>>>>> MSNoise at mailman-as.oma.be
>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>> _______________________________________________
>>>>> MSNoise mailing list
>>>>> MSNoise at mailman-as.oma.be
>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>> _______________________________________________
>>>> MSNoise mailing list
>>>> MSNoise at 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 at oma.be>
>>> To: msnoise at mailman-as.oma.be
>>> Subject: Re: [MSNoise] [EXTERNAL] Re: Possible to use data from
>>> ref2mseed?
>>> Message-ID: <52DFE34C.8020206 at 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 at 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 at 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 at 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 at 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 at mailman-as.oma.be
>>>>>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>>>>>> _______________________________________________
>>>>>>>>>>> MSNoise mailing list
>>>>>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>>>>> _______________________________________________
>>>>>>>>>> MSNoise mailing list
>>>>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>>>> _______________________________________________
>>>>>>>>> MSNoise mailing list
>>>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>>> _______________________________________________
>>>>>>>> MSNoise mailing list
>>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>> _______________________________________________
>>>>>>> MSNoise mailing list
>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>> _______________________________________________
>>>>>> MSNoise mailing list
>>>>>> MSNoise at mailman-as.oma.be
>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>> _______________________________________________
>>>>> MSNoise mailing list
>>>>> MSNoise at mailman-as.oma.be
>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>> _______________________________________________
>>>> MSNoise mailing list
>>>> MSNoise at 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 at sandia.gov>
>>> To: "Python Package for Monitoring Seismic Velocity Changes using
>>> Ambient Seismic Noise" <msnoise at mailman-as.oma.be>
>>> Subject: Re: [MSNoise] [EXTERNAL] Re: Possible to use data from
>>> ref2mseed?
>>> Message-ID: <CF053777.4BF8%reabbot at 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 at 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;samplerat
>>>>> e;
>>>>> 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 at 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 at 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 at 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 at 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 at mailman-as.oma.be
>>>>>>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>>>>>>> _______________________________________________
>>>>>>>>>>>> MSNoise mailing list
>>>>>>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>>>>>> _______________________________________________
>>>>>>>>>>> MSNoise mailing list
>>>>>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>>>>> _______________________________________________
>>>>>>>>>> MSNoise mailing list
>>>>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>>>> _______________________________________________
>>>>>>>>> MSNoise mailing list
>>>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>>> _______________________________________________
>>>>>>>> MSNoise mailing list
>>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>> _______________________________________________
>>>>>>> MSNoise mailing list
>>>>>>> MSNoise at mailman-as.oma.be
>>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>> _______________________________________________
>>>>>> MSNoise mailing list
>>>>>> MSNoise at mailman-as.oma.be
>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>> _______________________________________________
>>>>> MSNoise mailing list
>>>>> MSNoise at mailman-as.oma.be
>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>> _______________________________________________
>>>> MSNoise mailing list
>>>> MSNoise at mailman-as.oma.be
>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>
>>>
>>> ------------------------------
>>>
>>> _______________________________________________
>>> MSNoise mailing list
>>> MSNoise at mailman-as.oma.be
>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>
>>>
>>> End of MSNoise Digest, Vol 3, Issue 7
>>> *************************************
>>>
>> _______________________________________________
>> MSNoise mailing list
>> MSNoise at mailman-as.oma.be
>> http://mailman-as.oma.be/mailman/listinfo/msnoise
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise
More information about the MSNoise
mailing list