Commit 47cf534b authored by Reza  ANSARI's avatar Reza ANSARI
Browse files

Correction ds generation track-files (angles converti en degres et...

Correction ds generation track-files (angles converti en degres et DeltaLongitude), Reza 25/11/2018-
parent 3d98a78b
......@@ -72,10 +72,11 @@ Les satellites dont on telecharge les TLE sont listes dans les variables: naviga
.pour changer l'emplacement, ajouter a l'ordre "make ..." les arguments:
OBJ="la ou je veux mettre mes objets" EXE="la ou je veux mettre mes executables"
.Reza 25/11/2018 : Avec le build comme indiqu ci-dessus, j'ai du faire
.Reza 25/11/2018 : Avec le build comme indique ci-dessus, j'ai du faire
> make SGP4INC=${HERESGP4}/sgp4/libsgp4/ SGP4LIB=${HERESGP4}/libsgp4/
.Reza 25/11/2018 : Au CC-IN2P3, il faut ajouter l'option de compilation -std=c++11
J'ai donc modifie le Makefile en ajoutant un flag optionnel de compilation - pour faire make au CC
J'ai donc modifie le Makefile en ajoutant un flag optionnel de compilation - pour faire make au CC
> set HERESGP4 = /pbs/throng/baoradio/Library/SGP4_Code/
> make SGP4INC=${HERESGP4}/sgp4/libsgp4/ SGP4LIB=${HERESGP4}/libsgp4/ SGP4CCFLG="-std=c++11"
-----
......@@ -91,7 +92,7 @@ J'ai donc modifie le Makefile en ajoutant un flag optionnel de compilation - pou
> ./Objs/predictsatsgp4 -D 37847 satname YYYMMDD/galileo.txt [...]
.Reza 25/11/2018 : j'ai ajoute l'option -K au programme, pour creer un fichier avec la trace (track) du satellite
..../AnaPAON4/Satellites/Objs/predictsatsgp4 -T "2018/10/20 15:16:00" -H 0.,81. -T 41859 -p 2 TLE_20181123/*.txt
..../AnaPAON4/Satellites/Objs/predictsatsgp4 -T "2018/10/20 15:16:00" -H 0.,81. -K 41859 -p 2 galileo.txt
-----
5-/ Remarques
......
......@@ -379,11 +379,13 @@ void SaveTrack(string & satname, SGP4 & sgp4, Observer & obs, DateTime& datestar
string flnm = "trk_"+satname+".txt";
ofstream ofs(flnm);
ofs << "#### Track for satellite : " << satname << endl;
ofs << "## Date/Time Azimuth Elevation Latitude Longitude " << endl;
ofs << "## Date/Time UTC Azimuth Elevation Latitude Longitude DeltaLongitude (Angles in degree)" << endl;
if (prtlev>1) {
cout << "#### Track for satellite : " << satname << endl;
cout << "## Date/Time Azimuth Elevation Latitude Longitude " << endl;
cout << "## Date/Time UTC Azimuth Elevation Latitude Longitude DeltaLongitude (Angles in degree)" << endl;
}
CoordGeodetic obsgeo=obs.GetLocation();
DateTime datecur = datestart;
while(datecur <= dateend+tspaninc) {
Eci eci = sgp4.FindPosition(datecur);
......@@ -395,9 +397,15 @@ void SaveTrack(string & satname, SGP4 & sgp4, Observer & obs, DateTime& datestar
//...Convert satellite position to geodetic coordinates
CoordGeodetic geo = eci.ToGeodetic();
// Vector Vungeo = AzAlt2Vec(RadiansToDegrees(geo.longitude),RadiansToDegrees(geo.latitude));
ofs <<datecur<<" "<<topo.azimuth<<" "<<topo.elevation<<" "<<geo.latitude<<" "<<geo.longitude<<endl;
if (prtlev>1)
cout <<datecur<<" "<<topo.azimuth<<" "<<topo.elevation<<" "<<geo.latitude<<" "<<geo.longitude<<endl;
ofs <<datecur<<" "<<RadiansToDegrees(topo.azimuth)<<" "<<RadiansToDegrees(topo.elevation)
<<" "<<RadiansToDegrees(geo.latitude)<<" "<<RadiansToDegrees(geo.longitude)
<<" "<<RadiansToDegrees(geo.longitude-obsgeo.longitude)<<endl;
if (prtlev>1) {
cout <<datecur<<" "<<RadiansToDegrees(topo.azimuth)<<" "<<RadiansToDegrees(topo.elevation)
<<" "<<RadiansToDegrees(geo.latitude)<<" "<<RadiansToDegrees(geo.longitude)
<<" "<<RadiansToDegrees(geo.longitude-obsgeo.longitude)<<endl;
}
datecur = datecur + tspaninc;
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment