[teqc] helpful tip of week 1921

Lou Estey lou at unavco.org
Wed Nov 2 08:08:15 MDT 2016


This week's tip: using time filtering options in combination

You may have noticed that quite a few tips so far have involved time filtering options.
Let's review (links at the end of this email):

- setting GPS week with '-week' (tip of week 1904)
- ordinary input file splicing (tip of week 1914)
- explicit global start, delta, and/or end time windowing (tip of week 1919)
- GNSS observation decimation (tip of week 1915)
- met observation decimation (tip of week 1920)
- '-tbin' option for "time binning" (tip of week 1917)

This week's tip is a reminder that you can use these in essentially any combination
(at least with the beta version and the next official release, which will probably
be next week; caveats for past versions of teqc are at the end of this email).

As an example, let's suppose you have some 10-Hz hourly raw data files, starting at
an irregular time on Oct 18, but you need daily RINEX decimated to 15 seconds and want
the output to start on Oct 19.  And the first raw data file is messed up in some way
so that you need to tell teqc the starting GPS week and that is week 1920.  And let's
also suppose the raw data contain met observations with met sampling at 5 seconds with
a 1 second offset from the nominal times, and that you want daily RINEX met as well
but at 1 minute intervals.  (And let's suppose that teqc's auto-identification of the
raw data files is working to simply that part of the command.)  You might then have
a command that looks something like:

teqc -week 1920 -O.dec 15 -O.int 0.1 -M.dec 1m:1 -st 2016-10-19T00:00:00 +obs + +met + -tbin 1d site file*

here using UNXI/Linux filename globbing with 'file*' (see tip of week 1901) where file1,
file2, etc. are your raw data files listed in time order.  The end result from teqc
should be daily RINEX obs and met files with the first file starting on 2016 Oct 19
(2016-293), all output files with 'site' as the file prefix.

Using '-tbin' while doing qc is not possible, but you can use the other time filtering
options during qc.  For example, let's suppose you you want to qc the observation data
at 30-second sampling (say, to compare with some other 30-second site), and you want to
quickly see 2 days worth of qc staring on the 19th:

teqc +qc -week 1920 -O.dec 30 -O.int 0.1 -st 2016-10-19T00:00:00 +dh 48 file*

Remember from tip of week 1919 that the global time window options specify limits that
are inclusive, so the last epoch using the above command will be 2016-10-21T00:00:00 --
which means you'll probably have 5761 epochs in the qc -- so if you want to make sure
the qc only goes to 2016-10-20T23:59:30 with 5760 epochs, use:

teqc +qc -week 1920 -O.dec 30 -O.int 0.1 -st 2016-10-19T00:00:00 -e 2016-10-20T23:59:30 file*
equivalent to:
teqc +qc -week 1920 -O.dec 30 -O.int 0.1 -st 2016-10-19T00:00:00 +dm 2879.5 file*

A couple of notes on this:

Yes, you can do a qc on raw data, as long as teqc knows what the window length will be,
which it is getting from the global time windowing options in the above commands.

And because you are decimating to 30 seconds, you don't have to worry about adjustments
for 10-Hz high-rate qc; see tip of week 1895: http://postal.unavco.org/pipermail/teqc/2016/002071.html

---------

Caveats with past versions of teqc:

- depending on the version, the global time windowing option may not work with '-tbin'
when creating the RINEX observation files; this was fixed on 17 Oct 2016

- the optional offset for met observations with '-M.dec' was included on 26 Sept 2016

So both of these fixes are in the 2016Oct18 beta (see
http://postal.unavco.org/pipermail/teqc/2016/002160.html ) and will be in the next
official release.

Happy teqc-ing!

cheers,
--lou

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Louis H. Estey, Ph.D.              office:  [+001] 303-381-7456
UNAVCO, 6350 Nautilus Drive           FAX:  [+001] 303-381-7451
Boulder, CO  80301-5554            e-mail:  lou  unavco.org
      WWW:  http://www.unavco.org   http://jules.unavco.org

"If the universe is the answer, what is the question?"
                                                -- Leon Lederman
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

Past helpful tips:

week 1894: using teqc config files - http://postal.unavco.org/pipermail/teqc/2016/002067.html
week 1895: qc of high-rate data - http://postal.unavco.org/pipermail/teqc/2016/002071.html
week 1896: UNIX/Linux shells for Windows - http://postal.unavco.org/pipermail/teqc/2016/002072.html
week 1897: '-' vs. '+' teqc options - http://postal.unavco.org/pipermail/teqc/2016/002076.html
week 1898: auto-identification of formats - http://postal.unavco.org/pipermail/teqc/2016/002092.html
week 1899: auto-identification vs. format flags - http://postal.unavco.org/pipermail/teqc/2016/002096.html
week 1900: square brackets in options - http://postal.unavco.org/pipermail/teqc/2016/002105.html
week 1901: using option '+mds' - http://postal.unavco.org/pipermail/teqc/2016/002108.html
week 1902: qc results w/ problematic nav messages - http://postal.unavco.org/pipermail/teqc/2016/002113.html
week 1903: '-no_orb[it]' and '-no_pos[ition]' options - http://postal.unavco.org/pipermail/teqc/2016/002115.html
week 1904: '-week' option - http://postal.unavco.org/pipermail/teqc/2016/002117.html
week 1905: using '+bcf' for XYZ/geodetic conversion - http://postal.unavco.org/pipermail/teqc/2016/002126.html
week 1906: the '+v[erify]' option - http://postal.unavco.org/pipermail/teqc/2016/002128.html
week 1907: '+C2', '+L5', "+L6', '+L7', '+L8', and '+all' options - http://postal.unavco.org/pipermail/teqc/2016/002130.html
week 1908: no doppler shortcut; RINEX L2 - http://postal.unavco.org/pipermail/teqc/2016/002131.html
week 1909: using paths w/ file names - http://postal.unavco.org/pipermail/teqc/2016/002132.html
week 1910: the (un)importance of file names - http://postal.unavco.org/pipermail/teqc/2016/002133.html
week 1911: notices, warnings, and errors - http://postal.unavco.org/pipermail/teqc/2016/002134.html
week 1912: the '-max_rx_SVs' option - http://postal.unavco.org/pipermail/teqc/2016/002137.html
week 1913: the end of '++igs' and '+igs' - http://postal.unavco.org/pipermail/teqc/2016/002140.html
week 1914: splicing together RINEX files - http://postal.unavco.org/pipermail/teqc/2016/002144.html
week 1915: using '-O.int' and '-O.dec' - http://postal.unavco.org/pipermail/teqc/2016/002145.html
week 1916: '+doy' option - http://postal.unavco.org/pipermail/teqc/2016/002146.html
week 1917: '-tbin' and '-ast' options - http://postal.unavco.org/pipermail/teqc/2016/002152.html
week 1918: mp12 RMS before/after Oct 2013 - http://postal.unavco.org/pipermail/teqc/2016/002158.html
week 1919: the global windowing options - http://postal.unavco.org/pipermail/teqc/2016/002159.html
week 1920: '-M.dec' and '-N.dec' options - http://postal.unavco.org/pipermail/teqc/2016/002163.html




More information about the teqc mailing list