import matplotlib.pyplot as plt
# Monkey patch in fftn and ifftn from pyfftw.interfaces.scipy_fftpack
fftpack.fft2 = pyfftw.interfaces.scipy_fftpack.fft2
fftpack.ifft2 = pyfftw.interfaces.scipy_fftpack.ifft2
fftpack.fft = pyfftw.interfaces.scipy_fftpack.fft
fftpack.ifft = pyfftw.interfaces.scipy_fftpack.ifft
# Turn on the cache for optimum performance
pyfftw.interfaces.cache.enable()
#parameters
N = 256
floatType = np.complex
twoQuads = True
p = nt.nearestPrime(N)
p = N
#-------------------------------
#load kspace data
#from scipy.io import loadmat
#load Cartesian data
#Attention: You must ensure the kspace data is correctly centered or not centered.
image = shepp_logan(N)
kspace = fftpack.fft2(image)
print("kSpace Shape:", kspace.shape)
kMaxValue = np.max(kspace)
kMinValue = np.min(kspace)
print("k-Space Max Value:", kMaxValue)
print("k-Space Min Value:", kMinValue)
from __future__ import print_function # (at top of module)
import _libpath #add custom libs
import finitetransform.numbertheory as nt #local modules
import finitetransform.mojette as mojette
import finitetransform.radon as radon
import finitetransform.imageio as imageio #local module
import finitetransform.farey as farey #local module
import numpy as np
#parameters
N = 512
M = 2 * N
K = 1
twoQuads = True
print("N:", N, "M:", M)
p = nt.nearestPrime(M)
print("p:", p)
pDash = nt.nearestPrime(N)
print("p':", pDash)
#angles = mojette.angleSet_Finite(pDash, 2)
angles, lengths = mojette.angleSet_Symmetric(N, N, 1, True, K)
perpAngle = farey.farey(1, 0)
angles.append(perpAngle)
print("Number of Angles:", len(angles))
print("angles:", angles)
powerSpect = np.zeros((p, p))
#compute lines
print("Computing Finite lines...")
centered = True
# -*- coding: utf-8 -*- """ Test if nearest prime function works in NT Created on Tue Oct 28 10:48:41 2014 @author: uqscha22 """ import _libpath #add custom libs import finitetransform.numbertheory as nt n = 98 print "Nearest prime to", n, "is", nt.nearestPrime(n)