Skip to content

Instantly share code, notes, and snippets.

View mhvk's full-sized avatar

Marten van Kerkwijk mhvk

  • University of Toronto
  • Toronto, Ontario, Canada
View GitHub Profile
@mhvk
mhvk / drao_vdif.py
Created June 1, 2017 19:48
DRAO vdif conversion script (quite a hack)
from baseband import vdif
from baseband.vdif.header import VDIFHeader0, HeaderParser, VDIFHeader
import numpy as np
import astropy.units as u
files = list('{:07d}.dat'.format(no) for no in range(100000, 100001))
class DRAOVDIFHeader(VDIFHeader0):
"""Test of the scintillometry dedispersion routines on PSR B1937+21"""
import numpy as np
from numpy.lib.format import open_memmap
import astropy.units as u
from baseband import guppi
from scintillometry.dispersion import Dedisperse
from scintillometry.functions import Square
from scintillometry.fourier import get_fft_maker
from scintillometry.integration import Stack
import matplotlib.pylab as plt
@mhvk
mhvk / rdf25b.c
Created January 3, 2020 18:07
C code to convert RadioAstron RDF to Mark 5B
#include <math.h>
#include <string.h>
#include <stdlib.h>
#include <time.h>
#include <stdio.h>
#include <sys/types.h>
#define _FILE_OFFSET_BITS 64
//much of this taken from from Ari Mujunen vsib code
static unsigned char bytetable[256];
unsigned char crcdata[6];