pygnss 0.7.0__cp311-cp311-musllinux_1_2_x86_64.whl → 1.1.0__cp311-cp311-musllinux_1_2_x86_64.whl

pygnss/__init__.py

@@ -1 +1 @@
-__version__ = "0.7.0"
+__version__ = "1.1.0"

pygnss/geodetic.py

@@ -969,15 +969,23 @@ def haversine(lon1_deg, lat1_deg, lon2_deg, lat2_deg, r=6371):
 
     >>> haversine(-86.67, 36.12, -118.40, 33.94, r=6371.8)
     2886.8068907353736
+
+    Arrays can also be used (for a vectorized computation)
+    >>> haversine([2, 2, 2], [89, 40, 0], [3, 3, 3], [89, 40, 0])
+    array([  1.94059443,  85.17980895, 111.19492664])
     """
     # convert decimal degrees to radians
-    lon1, lat1, lon2, lat2 = map(math.radians, [lon1_deg, lat1_deg, lon2_deg, lat2_deg])
+    lon1, lat1, lon2, lat2 = map(np.radians, [lon1_deg, lat1_deg, lon2_deg, lat2_deg])
 
     # haversine formula
     dlon = lon2 - lon1
     dlat = lat2 - lat1
-    a = math.sin(dlat / 2) ** 2 + math.cos(lat1) * math.cos(lat2) * math.sin(dlon / 2) ** 2
-    c = 2 * math.asin(math.sqrt(a))
+    a = np.sin(dlat / 2) ** 2 + np.cos(lat1) * np.cos(lat2) * np.sin(dlon / 2) ** 2
+    c = 2 * np.arcsin(np.sqrt(a))
+
+    # Return a float if input is a single float (to avoid test failures due to numpy float64)
+    if np.ndim(lon1_deg) == 0:
+        c = float(c)
 
     return c * r
 

pygnss/iono/chapman.py

@@ -0,0 +1,35 @@
+"""
+Module to compute the Chapman ionospheric model.
+
+This module provides functions to compute the Chapman ionospheric model,
+which is used to estimate the electron density in the ionosphere based on
+solar zenith angle and other parameters.
+"""
+
+import numpy as np
+
+HEIGHTS_KM = np.arange(60, 1000, 10)
+
+def compute_profile(hmF2_km: float, NmF2: float, H_km: float, zenith_angle_deg: float, heights_km: float | np.ndarray =HEIGHTS_KM) -> float | np.ndarray:
+    """
+    Compute the Chapman ionospheric model.
+
+    Parameters:
+    zenith_angle (float or np.ndarray): Solar zenith angle, the angle between the vertical and the direction to the Sun.
+    hmF2_km (float): Height of the F2 layer in kilometers.
+    NmF2 (float): Peak electron density at hmF2 in electrons/m^3.
+    H_km (float): Scale height in kilometers (typical values between 30 and 50 km).
+    heights_km (np.ndarray, optional): Heights at which to compute the electron density.
+        If None, defaults to a range from 0 to 100 km.
+
+    Returns:
+    float or np.ndarray: Computed electron density based on the Chapman model.
+    """
+    # Convert zenith angle from degrees to radians
+    zenith_angle_rad = np.radians(zenith_angle_deg)
+
+    # Include scale height
+    z = (heights_km - hmF2_km) / H_km
+
+    # Compute the Chapman function
+    return NmF2 * np.exp( 0.5 * (1 - z - np.exp(-z) / np.cos(zenith_angle_rad)))

{pygnss-0.7.0.dist-info → pygnss-1.1.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pygnss
-Version: 0.7.0
+Version: 1.1.0
 Summary: Package with utilities and tools for GNSS data processing
 Author-email: Miquel Garcia-Fernandez <miquel@mgfernan.com>
 License: MIT

{pygnss-0.7.0.dist-info → pygnss-1.1.0.dist-info}/RECORD

@@ -1,10 +1,10 @@
-pygnss/__init__.py,sha256=RaANGbRu5e-vehwXI1-Qe2ggPPfs1TQaZj072JdbLk4,22
-pygnss/_c_ext.cpython-311-x86_64-linux-musl.so,sha256=GKssFLbFGJ82ZKuhZ8ME6KG3XocBiWLrLoZOf5_pE5E,82984
+pygnss/__init__.py,sha256=LGVQyDsWifdACo7qztwb8RWWHds1E7uQ-ZqD8SAjyw4,22
+pygnss/_c_ext.cpython-311-x86_64-linux-musl.so,sha256=d29BZDsvWFnjanzqPaI-HUc_EKNzWc-CGkPjepsC8WY,82976
 pygnss/cl.py,sha256=ISmd2RjikUMmj3nLPN0VSjvQLG5rLizp2X2ajeBkoDE,4509
 pygnss/constants.py,sha256=1hF6K92X6E6Ofo0rAuCBCgrwln9jxio26RV2a6vyURk,133
 pygnss/decorator.py,sha256=qB-0jl2GTEHJdvmDruNll5X3RrdySssv94u9Hok5-lA,1438
 pygnss/file.py,sha256=kkMBWjoTPkxJD1UgH0mXJT2fxnhU8u7_l2Ph5Xz2-hY,933
-pygnss/geodetic.py,sha256=3q8Rpl4b5CxGlhdn1nQRBHHSW1v-0PBFz54zOeVyO74,33633
+pygnss/geodetic.py,sha256=RCXD4_RlzQR2DxNIFGS4iODwZZcnxklMsCw9mohYcBg,33948
 pygnss/hatanaka.py,sha256=P9XG6bZwUzfAPYn--6-DXfFQIEefeimE7fMJm_DF5zE,1951
 pygnss/ionex.py,sha256=WcY6PVwjzATjqyC0Htc5Jh_cv2a7yz4TZIcD723TpsY,14325
 pygnss/logger.py,sha256=4kvcTWXPoiG-MlyP6B330l4Fu7MfCuDjuIlIiLA8f1Y,1479
@@ -27,14 +27,15 @@ pygnss/gnss/observables.py,sha256=0x0NLkTjxf8cO9F_f_Q1b-1hEeoNjWB2x-53ecUEv0M,16
 pygnss/gnss/residuals.py,sha256=8qKGNOYkrqxHGOSjIfH21K82PAqEh2068kf78j5usL8,1244
 pygnss/gnss/types.py,sha256=lmL15KRckRiTwVkYvGzF4c1BrojnQlegrYCXSz1hGaI,10377
 pygnss/iono/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pygnss/iono/chapman.py,sha256=eyq2RaCkLpQJRzuHt_DSNaLjto1w_W1PdT_ys8A3Bjc,1388
 pygnss/iono/gim.py,sha256=khKzClxLf46mVO6BCDk9u9DgS_50d-sbWEyoUQu7Jng,3543
 pygnss/orbit/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 pygnss/orbit/kepler.py,sha256=QORTgg5yBtsQXxLWSzoZ1pmh-CwPiZlFdIYqhQhv1a0,1745
 pygnss/orbit/tle.py,sha256=6CIEielPgui3DXNv46XxOGlig31ROIwjH42xLGaeE5M,5905
 pygnss/parsers/rtklib/stats.py,sha256=YV6yadxMeQMQYZvsUCaSf4ZTpK8Bbv3f2xgu0l4PekA,5449
-pygnss-0.7.0.dist-info/METADATA,sha256=UGioECEGQJ7eQSa5yaZTnf3aNKvXG2dsiJVrhu0sPk0,1666
-pygnss-0.7.0.dist-info/WHEEL,sha256=3KlS4fjXmSH_tZmtgXZWnBbgytgKYhwpPQ7mdxt5Bus,112
-pygnss-0.7.0.dist-info/entry_points.txt,sha256=awWUCMfX3uMIAPPZqZfSUuw7Kd7defip8L7GtthKqgg,292
-pygnss-0.7.0.dist-info/top_level.txt,sha256=oZRSR-qOv98VW2PRRMGCVNCJmewcJjyJYmxzxfeimtg,7
-pygnss-0.7.0.dist-info/RECORD,,
-pygnss-0.7.0.dist-info/licenses/LICENSE,sha256=Wwany6RAAZ9vVHjFLA9KBJ0HE77d52s2NOUA1CPAEug,1067
+pygnss-1.1.0.dist-info/METADATA,sha256=ddxRfkulzB_nsS6xbvrsLHzqzRmjNBiFZHtBoo6g0dk,1666
+pygnss-1.1.0.dist-info/WHEEL,sha256=kA_iIvT-cxTFNl4I8QDfFHN1DAyqZDYakVXCaObxeLo,112
+pygnss-1.1.0.dist-info/entry_points.txt,sha256=awWUCMfX3uMIAPPZqZfSUuw7Kd7defip8L7GtthKqgg,292
+pygnss-1.1.0.dist-info/top_level.txt,sha256=oZRSR-qOv98VW2PRRMGCVNCJmewcJjyJYmxzxfeimtg,7
+pygnss-1.1.0.dist-info/RECORD,,
+pygnss-1.1.0.dist-info/licenses/LICENSE,sha256=Wwany6RAAZ9vVHjFLA9KBJ0HE77d52s2NOUA1CPAEug,1067

{pygnss-0.7.0.dist-info → pygnss-1.1.0.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (77.0.3)
+Generator: setuptools (80.9.0)
 Root-Is-Purelib: false
 Tag: cp311-cp311-musllinux_1_2_x86_64