def smooth(arr, sigma=5):
filled = arr.copy()
nans = sp.isnan(arr)
filled[nans] = 0
smoothed = sp.ndimage.gaussian_filter(filled, sigma=sigma, truncate=10)
missing_coefs = sp.ndimage.gaussian_filter((~nans).astype(float), sigma=sigma, truncate=10)
corrected = smoothed/missing_coefs
return corrected
评论列表
文章目录