<font size=2 face="sans-serif">Dear proj-members,</font>
<br>
<br><font size=2 face="sans-serif">is there a way to read the vertical
deflection (xi, eta) directly with proj?</font>
<br>
<br><font size=2 face="sans-serif">My work-around is sampling the geoid-grid
shift:</font>
<br>
<br><font size=2 face="sans-serif">```python</font>
<br><font size=3 color=#000080 face="Courier New"><b>def </b></font><font size=3 face="Courier New">get_vertical_deflection(lon,
lat):<br>
    </font><font size=3 color=#808080 face="Courier New"><i>"""<br>
    take origin and 1 east + 1 north on GEOID, transform to topocentric
and check<br>
    </i><b>:param</b><i> lon: longitude of point of interest
on WGS84 (°)<br>
    </i><b>:param</b><i> lat: latitude of point of interest on
WGS84 (°)<br>
    </i><b>:return</b><i>: (xi, eta) deflection in meridional
and prime vertical (rad)<br>
    """<br>
    </i></font><font size=3 face="Courier New">eps = </font><font size=3 color=blue face="Courier New">0.0001
 </font><font size=3 color=#808080 face="Courier New"><i># °<br>
    </i></font><font size=3 face="Courier New">lons = np.array([lon,
lon + eps, lon])<br>
    lats = np.array([lat, lat,       lat + eps])<br>
    geoid_heights = get_geoid_height_egm96(lons, lats)<br>
<br>
    </font><font size=3 color=#808080 face="Courier New"><i>#
W. E. Featherstone, "The Use and Abuse of Vertical Deflection"
p. 4<br>
    </i></font><font size=3 face="Courier New">e = eccentricity(f_wgs84)<br>
    phi = np.deg2rad(lat)<br>
<br>
    </font><font size=3 color=#808080 face="Courier New"><i>#
meridional, north/south<br>
    # R_M, should be radius of curvature in the meridian at point
of interest<br>
    </i></font><font size=3 face="Courier New">rad_meridian =
get_radius_of_curvature_meridian(a_wgs84, e, phi)<br>
    </font><font size=3 color=#808080 face="Courier New"><i>#
d_North = R d_lat; in rad<br>
    </i></font><font size=3 face="Courier New">xi = -(geoid_heights[</font><font size=3 color=blue face="Courier New">2</font><font size=3 face="Courier New">]
- geoid_heights[</font><font size=3 color=blue face="Courier New">0</font><font size=3 face="Courier New">])/(rad_meridian
* np.deg2rad(eps))<br>
<br>
    </font><font size=3 color=#808080 face="Courier New"><i>#
prime vertical, east/west<br>
    # R_N, should be radius of curvature in the prime vertical
at point of interest<br>
    </i></font><font size=3 face="Courier New">rad_prime_vertical
= get_radius_of_curvature_prime_vertical(a_wgs84, e, phi)<br>
    </font><font size=3 color=#808080 face="Courier New"><i>#
d_East = R cos(lat) d_long; in rad<br>
    </i></font><font size=3 face="Courier New">eta = -(geoid_heights[</font><font size=3 color=blue face="Courier New">1</font><font size=3 face="Courier New">]
- geoid_heights[</font><font size=3 color=blue face="Courier New">0</font><font size=3 face="Courier New">])/(rad_prime_vertical
* np.deg2rad(eps) * np.cos(phi))<br>
    </font><font size=3 color=#000080 face="Courier New"><b>return
</b></font><font size=3 face="Courier New">xi, eta</font>
<br>
<br><font size=2 face="sans-serif">```</font>
<br><font size=2 face="sans-serif">Over which (arc) distance should I sample
the geoid gtx?</font>
<br>
<br><font size=2 face="Verdana">Kind regards</font>
<br>
<br><font size=2 face="Verdana">Fabian Gross<br>
Telefon +49 (711) 648 71-995</font><font size=1 face="Verdana"><br>
_________________________________________________</font><font size=5 face="Verdana"><b><br>
sbp</b> </font><font size=1 face="Verdana"><b><br>
schlaich<br>
bergermann partner</b><br>
<br>
Beratende Ingenieure<br>
für erneuerbare Energie<br>
<br>
Stuttgart . Berlin . New York<br>
São Paulo . Shanghai . Paris<br>
<b><br>
sbp sonne gmbh</b><br>
<br>
Markus Balz Dipl.Ing. (FH)<br>
Andreas Keil Dipl.Ing.<br>
<br>
Schwabstrasse 43<br>
70197 Stuttgart<br>
Telefon +49 (711) 648 71-0<br>
</font><font size=1 color=blue face="Verdana"><u><br>
</u></font><a href=http://www.sbp.de/><font size=1 face="Verdana"><u>www.sbp.de</u></font></a><font size=1 face="Verdana"><br>
_________________________________________________</font>