File: PyGetIntersect.py

package info (click to toggle)
geographiclib 2.7-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 7,572 kB
  • sloc: cpp: 27,765; sh: 5,463; makefile: 695; python: 12; ansic: 10
file content (16 lines) | stat: -rw-r--r-- 437 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
import ctypes

intersect_dso = ctypes.CDLL('./intersect_dso.so')
intersect_c = intersect_dso.intersect
intersect_c.argtypes = ([ctypes.c_double] * 6) + ([ctypes.POINTER(ctypes.c_double)] * 2)

def intersect(lat1, lon1, az1, lat2, lon2, az2):
    dX = ctypes.c_double()
    dY = ctypes.c_double()

    intersect_c(
        lat1, lon1, az1, lat2, lon2, az2,
        ctypes.byref(dX), ctypes.byref(dY),
    )

    return dX.value, dY.value