First commit of package

Author: geomsi

LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2019 Mustafa Serkan Işık and Volkan Özbey
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

MANIFEST.in

@@ -0,0 +1,2 @@
+include LICENSE
+include README.md

README.md

@@ -0,0 +1,67 @@
+# GNSSpy
+Python Toolkit for GNSS Data developed by Mustafa Serkan Isik (isikm@itu.edu.tr) and Volkan Ozbey (ozbeyv@itu.edu.tr). This project is still a work in progress. Send us your feedback if possible...
+
+# What is GNSSpy? 
+GNSSpy is a free and open source library for handling multi GNSS and different versions (2.X and 3.X) of RINEX files. It provides
+Single Point Positioning (SPP) solutions by least squares adjustment using pseudo-range observations using precise ephemeris and clock files. GNSSpy can be used for editing (slicing, decimating, merging) and quality checking (multipath,ionospheric delay, SNR) for RINEX files. Ionospheric delay can be calculated from GNSS atmospheric models of IGS for single frequency RINEX data or removed using dual frequency RINEX data. It can be used for visualizing GNSS data such as skyplot, azimuth-elevation,time-elevation, ground track and band plots. Additionally, this library can be used for basic geodetic computations such as geodetic positions on reference ellipsoid and projection computations.
+
+# How to install?
+Download the package and change directory of your terminal to gnsspy-master folder. Then, simply type
+```
+pip setup.py install
+```
+Or you can directly install package from GitHub via
+```
+pip install git+https://github.com/GNSSpy-Project/gnsspy
+```
+# How to use?
+A detailed version of manual will be released soon.
+
+## Read RINEX Observation File
+`read_obsFile` function reads RINEX 2.x/3.x observation files. If the station is IGS station, RINEX file does not necessarily need to be exist in working directory. In that case, the file is automatically downloaded to working directory. 
+```python
+import gnsspy as gp
+station = gp.read_obsFile("mate2440.17o")
+```
+`read_obsFile` function returns to a class instance. This instance
+```python
+# Epoch of RINEX file as datetime
+station.epoch
+# Pandas.DataFrame of observations
+station.observation
+# Approximate position [type:list-> x,y,z]
+station.approx_position
+# Antenna Type
+station.antenna_type
+# Observation interval(seconds)
+station.interval
+# Receiver clock error
+# (if available)
+station.receiver_clock
+# Receiver Type
+station.receiver_type
+# RINEX version
+staton.version
+# RINEX filename
+station.filename
+```
+
+## Interpolation of SP3 Final Products
+`sp3_interp` function interpolates final precise orbit coordinates of satellites at RINEX observation epochs. Default interpolation method is 16 degree polynomial interpolation. Degree of polynomial can be changed, though it is not recommended to use lower than 11 degree. Above 16 degree is not applicable for 15 minute intervals of precise orbit solution. GFZ orbit and clock files are default product names. Alternatives are IGS, WUM, ESA etc. Of course, each product provides solution for different satellite systems, hence number of satellites may vary for each product choice.
+```python
+orbit = gp.sp3_interp(station.epoch, interval=station.interval, poly_degree=16, sp3_product="gfz", clock_product="gfz")
+```
+## Single Point Positioning (SPP)
+In order to use `spp` function, station file must be read and SP3 interpolation must be done. Satellite system can be chosed via `system` argument as in the example. G: GPS - R: GLONASS - E: GALILEO - C: COMPASS - J: QZSS - I: IRNSS. It should be noted that the choice of station file and products used for the SP3 interpolation can constrain satellite system selection. Additionally, elevation mask angle (cut_off) is by default 7.0 degree. In next release, more options will be added to this function.
+```python
+spp_result = gp.spp(station, orbit, system="G", cut_off=7.0)
+```
+
+# Notes
+crx2rnx function is not pure python implimentation and depends on 
+RNXCMP software for compression/restoration of RINEX observation files 
+developed by Y. Hatanaka of GSI.
+
+Source: http://terras.gsi.go.jp/ja/crx2rnx.html
+
+Reference: Hatanaka, Y. (2008): A Compression Format and Tools for GNSS Observation Data, Bulletin of the Geographical Survey  Institute, 55, 21-30, available at http://www.gsi.go.jp/ENGLISH/Bulletin55.html

gnsspy/__init__.py

@@ -0,0 +1,37 @@
+"""
+Copyright (c) 2019 Mustafa Serkan Işık and Volkan Özbey
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.
+"""
+
+from gnsspy.io.readFile import (read_navFile, read_obsFile, read_sp3File,
+                                read_clockFile, read_ionFile)
+from gnsspy.io.manipulate import (rinex_merge, crx2rnx)
+from gnsspy.position.atmosphere import (tropospheric_delay)
+from gnsspy.position.interpolation import (sp3_interp, ionosphere_interp)
+from gnsspy.position.position import (spp, multipath)
+from gnsspy.geodesy import (coordinate, projection)
+from gnsspy.funcs.funcs import (gpsweekday, gpswdtodate, jday, julianday2date,
+                                doy, doy2date, datetime2doy)
+from gnsspy.download import (get_rinex, get_rinex3, get_navigation, 
+                             get_sp3, get_clock, get_ionosphere)
+from gnsspy import plot
+
+__version__ = '2019.08.02'
+__author__ = 'Mustafa Serkan Isik & Volkan Ozbey'

gnsspy/doc/IGS.py

@@ -0,0 +1,21 @@
+# ===========================================================
+# ========================= imports =========================
+import os
+import pandas as _pd
+# ===========================================================
+
+global _FILEPATH
+_FILEPATH = os.path.dirname(os.path.abspath(__file__))
+
+def is_IGS(siteName):
+    igs = _pd.read_table(_FILEPATH+"/IGSList.txt",sep="\t", names = ["CODE","SITE","COUNTRY","LATITUDE","LONGITUDE","HEIGHT"], index_col = "CODE")
+    if siteName.upper() not in igs.index.tolist():
+        return False
+    return True
+
+def IGS(siteName):
+    igs = _pd.read_table(_FILEPATH+"/IGSList.txt",sep="\t", names = ["CODE","SITE","COUNTRY","LATITUDE","LONGITUDE","HEIGHT"], index_col = "CODE")
+    if is_IGS(siteName) == True:
+        siteInfo = igs[igs.index==siteName.upper()]
+        return siteInfo
+