Gráfica de lecturas de un DHT11 en Python

Esta nota corta será para no olvidar como hacer una gráfica en Python de un conjunto de mediciones de temperatura a lo largo del tiempo. El sensor utilizado en un DHT11 (que muy tardé me di cuenta que tenía muy baja resolución y que es mejor el DHT22). Generé un archivo de texto a partir de comunicación serial entre un PIC y una Raspberry Pi dónde cada linea tiene la siguiente estructura:

2019-03-22 23:11:33,22,46

fecha, temperatura, humedad relativa 

Algo muy útilil que aprendí a hacer es trasferir archivos entre un host local y la Raspberry a través de la red con un solo comando de consola. El comando es scp y su sintaxis es la siguiente para copiar del host local a la Raspberry :

$ scp /un/dir/local/archivo usuario@ipremota:/un/dir/remoto/archivo

Y de la Raspberri al host local:

$ scp usuario@ipremota:/un/dir/remoto/archivo /un/dir/local/archivo

Ya con el archivo en mi computadora, utilicé el módulo Pandas para leer el archivo .txt dónde registre cerca de dos dias y medio de datos de temperatura y humedad relativa. Aplique un filtro de mediana con el modulo scipy.signal para quitar los picos abruptos en los datos. Finalmente realice los ajustes necesarios para mostrar únicamente la hora en el eje x. El programa completo es el siguiente:

	#!/usr/bin/env python3
	# -*- coding: utf-8 -*-
	"""
	Created on Mon Mar 25 13:22:51 2019
	@author: rodolfo
	"""
	import pandas as pd
	import matplotlib.pyplot as plt
	import matplotlib.dates as mdates
	from scipy import signal
	DATOS = pd.read_csv('file.txt',sep=',',header = None, parse_dates = [0])
	fecha = DATOS[0]
	Temp = DATOS[1]
	Temp_med = signal.medfilt(DATOS[1],kernel_size=11)
	RH = DATOS[2]
	plt.plot(fecha, Temp_med)
	plt.axvline(pd.to_datetime("2019-03-24 00:00:00"),color='g',linestyle='--')
	plt.axvline(pd.to_datetime("2019-03-25 00:00:00"),color='g',linestyle='--')
	plt.grid(True)
	plt.title('Lectura DHT11')
	plt.xlabel('Tiempo')
	plt.ylabel('Temperatura (°C)')
	# Etiquetas de hora en eje x
	plt.gcf().autofmt_xdate()
	myFmt = mdates.DateFormatter('%H:%M')
	plt.gca().xaxis.set_major_formatter(myFmt)

view raw DHT11_plots.py hosted with ❤ by GitHub

Gráfica: