ilt.py 文件源码

def mask_init(self):
        x = np.linspace(-10,10,21)
        X, Y = np.meshgrid(x,x)
        R = X**2 + Y**2
        O = np.exp(-R/2/(4**2))
        OO = O/np.sum(O)
        D = sg.fftconvolve(1.0*self.image.mask.data+0.0, OO,'same')
        # D = pyfftw.interfaces.scipy_fftpack.convolve(1.0*self.image.mask.data+0.0, OO,'same')

        self.target = copy.deepcopy(self.image.mask.data)
        self.maskdata = 0.99*D + 0.01
        AA = 2*self.maskdata - 1
        AA = np.complex64(AA)
        BB = np.arccos(AA)
        self.masktheta = BB.real

        self.image.mask.data = self.maskdata