[1853] | 1 | import sys |
---|
| 2 | import os |
---|
| 3 | import shutil |
---|
| 4 | import datetime |
---|
| 5 | from asap import scantable, selector, mask_not |
---|
| 6 | from asap.logging import asaplog |
---|
| 7 | asaplog.disable() |
---|
| 8 | |
---|
| 9 | from nose.tools import * |
---|
| 10 | |
---|
| 11 | def tempdir_setup(): |
---|
| 12 | os.makedirs("test_temp") |
---|
| 13 | |
---|
| 14 | def tempdir_teardown(): |
---|
| 15 | shutil.rmtree("test_temp", True) |
---|
| 16 | |
---|
| 17 | class TestNRO(object): |
---|
| 18 | def setup(self): |
---|
| 19 | s = scantable('data/B68test.nro', average=False) |
---|
| 20 | sel = selector() |
---|
| 21 | # make sure this order is always correct - it can be random |
---|
| 22 | sel.set_order(["SCANNO", "POLNO"]) |
---|
| 23 | s.set_selection(sel) |
---|
| 24 | self.st = s.copy() |
---|
| 25 | |
---|
| 26 | def test_init(self): |
---|
| 27 | assert_equal(self.st.nrow(), 36) |
---|
| 28 | assert_equal(self.st.get_fluxunit(), 'K') |
---|
| 29 | assert_equal(self.st.nchan(), 2048) |
---|
| 30 | assert_equal(self.st.nif(), 4) |
---|
| 31 | |
---|
| 32 | def test_spectrum(self): |
---|
| 33 | import numpy |
---|
| 34 | sp=numpy.array(self.st._getspectrum(0)) |
---|
| 35 | assert_equal(sp.max(),1.9524071216583252) |
---|
| 36 | assert_equal(sp.min(),-0.0071961274370551109) |
---|
| 37 | assert_equal(sp.argmax(),1648) |
---|
| 38 | assert_equal(sp.argmin(),1614) |
---|
| 39 | |
---|
| 40 | def test_stats(self): |
---|
| 41 | assert_equal(self.st.stats('rms')[0],1.1064267158508301) |
---|
| 42 | assert_equal(self.st.stats('stddev')[0],0.29552212357521057) |
---|
| 43 | assert_equal(self.st.stats('median')[0],1.0714811086654663) |
---|
| 44 | assert_equal(self.st.stats('mean')[0],1.0662506818771362) |
---|
| 45 | assert_equal(self.st.stats('sum')[0],2183.681396484375) |
---|
| 46 | |
---|
| 47 | def test_frequency(self): |
---|
| 48 | rf=self.st.get_restfreqs() |
---|
| 49 | assert_equal(len(rf),2) |
---|
| 50 | assert_equal(rf[0][0],85162157000.0) |
---|
| 51 | assert_equal(rf[1][0],86754330000.0) |
---|
| 52 | self.st.set_unit('GHz') |
---|
| 53 | abc=self.st._getabcissa(0) |
---|
| 54 | assert_equal(abc[0],85.17337639557438) |
---|