def test_process_reverse_polarity(): """algorithm_test.AdjustedAlgorithm_test.test_process() Check adjusted data processing versus files generated from original script. Tests reverse polarity martix. """ # load adjusted data transform matrix and pier correction a = adj( statefile="etc/adjusted/adjbou_state_HE_.json", inchannels=["H", "E"], outchannels=["H", "E"], ) # load boulder May 20 files from /etc/ directory with open("etc/adjusted/BOU202005vmin.min") as f: raw = i2.IAGA2002Factory().parse_string(f.read()) with open("etc/adjusted/BOU202005adj.min") as f: expected = i2.IAGA2002Factory().parse_string(f.read()) # process he(raw) channels with loaded transform adjusted = a.process(raw) # compare channels from adjusted and expected streams assert_almost_equal( actual=adjusted.select(channel="H")[0].data, desired=expected.select(channel="H")[0].data, decimal=2, ) assert_almost_equal( actual=adjusted.select(channel="E")[0].data, desired=expected.select(channel="E")[0].data, decimal=2, )
def test_construct(): """algorithm_test.AdjustedAlgorithm_test.test_construct()""" # load adjusted data transform matrix and pier correction a = adj(statefile="etc/adjusted/adjbou_state_.json") assert_almost_equal(actual=a.matrix[0, 0], desired=9.83427577e-01, decimal=6) assert_equal(actual=a.pier_correction, desired=-22)
def test_construct(): """algorithm_test.AdjustedAlgorithm_test.test_construct() """ matrix = None pier_correction = None # load adjusted data transform matrix and pier correction a = adj(matrix, pier_correction, "etc/adjusted/adjbou_state_.json") assert_almost_equal(a.matrix[0, 0], 9.83427577e-01, 6) assert_equal(a.pier_correction, -22)
def test_construct(): """algorithm_test.AdjustedAlgorithm_test.test_construct() """ matrix = None pier_correction = None # load adjusted data transform matrix and pier correction a = adj(matrix, pier_correction, 'etc/adjusted/adjbou_state_.json') assert_almost_equals(a.matrix[0, 0], 9.83427577e-01, 6) assert_equals(a.pier_correction, -22)
def test_construct(): """ path to """ matrix = None pier_correction = None # load adjusted data transform matrix and pier correction a = adj(matrix, pier_correction, 'etc/adjusted/adjbou_state_.json') assert_almost_equals(a.matrix[0, 0], 9.83427577e-01, 6) assert_equals(a.pier_correction, -22)
def test_process_XYZF(): """algorithm_test.AdjustedAlgorithm_test.test_process() Check adjusted data processing versus files generated from original script """ # load adjusted data transform matrix and pier correction a = adj(statefile="etc/adjusted/adjbou_state_.json") # load boulder Jan 16 files from /etc/ directory with open("etc/adjusted/BOU201601vmin.min") as f: raw = i2.IAGA2002Factory().parse_string(f.read()) with open("etc/adjusted/BOU201601adj.min") as f: expected = i2.IAGA2002Factory().parse_string(f.read()) # process hezf (raw) channels with loaded transform adjusted = a.process(raw) # compare channels from adjusted and expected streams assert_almost_equal( actual=adjusted.select(channel="X")[0].data, desired=expected.select(channel="X")[0].data, decimal=2, ) assert_almost_equal( actual=adjusted.select(channel="Y")[0].data, desired=expected.select(channel="Y")[0].data, decimal=2, ) assert_almost_equal( actual=adjusted.select(channel="Z")[0].data, desired=expected.select(channel="Z")[0].data, decimal=2, ) assert_almost_equal( actual=adjusted.select(channel="F")[0].data, desired=expected.select(channel="F")[0].data, decimal=2, )
def test_process(): """algorithm_test.AdjustedAlgorithm_test.test_process() Check adjusted data processing versus files generated from original script """ matrix = None pier_correction = None # load adjusted data transform matrix and pier correction a = adj(matrix, pier_correction, 'etc/adjusted/adjbou_state_.json') # load boulder Jan 16 files from /etc/ directory hezf_iaga2002_file = open('etc/adjusted/BOU201601vmin.min') hezf_iaga2002_string = hezf_iaga2002_file.read() xyzf_iaga2002_file = open('etc/adjusted/BOU201601adj.min') xyzf_iaga2002_string = xyzf_iaga2002_file.read() factory = i2.IAGA2002Factory() hezf = factory.parse_string(hezf_iaga2002_string) xyzf = factory.parse_string(xyzf_iaga2002_string) # process hezf (raw) channels with loaded transform adj_bou = a.process(hezf) # unpack channels from loaded adjusted data file x = xyzf.select(channel='X')[0] y = xyzf.select(channel='Y')[0] z = xyzf.select(channel='Z')[0] f = xyzf.select(channel='F')[0] # unpack channels from adjusted processing of raw data x_adj = adj_bou.select(channel='X')[0] y_adj = adj_bou.select(channel='Y')[0] z_adj = adj_bou.select(channel='Z')[0] f_adj = adj_bou.select(channel='F')[0] for r in range(hezf[0].data.size): assert_almost_equals(x.data[r], x_adj.data[r], 1) assert_almost_equals(y.data[r], y_adj.data[r], 1) assert_almost_equals(z.data[r], z_adj.data[r], 1) assert_almost_equals(f.data[r], f_adj.data[r], 1)