³ò dÄGJc@sÔdZddkZddkZddkZddkZddkZddkZddkZde fd„ƒYZ ddk l Z ddk lZddklZddklZdd klZdd klZdS( s> Generic routines for manipulating log files. 5 Dec 2008, Lev iÿÿÿÿNtTimeLogcBsøeZdZddd„Zd„Zd„Zd„Zd„Zd„Z d„Z d„Z d „Z d „Z d „Zd „Zddd „Zd„Zdd„Zd„Zd„Zd„Zd„Zdd„Zedded„ƒZdd„ZRS(s_ A base class for various logs. Data is stored in a 2-dimensional array 'self.data': First column ([:,0]) is time (sec, since 00:00 1 January 1970 GMT), the rest is log-specific 'self.device' contains the name of the logged device. The purpose is to separate logs of the same type, e.g. paroscientific and melting curve thermometry logs. Two logs can be concatenated with a '&' sign, one of the logs can be None - in this case the other is returned as a result of concatenation. The child classes should provide the following functionality: 1. Methods to access custom 'self.data' fields. 2. A 'fieldcount()' static method that returns number of columns in 'self.data' array. 3. A 'load(filename)' static method that should load a log from a file and return an object of the appropriate child class. cCs‰|djotid|iƒfƒ}nti|ƒ}t|iƒdjp|id|iƒjotSn||_||_ dS(s„ TimeLog constructor. 'log' is an N x fieldcount() array, containing the log. If the argument is None, the empty log is created. 'device' is the name of the device that was used for the log. Columns have the following meaning: 1 [:,0] Time (sec, since 00:00 1 January 1970 GMT) ... rest is log-specific iiiN( tNonetnumpytzerost fieldcounttarraytlentshapet ValueErrortdatatdevice(tselftlogR ((sc:\py\lib\timelog.pyt__init__(s 0 cCs t|iƒS(N(RR (R ((sc:\py\lib\timelog.pyt__len__>scCs|idd…dfS(s,time, seconds since 00:00 1 January 1970 GMTNi(R (R ((sc:\py\lib\timelog.pytt@scCstii|iƒƒS(s!time for pylab.plot_date function(t matplotlibtdatest epoch2numR(R ((sc:\py\lib\timelog.pytplottimeAscCs1titig|iƒd|iƒd fƒS(sn Return intervals between pulses, first element is NaN - the delay before the first pulse, so iiÿÿÿÿ(Rt concatenatetNaNR(R ((sc:\py\lib\timelog.pyt delay_beforeCscCst|ƒdjS(s?return True if there are no entries in the log, False otherwisei(R(R ((sc:\py\lib\timelog.pytisemptyIscCst|ƒ|iiƒ|iƒS(N(ttypeR tcopyR (R ((sc:\py\lib\timelog.pyRKscCsJt|ƒt|ƒjo1|iid|iidjo|i|ijS(sjReturn True if the two TimeLogs are compatible (of the same type and for the same device), otherwise Falsei(RR RR (R tother((sc:\py\lib\timelog.pyt compatibleMscCsZ|djo|Sn|i|ƒp t‚nt|ƒti|i|ifƒ|iƒS(N(RRRRRRR R (R R((sc:\py\lib\timelog.pyt__and__Ss   cCs|djo|Snt‚dS(N(RR(R R((sc:\py\lib\timelog.pyt__rand__^s cCsN|iƒodSn|i|idd…dfiƒdd…f|_dS(s7 Sort the log accending the timestamps Ni(RR targsort(R ((sc:\py\lib\timelog.pytsortds cCsLt|ƒt|ƒdjo|i|ntid|iƒfƒ|iƒS(sV Return a subset of the log consisting of entries listed in 'indices' i(RRR RRRR (R tindices((sc:\py\lib\timelog.pytsubsetlscCs¢|iƒo|Sn|djot|iƒƒ}nti|ƒ}|djot|iƒƒ}nti|ƒ}|i|iƒ|j|iƒ|j@ƒS(sô Return a subset of the log between the two dates, in C and UNIX format or free-form strings. If the start date is not specified, the earliest is assumed. If the end date is not specified, the latest is assumed. N(RRtminRtflibttmtmaxR!(R tstarttend((sc:\py\lib\timelog.pytperiodrs   cCsP|iƒo|Sn|i|iƒti|ƒj|iƒti|ƒjBƒS(su Return a subset of the log outside the specified period, in C and UNIX format or free-form strings. (RR!RR#R$(R R&R'((sc:\py\lib\timelog.pytnoperiodˆs cCs»ti|ƒ}|dj oti|ƒn|}|iƒ}t||jƒoti|||jƒ}nt||jƒoti|||jƒ}n|i||j||j@ƒS(s¶ Return a subset of the log encompassing a period of time defined by 'start'-'end' or a moment in time defined by 'start', if 'end' is not specified. N( R#R$RRtanyRR%R"R!(R R&R'R((sc:\py\lib\timelog.pytaround’s# cCsX|idjo!dt|ƒi|iƒfSn$d|it|ƒi|iƒfSdS(Ns%s %ss%s %s %s(R RRt__name__tdatespan(R ((sc:\py\lib\timelog.pyt__repr__£s!cCsA|iƒodSn)tit|iƒƒt|iƒƒƒSdS(sK Return string representation of the date range of the log s(empty)N(RR#R-R"RR%(R ((sc:\py\lib\timelog.pyR-©s cCs8|iƒ}|i|t||ti|ƒjƒjƒS(s¬ Return a single point log with the regord at the moment 't', or just before it. Time is accepted in any format recognised by 'flib.tm()'. (RR!R%R#R$(R Rttt((sc:\py\lib\timelog.pyt justbefore²s cCs8|iƒ}|i|t||ti|ƒjƒjƒS(s¬ Return a single point log with the regord at the moment 't', or just before it. Time is accepted in any format recognised by 'flib.tm()'. (RR!R"R#R$(R RR/((sc:\py\lib\timelog.pyt justafter»s c s~ˆdjo d„}n‡fd†}t|ƒt|ƒjotdƒ‚ntit|ƒgdtƒ}|iƒ}xûtt|ƒƒD]ç}xct|dddƒD]K} || |||joPn||| ||ƒpt||Òscst||ƒˆjS((tabs(R2R3(tallow_change_by(sc:\py\lib\timelog.pyR4Ôss"condition" is wrong sizetdtypeiiÿÿÿÿN( RRRRtonestboolRtxrangetFalse( R t conditiontbeforetafterR6tsametmaskRtitj((R6sc:\py\lib\timelog.pytsteadyÄs2       c Csh|dj otiti|ƒƒ}nd}|dj otiti|ƒƒ}nd}d}xÒti|ƒD]Á}d} t|ƒdjo‡|d |jov|d |joe|o<y"|itii ||ƒƒ} Wqt j oqXq|itii ||ƒƒ} n| dj o|| M}qtqtW|djo |ƒSn|i ƒ|i |ƒS(s4 Load log files from a given directory. The function walkes through a directory and treats files named 'YYYYMMDD.dat' as log files. If either boundary time (start, end) is specified, only files with relevant dates are read. Selected files are read using the load() static method of a child class and appended to the log using '&'. If no files were read, an empty log is created using the child class default constructor. The resulting log is returned. t00000000t99999999iN( RR#t dt2filenameR$tostlistdirRtloadtpathtjoint ExceptionRR(( tclstdirnameR&R't skipErrorststartfntendfnR tfilenametday((sc:\py\lib\timelog.pytloaddirðs.  5" "    i cCs(|iƒ}ti|i|ƒ|_|S(se Average 'Nwindow' point long sections of the log to reduce size and scatter (RR#taverageR (R tNwindowtl((sc:\py\lib\timelog.pyRU!s N(R,t __module__t__doc__RR RRRRRRRRRRR!R(R)R+R.R-R0R1RCt classmethodR;RTRU(((sc:\py\lib\timelog.pyRs0             ,0(tPLMLog(tMCTLog(tParoLog(tPIDLog(t KelvinoxLog(t OneKPotLog(RYRRGtos.pathttimetdatetimetmatplotlib.datesRR#tobjectRtplmlogR[tmctlogR\tparologR]tpidlogR^t kelvinoxlogR_t oneKpotlogR`(((sc:\py\lib\timelog.pyss       ÿ