pyBADA 0.1.0__py3-none-any.whl

pyBADA/configuration.py

@@ -0,0 +1,98 @@
+# -*- coding: utf-8 -*-
+"""
+pyBADA
+Common configuration module
+Developped @EUROCONTROL (EIH)
+2024
+"""
+
+__author__ = "Henrich Glaser-Opitz"
+__copyright__ = "Copyright 2024, EUROCONTROL (EIH)"
+__license__ = "BADA Eurocontrol"
+__version__ = "1.0.0"
+__maintainer__ = "Henrich Glaser-Opitz"
+__email__ = "henrich.glaser-opitz@eurocontrol.int"
+__status__ = "Development"
+__docformat__ = "reStructuredText"
+
+
+import os
+import sys
+from pathlib import Path
+import inspect
+
+# get current and parent directory
+currentDir = Path(os.path.dirname(os.path.realpath(__file__)))
+parentDir = currentDir.parent
+
+
+def getVersionsList(badaFamily):
+    # list file and directories
+    path = getBadaFamilyPath(badaFamily)
+    items = os.listdir(path)
+
+    versionsList = []
+    for item in items:
+        if os.path.isdir(os.path.join(path, item)):
+            versionsList.append(item)
+
+    return versionsList
+
+
+def getAircraftList(badaFamily, badaVersion):
+    path = getBadaVersionPath(badaFamily, badaVersion)
+
+    if not os.path.exists(path):
+        return []
+    else:
+        items = os.listdir(path)
+
+    # check if I have BADA3 xml or standard ACSII files
+    xml = True
+    for item in items:
+        if "OPF" in item:
+            xml = False
+            break
+
+    if badaFamily == "BADA4" or badaFamily == "BADAH" or badaFamily == "BADAE" or xml:
+        aircraftList = []
+        for item in items:
+            if os.path.isdir(os.path.join(path, item)):
+                aircraftList.append(item)
+
+    elif badaFamily == "BADA3":
+        aircraftList = []
+        for item in items:
+            if len(item.split(".")) == 2:
+                if item.split(".")[0].rstrip("_") not in aircraftList and (
+                    item.split(".")[1] == "PTD"
+                    or item.split(".")[1] == "PTF"
+                    or item.split(".")[1] == "OPF"
+                    or item.split(".")[1] == "APF"
+                ):
+                    aircraftList.append(item.split(".")[0].rstrip("_"))
+
+    return aircraftList
+
+
+def getBadaFamilyPath(badaFamily):
+    path = os.path.join(getAircraftPath(), badaFamily)
+    return path
+
+
+def getBadaVersionPath(badaFamily, badaVersion):
+    path = os.path.join(getAircraftPath(), badaFamily, badaVersion)
+    return path
+
+
+def getAircraftPath():
+    # path = os.path.join(parentDir, "aircraft")
+    package_dir = os.path.dirname(__file__)
+    path = os.path.join(package_dir, "aircraft")
+    return path
+
+
+def getDataPath():
+    package_dir = os.path.dirname(__file__)
+    path = os.path.join(package_dir, "data")
+    return path

pyBADA/constants.py

@@ -0,0 +1,64 @@
+# -*- coding: utf-8 -*-
+"""
+pyBADA
+Common constants module
+Developped @EUROCONTROL (EIH)
+2024
+"""
+
+__author__ = "Henrich Glaser-Opitz"
+__copyright__ = "Copyright 2024, EUROCONTROL (EIH)"
+__license__ = "BADA Eurocontrol"
+__version__ = "1.0.0"
+__maintainer__ = "Henrich Glaser-Opitz"
+__email__ = "henrich.glaser-opitz@eurocontrol.int"
+__status__ = "Development"
+__docformat__ = "reStructuredText"
+
+# Gravitational acceleration
+g = 9.80665
+
+# Real gas constant for air
+R = 287.05287
+
+# ISA temperature gradient with altitude below the tropopause
+temp_h = 0.0065
+
+# Standard atmospheric temperature at MSL
+temp_0 = 288.15
+
+# Standard atmospheric temperature at tropopause (11km)
+temp_11 = 216.65
+
+# Tropopause geopotatial pressure altitude [m]
+h_11 = 11000.0
+
+# Standard atmospheric pressure at tropopause (11km)
+p_11 = 22632.04
+
+# Standard atmospheric pressure at MSL
+p_0 = 101325.0
+
+# Standard atmospheric density at MSL
+rho_0 = 1.225
+
+# Speed of sound
+a_0 = 340.294
+
+# Adiabatic index of air
+Agamma = 1.4
+
+# Adiabatic index of air - ratio
+Amu = (Agamma - 1) / Agamma
+
+# Average Earth radius [km] - https://en.wikipedia.org/wiki/Earth_radius#Mean_radius
+AVG_EARTH_RADIUS_KM = 6371
+
+# semi-major axis (WGS-84) [m]
+a = 6378137.0
+
+# semi-minor axis (WGS-84) [m]
+b = 6356752.314245
+
+# flattening (WGS-84)
+f = (a - b) / a

pyBADA/conversions.py

@@ -0,0 +1,197 @@
+# -*- coding: utf-8 -*-
+"""
+pyBADA
+Common unit conversions module
+Developped @EUROCONTROL (EIH)
+2024
+"""
+
+__author__ = "Henrich Glaser-Opitz"
+__copyright__ = "Copyright 2024, EUROCONTROL (EIH)"
+__license__ = "BADA Eurocontrol"
+__version__ = "1.0.0"
+__maintainer__ = "Henrich Glaser-Opitz"
+__email__ = "henrich.glaser-opitz@eurocontrol.int"
+__status__ = "Development"
+__docformat__ = "reStructuredText"
+
+
+from math import pi
+from datetime import datetime
+from time import mktime
+
+
+def ft2m(val):
+    """
+    This function converts from ft to m s
+
+    :param val: value in ft
+    :returns: vaue in m
+    """
+    return round(float(val) * 0.3048, 10)
+
+
+def nm2m(val):
+    """
+    This function converts from nautical miles to m
+
+    :param val: value in nautical miles
+    :returns: vaue in m
+    """
+    return val * 1852.0
+
+
+def h2s(val):
+    """
+    This function converts from hours to m seconds
+
+    :param val: value in hours
+    :returns: vaue in seconds
+    """
+    return val * 3600.0
+
+
+def kt2ms(val):
+    """
+    This function converts from kt to m s^-1
+
+    :param val: value in kt
+    :returns: vaue in m s^-1
+    """
+
+    if val is None:
+        return None
+    else:
+        return round(float(val) * 0.514444, 10)
+
+
+def lb2kg(val):
+    """
+    This function converts from lb to kg
+
+    :param val: value in lb
+    :returns: vaue in kg
+    """
+    return val * 0.453592
+
+
+def deg2rad(val):
+    """
+    This function converts from decimal degrees to radians
+
+    :param val: value in decimal degrees
+    :returns: vaue in radians
+    """
+    return val * pi / 180.0
+
+
+def m2ft(val):
+    """
+    This function converts from meters to feets
+
+    :param val: value in meters
+    :returns: value in feets
+    """
+
+    return round(float(val) / 0.3048, 10)
+
+
+def m2nm(val):
+    """
+    This function converts from meters to nautical miles
+
+    :param val: value in meters
+    :returns: value in nautical miles
+    """
+    return val / 1852.0
+
+
+def s2h(val):
+    """
+    This function converts from seconds to hours
+
+    :param val: value in seconds
+    :returns: value in hours
+    """
+    return val / 3600.0
+
+
+def ms2kt(val):
+    """
+    This function converts from m s^-1 to kt
+
+    :param val: value in m s^-1
+    :returns: value in kt
+    """
+
+    if val is None:
+        return None
+    else:
+        return round(float(val) / 0.514444, 10)
+
+
+def kg2lb(val):
+    """
+    This function converts from kg to lb
+
+    :param val: value in kg
+    :returns: value in lb
+    """
+    return val / 0.453592
+
+
+def rad2deg(val):
+    """
+    This function converts from radians to decimal degrees
+
+    :param val: value in radians
+    :returns: value in decimal degrees
+    """
+    return val / pi * 180.0
+
+
+def hp2W(val):
+    """
+    This function converts from horsepower to watts
+
+    :param val: value in horsepower
+    :returns: value in watts
+    """
+    return val * 745.699872
+
+
+def date2posix(val):
+    """
+    This function converts a date format to posix
+
+    :param val: date in %Y-%m-%d %H:%M:%S format
+    :returns: posix time referenced to 01-01-1970 [s]
+    """
+    return mktime(datetime.strptime(val, "%Y-%m-%d %H:%M:%S").timetuple())
+
+
+def unix2date(val):
+    """
+    This function converts posix to date format
+
+    :param val: time referenced to 01-01-1970 [s]
+    :returns: date in %Y-%m-%d %H:%M:%S format
+    """
+    return datetime.fromtimestamp(int(val)).strftime("%Y-%m-%d %H:%M:%S")
+
+
+convertFrom = {
+    "unix": unix2date,
+    "ft": ft2m,
+    "nm": nm2m,
+    "h": h2s,
+    "kt": kt2ms,
+    "lb": lb2kg,
+    "deg": deg2rad,
+    "date": date2posix,
+    "rad": rad2deg,
+    "ms": ms2kt,
+    "m": m2ft,
+    "kg": kg2lb,
+    "s": s2h,
+}