test_skycomponent.py 文件源码

def test_insert_skycomponent_lanczos_bandwidth(self):
        dphasecentre = SkyCoord(ra=+181.0 * u.deg, dec=-58.0 * u.deg, frame='icrs', equinox='J2000')
        sc = create_skycomponent(direction=dphasecentre, flux=numpy.array([[1.0]]), frequency=self.frequency,
                                 polarisation_frame=PolarisationFrame('stokesI'))
        self.model.data *= 0.0
        insert_skycomponent(self.model, sc, insert_method='Lanczos', bandwidth=0.5)
        # These test a regression but are not known a priori to be correct
        self.assertAlmostEqual(self.model.data[0,0,151, 122], 0.24031092091707615, 7)
        self.assertAlmostEqual(self.model.data[0,0,152, 122], 0.18648989466050975, 7)