[MSNoise] [EXTERNAL] Re: Possible to use data from ref2mseed?

Thomas Lecocq thomas.lecocq at oma.be
Wed Jan 22 15:27:08 UTC 2014 

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

More information about the MSNoise mailing list