How do I use mapproject to project onto a line?

Shahar Barak <[email protected]>
Newsgroups gmane.comp.gis.gmt.user
Message-ID <CAL6UvMoXRxC=J4yFW=u1as8wBKYjcHiJeyoGQWcayGSgvPNLbg@mail.gmail.com>
1. I need to project some points onto a line, but "project" is not giving
me accurate results.
The documentation says:
"Calculations of specific great-circle and geodesic distances or for
back-azimuths or azimuths are better done using mapproject."

However, I don't see how to a project onto a line in "mapproject".

2. Also, how come reversing the end-points in "project" gives me different
results?

echo "-117.000001 33" | /usr/local/gmt32bit/bin/project -C-118/33 -E-117/33
-Lw -Fxypqrs

-117.000001 33 0.838666571653 -3.9861033073e-09 -117.000001 33.000000004

echo "-117.000001 33" | /usr/local/gmt32bit/bin/project -C-117/33 -E-118/33
-Lw -Fxypqrs

-117.000001 33 0 3.9861033073e-09 -117 33

How come p=~0.84 and not ~0.999999 in the first case?

To unsubscribe, send the message "signoff gmt-help" to [email protected]
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.