Beispiel #1
0
 def test_raises(self) :
     Data = self.blocks[0]
     rotate_pol.rotate(Data, (1,))
     self.assertRaises(ce.DataError, pol_cal.calibrate_pol, Data,self.mueler)
     Data = self.blocks[1]
     rotate_pol.rotate(Data, (1,2,3,4))
     Data.field['CRVAL4'][0] = 3
     Data.field['CRVAL4'][2] = 1
     self.assertRaises(ce.DataError, pol_cal.calibrate_pol, Data,
                       self.mueler)
Beispiel #2
0
    def setUp(self) :
        Reader = fitsGBT.Reader(test_file, feedback=0)
        blocks = Reader.read((),())
        for Data in blocks :
            rotate_pol.rotate(Data, (1,2,3,4))
        Writer = fitsGBT.Writer(blocks, feedback=0)
        Writer.write('test_rot.testout.fits')

        self.params = {'pc_file_middles' : ("test",),
                       'pc_input_end' : "_rot.testout.fits",
                       'pc_output_end' : "_polcal.testout.fits",
                       'pc_mueler_file' : test_mueler_file
                       }
    def setUp(self):
        Reader = fitsGBT.Reader(test_file, feedback=0)
        self.blocks = Reader.read((), ())
        b = self.blocks[0]
        b.calc_freq()
        map = sp.zeros((150, 10, 11))
        map = algebra.make_vect(map, axis_names=('freq', 'ra', 'dec'))
        map.set_axis_info('freq', 700.0e6, -2.0e6)
        map.set_axis_info('ra', 325.6, -0.2)
        map.set_axis_info('dec', 0, 0.2)
        self.map = map

        for Data in self.blocks:
            rotate_pol.rotate(Data, (1, 2, 3, 4))
    def setUp(self) :
        Reader = fitsGBT.Reader(test_file, feedback=0)
        self.blocks = Reader.read((),())
        b = self.blocks[0]
        b.calc_freq()
        map = sp.zeros((150, 10, 11))
        map = algebra.make_vect(map, axis_names=('freq', 'ra', 'dec'))
        map.set_axis_info('freq', 700.0e6, -2.0e6)
        map.set_axis_info('ra', 325.6, -0.2)
        map.set_axis_info('dec', 0, 0.2)
        self.map = map

        for Data in self.blocks :
            rotate_pol.rotate(Data, (1, 2, 3, 4))
    def setUp(self) :
        Reader = fitsGBT.Reader(test_file, feedback=0)
        blocks = Reader.read((),())
        for Data in blocks :
            rotate_pol.rotate(Data, (1,2,3,4))
        Writer = fitsGBT.Writer(blocks, feedback=0)
        Writer.write('test_rot.testout.fits')
        
        map = sp.zeros((150, 10, 11))
        map = algebra.make_vect(map, axis_names=('freq', 'ra', 'dec'))
        map.set_axis_info('freq', 700.0e6, -2.0e6)
        map.set_axis_info('ra', 325.6, -0.2)
        map.set_axis_info('dec', 0, 0.2)
        algebra.save('test_inmap.npy', map)

        self.params = {'sm_file_middles' : ("test",),
                       'sm_input_end' : "_rot.testout.fits",
                       'sm_output_end' : "_sub.testout.fits",
                       'sm_solve_for_gain' : True,
                       'sm_gain_output_end' : '_gain.pickle',
                       'sm_map_file' : 'test_inmap.npy'
                       }
Beispiel #6
0
 def test_runs(self) :
     for Data in self.blocks :
         rotate_pol.rotate(Data, (1,2,3,4))
         pol_cal.calibrate_pol(Data, self.mueler)
 def setUp(self) :
     Reader = fitsGBT.Reader(test_file, feedback=0)
     self.Data = Reader.read(1,0)
     rotate_pol.rotate(self.Data, (1,))
     self.DataSub = Reader.read(1,0)
     rotate_pol.rotate(self.DataSub, (1,))
Beispiel #8
0
 def setUp(self):
     Reader = fitsGBT.Reader(test_file, feedback=0)
     self.Data = Reader.read(1, 0)
     rotate_pol.rotate(self.Data, (1, ))
     self.DataSub = Reader.read(1, 0)
     rotate_pol.rotate(self.DataSub, (1, ))