def get_coefficients(self): log.log('get_coefficients function called') self.coefficients = {} converter = 8192 if self.mj in ('mj04', 'mj22', 'mj36'): converter = 4096 l = [] f = './temp/ini/base.ini' if not os.path.isfile(f): return i = IniParser() i.read(f) for s in i.options('AnalogInput_Signals'): if 'scale' in s.lower() and i.get('AnalogInput_Signals', s) != 'x': l.append((s, i.getint('AnalogInput_Signals', s))) #endfor f = './temp/ini/parameters.ini' if not os.path.isfile(f): return i = IniParser() i.read(f) for s, a in l: self.coefficients[s] = i.getfloat(s, 'Range') #endfor f = './temp/ini/settings.ini' if not os.path.isfile(f): return i = IniParser() i.read(f) for s, a in l: print s, a print self.coefficients[s] print i.getfloat('AnalogInput_Coefficients', 'Coefficient%i' % a) print type(self.coefficients[s]) print type(i.getfloat('AnalogInput_Coefficients', 'Coefficient%i' % a)) capacity = self.coefficients[s] difference = (i.getfloat('AnalogInput_Coefficients', 'Coefficient%i' % a) / (capacity/converter) - 1) * 100 offset = i.getint('AnalogInput_Offsets', 'Offset%i' % a) self.coefficients[s] = (capacity, difference, offset)
def get_little_amounts(self): log.info('little amounts function called') silo = {} self.little_amounts_fails = {} doba = None vaha = None navazeno = None analog = None adaptivita = None pozadavek = None zasobnik = 99 # can delete this line after three month (bacuse of error in info txt from Piestany) fn = 'control.log' for l in give_lines(fn, 3): if 'poddavkovani ' in l.lower(): zasobnik = int(l.split()[-1]) vaha = None navazeno = None pozadavek = None elif 'dobaotevreni ' in l.lower(): pom = float(l.split()[-1]) if pom > 3 or pom < 0: continue # filtrace demetskych hodnot v milisekundach a zapornych hodnot doba = pom elif 'mezistav1 ' in l.lower(): vaha = float(l.split()[-1]) elif 'NejblizsiCil2' in l and vaha: pozadavek = float(l.split()[-1]) - vaha elif 'PrumernaHodnotaAI[CisloAnalogu] - VychoziVaha' in l: analog = float(l.split()[-1]) elif 'frakcenamereno2 ' in l.lower() and vaha: navazeno = float(l.split()[-1]) - vaha #endif if doba and navazeno: log.info('%i, %0.2f, %i' % (navazeno, doba, zasobnik)) if navazeno < 30 or navazeno > 500: log.info('little amount fail: %i, %0.2f, %i' % (navazeno, doba, zasobnik)) if not zasobnik in self.little_amounts_fails: self.little_amounts_fails[zasobnik] = [] #endif self.little_amounts_fails[zasobnik].append('%i, %0.2f' % (navazeno, doba)) doba = None continue #endif if analog: adaptivita = navazeno - analog + vaha if not zasobnik in silo: silo[zasobnik] = [] silo[zasobnik].append((doba, navazeno, pozadavek, adaptivita)) doba = None vaha = None navazeno = None analog = None adaptivita = None #endif #endfor #log.info('%s' % silo) self.silo = {} self.little_amounts_error = False f = './temp/ini/adapt.ini' if not os.path.isfile(f): return adapt = IniParser() adapt.read(f) f = './temp/ini/parameters.ini' if not os.path.isfile(f): return parameters = IniParser() parameters.read(f) for k in silo: self.little_amounts_error = adapt.has_option('Aggregate', 'VelocityOfFilling1-%s' % k) param_velocity = 0 param_offset = 0 if parameters.has_option('Aggregate_Bin1_Silo%s' % k, 'Velocity_Base'): param_velocity = parameters.getint('Aggregate_Bin1_Silo%s' % k, 'Velocity_Base') param_offset = parameters.getfloat('Aggregate_Bin1_Silo%s' % k, 'Additional_Time') else: self.little_amounts_error = True #endif x = 0 sez = {} while x < 1.5: if not x in sez: sez[x] = [] for doba, navazeno, pozadavek, adaptivita in silo[k]: if doba-x == 0: continue #I dont know why sez[x].append(navazeno/(doba - x)) #endfor x = x + 0.01 #endwhile minimum = None for x in sorted(sez): sqr_sum = 0 avg = 0 dev = 0 avg = sum(sez[x]) * 1.0 / len(sez[x]) for velocity in sez[x]: sqr_sum += ((velocity - avg) ** 2) / len(sez[x]) #endfor dev = sqr_sum ** 0.5 if not minimum or dev < minimum: #print dev minimum = dev self.silo[k] = (x, avg, len(sez[x]), dev, param_velocity, param_offset)
def get_adaptivity(mj): log.log('get_adaptivity function called') class Adaptivity(): def __init__(self, silo, scale, coefficient): self.silo = silo #log.info('creating new class for %s' % self.silo) self.scale = scale self.coefficient = coefficient self.first = None #enddef def check_line(self, line): if not self.scale in line: return if not self.first and self.first != 0: self.first = int(line.split('=')[1]) #log.info('%s first %s' % (self.silo, self.first)) else: self.adaptivity = (int(line.split('=')[1]) - self.first) * self.coefficient self.time = get_datetime_from_log(line) return 'all done' #endif #enddef #enclass ret = {} # ret[silo] = [(time, adaptivity), ...] materials = {'AGG': 'Aggregate', 'CEM': 'Cement', 'ADM': 'Admixture', 'WTR': 'Water'} cache = {} # cache[silo] = {'Scale': scale, 'Coefficient': coefficient} incompletes = [] # incompletes = [Adativity(), ...] f = './temp/ini/base.ini' if not os.path.isfile(f): return ret base = IniParser() base.read(f) f = './temp/ini/settings.ini' if not os.path.isfile(f): return ret settings = IniParser() settings.read(f) fn = 'recordedsignals.log' for l in give_lines(fn, 3): for m in materials: if m in l and 'O_' in l and 'intoBin' in l and ')=0' in l: silo = '%s%s' % (m, l.split(m)[1].split('into')[0]) # for example AGG5 if not silo in cache: bin = 'Bin%s' % l.split('intoBin')[1].split(')=')[0] scale = '%s_Scale%s' % (materials[m], base.get('%s_%s' % (materials[m], bin), 'CorrespondingScale')) signal = base.get('AnalogInput_Signals', scale) coefficient = settings.getfloat('AnalogInput_Coefficients', 'Coefficient%s' % signal) cache[silo] = {} cache[silo]['Scale'] = scale cache[silo]['Coefficient'] = coefficient #endif if not silo in ret: ret[silo] = [] # Adding new adaptivities to list incompletes.append(Adaptivity(silo, cache[silo]['Scale'], cache[silo]['Coefficient'])) break #endif #endfor if 'Scale' in l: for i in incompletes: if i.check_line(l): ret[i.silo].append((i.time, i.adaptivity)) incompletes.remove(i) #endif #endfor #endif #endfor # read alternately adaptivity from info.txt fn = 'control.log' materials = {'aggregate': 'AGG', 'cement': 'CEM', 'admixture': 'ADM', 'water': 'WTR'} for l in give_lines(fn, 3): if not 'adaptivity for' in l: continue material_type = materials[l.split(' for ')[1].split(' bin ')[0]] silo_number = l.split(' silo ')[1].split(':')[0] adaptivity = l.split(': ')[-1].strip() time = get_datetime_from_log(l) silo = '%s%s' % (material_type, silo_number) if not silo in ret: ret[silo] = [] ret[silo].append((time, float(adaptivity))) #endfor # create output files for silo in ret: log.info('creating file for silo %s' % silo) a = '../public_html_incoming' if not os.path.isdir(a): a = './public_html_incoming' #endif f = open('%s/%s%s.csv' % (a, mj, silo), 'w') for x,y in ret[silo]: f.write('%s;%i\n' % (x.strftime('%Y-%m-%d %H:%M:%S'), y)) #endfor f.close() #endfor return ret