データのないピクセルのヘルピー マップをマスクしようとしていますが、でヘルピー マスクを適用するとbadval=-1.6375e+30
、マスクの周囲に境界線が表示されます。
import healpy as hp
import numpy as np
# load the mask
mask = hp.read_map('mask_nvss_S20-S1000_Ns64.fits')
# degrade the mask to NSIDE = 16
mask16 = np.round(hp.ud_grade(mask, 16))
# apply it to the map
masked_map = hp.ma(map1)
masked_map.mask = np.logical_not(mask16)
# masked map
nvss_map = np.round(masked_map)
この境界効果を取り除くのを手伝ってくれる人はいますか?