Incorporate dump1090-mutability. UI signal strength and traffic improvements.

pull/292/head
AvSquirrel 2016-02-25 05:53:35 +00:00
rodzic 27dc75ca0e
commit 5d082e3292
11 zmienionych plików z 203 dodań i 171 usunięć

@ -1 +1 @@
Subproject commit b360db74b99a408cb83833af802042d29cccffd9
Subproject commit 7eca082524fa135eb8efd78780f25739efb6a82c

Wyświetl plik

@ -131,7 +131,7 @@ void make_atan2_table(void) {
double scaled_ang = round(32768 * ang / M_PI);
double amp = sqrt(d_i * d_i + d_q * d_q);
uint16_t scaled_amp = amp * 1000.0 / 127.5 + .5;
uint16_t scaled_amp = amp * 1000.0 / 127.5;
u.iq[0] = i;
u.iq[1] = q;
@ -150,13 +150,13 @@ static void convert_to_phi(uint16_t *dest, uint16_t *src, int n) {
dest[i] = iqphase[src[i]];
}
static void calc_power(uint16_t *samples, int len) {
static void calc_power(uint16_t *samples, int len) { // sets signal_strength to scaled amplitude. 0 = no signal, 1000 = saturated receiver on all samples in measurement.
long avg = 0;
int n = len;
while (n--) {
avg += iqamplitude[*samples++];
}
signal_strength = (avg + len / 2) / len;
signal_strength = avg / len;
}
#ifndef BUILD_LIB

Wyświetl plik

@ -20,6 +20,7 @@ import (
"io"
"io/ioutil"
"log"
"math"
"os"
"os/signal"
"runtime"
@ -112,7 +113,8 @@ type msg struct {
TimeReceived time.Time
Data []byte
Products []uint32
Signal_strength int
Signal_amplitude int
Signal_strength float64
ADSBTowerID string // Index in the 'ADSBTowers' map, if this is a parseable uplink message.
}
@ -125,9 +127,10 @@ var timeStarted time.Time
type ADSBTower struct {
Lat float64
Lng float64
Signal_strength_last_minute int
signal_power_last_minute int64 // Over total messages.
Signal_strength_max int
Signal_strength_now float64 // Current RSSI (dB)
Signal_strength_max float64 // all-time peak RSSI (dB) observed for this tower
Energy_last_minute uint64 // Summation of power observed for this tower across all messages last minute
Signal_strength_last_minute float64 // Average RSSI (dB) observed for this tower last minute
Messages_last_minute uint64
Messages_total uint64
}
@ -594,7 +597,7 @@ func relayMessage(msgtype uint16, msg []byte) {
func heartBeatSender() {
timer := time.NewTicker(1 * time.Second)
timerMessageStats := time.NewTicker(5 * time.Second)
timerMessageStats := time.NewTicker(2 * time.Second)
for {
select {
case <-timer.C:
@ -641,7 +644,7 @@ func updateMessageStats() {
// Clear out ADSBTowers stats.
for t, tinf := range ADSBTowers {
tinf.Messages_last_minute = 0
tinf.Signal_strength_last_minute = 0
tinf.Energy_last_minute = 0
ADSBTowers[t] = tinf
}
@ -656,8 +659,8 @@ func updateMessageStats() {
if len(MsgLog[i].ADSBTowerID) > 0 { // Update tower stats.
tid := MsgLog[i].ADSBTowerID
twr := ADSBTowers[tid]
twr.Energy_last_minute += uint64((MsgLog[i].Signal_amplitude) * (MsgLog[i].Signal_amplitude))
twr.Messages_last_minute++
twr.signal_power_last_minute += int64(MsgLog[i].Signal_strength)
if MsgLog[i].Signal_strength > twr.Signal_strength_max { // Update alltime max signal strength.
twr.Signal_strength_max = MsgLog[i].Signal_strength
}
@ -684,9 +687,9 @@ func updateMessageStats() {
// Update average signal strength over last minute for all ADSB towers.
for t, tinf := range ADSBTowers {
if tinf.Messages_last_minute == 0 {
tinf.Signal_strength_last_minute = 0
tinf.Signal_strength_last_minute = -99
} else {
tinf.Signal_strength_last_minute = int(tinf.signal_power_last_minute / int64(tinf.Messages_last_minute))
tinf.Signal_strength_last_minute = 10 * (math.Log10(float64((tinf.Energy_last_minute / tinf.Messages_last_minute))) - 6)
}
ADSBTowers[t] = tinf
}
@ -868,10 +871,9 @@ func parseInput(buf string) ([]byte, uint16) {
}
if s[0] == '-' {
parseDownlinkReport(s, int32(thisSignalStrength))
parseDownlinkReport(s, int(thisSignalStrength))
}
s = s[1:]
msglen := len(s) / 2
@ -901,7 +903,8 @@ func parseInput(buf string) ([]byte, uint16) {
thisMsg.MessageClass = MSGCLASS_UAT
thisMsg.TimeReceived = stratuxClock.Time
thisMsg.Data = frame
thisMsg.Signal_strength = thisSignalStrength
thisMsg.Signal_amplitude = thisSignalStrength
thisMsg.Signal_strength = 10 * (math.Log10(float64((thisSignalStrength * thisSignalStrength))) - 6)
thisMsg.Products = make([]uint32, 0)
if msgtype == MSGTYPE_UPLINK {
// Parse the UAT message.
@ -914,10 +917,13 @@ func parseInput(buf string) ([]byte, uint16) {
var newTower ADSBTower
newTower.Lat = uatMsg.Lat
newTower.Lng = uatMsg.Lon
newTower.Signal_strength_now = thisMsg.Signal_strength
newTower.Signal_strength_max = -999 // dBmax = 0, so this needs to initialize below scale ( << -48 dB)
ADSBTowers[towerid] = newTower
}
twr := ADSBTowers[towerid]
twr.Messages_total++
twr.Signal_strength_now = thisMsg.Signal_strength
ADSBTowers[towerid] = twr
// Get all of the "product ids".
for _, f := range uatMsg.Frames {

Wyświetl plik

@ -294,8 +294,10 @@ func messageQueueSender() {
averageSendableQueueSize = averageSendableQueueSize / float64(len(outSockets)) // It's a total, not an average, up until this point.
pd = math.Max(float64(1.0/750.0), float64(1.0/(4.0*averageSendableQueueSize))) // Say 250ms is enough to get through the whole queue.
} else {
pd = float64(1.0 / 0.1) // 100ms.
pd = float64(0.1) // 100ms.
}
log.Printf("Average sendable queue is %v messages. Changing queue timer to %f seconds\n", averageSendableQueueSize, pd)
queueTimer.Stop()
queueTimer = time.NewTicker(time.Duration(pd*1000000000.0) * time.Nanosecond)
lastQueueTimeChange = stratuxClock.Time

Wyświetl plik

@ -64,7 +64,7 @@ func readToChan(fp io.ReadCloser, ch chan []byte) {
func (e *ES) read() {
defer es_wg.Done()
log.Println("Entered ES read() ...")
cmd := exec.Command("/usr/bin/dump1090", "--net", "--device-index", strconv.Itoa(e.indexID), "--ppm", strconv.Itoa(e.ppm))
cmd := exec.Command("/usr/bin/dump1090", "--oversample", "--net", "--device-index", strconv.Itoa(e.indexID), "--ppm", strconv.Itoa(e.ppm))
stdout, _ := cmd.StdoutPipe()
stderr, _ := cmd.StderrPipe()
@ -128,12 +128,12 @@ func (u *UAT) read() {
}
func getPPM(serial string) int {
r, err := regexp.Compile("str?a?t?u?x:\\d+:?(-?\\d*)");
r, err := regexp.Compile("str?a?t?u?x:\\d+:?(-?\\d*)")
if err != nil {
return globalSettings.PPM
}
arr := r.FindStringSubmatch(serial);
arr := r.FindStringSubmatch(serial)
if arr == nil {
return globalSettings.PPM
}
@ -403,7 +403,7 @@ func sdrWatcher() {
serial = ""
}
if (rES != nil) {
if rES != nil {
doSkip = rES.MatchString(serial)
} else {
doSkip = strings.Compare(serial, "stratux:1090") == 0
@ -452,7 +452,7 @@ func sdrWatcher() {
serial = ""
}
if (rUAT != nil) {
if rUAT != nil {
doSkip = rUAT.MatchString(serial)
} else {
doSkip = strings.Compare(serial, "stratux:978") == 0

Wyświetl plik

@ -65,11 +65,11 @@ const (
TARGET_TYPE_MODE_S = 0
TARGET_TYPE_ADSB = 1
TARGET_TYPE_ADSR = 2
TARGET_TYPE_TISB = 3
// Assign next type to UAT messages with address qualifier == 2
// (code corresponds to and UAT GBT targets with Mode S addresses..
// for we'll treat them as ADS-R in the UI. Might be able to assign as TYPE_ADSR if we see a proper emitter category and NIC > 7.
TARGET_TYPE_TISB_S = 4
// for we'll treat them as ADS-R in the UI. If we see a proper emitter category and NIC > 7, they'll be reassigned to TYPE_ADSR.
TARGET_TYPE_TISB_S = 3
TARGET_TYPE_TISB = 4
)
type TrafficInfo struct {
@ -79,7 +79,7 @@ type TrafficInfo struct {
OnGround bool // Air-ground status. On-ground is "true".
Addr_type uint8 // UAT address qualifier. Used by GDL90 format, so translations for ES TIS-B/ADS-R are needed.
TargetType uint8 // types decribed in const above
SignalLevel int32 // Arbitrary signal level reported by dump1090 and dump978.
SignalLevel float64 // Signal level, dB RSSI.
Position_valid bool // set when position report received. Unset after n seconds? (To-do)
Lat float32 // decimal degrees, north positive
Lng float32 // decimal degrees, east positive
@ -116,7 +116,7 @@ type dump1090Data struct {
TypeCode int // Mode S type code
SubtypeCode int // Mode S subtype code
SBS_MsgType int // type of SBS message (used in "old" 1090 parsing)
SignalLevel int // 0-255
SignalLevel float64 // Decimal RSSI (0-1 nominal) as reported by dump1090-mutability. Convert to dB RSSI before setting in TrafficInfo.
Tail *string
Squawk *int // 12-bit squawk code in octal format
Emitter_category *int
@ -152,14 +152,14 @@ func sendTrafficUpdates() {
defer trafficMutex.Unlock()
cleanupOldEntries()
var msg []byte
if ((stratuxClock.Time.Second() % 10) == 0) {
if (stratuxClock.Time.Second() % 30) == 0 {
log.Printf("List of all aircraft being tracked:\n")
log.Printf("==================================================================\n")
}
for icao, ti := range traffic { // TO-DO: Limit number of aircraft in traffic message. ForeFlight 7.5 chokes at ~1000-2000 messages depending on iDevice RAM. Practical limit likely around ~500 aircraft without filtering.
// DEBUG: Print the list of all tracked targets (with data) to the log every ten seconds
if ((stratuxClock.Time.Second() % 10) == 0) {
if (stratuxClock.Time.Second() % 10) == 0 {
s_out, err := json.Marshal(ti)
if err != nil {
log.Printf("Error generating output: %s\n", err.Error())
@ -266,7 +266,7 @@ func makeTrafficReportMsg(ti TrafficInfo) []byte {
}
// Position containment / navigational accuracy
msg[13] = ((byte(ti.NIC) << 4) & 0xF0) | ((byte(ti.NACp) & 0x0F))
msg[13] = ((byte(ti.NIC) << 4) & 0xF0) | (byte(ti.NACp) & 0x0F)
// Horizontal velocity (speed).
@ -296,7 +296,7 @@ func makeTrafficReportMsg(ti TrafficInfo) []byte {
return prepareMessage(msg)
}
func parseDownlinkReport(s string, signalLevel int32) {
func parseDownlinkReport(s string, signalLevel int) {
var ti TrafficInfo
s = s[1:]
@ -363,11 +363,15 @@ func parseDownlinkReport(s string, signalLevel int32) {
ti.NIC = int(frame[11] & 0x0F)
if ((msg_type == 1) || (msg_type == 3)) { // Since NACp is passed with normal UATreports, no need to use our ES hack.
if (msg_type == 1) || (msg_type == 3) { // Since NACp is passed with normal UATreports, no need to use our ES hack.
ti.NACp = int((frame[25] >> 4) & 0x0F)
}
ti.SignalLevel = signalLevel
power := 10 * (math.Log10(float64((signalLevel * signalLevel))) - 6) // reported amplitude is 0-1000. Subtract 60 (10*log 1e6) to set nominal full scale = 0.
//log.Printf("%s (%X) seen with amplitude of %d, corresponding to normalized power of %f.2 dB\n",ti.Tail,ti.Icao_addr,signalLevel,power)
ti.SignalLevel = power
if ti.Addr_type == 0 {
ti.TargetType = TARGET_TYPE_ADSB
@ -377,7 +381,7 @@ func parseDownlinkReport(s string, signalLevel int32) {
ti.TargetType = TARGET_TYPE_ADSR
} else if ti.Addr_type == 2 {
ti.TargetType = TARGET_TYPE_TISB_S
if ((ti.NIC >= 7) && (ti.Emitter_category > 0)) { // If NIC is sufficiently and emitter type is transmitted, we'll assume it's ADS-R.
if (ti.NIC >= 7) && (ti.Emitter_category > 0) { // If NIC is sufficiently and emitter type is transmitted, we'll assume it's ADS-R.
ti.TargetType = TARGET_TYPE_ADSR
}
}
@ -583,7 +587,7 @@ func esListen() {
ti.ExtrapolatedPosition = false
}
ti.SignalLevel = int32(newTi.SignalLevel)
ti.SignalLevel = 10 * math.Log10(newTi.SignalLevel)
// generate human readable summary of message types for debug
// TO-DO: Use for ES message statistics?
@ -692,7 +696,7 @@ func esListen() {
ti.Speed_valid = true
ti.Last_speed = stratuxClock.Time // only update "last seen" data on position updates
}
} else if (((newTi.DF == 17) || (newTi.DF == 18)) && (newTi.TypeCode == 19)) { // invalid speed on velocity message only
} else if ((newTi.DF == 17) || (newTi.DF == 18)) && (newTi.TypeCode == 19) { // invalid speed on velocity message only
ti.Speed_valid = false
}
@ -731,7 +735,7 @@ func esListen() {
}
ti.NIC = nic
if ((ti.NACp < 7) && (ti.NACp < ti.NIC)) {
if (ti.NACp < 7) && (ti.NACp < ti.NIC) {
ti.NACp = ti.NIC // initialize to NIC, since NIC is sent with every position report, and not all emitters report NACp.
}
}
@ -753,10 +757,10 @@ func esListen() {
if newTi.CA == 6 {
ti.TargetType = TARGET_TYPE_ADSR
ti.Addr_type = 2
} else if (newTi.CA == 2) { // 2 = TIS-B with ICAO address, 5 = TIS-B without ICAO address
} else if newTi.CA == 2 { // 2 = TIS-B with ICAO address, 5 = TIS-B without ICAO address
ti.TargetType = TARGET_TYPE_TISB
ti.Addr_type = 2
} else if (newTi.CA == 5) {
} else if newTi.CA == 5 {
ti.TargetType = TARGET_TYPE_TISB
ti.Addr_type = 3
}
@ -835,6 +839,17 @@ func updateDemoTraffic(icao uint32, tail string, relAlt float32, gs float64, off
if ti.Addr_type == 1 { //reserved value
ti.Addr_type = 0
}
if ti.Addr_type == 0 {
ti.TargetType = TARGET_TYPE_ADSB
}
if ti.Addr_type == 2 {
ti.TargetType = TARGET_TYPE_ADSR
}
if ti.Addr_type == 3 {
ti.TargetType = TARGET_TYPE_TISB
}
ti.Emitter_category = 1
ti.Lat = float32(lat + traffRelLat)
ti.Lng = float32(lng + traffRelLng)
@ -857,13 +872,12 @@ func updateDemoTraffic(icao uint32, tail string, relAlt float32, gs float64, off
ti.Last_seen = stratuxClock.Time
ti.Last_alt = stratuxClock.Time
ti.Last_speed = stratuxClock.Time
ti.TargetType = TARGET_TYPE_ADSB
ti.NACp = 8
ti.NIC = 8
//ti.Age = math.Floor(ti.Age) + hdg / 1000
ti.Last_source = 1
if icao%7 == 1 { // make some of the traffic look like it came from UAT
if icao%5 == 1 { // make some of the traffic look like it came from UAT
ti.Last_source = 2
}

Wyświetl plik

@ -160,7 +160,8 @@
}
.traffic-style12,
.traffic-style13 {
.traffic-style13,
.traffic-style14 {
color: #000000;
background-color: skyblue;
}
@ -174,7 +175,8 @@
}
.traffic-style22,
.traffic-style23 {
.traffic-style23,
.traffic-style24 {
color: #000000;
background-color: khaki
}

Wyświetl plik

@ -19,8 +19,9 @@ function TowersCtrl($rootScope, $scope, $state, $http, $interval) {
function setTower(obj, new_tower) {
new_tower.lat = dmsString(obj.Lat);
new_tower.lon = dmsString(obj.Lng);
new_tower.signal = obj.Signal_strength_last_minute;
// Signal_strength_max int
new_tower.power = obj.Signal_strength_now.toFixed(2);
new_tower.power_last_min = obj.Signal_strength_last_minute.toFixed(2);
new_tower.power_max = obj.Signal_strength_max.toFixed(2);
// Messages_last_minute uint64
new_tower.messages = obj.Messages_total;
}
@ -57,7 +58,7 @@ function TowersCtrl($rootScope, $scope, $state, $http, $interval) {
var updateTowers = $interval(function () {
// refresh tower list once every 5 seconds (aka polling)
getTowers();
}, (5 * 1000), 0, false);
}, (2 * 1000), 0, false);
$state.get('towers').onEnter = function () {
// everything gets handled correctly by the controller

Wyświetl plik

@ -35,10 +35,11 @@ function TrafficCtrl($rootScope, $scope, $state, $http, $interval) {
function setAircraft(obj, new_traffic) {
new_traffic.icao_int = obj.Icao_addr;
new_traffic.addr_type = obj.Addr_type;
new_traffic.targettype = obj.TargetType;
new_traffic.signal = obj.SignalLevel.toFixed(2);
new_traffic.addr_symb ='\u2708';
if (new_traffic.addr_type > 1) {
new_traffic.addr_symb ='\u{1F4E1}';
if (new_traffic.targettype > 2) {
new_traffic.addr_symb ='\ud83d\udce1';
}
new_traffic.icao = obj.Icao_addr.toString(16).toUpperCase();
new_traffic.tail = obj.Tail;

Wyświetl plik

@ -7,16 +7,20 @@
</div>
<div class="panel-body towers-page">
<div class="row">
<span class="col-xs-6"><strong>Location</strong></span>
<span class="col-xs-3 text-right">Signal</span>
<span class="col-xs-3 text-right">Msgs</span>
<span class="col-xs-4"><strong>Location</strong></span>
<span class="col-xs-2 text-right">Current Power (dB)</span>
<span class="col-xs-2 text-right">Avg Power Last Min(dB)</span>
<span class="col-xs-2 text-right">Max Power (dB)</span>
<span class="col-xs-2 text-right">Msgs Last Minute</span>
</div>
<div class="row" ng-repeat="tower in data_list | orderBy: messages">
<div class="separator"></div>
<span class="col-xs-6">{{tower.lat}} {{tower.lon}}</span>
<span class="col-xs-3 text-right">{{tower.signal}}</span>
<span class="col-xs-3 text-right">{{tower.messages}}</span>
<span class="col-xs-4">{{tower.lat}} {{tower.lon}}</span>
<span class="col-xs-2 text-right">{{tower.power}}</span>
<span class="col-xs-2 text-right">{{tower.power_last_min}}</span>
<span class="col-xs-2 text-right">{{tower.power_max}}</span>
<span class="col-xs-2 text-right">{{tower.messages}}</span>
</div>
</div>
</div>

Wyświetl plik

@ -15,7 +15,8 @@
</div>
<div class="col-sm-6">
<span class="col-xs-2">&nbsp;</span>
<span class="col-xs-7">Location</span>
<span class="col-xs-5">Location</span>
<span class="col-xs-2">Power (dB)</span>
<span class="col-xs-3 text-right">Age</span>
</div>
</div>
@ -24,8 +25,8 @@
<div class="separator"></div>
<div class="col-sm-6">
<span class="col-xs-3">
<span ng-show="aircraft.tail" ng-class="'label traffic-style'+aircraft.src+aircraft.addr_type"><strong>{{aircraft.addr_symb}}&nbsp;{{aircraft.tail}}</strong></span>
<span ng-hide="aircraft.tail" ng-class="'label traffic-style'+aircraft.src+aircraft.addr_type"><strong class="text-muted">{{aircraft.addr_symb}}&nbsp;{{aircraft.icao}}</strong></span>
<span ng-show="aircraft.tail" ng-class="'label traffic-style'+aircraft.src+aircraft.targettype"><strong>{{aircraft.addr_symb}}&nbsp;{{aircraft.tail}}</strong></span>
<span ng-hide="aircraft.tail" ng-class="'label traffic-style'+aircraft.src+aircraft.targettype"><strong class="text-muted">{{aircraft.addr_symb}}&nbsp;{{aircraft.icao}}</strong></span>
</span>
<span class="col-xs-3 text-right">{{aircraft.speed}} KTS</span>
@ -38,7 +39,8 @@
</div>
<div class="col-sm-6">
<span class="col-xs-2">&nbsp;</span>
<span class="col-xs-7">{{aircraft.lat}} {{aircraft.lon}}</span>
<span class="col-xs-5">{{aircraft.lat}} {{aircraft.lon}}</span>
<span class="col-xs-2">{{aircraft.signal}}</span>
<!--<span class="col-xs-3 text-right">{{aircraft.time}}</span>-->
<span class="col-xs-3 text-right">{{aircraft.age}}s</span>
</div>