[Proj] (no subject) Now: How to propagate a line by bearing & distance

Charles Karney charles.karney at sri.com
Mon Dec 16 18:11:32 PST 2013


geod computes the geodesic (see the man page).  It's probably not a good
idea to call it a great circle (unless you're specifying a sphere).  The
result Frank quotes is from proj 4.9.0 which provides a more accurate
result than earlier versions of proj.4.

   --Charles

On 12/16/2013 08:26 PM, Brent Wood wrote:
> Quick question-
>
> having had a very similar discussion on the Postgis list regarding
> ST_Project()
> does geod propagate a straight line at constant bearing, or a great
> circle arc, to locate the new point?
>
> looking at these examples (the first gives a slightly different return
> azimuth to Frank's example for some reason):
>
> echo "33 -117 15 10000" | geod +ellps=WGS84
> 33d5'12.363"N    116d58'20.1"W    -164d59'5.527"
>
> echo "33 -117 90 10000" | geod +ellps=WGS84
> 32d59'59.835"N 116d53'34.398"W    -89d56'29.986"
>
> echo "33 -117 90 1000000" | geod +ellps=WGS84
> 32d32'41.22"N    106d19'31.529"W    -84d12'35.763"
>
> given the changes in return azimuth & latitude I'm assuming the azimuth
> given on the command line is not the constant bearing to the point, but
> the start bearing of a great circle arc crossing that point.
>
>
> Thanks
>
>    Brent Wood
>
>
> ------------------------------------------------------------------------
> *From:* Frank Warmerdam <warmerdam at pobox.com>
> *To:* PROJ.4 and general Projections Discussions <proj at lists.maptools.org>
> *Sent:* Tuesday, December 17, 2013 11:16 AM
> *Subject:* Re: [Proj] (no subject)
>
> Siva,
>
> I believe you want to use the "geod" program distributed with PROJ.4.
>   The man page, possibly a bit out of date is at:
>
> http://trac.osgeo.org/proj/wiki/man_geod
>
> Example use:
> $ geod +ellps=WGS84
>
> Input:
> 33 -117 15 10000
>
> Output:
> 33d5'13.53"N116d58'20.2"W-164d59'5.581"
>
> So if you gave it an initial location of 33N 117W, an azimuth(bearing)
> of 15 degrees and a distance 10000m you would arrive at the output
> location given with the third value being the azimuth back to the
> original point.
>
> You should be able to construct what you want from this.
>
> There is also a corresponding C API.
>
> Best regards,
> Frank
>
>
>
>
> On Sun, Dec 15, 2013 at 10:56 PM, SIVA RAMA KRISHNA
> <s.r.kriishna at gmail.com <mailto:s.r.kriishna at gmail.com>> wrote:
>
>     Hello,
>
>     Thanks in advance
>
>     I have a *lat, lon* coordinate with a known projection system I
>     assume it a initial point. I want to generate a next Position with
>     the an
>
>     *angle* with  (horizontal/vertical)  and *bearing length*  in
>     lat,lon  coordinates and complete it with angles and bearing lengths
>     to form a
>
>     polygon
>




More information about the Proj mailing list