Анализ данных с электросчетчика Eastron SDM220 средствами ThingSpeak09.09.2016 17:03

apt-get -y install python-pip
apt-get install python2.7-dev

pip install -U pymodbus

def read_float_reg (client, basereg, unit=1):
resp = client.read_input_registers (basereg,2, unit=1)
if resp == None:
return None
# according to spec, each pair of registers returned
# encodes a IEEE754 float where the first register carries
# the most significant 16 bits, the second register carries the
# least significant 16 bits.
return struct.unpack ('>f', struct.pack ('>HH',*resp.registers))

def fmt_or_dummy (regfmt, val):
if val is None:
return '.'*len (regfmt[2]%(0))
return regfmt[2]%(val)

def main ():
regs = [
# Symbol Reg# Format
('V:', 0×00, '%6.2f'), # Voltage [V]
('Curr:', 0×06, '%6.2f'), # Current [A]
('Pact:', 0×0c, '%6.0f'), # Active Power («Wirkleistung») [W]
('Papp:', 0×12, '%6.0f'), # Apparent Power («Scheinl.») [W]
('Prea:', 0×18, '%6.0f'), # Reactive Power («Blindl.») [W]
('PF:', 0×1e, '%6.3f'), # Power Factor [1]
('Phi:', 0×24, '%6.1f'), # cos (Phi)? [1]
('Freq:', 0×46, '%6.2f'), # Line Frequency [Hz]
('Wact:', 0×0156, '%6.2f'), # Energy [kWh]
('Wrea:', 0×0158, '%6.2f'), # Energy react [kvarh]
]

cl = pymodbus.client.sync.ModbusSerialClient ('rtu',
port='/dev/ttyUSB0', baudrate=9600, parity='N', stopbits=1,
timeout=0.8)

values = [ read_float_reg (cl, reg[1], unit=1) for reg in regs ]
outvals = list ((' '.join ([fmt_or_dummy (*t) for t in zip (regs, values)])).split ())
params = urllib.urlencode ({'key': 'xxxxxxxxxxxxxxxx', 'field1': outvals[0], 'field2': outvals[1], 'field3': outvals[2], 'field4': outvals[8]})
f = urllib.urlopen («api.thingspeak.com/update», data=params)
print (outvals)
sys.stdout.flush ()

if __name__ == '__main__':
main ()

% ID канала для записи обработанных данных (PowerStatistic)
writeChannelID = 157182;
writeAPIKey = 'xxxxxxxxxxxxxxxxxx';

%% Read Data %%
data1 = thingSpeakRead (readChannelID, 'Fields', 4);
data2 = thingSpeakRead (readChannelID, 'Fields', 4, 'NumMinutes', 60);
value = data1-data2(1);

%в data1 читаем последнее значение из 4-го поля (Energy)
%чтение в data2 сделано криво, но по другому я пока не придумал, читаются последние 60 значений (за час) и берется первое из них

%% Analyze Data %%
% пришлось добавить округление до двух знаков, так как несмотря на то, что в исходном массиве значения лежат округленные до двух знаков после запятой, результат вычитания почему-то выглядел так: 0.2700000000000031
analyzedData = round (value,2);

%disp можно использовать для отладки данных
%disp (analyzedData);

Если в analyzedData одно значение, а не вектор, оно будет записано в первое поле канала. Для записи в другие поля нужно добавить 'Fields',2, например.
%% Write Data %%
thingSpeakWrite (writeChannelID, analyzedData, 'WriteKey', writeAPIKey);

writeChannelID = 157182;
writeAPIKey = 'xxxxxxxxxxxxxxxxxxxx';

%% Read Data %%
% читаем последние 24 значения из первого поля канала
data = thingSpeakRead (readChannelID, 'ReadKey', readAPIKey, 'Fields',1, 'NumPoints',24);

%% Analyze Data %%
%суммируем
analyzedData = sum (data);

%disp (analyzedData);

%% Write Data %%
thingSpeakWrite (writeChannelID, analyzedData, 'WriteKey', writeAPIKey, 'Fields',2);