|
#include <stdio.h>
|
|
#include <stdlib.h>
|
|
#include <strings.h>
|
|
#include "trame.h"
|
|
#include <math.h>
|
|
|
|
/***********************************************
|
|
*
|
|
* D?finition des structures utiles au programme
|
|
*
|
|
************************************************/
|
|
|
|
typedef struct {
|
|
float latitude;
|
|
float longitude;
|
|
} Position ;
|
|
|
|
typedef struct {
|
|
Position rpos;
|
|
float vitmax;
|
|
} Zone ;
|
|
|
|
//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",
|
|
|
|
0};
|
|
|
|
|
|
/**Fonction qui compare une trame et renvoie si elle est conforme ou non **/
|
|
|
|
int trame_cmp(char * trame, char * type)
|
|
{
|
|
int i=0;
|
|
int b=0;
|
|
int c=0;
|
|
int ok;
|
|
|
|
while (trame[b]!='\0')
|
|
{
|
|
b=b+1;
|
|
}
|
|
while (type[c]!='\0')
|
|
{
|
|
c=c+1;
|
|
}
|
|
for(i=0; i<=c-1; i++)
|
|
{
|
|
|
|
{
|
|
if (trame[i+1]==type[i])
|
|
{
|
|
ok=1; /** la trame est valide **/
|
|
}
|
|
else
|
|
{
|
|
return 0;
|
|
}
|
|
}
|
|
}
|
|
return ok;
|
|
}
|
|
|
|
/** Fonction qui prends en argument une cha?ne de caract?res (correspondant ? un chiffre) et qui le caste en entier (int) **/
|
|
int decode_int (char c)
|
|
{
|
|
int i=c-48;
|
|
int ok=-1;
|
|
if (i>=0 && i<=9){
|
|
ok=i;
|
|
}
|
|
return ok;
|
|
}
|
|
/** M?me principe que decode_int mais avec une cha?ne de caract?re repr?sentant un nombre **/
|
|
int decode_nombre(char *ch,int n)
|
|
{
|
|
int a=0;
|
|
int nb=0; /**contient le nombre qui ?tait initialement une ch?ine de caract?res **/
|
|
int i; /** compteur pour la boucle for **/
|
|
for (i=0; i<=n-1; i++){
|
|
|
|
a=decode_int(ch[i]);
|
|
nb=nb*10+a; /** on multiplie par 10 pour d?caler d'une dizaine **/
|
|
}
|
|
|
|
return nb;
|
|
}
|
|
/**************************************************************************
|
|
*
|
|
* Fonctions qui convertisseurs la latitude puis la longitude en flottants
|
|
*
|
|
***************************************************************************/
|
|
|
|
float convertisseur_latitude (char *ch)
|
|
{
|
|
|
|
float latitude=0.0;
|
|
float degre=(float)decode_nombre(ch,2);
|
|
float decim1=(float)decode_nombre(ch,4)-degre*pow(10,2);
|
|
|
|
char chaine[50]; /** cha?ne qui va nous permettre de stocker les derniers caract?res correspondant au derniers chiffres de la latitude **/
|
|
int i;
|
|
for (i=5; i<9; i++)
|
|
{
|
|
chaine[i-5]=ch[i];
|
|
}
|
|
float decim2=(float)decode_nombre(chaine, 4); /** on r?cup?re la fin de la cha?ne pour ensuite la convertir **/
|
|
decim1=decim1/60.0;
|
|
decim2=decim2/600000.0;
|
|
latitude=degre+decim1+decim2;
|
|
|
|
return latitude;
|
|
}
|
|
|
|
float convertisseur_longitude (char *ch)
|
|
{
|
|
float longitude=0.0;
|
|
float degre=(float)decode_nombre(ch,3); /** on r?cup?re dans un premier temps les degr?s car ils ne subissent pas d'op?ration **/
|
|
float decim1=(float)decode_nombre(ch,5)-degre*pow(10,2); /** on r?cup?re ici les 3 premiers chiffres de la partie d?cimale **/
|
|
|
|
char chaine[50]; /** cha?ne qui va nous permettre de stocker les derniers caract?res correspondant au derniers chiffres de la longitude **/
|
|
int i;
|
|
for (i=6; i<10; i++)
|
|
{
|
|
chaine[i-6]=ch[i]; /** on r?cup?re la fin de la cha?ne pour ensuite la convertir **/
|
|
}
|
|
float decim2=(float)decode_nombre(chaine, 4);
|
|
decim1=decim1/60.0;
|
|
decim2=decim2/600000.0;
|
|
longitude=degre+decim1+decim2;
|
|
|
|
return longitude;
|
|
}
|
|
|
|
/** Fonction qui calcule indiff?remment la latitude ou la longitude **/
|
|
float latoulong (char *ch)
|
|
{
|
|
int i;
|
|
i=strlen(ch);
|
|
if (i==9)
|
|
{
|
|
return convertisseur_latitude(ch);
|
|
}
|
|
else
|
|
{
|
|
return convertisseur_longitude(ch);
|
|
}
|
|
}
|
|
|
|
/********************************************************************************************
|
|
*
|
|
* Fonctions qui calcule la distance en km entre deux trames et calcule de la vitesse en km/h
|
|
*
|
|
*********************************************************************************************/
|
|
|
|
|
|
float calcul_distance(Position p_1, Position p_2)
|
|
{
|
|
float S;
|
|
float R=6371.00; /**rayon de la Terre**/
|
|
float P=2*M_PI*R;
|
|
float c=P/360.00;
|
|
|
|
S=c*pow((pow((p_2.latitude-p_1.latitude),2)+pow((p_2.longitude-p_1.longitude),2)),0.5);
|
|
return S;
|
|
}
|
|
|
|
float calcul_vitesse(Position p_1, Position p_2)
|
|
{
|
|
float vitesse;
|
|
vitesse=calcul_distance(p_1, p_2)*3600; /** on multiplie par 3600 car on donne un temps entre deux trames ?quivalent ? 1s **/
|
|
return vitesse;
|
|
}
|
|
|
|
/********************************************************************************************
|
|
*
|
|
* Fonction qui d?code la trame et r?cup?rer la latitude et la longitude
|
|
*
|
|
*********************************************************************************************/
|
|
|
|
|
|
int decode_trame(char* t, Position* p){
|
|
int ok,i;
|
|
char latitudecopie[9]; /** chaine qui va contenir la latitude extraite de la trame **/
|
|
char longitudecopie[10]; /** chaine qui va contenir la longitude extraite de la trame **/
|
|
|
|
ok=trame_cmp(t,"GPGGA"); /** comparaison pour voir si la trame est valide ? l'aide de la fonction faite au TP1 **/
|
|
if(ok==1)
|
|
{
|
|
for(i=0;i<9;++i)
|
|
{
|
|
latitudecopie[i]=t[17+i]; /** on a un tableau contenant la latitude en cha?ne de caract?res **/
|
|
}
|
|
(*p).latitude=convertisseur_latitude(latitudecopie);
|
|
|
|
for(i=0;i<10;++i)
|
|
{
|
|
longitudecopie[i]=t[29+i]; /** on a un tableau contenant la longitude en cha?ne de caract?res **/
|
|
|
|
}
|
|
(*p).longitude=convertisseur_longitude(longitudecopie); /** on convertie la longitude en flottant pour l'utiliser ensuite si besoin est **/
|
|
|
|
}else{
|
|
exit(-1);
|
|
} }
|
|
|
|
/***********************************************************************************************
|
|
*
|
|
* Fonction qui rep?re la zone dangereuse la plus proche de l? o? se trouve l'utilisateur du GPS
|
|
*
|
|
************************************************************************************************/
|
|
|
|
int distance_a_la_plus_proche_zone(Position p, Zone r[],int nb_zones, float *d)
|
|
{ int i;
|
|
float distance;
|
|
distance=calcul_distance(p, r[0].rpos); /** Initialise la premiere distance ? comparer **/
|
|
int index; /** Variable que l'on retourne ? la fin de la fonction **/
|
|
|
|
for (i=0; i<nb_zones; i++)
|
|
|
|
{
|
|
*d=calcul_distance(p, r[i].rpos);
|
|
if (*d<=distance) /**si la distance entre la position p et celle du radar ? la i-?me position recens?e est inf?rieur ? la premi?re distance calcul?e, distance prends cette valeur. on retrouvera ainsi de suite la zone dangereuse la plus pr?s de la position de la voiture! **/
|
|
{
|
|
distance=*d;
|
|
index=i;
|
|
}
|
|
|
|
}
|
|
return index;
|
|
}
|
|
//Fonction ? modifier !!!!!
|
|
void traitement(char * trame)
|
|
{
|
|
static int cpt=0;
|
|
cpt++;
|
|
|
|
if((trame[5]==65)&&(trame[4]==71)&&(trame[3]==71))
|
|
{
|
|
printf ("> %s\n",trame);
|
|
}
|
|
|
|
}
|
|
|
|
//Ajouter vos tests unitaires dans cette fonction.
|
|
void tests_unitaires(void){
|
|
|
|
|
|
Position a={0.0, 1.0};
|
|
Position *p;
|
|
p=&a;
|
|
|
|
/******************************************
|
|
*
|
|
* Tests unitaires de diff?rentes fonctions
|
|
*
|
|
*******************************************/
|
|
|
|
if (5!=5){
|
|
printf ("Erreur Test unitaire basique.\n");
|
|
exit(-1);
|
|
}
|
|
if (trame_cmp("$GPGGA,141922.00,4545.1810,N,00306.6046,E,1,05,3.4,500.6,M,,M,,*74","GPGGA")!=1){
|
|
printf ("Erreur 1 Test unitaire trame_cmp.\n");
|
|
exit(-1);
|
|
}
|
|
if (trame_cmp("$GPRMC suite chaine","GPGGA")!=0){
|
|
printf ("Erreur 2 Test unitaire trame_cmp.\n");
|
|
exit(-1);
|
|
}
|
|
if (trame_cmp("$GPRMC... ", "GPRMC" )!=1){
|
|
printf ("Erreur 3 Test unitaire trame_cmp.\n");
|
|
exit(-1);
|
|
}
|
|
if (trame_cmp("$APRMC...", "GPGGA")!=0){
|
|
printf ("Erreur 4 Test unitaire trame_cmp.\n");
|
|
exit(-1);
|
|
}
|
|
if(decode_int('G')!=-1){
|
|
printf("Erreur 1 test unitaire decode_int.\n");
|
|
}
|
|
if(decode_int('6')!=6){
|
|
printf("Erreur 2 test unitaire decode_int.\n");
|
|
}
|
|
if(decode_int('B')!=(-1)){
|
|
printf("Erreur 3 test unitaire decode_int.\n");
|
|
}
|
|
if(decode_int('4')!=4){
|
|
printf("Erreur 4 test unitaire decode_int.\n");
|
|
}
|
|
if(decode_nombre("789654",3)!=789){
|
|
printf("Erreur sur decode_nombre.\n");
|
|
}
|
|
if(decode_nombre("789654",4)!=7896){
|
|
printf("Erreur sur decode_nombre.\n");
|
|
}
|
|
if(convertisseur_latitude("4545.1810")<45.7000 || convertisseur_latitude("4545.1810")>45.8000 ) /**comme c'est des float on ne peut que encadrer la valeur **/
|
|
{
|
|
printf("Erreur sur le convertisseur latitude");
|
|
}
|
|
if(convertisseur_longitude("00306.6046")<03.0000 || convertisseur_latitude("00306.6046")>03.4000 )
|
|
{
|
|
printf("Erreur sur le convertisseur longitude");
|
|
}
|
|
if (decode_trame("$GPGGA,141922.00,4545.1810,N,00306.6046,E,1,05,3.4,500.6,M,,M,,*74",p)==1) {
|
|
printf("Erreur sur decode_trame");
|
|
}
|
|
|
|
/**********************************************************************************************************************************
|
|
*
|
|
* Tests de la fonction decode_trame et r?cup?ration de la latitude et la longitude des deux tests du sujet (TP2 question 5)
|
|
*
|
|
***********************************************************************************************************************************/
|
|
|
|
Position pos;
|
|
pos.latitude=1;
|
|
pos.longitude=1;
|
|
Position *posi;
|
|
posi=&pos;
|
|
|
|
Position position9, position8; /** D?claration des positions qui contiendront les coordon?es GPS des trames du sujet **/
|
|
|
|
|
|
decode_trame("$GPGGA,141914.00,4545.0000,N,00306.6036,E,1,05,3.4,499.3,M,,M,,*7D", &position8);
|
|
decode_trame("$GPGGA,141915.00,4545.0242,N,00306.6039,E,1,05,3.4,499.5,M,,M,,*72", &position9);
|
|
//printf("position latitude trame 1: %f \n", position8.latitude); //on r?cup?re les donn?es des positions de la latitude de la trame 2 question 5
|
|
//printf("position longtitude trame 1: %f \n", position8.longitude);//on r?cup?re les donn?es des positions de la longitude de la trame 2 question 5
|
|
//printf("position latitude trame 2: %f \n", position9.latitude); //on r?cup?re les donn?es des positions de la latitude de la trame 2 question 5
|
|
//printf("position longtitude trame 2: %f \n", position9.longitude); //on r?cup?re les donn?es des positions de la longitude de la trame 2 question 5
|
|
|
|
/**********************************************************************************************************************************
|
|
*
|
|
* Initialisation des coordonn?es GPS de diff?rentes villes de france pour tester les fonctions calcule_distance et calcule_vitesse
|
|
*
|
|
***********************************************************************************************************************************/
|
|
|
|
Position position1, position2, position3, position4, position5,position6, position7; //definitions des variables de type Position qui vont r?pertorier les diff?rentes coordonn?es des diff?rentes villes
|
|
|
|
/** Paris **/
|
|
position2.latitude=48.85666666666667;
|
|
position2.longitude=2.3519444444444444;
|
|
|
|
/** Rouen **/
|
|
position1.latitude=49.443888888888885;
|
|
position1.longitude=1.1033333333333335;
|
|
|
|
/** Lille **/
|
|
position3.latitude=50.63722222222222;
|
|
position3.longitude=3.063333333333333;
|
|
|
|
/** Bordeaux **/
|
|
position4.latitude=44.8404400;
|
|
position4.longitude=-0.5805000;
|
|
|
|
/** Orl?ans **/
|
|
position5.latitude=47.9028900;
|
|
position5.longitude=1.9038900;
|
|
|
|
/** Polytech Clermont **/
|
|
position6.latitude=45.77140808105469;
|
|
position6.longitude=3.0903782844543457;
|
|
|
|
/** Isima Clermont **/
|
|
position7.latitude=45.75858688354492;
|
|
position7.longitude=3.1108570098876953;
|
|
|
|
|
|
//printf("**** Test des distances entre plusieurs villes : ***** \n\n");
|
|
//printf("Distance Paris-Rouen=%f km\n", calcul_distance(position2, position1)); /**Paris-Rouen = 135 km**/
|
|
//printf("Distance Paris-Lille=%f km\n", calcul_distance(position2, position3)); /**Paris-Lille = 225 km**/
|
|
//printf("Distance Paris-Bordeaux=%f km\n", calcul_distance(position2, position4)); /**Paris-Bordeaux = 584 km**/
|
|
//printf("Distance Paris-Orleans=%f km\n\n", calcul_distance(position2, position5)); /**Paris-Orl?ans = 132 km**/
|
|
//printf("Distance entre les trames du sujet=%f km soit approximativement %.1f. m\n\n", calcul_distance(position8, position9), calcul_distance(position8, position9)*1000); /**distance entre les deux trames du sujet question 5**/
|
|
|
|
|
|
Zone zones[] = {{{44.7887762, -3.012}, 50},{{44.7891220, -3,013}, 70},{{45.7604, 3.120065}, 130}}; /**tableaux des zones dangereuses donn? dans l'?nonc? + une autre pour les tests**/
|
|
float dist;
|
|
int index;
|
|
|
|
/** Test de la fonction distance_a_la_plus_proche_zone **/
|
|
index=distance_a_la_plus_proche_zone(position9, zones, 3, &dist);
|
|
printf("**** Tests sur les zones dangereuses ***** \n\n");
|
|
|
|
printf("Index de la zone dangereuse la plus proche : %i.\n\n", index);
|
|
float vitesse;
|
|
int alarme;
|
|
|
|
vitesse=calcul_vitesse(position8, position9);
|
|
|
|
/** Affichage de la vitesse en km/h **/
|
|
|
|
printf("La vitesse du conducteur ou de la conductrice est de : %.0f km/h et la vitesse maximale autorisee est de : %.0f km/h.\n", vitesse, zones[index].vitmax);
|
|
if (vitesse>zones[index].vitmax) /** si la vitesse est sup?rieure au seuil autoris? alors l'alarme va s'activer car le conducteur ou la conductrice roule trop vite **/
|
|
{
|
|
printf("Alarme : active. Vous ?tes au dessus de la limitation. \n\n");
|
|
alarme=1;
|
|
}
|
|
else
|
|
{
|
|
printf("Alarme : non active. \n\n"); /** le conducteur est en dessous de la vitesse maximale autoris?e : pas besoin de d?clancher l'alarme **/
|
|
alarme=0;
|
|
}
|
|
printf("La vitesse entre l'utilisateur du GPS et la prochaine zone dangereuse est de : %.2f km soit %.0f m.", dist, dist*1000);
|
|
}
|
|
|
|
|
|
|
|
// Ne pas modifier cette fonction
|
|
int main(int argc,char ** argv)
|
|
{
|
|
|
|
tests_unitaires();
|
|
|
|
// 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;
|
|
}
|