root/branch/CELIK/sp4a12/main.c @ 496
| 494 | racelik | #include <stdio.h>
|
|
#include <stdlib.h>
|
|||
#include <strings.h>
|
|||
| 317 | racelik | #include "trame.h"
|
|
| 494 | racelik | #include<math.h>
|
|
| 317 | racelik | ||
| 494 | racelik | /** definition des constantes **/
|
|
#define rayon_terre 6370
|
|||
#define PI 3.14159265
|
|||
| 317 | racelik | ||
| 494 | racelik | /** definition des structures **/
|
|
| 317 | racelik | typedef struct{
|
|
| 494 | racelik | float latitude;
|
|
float longitude;
|
|||
} position;
|
|||
| 161 | racelik | ||
| 494 | racelik | typedef struct{
|
|
position rpos;
|
|||
float vitmax;
|
|||
} zone;
|
|||
| 137 | flgamel | ||
| 494 | racelik | /*Trames de tests ? modifier si n?cessaire.*/ /** tableau des differentes trames **/
|
|
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",
|
|||
0};
|
|||
/** description de toutes les zones **/
|
|||
zone zones[]= {
|
|||
{{44.7887762, -3.012}, 50}, /* Descripteur de la premiere zone */
|
|||
{{44.7891220, -3,013}, 70},
|
|||
{{45.7891220, 3,013}, 90},
|
|||
};
|
|||
position pos, pos_prec;
|
|||
| 161 | racelik | int trame_cmp(char *trame, char *type)
|
|
| 137 | flgamel | {
|
|
| 494 | racelik | int i; /* variable qui permet de parcourir le tableau de caractere trame et type */
|
|
| 317 | racelik | int test=1; /* test=1, si les 2 trames en question sont les memes */
|
|
| 161 | racelik | ||
| 494 | racelik | for (i=0; i<4; i++) /* on parcourt 4 fois car il y a 4 lettres sur les debuts de trame en prenant pas en compte le caractere $ du debut de trame */
|
|
| 137 | flgamel | {
|
|
| 317 | racelik | if (trame[i+1] != type[i]) /* verifie que le debut de la trame correspond au caractere type */
|
|
| 137 | flgamel | {
|
|
| 494 | racelik | test=0; /* test=0, si les 2 trames sont differentes */
|
|
| 137 | flgamel | }
|
|
}
|
|||
| 161 | racelik | return test;
|
|
| 137 | flgamel | }
|
|
| 494 | racelik | ||
| 317 | racelik | int decode_int(char c) /* fonction qui permet de passer d'un caractere a un entier */
|
|
| 161 | racelik | {
|
|
int res;
|
|||
| 137 | flgamel | ||
| 494 | racelik | if ((c>=48) && (c<=58)) /* comparaison si l'argument en entree est compris entre 0 et 9 */
|
|
| 161 | racelik | {
|
|
| 494 | racelik | res=c-48; /* on soustrait 48 car on veut afficher la bonne valeur */
|
|
| 161 | racelik | }
|
|
else
|
|||
{
|
|||
| 494 | racelik | res=-1; /* sinon le res vaut -1 s'il n'est pas compris entre 0 et 9 */
|
|
| 161 | racelik | }
|
|
return res;
|
|||
}
|
|||
| 168 | racelik | ||
| 494 | racelik | int decode_nombre(char *ch, int i, int n) /* fonction qui permet de passer d'une chaine de caractere a un entier */
|
|
| 168 | racelik | {
|
|
| 494 | racelik | int res=0; /* initialisation de res=0 */
|
|
int j; /* iteration du tableau de caractere ch */
|
|||
| 168 | racelik | ||
| 494 | racelik | for (j=i; j<i+n; j++) /* parcourt la chaine jusqu'au nombre de caractere souhaite */
|
|
| 168 | racelik | {
|
|
| 494 | racelik | res=(10*res) + decode_int(ch[j]); /* permet de reconstituer le nombre terme a terme par puissance de 10 */
|
|
| 168 | racelik | }
|
|
| 494 | racelik | ||
| 317 | racelik | return res;
|
|
| 168 | racelik | }
|
|
| 280 | racelik | ||
| 494 | racelik | float decode_latitude(char *ch) /* fonction convertissant la cha?ne de caract?re contenant la latitude en nombre flottant */
|
|
| 280 | racelik | {
|
|
| 494 | racelik | float degres, minute, flottant, res;
|
|
| 280 | racelik | ||
| 494 | racelik | degres = decode_nombre(ch, 17, 2); /* dans la trame GPGGA, retrouver les degres en sexagesimal */
|
|
minute = decode_nombre(ch, 19, 2); /* dans la trame GPGGA, retrouver les minutes en sexagesimal */
|
|||
flottant = decode_nombre(ch, 22, 4); /* dans la trame GPGGA, retrouver les secondes (flottant) en sexagesimal */
|
|||
| 280 | racelik | ||
| 494 | racelik | res = degres + (minute+flottant/10000)/60; /* conversion sexagesimal => degres */
|
|
| 317 | racelik | ||
| 494 | racelik | if (ch[27]=='S') /* identifier dans la trame GPGGA, si la latitude est nord ou bien le sud */
|
|
{
|
|||
res=-res; /* Si c'est au sud alors inverser les resultats, inversion de polarite */
|
|||
}
|
|||
return res;
|
|||
| 280 | racelik | }
|
|
| 494 | racelik | float decode_longitude(char *ch) /* fonction convertissant la cha?ne de caract?re contenant la longitude en nombre flottant, meme principe que decode_latitude */
|
|
| 280 | racelik | {
|
|
| 494 | racelik | float degres, minute, flottant, res;
|
|
| 280 | racelik | ||
| 494 | racelik | degres = decode_nombre(ch, 29, 3); /* dans la trame GPGGA, retrouver les degres en sexagesimal */
|
|
minute = decode_nombre(ch, 32, 2); /* dans la trame GPGGA, retrouver les minutes en sexagesimal */
|
|||
flottant = decode_nombre(ch, 35, 4); /* dans la trame GPGGA, retrouver les secondes (flottant) en sexagesimal */
|
|||
res = degres + (minute+flottant/10000)/60; /* conversion sexagesimal => degres */
|
|||
if (ch[40]=='W') /* identifier dans la trame GPGGA, si la longitude est ouest ou bien est */
|
|||
| 280 | racelik | {
|
|
| 494 | racelik | res=-res; /* Si c'est ouest alors inverser les resultats, inversion de polarite */
|
|
| 280 | racelik | }
|
|
| 494 | racelik | return res;
|
|
}
|
|||
| 280 | racelik | ||
| 494 | racelik | void decode_trame(char *ch, position *p) /* fonction qui extrait la position (latitude et longitude) pour les trames valides */
|
|
{
|
|||
/* recuperer la latitude ainsi que la longitude de la trame en question */
|
|||
p->latitude = decode_latitude(ch);
|
|||
p->longitude = decode_longitude(ch);
|
|||
| 280 | racelik | }
|
|
| 494 | racelik | float calcul_distance(position *p1, position *p2)/* fonction qui calcule la distance (en km) entre 2 points reperes par leur latitude et longitude */
|
|
| 280 | racelik | {
|
|
| 494 | racelik | float res_distance;
|
|
| 280 | racelik | ||
| 496 | racelik | /* formule de calcul d'une distance sur une sphere, d'apres les recherches sur internet */
|
|
res_distance = sqrt(pow((p2->latitude-p1->latitude)*PI*rayon_terre/180.0,2)+pow(((p2->longitude-p1->longitude)*(cos(((p2->latitude+p1->latitude)/2.0)*PI/180.0)*PI*rayon_terre/180.0)),2));
|
|||
| 494 | racelik | return res_distance;
|
|
| 280 | racelik | }
|
|
| 317 | racelik | ||
| 494 | racelik | float calcul_vitesse(position p1, position p2) /* fonction qui calcule la vitesse moyenne (en km/h) entre 2 releves GPS successifs */
|
|
| 317 | racelik | {
|
|
| 494 | racelik | float res_vitesse;
|
|
| 496 | racelik | res_vitesse = calcul_distance(&p1, &p2)*3600.0; /* v=d/t sachant que le tps est de 1s on a donc la vitesse en km/s, soit x3600 pour l'avoir en km/h */
|
|
| 494 | racelik | return res_vitesse;
|
|
}
|
|||
| 317 | racelik | ||
| 494 | racelik | int distance_a_la_plus_proche_zone(position p, zone r[], int nb_zones, float *d) /* fonction qui retourne le numero de la zone dont la distance est la plus proche */
|
|
{
|
|||
float reference;
|
|||
int i, j; /* indice j */
|
|||
| 317 | racelik | ||
| 494 | racelik | reference = calcul_distance(&p, &r[0].rpos); /* recupere la distance entre la position p et la position initiale de la zone dangereuse */
|
|
| 317 | racelik | ||
| 494 | racelik | for (i=0; i<nb_zones; i++) /* parcourt des zones dangereuses */
|
|
{
|
|||
*d = calcul_distance(&p, &r[i].rpos); /* recalcule la distance de la zone actuelle avec celle de la position */
|
|||
if (*d<=reference) /* comparaison entre la distance et la premiere zone dangereuse */
|
|||
{
|
|||
j=i; /* changement de l'indice de la zone la plus proche de la position courante */
|
|||
reference=*d; /* la reference est incremente pour recuperer le numero de la zone dangereuse */
|
|||
}
|
|||
}
|
|||
return j;
|
|||
| 317 | racelik | }
|
|
| 494 | racelik | //Fonction ? modifier !!!!!
|
|
void traitement(char * trame)
|
|||
| 317 | racelik | {
|
|
| 494 | racelik | /** initialisation des variables locales **/
|
|
static int cpt=0;
|
|||
float d;
|
|||
float v;
|
|||
int nbzone, trame_ok;
|
|||
int alarme=0;
|
|||
int seuil_d=100; /* distance max pour que l'alarme soit active */
|
|||
| 317 | racelik | ||
| 494 | racelik | trame_ok=1;
|
|
cpt++;
|
|||
| 317 | racelik | ||
| 494 | racelik | if (trame_cmp(trame,"GPGGA") == 1) /* test si les 2 trames sont equivalentes */
|
|
{
|
|||
printf ("> %s\n", trame); /* affiche uniquement les trames GPGGA */
|
|||
if (trame_ok==1) /* trame disponible */
|
|||
{
|
|||
trame_ok=0;
|
|||
decode_trame(trame, &pos); /* decoder la trame qui arrive pour obtenir la latitude et longitude */
|
|||
if (&pos) /* trame valide donc elle existe */
|
|||
{
|
|||
| 496 | racelik | v = calcul_vitesse(pos, pos_prec); /* calcule la vitesse entre la position de la voiture */
|
|
| 494 | racelik | nbzone=distance_a_la_plus_proche_zone(pos, zones, sizeof(zone), &d); /*definit la zone la plus proche de la position actuelle */
|
|
if ( (d<seuil_d) && (v>zones[nbzone].vitmax)) /*verification de la voiture dans une zone proche de la zone et si la vitesse est trop grande pour cette zone*/
|
|||
{
|
|||
alarme=1; /* alarme est active ==> ON */
|
|||
}
|
|||
if(alarme==1)
|
|||
{
|
|||
printf("alarme ON\n");
|
|||
}
|
|||
else
|
|||
{
|
|||
printf("alarme OFF\n");
|
|||
}
|
|||
pos_prec=pos; /* la position pos donnee par la trame traitee devient la position pos_prec precedente pour le calcul de vitesse et distance */
|
|||
}
|
|||
printf("Vitesse : %.3f Distance : %.3f\n\n", v, d); /* afficher la vitesse et la distance */
|
|||
}
|
|||
}
|
|||
| 317 | racelik | }
|
|
| 494 | racelik | //Ajouter vos tests unitaires dans cette fonction.
|
|
| 161 | racelik | void tests_unitaires(void)
|
|
| 494 | racelik | {
|
|
| 161 | racelik | if (5!=5)
|
|
| 494 | racelik | {
|
|
printf ("Erreur Test unitaire basique.\n");
|
|||
exit(-1);
|
|||
| 137 | flgamel | }
|
|
| 494 | racelik | ||
/** Tests pour voir "GPGGA" en debut de chaine pour la fonction trame_cmp **/
|
|||
| 161 | racelik | if (trame_cmp("$GPRMC suite chaine","GPGGA") != 0)
|
|
| 494 | racelik | {
|
|
printf ("Erreur Test unitaire trame_cmp_1\n");
|
|||
exit(-1);
|
|||
| 161 | racelik | }
|
|
if (trame_cmp("$GPGGA suite chaine","GPGGA") != 1)
|
|||
| 494 | racelik | {
|
|
printf ("Erreur Test unitaire trame_cmp_2\n");
|
|||
exit(-1);
|
|||
| 161 | racelik | }
|
|
| 280 | racelik | if (trame_cmp("$GPRMC suite chaine","GPRMC") != 1)
|
|
| 494 | racelik | {
|
|
printf ("Erreur Test unitaire trame_cmp_3\n");
|
|||
exit(-1);
|
|||
| 161 | racelik | }
|
|
| 280 | racelik | if (trame_cmp("$APRMC suite chaine","GPGGA") != 0)
|
|
| 494 | racelik | {
|
|
printf ("Erreur Test unitaire trame_cmp_4\n");
|
|||
exit(-1);
|
|||
| 280 | racelik | }
|
|
| 161 | racelik | ||
| 494 | racelik | /** test unitaire pour la fonction decode_int **/
|
|
| 161 | racelik | if (decode_int('0') != 0)
|
|
{
|
|||
| 280 | racelik | printf ("erreur test unitaire decode_int_10\n");
|
|
exit(-1);
|
|||
| 137 | flgamel | }
|
|
| 161 | racelik | ||
if (decode_int('1') != 1)
|
|||
{
|
|||
| 280 | racelik | printf ("erreur test unitaire decode_int_11\n");
|
|
exit(-1);
|
|||
| 137 | flgamel | }
|
|
| 161 | racelik | ||
if (decode_int('2') != 2)
|
|||
{
|
|||
| 280 | racelik | printf ("erreur test unitaire decode_int_12\n");
|
|
exit(-1);
|
|||
| 137 | flgamel | }
|
|
| 161 | racelik | ||
if (decode_int('3') != 3)
|
|||
{
|
|||
| 280 | racelik | printf ("erreur test unitaire decode_int_13\n");
|
|
exit(-1);
|
|||
| 137 | flgamel | }
|
|
| 161 | racelik | ||
if (decode_int('4') != 4)
|
|||
{
|
|||
| 280 | racelik | printf ("erreur test unitaire decode_int_14\n");
|
|
exit(-1);
|
|||
| 137 | flgamel | }
|
|
| 161 | racelik | ||
if (decode_int('5') != 5)
|
|||
{
|
|||
| 280 | racelik | printf ("erreur test unitaire decode_int_15\n");
|
|
exit(-1);
|
|||
| 137 | flgamel | }
|
|
| 161 | racelik | ||
if (decode_int('6') != 6)
|
|||
{
|
|||
| 280 | racelik | printf ("erreur test unitaire decode_int_16\n");
|
|
exit(-1);
|
|||
| 161 | racelik | }
|
|
if (decode_int('7') != 7)
|
|||
{
|
|||
| 280 | racelik | printf ("erreur test unitaire decode_int_17\n");
|
|
exit(-1);
|
|||
| 161 | racelik | }
|
|
if (decode_int('8') != 8)
|
|||
{
|
|||
| 280 | racelik | printf ("erreur test unitaire decode_int_18\n");
|
|
exit(-1);
|
|||
| 161 | racelik | }
|
|
if (decode_int('9') != 9)
|
|||
{
|
|||
| 280 | racelik | printf ("erreur test unitaire decode_int_19\n");
|
|
exit(-1);
|
|||
| 161 | racelik | }
|
|
if (decode_int('A') != -1)
|
|||
{
|
|||
| 494 | racelik | printf ("erreur test unitaire decode_int_20\n");
|
|
| 280 | racelik | exit(-1);
|
|
| 161 | racelik | }
|
|
| 280 | racelik | if (decode_int('Z') != -1)
|
|
{
|
|||
| 494 | racelik | printf ("erreur test unitaire decode_int_30\n");
|
|
| 280 | racelik | exit(-1);
|
|
}
|
|||
| 494 | racelik | /** test unitaire pour la fonction decode_nombre **/
|
|
if (decode_nombre("1234", 0, 2) != 12)
|
|||
| 280 | racelik | {
|
|
printf ("Erreur test unitaire decode_nombre_1\n");
|
|||
exit(-1);
|
|||
}
|
|||
| 494 | racelik | if (decode_nombre("6789", 1, 3) != 789)
|
|
| 280 | racelik | {
|
|
printf ("Erreur test unitaire decode_nombre_2\n");
|
|||
exit(-1);
|
|||
}
|
|||
| 494 | racelik | /** test unitaire pour la fonction decode_latitude, avec une precision de 10e-3 **/
|
|
if (((decode_latitude("$GPGGA,141914.00,4545.0000,N,00306.6036,E,1,05,3.4,499.3,M,,M,,*7D"))-45.75) > 10-3) /* test si la latitude est coherente */
|
|||
| 280 | racelik | {
|
|
printf ("Erreur test unitaire decode_latitude\n");
|
|||
exit(-1);
|
|||
}
|
|||
| 494 | racelik | /** test unitaire pour la fonction decode_longitude, avec une precision de 10e-3 **/
|
|
if (((decode_longitude("$GPGGA,141914.00,4545.0000,N,00306.6036,E,1,05,3.4,499.3,M,,M,,*7D"))-3.11006) > 10-3) /* test si la longitude est coherente */
|
|||
| 280 | racelik | {
|
|
printf ("Erreur test unitaire decode_longitude\n");
|
|||
exit(-1);
|
|||
}
|
|||
| 317 | racelik | ||
| 494 | racelik | /** test unitaire pour la fonction decode_trame qui renvoie la latitude et la longitude, avec precision de 10e-3 **/
|
|
position p1;
|
|||
decode_trame("$GPGGA,141916.00,4545.0484,N,00306.6037,E,1,05,3.4,500.0,M,,M,,*7A", &p1); /* exemple avec cette trame GPGGA pour la latitude et longitude */
|
|||
| 317 | racelik | ||
| 494 | racelik | if ((p1.latitude)-45.750806 > 10-3)
|
|
| 317 | racelik | {
|
|
| 494 | racelik | printf ("Erreur test unitaire decode_trame coordonnees (latitude)\n");
|
|
| 317 | racelik | exit(-1);
|
|
}
|
|||
| 494 | racelik | if ((p1.longitude)-3.110061 > 10-3)
|
|
| 317 | racelik | {
|
|
| 494 | racelik | printf ("Erreur test unitaire decode_trame coordonnees (longitude)\n");
|
|
| 317 | racelik | exit(-1);
|
|
}
|
|||
| 494 | racelik | /** declaration pts coordonnees de 2 villes pour les tests unitaires de calcul_distance (Paris et Lyon) **/
|
|
/*Paris*/
|
|||
position p2;
|
|||
p2.latitude = 48.8588897;
|
|||
p2.longitude = 2.320041;
|
|||
| 317 | racelik | ||
| 494 | racelik | /*Lyon*/
|
|
position p3;
|
|||
p3.latitude = 45.7578137;
|
|||
p3.longitude = 4.8320114;
|
|||
/* test calcul_distance entre Paris et Lyon pour voir si la distance est coherente */
|
|||
if((calcul_distance(&p2,&p3)>392) && (calcul_distance(&p2,&p3)<394))
|
|||
{}
|
|||
else
|
|||
| 317 | racelik | {
|
|
| 494 | racelik | printf ("Erreur test unitaire calcul_distance\n");
|
|
| 317 | racelik | exit(-1);
|
|
}
|
|||
| 494 | racelik | /** declaration 2 points coordonnes aleatoirement (latitude et longitude) pour simuler une diff?rence de position de 1s **/
|
|
position p4; /* premier point coordonnee */
|
|||
p4.latitude = 45.75040333;
|
|||
p4.longitude = 3.110065;
|
|||
position p5; /* second point coordonnee */
|
|||
p5.latitude = 45.7500;
|
|||
p5.longitude = 3.11006;
|
|||
if((calcul_vitesse(p4,p5)-1613.443) > 10e-3) /* precision de 10e-3 */
|
|||
| 317 | racelik | {
|
|
| 494 | racelik | printf ("Erreur test unitaire calcul_vitesse\n");
|
|
| 317 | racelik | exit(-1);
|
|
}
|
|||
| 494 | racelik | /*test fonction distance_a_la_proche_zone avec l'indice de la zone la plus proche dans le test, la zone 3 correspond a la 4eme zone du tableau zones[]*/
|
|
float distance = 0.0;
|
|||
int nb_zone = 3; /* variable concernant le test unitaire distance_a_la_plus_proche_zone */
|
|||
if(distance_a_la_plus_proche_zone(p4,zones,nb_zone,&distance) != 2)
|
|||
| 317 | racelik | {
|
|
| 494 | racelik | printf ("Erreur test unitaire calcul distance_a_la_plus_proche_zone\n");
|
|
| 317 | racelik | exit(-1);
|
|
}
|
|||
| 494 | racelik | }
|
|
| 317 | racelik | ||
| 494 | racelik | /* Ne pas modifier cette fonction*/
|
|
| 161 | racelik | int main(int argc,char ** argv)
|
|
{
|
|||
tests_unitaires();
|
|||
| 494 | racelik | /* Affichage des trames definies dans la table trames.*/
|
|
| 161 | racelik | printf ("Trames de tests tableau trames:\n");
|
|
int i=0;
|
|||
while (trames[i])
|
|||
traitement(trames[i++]);
|
|||
if (!trame_init())
|
|||
exit(-1);
|
|||
| 494 | racelik | /* Affichage des trames du fichier gps.log*/
|
|
| 161 | racelik | char *trame;
|
|
printf ("Trames de tests du fichier gps.log\n");
|
|||
while ((trame = trame_suivante()))
|
|||
traitement(trame);
|
|||
| 494 | racelik | return 0;
|
|
}
|