3. Procesa datos. Modelo de pérdidas en propagación LoRa

Para el modelo de perdidas de propagación, en cada punto se registra en el archivo las mediciones de Rssi y SNR. Cada archivo de datos procesan, tabulando y ordenanto los valores representativos del comportamiento del RSSI y SNR para revisar sus descriptores de estadística.

Para realizar el procesamiento de los datos, se crearon algunas funciones y procedimientos para simplificar la escritura de instrucciones, las que se resumen en el archivo girni_lora_libreria.

Las coordenadas geográficas de dada punto se registraron con un GPS diferencial usando el formato UTM en un archivo tipo texto.

El procesamiento de los datos ser realiza en varios pasos donde se revisan los resultados parciales.

El primero de ellos consiste en tabular los datos de Rssi y SNR de cada punto en un solo archivo, luego se añaden las coordenadas y distancias cada punto medido, para finalmente integrar ambos resultados en un solo archivo con RSSI, distancias, coordenadas de cada punto.

Cada sección permite disponer de archivos intermedios que pueden ser usados para observar y procesar resultados que permitan realizar observaciones y mejoras a los modelos planteados. Entre los pasos intemedios está por ejemplo: observar en gráficas las ubicaciones de los puntos en el plano XY usando sus coordenadas, o en otro caso observar los valores de Rssi distribuidos en el espacio formado por el plano del ejemplo anterior y en el eje Z los valores promedios RSSI.

Las siguientes secciones describen lo realizado para:

Procesa datos. Resumen de archivos.txt por punto, Rssi o SNR

Funciones girni_lora_libreria

Procesa datos. Coordenadas GPS y distancias al archivo.txt

Procesa datos. Ubica los puntos en gráfica 2D o 3D

Procesa datos. Integra tablas Rssi, distancias y coordenadas

Procesa Datos. Un punto – revisa descriptores estadísticos

2. Captura datos. Modelo de pérdidas en propagación

El esquema prototipo para el modelo de pérdidas de señal en propagación, se registra los niveles de señal Rssi a diferentes distancias, condiciones y escenarios.

El esquema básico consta de tres balizas y un dispositivo de captura de datos. Cada baliza emite una señal a intervalos de tiempo que al ser captadas por el dispositivo permite obtener los valores Rssi y SNR. Dado que el objetivo final es revisar un modelo de localización, se requieren al menos tres balizas en el esquema.

El modelo de propagación require realizar mediciones en varios puntos a diferentes distancias que en campus presenta diferentes entornos.

Las mediciones de Rssi en un mismo punto no son constantes, por lo que se registran «muchos» valores, al menos 100 lecturas para estimar el comportamiento y usar valores representativos como el promedio y la desviación estándar..

La captura de los datos  se simplifica usando un dispositivo LoRa y un computador portátil que almacena los datos en un archivo txt.

El dispositivo LoRa se programa para tomar lecturas de Rssi y SNR de cada baliza e inmediatamente enviarlos por la conexión serial/USB, al mismo tiempo que el dispositivo se alimenta de energía por el puerto USB.

Desde un computador portátil se leen los datos capturados por el dispositivo y leídos por la conexión USB  y se almacenan en un archivo.txt. La lectura de datos serial/USB se realiza con un algoritmo en Python, que secuencialmente añade datos al archivo.txt para posterior procesamiento.

Las siguientes secciones describen lo realizado para:

Captura datos. Balizas y dispositivo – configuración algoritmo.ino

Captura datos. USB-Serial a un archivo.txt con Python

3.1 Procesa datos. Resumen de archivos.txt por punto, Rssi o SNR

Los archivos de las capturas de datos por USB-Serial que se obtienen en las campañas de medición en cada punto, se requieren: procesar, ordenar y obtener los valores representativos del comportamiento del RSSI y SNR.

Cada archivo contiene al menos 100 lecturas a cada baliza. Los archivos se encuentran en una carpeta desde donde se procesarán y se obtendrá un archivo «resumen».

Para procesar los archivos de deben indicar los siguientes parámetros:

  • medida a observar: rssi, snr
  • descriptor estadistico a usar: mean, std, count, max, min, …
  • trama: rx
  • una ‘carpeta’ donde se encuentran los archivos a procesar
  • nombre del archivo de resultados

Un ejemplo de los resultados en el resumen por punto del promedio Rssi es:

resumen trama: rx  medida: rssi
         rssi_rx_d1  rssi_rx_d2  rssi_rx_d3  rssi_tx_d1  rssi_tx_d2  rssi_tx_d3
FCNM101 -124.396907 -129.132184  -81.620000  -98.000000 -126.632184  -78.040000
FCNM102 -114.345455 -118.917526  -92.053571  -91.000000 -118.474227  -80.482143
FCNM103 -115.459854 -124.634615  -95.453901  -91.000000 -123.173077  -76.617021
FCNM104 -115.922414 -123.678899 -110.090909  -91.000000 -122.844037 -107.933884
FCNM105 -108.239726 -124.893204 -102.007042  -91.000000 -123.582524  -99.281690
FCNM106 -128.383929 -127.910000 -117.589041  -91.000000 -127.290000 -110.082192
FCNM107 -124.429907 -123.009009 -105.078125  -91.000000 -123.225225 -102.343750

Carpeta con archivos a procesar: lectura_puntos

Archivos resultado:

para Rssi: resumen_rssimean01.txt

Para SNR: resumen_snrmean01.txt


Algortimo en Python

Para simplificar el algoritmo, se realizó una libreria de funciones, «girni_lora_libreria.py»,  que se debe de ubicar en la misma carpeta del algoritmo que lo usa.

Para almacenar los datos tabulados, se crea tambien el archivo de ésta tabla en formato .json, en caso de que se requiera realizar otro tipo de análisis a los datos.

'''
Procesa Archivo.txt de LoRa Rssi y SNR
de muestras desde un puerto Serial
Girni 2020-10-07 propuesta: edelros@espol.edu.ec
medida: 'rssi', 'snr'
unestadistico: 'count', 'mean', 'std',
      'min', '25%', '50%', '75%', 'max'
modo: rx baliza a dispositivo
      tx difusion hacia balizas
'''

import os
import numpy as np
import pandas as pd
import girni_lora_libreria as girni

# INGRESO
# revisar parametros al inicio
medida     = 'rssi'
descriptor = 'mean'
trama      = 'rx'

# carpeta de Archivos entrada
carpeta = 'Lecturas_dispositivo'

# Archivos para resultados
arch_rsmUnestadistico = 'resumen_'+medida+descriptor+'02.txt'
arch_detalle = 'tabla_puntosdatos2.json'

# PROCEDIMIENTO
# lista los archivos entrada en la carpeta
archivosPuntos = os.listdir(carpeta)

# Lee cada archivo de un punto y tabula
# las lecturas por cada trama: rx, tx
# cada remitente de paquete: baliza
tabula = {}
for unarchivo in archivosPuntos:
    unpunto  = girni.tabulaPunto(unarchivo,carpeta)
    p_nombre = unpunto['nombre']
    tabula[p_nombre] = unpunto
tabula = pd.DataFrame(tabula)
tabula = tabula.drop(['nombre'])
tabula = tabula.T

# Estadistica descriptiva de cada punto tabulado
descrito = {}
for cadapunto in tabula.index:
    punto = tabula.loc[cadapunto]
    descrito[cadapunto] = girni.describePunto(punto)

# resumen de una medida: rssi, snr
rsm_medida = pd.DataFrame()
for cadapunto in descrito.keys():
    punto = descrito[cadapunto][trama]
    unafila = girni.resumen_medida(punto,medida,descriptor)
    rsm_medida[cadapunto] = unafila
rsm_medida = rsm_medida.T
rsm_medida = rsm_medida.sort_index(axis=1)

# SALIDA
# muestra en pantalla solo un ejemplo.head()
print('resumen trama: ', trama,', medida: ', medida)
print(rsm_medida.head(n=10))

# Se escribe todo el archivo a csv
rsm_medida.to_csv(arch_rsmUnestadistico)

# Para grabar en formato .json
tabula.to_json(arch_detalle)

2.2 Captura datos. Dispositivo con USB-Serial a un archivo.txt usando Python

La relación de valores de Rssi y distancia hacia cada punto de referencia (baliza), se estima a partir de las mediciones en varios puntos establecidos y ubicados en un mapa.

En cada punto, los valores se estiman considerando usar al menos 100 lecturas hacia cada baliza en cada punto.

En éste proceso usa un computador portátil que captura las tramas de cada paquete recibido por el dispositivo desde el puerto USB-Serial a un archivo.txt.

Cada baliza periodicamente emite un paquete que permite al dispositivo registar los valores de Rssi y SNR. Cada captura se registra en una línea de texto como un paquete recibido y etiquetado como «rx».

Para registrar el Rssi y SNR desde el otro extremo, la baliza, se emite un paquete de difusión desde el dispositivo, con el que la baliza al recibirlo toma los valores de Rssi y SNR y los envia al dispositivo en el próximo envío. El paquete de difusión del dispositivo se etiqueta como ‘tx’ en el archivo txt.

En la pantalla del computador portátil se muesta el estado del registro de datos en cada evento de recepción.

Configuración del algoritmo en Python

La captura de datos desde USB se realiza usando instrucciones en Python, con la que se pre-procesa cada paquete y se almacena en un archivo.

En el bloque de INGRESO es necesario indicar:

  • el puerto ‘com’ y los baudios del dispositivo conectado por USB
  • el nombre del archivo de texto que registra la captura de datos
  • las direcciones de las balizas, identificadores indicados en la sección anterior
  • la cantidad mínima de lecturas desde cada baliza

El puerto ‘com’ se puede verificar antes de iniciar la captura, usando el IDE de arduino, al conectar lo muestra en el menu de herramientas/puerto.

El proceso de captura se ejecuta en un lazo, en caso de ser necesario deterner se pueden usar las teclas «Ctrl-c».

En cada punto de referencia se requiere asignar un nombre a cada archivo, iniciando el nombre con «multipunto» añadiendo letras para identificar el sector y un número como identificador del punto.

En los puntos de prueba el nombre del archivo inicia con «mpcircuito». Observe que ambas palabras de inicio contienen la misma cantidad de caracteres, pues facilita el procesamiento de los archivos.

Los archivos se almacenan en una carpeta de donde se puedan leer por lotes en el siguiente proceso de los datos.


Algoritmo en Python

''' Rssi y SNR LoRa punto a punto
paquete de datos desde puerto Serial
generados desde dispositivo LoRa
modelo Heltec Lora 32 v.2
Girni 2020-10-07 propuesta: edelros@espol.edu.ec
'''
import numpy as np
import matplotlib.pyplot as plt
import serial, time

# INGRESO
# Puerto de captura de datos USB-Serial
puerto = 'com3'
baudios = 115200

# Archivo para el registro de cada evento rx,tx
nombrearchivo = 'multipuntoFIEC101.txt'

# identificadores de balizas
baliza = ['d1','d2','d3']

# Mínimo de lecturas por baliza
lecturasMin = 100

# PROCEDIMIENTO
# Crea registro de lecturas
n = len (baliza)
registro = {}
for elemento in baliza:
    registro[elemento] = {'cuenta': 0 ,
                          'rssi': np.array([0.0,0.0]),
                          'snr':  np.array([0.0,0.0]),
                          'SumaRssi': np.array([0.0,0.0]),
                          'SumaSnr':  np.array([0.0,0.0]),
                          'minmaxRssi': np.array([0.0,-200.0,0.0,-200.0]),
                          'minmaxSnr': np.array([100.0,0.0,100.0,0.0])
                          }

# inicializa archivo.txt a vacio
archivo = open(nombrearchivo,'w')
archivo.close()  # Cierra el archivo

# Abre puerto Serial
arduino = serial.Serial(puerto, baudios)
arduino.setDTR(False)
time.sleep(0.3)

# limpia buffer de datos anteriores
arduino.flushInput()  
arduino.setDTR()  
time.sleep(0.3)
print('\nEstado del puerto: ',arduino.isOpen())
print('Nombre del dispositivo conectado: ', arduino.name)
print('Dump de la configuración:\n ',arduino)
print('\n###############################################\n')

# Lectura de datos
np.set_printoptions(precision=2)
conteo = 0
difunde = 0
while conteo<lecturasMin:
    #esperar hasta recibir un paquete
    while (arduino.inWaiting()==0):
        pass 

    # leer linea desde puerto serial
    linea = arduino.readline()
    # binario a texto, elimina /r/n
    texto = linea.decode()
    texto = linea.strip()
    
    # identificar la trama como rx, tx
    cond1 = texto.startswith('tx')
    cond2 = texto.startswith('rx')
    if cond1 or cond2:
        archivo = open(nombrearchivo,'a')
        archivo.write(texto+'\n')
        archivo.close()
        
        if (texto.startswith('tx')):
            difunde = difunde + 1
        if (texto.startswith('rx')):
            texto = texto.split(',')
            tipo  = texto[0]
            dir_remite = texto[2]
            paqrcbvID  = texto[3]
            rssi_tx    = float(texto[4])
            snr_tx     = float(texto[5])
            rssi_rx    = float(texto[6])
            snr_rx     = float(texto[7])
            if tipo == "rx":
                # conteo de lecturas
                cual = dir_remite
                registro[cual]['cuenta']=registro[cual]['cuenta']+1
                if registro[cual]['cuenta']<conteo:
                    conteo = registro[cual]['cuenta']

                # acumulado
                registro[cual]['SumaRssi'][0] = registro[cual]['SumaRssi'][0]+rssi_tx
                registro[cual]['SumaSnr'][0]  = registro[cual]['SumaSnr'][0]+snr_tx
                registro[cual]['SumaRssi'][1] = registro[cual]['SumaRssi'][1]+rssi_rx
                registro[cual]['SumaSnr'][1]  = registro[cual]['SumaSnr'][1]+snr_rx

                # promedios
                cuantos = registro[cual]['cuenta']
                registro[cual]['rssi'] = registro[cual]['SumaRssi']/cuantos
                registro[cual]['snr']  = registro[cual]['SumaSnr']/cuantos

                # minimos y maximos
                registro[cual]['minmaxRssi'][0] = np.min([rssi_tx,registro[cual]['minmaxRssi'][0]])
                registro[cual]['minmaxRssi'][1] = np.max([rssi_tx,registro[cual]['minmaxRssi'][1]])
                registro[cual]['minmaxRssi'][2] = np.min([rssi_rx,registro[cual]['minmaxRssi'][2]])
                registro[cual]['minmaxRssi'][3] = np.max([rssi_rx,registro[cual]['minmaxRssi'][3]])

                registro[cual]['minmaxSnr'][0] = np.min([snr_tx,registro[cual]['minmaxSnr'][0]])
                registro[cual]['minmaxSnr'][1] = np.max([snr_tx,registro[cual]['minmaxSnr'][1]])
                registro[cual]['minmaxSnr'][2] = np.min([snr_rx,registro[cual]['minmaxSnr'][2]])
                registro[cual]['minmaxSnr'][3] = np.max([snr_rx,registro[cual]['minmaxSnr'][3]])

            # Muestra en pantalla el estado de recepción
            print('\n difusion: ',difunde)
            print(texto)
            for elemento in baliza:
                print(elemento,registro[elemento]['cuenta'],
                      '\tprom Rssi[tx,rx] \t   Snr[tx,rx]')
                print("prom   :",registro[elemento]['rssi'],"\t  ",
                      registro[elemento]['snr'])
                print("min,max:",registro[elemento]['minmaxRssi'],
                      registro[elemento]['minmaxSnr'])

# Cerrar el puerto serial.
serial.Serial.close

2.1 Captura datos. Balizas y dispositivo – configuración algoritmo.ino

Instrucciones para el dispositivo en Arduino IDE

El esquema para baliza y dispositivo propuesto para el modelo de localización usa el modelo multipunto de LoRa que usa direcciones para identificar la función de cada elemento:

  • Baliza: d1,d2,d3
  • Dispositivo: c1

El esquema se muestra simplificado para un solo punto, se plantea usar el dispositivo ‘C1’ como punto de captura de datos hacia un archivo tipo texto mediante la conexión USB a un computador portatil.

En el dispositivo, por cada paquete recibido se toman los datos del paquete y valores medidos, creando un registro de trama con: el tipo rx/tx, id receptor, remitente, los datos recibidos y lecturas de RSSI y SNR del paquete recibido desde el remitente que es cada baliza.

    rssi_lorarx = LoRa.packetRssi();
    snr_lorarx = LoRa.packetSnr();

Registro de la trama de datos

El registro de la trama de datos recibida en el computador se compone de la parte de datos recibido y las lecturas de Rssi-SNR del paquete. El orden de los datos se muestra en la tabla:

esquema de trama de datos
 tipo receptor remitente dato1 dato2 dato3 Rssi recibido SNR recibido
rx, tx c1, ff d1, d2, d3 medido en el dispositivo medido en el dispositivo

Ejemplo de trama recibida:

rx,c1,d3,36,-62,6.00,-65,6.00
tx,ff,c1,2,-65,6.00,-65,6.00
rx,c1,d2,21,-130,6.00,-130,6.00
rx,c1,d1,138,-98,18.75,-131,10.25

Para obtener los valores de Rssi y SNR de recepción desde cada baliza, se emite periódicamente un paquete de difusión (‘ff’) desde el dispositivo c1. Los valores recibidos en la baliza se replican en los datos del próximo paquete de la baliza hacia el dispositivo c1, que corresponden a dato2 y dato3.

Para revisar la secuencia del paquete se registra un contador desde el emisor, enviado en la posición de dato1.

Configurar la función de baliza o dispositivo

se  establece mediante las direcciones de envío y recepción descritas en la sección principal del algoritmo.

Para el dispositivo la dirección local es C1 , para las balizas serán D1, D2, D3.

En el dispositivo se usa la dirección local ‘C1’, y como destino una trama de difusión ‘FF’ hacia todas las balizas de tal manera que cada baliza pueda medir el nivel de recepción de la señal emitida desde el dispositivo.

// Direcciones Enviar mensaje
byte dir_local   = 0xC1; // Concentrador 1
byte dir_destino = 0xFF; // Broadcast para medir Rssi
byte msjtxContador = 0;

En las balizas,  la dirección local de la baliza es D1, D2, D3, y la dirección destino es el dispositivo.

// Direcciones Enviar mensaje
byte dir_local   = 0xD1; // Baliza 1
byte dir_destino = 0xC1; // Dispositivo

Otros parámetros como tiempo entre paquetes, se configuran de acuerdo a las lecturas que se tomarán en un intervalo de tiempo. El parámetro t_ciclotx establece el intervalo de tiempo en que cada dispositivo o baliza emite un paquete.

Para LoRa, los tiempos de envio de paquetes son bastante largos de acuerdo al diseño básico de la red estándar, sin embargo durante el experimiento para tomar las lecturas más rápido en cada ubicación de prueba se puede reducir a por ejemplo, 10 segundos (10000ms)

int  t_ciclotx =10000;

otros detalles de configuración del dispositivo se pueden revisar en la sección de referencias al final de la hoja.


Módulos LoRa

El prototipo se implentó con los módulos de desarrollo Heltec Lora 32 V2 disponibles al momento de las pruebas. Los módulos son de 915MHz, y se usaron también con las pruebas de conexión multipunto descritas en LoRa Multipunto – Esquema


Instrucciones en Arduino

Las instrucciones del dispositivo escritas usando el IDE Arduino se presentan en secciones o bloques:

  • Principal: LoRa_Listener02
  • Receptor: LoRaReceptor
  • Transmisor: LoRaTransmisor
  • Sensor: Sensores


Sección principal

/* Pruebas de Rssi y SNR multipunto, esquema de triangulación
  Nodo- Receptor. http://blog.espol.edu.ec/girni/
  Septiembre-2020 */
#include "heltec.h"

// SENSOR Parpadeo
String sensorBlink = "ON"; // inicializa on/off: 1/0

// Direcciones Enviar mensaje
byte dir_local   = 0xC1; // Concentrador 1
byte dir_destino = 0xFF; // Broadcast para medir Rssi
byte msjtxContador = 0;
int  t_ciclotx = 10000;
long t_anterior = 0;
int  t_intervalo = 3000;

// Direcciones Mensaje Recibido
byte dir_envio = 0xC1;
int dir_remite = 0xFF; 
String paqueteRcb = "";
byte   paqrcbID = 0;
byte   paqrcbEstado = 0;
       // 0:vacio, 1: nuevo, 2:incompleto
       // 3:otro destinatario, 4:Broadcast

int  rssi_lorarx = 0;
float snr_lorarx = 0;

// Mensajes por Puerto Serial
volatile boolean serial_msj = true;

#define BAND    915E6  // 868E6,915E6

void setup() {
  Heltec.begin(false /*DisplayEnable Enable*/,
               true /*Heltec.LoRa Disable*/,
               serial_msj /*Serial Enable*/,
               true /*PABOOST Enable*/,
               BAND /*long BAND*/);
  LoRa.setSpreadingFactor(8);
  LoRa.receive();
}

void loop() {

  // LoRa enviar mensajes entre intervalos
  long t_ahora = millis();
  long t_transcurrido = t_ahora - t_anterior;
 
  if (t_transcurrido >= t_intervalo){
    t_anterior = t_ahora;
    t_intervalo = t_ciclotx + random(1000);
    
    // sensorParpadea(); //actualiza sensor
    sensorBlink = String(rssi_lorarx)+","+String(snr_lorarx);
    
    String paqueteEnv = String(sensorBlink).c_str() ;
    enviarlora(dir_destino, dir_local,
               msjtxContador, paqueteEnv);
    yield();

    // mensaje a serial
    if (serial_msj==true){
      // muestra todos paqrcbEstado
      Serial.print("tx,");
      Serial.print(String(dir_destino,HEX)+",");
      Serial.print(String(dir_local,HEX)+",");
      Serial.print(String(msjtxContador)+",");
      Serial.print(String(paqueteEnv)+",");
      Serial.print(String(rssi_lorarx)+",");
      Serial.println(snr_lorarx);
    }
    
    msjtxContador = msjtxContador + 1;
    
    // LED parpadea envio LoRa
    digitalWrite(LED, HIGH); delay(100);
    digitalWrite(LED, LOW);  delay(100);
    yield(); // procesa wifi
  }

  // LoRa revisar mensajes entrantes
  int msjRcbLoRa = LoRa.parsePacket();
  if (msjRcbLoRa !=0){ //¿mensaje no vacio?
    recibirlora(msjRcbLoRa);
    rssi_lorarx = LoRa.packetRssi();
    snr_lorarx = LoRa.packetSnr();
    
    if (serial_msj==true){
      // muestra todos los tipos
      //Serial.print("Paquete recibido, Estado: ");
      //Serial.println(paqrcbEstado);
      Serial.print("rx,");
      Serial.print(String(dir_envio,HEX)+",");
      Serial.print(String(dir_remite,HEX)+",");
      Serial.print(String(paqrcbID)+",");
      Serial.print(String(paqueteRcb)+",");
      Serial.print(String(rssi_lorarx)+",");
      Serial.println(String(snr_lorarx));
    }
    yield(); // procesa wifi
  }
}

Sección Receptor

En esta sección se revisa la trama, separando cada una de las partes del paquete recibido: dirección de envío, dirección de remitente, identificador del paquete generado por un contador de secuencia.

Se verifica que el destinatario del paquete, si el tamaño de la trama indica que el mensaje está completo, si la trama es de difusión (broadcast) usada para registrar los niveles de recepción.

void recibirlora(int tamano){
  if (tamano == 0){ 
    paqrcbEstado = 0; //vacio
    return;
    }
    
  // lectura de paquete
  paqueteRcb = "";
  dir_envio = LoRa.read();
  dir_remite  = LoRa.read();
  paqrcbID = LoRa.read();
  byte paqrcbTamano = LoRa.read();
  while(LoRa.available()){
    paqueteRcb += (char)LoRa.read();
  }
  if (paqrcbTamano != paqueteRcb.length()){
    paqrcbEstado = 2; // Tamaño incompleto
    return;
  }
  if (dir_envio != dir_local){
    if (dir_envio == 0xFF) {
      paqrcbEstado = 4; // Broadcast
    }else{
      paqrcbEstado = 3; // otro destino
    }
    return;
  }
  paqrcbEstado = 1;  // mensaje Nuevo
}

Sección de transmisor

En esta sección se verifica que el canal no se encuentre ocupado o que el dispositivo este listo para transmitir.

Se arma el paquete a enviar y se transmite.

void enviarlora(byte destino, byte remite,
                byte paqueteID, String paquete) {
  // espera que el radio esté listo
  // para enviar un paquete
  while(LoRa.beginPacket() == 0){
    if (serial_msj==true){
      Serial.println("Esperando radio disponible...");
    }
    yield(); // procesa wifi
    delay(100);
  }
  // envio del mensaje LoRa
  LoRa.beginPacket();
  LoRa.write(destino);
  LoRa.write(remite);
  LoRa.write(paqueteID);
  LoRa.write(paquete.length());
  LoRa.print(paquete);
  LoRa.endPacket();
}

Sección Sensor

La sección del sensor no es requerida para el registro de Rssi y SNR, por lo que se simula con un parpadeo de LED.

// Sensor Simulado
void sensorParpadea(){
    if (sensorBlink == "ON"){
      sensorBlink = "OFF";
    }else{
      sensorBlink = "ON";
    }
   return;
}

Referencias: LoRa Multipunto – Esquema

3.5 LoRa multipunto – Temperatura, Humedad: MQTT-HA

Esta es la última sección a realizar, pues se supone que tiene listo el dispositivo, construido y operativas las partes: Dispositivo, Gateway y Broker.

la visualización de los valores requiere declarar los dispositivos en Home Assistant en el archivo configuration.yaml

1. Incorporar el dispositivo en Home Assistant

Se requiere modificar el archivo configuration.yaml en el raspberry.

Se puede realizar en forma local desde el raspberry que tiene monitor, teclado y mouse conectado, y editar el archivo que se encuentra en el directorio:

 /home/homeassistant/.homeassistant/configuration.yaml

las líneas a añadir en la sección sensor:

sensor: 
  - platform: mqtt
    name: 'inv_D1_temperatura'
    unit_of_measurement: '°C'
    state_topic: 'invernadero/loraD1/temperatura'

  - platform: mqtt
    name: 'inv_D1_humedad'
    state_topic: 'invernadero/loraD1/humedad'
    unit_of_measurement: '%'

  - platform: mqtt
    name: 'inv_D1_bateria'
    state_topic: 'invernadero/loraD1/voltaje'
    unit_of_measurement: 'volt'

La configuración para añadir una tarjeta en la página es:

o añadiendo las intrucciones en el editor:

entities:
  - entity: sensor.inv_d1_temperatura
  - entity: sensor.inv_d1_humedad
  - entity: sensor.inv_d1_bateria
show_header_toggle: false
theme: default
title: Invernadero Dispositivo 1
type: entities

3.4 LoRa multipunto – Temperatura, Humedad: Gateway Archivo.ino

1. Instrucciones en Arduino IDE

Para el envío de los mensajes hacia el broker MQTT y Home Assistant se usa un dispositivo configurado como gateway.

Se usa un dispositivo en lugar de un concentrador en la etapa de prototipo considerando los costos involucrados. La próxima tarea es desarrollar el gateway usando un concentrador, cuyo valor es  más alto.

/* Gateway LoRa
  Lora/Wifi/MQTT/Home-Assistant
  Envia información por red Lora/Gateway WiFi
  hacia un broker MQTT y gestionar datos en Home-Assistant
  edelros@espol.edu.ec
  http://blog.espol.edu.ec/edelros/
  
  Referencia: Aaron.Lee www.heltec.cn
  https://github.com/Heltec-Aaron-Lee/WiFi_Kit_series
*/
#include "heltec.h"
#include <WiFi.h>
#include <PubSubClient.h>

// DISPOSITIVO LORA Banda ISM en Región 915Mhz
#define BAND  915E6 //433E6,868E6,915E6
// ranges from 6-12,default 7 see API docs
byte spread_factor = 8;

// LoRa Mensaje a enviar por direcciones
String paqueteEnv = "";
byte dir_local   = 0xC1; // Concentrador 1
byte dir_destino = 0xD1; // Dispositivo 1
byte msjContador = 0; // identificador de mensaje
// tiempo entre lecturas
long t_anterior = 0;
int  t_intervalo = 4000;

// LoRa Mensaje Recibido
byte dir_envio = 0xC1; // Concentrador 1
int dir_remite = 0xD0; // Inicia Remitente
String paqueteRcb = "";
byte   paqrcbID = 0;
byte   paqrcbEstado = 0;
  // 0:vacio, 1: nuevo, 2:incompleto
  // 3:otro destinatario, 4:Broadcast

 // Mensajes por Puerto Serial
volatile boolean serial_msj = true;

// WIFI: conexión a Router
char* ssid = "giotirni20";
char* password = "Anera2020@";

// MQTT: Servidor
char* MQTT_IP = "192.168.10.50";
uint16_t MQTT_puerto = 1883;
char* MQTT_usuario = "usuarioprueba";
char* MQTT_contrasena = "usuarioclave";

// MQTT: Dispositivo Sensor
char* MQTT_ID = "LoraGatewayC1";
char MQTT_TOPIC_T[50] = "invernadero/loraD1/temperatura";
char MQTT_TOPIC_H[50] = "invernadero/loraD1/humedad";
char MQTT_TOPIC_V[50] = "invernadero/loraD1/voltaje";
char MQTT_SensorEstado[10] = "OFF";
volatile boolean mqtt_desconectado = true;
// MQTT: Dispositivo Actuador
char* MQTT_COMMAND = "invernadero/loraD1/cambia";
char MQTT_ActuadorEstado[10] = "OFF";
volatile boolean actuador_estado = false;
volatile boolean actuador_bandera = false;
char temperatura[10]  = "00.00";
char humedad[10] = "00.00";
char voltaje[10] = "00.00";

// Clientes WiFi y MQTT
WiFiClient wificlient;
PubSubClient mqttclient(wificlient);

void setup(){
  Heltec.begin(false /*DisplayEnable Enable*/,
    true /*Heltec.Heltec.Heltec.LoRa Disable*/,
    serial_msj /*Serial Enable*/,
    true /*PABOOST Enable*/,
    BAND /*long BAND*/);
  // ranges from 6-12,default 7 see API docs
  LoRa.setSpreadingFactor(spread_factor);
  
  //LoRa.onReceive(cbk);
  LoRa.receive();
  
  // conexión WIFI y MQTT
  inicia_wifi();
  if (WiFi.status() == WL_CONNECTED){
    inicia_mqtt();
    }
}

void loop(){
  // parametros de recepción
  int rssi_lora = 0;
  int snr_lora = 0;
  
  // Revisa mensajes LoRa entrantes
  int msjRcbLoRa = LoRa.parsePacket();
  if (msjRcbLoRa !=0){
    
    recibirlora(msjRcbLoRa);
    rssi_lora = LoRa.packetRssi();
    snr_lora = LoRa.packetSnr();
       
    if (serial_msj==true){
      Serial.println("remite,msjID,mensaje,estado,Rssi,Snr");
      Serial.print(String(dir_remite, HEX)); Serial.print(",");
      Serial.print(paqrcbID); Serial.print(",");
      Serial.print(paqueteRcb); Serial.print(",");
      Serial.print(paqrcbEstado); Serial.print(",");
      Serial.print(rssi_lora); Serial.print(",");
      Serial.println(snr_lora);
    }
    yield(); // procesa wifi
    
    // LED parpadea Rebibido Lora
    digitalWrite(LED, HIGH); delay(50);
    digitalWrite(LED, LOW); delay(50);
    digitalWrite(LED, HIGH); delay(50);
    digitalWrite(LED, LOW);
    yield(); // procesa wifi
    delay(100);
  }

  // Procesa a MQTT mensaje completo
  if (msjRcbLoRa !=0 && paqrcbEstado == 1){
    // Separa parámetros
    String t = paqueteRcb.substring(1,6);
    String h = paqueteRcb.substring(8,10);
    String v = paqueteRcb.substring(12);
    
    // procesa tópico MQTT
    
    // añade dispositivo
    String topico = "invernadero/lora";
    String remite = String(dir_remite, HEX);
    remite.toUpperCase();
    topico = topico + remite;
    
    // procesa topico
    String topicot = topico + "/temperatura";
    String topicoh = topico + "/humedad";
    String topicov = topico + "/voltaje";
    topicot.toCharArray(MQTT_TOPIC_T,topicot.length()+1);
    topicoh.toCharArray(MQTT_TOPIC_H,topicoh.length()+1);
    topicov.toCharArray(MQTT_TOPIC_V,topicov.length()+1);

    t.toCharArray(temperatura,t.length()+1);
    h.toCharArray(humedad,h.length()+1);
    v.toCharArray(voltaje,v.length()+1);
    Serial.println(topicov);

    
   publica_estado();
  }
  yield(); // procesa wifi
  delay(20);

  // reenviar a dispositivo
  if (actuador_bandera == true){
    msjContador = msjContador +1;
    enviarlora(dir_destino, dir_local, 
               msjContador, paqueteEnv);
    actuador_bandera = false;
  }
  yield(); // procesa wifi
  delay(20);
  
  if (WiFi.status() != WL_CONNECTED){
    inicia_wifi();
  }else{
    if (mqttclient.connected()==false){
      mqtt_desconectado = true;
      inicia_mqtt(); // reintento
    }
    if (mqttclient.connected()==true){
      if (mqtt_desconectado==true){
        publica_estado();
        mqtt_desconectado=false;
      }
      mqttclient.loop();
    }
  }
  yield(); // procesa wifi
}


void enviarlora(byte destino, byte remite,
                byte paqueteID, String paquete){
  // espera que el radio esté listo
  // para enviar un paquete
  while(LoRa.beginPacket() == 0){
    if (serial_msj==true){
      Serial.println("Esperando radio disponible...");
    }
    yield(); // procesa wifi
    delay(100);
  }
  // envio del mensaje LoRa
  LoRa.beginPacket();
  LoRa.write(destino);
  LoRa.write(remite);
  LoRa.write(paqueteID);
  LoRa.write(paquete.length());
  LoRa.print(paquete);
  LoRa.endPacket();
}

void recibirlora(int tamano){
  if (tamano == 0){ 
    paqrcbEstado = 0; //vacio
    return;
  }
    
  // lectura de paquete
  paqueteRcb = "";
  dir_envio = LoRa.read();
  dir_remite  = LoRa.read();
  paqrcbID = LoRa.read();
  byte paqrcbTamano = LoRa.read();
  while(LoRa.available()){
    paqueteRcb += (char)LoRa.read();
  }
  
  if (paqrcbTamano != paqueteRcb.length()){
    paqrcbEstado = 2; // Tamaño incompleto
    return;
  }
  if (dir_envio != dir_local){
    paqrcbEstado = 3; // otro destino
    return;
  }
  if (dir_envio == 0xFF) {
    paqrcbEstado = 4; // Broadcast
    return;
  }
  paqrcbEstado = 1;  // mensaje Nuevo
}

void inicia_mqtt(void){
  int intentosmqtt = 5;
  int cuentamqtt = 0;
  
  if (serial_msj){
    Serial.print(" MQTT Conectando a ");
    Serial.println(MQTT_IP);
    }

  mqttclient.setServer(MQTT_IP, MQTT_puerto);
  mqttclient.connect(MQTT_ID, MQTT_usuario, MQTT_contrasena);
  mqttclient.setCallback(recibirmqtt);
  
  while (!mqttclient.connected() && (cuentamqtt<=intentosmqtt)){
    if (serial_msj){
      Serial.print(".");
    }
    cuentamqtt = cuentamqtt + 1;
    // LED Monitor parpadeo MQTT
    digitalWrite(LED, HIGH); delay(200);
    digitalWrite(LED, LOW); delay(200);
  }
  if (mqttclient.connected()){
      publica_estado();
  }
  if (serial_msj){
    //Fin de "...."
    Serial.println();
    Serial.print(" MQTT Conectado: ");
    Serial.print(mqttclient.connected());
    Serial.print("\t MQTT Estado: ");
    Serial.println(mqttclient.state());
  }
}

void publica_estado() { 

  if (mqttclient.connected()==true){
    mqttclient.publish(MQTT_TOPIC_T,temperatura,true);
    mqttclient.publish(MQTT_TOPIC_H,humedad,true);
    mqttclient.publish(MQTT_TOPIC_V,voltaje,true);
    mqttclient.subscribe(MQTT_COMMAND);
  }else{
    mqtt_desconectado = true;
  }
}

// llega mensaje MQTT, callback mqtt
void recibirmqtt(char* p_topic, byte* p_payload,
                  unsigned int p_length) {
  Serial.println("un mensaje mqtt");
  Serial.println(p_topic);
  // convierte a texto
  String payload;
  for (uint8_t i = 0; i < p_length; i++) {
    payload.concat((char)p_payload[i]);
    }
  // String dispositivo = p_topic[16] + String(p_topic[17]);
  char dispositivo[3] = "D0";
  dispositivo[1]=p_topic[17];
  dir_destino = (int) strtol(dispositivo,NULL,16);
  paqueteEnv = payload;
  actuador_bandera = true;

  if (mqttclient.connected()==true){
    mqttclient.subscribe(MQTT_COMMAND);
  }else{
    mqtt_desconectado = true;
  }
}

void inicia_wifi(void) {
  int intentoswifi = 10;
  int cuentawifi = 0;
  
  if (serial_msj){
    Serial.print(" WiFi Conectando a ");
    Serial.println(ssid);
    }
  
  WiFi.disconnect(true);
  delay(1000);
  WiFi.mode(WIFI_STA);
  WiFi.setAutoConnect(true);
  WiFi.begin(ssid,password);
  delay(100);
  
  while(WiFi.status() != WL_CONNECTED && 
             cuentawifi < intentoswifi){
    if (serial_msj){
      Serial.print(".");
      }
    cuentawifi = cuentawifi + 1;
    // Parpadeo de Monitor Wifi
    digitalWrite(LED, HIGH);delay(300);
    digitalWrite(LED, LOW);delay(200);
  }
  if (serial_msj){
    // mensaje a serial
    Serial.println();
    if (WiFi.status() == WL_CONNECTED){
      Serial.print(" Estado: ");
      Serial.println(WiFi.status());
      Serial.print(" MAC: ");
      Serial.println(WiFi.macAddress());
      Serial.print(" IP: ");
      Serial.println(WiFi.localIP());
      Serial.print(" RSSI: ");
      Serial.println(WiFi.RSSI());
      Serial.println();
    }
    if (WiFi.status() != WL_CONNECTED){
        WiFi.printDiag(Serial);
        Serial.println();
    }
  }
}

3.3 LoRa multipunto – Temperatura, Humedad: Dispositivo Archivo.ino

1. Instrucciones en Arduino IDE

Para facilitar la programación, se separan en funciones las acciones para enviar y recibir mensajes Lora y las de manejo de sensor/actuador.

El sensor DTH-11 se conecta al pin 13

El sensor de bateria se conecta al pin 36 como entrada analógica.

/*
  Dispositivo Sensor Temperatura y Humedad con DHT11
  Broker: MQTT/Home-Assistant
  Red ruta: LoRa/WiFi/Ethernet
  edelros@espol.edu.ec
  http://blog.espol.edu.ec/edelros/
  Referencia: Ejemplos de Aaron.Lee www.heltec.cn
*/
#include "heltec.h"
#include "DHT.h"

// Sensor de Temperatura&Humedad
#define DHTPIN 13
#define DHTTYPE DHT11 
DHT dht(DHTPIN, DHTTYPE);
String temperatura = "";
String humedad = ""; 

// Sensor de Bateria
# define BattPIN 36
String battNivel;

//Banda LoRa 915Mhz - ISM en Región 
#define BAND  915E6 // 433E6,868E6,915E6
// ranges from 6-12,default 7 see API docs
byte spread_factor = 8;

// Mensaje a enviar por direcciones
byte dir_local   = 0xD1; // Dispositivo  1
byte dir_destino = 0xC1; // Concentrador 1
// identificador de mensaje
byte msjContador = 0;
// tiempo entre lecturas
long t_anterior = 0;
int  t_intervalo = 5000;

// Mensaje Recibido
byte dir_envio = 0xC1; // Concentrador 1
int dir_remite = 0xD0; // Inicia Remitente
String paqueteRcb = "";
byte   paqrcbvID = 0;
byte   paqrcbEstado = 0;
  // 0:vacio, 1: nuevo, 2:incompleto
  // 3:otro destinatario, 4:Broadcast

 // Mensajes por Puerto Serial
volatile boolean serial_msj = true;

void setup(){
  Heltec.begin(false /*DisplayEnable Enable*/,
    true /*Heltec.Heltec.Heltec.LoRa Disable*/,
    serial_msj /*Serial Enable*/,
    true /*PABOOST Enable*/,
    BAND /*long BAND*/);
    
  // ranges from 6-12,default 7 see API docs
  LoRa.setSpreadingFactor(spread_factor);
  //LoRa.onReceive(cbk);
  LoRa.receive();

  //inicializa sensores
  pinMode(DHTPIN, INPUT);
  dht.begin();
  pinMode(BattPIN, INPUT);
  }

void loop(){
  // Enviar mensajes entre intervalos
  long t_ahora = millis();
  long t_transcurrido = t_ahora - t_anterior;

  // parametros de recepción
  int rssi_lora = 0;
  int snr_lora = 0;
  
  if (t_transcurrido >= t_intervalo){
    sensorLeeDHT(); //actualiza estado del sensor
    sensorBateria(); // actualiza estado de bateria
    
    // Construye paquete a enviar
    String paqueteEnv = "";
    paqueteEnv = paqueteEnv + "t" + temperatura;
    paqueteEnv = paqueteEnv + "|";
    paqueteEnv = paqueteEnv + "h" + humedad;
    paqueteEnv = paqueteEnv + "|";
    paqueteEnv = paqueteEnv + "v" + battNivel;
    
    enviarlora(dir_destino, dir_local,
               msjContador, paqueteEnv);
    
    msjContador = msjContador + 1;
    
    // parametros de recepción
    rssi_lora = LoRa.packetRssi();
    snr_lora = LoRa.packetSnr();
    yield();
    
    // mensaje a serial
    if (serial_msj==true){
      Serial.print(String(dir_destino,HEX));
      Serial.print(",");
      Serial.print(String(dir_local,HEX));
      Serial.print(",");
      Serial.print(msjContador);Serial.print(",");
      Serial.print(paqueteEnv.length());Serial.print(",");
      Serial.print(paqueteEnv);Serial.print(",");
      Serial.print(rssi_lora);Serial.print(",");
      Serial.println(snr_lora);
    }
    
    t_anterior = millis();
    t_intervalo = 3000 + random(2000);
    
    // LED parpadea envio lora
    digitalWrite(LED, HIGH); delay(100);
    digitalWrite(LED, LOW);  delay(100);
    yield(); // procesa wifi
  }

  // Revisar mensajes LoRa entrantes
  int msjRcbLoRa = LoRa.parsePacket();
  if (msjRcbLoRa !=0){
    recibirlora(msjRcbLoRa);
    rssi_lora = LoRa.packetRssi();
    snr_lora = LoRa.packetSnr();
    
    if (serial_msj==true){
      if (paqrcbEstado == 1){
        Serial.println("Mensaje: " + paqueteRcb);
        Serial.println("RSSI: " + String(rssi_lora));
        Serial.println("Snr: " + String(snr_lora));
        Serial.println();
      }else{
        Serial.print("Paquete recibido Estado: ");
        Serial.println(paqrcbEstado);
      }
    }
    yield(); // procesa wifi
    
    // LED parpadea Recibido Lora
    digitalWrite(LED, HIGH); delay(50);
    digitalWrite(LED, LOW); delay(50);
    digitalWrite(LED, HIGH); delay(50);
    digitalWrite(LED, LOW);
  }
  delay(100);
  yield(); // procesa wifi
}

void enviarlora(byte destino, byte remite,
                byte paqueteID, String paquete){
  // espera que el radio esté listo
  // para enviar un paquete
  while(LoRa.beginPacket() == 0){
    if (serial_msj==true){
      Serial.println("Esperando radio disponible...");
    }
    yield(); // procesa wifi
    delay(100);
  }
  // envio del mensaje LoRa
  LoRa.beginPacket();
  LoRa.write(destino);
  LoRa.write(remite);
  LoRa.write(paqueteID);
  LoRa.write(paquete.length());
  LoRa.print(paquete);
  LoRa.endPacket();
}

void recibirlora(int tamano){
  if (tamano == 0){ 
    paqrcbEstado = 0; //vacio
    return;
    }
  // lectura de paquete
  paqueteRcb = "";
  dir_envio = LoRa.read();
  dir_remite  = LoRa.read();
  paqrcbvID = LoRa.read();
  byte paqrcbTamano = LoRa.read();
  while(LoRa.available()){
    paqueteRcb += (char)LoRa.read();
  }
  if (paqrcbTamano != paqueteRcb.length()){
    paqrcbEstado = 2; // Tamaño incompleto
    return;
  }
  if (dir_envio != dir_local){
    paqrcbEstado = 3; // otro destino
    return;
  }
  if (dir_envio == 0xFF) {
    paqrcbEstado = 4; // Broadcast
    return;
  }
  paqrcbEstado = 1;  // mensaje Nuevo
}

// Sensor lecturas
void sensorLeeDHT(){
  humedad = String(int(dht.readHumidity())).c_str();
  temperatura = String(dht.readTemperature()).c_str();
}

void sensorBateria(){
  int lectura = analogRead(BattPIN);
  // convierte a equivalente en voltios
  float voltaje = (float(lectura)/4096.0)*4.2*3.0/2.0;
  battNivel = String(voltaje).c_str();
}

3.2 LoRa multipunto – Temperatura, Humedad: Esquemático LoRa32

El primer prototipo se realiza usando una placa de desarrollo Heltec Lora 32. EL punto de partida es LoRa multipunto añadiendo los siguientes componentes:

  • Sensor de Temperatura y Humedad DHT-11, puede ser actualizado a DHT-22 para mayor precisión. Se empieza con DHT-11 por tenerlo disponible al inicio.
  • Bateria para añadir portabilidad
  • Panel solar para alimentación, usando un módulo de carga de batería.

3.1 LoRa multipunto – Temperatura, Humedad

Presentación

El sensor numérico más sencillo de implementar es del de temperatura-Humedad con el sensor DHT-11 o DHT-22. Un requerimiento para el  dispositivo es operar a batería, con opción de carga con un panel solar.

El uso de la batería limita el consumo de energía, los módulos LoRa al ser de bajo consumo son los seleccionados.

El punto de partida es la configuración LoRa Multipunto, que modificando e valor del sensor enviado y los elementos correspondientes en MQTT, permite visualizar en la página del broker Home-assistant el valor del sensor.