Simple test¶
Ensure your device works with this simple test.
examples/mag_cal_simpletest.py¶
1# SPDX-FileCopyrightText: Copyright (c) 2022 Phil Underwood for Underwood Underground
2#
3# SPDX-License-Identifier: Unlicense
4from mag_cal.calibration import Calibration
5from mag_cal.utils import read_fixture
6
7PATH = "../tests/fixtures/cal_data/hj2.json"
8
9with open(PATH) as f:
10 aligned, grav, mag = read_fixture(f.read())
11
12calib = Calibration()
13calib.fit_ellipsoid(mag, grav)
14calib.fit_to_axis(aligned)
15calib.fit_non_linear_quick(aligned, param_count=5)
16
17# calib.fit_non_linear(aligned, param_count=3)
18for m, g in zip(mag, grav):
19 azimuth, inclination, roll = calib.get_angles(m, g)
20 print(f"{azimuth:05.1f}° {inclination:+05.1f}° {roll:+04.0f}°")