Skip to content

Instantly share code, notes, and snippets.

View cbassa's full-sized avatar

Cees Bassa cbassa

  • ASTRON Netherlands Institute for Radio Astronomy
  • Netherlands
View GitHub Profile
@cbassa
cbassa / plot_waterfall.py
Created November 27, 2020 08:20
Waterfall Plotting
#!/usr/bin/env python3
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.gridspec as gridspec
from matplotlib.colorbar import Colorbar
import sys
import warnings
with warnings.catch_warnings():
warnings.filterwarnings("ignore",category=FutureWarning)
import h5py
@cbassa
cbassa / preobs_uhf.sh
Created November 29, 2020 21:15
SatNOGS STRF pre-obs script
#!/usr/bin/bash
echo $1 $2 $3 $4 $5 $6
# Run strf
echo "Killing left over netcat"
ps aux | grep -e "nc" | grep -e "45360" | awk '{printf("kill -9 %s\n",$2)}' | sh
freq=`echo $3 | awk '{printf("%d\n",$1-100000)}'`
fmin=`echo $3 | awk '{printf("%d\n",$1-100000)}'`
fmax=`echo $3 | awk '{printf("%d\n",$1+100000)}'`
@cbassa
cbassa / parse_tle.py
Created July 24, 2021 19:56
Simple TLE parsing script
#!/usr/bin/env python3
from sgp4.api import Satrec
from sgp4.earth_gravity import wgs84
if __name__ == "__main__":
with open("iss.txt", "r") as fp:
lines = fp.readlines()
sats = []
for i in range(len(lines) - 1):
@cbassa
cbassa / plot_hydrogen.py
Created July 27, 2021 15:27
Hydrogen plotting script
#!/usr/bin/env python3
import os
import re
import glob
import tqdm
from datetime import datetime
import numpy as np
import matplotlib.pyplot as plt
#!/usr/bin/env python3
import os
import json
from spacetrack import SpaceTrackClient
from sgp4.earth_gravity import wgs84
from sgp4.io import twoline2rv
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.dates as mdates
from astropy.time import Time
@cbassa
cbassa / tle_swap.py
Created November 16, 2021 20:42
Plot TLE parameters to find swapped objects
#!/usr/bin/env python3
import os
import json
from spacetrack import SpaceTrackClient
from sgp4.earth_gravity import wgs84
from sgp4.io import twoline2rv
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.dates as mdates
from astropy.time import Time
@cbassa
cbassa / sanitize_csv.py
Created April 8, 2022 13:58
Sanitize Bobcat-1 GPS measurements
#!/usr/bin/env python3
import sys
import numpy as np
import matplotlib.pyplot as plt
from astropy.io import ascii
if __name__ == "__main__":
d = ascii.read(sys.argv[1])
@cbassa
cbassa / get.sh
Created August 11, 2022 20:16
gphoto2 DSLR image capture and FITS conversion
#!/bin/bash
# Wait for nearest full second
waitfor
# Log date
timestamp=`date -u +%FT%T`
echo $timestamp >>log.txt
# Capture
@cbassa
cbassa / find.sh
Created August 11, 2022 20:17
Star finding script
#!/bin/bash
# Copy first file
if [ ! -e "test.fits" ]; then
ls -1 2*.fits | head -n1 | awk '{printf("cp %s test.fits\n",$1)}' | sh
fi
# Loop over files
for file in *.fits; do
@cbassa
cbassa / id.sh
Last active August 11, 2022 20:18
Astrometric calibration and prediction
wcsfit
for file in 2*.fits; do
if [ ! -e $file.cal ]; then
addwcs -f $file -r test.fits -m 9 -R 20
satid $file $file.png/png 2>/dev/null
fi
done