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
3 Sep 30, 2022
7 Sep 24, 2022
3 Dec 20, 2021
7 Dec 01, 2022
17 Nov 13, 2022
16 Nov 30, 2022
4 Dec 09, 2022
5 Feb 19, 2022
1 Dec 28, 2021
27 Nov 01, 2022
16 Nov 26, 2022
258 Sep 20, 2021
13 Jul 11, 2022
4 Jan 15, 2022
2 Jan 28, 2022
7 Jul 31, 2021
30 Dec 04, 2022
68 Jan 03, 2023
33 Nov 07, 2022
49 Feb 10, 2022