[TLE] New tle api to load TLE files only once
This is a wrapper over the existing code, underlying parsing is not changed.pull/36/head
parent
289eaf14dd
commit
7177ad68bd
|
@ -19,8 +19,8 @@ all:
|
||||||
rffit: rffit.o sgdp4.o satutl.o deep.o ferror.o dsmin.o simplex.o versafit.o
|
rffit: rffit.o sgdp4.o satutl.o deep.o ferror.o dsmin.o simplex.o versafit.o
|
||||||
gfortran -o rffit rffit.o sgdp4.o satutl.o deep.o ferror.o dsmin.o simplex.o versafit.o $(LFLAGS)
|
gfortran -o rffit rffit.o sgdp4.o satutl.o deep.o ferror.o dsmin.o simplex.o versafit.o $(LFLAGS)
|
||||||
|
|
||||||
rfpng: rfpng.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o
|
rfpng: rfpng.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o rftles.o
|
||||||
gfortran -o rfpng rfpng.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o $(LFLAGS)
|
gfortran -o rfpng rfpng.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o rftles.o $(LFLAGS)
|
||||||
|
|
||||||
rfedit: rfedit.o rfio.o rftime.o
|
rfedit: rfedit.o rfio.o rftime.o
|
||||||
$(CC) -o rfedit rfedit.o rfio.o rftime.o -lm
|
$(CC) -o rfedit rfedit.o rfio.o rftime.o -lm
|
||||||
|
@ -31,8 +31,8 @@ rffind: rffind.o rfio.o rftime.o
|
||||||
rftrack: rftrack.o rfio.o rftime.o rftrace.o sgdp4.o satutl.o deep.o ferror.o
|
rftrack: rftrack.o rfio.o rftime.o rftrace.o sgdp4.o satutl.o deep.o ferror.o
|
||||||
$(CC) -o rftrack rftrack.o rfio.o rftime.o rftrace.o sgdp4.o satutl.o deep.o ferror.o -lm
|
$(CC) -o rftrack rftrack.o rfio.o rftime.o rftrace.o sgdp4.o satutl.o deep.o ferror.o -lm
|
||||||
|
|
||||||
rfplot: rfplot.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o
|
rfplot: rfplot.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o rftles.o
|
||||||
gfortran -o rfplot rfplot.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o $(LFLAGS)
|
gfortran -o rfplot rfplot.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o rftles.o $(LFLAGS)
|
||||||
|
|
||||||
rffft: rffft.o rftime.o
|
rffft: rffft.o rftime.o
|
||||||
$(CC) -o rffft rffft.o rftime.o -lfftw3f -lm
|
$(CC) -o rffft rffft.o rftime.o -lfftw3f -lm
|
||||||
|
|
|
@ -33,8 +33,8 @@ all:
|
||||||
rffit: rffit.o sgdp4.o satutl.o deep.o ferror.o dsmin.o simplex.o versafit.o
|
rffit: rffit.o sgdp4.o satutl.o deep.o ferror.o dsmin.o simplex.o versafit.o
|
||||||
$(CC) -o rffit rffit.o sgdp4.o satutl.o deep.o ferror.o dsmin.o simplex.o versafit.o $(LFLAGS)
|
$(CC) -o rffit rffit.o sgdp4.o satutl.o deep.o ferror.o dsmin.o simplex.o versafit.o $(LFLAGS)
|
||||||
|
|
||||||
rfpng: rfpng.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o
|
rfpng: rfpng.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o rftles.o
|
||||||
$(CC) -o rfpng rfpng.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o $(LFLAGS)
|
$(CC) -o rfpng rfpng.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o rftles.o $(LFLAGS)
|
||||||
|
|
||||||
rfedit: rfedit.o rfio.o rftime.o
|
rfedit: rfedit.o rfio.o rftime.o
|
||||||
$(CC) -o rfedit rfedit.o rfio.o rftime.o -lm
|
$(CC) -o rfedit rfedit.o rfio.o rftime.o -lm
|
||||||
|
@ -45,8 +45,8 @@ rffind: rffind.o rfio.o rftime.o
|
||||||
rftrack: rftrack.o rfio.o rftime.o rftrace.o sgdp4.o satutl.o deep.o ferror.o
|
rftrack: rftrack.o rfio.o rftime.o rftrace.o sgdp4.o satutl.o deep.o ferror.o
|
||||||
$(CC) -o rftrack rftrack.o rfio.o rftime.o rftrace.o sgdp4.o satutl.o deep.o ferror.o -lm
|
$(CC) -o rftrack rftrack.o rfio.o rftime.o rftrace.o sgdp4.o satutl.o deep.o ferror.o -lm
|
||||||
|
|
||||||
rfplot: rfplot.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o versafit.o dsmin.o simplex.o
|
rfplot: rfplot.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o versafit.o dsmin.o simplex.o rftles.o
|
||||||
$(CC) -o rfplot rfplot.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o versafit.o dsmin.o simplex.o $(LFLAGS)
|
$(CC) -o rfplot rfplot.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o versafit.o dsmin.o simplex.o rftles.o $(LFLAGS)
|
||||||
|
|
||||||
rffft: rffft.o rftime.o
|
rffft: rffft.o rftime.o
|
||||||
$(CC) -o rffft rffft.o rftime.o -lfftw3f -lm $(LFLAGS)
|
$(CC) -o rffft rffft.o rftime.o -lfftw3f -lm $(LFLAGS)
|
||||||
|
|
12
makefile
12
makefile
|
@ -19,11 +19,11 @@ all:
|
||||||
rffit: rffit.o sgdp4.o satutl.o deep.o ferror.o dsmin.o simplex.o versafit.o
|
rffit: rffit.o sgdp4.o satutl.o deep.o ferror.o dsmin.o simplex.o versafit.o
|
||||||
gfortran -o rffit rffit.o sgdp4.o satutl.o deep.o ferror.o dsmin.o simplex.o versafit.o $(LFLAGS)
|
gfortran -o rffit rffit.o sgdp4.o satutl.o deep.o ferror.o dsmin.o simplex.o versafit.o $(LFLAGS)
|
||||||
|
|
||||||
rfpng: rfpng.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o
|
rfpng: rfpng.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o rftles.o
|
||||||
gfortran -o rfpng rfpng.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o $(LFLAGS)
|
gfortran -o rfpng rfpng.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o rftles.o $(LFLAGS)
|
||||||
|
|
||||||
rfdop: rfdop.o rftrace.o rfio.o rftime.o sgdp4.o satutl.o deep.o ferror.o
|
rfdop: rfdop.o rftrace.o rfio.o rftime.o sgdp4.o satutl.o deep.o ferror.o rftles.o
|
||||||
$(CC) -o rfdop rfdop.o rftrace.o rfio.o rftime.o sgdp4.o satutl.o deep.o ferror.o -lm
|
$(CC) -o rfdop rfdop.o rftrace.o rfio.o rftime.o sgdp4.o satutl.o deep.o ferror.o rftles.o -lm
|
||||||
|
|
||||||
rfedit: rfedit.o rfio.o rftime.o
|
rfedit: rfedit.o rfio.o rftime.o
|
||||||
$(CC) -o rfedit rfedit.o rfio.o rftime.o -lm
|
$(CC) -o rfedit rfedit.o rfio.o rftime.o -lm
|
||||||
|
@ -34,8 +34,8 @@ rffind: rffind.o rfio.o rftime.o
|
||||||
rftrack: rftrack.o rfio.o rftime.o rftrace.o sgdp4.o satutl.o deep.o ferror.o
|
rftrack: rftrack.o rfio.o rftime.o rftrace.o sgdp4.o satutl.o deep.o ferror.o
|
||||||
$(CC) -o rftrack rftrack.o rfio.o rftime.o rftrace.o sgdp4.o satutl.o deep.o ferror.o -lm
|
$(CC) -o rftrack rftrack.o rfio.o rftime.o rftrace.o sgdp4.o satutl.o deep.o ferror.o -lm
|
||||||
|
|
||||||
rfplot: rfplot.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o versafit.o dsmin.o simplex.o
|
rfplot: rfplot.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o versafit.o dsmin.o simplex.o rftles.o
|
||||||
gfortran -o rfplot rfplot.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o versafit.o dsmin.o simplex.o $(LFLAGS)
|
gfortran -o rfplot rfplot.o rftime.o rfio.o rftrace.o sgdp4.o satutl.o deep.o ferror.o versafit.o dsmin.o simplex.o rftles.o $(LFLAGS)
|
||||||
|
|
||||||
rffft: rffft.o rftime.o
|
rffft: rffft.o rftime.o
|
||||||
$(CC) -o rffft rffft.o rftime.o -lfftw3f -lm
|
$(CC) -o rffft rffft.o rftime.o -lfftw3f -lm
|
||||||
|
|
|
@ -0,0 +1,92 @@
|
||||||
|
#include "rftles.h"
|
||||||
|
|
||||||
|
#include "satutl.h"
|
||||||
|
|
||||||
|
#include <stdio.h>
|
||||||
|
#include <stdlib.h>
|
||||||
|
#include <string.h>
|
||||||
|
|
||||||
|
tles_t load_tles(char *tlefile) {
|
||||||
|
tles_t tles;
|
||||||
|
|
||||||
|
tles.orbits = NULL;
|
||||||
|
tles.number_of_elements = 0;
|
||||||
|
|
||||||
|
char filename[1024];
|
||||||
|
|
||||||
|
if (tlefile) {
|
||||||
|
strncpy(filename, tlefile, sizeof(filename));
|
||||||
|
} else {
|
||||||
|
char * env = getenv("ST_TLEDIR");
|
||||||
|
|
||||||
|
if (env == NULL || strlen(env) == 0) {
|
||||||
|
env=".";
|
||||||
|
}
|
||||||
|
|
||||||
|
sprintf(filename, "%s/bulk.tle", env);
|
||||||
|
}
|
||||||
|
|
||||||
|
FILE * file = fopen(filename, "r");
|
||||||
|
|
||||||
|
if (file == NULL) {
|
||||||
|
fprintf(stderr, "TLE file %s not found\n", filename);
|
||||||
|
|
||||||
|
return tles;
|
||||||
|
}
|
||||||
|
|
||||||
|
size_t linesize = 256;
|
||||||
|
char * line = malloc(linesize);
|
||||||
|
ssize_t read;
|
||||||
|
|
||||||
|
// Count number of entries
|
||||||
|
long num_elements = 0;
|
||||||
|
|
||||||
|
while ((read = getline(&line, &linesize, file)) != -1) {
|
||||||
|
if (read > 0 && strncmp(line, "1 ", 2) == 0) {
|
||||||
|
num_elements++;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
tles.orbits = (orbit_t *)calloc(num_elements, sizeof(orbit_t));
|
||||||
|
|
||||||
|
// Rewind and parse file
|
||||||
|
rewind(file);
|
||||||
|
|
||||||
|
while (read_twoline(file, 0, &(tles.orbits[tles.number_of_elements])) == 0) {
|
||||||
|
tles.number_of_elements++;
|
||||||
|
}
|
||||||
|
|
||||||
|
free(line);
|
||||||
|
fclose(file);
|
||||||
|
|
||||||
|
printf("Loaded %ld orbits\n", tles.number_of_elements);
|
||||||
|
|
||||||
|
return tles;
|
||||||
|
}
|
||||||
|
|
||||||
|
void free_tles(tles_t *tles) {
|
||||||
|
if (tles) {
|
||||||
|
free(tles->orbits);
|
||||||
|
tles->number_of_elements = 0;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
orbit_t *get_orbit_by_index(tles_t *tles, long index) {
|
||||||
|
if (tles) {
|
||||||
|
return &(tles->orbits[index]);
|
||||||
|
}
|
||||||
|
|
||||||
|
return NULL;
|
||||||
|
}
|
||||||
|
|
||||||
|
orbit_t *get_orbit_by_catalog_id(tles_t *tles, long satno) {
|
||||||
|
if (tles) {
|
||||||
|
for (long i = 0; i < tles->number_of_elements; i++) {
|
||||||
|
if (tles->orbits[i].satno == satno) {
|
||||||
|
return &(tles->orbits[i]);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
return NULL;
|
||||||
|
}
|
|
@ -0,0 +1,24 @@
|
||||||
|
#ifndef _RFTLES_H
|
||||||
|
#define _RFTLES_H
|
||||||
|
|
||||||
|
#include "sgdp4h.h"
|
||||||
|
|
||||||
|
#ifdef __cplusplus
|
||||||
|
extern "C" {
|
||||||
|
#endif
|
||||||
|
|
||||||
|
typedef struct tles {
|
||||||
|
long number_of_elements;
|
||||||
|
orbit_t *orbits;
|
||||||
|
} tles_t;
|
||||||
|
|
||||||
|
tles_t load_tles(char *tlefile);
|
||||||
|
void free_tles(tles_t *tles);
|
||||||
|
orbit_t *get_orbit_by_index(tles_t *tles, long index);
|
||||||
|
orbit_t *get_orbit_by_catalog_id(tles_t *tles, long satno);
|
||||||
|
|
||||||
|
#ifdef __cplusplus
|
||||||
|
}
|
||||||
|
#endif
|
||||||
|
|
||||||
|
#endif /* _RFTLES_H */
|
217
rftrace.c
217
rftrace.c
|
@ -9,6 +9,9 @@
|
||||||
#include <sys/time.h>
|
#include <sys/time.h>
|
||||||
#include <time.h>
|
#include <time.h>
|
||||||
|
|
||||||
|
|
||||||
|
#include "rftles.h"
|
||||||
|
|
||||||
#define LIM 80
|
#define LIM 80
|
||||||
#define D2R M_PI/180.0
|
#define D2R M_PI/180.0
|
||||||
#define R2D 180.0/M_PI
|
#define R2D 180.0/M_PI
|
||||||
|
@ -41,7 +44,7 @@ double modulo(double x,double y)
|
||||||
int fgetline(FILE *file,char *s,int lim)
|
int fgetline(FILE *file,char *s,int lim)
|
||||||
{
|
{
|
||||||
int c,i=0;
|
int c,i=0;
|
||||||
|
|
||||||
while (--lim > 0 && (c=fgetc(file)) != EOF && c != '\n')
|
while (--lim > 0 && (c=fgetc(file)) != EOF && c != '\n')
|
||||||
s[i++] = c;
|
s[i++] = c;
|
||||||
// if (c == '\n')
|
// if (c == '\n')
|
||||||
|
@ -88,10 +91,10 @@ void obspos_xyz(double mjd,double lng,double lat,float alt,xyz_t *pos,xyz_t *vel
|
||||||
dtheta=dgmst(mjd)*D2R/86400;
|
dtheta=dgmst(mjd)*D2R/86400;
|
||||||
|
|
||||||
pos->x=gc*cos(lat*D2R)*cos(theta*D2R)*XKMPER;
|
pos->x=gc*cos(lat*D2R)*cos(theta*D2R)*XKMPER;
|
||||||
pos->y=gc*cos(lat*D2R)*sin(theta*D2R)*XKMPER;
|
pos->y=gc*cos(lat*D2R)*sin(theta*D2R)*XKMPER;
|
||||||
pos->z=gs*sin(lat*D2R)*XKMPER;
|
pos->z=gs*sin(lat*D2R)*XKMPER;
|
||||||
vel->x=-gc*cos(lat*D2R)*sin(theta*D2R)*XKMPER*dtheta;
|
vel->x=-gc*cos(lat*D2R)*sin(theta*D2R)*XKMPER*dtheta;
|
||||||
vel->y=gc*cos(lat*D2R)*cos(theta*D2R)*XKMPER*dtheta;
|
vel->y=gc*cos(lat*D2R)*cos(theta*D2R)*XKMPER*dtheta;
|
||||||
vel->z=0.0;
|
vel->z=0.0;
|
||||||
|
|
||||||
return;
|
return;
|
||||||
|
@ -103,7 +106,7 @@ void equatorial2horizontal(double mjd,double ra,double de,double lng,double lat,
|
||||||
double h;
|
double h;
|
||||||
|
|
||||||
h=gmst(mjd)+lng-ra;
|
h=gmst(mjd)+lng-ra;
|
||||||
|
|
||||||
*azi=modulo(atan2(sin(h*D2R),cos(h*D2R)*sin(lat*D2R)-tan(de*D2R)*cos(lat*D2R))*R2D,360.0);
|
*azi=modulo(atan2(sin(h*D2R),cos(h*D2R)*sin(lat*D2R)-tan(de*D2R)*cos(lat*D2R))*R2D,360.0);
|
||||||
*alt=asin(sin(lat*D2R)*sin(de*D2R)+cos(lat*D2R)*cos(de*D2R)*cos(h*D2R))*R2D;
|
*alt=asin(sin(lat*D2R)*sin(de*D2R)+cos(lat*D2R)*cos(de*D2R)*cos(h*D2R))*R2D;
|
||||||
|
|
||||||
|
@ -148,7 +151,7 @@ struct site get_site(int site_id)
|
||||||
|
|
||||||
// Change to km
|
// Change to km
|
||||||
alt/=1000.0;
|
alt/=1000.0;
|
||||||
|
|
||||||
// Copy site
|
// Copy site
|
||||||
if (id==site_id) {
|
if (id==site_id) {
|
||||||
s.lat=lat;
|
s.lat=lat;
|
||||||
|
@ -171,7 +174,7 @@ void identify_trace_graves(char *tlefile,struct trace t,int satno,char *freqlist
|
||||||
struct point *p;
|
struct point *p;
|
||||||
struct site s,sg;
|
struct site s,sg;
|
||||||
double *v,*vg;
|
double *v,*vg;
|
||||||
orbit_t orb;
|
orbit_t *orb;
|
||||||
xyz_t satpos,satvel;
|
xyz_t satpos,satvel;
|
||||||
FILE *file;
|
FILE *file;
|
||||||
double dx,dy,dz,dvx,dvy,dvz,r,za;
|
double dx,dy,dz,dvx,dvy,dvz,r,za;
|
||||||
|
@ -204,26 +207,31 @@ void identify_trace_graves(char *tlefile,struct trace t,int satno,char *freqlist
|
||||||
// Mid point
|
// Mid point
|
||||||
imid=t.n/2;
|
imid=t.n/2;
|
||||||
|
|
||||||
// Loop over TLEs
|
// Load TLEs
|
||||||
file=fopen(tlefile,"r");
|
tles_t twolines = load_tles(tlefile);
|
||||||
if (file==NULL) {
|
|
||||||
fprintf(stderr,"TLE file %s not found\n",tlefile);
|
if (twolines.number_of_elements == 0) {
|
||||||
|
fprintf(stderr,"TLE file %s not found or empty\n", tlefile);
|
||||||
return;
|
return;
|
||||||
}
|
}
|
||||||
while (read_twoline(file,satno,&orb)==0) {
|
|
||||||
|
for (long elem = 0; elem < twolines.number_of_elements; elem++) {
|
||||||
|
// Get TLE
|
||||||
|
orb = get_orbit_by_index(&twolines, elem);
|
||||||
|
|
||||||
// Initialize
|
// Initialize
|
||||||
imode=init_sgdp4(&orb);
|
imode=init_sgdp4(orb);
|
||||||
if (imode==SGDP4_ERROR) {
|
if (imode==SGDP4_ERROR) {
|
||||||
printf("Error with %d, skipping\n",orb.satno);
|
printf("Error with %d, skipping\n",orb->satno);
|
||||||
continue;
|
continue;
|
||||||
}
|
}
|
||||||
|
|
||||||
// Loop over points
|
// Loop over points
|
||||||
for (i=0,sum1=0.0,sum2=0.0;i<t.n;i++) {
|
for (i=0,sum1=0.0,sum2=0.0;i<t.n;i++) {
|
||||||
// Get satellite position
|
// Get satellite position
|
||||||
satpos_xyz(t.mjd[i]+2400000.5,&satpos,&satvel);
|
satpos_xyz(t.mjd[i]+2400000.5,&satpos,&satvel);
|
||||||
|
|
||||||
dx=satpos.x-p[i].obspos.x;
|
dx=satpos.x-p[i].obspos.x;
|
||||||
dy=satpos.y-p[i].obspos.y;
|
dy=satpos.y-p[i].obspos.y;
|
||||||
dz=satpos.z-p[i].obspos.z;
|
dz=satpos.z-p[i].obspos.z;
|
||||||
dvx=satvel.x-p[i].obsvel.x;
|
dvx=satvel.x-p[i].obsvel.x;
|
||||||
|
@ -237,7 +245,7 @@ void identify_trace_graves(char *tlefile,struct trace t,int satno,char *freqlist
|
||||||
de=asin(dz/r)*R2D;
|
de=asin(dz/r)*R2D;
|
||||||
equatorial2horizontal(t.mjd[i],ra,de,s.lng,s.lat,&azi,&alt);
|
equatorial2horizontal(t.mjd[i],ra,de,s.lng,s.lat,&azi,&alt);
|
||||||
}
|
}
|
||||||
dx=satpos.x-p[i].grpos.x;
|
dx=satpos.x-p[i].grpos.x;
|
||||||
dy=satpos.y-p[i].grpos.y;
|
dy=satpos.y-p[i].grpos.y;
|
||||||
dz=satpos.z-p[i].grpos.z;
|
dz=satpos.z-p[i].grpos.z;
|
||||||
dvx=satvel.x-p[i].grvel.x;
|
dvx=satvel.x-p[i].grvel.x;
|
||||||
|
@ -254,7 +262,7 @@ void identify_trace_graves(char *tlefile,struct trace t,int satno,char *freqlist
|
||||||
freq0=143050000.0;
|
freq0=143050000.0;
|
||||||
|
|
||||||
// Compute residuals
|
// Compute residuals
|
||||||
for (i=0,rms=0.0;i<t.n;i++)
|
for (i=0,rms=0.0;i<t.n;i++)
|
||||||
rms+=pow(t.freq[i]-(1.0-v[i]/C)*(1.0-vg[i]/C)*freq0,2);
|
rms+=pow(t.freq[i]-(1.0-v[i]/C)*(1.0-vg[i]/C)*freq0,2);
|
||||||
rms=sqrt(rms/(double) t.n);
|
rms=sqrt(rms/(double) t.n);
|
||||||
|
|
||||||
|
@ -262,18 +270,18 @@ void identify_trace_graves(char *tlefile,struct trace t,int satno,char *freqlist
|
||||||
for (i=1,mjd0=0.0;i<t.n;i++)
|
for (i=1,mjd0=0.0;i<t.n;i++)
|
||||||
if (v[i]*v[i-1]<0.0)
|
if (v[i]*v[i-1]<0.0)
|
||||||
mjd0=t.mjd[i];
|
mjd0=t.mjd[i];
|
||||||
|
|
||||||
if (mjd0>0.0)
|
if (mjd0>0.0)
|
||||||
mjd2nfd(mjd0,nfd);
|
mjd2nfd(mjd0,nfd);
|
||||||
else
|
else
|
||||||
strcpy(nfd,"0000-00-00T00:00:00");
|
strcpy(nfd,"0000-00-00T00:00:00");
|
||||||
|
|
||||||
if (rms<1000) {
|
if (rms<1000) {
|
||||||
if (rms<50.0)
|
if (rms<50.0)
|
||||||
printf("%05d: %s %8.1f Hz (%.1f,%.1f)\n",orb.satno,nfd,rms,modulo(azi+180.0,360.0),alt);
|
printf("%05d: %s %8.1f Hz (%.1f,%.1f)\n",orb->satno,nfd,rms,modulo(azi+180.0,360.0),alt);
|
||||||
// printf("%05d: %s %8.3f MHz %8.3f kHz\n",orb.satno,nfd,1e-6*freq0,1e-3*rms);
|
// printf("%05d: %s %8.3f MHz %8.3f kHz\n",orb->satno,nfd,1e-6*freq0,1e-3*rms);
|
||||||
if (flag==0 || rms<rmsmin) {
|
if (flag==0 || rms<rmsmin) {
|
||||||
satnomin=orb.satno;
|
satnomin=orb->satno;
|
||||||
strcpy(nfdmin,nfd);
|
strcpy(nfdmin,nfd);
|
||||||
freqmin=freq0;
|
freqmin=freq0;
|
||||||
rmsmin=rms;
|
rmsmin=rms;
|
||||||
|
@ -283,7 +291,6 @@ void identify_trace_graves(char *tlefile,struct trace t,int satno,char *freqlist
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
fclose(file);
|
|
||||||
fclose(stderr);
|
fclose(stderr);
|
||||||
|
|
||||||
if (flag==1) {
|
if (flag==1) {
|
||||||
|
@ -305,6 +312,7 @@ void identify_trace_graves(char *tlefile,struct trace t,int satno,char *freqlist
|
||||||
}
|
}
|
||||||
|
|
||||||
// Free
|
// Free
|
||||||
|
free_tles(&twolines);
|
||||||
free(p);
|
free(p);
|
||||||
free(v);
|
free(v);
|
||||||
free(vg);
|
free(vg);
|
||||||
|
@ -319,7 +327,7 @@ void identify_trace(char *tlefile,struct trace t,int satno,char *freqlist)
|
||||||
struct point *p;
|
struct point *p;
|
||||||
struct site s;
|
struct site s;
|
||||||
double *v;
|
double *v;
|
||||||
orbit_t orb;
|
orbit_t *orb;
|
||||||
xyz_t satpos,satvel;
|
xyz_t satpos,satvel;
|
||||||
FILE *file;
|
FILE *file;
|
||||||
double dx,dy,dz,dvx,dvy,dvz,r,za;
|
double dx,dy,dz,dvx,dvy,dvz,r,za;
|
||||||
|
@ -329,7 +337,7 @@ void identify_trace(char *tlefile,struct trace t,int satno,char *freqlist)
|
||||||
double rmsmin,freqmin;
|
double rmsmin,freqmin;
|
||||||
struct timeval tv;
|
struct timeval tv;
|
||||||
char tbuf[30];
|
char tbuf[30];
|
||||||
|
|
||||||
// Reloop stderr
|
// Reloop stderr
|
||||||
if (freopen("/tmp/stderr.txt","w",stderr)==NULL)
|
if (freopen("/tmp/stderr.txt","w",stderr)==NULL)
|
||||||
fprintf(stderr,"Failed to redirect stderr\n");
|
fprintf(stderr,"Failed to redirect stderr\n");
|
||||||
|
@ -342,31 +350,36 @@ void identify_trace(char *tlefile,struct trace t,int satno,char *freqlist)
|
||||||
v=(double *) malloc(sizeof(double)*t.n);
|
v=(double *) malloc(sizeof(double)*t.n);
|
||||||
|
|
||||||
// Get observer position
|
// Get observer position
|
||||||
for (i=0;i<t.n;i++)
|
for (i=0;i<t.n;i++)
|
||||||
obspos_xyz(t.mjd[i],s.lng,s.lat,s.alt,&p[i].obspos,&p[i].obsvel);
|
obspos_xyz(t.mjd[i],s.lng,s.lat,s.alt,&p[i].obspos,&p[i].obsvel);
|
||||||
|
|
||||||
printf("Fitting trace:\n");
|
printf("Fitting trace:\n");
|
||||||
|
|
||||||
// Loop over TLEs
|
// Load TLEs
|
||||||
file=fopen(tlefile,"r");
|
tles_t twolines = load_tles(tlefile);
|
||||||
if (file==NULL) {
|
|
||||||
fprintf(stderr,"TLE file %s not found\n",tlefile);
|
if (twolines.number_of_elements == 0) {
|
||||||
|
fprintf(stderr,"TLE file %s not found or empty\n", tlefile);
|
||||||
return;
|
return;
|
||||||
}
|
}
|
||||||
while (read_twoline(file,satno,&orb)==0) {
|
|
||||||
|
for (long elem = 0; elem < twolines.number_of_elements; elem++) {
|
||||||
|
// Get TLE
|
||||||
|
orb = get_orbit_by_index(&twolines, elem);
|
||||||
|
|
||||||
// Initialize
|
// Initialize
|
||||||
imode=init_sgdp4(&orb);
|
imode=init_sgdp4(orb);
|
||||||
if (imode==SGDP4_ERROR) {
|
if (imode==SGDP4_ERROR) {
|
||||||
printf("Error with %d, skipping\n",orb.satno);
|
printf("Error with %d, skipping\n",orb->satno);
|
||||||
continue;
|
continue;
|
||||||
}
|
}
|
||||||
|
|
||||||
// Loop over points
|
// Loop over points
|
||||||
for (i=0,sum1=0.0,sum2=0.0;i<t.n;i++) {
|
for (i=0,sum1=0.0,sum2=0.0;i<t.n;i++) {
|
||||||
// Get satellite position
|
// Get satellite position
|
||||||
satpos_xyz(t.mjd[i]+2400000.5,&satpos,&satvel);
|
satpos_xyz(t.mjd[i]+2400000.5,&satpos,&satvel);
|
||||||
|
|
||||||
dx=satpos.x-p[i].obspos.x;
|
dx=satpos.x-p[i].obspos.x;
|
||||||
dy=satpos.y-p[i].obspos.y;
|
dy=satpos.y-p[i].obspos.y;
|
||||||
dz=satpos.z-p[i].obspos.z;
|
dz=satpos.z-p[i].obspos.z;
|
||||||
dvx=satvel.x-p[i].obsvel.x;
|
dvx=satvel.x-p[i].obsvel.x;
|
||||||
|
@ -383,7 +396,7 @@ void identify_trace(char *tlefile,struct trace t,int satno,char *freqlist)
|
||||||
freq0=sum1/sum2;
|
freq0=sum1/sum2;
|
||||||
|
|
||||||
// Compute residuals
|
// Compute residuals
|
||||||
for (i=0,rms=0.0;i<t.n;i++)
|
for (i=0,rms=0.0;i<t.n;i++)
|
||||||
rms+=pow(t.freq[i]-(1.0-v[i]/C)*freq0,2);
|
rms+=pow(t.freq[i]-(1.0-v[i]/C)*freq0,2);
|
||||||
rms=sqrt(rms/(double) t.n);
|
rms=sqrt(rms/(double) t.n);
|
||||||
|
|
||||||
|
@ -391,16 +404,16 @@ void identify_trace(char *tlefile,struct trace t,int satno,char *freqlist)
|
||||||
for (i=1,mjd0=0.0;i<t.n;i++)
|
for (i=1,mjd0=0.0;i<t.n;i++)
|
||||||
if (v[i]*v[i-1]<0.0)
|
if (v[i]*v[i-1]<0.0)
|
||||||
mjd0=t.mjd[i];
|
mjd0=t.mjd[i];
|
||||||
|
|
||||||
if (mjd0>0.0)
|
if (mjd0>0.0)
|
||||||
mjd2nfd(mjd0,nfd);
|
mjd2nfd(mjd0,nfd);
|
||||||
else
|
else
|
||||||
strcpy(nfd,"0000-00-00T00:00:00");
|
strcpy(nfd,"0000-00-00T00:00:00");
|
||||||
|
|
||||||
if (rms<1000) {
|
if (rms<1000) {
|
||||||
printf("%05d: %s %8.3f MHz %8.3f kHz\n",orb.satno,nfd,1e-6*freq0,1e-3*rms);
|
printf("%05d: %s %8.3f MHz %8.3f kHz\n",orb->satno,nfd,1e-6*freq0,1e-3*rms);
|
||||||
if (flag==0 || rms<rmsmin) {
|
if (flag==0 || rms<rmsmin) {
|
||||||
satnomin=orb.satno;
|
satnomin=orb->satno;
|
||||||
strcpy(nfdmin,nfd);
|
strcpy(nfdmin,nfd);
|
||||||
freqmin=freq0;
|
freqmin=freq0;
|
||||||
rmsmin=rms;
|
rmsmin=rms;
|
||||||
|
@ -408,7 +421,6 @@ void identify_trace(char *tlefile,struct trace t,int satno,char *freqlist)
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
fclose(file);
|
|
||||||
fclose(stderr);
|
fclose(stderr);
|
||||||
|
|
||||||
if (flag==1) {
|
if (flag==1) {
|
||||||
|
@ -432,6 +444,7 @@ void identify_trace(char *tlefile,struct trace t,int satno,char *freqlist)
|
||||||
}
|
}
|
||||||
|
|
||||||
// Free
|
// Free
|
||||||
|
free_tles(&twolines);
|
||||||
free(p);
|
free(p);
|
||||||
free(v);
|
free(v);
|
||||||
|
|
||||||
|
@ -444,12 +457,12 @@ int is_classified(int satno)
|
||||||
int flag=0,no;
|
int flag=0,no;
|
||||||
char *env,tlefile[128],line[LIM];
|
char *env,tlefile[128],line[LIM];
|
||||||
FILE *file;
|
FILE *file;
|
||||||
|
|
||||||
// Get classfd.tle path
|
// Get classfd.tle path
|
||||||
env=getenv("ST_TLEDIR");
|
env=getenv("ST_TLEDIR");
|
||||||
if(env==NULL||strlen(env)==0)
|
if(env==NULL||strlen(env)==0)
|
||||||
env=".";
|
env=".";
|
||||||
sprintf(tlefile,"%s/classfd.tle",env);
|
sprintf(tlefile,"%s/classfd.tle",env);
|
||||||
|
|
||||||
// Does it exist
|
// Does it exist
|
||||||
file=fopen(tlefile,"r");
|
file=fopen(tlefile,"r");
|
||||||
|
@ -467,7 +480,7 @@ int is_classified(int satno)
|
||||||
}
|
}
|
||||||
fclose(file);
|
fclose(file);
|
||||||
}
|
}
|
||||||
|
|
||||||
return flag;
|
return flag;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
@ -478,7 +491,7 @@ struct trace *compute_trace(char *tlefile,double *mjd,int n,int site_id,float fr
|
||||||
struct point *p;
|
struct point *p;
|
||||||
struct site s,sg;
|
struct site s,sg;
|
||||||
FILE *file,*infile;
|
FILE *file,*infile;
|
||||||
orbit_t orb;
|
orbit_t *orb;
|
||||||
xyz_t satpos,satvel;
|
xyz_t satpos,satvel;
|
||||||
double dx,dy,dz,dvx,dvy,dvz,r,v,za,vg;
|
double dx,dy,dz,dvx,dvy,dvz,r,v,za,vg;
|
||||||
double freq0;
|
double freq0;
|
||||||
|
@ -494,7 +507,7 @@ struct trace *compute_trace(char *tlefile,double *mjd,int n,int site_id,float fr
|
||||||
// Reloop stderr
|
// Reloop stderr
|
||||||
if (freopen("/tmp/stderr.txt","w",stderr)==NULL)
|
if (freopen("/tmp/stderr.txt","w",stderr)==NULL)
|
||||||
fprintf(stderr,"Failed to redirect stderr\n");
|
fprintf(stderr,"Failed to redirect stderr\n");
|
||||||
|
|
||||||
// Find number of satellites in frequency range
|
// Find number of satellites in frequency range
|
||||||
infile=fopen(freqlist,"r");
|
infile=fopen(freqlist,"r");
|
||||||
if (infile==NULL) {
|
if (infile==NULL) {
|
||||||
|
@ -503,7 +516,7 @@ struct trace *compute_trace(char *tlefile,double *mjd,int n,int site_id,float fr
|
||||||
return NULL;
|
return NULL;
|
||||||
} else {
|
} else {
|
||||||
for (i=0;;) {
|
for (i=0;;) {
|
||||||
if (fgetline(infile,line,LIM)<=0)
|
if (fgetline(infile,line,LIM)<=0)
|
||||||
break;
|
break;
|
||||||
if (line[0]=='#')
|
if (line[0]=='#')
|
||||||
continue;
|
continue;
|
||||||
|
@ -513,7 +526,7 @@ struct trace *compute_trace(char *tlefile,double *mjd,int n,int site_id,float fr
|
||||||
i++;
|
i++;
|
||||||
else if (freq0>=fmin && freq0<=fmax && graves==0)
|
else if (freq0>=fmin && freq0<=fmax && graves==0)
|
||||||
i++;
|
i++;
|
||||||
|
|
||||||
}
|
}
|
||||||
fclose(infile);
|
fclose(infile);
|
||||||
*nsat=i;
|
*nsat=i;
|
||||||
|
@ -541,16 +554,24 @@ struct trace *compute_trace(char *tlefile,double *mjd,int n,int site_id,float fr
|
||||||
p=(struct point *) malloc(sizeof(struct point)*m);
|
p=(struct point *) malloc(sizeof(struct point)*m);
|
||||||
|
|
||||||
// Get observer position
|
// Get observer position
|
||||||
for (i=0;i<m;i++)
|
for (i=0;i<m;i++)
|
||||||
obspos_xyz(mjd[i],s.lng,s.lat,s.alt,&p[i].obspos,&p[i].obsvel);
|
obspos_xyz(mjd[i],s.lng,s.lat,s.alt,&p[i].obspos,&p[i].obsvel);
|
||||||
|
|
||||||
// Compute Graves positions
|
// Compute Graves positions
|
||||||
if (graves==1) {
|
if (graves==1) {
|
||||||
sg=get_site(9999);
|
sg=get_site(9999);
|
||||||
for (i=0;i<m;i++)
|
for (i=0;i<m;i++)
|
||||||
obspos_xyz(mjd[i],sg.lng,sg.lat,sg.alt,&p[i].grpos,&p[i].grvel);
|
obspos_xyz(mjd[i],sg.lng,sg.lat,sg.alt,&p[i].grpos,&p[i].grvel);
|
||||||
}
|
}
|
||||||
|
|
||||||
|
// Load TLEs
|
||||||
|
tles_t twolines = load_tles(tlefile);
|
||||||
|
|
||||||
|
if (twolines.number_of_elements == 0) {
|
||||||
|
fprintf(stderr,"TLE file %s not found or empty\n", tlefile);
|
||||||
|
return NULL;
|
||||||
|
}
|
||||||
|
|
||||||
infile=fopen(freqlist,"r");
|
infile=fopen(freqlist,"r");
|
||||||
for (j=0;;) {
|
for (j=0;;) {
|
||||||
if (fgetline(infile,line,LIM)<=0)
|
if (fgetline(infile,line,LIM)<=0)
|
||||||
|
@ -562,10 +583,10 @@ struct trace *compute_trace(char *tlefile,double *mjd,int n,int site_id,float fr
|
||||||
flag=1;
|
flag=1;
|
||||||
else if (freq0>=fmin && freq0<=fmax && graves==0)
|
else if (freq0>=fmin && freq0<=fmax && graves==0)
|
||||||
flag=1;
|
flag=1;
|
||||||
|
|
||||||
if (flag==0)
|
if (flag==0)
|
||||||
continue;
|
continue;
|
||||||
|
|
||||||
// Allocate
|
// Allocate
|
||||||
t[j].satno=satno;
|
t[j].satno=satno;
|
||||||
t[j].site=site_id;
|
t[j].site=site_id;
|
||||||
|
@ -576,27 +597,24 @@ struct trace *compute_trace(char *tlefile,double *mjd,int n,int site_id,float fr
|
||||||
t[j].za=(float *) malloc(sizeof(float)*m);
|
t[j].za=(float *) malloc(sizeof(float)*m);
|
||||||
t[j].classfd=is_classified(t[j].satno);
|
t[j].classfd=is_classified(t[j].satno);
|
||||||
t[j].graves=graves;
|
t[j].graves=graves;
|
||||||
|
|
||||||
// Loop over TLEs
|
|
||||||
hastle=0;
|
|
||||||
file=fopen(tlefile,"r");
|
|
||||||
while (read_twoline(file,satno,&orb)==0) {
|
|
||||||
if (orb.satno==satno)
|
|
||||||
hastle=1;
|
|
||||||
|
|
||||||
|
// Get TLE
|
||||||
|
orb = get_orbit_by_catalog_id(&twolines, satno);
|
||||||
|
|
||||||
|
if (orb) {
|
||||||
// Initialize
|
// Initialize
|
||||||
imode=init_sgdp4(&orb);
|
imode=init_sgdp4(orb);
|
||||||
if (imode==SGDP4_ERROR) {
|
if (imode==SGDP4_ERROR) {
|
||||||
printf("Error with %d, skipping\n",orb.satno);
|
printf("Error with %d, skipping\n",orb->satno);
|
||||||
continue;
|
continue;
|
||||||
}
|
}
|
||||||
|
|
||||||
// Loop over points
|
// Loop over points
|
||||||
for (i=0,flag=0,tflag=0;i<m;i++) {
|
for (i=0,flag=0,tflag=0;i<m;i++) {
|
||||||
// Get satellite position
|
// Get satellite position
|
||||||
satpos_xyz(mjd[i]+2400000.5,&satpos,&satvel);
|
satpos_xyz(mjd[i]+2400000.5,&satpos,&satvel);
|
||||||
|
|
||||||
dx=satpos.x-p[i].obspos.x;
|
dx=satpos.x-p[i].obspos.x;
|
||||||
dy=satpos.y-p[i].obspos.y;
|
dy=satpos.y-p[i].obspos.y;
|
||||||
dz=satpos.z-p[i].obspos.z;
|
dz=satpos.z-p[i].obspos.z;
|
||||||
dvx=satvel.x-p[i].obsvel.x;
|
dvx=satvel.x-p[i].obsvel.x;
|
||||||
|
@ -605,7 +623,7 @@ struct trace *compute_trace(char *tlefile,double *mjd,int n,int site_id,float fr
|
||||||
r=sqrt(dx*dx+dy*dy+dz*dz);
|
r=sqrt(dx*dx+dy*dy+dz*dz);
|
||||||
v=(dvx*dx+dvy*dy+dvz*dz)/r;
|
v=(dvx*dx+dvy*dy+dvz*dz)/r;
|
||||||
za=acos((p[i].obspos.x*dx+p[i].obspos.y*dy+p[i].obspos.z*dz)/(r*XKMPER))*R2D;
|
za=acos((p[i].obspos.x*dx+p[i].obspos.y*dy+p[i].obspos.z*dz)/(r*XKMPER))*R2D;
|
||||||
|
|
||||||
// Store
|
// Store
|
||||||
t[j].mjd[i]=mjd[i];
|
t[j].mjd[i]=mjd[i];
|
||||||
t[j].freq[i]=(1.0-v/C)*freq0;
|
t[j].freq[i]=(1.0-v/C)*freq0;
|
||||||
|
@ -613,7 +631,7 @@ struct trace *compute_trace(char *tlefile,double *mjd,int n,int site_id,float fr
|
||||||
|
|
||||||
// Compute Graves velocity/frequency
|
// Compute Graves velocity/frequency
|
||||||
if (graves==1) {
|
if (graves==1) {
|
||||||
dx=satpos.x-p[i].grpos.x;
|
dx=satpos.x-p[i].grpos.x;
|
||||||
dy=satpos.y-p[i].grpos.y;
|
dy=satpos.y-p[i].grpos.y;
|
||||||
dz=satpos.z-p[i].grpos.z;
|
dz=satpos.z-p[i].grpos.z;
|
||||||
dvx=satvel.x-p[i].grvel.x;
|
dvx=satvel.x-p[i].grvel.x;
|
||||||
|
@ -630,22 +648,22 @@ struct trace *compute_trace(char *tlefile,double *mjd,int n,int site_id,float fr
|
||||||
t[j].za[i]=100.0;
|
t[j].za[i]=100.0;
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
}
|
|
||||||
fclose(file);
|
|
||||||
|
|
||||||
// Increment
|
// Increment
|
||||||
if (hastle==1)
|
|
||||||
j++;
|
j++;
|
||||||
|
}
|
||||||
}
|
}
|
||||||
fclose(infile);
|
fclose(infile);
|
||||||
fclose(stderr);
|
fclose(stderr);
|
||||||
|
|
||||||
// Free
|
// Free
|
||||||
|
|
||||||
|
free_tles(&twolines);
|
||||||
free(p);
|
free(p);
|
||||||
|
|
||||||
// Update counter
|
// Update counter
|
||||||
*nsat=j;
|
*nsat=j;
|
||||||
|
|
||||||
return t;
|
return t;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
@ -655,8 +673,8 @@ void compute_doppler(char *tlefile,double *mjd,int n,int site_id,int satno,int g
|
||||||
int i,j,imode,flag,tflag,m,status;
|
int i,j,imode,flag,tflag,m,status;
|
||||||
struct point *p;
|
struct point *p;
|
||||||
struct site s,sg;
|
struct site s,sg;
|
||||||
FILE *file,*outfile;
|
FILE *outfile;
|
||||||
orbit_t orb;
|
orbit_t *orb;
|
||||||
xyz_t satpos,satvel;
|
xyz_t satpos,satvel;
|
||||||
double dx,dy,dz,dvx,dvy,dvz,r,v,rg,vg;
|
double dx,dy,dz,dvx,dvy,dvz,r,v,rg,vg;
|
||||||
double freq0;
|
double freq0;
|
||||||
|
@ -669,7 +687,7 @@ void compute_doppler(char *tlefile,double *mjd,int n,int site_id,int satno,int g
|
||||||
// Reloop stderr
|
// Reloop stderr
|
||||||
if (freopen("/tmp/stderr.txt","w",stderr)==NULL)
|
if (freopen("/tmp/stderr.txt","w",stderr)==NULL)
|
||||||
fprintf(stderr,"Failed to redirect stderr\n");
|
fprintf(stderr,"Failed to redirect stderr\n");
|
||||||
|
|
||||||
// Get site
|
// Get site
|
||||||
s=get_site(site_id);
|
s=get_site(site_id);
|
||||||
|
|
||||||
|
@ -677,13 +695,13 @@ void compute_doppler(char *tlefile,double *mjd,int n,int site_id,int satno,int g
|
||||||
p=(struct point *) malloc(sizeof(struct point)*n);
|
p=(struct point *) malloc(sizeof(struct point)*n);
|
||||||
|
|
||||||
// Get observer position
|
// Get observer position
|
||||||
for (i=0;i<n;i++)
|
for (i=0;i<n;i++)
|
||||||
obspos_xyz(mjd[i],s.lng,s.lat,s.alt,&p[i].obspos,&p[i].obsvel);
|
obspos_xyz(mjd[i],s.lng,s.lat,s.alt,&p[i].obspos,&p[i].obsvel);
|
||||||
|
|
||||||
// Compute Graves positions
|
// Compute Graves positions
|
||||||
if (graves==1) {
|
if (graves==1) {
|
||||||
sg=get_site(9999);
|
sg=get_site(9999);
|
||||||
for (i=0;i<n;i++)
|
for (i=0;i<n;i++)
|
||||||
obspos_xyz(mjd[i],sg.lng,sg.lat,sg.alt,&p[i].grpos,&p[i].grvel);
|
obspos_xyz(mjd[i],sg.lng,sg.lat,sg.alt,&p[i].grpos,&p[i].grvel);
|
||||||
}
|
}
|
||||||
|
|
||||||
|
@ -695,27 +713,32 @@ void compute_doppler(char *tlefile,double *mjd,int n,int site_id,int satno,int g
|
||||||
fprintf(outfile, "# satno mjd r v azi alt rg vg azig altg\n");
|
fprintf(outfile, "# satno mjd r v azi alt rg vg azig altg\n");
|
||||||
else
|
else
|
||||||
fprintf(outfile, "# satno mjd r v azi alt\n");
|
fprintf(outfile, "# satno mjd r v azi alt\n");
|
||||||
|
|
||||||
// Loop over TLEs
|
|
||||||
file=fopen(tlefile,"r");
|
|
||||||
while (read_twoline(file,satno,&orb)==0) {
|
|
||||||
// Initialize
|
|
||||||
imode=init_sgdp4(&orb);
|
|
||||||
if (imode==SGDP4_ERROR) {
|
|
||||||
printf("Error with %d, skipping\n",orb.satno);
|
|
||||||
continue;
|
|
||||||
}
|
|
||||||
|
|
||||||
// Skip high satellites
|
// Load TLEs
|
||||||
if (skiphigh==1 && orb.rev<10.0)
|
tles_t twolines = load_tles(tlefile);
|
||||||
continue;
|
|
||||||
|
if (twolines.number_of_elements == 0) {
|
||||||
|
fprintf(stderr,"TLE file %s not found or empty\n", tlefile);
|
||||||
|
return;
|
||||||
|
}
|
||||||
|
|
||||||
|
// Get TLE
|
||||||
|
orb = get_orbit_by_catalog_id(&twolines, satno);
|
||||||
|
|
||||||
|
// Skip high satellites
|
||||||
|
if (orb && !(skiphigh == 1 && orb->rev < 10.0)) {
|
||||||
|
// Initialize
|
||||||
|
imode=init_sgdp4(orb);
|
||||||
|
if (imode==SGDP4_ERROR) {
|
||||||
|
printf("Error with %d, skipping\n",orb->satno);
|
||||||
|
}
|
||||||
|
|
||||||
// Loop over points
|
// Loop over points
|
||||||
for (i=0,flag=0,tflag=0;i<n;i++) {
|
for (i=0,flag=0,tflag=0;i<n;i++) {
|
||||||
// Get satellite position
|
// Get satellite position
|
||||||
satpos_xyz(mjd[i]+2400000.5,&satpos,&satvel);
|
satpos_xyz(mjd[i]+2400000.5,&satpos,&satvel);
|
||||||
|
|
||||||
dx=satpos.x-p[i].obspos.x;
|
dx=satpos.x-p[i].obspos.x;
|
||||||
dy=satpos.y-p[i].obspos.y;
|
dy=satpos.y-p[i].obspos.y;
|
||||||
dz=satpos.z-p[i].obspos.z;
|
dz=satpos.z-p[i].obspos.z;
|
||||||
dvx=satvel.x-p[i].obsvel.x;
|
dvx=satvel.x-p[i].obsvel.x;
|
||||||
|
@ -726,10 +749,10 @@ void compute_doppler(char *tlefile,double *mjd,int n,int site_id,int satno,int g
|
||||||
ra=modulo(atan2(dy,dx)*R2D,360.0);
|
ra=modulo(atan2(dy,dx)*R2D,360.0);
|
||||||
de=asin(dz/r)*R2D;
|
de=asin(dz/r)*R2D;
|
||||||
equatorial2horizontal(mjd[i],ra,de,s.lng,s.lat,&azi,&alt);
|
equatorial2horizontal(mjd[i],ra,de,s.lng,s.lat,&azi,&alt);
|
||||||
|
|
||||||
// Compute Graves velocity/frequency
|
// Compute Graves velocity/frequency
|
||||||
if (graves==1) {
|
if (graves==1) {
|
||||||
dx=satpos.x-p[i].grpos.x;
|
dx=satpos.x-p[i].grpos.x;
|
||||||
dy=satpos.y-p[i].grpos.y;
|
dy=satpos.y-p[i].grpos.y;
|
||||||
dz=satpos.z-p[i].grpos.z;
|
dz=satpos.z-p[i].grpos.z;
|
||||||
dvx=satvel.x-p[i].grvel.x;
|
dvx=satvel.x-p[i].grvel.x;
|
||||||
|
@ -740,18 +763,18 @@ void compute_doppler(char *tlefile,double *mjd,int n,int site_id,int satno,int g
|
||||||
rag=modulo(atan2(dy,dx)*R2D,360.0);
|
rag=modulo(atan2(dy,dx)*R2D,360.0);
|
||||||
deg=asin(dz/rg)*R2D;
|
deg=asin(dz/rg)*R2D;
|
||||||
equatorial2horizontal(mjd[i],rag,deg,sg.lng,sg.lat,&azig,&altg);
|
equatorial2horizontal(mjd[i],rag,deg,sg.lng,sg.lat,&azig,&altg);
|
||||||
fprintf(outfile,"%05d %14.8lf %f %f %f %f %f %f %f %f\n",orb.satno,mjd[i],r,v,azi,alt,rg,vg,azig,altg);
|
fprintf(outfile,"%05d %14.8lf %f %f %f %f %f %f %f %f\n",orb->satno,mjd[i],r,v,azi,alt,rg,vg,azig,altg);
|
||||||
} else {
|
} else {
|
||||||
fprintf(outfile,"%05d %14.8lf %f %f %f %f\n",orb.satno,mjd[i],r,v,azi,alt);
|
fprintf(outfile,"%05d %14.8lf %f %f %f %f\n",orb->satno,mjd[i],r,v,azi,alt);
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
fclose(file);
|
|
||||||
fclose(outfile);
|
fclose(outfile);
|
||||||
|
|
||||||
fclose(stderr);
|
fclose(stderr);
|
||||||
|
|
||||||
// Free
|
// Free
|
||||||
|
free_tles(&twolines);
|
||||||
free(p);
|
free(p);
|
||||||
|
|
||||||
return;
|
return;
|
||||||
|
|
Loading…
Reference in New Issue