Skip to content

Commit 5cb6970

Browse files
committed
Completed vertical datums tutorial
1 parent d2e1052 commit 5cb6970

File tree

4 files changed

+269
-6
lines changed

4 files changed

+269
-6
lines changed

docs/features/height.rst

Lines changed: 43 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,4 +1,46 @@
11
.. _features/height:
22

33
Vertical Datums
4-
=================
4+
=================
5+
6+
This module includes functions for dealing with height datums and geoids.
7+
8+
.. note:: GDAL is needed for this module. Please consult the GDAL pypi for instructions on installing GDAL
9+
10+
Converting Between Height Datums
11+
--------------------------------
12+
13+
The following functions all convert between different height datums with a particular focus on Australian datums.
14+
15+
.. autofunction:: geodepy.height.GPS_to_AHD
16+
.. autofunction:: geodepy.height.AHD_to_GPS
17+
.. autofunction:: geodepy.height.GPS_to_AVWS
18+
.. autofunction:: geodepy.height.AVWS_to_GPS
19+
.. autofunction:: geodepy.height.AHD_to_AVWS
20+
.. autofunction:: geodepy.height.AVWS_to_AHD
21+
.. autofunction:: geodepy.height.GPS_to_AUSGeoid98
22+
.. autofunction:: geodepy.height.AUSGeoid98_to_GPS
23+
.. autofunction:: geodepy.height.GPS_to_AUSGeoid09
24+
.. autofunction:: geodepy.height.AUSGeoid09_to_GPS
25+
.. autofunction:: geodepy.height.DOV
26+
.. autofunction:: geodepy.height.DOV_09
27+
.. autofunction:: geodepy.height.DOV_98
28+
29+
Calculating Gravity
30+
-------------------
31+
32+
The functions below can be used to calculate different components of gravity.
33+
34+
.. autofunction:: geodepy.height.mean_surface_grav
35+
.. autofunction:: geodepy.height.interp_grav
36+
.. autofunction:: geodepy.height.normal_grav
37+
.. autofunction:: geodepy.height.mean_normal_grav
38+
.. autofunction:: geodepy.height.normal_correction
39+
.. autofunction:: geodepy.height.normal_orthometric_correction
40+
41+
Auxilary Function
42+
-----------------
43+
44+
Functions used to enable other functions in the module.
45+
46+
.. autofunction:: geodepy.height.interp_file

docs/requirements.txt

Lines changed: 4 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,3 +1,6 @@
11
furo
22
pandas
3-
sphinxext-opengraph
3+
sphinxext-opengraph
4+
gdal
5+
numpy>1.0.0 wheel setuptools>=67
6+
gdal[numpy]=="$(gdal-config --version).*"
Lines changed: 65 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,4 +1,68 @@
11
.. _tutorials/vert:
22

33
Vertical Datums
4-
===============
4+
===============
5+
6+
This tutorial will work through converting between different vertical datums and working with gravity.
7+
8+
Converting Between Vertical Datums
9+
----------------------------------
10+
11+
Geodepy allows you to convert between different vertical datums within Australia. Below we will convert from
12+
an ellipsoidal height to AHD and then to AVWS.
13+
14+
.. code:: python
15+
16+
import geodepy.height
17+
18+
First we will convert from ellipsoidal height to AHD.
19+
20+
.. code:: python
21+
22+
rl = geodepy.height.GPS_to_AHD(-35.34405212, 149.15847673, 594.495)
23+
24+
print(rl[0])
25+
26+
>>[575.176]
27+
28+
This is the AHD height of this point. Now we can convert to AVWS.
29+
30+
.. code:: python
31+
32+
rl = geodepy.height.AHD_to_AVWS(-35.34405212, 149.15847673,rl[0])
33+
34+
print(rl[0])
35+
36+
>>[575.303]
37+
38+
Finding Gravity Values
39+
----------------------
40+
41+
Geodepy can also find gravity values at any point around Australia. This can be seen below.
42+
43+
.. code:: python
44+
45+
import geodepy.height
46+
47+
To find the gravity at any point the geodepy.height.mean_surface_grav function can be used.
48+
Normally this function is used to find the mean surface gravty between two points however
49+
if the same lat and long are used for both point A and B then the gravity at a set location can be found.
50+
51+
.. code:: python
52+
53+
grav = geodepy.height.mean_surface_grav(-35.34405212, 149.15847673, 575.176, -35.34405212, 149.15847673, 575.176,)
54+
55+
print(grav)
56+
57+
>>[9.79603591]
58+
59+
To calculate the average gravity between the ellipsoid and a certain height the following function can be used.
60+
61+
.. code:: python
62+
63+
grav = geodepy.height.mean_normal_grav(-35.34405212, 594.495)
64+
65+
print(grav)
66+
67+
>>9.79671297
68+

0 commit comments

Comments
 (0)