39 from scipy.ndimage
import interpolation
as interp
41 from .
import utilities
as util
58 print(
"Opening apodizer")
59 print(
"reading file:", filename)
60 pup = np.load(filename)
64 raise ValueError(
"Apodizer dimensions must be smaller.")
68 print(
"TODO pup=bilinear(pup,pupd,pupd)")
72 print(
"TODO pup=rotate2(pup,angle)")
73 pup = interp.rotate(pup, angle, reshape=
False, order=2)
75 reg = np.where(util.dist(pupd) > pupd / 2.)
78 pupf = np.zeros((dim, dim), dtype=np.float32)
81 if ((dim - pupd) % 2 != 0):
82 pupf[(dim - pupd + 1) / 2:(dim + pupd + 1) / 2, (dim - pupd + 1) /
83 2:(dim + pupd + 1) / 2] = pup
86 pupf[(dim - pupd) / 2:(dim + pupd) / 2, (dim - pupd) / 2:(dim + pupd) /
92 pupf = np.abs(pupf).astype(np.float32)