Subversion Repositories Projects

#include <inttypes.h>

#include "ubx.h"

#include "timer0.h"

#include "uart0.h"

// ------------------------------------------------------------------------------------------------

// defines

#define DAYS_FROM_JAN01YEAR0001_TO_JAN6_1980 722819 // the year 0 does not exist!

#define DAYS_PER_YEAR           365

#define DAYS_PER_LEAPYEAR       366

#define DAYS_PER_4YEARS         1461    //((3 * DAYS_PER_YEAR) + DAYS_PER_LEAPYEAR) // years dividable by 4 are leap years

#define DAYS_PER_100YEARS       36524   //((25 * DAYS_PER_4YEARS) - 1) // years dividable by 100 are no leap years

#define DAYS_PER_400YEARS       146097  //((4 * DAYS_PER_100YEARS) + 1L) // but years dividable by 400 are leap years

#define SECONDS_PER_MINUTE      60

#define MINUTES_PER_HOUR        60

#define HOURS_PER_DAY           24

#define DAYS_PER_WEEK           7

#define SECONDS_PER_HOUR        3600    //(SECONDS_PER_MINUTE * MINUTES_PER_HOUR)

#define SECONDS_PER_DAY         86400   //(SECONDS_PER_HOUR * HOURS_PER_DAY)

#define SECONDS_PER_WEEK        604800  //(SECONDS_PER_DAY * DAYS_PER_WEEK)

// days per month in normal and leap years

const uint32_t   Leap[ 13 ] = { 0,  31,  60,  91, 121, 152, 182, 213, 244, 274, 305, 335, 366 };

const uint32_t Normal[ 13 ]     = { 0,  31,  59,  90, 120, 151, 181, 212, 243, 273, 304, 334, 365 };

#define LEAP_SECONDS_FROM_1980  15 // the last one was on the Dec 31th 2008

// message sync bytes

#define UBX_SYNC1_CHAR  0xB5

#define UBX_SYNC2_CHAR  0x62

// protocoll identifier

#define UBX_CLASS_NAV   0x01

// message id

#define UBX_ID_POSLLH   0x02

#define UBX_ID_SOL              0x06

#define UBX_ID_VELNED   0x12

// ------------------------------------------------------------------------------------------------

// typedefs

// ubx parser state

typedef enum

{

        UBXSTATE_IDLE,

        UBXSTATE_SYNC1,

        UBXSTATE_SYNC2,

        UBXSTATE_CLASS,

        UBXSTATE_LEN1,

        UBXSTATE_LEN2,

        UBXSTATE_DATA,

        UBXSTATE_CKA,

        UBXSTATE_CKB

} ubxState_t;

typedef struct

{

        uint32_t        itow;           // ms GPS Millisecond Time of Week

        int32_t         frac;           // ns remainder of rounded ms above

        int16_t         week;           // GPS week

        uint8_t         GPSfix;         // GPSfix Type, range 0..6

        uint8_t         Flags;          // Navigation Status Flags

        int32_t         ECEF_X;         // cm ECEF X coordinate

        int32_t         ECEF_Y;         // cm ECEF Y coordinate

        int32_t         ECEF_Z;         // cm ECEF Z coordinate

        int32_t         PAcc;           // cm 3D Position Accuracy Estimate

        int32_t         ECEFVX;         // cm/s ECEF X velocity

        int32_t         ECEFVY;         // cm/s ECEF Y velocity

        int32_t         ECEFVZ;         // cm/s ECEF Z velocity

        uint32_t        SAcc;           // cm/s Speed Accuracy Estimate

        uint16_t        PDOP;           // 0.01 Position DOP

        uint8_t         res1;           // reserved

        uint8_t         numSV;          // Number of SVs used in navigation solution

        uint32_t        res2;           // reserved

        uint8_t         Status;     // invalid/newdata/processed

} __attribute__((packed)) ubx_nav_sol_t;

typedef struct

{

        uint32_t        itow;           // ms  GPS Millisecond Time of Week

        int32_t         VEL_N;          // cm/s  NED north velocity

        int32_t         VEL_E;          // cm/s  NED east velocity

        int32_t         VEL_D;          // cm/s  NED down velocity

        int32_t         Speed;          // cm/s  Speed (3-D)

        int32_t         GSpeed;         // cm/s  Ground Speed (2-D)

        int32_t         Heading;        // 1e-05 deg  Heading 2-D

        uint32_t        SAcc;           // cm/s  Speed Accuracy Estimate

        uint32_t        CAcc;           // deg  Course / Heading Accuracy Estimate

        uint8_t         Status;         // invalid/newdata/processed

} __attribute__((packed)) ubx_nav_velned_t;

typedef struct

{

        uint32_t        itow;           // ms GPS Millisecond Time of Week

        int32_t         LON;            // 1e-07 deg Longitude

        int32_t         LAT;            // 1e-07 deg Latitude

        int32_t         HEIGHT;         // mm Height above Ellipsoid

        int32_t         HMSL;           // mm Height above mean sea level

        uint32_t        Hacc;           // mm Horizontal Accuracy Estimate

        uint32_t        Vacc;           // mm Vertical Accuracy Estimate

        uint8_t         Status;         // invalid/newdata/processed

} __attribute__((packed)) ubx_nav_posllh_t;

//------------------------------------------------------------------------------------

// global variables

// local buffers for the incomming ubx messages

volatile ubx_nav_sol_t          UbxSol    = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, INVALID};

volatile ubx_nav_posllh_t       UbxPosLlh = {0,0,0,0,0,0,0, INVALID};

volatile ubx_nav_velned_t       UbxVelNed = {0,0,0,0,0,0,0,0,0, INVALID};

uint16_t CheckGPSOkay = 0;

// shared buffer

gps_data_t              GPSData = {{0,0,0,INVALID},0,0,0,0,0,0,0, INVALID};

//------------------------------------------------------------------------------------

// functions

uint8_t IsLeapYear(uint16_t year)

{

        if((year%400 == 0) || ( (year%4 == 0) && (year%100 != 0) ) ) return 1;

        else return 0;

}

/********************************************************/

/*  Calculates the UTC Time from the GPS week and tow   */

/********************************************************/

void SetGPSTime(DateTime_t * pTimeStruct)

{

        uint32_t Days, Seconds, Week;

        uint16_t YearPart;

        uint32_t * MonthDayTab = 0;

        uint8_t  i;

        // if GPS data show valid time data

        if((UbxSol.Status != INVALID) && (UbxSol.Flags & FLAG_WKNSET) && (UbxSol.Flags & FLAG_TOWSET) )

        {

                Seconds = UbxSol.itow / 1000L;

                Week = (uint32_t)UbxSol.week;

                // correct leap seconds since 1980

                if(Seconds < LEAP_SECONDS_FROM_1980)

                {

                        Week--;

                        Seconds = SECONDS_PER_WEEK - LEAP_SECONDS_FROM_1980 + Seconds;

                }

                else Seconds -= LEAP_SECONDS_FROM_1980;

                Days = DAYS_FROM_JAN01YEAR0001_TO_JAN6_1980;

                Days += (Week * DAYS_PER_WEEK);

                Days += Seconds / SECONDS_PER_DAY; // seperate days from GPS seconds of week

                pTimeStruct->Year = 1;

                YearPart = (uint16_t)(Days / DAYS_PER_400YEARS);

                pTimeStruct->Year += YearPart * 400;

                Days = Days % DAYS_PER_400YEARS;

                YearPart = (uint16_t)(Days / DAYS_PER_100YEARS);

                pTimeStruct->Year += YearPart * 100;

                Days = Days % DAYS_PER_100YEARS;

                YearPart = (uint16_t)(Days / DAYS_PER_4YEARS);

                pTimeStruct->Year += YearPart * 4;

                Days = Days % DAYS_PER_4YEARS;

                if(Days < (3* DAYS_PER_YEAR)) YearPart = (uint16_t)(Days / DAYS_PER_YEAR);

                else YearPart = 3;

                pTimeStruct->Year += YearPart;

                // calculate remaining days of year

                Days -= (uint32_t)(YearPart *  DAYS_PER_YEAR);

                Days += 1;

                // check if current year is a leap year

                if(IsLeapYear(pTimeStruct->Year)) MonthDayTab = (uint32_t*)Leap;

                else MonthDayTab = (uint32_t*)Normal;

            // seperate month and day from days of year

                for ( i = 0; i < 12; i++ )

                {

                        if ( (MonthDayTab[i]< Days) && (Days <= MonthDayTab[i+1]) )

                        {

                                pTimeStruct->Month = i+1;

                                pTimeStruct->Day = Days - MonthDayTab[i];

                                i = 12;

                        }

                }

                Seconds = Seconds % SECONDS_PER_DAY; // remaining seconds of current day

                pTimeStruct->Hour = (uint8_t)(Seconds / SECONDS_PER_HOUR);

                Seconds = Seconds % SECONDS_PER_HOUR; // remaining seconds of current hour

                pTimeStruct->Min = (uint8_t)(Seconds / SECONDS_PER_MINUTE);

                Seconds = Seconds % SECONDS_PER_MINUTE; // remaining seconds of current minute

                pTimeStruct->Sec = (uint8_t)(Seconds);

                pTimeStruct->mSec  = (uint16_t)(UbxSol.itow % 1000L);

                pTimeStruct->Valid = 1;

        }

        else

        {

                pTimeStruct->Valid = 0;

        }

}

/********************************************************/

/*                  Initialize UBX Parser               */

/********************************************************/

void UBX_Init(void)

{

        // mark msg buffers invalid

        UbxSol.Status = INVALID;

        UbxPosLlh.Status = INVALID;

        UbxVelNed.Status = INVALID;

        GPSData.Status = INVALID;

}

/********************************************************/

/*            Upate GPS data stcructure                 */

/********************************************************/

void Update_GPSData (void)

{

        static uint16_t Ubx_Timeout = 0;

        static uint8_t  Msg_Count = 0;

        // the timeout is used to detect the delay between two message sets

        // and is used for synchronisation so that always a set is collected

        // that belongs together

        // _______NAVSOL|POSLLH|VELNED|___________________NAVSOL|POSLLH|VELNED|_____________

        //              |  8ms | 8ms  |         184 ms          |      |      |

        // msg_count:   0      1      2                         0      1      2

        if(CheckDelay(Ubx_Timeout))     Msg_Count = 0;

        else Msg_Count++;

        Ubx_Timeout = SetDelay(100); // reset ubx msg timeout

        // if a new set of ubx messages was collected

        if((Msg_Count >= 2))

        {       // if set is complete

                if((UbxSol.Status == NEWDATA) && (UbxPosLlh.Status == NEWDATA) && (UbxVelNed.Status == NEWDATA))

                {

                        CheckGPSOkay++;

                        // update GPS data only if the status is INVALID or PROCESSED  and the last ubx message was received within less than 100 ms

                        if(GPSData.Status != NEWDATA) // if last data were processed

                        { // wait for new data at all neccesary ubx messages

                                GPSData.Status = INVALID;

                                // NAV SOL

                                GPSData.Flags =                                 UbxSol.Flags;

                                GPSData.NumOfSats =                     UbxSol.numSV;

                                GPSData.SatFix =                                UbxSol.GPSfix;

                                GPSData.Position_Accuracy =             UbxSol.PAcc;

                                GPSData.Speed_Accuracy =                UbxSol.SAcc;

                                SetGPSTime(&SystemTime); // update system time

                                // NAV POSLLH

                                GPSData.Position.Status =               INVALID;

                                GPSData.Position.Longitude =    UbxPosLlh.LON;

                                GPSData.Position.Latitude =     UbxPosLlh.LAT;

                                GPSData.Position.Altitude =     UbxPosLlh.HMSL;

                                GPSData.Position.Status =               NEWDATA;

                                // NAV VELNED

                                GPSData.Speed_East =                    UbxVelNed.VEL_E;

                                GPSData.Speed_North =                   UbxVelNed.VEL_N;

                                GPSData.Speed_Top       =                       -UbxVelNed.VEL_D;

                                GPSData.Speed_Ground =                  UbxVelNed.GSpeed;

                                GPSData.Heading =                               UbxVelNed.Heading;

                                GPSData.Status = NEWDATA; // new data available

                        } // EOF if(GPSData.Status != NEWDATA)

                } // EOF all ubx messages received

                // set state to collect new data

                UbxSol.Status =                                 PROCESSED;      // ready for new data

                UbxPosLlh.Status =                              PROCESSED;      // ready for new data

                UbxVelNed.Status =                              PROCESSED;      // ready for new data

        }

}

/********************************************************/

/*                   UBX Parser                         */

/********************************************************/

void UBX_Parser(uint8_t c)

{

        static ubxState_t ubxState = UBXSTATE_IDLE;

        static uint16_t msglen;

        static uint8_t cka, ckb;

        static uint8_t *ubxP, *ubxEp, *ubxSp; // pointers to data currently transfered

        //state machine

        switch (ubxState)       // ubx message parser

        {

                case UBXSTATE_IDLE: // check 1st sync byte

                        if (c == UBX_SYNC1_CHAR) ubxState = UBXSTATE_SYNC1;

                        else ubxState = UBXSTATE_IDLE; // out of synchronization

                        break;

                case UBXSTATE_SYNC1: // check 2nd sync byte

                        if (c == UBX_SYNC2_CHAR) ubxState = UBXSTATE_SYNC2;

                        else ubxState = UBXSTATE_IDLE; // out of synchronization

                        break;

                case UBXSTATE_SYNC2: // check msg class to be NAV

                        if (c == UBX_CLASS_NAV) ubxState = UBXSTATE_CLASS;

                        else ubxState = UBXSTATE_IDLE; // unsupported message class

                        break;

                case UBXSTATE_CLASS: // check message identifier

                        switch(c)

                        {

                                case UBX_ID_POSLLH: // geodetic position

                                        ubxP =  (uint8_t *)&UbxPosLlh; // data start pointer

                                        ubxEp = (uint8_t *)(&UbxPosLlh + 1); // data end pointer

                                        ubxSp = (uint8_t *)&UbxPosLlh.Status; // status pointer

                                        break;

                                case UBX_ID_SOL: // navigation solution

                                        ubxP =  (uint8_t *)&UbxSol; // data start pointer

                                        ubxEp = (uint8_t *)(&UbxSol + 1); // data end pointer

                                        ubxSp = (uint8_t *)&UbxSol.Status; // status pointer

                                        break;

                                case UBX_ID_VELNED: // velocity vector in tangent plane

                                        ubxP =  (uint8_t *)&UbxVelNed; // data start pointer

                                        ubxEp = (uint8_t *)(&UbxVelNed + 1); // data end pointer

                                        ubxSp = (uint8_t *)&UbxVelNed.Status; // status pointer

                                        break;

                                default:                        // unsupported identifier

                                        ubxState = UBXSTATE_IDLE;

                                        break;

                        }

                        if (ubxState != UBXSTATE_IDLE)

                        {

                                ubxState = UBXSTATE_LEN1;

                                cka = UBX_CLASS_NAV + c;

                                ckb = UBX_CLASS_NAV + cka;

                        }

                        break;

                case UBXSTATE_LEN1: // 1st message length byte

                        msglen = (uint16_t)c; // lowbyte first

                        cka += c;

                        ckb += cka;

                        ubxState = UBXSTATE_LEN2;

                        break;

                case UBXSTATE_LEN2: // 2nd message length byte

                        msglen += ((uint16_t)c)<<8; // high byte last

                        cka += c;

                        ckb += cka;

                        // if the old data are not processed so far then break parsing now

                        // to avoid writing new data in ISR during reading by another function

                        if ( *ubxSp == NEWDATA )

                        {

                                ubxState = UBXSTATE_IDLE;

                                Update_GPSData(); //update GPS info respectively

                        }

                        else // data invalid or allready processd

                        {

                                *ubxSp = INVALID; // mark invalid during buffer filling

                                ubxState = UBXSTATE_DATA;

                        }

                        break;

                case UBXSTATE_DATA: // collecting data

                        if (ubxP < ubxEp)

                        {

                                *ubxP++ = c; // copy curent data byte if any space is left

                                cka += c;

                                ckb += cka;

                                if (--msglen == 0)      ubxState = UBXSTATE_CKA; // switch to next state if all data was read

                        }

                        else // rx buffer overrun

                        {

                                ubxState = UBXSTATE_IDLE;

                        }

                        break;

                case UBXSTATE_CKA:

                        if (c == cka) ubxState = UBXSTATE_CKB;

                        else

                        {

                                *ubxSp = INVALID;

                                ubxState = UBXSTATE_IDLE;

                        }

                        break;

                case UBXSTATE_CKB:

                        if (c == ckb)

                        {

                                *ubxSp = NEWDATA; // new data are valid

                                Update_GPSData(); //update GPS info respectively

                        }

                        else

                        {       // if checksum not match then set data invalid

                                *ubxSp = INVALID;

                        }

                        ubxState = UBXSTATE_IDLE; // ready to parse new data

                        break;

                default: // unknown ubx state

                        ubxState = UBXSTATE_IDLE;

                        break;

        }

        DebugOut.Analog[9] = GPSData.Speed_North;

        DebugOut.Analog[10] = GPSData.Speed_East;

        DebugOut.Analog[11] = GPSData.Speed_Top;

        DebugOut.Analog[12] = GPSData.NumOfSats;

        DebugOut.Analog[13] = GPSData.Position.Longitude;

        DebugOut.Analog[14] = GPSData.Position.Latitude;

        DebugOut.Analog[15] = GPSData.Position.Altitude;

}