@@ -0,0 +1,934 @@
+#!/usr/bin/env python%0Aimport csv%0Aimport numpy as np%0Afrom scipy import signal%0Afrom scipy.fftpack import fft, ifft%0Afrom scipy import fftpack%0Aimport matplotlib.pyplot as plt%0A%0Anote_name = '77726_MS-DAR-00205-00001-000-00096'%0Abase_dir = '/home/ibanez/data/amnh/darwin_notes/'%0Aimages_dir = base_dir + 'images/'%0Acurves_dir = base_dir + 'image_csvs/'%0Aimage_name = images_dir + note_name + '.jpg'%0Anorth_curve_name = curves_dir + note_name + '_north.csv'%0Asouth_curve_name = curves_dir + note_name + '_south.csv'%0A%0Ainput_north_curve = list(csv.reader(open(north_curve_name)))%0Ainput_south_curve = list(csv.reader(open(south_curve_name)))%0A%0AN = 4096%0Anorth_curve = signal.resample(input_north_curve, N)%0Asouth_curve = signal.resample(input_south_curve, N)%0A%0Ax_signal = %5B%5D%0Ay_signal = %5B%5D%0Afor x, y in north_curve:%0A  y_signal.append(y)%0A  x_signal.append(x)%0A%0Ayfft = fftpack.fftshift(fft(y_signal))%0A%0Aplt.semilogy(2.0/N * np.abs(yfft))%0Aplt.grid()%0Aplt.show()%0A%0A

Add python script to compute fast fourier transforms.