median_fiter_hpss.py 文件源码

def process(self, S):
        mag = S.magnitude()        
        H = np.empty_like(mag)
        H[:] = signal.medfilt(mag, self.h_kernel)
        P = np.empty_like(mag)
        P[:] = signal.medfilt(mag, self.p_kernel)
        h_mask = self.mask_class(H, P, self.mask_exp)
        p_mask = self.mask_class(P, H, self.mask_exp)
        return (S * h_mask, S * p_mask)