1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
|
#!/usr/bin/env python
from pymap3d.vincenty import vreckon
from argparse import ArgumentParser
def main():
p = ArgumentParser(
description="Given starting latitude, longitude: find final lat,lon for distance and azimuth"
)
p.add_argument("lat", help="latitude WGS-84 [degrees]", type=float)
p.add_argument("lon", help="longitude WGS-84 [degrees]", type=float)
p.add_argument("range", help="range from start point [meters]", type=float)
p.add_argument("azimuth", help="clockwise from north: azimuth to start [degrees]", type=float)
P = p.parse_args()
lat2, lon2 = vreckon(P.lat, P.lon, P.range, P.azimuth)
print("lat, lon = ({:.4f}, {:.4f})".format(lat2, lon2))
if __name__ == "__main__": # pragma: no cover
main()
|