Sensly/Sensly_Calibration.py at master · Altitude-Tech/Sensly · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
# File Description: puts into Sensly_Calibration_<timeStamp>.csv values from Sensors
from time import *

import RPi.GPIO as GPIO
from Sensors import Sensor, Gas
from bme_combo import *
import logging
import sys

# Sensly Constants
R0 = [0,0,0]
RSAir = [9.5,27,3.62]
RLOAD = 10000               # In Ohms
SAMPLES_COUNT_FOR_VALUES_TO_STORE = 300 # Times the loop getting data will turn to get an average of data before writting into log file (Setted by AltitudeTech to 300)
SECONDS_TO_WAIT_TILL_HAT_IS_HOT = 600 # Time in seconds the script will wait till the HAT is @ working temperature (Setted by Altitude to 600)

MQ2 = Sensor('MQ2',R0[0],RSAir[0]) # name, max adc value, Calibrated R0 value, load resistance
MQ7 = Sensor('MQ7',R0[1],RSAir[1])
MQ135 = Sensor('MQ135',R0[2],RSAir[2])
ADCMax = 4095


def initialize():
	"""Sets the logging variable & Resets GPIO"""
	logging.basicConfig(stream=sys.stdout, level=logging.DEBUG, format='%(asctime)s.%(msecs)03d %(levelname)s %(module)s - %(funcName)s: %(message)s', datefmt="%Y-%m-%d %H:%M:%S")
	# Set up GPIO ports
	GPIO.setmode(GPIO.BCM)      # Use BCM numbering
	GPIO.setup(23, GPIO.OUT)
	# Reset Sensly HAT
	GPIO.output(23, False)      # Set GPIO Pin 4 to low
	time.sleep(1)
	GPIO.output(23, True)       # Set GPIO Pin 4 to low
	# Clean up GPIO
	GPIO.cleanup()

def gpioReset():
	"""Resets GPIO & HAT"""
	GPIO.setmode(GPIO.BCM)      # Use BCM numbering
	GPIO.setup(23, GPIO.OUT)
	# Reset Sensly HAT
	GPIO.output(23, False)      # Set GPIO Pin 4 to low
	#GPIO.output(25, True)      # Set GPIO Pin 4 to low
	time.sleep(0.5)
	GPIO.output(23, True)       # Set GPIO Pin 4 to low
	time.sleep(0.5)


initialize()
datafile = time.strftime('./Sensly_Calibration_%d-%m-%Y_%H_%M_%S.csv')
logging.debug("Writting data to " + datafile)
with open(datafile, 'w+') as f1:
    f1.write('Time, MQ2RZero, MQ7RZero, MQ135RZero\n')
try:
    # Set commands for getting data from snesors
    MQ2cmd = 0x01
    MQ7cmd = 0x02
    MQ135cmd = 0x03
    # Initialise the last stored R0 value used for the running average
    L_AvRs_MQ2 = 0
    L_AvRs_MQ7 = 0
    L_AvRs_MQ135 = 0
    logging.info("Waiting for " + str(SAMPLES_COUNT_FOR_VALUES_TO_STORE) + " seconds the HAT to be @ working temperature...")
    sleep(SECONDS_TO_WAIT_TILL_HAT_IS_HOT) # 600 seconds to get the HAT in Working temperature
    logging.info("Enterring to a True loop, Ctrl+c to stop it...")
    callIndex = 0
    while True:
        data = []
        # Get current time and add data array
        data.append(time.strftime('%H:%M:%S'))

        # Initalise Average resistance values
        AvRs_MQ2 = 0
        AvRs_MQ7 = 0
        AvRs_MQ135 = 0

        logging.debug("Starting calibration for an average of " + str(SAMPLES_COUNT_FOR_VALUES_TO_STORE))
        # Run calibration
        for x in range(SAMPLES_COUNT_FOR_VALUES_TO_STORE):
	    callIndex = callIndex + 1
	    gpioReset()
	    logging.debug("Trying to get Data from MQ2 Sensor 1")
            AvRs_MQ2 = AvRs_MQ2 + MQ2.Get_RS(MQ2cmd)
            logging.debug("Got first data AvRs_MQ2:" + str(AvRs_MQ2))
            sleep(0.01)
	    gpioReset()
	    logging.debug("Trying to get Data from MQ7 Sensor 2")
            AvRs_MQ7 = AvRs_MQ7 + MQ7.Get_RS(MQ7cmd)
            logging.debug("Got data AvRs_MQ7:" + str(AvRs_MQ7))
            sleep(0.01)
	    gpioReset()
	    logging.debug("Trying to get Data from MQ135 Sensor 3")
            AvRs_MQ135 = AvRs_MQ135 + MQ135.Get_RS(MQ135cmd)
            logging.debug("Got data MQ135:" + str(AvRs_MQ135))
            sleep(0.01)
	    logging.debug("----> Called " + str(callIndex) + " times")

	# Setting all averages
        AvRs_MQ2 = AvRs_MQ2/SAMPLES_COUNT_FOR_VALUES_TO_STORE
        AvRs_MQ7 = AvRs_MQ7/SAMPLES_COUNT_FOR_VALUES_TO_STORE
        AvRs_MQ135 = AvRs_MQ135/SAMPLES_COUNT_FOR_VALUES_TO_STORE

        # Find R0
        MQ2.R0 = ((AvRs_MQ2/MQ2.RSAIR) + L_AvRs_MQ2)/2
        MQ7.R0 = ((AvRs_MQ7/MQ7.RSAIR) + L_AvRs_MQ7)/2
        MQ135.R0 = ((AvRs_MQ135/MQ135.RSAIR) + L_AvRs_MQ135)/2

        L_AvRs_MQ2 =  MQ2.R0
        L_AvRs_MQ7 = MQ7.R0
        L_AvRs_MQ135 = MQ135.R0

        # Store the values in a file
        data.append(MQ2.R0)
        data.append(MQ7.R0)
        data.append(MQ135.R0)

        with open(datafile, 'a') as f2:
	    toWrite = ','.join(str(d) for d in data) + '\n'
            f2.write(toWrite)
	    logging.info("Writting to datafile:" + toWrite)

except KeyboardInterrupt:
    print "Bye"
finally:
    GPIO.cleanup()