In the PROJ package, the geod utility is included - you could do echo lat0 lon0 lat1 lon1 | geod -I +ellps=GRS80
Where the "-I" indicates that you want to calculate the "inverse geodetic problem", i.e. you know where you are, and where you want to go, but you need to know how far to go and in which direction. e.g. $ echo 55 12 56 12 |geod -I +ellps=GRS80 > 0d -180d 111332.699 Where the distance (in m) you look for is the third number on the output line. The two first are the forward and return azimuths respectively. Otherwise, you could use the "proj_lp_dist()" entry of the PROJ API - cf. https://proj.org/development/reference/functions.html (and probably you should continue this thread on the PROJ mailing list if you need additional assistance) /Thomas Den lør. 15. jun. 2019 kl. 06.43 skrev Nicolas Cadieux < [email protected]>: > Thanks, > Could work but I think this will be too slow. I wonder how QGIS does it? I > guess they use code from Proj.4. If anyone has an other idea, shoot! > Cheers > Nicolas > > Le 14 juin 2019 à 13:20, Patrick Young <[email protected]> > a écrit : > > Not exactly what you want, but you can do this with PostGIS by casting > your geometry to the geography type: > > https://postgis.net/workshops/postgis-intro/geography.html > > On Fri, Jun 14, 2019 at 10:26 AM Atle Frenvik Sveen <[email protected]> > wrote: > >> Hi! >> >> Not sure if would want to use gdal for this task*, but take a look at >> this blog post: >> >> https://janakiev.com/blog/gps-points-distance-python/ >> >> >> *or if it's doable, i guess not, since the scope of gdal is >> reading/writing geospatial formats >> >> -a >> >> - >> Atle Frenvik Sveen >> [email protected] >> 45278689 >> atlefren.net >> >> On Fri, Jun 14, 2019, at 17:32, Nicolas Cadieux wrote: >> > Hi, >> > >> > I am trying to get the length of a line in python. (Not just the >> > straight length between the first and last nodes). Using geopandas, >> > (therefore the Shapely lib) I am getting the euclidien distance even >> > though the dataframe holdings the line geometries has a CRS (WGS84, >> > zone UTM 18 S). Obviously, the WGS84 Ellipsoid is not taken into >> > account. >> > >> > Can I do this with gdal/ogr? >> > >> > Thanks for the help >> > Nicolas >> > _______________________________________________ >> > gdal-dev mailing list >> > [email protected] >> > https://lists.osgeo.org/mailman/listinfo/gdal-dev >> _______________________________________________ >> gdal-dev mailing list >> [email protected] >> https://lists.osgeo.org/mailman/listinfo/gdal-dev > > _______________________________________________ > gdal-dev mailing list > [email protected] > https://lists.osgeo.org/mailman/listinfo/gdal-dev
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
