Magnetic field calculator
Python API for British Geological Survey magnetic field calculator.
Description
This project magnetic field calculator. It uses is British Geological Survey (BGS) API web service for calculation.
The web service makes the World Magnetic Model (WMM), the International Geomagnetic Reference Field (IGRF) and the BGS Global Geomagnetic Model (BGGM) available as a web service. The IGRF and WMM have no restrictions on use, the BGGM is only available to subscribers. The API provides options to select which model and revision to use. Values of the magnetic field at any point around the world can be retrieved for a given date.
The project is not part of the BGS. It is just API client for it, made by community. For any information about BGS, you should use the official website.
Usage
First, you need to import the calculator class:
from magnetic_field_calculator import MagneticFieldCalculator
Then you need to init the calculator client, with default WMM model and it's latest available version:
calculator = MagneticFieldCalculator()
Model, it's revisions and custom URS can be changed with parameters:
calculator = MagneticFieldCalculator(
model='wmm',
revision='2020',
sub_revision='2',
custom_url='https://example.com'
)
You can then get the calculations for given location:
result = calculator.calculate(
latitude=-80,
longitude=140,
altitude=200,
date='2028-12-31'
)
The only required parameters are latitude and longitude, but it is reccomended to also set others for best results. Some of them are exclusive and you should not use them at same time. You can also set username and password for HTTP auth for protected models. For more details about parameters, see the official documentation.
Output is returned as dictionary:
field_value = result['field-value']
declination = field_value['declination']
inclination = field_value['inclination']
total_intensity = field_value['total-intensity']
north_intensity = field_value['north-intensity']
east_intensity = field_value['east-intensity']
vertical_intensity = field_value['vertical-intensity']
horizontal_intensity = field_value['horizontal-intensity']
Some output properties will have units and value sub properties. Output format will depend on used input, specially for height and date.
You can also look to example file for more examples.
Versioning
This library uses SemVer for versioning. For the versions available, see the tags on this repository.
License
This library is licensed under the GPLv3+ license. See the LICENSE file for details.
168 May 17, 2022
904 Jan 2, 2023
20 Jun 18, 2022
16 Mar 29, 2022
125 Dec 26, 2022
11 Nov 22, 2022
46 Sep 1, 2022
1 Nov 19, 2022
21 May 31, 2022
14 Oct 21, 2022
6 Nov 24, 2022
1 Dec 28, 2021
7 Jul 04, 2022
9 Feb 17, 2022
36 Dec 09, 2022
3 Jul 03, 2022
8 Jan 30, 2022
8 Dec 31, 2022
2 Jan 20, 2022
34 Dec 25, 2022
14 Oct 06, 2022
4 Jan 27, 2022
12 Dec 27, 2022
1 Dec 13, 2021
22 Nov 12, 2022
32 Oct 01, 2022
65 Nov 09, 2022
2 Mar 18, 2022
4 Jan 30, 2022