Reading RSF (Madagascar) files into numpy arrays.

read_rsf.py

# Example modified from here: https://github.com/ahay/src/issues/130
from subprocess import run
 
import m8r
import numpy as np

# create example data
run('sfspike n1=5 n2=5 k1=3 > spike.rsf', shell=True)
 
# Read rsf file
rsf = m8r.Input('spike.rsf')
 
# Init output array of appropriate size
output = np.zeros(rsf.shape(), dtype='float32')

# Populate output array with contents of rsf file
rsf.read(output)

Actually, I just noticed m8r.Input has a read method, so this can be simplified:

# Example modified from here: https://github.com/ahay/src/issues/130
from subprocess import run
 
import m8r

# create example data
run('sfspike n1=5 n2=5 k1=3 > spike.rsf', shell=True)
 
# Read rsf file
output = m8r.Input('spike.rsf').read()