cbassa

# Open file
f=open(fname)

# Read number of channels
nchan=int(np.fromfile(f,dtype='float32',count=1)[0])

# Read channel frequencies
freq=np.fromfile(f,dtype='float32',count=nchan)/1000.0

# Read entire file

cbassa

#include <stdio.h>
#include <cpgplot.h>

int main(int argc,char *argv[])
{

  cpgopen("?");
  cpgend();

  return 0;

cbassa

import numpy as np
from astropy.time import Time
import matplotlib.pyplot as plt

def parse_line(line):
    # Strip line and split
    line = [part.strip() for part in line.strip().split(",")]

    # Get MJD
    t = Time("%sT%s"%(line[0], line[1]), format="isot", scale="utc")

cbassa

20442 437.126562
24278 435.795813
25397 435.225625
25544 437.005969
25544 437.049719
27607 436.793406
27844 436.836844
27848 436.849562
27939 435.354062
28895 437.465969

cbassa

#!/usr/bin/env python3
import numpy as np
import h5py
import argparse
import os

if __name__ == "__main__":
    # Parse arguments
    parser = argparse.ArgumentParser(
        description="Convert a SatNOGS waterfall file to HDF5 format")

cbassa

Object 80000
1 80000U 20001A   20180.00000000  .00000000  00000-0  50000-4 0    04
2 80000  53.0000   0.0000 0001000   0.0000   0.0000 15.05490756    01
Object 80001
1 80001U 20001A   20180.00000000  .00000000  00000-0  50000-4 0    05
2 80001  53.0000   0.0000 0001000   0.0000  16.3636 15.05490756    07
Object 80002
1 80002U 20001A   20180.00000000  .00000000  00000-0  50000-4 0    06
2 80002  53.0000   0.0000 0001000   0.0000  32.7273 15.05490756    07
Object 80003

cbassa

#!/usr/bin/env python3
import json
import datetime
from spacetrack import SpaceTrackClient
import matplotlib.pyplot as plt
import matplotlib.dates as mdates

if __name__ == "__main__":
    # Login
    st = SpaceTrackClient("username", "password")

cbassa

#!/usr/bin/env python3
"""
This test will compare the output of ephem, skyfield and sgp4

Reference information from cbassa/sattools skymap
Time: 2019-10-31T18:22:00
Location: 52.8344N, 6.3785E, 10m
TLE:
USA 276
1 42689U 17022A   19300.77210760 0.00002000  00000-0  26008-4 0    09

cbassa

#!/usr/bin/env python3
import sys
import numpy as np
from astropy.wcs import wcs

w = wcs.WCS(sys.argv[1])

sx = np.sqrt(w.wcs.cd[0, 0]**2+w.wcs.cd[0, 1]**2)*3600.0
sy = np.sqrt(w.wcs.cd[1, 0]**2+w.wcs.cd[1, 1]**2)*3600.0
ang1 = np.arctan2(w.wcs.cd[0, 1], w.wcs.cd[0, 0])*180.0/np.pi

cbassa

#!/usr/bin/env python3
import numpy as np
from sgp4.earth_gravity import wgs84
from sgp4.io import twoline2rv
from astropy.time import Time

class TwoLineElement:
    """TLE class"""

    def __init__(self, tle0, tle1, tle2):