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

Abbott, Robert E reabbot at sandia.gov
Wed Jan 22 17:27:38 UTC 2014


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



More information about the MSNoise mailing list