« ctf2st | Software | StockwellDs »

The Stockwell Transform

This is a C implementation of the Stockwell Transform. The ST produces a time-frequency representation, similar to a wavelet decomposition.

• st.c — The Stockwell Transform
• stmodule.c — Python wrapper for st()
• smt.py — Python/numpy code for multitaper Stockwell transforms
• sinemodule.c — Sine tapers for smt.py

The code here has been updated to work with fftw3.

To compile the Python wrappers, you'll want to set up a Makefile similar to this:

PYINCDIR = /usr/include/python2.4
NUMPYDIR = /usr/lib/python2.4/site-packages/numpy/core/include/numpy
CFLAGS = -O3 -I$(PYINCDIR) -I$(NUMPYDIR) -fPIC

all: sinemodule.so stmodule.so

sinemodule.so: sinemodule.o
        $(CC) -shared -o $@ sinemodule.o

stmodule.so: stmodule.o st.o
        $(CC) -shared -o $@ stmodule.o st.o -lfftw3