cssrlib 1.2.0__tar.gz → 1.3.2__tar.gz

{cssrlib-1.2.0/src/cssrlib.egg-info → cssrlib-1.3.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: cssrlib
-Version: 1.2.0
+Version: 1.3.2
 Summary: Python Toolkit for PPP/PPP-RTK positioning
 Home-page: https://github.com/hirokawa/cssrlib
 Author: Rui Hirokawa
@@ -16,7 +16,7 @@ Requires-Dist: notebook
 Requires-Dist: numpy
 Requires-Dist: matplotlib
 Requires-Dist: cartopy
-Requires-Dist: pysolid
+Requires-Dist: cryptography
 Dynamic: license-file
 
 CSSRlib - Toolkit for PPP-RTK/RTK in Python using Compact SSR

{cssrlib-1.2.0 → cssrlib-1.3.2}/pyproject.toml

@@ -7,7 +7,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "cssrlib"
-version = "1.2.0"
+version = "1.3.2"
 description = "Python Toolkit for PPP/PPP-RTK positioning"
 authors = [
     {name="Rui Hirokawa", email="rui.hirokawa@gmail.com"},
@@ -23,5 +23,5 @@ dependencies = [
     "numpy",
     "matplotlib",
     "cartopy",
-    "pysolid",
+    "cryptography",
 ]

{cssrlib-1.2.0 → cssrlib-1.3.2}/setup.cfg

@@ -1,6 +1,6 @@
 [metadata]
 name = cssrlib
-version = 1.2.0
+version = 1.3.2
 author = Rui Hirokawa
 author_email = rui.hirokawa@gmail.com
 description = Python Toolkit for PPP/PPP-RTK positioning

{cssrlib-1.2.0 → cssrlib-1.3.2}/src/cssrlib/__init__.py

@@ -1,3 +1,5 @@
+"""CSSRLIB package for GNSS SSR and PPP/PPP-RTK processing."""
+
 
 from cssrlib.gnss import (gpst2time,epoch2time,ecef2pos,
                           prn2sat)

{cssrlib-1.2.0 → cssrlib-1.3.2}/src/cssrlib/cssr_bds.py

@@ -1,7 +1,6 @@
-"""
-BDS PPP correction data decoder
+"""BDS PPP correction data decoder.
 
-[1] BeiDou Navigation Satellite System Signal In Space Interface Control Document
+[1] BeiDou Navigation Satellite System Signal In Space ICD
 Precise Point Positioning Service Signal PPP-B2b (Version 1.0), 2020
 """
 
@@ -12,7 +11,10 @@ from cssrlib.gnss import bdt2time, bdt2gpst, uGNSS, uSIG, uTYP, rSigRnx
 
 
 class cssr_bds(cssr):
+    """Class to decode BDS PPP (B2b) Compact SSR messages."""
+
     def __init__(self, foutname=None):
+        """Initialize BDS PPP SSR decoder."""
         super().__init__(foutname)
         self.MAXNET = 1
         self.cssrmode = sCSSRTYPE.BDS_PPP
@@ -23,6 +25,7 @@ class cssr_bds(cssr):
         self.cb_scl = 0.017
 
     def ssig2rsig(self, sys: sGNSS, utyp: uTYP, ssig):
+        """Convert SSR signal index to RINEX signal code."""
         gps_tbl = {
             0: uSIG.L1C,
             1: uSIG.L1P,
@@ -72,13 +75,14 @@ class cssr_bds(cssr):
         return rSigRnx(sys, utyp, usig_tbl[ssig])
 
     def sval(self, u, n, scl):
-        """ calculate signed value based on n-bit int, lsb """
+        """Calculate signed value based on n-bit int, lsb."""
         invalid = -2**(n-1)
         dnu = -(2**(n-1)-1)  # this value seems to be invalid
         y = np.nan if u == invalid or u == dnu else u*scl
         return y
 
     def slot2prn(self, slot):
+        """Convert BDS B2b satellite slot index to system and PRN."""
         prn = 0
         sys = uGNSS.NONE
         if slot >= 1 and slot <= 63:
@@ -96,6 +100,7 @@ class cssr_bds(cssr):
         return sys, prn
 
     def decode_head(self, msg, i, st=-1):
+        """Decode BDS PPP message header."""
         self.tod, _, iodssr = bs.unpack_from('u17u4u2', msg, i)
         i += 23
 
@@ -107,10 +112,11 @@ class cssr_bds(cssr):
             if self.week >= 0:
                 self.time = bdt2gpst(bdt2time(self.week, self.tow))
 
-        head = {'uint': 0, 'mi': 0, 'iodssr': iodssr}
+        head = {'udi': 0, 'mi': 0, 'iodssr': iodssr}
         return head, i
 
     def add_gnss(self, mask, blen, gnss):
+        """Add decoded GNSS satellite mask entries."""
         prn, nsat = self.decode_mask(mask, blen)
         self.nsat_g[gnss] = nsat
         self.nsat_n += nsat
@@ -124,7 +130,7 @@ class cssr_bds(cssr):
             self.gnss_n.append(gnss)
 
     def decode_cssr_mask(self, msg, i):
-        """decode MT1 Mask message """
+        """Decode MT1 Mask message."""
         head, i = self.decode_head(msg, i, sCSSR.MASK)
 
         self.iodp = bs.unpack_from('u4', msg, i)[0]
@@ -179,6 +185,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_orb_sat(self, msg, i, inet, sat_n):
+        """Decode satellite orbit correction for BDS PPP."""
         slot, iodn, iodc = bs.unpack_from('u9u10u3', msg, i)
         i += 22
 
@@ -216,7 +223,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_orb(self, msg, i, inet=0):
-        """decode MT2 orbit + URA message """
+        """Decode MT2 orbit + URA message."""
         head, i = self.decode_head(msg, i)
         sat_n = np.array(self.sat_n)
 
@@ -233,7 +240,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_cbias(self, msg, i, inet=0):
-        """decode MT3 Code Bias Correction message """
+        """Decode MT3 Code Bias Correction message."""
         head, i = self.decode_head(msg, i)
         nsat = self.nsat_n
         self.flg_net = False
@@ -270,7 +277,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_clk_sat(self, msg, i, inet, sat):
-        """ decode clock correction for satellite """
+        """Decode clock correction for satellite."""
         iodc, dclk = bs.unpack_from('u3s{:d}'.format(self.dclk_blen), msg, i)
         i += 3+self.dclk_blen
 
@@ -285,7 +292,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_clk(self, msg, i, inet=0):
-        """decode MT4 Clock Correction message """
+        """Decode MT4 Clock Correction message."""
         head, i = self.decode_head(msg, i)
         if self.iodssr != head['iodssr']:
             return -1
@@ -311,7 +318,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_ura(self, msg, i):
-        """decode MT5 URA message """
+        """Decode MT5 URA message."""
         head, i = self.decode_head(msg, i)
         if self.iodssr != head['iodssr']:
             return -1
@@ -332,7 +339,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_comb1(self, msg, i, inet=0):
-        """decode MT6 combined message #1 """
+        """Decode MT6 combined message #1."""
         numc, numo = bs.unpack_from('u5u3', msg, i)
         sat_n = np.array(self.sat_n)
 
@@ -354,7 +361,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr_comb2(self, msg, i, inet=0):
-        """decode MT7 combined message #2 """
+        """Decode MT7 combined message #2."""
         numc, numo = bs.unpack_from('u5u3', msg, i)
         sat_n = np.array(self.sat_n)
 
@@ -377,6 +384,7 @@ class cssr_bds(cssr):
         return i
 
     def decode_cssr(self, msg, i=0):
+        """Decode BDS PPP B2b message."""
         mt = bs.unpack_from('u6', msg, i)[0]
         i += 6
 

{cssrlib-1.2.0 → cssrlib-1.3.2}/src/cssrlib/cssr_has.py

@@ -1,5 +1,4 @@
-"""
-Galileo HAS correction data decoder
+"""Galileo HAS correction data decoder.
 
 [1] Galileo High Accuracy Service Signal-in-Space
   Interface Control Document (HAS SIS ICD), Issue 1.0, May 2022
@@ -15,7 +14,10 @@ from binascii import unhexlify
 
 
 class cssr_has(cssr):
+    """Class to decode Galileo HAS Compact SSR messages."""
+
     def __init__(self, foutname=None):
+        """Initialize Galileo HAS SSR decoder."""
         super().__init__(foutname)
         self.MAXNET = 1
         self.cssrmode = sCSSRTYPE.GAL_HAS_SIS
@@ -29,13 +31,14 @@ class cssr_has(cssr):
         self.pb_scl = 0.01  # cycles
 
     def sval(self, u, n, scl):
-        """ calculate signed value based on n-bit int, lsb """
+        """Calculate signed value based on n-bit int, lsb."""
         invalid = -2**(n-1)
         dnu = 2**(n-1)-1
         y = np.nan if u == invalid or u == dnu else u*scl
         return y
 
     def decode_cssr_clk_sub(self, msg, i=0):
+        """Decode Galileo HAS clock subset correction message."""
         head, i = self.decode_head(msg, i)
         self.flg_net = False
         if self.iodssr != head['iodssr']:
@@ -58,12 +61,11 @@ class cssr_has(cssr):
         return i
 
     def decode_head(self, msg, i, st=-1):
-        """ decode header part of HAS messages """
-
+        """Decode header part of HAS messages."""
         if st == sCSSR.MASK:
-            ui = 0
+            udi = 0
         else:
-            ui = bs.unpack_from('u4', msg, i)[0]
+            udi = bs.unpack_from('u4', msg, i)[0]
             i += 4
 
         if self.tow0 >= 0:
@@ -71,13 +73,12 @@ class cssr_has(cssr):
             if self.week >= 0:
                 self.time = gpst2time(self.week, self.tow)
 
-        head = {'uint': ui, 'mi': 0, 'iodssr': self.iodssr}
+        head = {'udi': udi, 'mi': 0, 'iodssr': self.iodssr}
         return head, i
 
     def decode_cssr(self, msg, i=0):
-        """ decode HAS messages """
-
-        # Galileo HAS-SIS only defines MT=1, MT=4073 is CSSR message used in the
+        """Decode HAS messages."""
+        # Galileo HAS-SIS only defines MT=1; MT=4073 is CSSR used in the
         # Galileo HAS test dataset
         #
         if self.msgtype != 1 and self.msgtype != 4073:
@@ -131,9 +132,10 @@ class cssr_has(cssr):
 
 
 class cnav_msg():
-    """ class to handle Galileo E6 (CNAV) message """
+    """Class to handle Galileo E6 (CNAV) message."""
 
     def __init__(self):
+        """Initialize Galileo C/NAV HAS message decoder."""
         self.mid_decoded = []
         self.mid_ = -1
         self.ms_ = -1
@@ -148,10 +150,11 @@ class cnav_msg():
         self.monlevel = 0
 
     def load_gmat(self, file_gm):
+        """Load Reed-Solomon generator matrix for Galileo HAS decoding."""
         self.gMat = np.genfromtxt(file_gm, dtype="u1", delimiter=",")
 
     def decode_has_header(self, buff, i):
-        """ decode header of HAS pages (obsoleted) """
+        """Decode header of HAS pages (obsoleted)."""
         if bs.unpack_from('u24', buff, i)[0] == 0xaf3bc3:
             return 0, 0, 0, 0, 0
 
@@ -160,7 +163,7 @@ class cnav_msg():
         return hass, mt, mid, ms, pid
 
     def decode_has_page(self, idx, has_pages, gMat, ms):
-        """ HPVRS decoding for RS(255,32,224) """
+        """Decode HPVRS Reed-Solomon RS(255,32,224) page."""
         HASmsg = bytes()
         k = len(idx)
         if k >= ms:
@@ -172,8 +175,7 @@ class cnav_msg():
         return HASmsg
 
     def decode_cnav(self, tow, vi):
-        """ decode Galileo CNAV message """
-
+        """Decode Galileo CNAV message."""
         HASmsg = None
 
         for vn in vi:
@@ -183,8 +185,13 @@ class cnav_msg():
                 continue
             hass, res = bs.unpack_from('u2u2', buff, i)
             i += 4
+            # HAS status
             if hass >= 2:  # 0:test,1:operational,2:res,3:dnu
                 continue
+            # mt: type of message
+            # mid: id of the message
+            # ms: size of non-encoded message in number of pages 0:1,...
+            # pid: id of the transmitted HAS encoded page
             mt, mid, ms, pid = bs.unpack_from('u2u5u5u8', buff, i)
 
             self.msgtype = mt
@@ -202,7 +209,7 @@ class cnav_msg():
             if self.monlevel >= 2:
                 print(f"{mt} {mid} {ms} {pid}")
 
-        if len(self.rec) >= self.ms_:
+        if self.ms_ > 0 and len(self.rec) >= self.ms_:
             if self.monlevel >= 2:
                 print(" data collected mid={:2d} ms={:2d} tow={:.0f}"
                       .format(self.mid_, self.ms_, tow))

{cssrlib-1.2.0 → cssrlib-1.3.2}/src/cssrlib/cssr_mdc.py

@@ -1,20 +1,23 @@
-"""
-QZSS MADOCA-PPP correction data decoder
+"""QZSS MADOCA-PPP correction data decoder.
 
 [1] Quasi-Zenith Satellite System Interface Specification Multi-GNSS
     Advanced Orbit and Clock Augmentation - Precise Point Positioning
-    (IS-QZSS-MDC-004), 2025
+    (IS-QZSS-MDC-005draft), October, 2026
 
 """
 
 import numpy as np
 import bitstruct as bs
 from cssrlib.cssrlib import cssr, sCSSRTYPE, sCSSR, local_corr, sCType
-from cssrlib.gnss import gpst2time, uGNSS, prn2sat, rCST
+from cssrlib.gnss import gpst2time, time2str, uGNSS, prn2sat, rCST, sat2id
+from cssrlib.gnss import pos2ecef, gtime_t
 
 
 class areaInfo():
+    """Class to store MADOCA-PPP service area information."""
+
     def __init__(self, sid, latr, lonr, p1=0, p2=0):
+        """Initialize MADOCA-PPP area information."""
         self.sid = sid
         self.latr = latr
         self.lonr = lonr
@@ -26,7 +29,10 @@ class areaInfo():
 
 
 class ionoCorr():
+    """Class to store MADOCA-PPP regional ionospheric corrections."""
+
     def __init__(self):
+        """Initialize MADOCA-PPP ionospheric correction data."""
         self.t0 = None
         self.qi = []
         self.iodssr = 0
@@ -36,7 +42,10 @@ class ionoCorr():
 
 
 class cssr_mdc(cssr):
+    """Class to decode QZSS MADOCA-PPP L6D/L6E messages."""
+
     def __init__(self, foutname=None):
+        """Initialize MADOCA-PPP L6 message decoder."""
         super().__init__(foutname)
 
         self.cssrmode = sCSSRTYPE.QZS_MADOCA
@@ -44,11 +53,19 @@ class cssr_mdc(cssr):
 
         self.pnt = {}
         self.ci = {}
+        self.reg_alrt = {}
+        self.reg_iod = {}
 
         self.area = -1
+        self.reg = -1
+        self.alrt = 0
         self.narea_t = {1: 8, 2: 16, 3: 5, 4: 1, 5: 8}
         self.MAXNET = np.sum(list(self.narea_t.values()))
 
+        self.tmax = {sCType.CLOCK: 60.0, sCType.ORBIT: 60.0,
+                     sCType.CBIAS: 60.0, sCType.PBIAS: 60.0,
+                     sCType.TROP: 120.0, sCType.STEC: 300.0}
+
         for inet in range(self.MAXNET+1):
             self.lc.append(local_corr())
             self.lc[inet].inet = inet
@@ -58,35 +75,40 @@ class cssr_mdc(cssr):
             self.lc[inet].t0 = {}
 
     def find_grid_index(self, pos):
-        """ find region/area   """
-
+        """Find region/area."""
         latd, lond = np.rad2deg(pos[0:2])
+        rr = pos2ecef(pos)
         self.reg = self.area = -1
+        min_dist = 1e5
 
-        for reg in self.pnt.key():
-            for area in self.pnt[reg].key():
+        for reg in self.pnt.keys():
+            if self.reg_alrt.get(reg, 0) != 0:
+                continue
+            for area in self.pnt[reg].keys():
                 p = self.pnt[reg][area]
+                ref_ecef = pos2ecef(
+                    np.array([p.latr*rCST.D2R, p.lonr*rCST.D2R, 0.0]))
+                dist = np.linalg.norm(rr - ref_ecef)*1e-3  # [km]
+                if dist > min_dist:
+                    continue
                 if p.sid == 0:  # rectangle shape:
                     if p.latr-p.lats <= latd and latd <= p.latr+p.lats and \
                             p.lonr-p.lons <= lond and lond <= p.lonr+p.lons:
                         self.reg = reg
                         self.area = area
-                        break
+                        min_dist = dist
                 elif p.sid == 1:  # circle
-                    dn = (latd - p.latr)*rCST.D2R*rCST.RE_WGS84
-                    de = (lond - p.lonr)*rCST.D2R * \
-                        rCST.RE_WGS84*np.cos(p.lonr*rCST.D2R)
-                    r = np.sqrt(dn**2+de**2)
-                    if r <= p.rng:
+                    if dist <= p.rng:
                         self.reg = reg
                         self.area = area
-                        break
+                        min_dist = dist
 
         inet = self.get_inet(self.reg, self.area)
+        self.inet = inet
         return inet
 
     def get_dpos(self, pos):
-        """ calculate position offset from reference """
+        """Calculate position offset from reference."""
         if self.reg < 0 or self.area < 0:
             print("get_dpos: region, area not defined.")
             return 0, 0
@@ -99,28 +121,37 @@ class cssr_mdc(cssr):
         return dlat, dlon
 
     def get_stec(self, dlat=0.0, dlon=0.0):
-        """ calculate STEC correction by interporation """
+        """Calculate STEC correction by interporation."""
         if self.inet < 0:
             print("get_stec: region, area not defined.")
-            return 0.0
+            return np.zeros(self.nsat_n)
 
         p = self.lc[self.inet]
         # if p.inet_ref != self.iodssr:
         #    return 0.0
-        stec = np.zeros(self.nsat_n)
+        stec = np.zeros(p.nsat_n)
+        v = np.array([1, dlat, dlon, dlat*dlon, dlat**2, dlon**2])
 
-        for i, sat in enumerate(self.sat_n):
-            stec[i] = [1, dlat, dlon, dlat*dlon, dlat**2, dlon**2]@p.ci[sat]
+        for i, sat in enumerate(p.sat_n):
+            stec[i] = v@p.ci[sat]
 
         return stec
 
     def decode_mdc_stec_area(self, buff, i=0):
-        """ decoder for MT1 - STEC Coverage Message """
+        """Decode MT1 STEC coverage message."""
         tow, uid, mi, iodssr = bs.unpack_from('u20u4u1u4', buff, i)
         i += 29
+        self.tow = tow
         self.tow0 = tow//3600*3600
+        if self.week >= 0:
+            self.time = gpst2time(self.week, self.tow)
         reg, alrt, len_, narea = bs.unpack_from('u8u1u16u5', buff, i)
         i += 30
+        self.reg = reg
+        self.alrt = alrt
+        self.reg_alrt[reg] = alrt
+        self.reg_iod[reg] = iodssr
+
         if reg not in self.pnt:
             self.pnt[reg] = {}
 
@@ -130,14 +161,14 @@ class cssr_mdc(cssr):
 
             if sid == 0:  # rectangle shape
                 latr, lonr, lats, lons = bs.unpack_from('s11u12u8u8', buff, i)
-                if self.monlevel > 2:
+                if self.monlevel >= 2:
                     print(f"{reg} {area:2d} {sid} {latr*0.1:5.1f} "
                           f"{lonr*0.1:5.1f} {lats*0.1:3.1f} {lons*0.1:3.1f}")
                 self.pnt[reg][area] = areaInfo(
                     sid, latr*0.1, lonr*0.1, lats*0.1, lons*0.1)
             else:  # circle range
                 latr, lonr, rng = bs.unpack_from('s15u16u8', buff, i)
-                if self.monlevel > 2:
+                if self.monlevel >= 2:
                     print(f"{reg} {area:2d} {sid} {latr*0.01:6.2f} "
                           f"{lonr*0.01:6.2f} {rng*10}")
                 self.pnt[reg][area] = areaInfo(
@@ -146,7 +177,9 @@ class cssr_mdc(cssr):
         return i
 
     def get_inet(self, reg, area):
-        """ region, area to inet conversion """
+        """Convert region and area ID to network index."""
+        if reg < 0 or area < 0:
+            return -1
 
         inet = 0
         for r in self.narea_t.keys():
@@ -159,12 +192,23 @@ class cssr_mdc(cssr):
         return inet
 
     def decode_mdc_stec_corr(self, buff, i=0):
-        """ decoder for MT2 - STEC Correction Message """
+        """Decode MT2 STEC correction message."""
         dtow, uid, mi, iodssr = bs.unpack_from('u12u4u1u4', buff, i)
         i += 21
+        if self.tow0 >= 0:
+            tow = self.tow0 + dtow
+            if self.tow >= 0 and tow < self.tow - 1800:
+                tow += 3600
+            self.tow = tow
+        t0 = gpst2time(self.week, self.tow) if self.week >= 0 else gtime_t()
+        if self.week >= 0:
+            self.time = t0
         reg, area, stype_ = bs.unpack_from('u8u5u2', buff, i)
         i += 15
 
+        self.reg = reg
+        self.area = area
+
         nsat = bs.unpack_from('u5u5u5u5u5', buff, i)
         i += 25
         # gps, glo, gal, bds, qzss
@@ -172,7 +216,6 @@ class cssr_mdc(cssr):
 
         inet = self.get_inet(reg, area)
 
-        self.lc[inet].nsat_n = np.sum(nsat)
         ci = {}
         qi = {}
         stype = {}
@@ -181,39 +224,29 @@ class cssr_mdc(cssr):
         for gnss in range(5):
             sys = sys_t[gnss]
             for k in range(nsat[gnss]):
-                c01 = c10 = c11 = c02 = c20 = 0.0
-
-                prn, qi_, c00 = bs.unpack_from('u6u6s14', buff, i)
-                i += 26
+                prn, cls_, val = bs.unpack_from('u6u3u3', buff, i)
+                i += 12
 
                 if sys == uGNSS.QZS:
                     prn += 192
                 sat = prn2sat(sys, prn)
-                qi[sat] = qi_
+                ci_sat, i = self.decode_cssr_stec_coeff(buff, stype_, i)
+                if np.isnan(ci_sat).any():
+                    continue
+                qi[sat] = self.quality_idx(cls_, val)
                 stype[sat] = stype_
-
-                if stype_ > 0:
-                    c01, c10 = bs.unpack_from('s12s12', buff, i)
-                    i += 24
-                if stype_ > 1:
-                    c11 = bs.unpack_from('s10', buff, i)[0]
-                    i += 10
-                if stype_ > 2:
-                    c02, c20 = bs.unpack_from('s8s8', buff, i)
-                    i += 16
-
                 sat_ += [sat]
-                ci[sat] = np.array(
-                    [c00*0.05, c01*0.02, c10*0.02, c11*0.02,
-                     c02*5e-3, c20*5e-3])
+                ci[sat] = ci_sat
+                self.set_t0(inet, sat, sCType.STEC, t0)
                 j += 1
 
+        self.lc[inet].nsat_n = len(sat_)
+        self.lc[inet].flg_stec = 2
         self.lc[inet].stype = stype
         self.lc[inet].sat_n = sat_
         self.lc[inet].ci = ci
         self.lc[inet].stec_quality = qi
 
-        t0 = gpst2time(self.week, self.tow0+dtow)
         self.set_t0(inet, 0, sCType.STEC, t0)
         self.lc[inet].cstat |= (1 << sCType.STEC)
 
@@ -221,15 +254,60 @@ class cssr_mdc(cssr):
 
         return i
 
+    def out_log(self):
+        # if self.msgtype not in (1, 2):
+        #    return super(cssr_mdc, self).out_log()
+        """Write decoded MADOCA-PPP STEC messages to log file."""
+        sz_t = [1, 3, 4, 6]
+
+        if self.time == -1:
+            return
+
+        self.fh.write("{:4d}\t{:s}\n".format(self.msgtype,
+                                             time2str(self.time)))
+
+        if self.msgtype == 1:
+
+            self.fh.write(f"Reg\tArea\tsid\t"
+                          f"latr\tlonr\tlats\tlons\n")
+
+            for area in self.pnt[self.reg].keys():
+                p = self.pnt[self.reg][area]
+                self.fh.write(f"{self.reg}\t{area:2d}\t{p.sid}\t"
+                              f"{p.latr:3.1f}\t{p.lonr:4.1f}\t")
+
+                if p.sid == 0:
+                    self.fh.write(f"{p.lats:3.1f}\t{p.lons:3.1f}\n")
+                else:
+                    self.fh.write(f"{p.rng:3.1f}\n")
+
+        elif self.msgtype == 2:
+            inet = self.get_inet(self.reg, self.area)
+
+            self.fh.write(f"Reg:{self.reg}\tArea:{self.area:2d}\n")
+            self.fh.write("Sat,stype,c00,c01,c10,c11,c20,c02\n")
+
+            for sat in self.lc[inet].sat_n:
+                stype = self.lc[inet].stype[sat]
+                ci = self.lc[inet].ci[sat]
+                self.fh.write(f"{sat2id(sat)}\t{stype}")
+                for k in range(sz_t[stype]):
+                    self.fh.write(f"\t{ci[k]:6.2f}")
+                self.fh.write("\n")
+        self.fh.flush()
+
     def decode_cssr(self, msg, i=0):
-        """decode Compact SSR message with MADOCA-PPP extension """
+        """Decode Compact SSR message with MADOCA-PPP extension."""
         df = {'msgtype': 4073}
         while df['msgtype'] in [1, 2, 4073]:
+            if len(msg)*8 < i + 16:
+                return 0
             df = bs.unpack_from_dict('u12u4', ['msgtype', 'subtype'], msg, i)
             i += 16
             if df['msgtype'] not in [1, 2, 4073]:
                 return -1
 
+            self.msgtype = df['msgtype']
             self.subtype = df['subtype']
             if df['msgtype'] == 4073:  # Compact SSR
                 if self.subtype == sCSSR.MASK: