duanemalcolm

import numpy as np
import shapefile
import sys
from scipy.spatial import cKDTree
import scipy.optimize


def ellipse(x, theta):
    phi = x[4]
    r = 1 / np.sqrt((np.cos(theta)) ** 2 + (np.sin(theta)) ** 2)