pygnss 0.2.0__cp312-cp312-musllinux_1_2_x86_64.whl → 0.3.0__cp312-cp312-musllinux_1_2_x86_64.whl

pygnss/__init__.py

@@ -1 +1 @@
-__version__ = "0.2.0"
+__version__ = "0.3.0"

pygnss/geodetic.py

@@ -701,7 +701,7 @@ def body_to_enu(yaw_deg, pitch_deg, roll_deg, x, y, z, matrix=None):
     BODY_TO_ENU_MATRIX IS ACTUALLY IMPLEMENTING BODY2NED_MATRIX!!!!!!
 
     >>> body_to_enu(0, 0, 0, 1, 1, 1)
-    (np.float64(1.0), np.float64(1.0), np.float64(-1.0))
+    (1.0, 1.0, -1.0)
 
     >>> enu = body_to_enu(90, 0, 0, 1, 1, 1)
     >>> np.round(enu)
@@ -809,11 +809,11 @@ def enu_to_ecef(longitude_deg, latitude_deg, e, n, u, matrix=None):
 
     >>> enu = (0.5, 0.5, 1.0)
     >>> enu_to_ecef(0, 0, *enu)
-    (np.float64(1.0), np.float64(0.5), np.float64(0.5))
+    (1.0, 0.5, 0.5)
 
     >>> enu = (0.0, 0.0, 1.0)
     >>> enu_to_ecef(0, 0, *enu)
-    (np.float64(1.0), np.float64(0.0), np.float64(0.0))
+    (1.0, 0.0, 0.0)
 
     >>> lons = [90, 180, 270]
     >>> lats = [0, 0, 0]
@@ -822,23 +822,23 @@ def enu_to_ecef(longitude_deg, latitude_deg, e, n, u, matrix=None):
     >>> us = [1, -1, -1]
     >>> dxyz = enu_to_ecef(lons, lats, es, ns, us)
     >>> round(dxyz[0][0], 8)
-    np.float64(1.0)
+    1.0
     >>> round(dxyz[0][1], 8)
-    np.float64(1.0)
+    1.0
     >>> round(dxyz[0][2], 8)
-    np.float64(1.0)
+    1.0
     >>> round(dxyz[1][0], 8)
-    np.float64(1.0)
+    1.0
     >>> round(dxyz[1][1], 8)
-    np.float64(1.0)
+    1.0
     >>> round(dxyz[1][2], 8)
-    np.float64(1.0)
+    1.0
     >>> round(dxyz[2][0], 8)
-    np.float64(1.0)
+    1.0
     >>> round(dxyz[2][1], 8)
-    np.float64(1.0)
+    1.0
     >>> round(dxyz[2][2], 8)
-    np.float64(1.0)
+    1.0
     """
 
     if matrix is None:
@@ -860,12 +860,12 @@ def xyz_to_enu(ref_pos, x, y, z, a=WGS84_A, e=WGS84_E):
     >>> dxyz = (1.0, 0.0, 0.0)      # Deviation relative to reference position
     >>> enu = xyz_to_enu(xyz, *dxyz)  # Conversion to ENU at reference position
     >>> enu
-    (np.float64(0.0), np.float64(0.0), np.float64(1.0))
+    (0.0, 0.0, 1.0)
 
     >>> enu = (0.5, 0.5, 1.0)
     >>> dxyz = enu_to_ecef(0, 0, *enu)
     >>> dxyz
-    (np.float64(1.0), np.float64(0.5), np.float64(0.5))
+    (1.0, 0.5, 0.5)
     >>> xyz_to_enu(xyz, *dxyz) == enu
     True
     """
@@ -886,16 +886,16 @@ def ecef_to_enu(longitude_deg, latitude_deg, x, y, z, matrix=None):
 
     >>> dxyz = (1, 0, 0)
     >>> ecef_to_enu(0, 0, *dxyz)
-    (np.float64(0.0), np.float64(0.0), np.float64(1.0))
+    (0.0, 0.0, 1.0)
 
     >>> dxyz = (1, 0.5, 0.5)
     >>> ecef_to_enu(0, 0, *dxyz)
-    (np.float64(0.5), np.float64(0.5), np.float64(1.0))
+    (0.5, 0.5, 1.0)
 
     >>> enu = (0.5, 0.5, 1)
     >>> dxyz = enu_to_ecef(0, 0, *enu)
     >>> dxyz
-    (np.float64(1.0), np.float64(0.5), np.float64(0.5))
+    (1.0, 0.5, 0.5)
     >>> ecef_to_enu(0, 0, *dxyz) == enu
     True
     """

pygnss/nequick.py

@@ -0,0 +1,133 @@
+import datetime
+import math
+from typing import List
+
+import numpy as np
+
+import nequick
+
+
+class GimIonexHandler(nequick.GimHandler):
+    """
+    A handler that accumulates GIMs and then generates an IONEX file
+    """
+
+    def __init__(self, coeffs: nequick.Coefficients):
+        self._coeffs = coeffs
+        self._gims: List[nequick.gim.Gim] = []
+
+    def process(self, gim: nequick.Gim):
+        """
+        Store the incoming gim for later process
+        """
+
+        # Check that the latitude and longitude values are
+        # the same as the last appended gim
+        if len(self._gims) > 0:
+            last_gim = self._gims[-1]
+            if np.array_equal(last_gim.latitudes, gim.latitudes) == False:
+                raise ValueError("Latitude values do not match")
+            if np.array_equal(last_gim.longitudes, gim.longitudes) == False:
+                raise ValueError("Longitude values do not match")
+
+        self._gims.append(gim)
+
+    def to_ionex(self, pgm: str = "pygnss", runby: str = "pygnss") -> str:
+
+        EXPONENT = -1
+
+        # Sort the IONEX files by epoch
+        self._gims.sort(key=lambda gim: gim.epoch)
+
+        first_epoch = self._gims[0].epoch
+        last_epoch = self._gims[-1].epoch
+        n_maps = len(self._gims)
+
+        lat_0 = self._gims[0].latitudes[0]
+        lat_1 = self._gims[0].latitudes[-1]
+        dlat = self._gims[0].latitudes[1] - self._gims[0].latitudes[0]
+
+        # We will print the map from North to South, therefore check if the
+        # latitudes need to be reversed
+        latitude_reversal = lat_0 < lat_1
+        if latitude_reversal:
+            lat_0 = self._gims[0].latitudes[-1]
+            lat_1 = self._gims[0].latitudes[0]
+            dlat = self._gims[0].latitudes[0] - self._gims[0].latitudes[1]
+
+        lon_0 = self._gims[0].longitudes[0]
+        lon_1 = self._gims[0].longitudes[-1]
+        dlon = self._gims[0].longitudes[1] - self._gims[0].longitudes[0]
+
+        doc = ""
+
+        # Header
+        today = datetime.datetime.now()
+        epoch_str = today.strftime('%d-%b-%y %H:%M')
+
+        doc +="     1.0            IONOSPHERE MAPS     NEQUICK             IONEX VERSION / TYPE\n"
+        doc +=f"{pgm[:20]:<20}{runby[:20]:<20}{epoch_str[:20]:<20}PGM / RUN BY / DATE\n"
+        doc +="Maps computed using the NeQuick model with the following    COMMENT\n"
+        doc +="coefficients:                                               COMMENT\n"
+        doc += f"{EXPONENT:>6}                                                      EXPONENT\n"
+        doc +=f"a0={self._coeffs.a0:<17.6f}a1={self._coeffs.a1:<17.8f}a2={self._coeffs.a2:<17.11f}COMMENT\n"
+        doc += first_epoch.strftime("  %Y    %m    %d    %H    %M    %S                        EPOCH OF FIRST MAP\n")
+        doc += last_epoch.strftime("  %Y    %m    %d    %H    %M    %S                        EPOCH OF LAST MAP\n")
+        doc += "     0                                                      INTERVAL\n"
+        doc += f"{n_maps:>6}                                                      # OF MAPS IN FILE\n"
+        doc += "  NONE                                                      MAPPING FUNCTION\n"
+        doc += "     0.0                                                    ELEVATION CUTOFF\n"
+        doc += "                                                            OBSERVABLES USED\n"
+        doc += "  6371.0                                                    BASE RADIUS\n"
+        doc += "     2                                                      MAP DIMENSION\n"
+        doc += "   450.0 450.0   0.0                                        HGT1 / HGT2 / DHGT\n"
+        doc += f"  {lat_0:6.1f}{lat_1:6.1f}{dlat:6.1f}                                        LAT1 / LAT2 / DLAT\n"
+        doc += f"  {lon_0:6.1f}{lon_1:6.1f}{dlon:6.1f}                                        LON1 / LON2 / DLON\n"
+        doc += "                                                            END OF HEADER\n"
+
+        # Body: For each GIM file, write the VTEC values
+        for i, gim in enumerate(self._gims):
+
+            doc += f"{i+1:>6}                                                      START OF TEC MAP\n"
+            doc += gim.epoch.strftime("  %Y    %m    %d    %H    %M    %S                        EPOCH OF CURRENT MAP\n")
+
+            n_latitudes = len(gim.latitudes)
+
+            for i_lat, lat in enumerate(gim.latitudes):
+
+                if latitude_reversal:
+                    i_lat = n_latitudes - 1 - i_lat
+
+                lat = gim.latitudes[i_lat]
+                doc += f"  {lat:6.1f}{lon_0:6.1f}{lon_1:6.1f}{dlon:6.1f} 450.0                            LAT/LON1/LON2/DLON/H"
+
+                lat_row = gim.vtec_values[i_lat]
+                for i_lon, _ in enumerate(gim.longitudes):
+
+                    if i_lon % 16 == 0:
+                        doc += "\n"
+
+                    vtec = lat_row[i_lon] / math.pow(10, EXPONENT)
+                    doc += f"{int(vtec):>5d}"
+
+                doc += "\n"
+
+            doc += f"{i+1:>6}                                                      END OF TEC MAP\n"
+
+        # Tail
+        doc += "                                                            END OF FILE\n"
+
+        return doc
+
+def to_ionex(coeffs: nequick.Coefficients, dates: List[datetime.datetime]) -> str:
+
+    doc = None
+
+    gim_handler = GimIonexHandler(coeffs)
+
+    for date in dates:
+        nequick.to_gim(coeffs, date, gim_handler=gim_handler)
+
+    doc = gim_handler.to_ionex()
+
+    return doc

pygnss/stats.py

@@ -70,6 +70,6 @@ def rms(values: Iterable) -> float:
 
     >>> array = [1, 2, 3, 4, 5]
     >>> rms(array)
-    np.float64(3.3166247903554)
+    3.3166247903554
     """
     return np.sqrt(np.mean(np.square(values)))

{pygnss-0.2.0.dist-info → pygnss-0.3.0.dist-info}/METADATA

@@ -1,12 +1,13 @@
-Metadata-Version: 2.2
+Metadata-Version: 2.4
 Name: pygnss
-Version: 0.2.0
+Version: 0.3.0
 Summary: Package with utilities and tools for GNSS data processing
 Author-email: Miquel Garcia-Fernandez <miquel@mgfernan.com>
 License: MIT
 Requires-Python: >=3.10
 Description-Content-Type: text/markdown
 License-File: LICENSE
+Requires-Dist: nequick
 Requires-Dist: numpy
 Requires-Dist: pandas>=2.2
 Requires-Dist: pyarrow>=18.0.0
@@ -17,6 +18,7 @@ Requires-Dist: pytest-mocha>=0.4.0; extra == "test"
 Requires-Dist: flake8>=7.0.0; extra == "test"
 Provides-Extra: release
 Requires-Dist: python-semantic-release>=9.4.0; extra == "release"
+Dynamic: license-file
 
 # GNSS and Navigation modules
 

{pygnss-0.2.0.dist-info → pygnss-0.3.0.dist-info}/RECORD

@@ -1,36 +1,37 @@
-pygnss-0.2.0.dist-info/LICENSE,sha256=Wwany6RAAZ9vVHjFLA9KBJ0HE77d52s2NOUA1CPAEug,1067
-pygnss-0.2.0.dist-info/WHEEL,sha256=c9GU1jhDUn7S-APtjPhve05mT2FAOW_4NMGdpCIyIBM,112
-pygnss-0.2.0.dist-info/METADATA,sha256=Hak3b6_DXaFefbHOylzU0BfjjErmA7MxzZhbTqzSdks,1614
-pygnss-0.2.0.dist-info/top_level.txt,sha256=oZRSR-qOv98VW2PRRMGCVNCJmewcJjyJYmxzxfeimtg,7
-pygnss-0.2.0.dist-info/RECORD,,
-pygnss-0.2.0.dist-info/entry_points.txt,sha256=mCuKrljB_wh9ZQVROiId9m68EDbTiY1oef_L1N3IDDA,262
-pygnss/sinex.py,sha256=nErOmGCFFmGSnmWGNTJhaj3yZ6IIB8GgtW5WPypJc6U,3057
-pygnss/file.py,sha256=kkMBWjoTPkxJD1UgH0mXJT2fxnhU8u7_l2Ph5Xz2-hY,933
-pygnss/__init__.py,sha256=Zn1KFblwuFHiDRdRAiRnDBRkbPttWh44jKa5zG2ov0E,22
-pygnss/geodetic.py,sha256=gfVsOeEKLn2RaJYpaCk0OrQpYz6QiDPMX6PoJHEaP9Q,34029
+pygnss/__init__.py,sha256=VrXpHDu3erkzwl_WXrqINBm9xWkcyUy53IQOj042dOs,22
+pygnss/_c_ext.cpython-312-x86_64-linux-musl.so,sha256=ctqHTTBBdSAmNdjEjiY4tm7dgr5yDsZfDqepX3tb1Lc,83816
 pygnss/cl.py,sha256=ISmd2RjikUMmj3nLPN0VSjvQLG5rLizp2X2ajeBkoDE,4509
 pygnss/constants.py,sha256=1hF6K92X6E6Ofo0rAuCBCgrwln9jxio26RV2a6vyURk,133
 pygnss/decorator.py,sha256=ldlZuvwuIlJf2pkoWteyXyp5tLds8KRkphrPsrURw9U,491
-pygnss/tensorial.py,sha256=aA0-0WK2MXhDUg0_8HMbECOt9cXmp3EnKFQXjdYMBXA,1598
-pygnss/time.py,sha256=YdMNs2xA43LrSgEOgB7jpEq0dCWv89fUBF5syDLjbu0,11178
+pygnss/file.py,sha256=kkMBWjoTPkxJD1UgH0mXJT2fxnhU8u7_l2Ph5Xz2-hY,933
+pygnss/geodetic.py,sha256=3q8Rpl4b5CxGlhdn1nQRBHHSW1v-0PBFz54zOeVyO74,33633
 pygnss/hatanaka.py,sha256=P9XG6bZwUzfAPYn--6-DXfFQIEefeimE7fMJm_DF5zE,1951
-pygnss/_c_ext.cpython-312-x86_64-linux-musl.so,sha256=HIZmPjMg3D5ZISMjeWWIsqK-MCepnPHQLwPJnAYO82Q,83184
-pygnss/stats.py,sha256=mDiY0K-VTndlFEkbxTzq9PYxCOjYDYsY3ZQV0PuMREM,1924
-pygnss/rinex.py,sha256=LsOOh3Fc263kkM8KOUBNeMeIAmbOn2ASSBO4rAUJWj8,68783
 pygnss/logger.py,sha256=4kvcTWXPoiG-MlyP6B330l4Fu7MfCuDjuIlIiLA8f1Y,1479
-pygnss/gnss/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-pygnss/gnss/types.py,sha256=lmL15KRckRiTwVkYvGzF4c1BrojnQlegrYCXSz1hGaI,10377
-pygnss/gnss/observables.py,sha256=0x0NLkTjxf8cO9F_f_Q1b-1hEeoNjWB2x-53ecUEv0M,1656
-pygnss/gnss/edit.py,sha256=T1r0WbJmt8tLJpG_IIsy4Atej6cy0IStBaSGxw0S5ho,1884
-pygnss/gnss/residuals.py,sha256=8qKGNOYkrqxHGOSjIfH21K82PAqEh2068kf78j5usL8,1244
-pygnss/filter/ukf.py,sha256=koaal4aLyYtYwPhVIvjn4uXsmSDDOQnnFEv4xenryqo,10646
+pygnss/nequick.py,sha256=ewtccKH1tdh8z6iuM6yEvDCakOKbzgkMnEWh3jWg6R4,5490
+pygnss/rinex.py,sha256=LsOOh3Fc263kkM8KOUBNeMeIAmbOn2ASSBO4rAUJWj8,68783
+pygnss/sinex.py,sha256=nErOmGCFFmGSnmWGNTJhaj3yZ6IIB8GgtW5WPypJc6U,3057
+pygnss/stats.py,sha256=GYZfcyDvbM9xamWIyVlqyN5-DPJzTLJrybRrcNV6Z6o,1912
+pygnss/tensorial.py,sha256=aA0-0WK2MXhDUg0_8HMbECOt9cXmp3EnKFQXjdYMBXA,1598
+pygnss/time.py,sha256=YdMNs2xA43LrSgEOgB7jpEq0dCWv89fUBF5syDLjbu0,11178
+pygnss/_c_ext/src/hatanaka.c,sha256=YNWaMzQQQnTNls5J6TMNuyhlq505NGDfzU-MJAHab8Q,2520
+pygnss/_c_ext/src/helpers.c,sha256=gINr73ktRgox_S7fYdFR58lLqAUACRpJfog4M5BW1-Q,364
+pygnss/_c_ext/src/mtable_init.c,sha256=5w869E6PX-ca9UHhKBxLFRW694-VaNwGlMs0I5v99mk,1132
 pygnss/filter/__init__.py,sha256=0gGBaZZgaFSiTxt0Mycn0oGZjJyATo8AAAgXhxh9k6c,1850
 pygnss/filter/ekf.py,sha256=D-qtalFRguVj4HcOtCtE0lGGXtVibC7qg-AZhKjmwgM,2017
 pygnss/filter/models.py,sha256=gXq7-YBcAoDq4-7Wr0ChNWxwXr9m1EEhUnlLtKVlsAQ,2165
+pygnss/filter/ukf.py,sha256=koaal4aLyYtYwPhVIvjn4uXsmSDDOQnnFEv4xenryqo,10646
+pygnss/gnss/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pygnss/gnss/edit.py,sha256=T1r0WbJmt8tLJpG_IIsy4Atej6cy0IStBaSGxw0S5ho,1884
+pygnss/gnss/observables.py,sha256=0x0NLkTjxf8cO9F_f_Q1b-1hEeoNjWB2x-53ecUEv0M,1656
+pygnss/gnss/residuals.py,sha256=8qKGNOYkrqxHGOSjIfH21K82PAqEh2068kf78j5usL8,1244
+pygnss/gnss/types.py,sha256=lmL15KRckRiTwVkYvGzF4c1BrojnQlegrYCXSz1hGaI,10377
 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/_c_ext/src/hatanaka.c,sha256=YNWaMzQQQnTNls5J6TMNuyhlq505NGDfzU-MJAHab8Q,2520
-pygnss/_c_ext/src/helpers.c,sha256=gINr73ktRgox_S7fYdFR58lLqAUACRpJfog4M5BW1-Q,364
-pygnss/_c_ext/src/mtable_init.c,sha256=5w869E6PX-ca9UHhKBxLFRW694-VaNwGlMs0I5v99mk,1132
 pygnss/parsers/rtklib/stats.py,sha256=YV6yadxMeQMQYZvsUCaSf4ZTpK8Bbv3f2xgu0l4PekA,5449
+pygnss-0.3.0.dist-info/METADATA,sha256=kDz5U4Y1CzGww42HgE5Fkj05KGNXZOs1SkYhuiHmPL8,1659
+pygnss-0.3.0.dist-info/WHEEL,sha256=ENczIf-GYPIfQ0_PXY8w5mFbOCUfjnGClbW6hbN_vLc,112
+pygnss-0.3.0.dist-info/entry_points.txt,sha256=mCuKrljB_wh9ZQVROiId9m68EDbTiY1oef_L1N3IDDA,262
+pygnss-0.3.0.dist-info/top_level.txt,sha256=oZRSR-qOv98VW2PRRMGCVNCJmewcJjyJYmxzxfeimtg,7
+pygnss-0.3.0.dist-info/RECORD,,
+pygnss-0.3.0.dist-info/licenses/LICENSE,sha256=Wwany6RAAZ9vVHjFLA9KBJ0HE77d52s2NOUA1CPAEug,1067

{pygnss-0.2.0.dist-info → pygnss-0.3.0.dist-info}/WHEEL

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