Add files via upload

Added BME support and persistent USB names.

apsetup.sh

@@ -22,8 +22,8 @@ wpa=2
 wpa_key_mgmt=WPA-PSK
 wpa_pairwise=CCMP
 rsn_pairwise=CCMP
-wpa_passphrase=PASSWORD
-ssid=NAME
+wpa_passphrase=lighthouse
+ssid=freyja
 EOF
 
 sed -i -- 's/exit 0/ /g' /etc/rc.local

bme.py

@@ -0,0 +1,221 @@
+#coding: utf-8
+import math
+import smbus
+import time
+import socket
+import operator
+
+bus_number  = 1
+i2c_address = 0x77
+
+bus = smbus.SMBus(bus_number)
+
+digT = []
+digP = []
+digH = []
+
+t_fine = 0.0
+
+
+def writeReg(reg_address, data):
+	bus.write_byte_data(i2c_address,reg_address,data)
+
+def get_calib_param():
+	calib = []
+
+	for i in range (0x88,0x88+24):
+		calib.append(bus.read_byte_data(i2c_address,i))
+	calib.append(bus.read_byte_data(i2c_address,0xA1))
+	for i in range (0xE1,0xE1+7):
+		calib.append(bus.read_byte_data(i2c_address,i))
+
+	digT.append((calib[1] << 8) | calib[0])
+	digT.append((calib[3] << 8) | calib[2])
+	digT.append((calib[5] << 8) | calib[4])
+	digP.append((calib[7] << 8) | calib[6])
+	digP.append((calib[9] << 8) | calib[8])
+	digP.append((calib[11]<< 8) | calib[10])
+	digP.append((calib[13]<< 8) | calib[12])
+	digP.append((calib[15]<< 8) | calib[14])
+	digP.append((calib[17]<< 8) | calib[16])
+	digP.append((calib[19]<< 8) | calib[18])
+	digP.append((calib[21]<< 8) | calib[20])
+	digP.append((calib[23]<< 8) | calib[22])
+	digH.append( calib[24] )
+	digH.append((calib[26]<< 8) | calib[25])
+	digH.append( calib[27] )
+	digH.append((calib[28]<< 4) | (0x0F & calib[29]))
+	digH.append((calib[30]<< 4) | ((calib[29] >> 4) & 0x0F))
+	digH.append( calib[31] )
+
+	for i in range(1,2):
+		if digT[i] & 0x8000:
+			digT[i] = (-digT[i] ^ 0xFFFF) + 1
+
+	for i in range(1,8):
+		if digP[i] & 0x8000:
+			digP[i] = (-digP[i] ^ 0xFFFF) + 1
+
+	for i in range(0,6):
+		if digH[i] & 0x8000:
+			digH[i] = (-digH[i] ^ 0xFFFF) + 1  
+
+def readData():
+	data = []
+	for i in range (0xF7, 0xF7+8):
+		data.append(bus.read_byte_data(i2c_address,i))
+	pres_raw = (data[0] << 12) | (data[1] << 4) | (data[2] >> 4)
+	temp_raw = (data[3] << 12) | (data[4] << 4) | (data[5] >> 4)
+	hum_raw  = (data[6] << 8)  |  data[7]
+
+	temperature = compensate_T(temp_raw)
+	pressure = compensate_P(pres_raw)
+	humidity = compensate_H(hum_raw)
+        return dict(temperature=temperature,
+                    pressure=pressure,
+                    humidity=humidity)
+
+def compensate_P(adc_P):
+	global  t_fine
+	pressure = 0.0
+
+	v1 = (t_fine / 2.0) - 64000.0
+	v2 = (((v1 / 4.0) * (v1 / 4.0)) / 2048) * digP[5]
+	v2 = v2 + ((v1 * digP[4]) * 2.0)
+	v2 = (v2 / 4.0) + (digP[3] * 65536.0)
+	v1 = (((digP[2] * (((v1 / 4.0) * (v1 / 4.0)) / 8192)) / 8)  + ((digP[1] * v1) / 2.0)) / 262144
+	v1 = ((32768 + v1) * digP[0]) / 32768
+
+	if v1 == 0:
+		return 0
+	pressure = ((1048576 - adc_P) - (v2 / 4096)) * 3125
+	if pressure < 0x80000000:
+		pressure = (pressure * 2.0) / v1
+	else:
+		pressure = (pressure / v1) * 2
+	v1 = (digP[8] * (((pressure / 8.0) * (pressure / 8.0)) / 8192.0)) / 4096
+	v2 = ((pressure / 4.0) * digP[7]) / 8192.0
+	pressure = pressure + ((v1 + v2 + digP[6]) / 16.0)  
+
+	#print "pressure : %7.2f hPa" % (pressure/100)
+        return "{:.2f}".format(pressure/100)
+
+def compensate_T(adc_T):
+	global t_fine
+	v1 = (adc_T / 16384.0 - digT[0] / 1024.0) * digT[1]
+	v2 = (adc_T / 131072.0 - digT[0] / 8192.0) * (adc_T / 131072.0 - digT[0] / 8192.0) * digT[2]
+	t_fine = v1 + v2
+	temperature = t_fine / 5120.0
+	#print "temp : %-6.2f ℃" % (temperature) 
+        return "{:.2f}".format(temperature)
+
+def compensate_H(adc_H):
+	global t_fine
+	var_h = t_fine - 76800.0
+	if var_h != 0:
+		var_h = (adc_H - (digH[3] * 64.0 + digH[4]/16384.0 * var_h)) * (digH[1] / 65536.0 * (1.0 + digH[5] / 67108864.0 * var_h * (1.0 + digH[2] / 67108864.0 * var_h)))
+	else:
+		return 0
+	var_h = var_h * (1.0 - digH[0] * var_h / 524288.0)
+	if var_h > 100.0:
+		var_h = 100.0
+	elif var_h < 0.0:
+		var_h = 0.0
+	#print "hum : %6.2f ％" % (var_h)
+        return "{:.2f}".format(var_h)
+
+def setup():
+	osrs_t = 1			#Temperature oversampling x 1
+	osrs_p = 1			#Pressure oversampling x 1
+	osrs_h = 1			#Humidity oversampling x 1
+	mode   = 3			#Normal mode
+	t_sb   = 5			#Tstandby 1000ms
+	filter = 0			#Filter off
+	spi3w_en = 0			#3-wire SPI Disable
+
+	ctrl_meas_reg = (osrs_t << 5) | (osrs_p << 2) | mode
+	config_reg    = (t_sb << 5) | (filter << 2) | spi3w_en
+	ctrl_hum_reg  = osrs_h
+
+	writeReg(0xF2,ctrl_hum_reg)
+	writeReg(0xF4,ctrl_meas_reg)
+	writeReg(0xF5,config_reg)
+
+setup()
+get_calib_param()
+
+BME_IP = "127.0.0.7"
+BME_PORT = 5005
+
+MON_IP = "127.0.0.8"
+MON_PORT = 5005
+mon = 0
+
+t_print = 0
+
+while True:
+
+    # health monitor
+    hack = time.time()
+    if hack - mon > .5:
+        sock = socket.socket(socket.AF_INET, socket.SOCK_DGRAM)
+        sock.sendto(str(hack), (MON_IP, MON_PORT))
+        mon = hack
+
+    if hack - t_print > 10:
+	try:
+        	f = open('dst_bus', 'r')
+                mtw = f.readline()
+                f.close()
+	except ValueError:
+                time.sleep(.03)
+                f = open('dst_bus', 'r')
+                mtw = f.readline()
+                f.close()
+
+        data = readData()
+        temp = round(float(data['temperature']), 1)
+        bars = round(float(data['pressure']) * .001, 5)
+        inch = round(bars * 29.53, 2)
+        humi = round(float(data['humidity']), 1)
+        dewp = round((temp + 273.15) * math.sqrt(1 + ((temp + 273.15) * math.log(humi / 100)) / ((3.167485 - .00244 * (temp + 273.15)) * 1000000 / 461.5)) - 273.15, 1)
+
+        # MDA
+        mda = "IIMDA," + str(inch) + ",I," + str(round(bars, 2)) + ",B," + str(temp) + ",C," + mtw + ",C," + str(humi) + ",," + str(dewp) + ",C,,,,,,,,"
+        mdacs = format(reduce(operator.xor,map(ord,mda),0),'X')
+        if len(mdacs) == 1:
+            mdacs = "0" + mdacs
+        iimda = "$" + mda + "*" + mdacs + "\r\n"
+
+        # Simulate NKE Instruments
+	xdrb = "IIXDR,P," + str(bars) + ",B,Barometer"
+	xdrbcs = format(reduce(operator.xor,map(ord,xdrb),0),'X')
+        if len(xdrbcs) == 1:
+            xdrbcs = "0" + xdrbcs
+	iixdrb = "$" + xdrb + "*" + xdrbcs + "\r\n"
+
+	xdrt = "IIXDR,C," + str(temp) + ",C,AirTemp"
+	xdrtcs = format(reduce(operator.xor,map(ord,xdrt),0),'X')
+        if len(xdrtcs) == 1:
+            xdrtcs = "0" + xdrtcs
+        iixdrt = "$" + xdrt + "*" + xdrtcs + "\r\n"
+
+	mmb = "IIMMB," + str(inch) + ",I," + str(bars) + ",B"
+	mmbcs = format(reduce(operator.xor,map(ord,mmb),0),'X')
+        if len(mmbcs) == 1:
+            mmbcs = "0" + mmbcs
+        iimmb = "$" + mmb + "*" + mmbcs + "\r\n"
+
+	mta = "IIMTA," + str(temp) + ",C"
+        mtacs = format(reduce(operator.xor,map(ord,mta),0),'X')
+        if len(mtacs) == 1:
+            mtacs = "0" + mtacs
+        iimta = "$" + mta + "*" + mtacs + "\r\n"
+
+	atm_sentence = iimda + iixdrb + iixdrt + iimmb + iimta
+
+        # to kplex
+        sock = socket.socket(socket.AF_INET, socket.SOCK_DGRAM)
+        sock.sendto(atm_sentence, (BME_IP, BME_PORT))
+
+        t_print = hack

dst.py

@@ -13,7 +13,6 @@
 MON_PORT = 5005
 
 mon = 0
-log = time.time()
 
 # initialize for VDR
 five = 0
@@ -28,14 +27,13 @@
 vlwfirst = 1
 vlwinit = 0.0
 
-# determine dst port
-f = open('gps_port', 'r')
-gps_port = f.readline()
+# initialize log
+log = time.time()
+f = open('dstraw', 'w')
+f.write("Last 60 Seconds DST Raw Input:\r\n")
 f.close()
-dst_port = "1"
-if gps_port == "1":
-    dst_port = "0"
-ser = serial.Serial('/dev/ttyUSB' + dst_port, 4800, timeout=5)
+
+ser = serial.Serial('/dev/dst', 4800, timeout=5)
 
 while True:
 
@@ -54,7 +52,7 @@
     if hack - log > 60:
         os.remove('dstraw')
         f = open('dstraw', 'w')
-        f.write("Serial port " + dst_port + ": " + time.asctime( time.localtime(time.time()) ) + "\r\n")
+        f.write("Last 60 Seconds DST Raw Input:\r\n")
         f.close()
         log = hack
 
@@ -113,21 +111,20 @@
                 sock.sendto(sddpt, (DST_IP, DST_PORT))
 
             # mean water temp sentence
-            if title == "YXMTW":
+            elif title == "YXMTW":
 
-                # convert to fahrenheit
-                mtw = dst_vars[0] + "," + str(int(float(dst_vars[1]) * 9 / 5 + 32)) + ",F"
-                cs = format(reduce(operator.xor,map(ord,mtw),0),'X')
-                if len(cs) == 1:
-                    cs = "0" + cs
-                yxmtw = "$" + mtw + "*" + cs + "\r\n"
+                # write to bus
+                mtw = dst_vars[1]
+		f = open('dst_bus', 'w')
+		f.write(mtw)
+		f.close()
 
                 # to kplex
                 sock = socket.socket(socket.AF_INET, socket.SOCK_DGRAM)
-                sock.sendto(yxmtw, (DST_IP, DST_PORT))
+                sock.sendto(dst_raw, (DST_IP, DST_PORT))
 
             # vessel waterspeed sentence and current set and drift
-            if title == "VWVHW":
+            elif title == "VWVHW":
 
                 # heading and roll from imu
                 try:
@@ -247,19 +244,12 @@
                         vdrcs = "0" + vdrcs
                     iivdr = "$" + vdr + "*" + vdrcs + "\r\n"
 
-                    # ghost MWV sentence to show set/drift instead of wind
-                    mwv = "IIMWV," + str(int(set_apparent)) + ",R," + str(round(drift,1)) + ",N,A"
-                    mwvcs = format(reduce(operator.xor,map(ord,mwv),0),'X')
-                    if len(mwvcs) == 1:
-                        mwvcs = "0" + mwvcs
-                    iimwv = "$" + mwv + "*" + mwvcs + "\r\n"
-
                     # to kplex
                     sock = socket.socket(socket.AF_INET, socket.SOCK_DGRAM)
-                    sock.sendto(iivdr + iimwv, (DST_IP, DST_PORT))
+                    sock.sendto(iivdr, (DST_IP, DST_PORT))
 
             # voyage log sentence
-            if title == "VWVLW":
+            elif title == "VWVLW":
 
                 # calculate present trip total vs overall total
                 if vlwfirst == 1:
@@ -275,3 +265,10 @@
                 # to kplex
                 sock = socket.socket(socket.AF_INET, socket.SOCK_DGRAM)
                 sock.sendto(vwvlw, (DST_IP, DST_PORT))
+
+	    # if it's any other valid NMEA sentence
+            else:
+
+                # to kplex
+                sock = socket.socket(socket.AF_INET, socket.SOCK_DGRAM)
+                sock.sendto(dst_raw, (DST_IP, DST_PORT))