Projet

Général

Profil

5 jalaffon
#include <stdio.h>
#include <stdlib.h>
10 albest
#include <strings.h>
#include <math.h>
5 jalaffon
#include "trame.h"

//Trames de tests ? modifier si n?cessaire.
char * trames[]= {"$GPGSV,3,2,10,15,03,077,,18,04,041,42,19,85,271,,20,08,214,*7C",
"$GPGSV,3,3,10,22,39,053,50,28,15,320,*7E",
"$GPRMC,141914.00,A,4545.6424,N,00306.6036,E,0.4,99.4,010206,,*0C",
"$GPGLL,4545.6424,N,00306.6036,E,141914.00,A*0E",
"$GPGGA,141914.00,4545.0000,N,00306.6036,E,1,05,3.4,499.3,M,,M,,*7D",
"$GPGSA,A,3,,03,,22,14,,01,,18,,,,3.9,3.4,1.9*39",
"$GPVTG,99.4,T,,M,0.4,N,0.7,K*57",
"$GPZDA,141914.00,01,02,2006,00,00*69",
10 albest
0};
5 jalaffon
233 albest
int trame_cmp(char* trame,char* type){ //Cette fonction renvoie 1 si la trame est de la forme "GPGGA".
10 albest
int i = 0;
while (type[i]!= '\0'){
if (trame[i+1] != type[i]){
return 0;
}
i++;
}
return 1;
5 jalaffon
}
10 albest
233 albest
int decode_int(char c){ //Cette fonction renvoie la valeur decimale associ?e ? un caract?re donn? en param?tre compris entre '0' et '9' et renvoie -1 sinon.
10 albest
int val = c - '0';
580 albest
if (val<0 || val>9){ //**En premier faire le test pour ne pas faire op?ration pour rien
10 albest
val = -1;
}
return val;
}

233 albest
int decode_nombre(char * ch,int n){ //Renvoie la valeur d?cimale des n premiers caract?res de la cha?ne ch.
10 albest
int res = 0;
for (int i =0;i < n;i++){
214 albest
res *= 10;
res += decode_int(ch[i]);
10 albest
}
return res;
}

233 albest
float conv_lat(char * latitude){ //La fonction convertie la latitude en sexa en degr?.
580 albest
float val = decode_nombre(latitude,4); //* bof
10 albest
float res = 0;
17 albest
for (int i = 5; i < 9;i++){
233 albest
val += (decode_int(latitude[i]))*pow(10,-i+4); //On replace chaque chiffre au bonne endroit (unit?, dizaine ...).
17 albest
}
580 albest
res = (int)(val/100); //****
17 albest
res += (val - res*100)/60;
580 albest
return res;
}
10 albest
233 albest
float conv_long(char * longitude){ //La fonction convertie la longitude en sexa en degr?.
580 albest
float val = decode_nombre(longitude,5); //**
17 albest
float res = 0;
for (int i = 6; i < 10;i++){
val += (decode_int(longitude[i]))*pow(10,-i+5);
}
res = (int)(val/100);
res += (val - res*100)/60;
return res;}

233 albest
float conversion(char * chaine){ //Cette fonction convertie soit une latitude soit une longitude en degr?.
17 albest
int cpt = 0;
580 albest
while (chaine[cpt]!='\0'){ //***, .
17 albest
cpt++;
}
if (cpt == 9){
return conv_lat(chaine);
}
else if(cpt==10){
return conv_long(chaine);
}
else{
25 albest
return 1000.0;
17 albest
}
}

580 albest
typedef struct { //*
214 albest
float latitude;
float longitude;
} Position;

233 albest
int decode_trame(char * trame, Position *p){
580 albest
char longi[10]={'\0'}; //*** trame -> @ lati -> @
233 albest
char lati[9];
214 albest
int cpt = 0;
int i = 0;
int j = 0;
int y = 0;
if (trame_cmp(trame,"GPGGA")==1){
while (trame[i] != '\0'){
233 albest
if (trame[i] == ','){ //On compte le nombre de virgule.
214 albest
cpt += 1;
580 albest
i += 1;
214 albest
}
233 albest
if (cpt == 2) { //On regarde si le nombre de virgule correspond ? 2, ie la latitude
lati[j] = trame[i]; //On stocke caract?re par caract?re pour isoler la latitude afin de la convertir.
214 albest
j ++;
}
233 albest
if (cpt == 4) { //On regarde si le nombre de virgule correspond ? 4, ie la longitude
580 albest
longi[y] = trame[i];
214 albest
y ++;
}
580 albest
i += 1;
214 albest
}
580 albest
p->latitude = conv_lat(lati); //***************** vous avez de la chance &trame[15], trame+15

233 albest
p->longitude = conv_long(longi); //On stocke la latitude et la longitude convertie dans structure position.
return 1;
214 albest
}
233 albest
return 0;
25 albest
}
233 albest
#define PI 3.14159
float calcule_distance(Position p_1, Position p_2){
float rayon = 6378.14; //Rayon de la Terre (km).
float DistAng = 0;
DistAng = (PI/180)*acos((sin(p_1.latitude)*sin(p_2.latitude)+cos(p_1.latitude)*cos(p_2.latitude)*cos(p_2.longitude-p_1.longitude))); //Calcul distance angulaire
return DistAng*rayon; //On renvoi le calcul de la distance : D (km) = DistanceAngulaire*R.
}

float calcule_vitesse(Position p_1, Position p_2){
580 albest
return calcule_distance(p_1,p_2)*3600; //On calcul la vitesse (en km) avec la formule v=d/t avec t=1s.
233 albest
}

typedef struct{
Position rpos;
float vitmax;
} Zone;

Zone zones[] = {
{{44.7887762, -3.012}, 50},
{{44.7891220, -3.013}, 70},
};

int distance_a_la_plus_poche_zone(Position p, Zone r[], int nb_zones, float *d){
int index = 0;
580 albest
*d = calcule_distance(p, r[0].rpos); //On initialise une distance pour pouvoir ensuite la comparer et chercher le minimum.
233 albest
for (int i = 1; i<nb_zones; i++){ //On parcours la table de zone dangereuse.
580 albest
if (calcule_distance(p,r[i].rpos) < *d){
233 albest
index = i;
580 albest
*d = calcule_distance(p,r[i].rpos);
233 albest
}
}
return index;
}

573 albest
Position pos_prec;
580 albest
10 albest
void traitement(char * trame){
573 albest
Position pos;
295 albest
int alarme;
233 albest
float vitesse;
295 albest
float seuil = 0.5;
int index;
573 albest
float distance;

233 albest
if (trame_cmp(trame,"GPGGA")==1){
295 albest
printf ("> %s\n",trame);
573 albest
if (decode_trame(trame,&pos)==1){
vitesse = calcule_vitesse(pos,pos_prec);
index = distance_a_la_plus_poche_zone(pos,zones,2,&distance);
if ((distance <= seuil) && (vitesse > zones[index].vitmax)){
295 albest
alarme = 1;
573 albest
printf("Alarme on \n");
233 albest
}
295 albest
alarme = 0;
573 albest
printf("Alarme off \n");
233 albest
}
573 albest
pos_prec.latitude = pos.latitude;
pos_prec.longitude = pos.longitude;
295 albest
}
}
5 jalaffon
295 albest
5 jalaffon
//Ajouter vos tests unitaires dans cette fonction.
void tests_unitaires(void){
if (5!=5){
printf ("Erreur Test unitaire basique.\n");
exit(-1);
}
if (trame_cmp("$GPGGA suite chaine","GPGGA")!=1){
printf ("Erreur Test unitaire trame_cmp.\n");
exit(-1);
}
if (trame_cmp("$GPRMC suite chaine","GPGGA")!=0){
printf ("Erreur Test unitaire trame_cmp.\n");
exit(-1);
}
if (trame_cmp("$GPRMC... ", "GPRMC" )!=1){
printf ("Erreur Test unitaire trame_cmp.\n");
exit(-1);
}
if (trame_cmp("$APRMC...", "GPGGA")!=0){
printf ("Erreur Test unitaire trame_cmp.\n");
exit(-1);
10 albest
}
if (decode_int('5')!=5){
printf ("Erreur Test unitaire decode_int.\n");
exit(-1);
}
if (decode_int('A')!=-1){
printf ("Erreur Test unitaire decode_int.\n");
exit(-1);
}
if (decode_nombre("123",3)!=123){
printf ("Erreur Test unitaire decode_nombre.\n");
exit(-1);
}
if (decode_nombre("987654321",2)!=98){
printf ("Erreur Test unitaire decode_nombre.\n");
exit(-1);
}
17 albest
if (fabs(conv_lat("3723.2475")-37.387458)>= pow(10,-6)){
10 albest
printf ("Erreur Test unitaire conv_lat.\n");
exit(-1);
17 albest
}
if (fabs(conversion("00306.6036")-3.11006)>= pow(10,-6)){
printf ("Erreur Test unitaire conversion.\n");
exit(-1);
}
if (fabs(conversion("3723.2475")-37.387458)>= pow(10,-6)){
printf ("Erreur Test unitaire conversion.\n");
exit(-1);
233 albest
}
//Position pos1,pos2;
//decode_trame("$GPGGA,141914.00,4545.0000,N,00306.6036,E,1,05,3.4,499.3,M,,M,,*7D", &pos1);
//decode_trame("$GPGGA,141914.00,4545.0242,N,00306.6039,E,1,05,3.4,499.3,M,,M,,*7D", &pos2);
//printf("%f\n",calcule_distance(pos1, pos2));
//printf("%f\n",calcule_vitesse(pos1, pos2));
//float distance;
//printf("%d\n",distance_a_la_plus_poche_zone(pos1 , zones, 2, &distance));
//printf("%f\n", distance);
5 jalaffon
}

// Ne pas modifier cette fonction
int main(int argc,char ** argv)
214 albest
{
17 albest
tests_unitaires();
5 jalaffon
// Affichage des trames definies dans la table trames.
printf ("Trames de tests tableau trames:\n");
int i=0;
while (trames[i])
traitement(trames[i++]);

if (!trame_init())
exit(-1);
// Affichage des trames du fichier gps.log
char *trame;
printf ("Trames de tests du fichier gps.log\n");
while ((trame = trame_suivante()))
traitement(trame);

return 0;
}