python将缺失值用0代替_平滑而不使用零填充缺失值

此问题与以下问题和答案有关(免责声明：来自我)：

它可以按如下方式转移到您的案例中：import numpy as np

import healpy as hp

nside = 128

npix = hp.nside2npix(nside)

# using random numbers here to see the actual smoothing

arr = np.random.rand(npix)

mask = np.zeros(npix, dtype=bool)

mask[:mask.size//2] = True

def masked_smoothing(U, rad=5.0):

V=U.copy()

V[U!=U]=0

VV=hp.smoothing(V, fwhm=np.radians(rad))

W=0*U.copy()+1

W[U!=U]=0

WW=hp.smoothing(W, fwhm=np.radians(rad))

return VV/WW

# setting array to np.nan to exclude the pixels in the smoothing

arr[~mask] = np.nan

arr_sm = masked_smoothing(arr)

arr_sm[~mask] = hp.UNSEEN

hp.mollview(arr, title='Input array')

hp.mollview(arr_sm, title='Smoothed array')