def obsheader2(f: TextIO, useindicators: bool = False, meas: Sequence[str] = None) -> Dict[str, Any]: """ End users should use rinexheader() """ if isinstance(f, (str, Path)): with opener(f, header=True) as h: return obsheader2(h, useindicators, meas) f.seek(0) # %% selection if isinstance(meas, str): meas = [meas] if not meas or not meas[0].strip(): meas = None hdr = rinexinfo(f) Nobs = 0 # not None due to type checking for ln in f: if "END OF HEADER" in ln: break h = ln[60:80].strip() c = ln[:60] # %% measurement types if '# / TYPES OF OBSERV' in h: if Nobs == 0: Nobs = int(c[:6]) hdr[h] = c[6:].split() else: hdr[h] += c[6:].split() elif h not in hdr: # Header label hdr[h] = c # string with info else: # concatenate hdr[h] += " " + c # %% useful values try: hdr['systems'] = hdr['RINEX VERSION / TYPE'][40] except KeyError: pass hdr['Nobs'] = Nobs # 5 observations per line (incorporating LLI, SSI) hdr['Nl_sv'] = ceil(hdr['Nobs'] / 5) # %% list with receiver location in x,y,z cartesian ECEF (OPTIONAL) try: hdr['position'] = [float(j) for j in hdr['APPROX POSITION XYZ'].split()] if ecef2geodetic is not None: hdr['position_geodetic'] = ecef2geodetic(*hdr['position']) except (KeyError, ValueError): pass # %% observation types try: hdr['fields'] = hdr['# / TYPES OF OBSERV'] if Nobs != len(hdr['fields']): raise ValueError(f'{f.name} header read incorrectly') if isinstance(meas, (tuple, list, np.ndarray)): ind = np.zeros(len(hdr['fields']), dtype=bool) for m in meas: for i, f in enumerate(hdr['fields']): if f.startswith(m): ind[i] = True hdr['fields_ind'] = np.nonzero(ind)[0] else: ind = slice(None) hdr['fields_ind'] = np.arange(Nobs) hdr['fields'] = np.array(hdr['fields'])[ind].tolist() except KeyError: pass hdr['Nobsused'] = hdr['Nobs'] if useindicators: hdr['Nobsused'] *= 3 # %% try: hdr['# OF SATELLITES'] = int(hdr['# OF SATELLITES'][:6]) except (KeyError, ValueError): pass # %% time try: hdr['t0'] = _timehdr(hdr['TIME OF FIRST OBS']) except (KeyError, ValueError): pass try: hdr['t1'] = _timehdr(hdr['TIME OF LAST OBS']) except (KeyError, ValueError): pass try: # This key is OPTIONAL hdr['interval'] = float(hdr['INTERVAL'][:10]) except (KeyError, ValueError): pass return hdr
def obsheader2(f: TextIO, useindicators: bool = False, meas: Sequence[str] = None) -> Dict[str, Any]: """ End users should use rinexheader() """ if isinstance(f, (str, Path)): with opener(f, header=True) as h: return obsheader2(h, useindicators, meas) f.seek(0) # %% selection if isinstance(meas, str): meas = [meas] if not meas or not meas[0].strip(): meas = None hdr = rinexinfo(f) Nobs = 0 # not None due to type checking for ln in f: if "END OF HEADER" in ln: break h = ln[60:80].strip() c = ln[:60] # %% measurement types if "# / TYPES OF OBSERV" in h: if Nobs == 0: Nobs = int(c[:6]) hdr[h] = c[6:].split() else: hdr[h] += c[6:].split() elif h not in hdr: # Header label hdr[h] = c # string with info else: # concatenate hdr[h] += " " + c # %% useful values try: hdr["systems"] = hdr["RINEX VERSION / TYPE"][40] except KeyError: pass hdr["Nobs"] = Nobs # 5 observations per line (incorporating LLI, SSI) hdr["Nl_sv"] = ceil(hdr["Nobs"] / 5) # %% list with receiver location in x,y,z cartesian ECEF (OPTIONAL) try: hdr["position"] = [ float(j) for j in hdr["APPROX POSITION XYZ"].split() ] if ecef2geodetic is not None: hdr["position_geodetic"] = ecef2geodetic(*hdr["position"]) except (KeyError, ValueError): pass # %% observation types try: hdr["fields"] = hdr["# / TYPES OF OBSERV"] if hdr["Nobs"] != len(hdr["fields"]): logging.error( f"{f.name} number of observations declared in header does not match fields" ) hdr["Nobs"] = len(hdr["fields"]) if isinstance(meas, (tuple, list, np.ndarray)): ind = np.zeros(len(hdr["fields"]), dtype=bool) for m in meas: for i, f in enumerate(hdr["fields"]): if f.startswith(m): ind[i] = True hdr["fields_ind"] = np.nonzero(ind)[0] else: ind = np.s_[:] hdr["fields_ind"] = np.arange(hdr["Nobs"]) hdr["fields"] = np.array(hdr["fields"])[ind].tolist() except KeyError: pass hdr["Nobsused"] = hdr["Nobs"] if useindicators: hdr["Nobsused"] *= 3 # %% try: hdr["# OF SATELLITES"] = int(hdr["# OF SATELLITES"][:6]) except (KeyError, ValueError): pass # %% time try: hdr["t0"] = _timehdr(hdr["TIME OF FIRST OBS"]) except (KeyError, ValueError): pass try: hdr["t1"] = _timehdr(hdr["TIME OF LAST OBS"]) except (KeyError, ValueError): pass try: # This key is OPTIONAL hdr["interval"] = float(hdr["INTERVAL"][:10]) except (KeyError, ValueError): pass try: s = " " hdr["rxmodel"] = s.join(hdr["REC # / TYPE / VERS"].split()[1:-1]) except (KeyError, ValueError): pass return hdr
def obsheader3(f: TextIO, use: Sequence[str] = None, meas: Sequence[str] = None) -> Dict[str, Any]: """ get RINEX 3 OBS types, for each system type optionally, select system type and/or measurement type to greatly speed reading and save memory (RAM, disk) """ if isinstance(f, (str, Path)): with opener(f, header=True) as h: return obsheader3(h, use, meas) fields = {} Fmax = 0 # %% first line hdr = rinexinfo(f) for ln in f: if "END OF HEADER" in ln: break h = ln[60:80] c = ln[:60] if 'SYS / # / OBS TYPES' in h: k = c[0] fields[k] = c[6:60].split() N = int(c[3:6]) # %% maximum number of fields in a file, to allow fast Numpy parse. Fmax = max(N, Fmax) n = N-13 while n > 0: # Rinex 3.03, pg. A6, A7 ln = f.readline() assert 'SYS / # / OBS TYPES' in ln[60:] fields[k] += ln[6:60].split() n -= 13 assert len(fields[k]) == N continue if h.strip() not in hdr: # Header label hdr[h.strip()] = c # don't strip for fixed-width parsers # string with info else: # concatenate to the existing string hdr[h.strip()] += " " + c # %% list with x,y,z cartesian (OPTIONAL) try: hdr['position'] = [float(j) for j in hdr['APPROX POSITION XYZ'].split()] if ecef2geodetic is not None: hdr['position_geodetic'] = ecef2geodetic(*hdr['position']) except (KeyError, ValueError): pass # %% time try: t0s = hdr['TIME OF FIRST OBS'] # NOTE: must do second=int(float()) due to non-conforming files hdr['t0'] = datetime(year=int(t0s[:6]), month=int(t0s[6:12]), day=int(t0s[12:18]), hour=int(t0s[18:24]), minute=int(t0s[24:30]), second=int(float(t0s[30:36])), microsecond=int(float(t0s[30:43]) % 1 * 1000000)) except (KeyError, ValueError): pass try: hdr['interval'] = float(hdr['INTERVAL'][:10]) except (KeyError, ValueError): pass # %% select specific satellite systems only (optional) if use is not None: if not set(fields.keys()).intersection(use): raise KeyError(f'system type {use} not found in RINEX file') fields = {k: fields[k] for k in use if k in fields} # perhaps this could be done more efficiently, but it's probably low impact on overall program. # simple set and frozenset operations do NOT preserve order, which would completely mess up reading! sysind = {} if isinstance(meas, (tuple, list, np.ndarray)): for sk in fields: # iterate over each system # ind = np.isin(fields[sk], meas) # boolean vector ind = np.zeros(len(fields[sk]), dtype=bool) for m in meas: for i, f in enumerate(fields[sk]): if f.startswith(m): ind[i] = True fields[sk] = np.array(fields[sk])[ind].tolist() sysind[sk] = np.empty(Fmax*3, dtype=bool) # *3 due to LLI, SSI for j, i in enumerate(ind): sysind[sk][j*3:j*3+3] = i else: sysind = {k: slice(None) for k in fields} hdr['fields'] = fields hdr['fields_ind'] = sysind hdr['Fmax'] = Fmax return hdr
def obsheader3(f: TextIO, use: Sequence[str] = None, meas: Sequence[str] = None) -> Dict[str, Any]: """ get RINEX 3 OBS types, for each system type optionally, select system type and/or measurement type to greatly speed reading and save memory (RAM, disk) """ if isinstance(f, (str, Path)): with opener(f, header=True) as h: return obsheader3(h, use, meas) fields = {} Fmax = 0 # %% first line hdr = rinexinfo(f) for ln in f: if "END OF HEADER" in ln: break h = ln[60:80] c = ln[:60] if 'SYS / # / OBS TYPES' in h: k = c[0] fields[k] = c[6:60].split() N = int(c[3:6]) # %% maximum number of fields in a file, to allow fast Numpy parse. Fmax = max(N, Fmax) n = N - 13 while n > 0: # Rinex 3.03, pg. A6, A7 ln = f.readline() assert 'SYS / # / OBS TYPES' in ln[60:] fields[k] += ln[6:60].split() n -= 13 assert len(fields[k]) == N continue if h.strip() not in hdr: # Header label hdr[h.strip()] = c # don't strip for fixed-width parsers # string with info else: # concatenate to the existing string hdr[h.strip()] += " " + c # %% list with x,y,z cartesian (OPTIONAL) try: hdr['position'] = [ float(j) for j in hdr['APPROX POSITION XYZ'].split() ] if ecef2geodetic is not None: hdr['position_geodetic'] = ecef2geodetic(*hdr['position']) except (KeyError, ValueError): pass # %% time try: t0s = hdr['TIME OF FIRST OBS'] # NOTE: must do second=int(float()) due to non-conforming files hdr['t0'] = datetime(year=int(t0s[:6]), month=int(t0s[6:12]), day=int(t0s[12:18]), hour=int(t0s[18:24]), minute=int(t0s[24:30]), second=int(float(t0s[30:36])), microsecond=int(float(t0s[30:43]) % 1 * 1000000)) except (KeyError, ValueError): pass try: hdr['interval'] = float(hdr['INTERVAL'][:10]) except (KeyError, ValueError): pass # %% select specific satellite systems only (optional) if use is not None: if not set(fields.keys()).intersection(use): raise KeyError(f'system type {use} not found in RINEX file') fields = {k: fields[k] for k in use if k in fields} # perhaps this could be done more efficiently, but it's probably low impact on overall program. # simple set and frozenset operations do NOT preserve order, which would completely mess up reading! sysind = {} if isinstance(meas, (tuple, list, np.ndarray)): for sk in fields: # iterate over each system # ind = np.isin(fields[sk], meas) # boolean vector ind = np.zeros(len(fields[sk]), dtype=bool) for m in meas: for i, f in enumerate(fields[sk]): if f.startswith(m): ind[i] = True fields[sk] = np.array(fields[sk])[ind].tolist() sysind[sk] = np.empty(Fmax * 3, dtype=bool) # *3 due to LLI, SSI for j, i in enumerate(ind): sysind[sk][j * 3:j * 3 + 3] = i else: sysind = {k: slice(None) for k in fields} hdr['fields'] = fields hdr['fields_ind'] = sysind hdr['Fmax'] = Fmax return hdr
def obsheader2(f: TextIO, useindicators: bool = False, meas: Sequence[str] = None) -> Dict[str, Any]: """ End users should use rinexheader() """ if isinstance(f, (str, Path)): with opener(f, header=True) as h: return obsheader2(h, useindicators, meas) f.seek(0) # %% selection if isinstance(meas, str): meas = [meas] if not meas or not meas[0].strip(): meas = None hdr = rinexinfo(f) Nobs = 0 # not None due to type checking for ln in f: if "END OF HEADER" in ln: break h = ln[60:80].strip() c = ln[:60] # %% measurement types if '# / TYPES OF OBSERV' in h: if Nobs == 0: Nobs = int(c[:6]) hdr[h] = c[6:].split() else: hdr[h] += c[6:].split() elif h not in hdr: # Header label hdr[h] = c # string with info else: # concatenate hdr[h] += " " + c # %% useful values try: hdr['systems'] = hdr['RINEX VERSION / TYPE'][40] except KeyError: pass hdr['Nobs'] = Nobs # 5 observations per line (incorporating LLI, SSI) hdr['Nl_sv'] = ceil(hdr['Nobs'] / 5) # %% list with receiver location in x,y,z cartesian ECEF (OPTIONAL) try: hdr['position'] = [float(j) for j in hdr['APPROX POSITION XYZ'].split()] if ecef2geodetic is not None: hdr['position_geodetic'] = ecef2geodetic(*hdr['position']) except (KeyError, ValueError): pass # %% observation types try: hdr['fields'] = hdr['# / TYPES OF OBSERV'] if Nobs != len(hdr['fields']): raise ValueError(f'{f.name} header read incorrectly') if isinstance(meas, (tuple, list, np.ndarray)): ind = np.zeros(len(hdr['fields']), dtype=bool) for m in meas: for i, f in enumerate(hdr['fields']): if f.startswith(m): ind[i] = True hdr['fields_ind'] = np.nonzero(ind)[0] else: ind = slice(None) hdr['fields_ind'] = np.arange(Nobs) hdr['fields'] = np.array(hdr['fields'])[ind].tolist() except KeyError: pass hdr['Nobsused'] = hdr['Nobs'] if useindicators: hdr['Nobsused'] *= 3 # %% try: hdr['# OF SATELLITES'] = int(hdr['# OF SATELLITES'][:6]) except (KeyError, ValueError): pass # %% time try: hdr['t0'] = _timehdr(hdr['TIME OF FIRST OBS']) except (KeyError, ValueError): pass try: hdr['t1'] = _timehdr(hdr['TIME OF LAST OBS']) except (KeyError, ValueError): pass try: # This key is OPTIONAL hdr['interval'] = float(hdr['INTERVAL'][:10]) except (KeyError, ValueError): pass return hdr