geoscienceaustralia/geodepy

A toolkit for Geodesy and Surveying in Python

This is a package of tools for manipulating geospatial datasets using Python and tested in Python 3.6.4.

Tutorials

See here for worked examples of common GeodePy functions and routines.

Dependencies

This package requires the following PyPI Packages installed:

NumPy
SciPy

Additionally, the geodepy.height module requires the GDAL library (tested using GDAL 3.0.4). For more information, see here for information about GDAL, here for Anaconda support for GDAL and here for GDAL Binaries for Windows.

Testing

Run: python -m unittest discover geodepy/tests/ --verbose

API

cd api/
virtualenv env
source env/bin/activate
pip install -r requirements.txt
zappa deploy dev

For subsequent updating run: zappa update dev

Authors

• Craig Harrison - Project Management - harry093
• Josh Batchelor - Initial Work, Geodesy and Surveying - BatchelorJ
• Jonathan Mettes - Testing, Integration and Deployment - jmettes
• Jack McCubbine - Height Module - JackMcCubbineGA

See also the list of contributors who participated in this project.

License

Copyright 2018-2020 Geoscience Australia

Licensed under the Apache License, Version 2.0 (the "License"); you may not use this file except in compliance with the License. You may obtain a copy of the License at

http://www.apache.org/licenses/LICENSE-2.0

Unless required by applicable law or agreed to in writing, software distributed under the License is distributed on an "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. See the License for the specific language governing permissions and limitations under the License.

geodepy's Issues

gpscal

dms2dd bug

dms2dd is outputting incorrect results in limited examples.

Found with input 133deg 48min 00sec, giving incorrect output of 133.81111111111107decdeg instead of correct output of 133.8decdeg.

I believe this is an issue with rounding in the output of modf parts of the function used to separate out the degrees, minutes and seconds components of the input.

rdsinex

Read & Search Geoid Data

Currently geoid data is stored in ASCII .dat and NTv2 .gsb formats. This module requires a way to read this data into a Python data structure (e.g. a NumPy array or PyTable object?) and search through to find specific latitude and longitude values.

MGA94 <-> MGA2020

It would be good to create a function that strings together all the other necessary functions to do this. Emails to geodesy@ga are making it obvious that people think this is a simple transformation. Required steps:

MGA94<->GDA94 lat/lon<->GDA94 XYZ<->GDA2020 XYZ<->GDA2020 lat/lon<->MGA2020

Vectors: Local to Cartesian

• xyz2enu
• enu2xyz

Add test coverage to statistics module

Move contents of ntv2reader module to the geoid module

Typos in convert.py

Non-issue

----

We need to develop this

Add ellipsoid parameters

These can be stored as dictionaries. Would require rewriting dependencies from grs80 in multiple pieces of code

Shift tests from test_tranform.py to test_convert.py

Functions grid2geo, geo2grid, llh2xyz and xyz2llh were shifted from transform to convert. Their tests need to be moved from test_transform to test_convert

Reference frame transformations

See Table 4.2 in the DynaNet User's Guide v3.3

Other geodesy related python packages

A list of existing geodesy (or related) python packages:

• https://geographiclib.sourceforge.io/

Bilinear interpolator

Move ntv2reader module tests

At the bottom of the ntv2reader module there is a test. This should be moved to a file in the geodepy/tests/ directory.

The existing test has a hardcoded path in it. Should this be removed?

Allow epochs in decimal years

Currently, only datetime.date objects are accepted to function that require an epoch. Extend code to accept decimal years

Grid Interpolation

Requires #16
Uses data from transformation grid files to interpolate (bilinear and/or bicubic) between grid points and return a set of values at a specified point.

vincdir: Geodesics crossing 180th meridian

Currently, where the geodesic path crosses the 180th meridian, the resulting long2 will be > 180°. Possibly this issue may affect geodesics crossing the prime meridian resulting in long2 < 0°.

sinex2epoch

Integer Negative Zero issue in DMSAngle and DDMAngle classes

When entering the degree component of a DMSAngle(degree, minute, second) or DDMAngle(degree, minute) object, if the degree component is between -1 and 0 (i.e. DMSAngle(-0,12,14)) this will incorrectly convert this to a positive value. This is because python interprets the input '-0' as an int with value 0 and does not retain the sign of this variable.

Tidy up convert.py

• Move dec2hp and hp2dec to below the class definitions
• Remove legacy functions (use PyCharm to find usages, if any, first)
• Ask Jack about dd2dms_v and dms2dd_v
• Move class DNACoord up with the other classes

Conformal Transformations: 14 Parameter

• conform14 - Functionality
• Validation and Unit Testing

Surveying observation equations

• Slope distance (S)
• Ellipsoid arc distance (E)
• Ellipsoid chord distance (C)
• Mean sea level arc distance (M)
• Ellipsoid heights (R) and height differences
• Orthometric heights (H)
• Orthometric height differences (L)
• Cartesian coordinates and GNSS point clusters (Y)
• GNSS baselines (G)
• GNSS baseline clusters (X)
• Geodetic latitude (P) and longitude (Q)
• Geodetic azimuths and horizontal bearings (B)
• Horizontal angles (A)
• Horizontal direction sets (D)
• Zenith distances (V)
• Vertical angles (Z)
• Astronomic latitude (I) and longitude (J)
• Astronomic (Laplace) azimuth (K)

Cartesian to Geographic

• xyz2llh
• llh2xyz

Review GDAL Dependency

The package dependency on GDAL is making installation with pip a bit of a nightmare. Is GDAL required, or is there an alternative pure-python package we could use?

Survey Conversions: Leica GSI to GeodesyML

• gsi2geodesyml
  Tool to convert Leica GSI files to GeodesyML format ready for import and processing in DynaNet.

Add methods from old tech manuals

For example, Lauff's Method

Polar to local conversion

Provide test coverage for the statistics module

• rotation_matrix()
• vcv_cart2local()
• vcv_local2cart2()
• error_ellipse()
• circ_hz_pu()
• k_val95()

Add ATRF to GDA2020 parameters

Can you please updates the Constant.py script to include ATRF to GDA2020 parameters. This should be a copy and paste of the ITRF2014 to GDA2020 parameters?

Please also include the comment that the parameters (for both) are derived from the Australian Plate Motion Model based on 109 reference stations in Australia.

Vincenty's Formulae

• vincdir
• vincinv
  .csv i/o functions:
• vincdirio
• vincinvio

Vectors: Local to Cartesian

• xyz2enu
• enu2xyz

Tidy up transform module

Remove grid2geoio() & geo2gridio() (use PyCharm to find usages, if any, first)

Geographic to Grid Conversions

• geo2grid
• grid2geo

.csv i/o functions:

• geo2gridio
• grid2geoio

Intersection of two lines

Given two lines, defined by a pair of lat/lons, calculate the point of intersection

Bicubic interpolator

Tidy up constants.py

Remove c_vac and k_0 - these were only ever added to give the initial version of the module some content :-)

Conformal Transformations: 7 Parameter

• conform7

----

• Conform7

Survey Conversion: Liscad Field File to GeodesyML

• lisfield2geodesyml
  Tool to convert Liscad Field file to GeodesyML format ready for import and processing in DynaNet.

Tidy up geodesy.py

Remove vindirio() & vininvio() (use PyCharm to find usages, if any, first)

Simplify Angular Conversions

There are currently 8 different functions for converting between angular data formats. Some of these duplicate functionality. These functions should be rewritten and renamed to avoid confusion.

Local to polar conversion

Add test coverage of gnss module

----

Github is shit

Variance reference frame transformations

• Cartesian <-> local
• Polar <-> local
• Geographical <-> Cartesian

Modify vcv_cart2local and vcv_local2cart to return a 3x1 matrix when necessary

vcv_cart2local and vcv_local2cart were modified to accept a 3x1 matrix of variances as well as the full 3x3 VCV. However, the code only returns the 3x3 matrix

Add Helmert Transformation Parameters

These are being added into the Constants module as a Transformation Class object

• GDA94 to GDA2020
• ITRF14 to GDA2020
  From Dawson and Woods
• ITRF08 to GDA94
• ITRF05 to GDA94
• ITRF00 to GDA94
• ITRF97 to GDA94
• ITRF96 to GDA94
  From GDA94 Tech Manual
• AGD66/84 to GDA94
  From IERS
• All current ITRFxx to ITRFyy

Shift llh2xyz, xyz2llh, grid2geo, and geo2grid from transform module to convert module

These four functions are coordinate conversions within the same datum and are better placed in the convert module