fubel

'''
This script is an implementation of Behrend's algorithm for set coloring.
Here we use it for communication complexity:

The question is, if the numbers [x,y,z] add up to n.
You can set n as well as the numbers
'''

import math
import numpy as np 

fubel

import numpy as np
import matplotlib.pyplot as plt


def fourier_polynomial(f: 'function', n: int, r: int, M: int) -> np.array:
    """ Computes the Fourier polynomial of degree n.
    
    Args:
        f: Any Python function that returns a float
        n: Degree of the Fourier polynomial

fubel

""" Radon Transform as described in Birkfellner, Wolfgang. Applied Medical Image Processing: A Basic Course. [p. 344] """
from scipy import misc
import numpy as np
import matplotlib.pyplot as plt

def discrete_radon_transform(image, steps):
    R = np.zeros((steps, len(image)), dtype='float64')
    for s in range(steps):
        rotation = misc.imrotate(image, -s*180/steps).astype('float64')
        R[:,s] = sum(rotation)