<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body dir="auto"><div dir="ltr"></div><div dir="ltr">Thanks every one for your help,</div><div dir="ltr"><br></div><div dir="ltr">I found good packages out there. I will probably go with Proj. I’am just surprised it’s not part standard libraries like Shapely or OGR since both can calculate line length. Most packages look like they are made to calculate only two point at a time and not, for example, a WKT string with multiple nodes. I guess calculating distance between each vertex would be the way to go.</div><div dir="ltr"><br></div><div dir="ltr">Cheer, </div><div dir="ltr">Nicolas</div><div dir="ltr"><br>Le 15 juin 2019 à 07:02, Thomas Knudsen <<a href="mailto:knudsen.thomas@gmail.com">knudsen.thomas@gmail.com</a>> a écrit :<br><br></div><blockquote type="cite"><div dir="ltr"><meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="ltr"><div>In the PROJ package, the geod utility is included - you could do</div><div><br></div><div>echo lat0 lon0 lat1 lon1 | geod -I +ellps=GRS80</div><div><br></div><div>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.</div><div><br></div><div>e.g.</div><div><br></div><div>$ echo 55 12 56 12 |geod -I +ellps=GRS80</div><div>> 0d -180d 111332.699 <br></div><div><br></div><div>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. <br></div><div></div><div>Otherwise, you could use the "proj_lp_dist()" entry of the PROJ API - cf. <a href="https://proj.org/development/reference/functions.html">https://proj.org/development/reference/functions.html</a></div><div><br></div><div>(and probably you should continue this thread on the PROJ mailing list if you need additional assistance)<br></div><div><br></div><div>/Thomas<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Den lør. 15. jun. 2019 kl. 06.43 skrev Nicolas Cadieux <<a href="mailto:nicolas.cadieux@archeotec.ca">nicolas.cadieux@archeotec.ca</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="auto"><div dir="ltr"></div><div dir="ltr">Thanks,</div><div dir="ltr">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!</div><div dir="ltr">Cheers </div><div dir="ltr">Nicolas </div><div dir="ltr"><br>Le 14 juin 2019 à 13:20, Patrick Young <<a href="mailto:patrick.mckendree.young@gmail.com" target="_blank">patrick.mckendree.young@gmail.com</a>> a écrit :<br><br></div><blockquote type="cite"><div dir="ltr"><div dir="ltr">Not exactly what you want, but you can do this with PostGIS by casting your geometry to the geography type:<div><br></div><div><a href="https://postgis.net/workshops/postgis-intro/geography.html" target="_blank">https://postgis.net/workshops/postgis-intro/geography.html</a></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jun 14, 2019 at 10:26 AM Atle Frenvik Sveen <<a href="mailto:atle@frenviksveen.net" target="_blank">atle@frenviksveen.net</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi! <br>
<br>
Not sure if would want to use gdal for this task*, but take a look at this blog post: <br>
<br>
<a href="https://janakiev.com/blog/gps-points-distance-python/" rel="noreferrer" target="_blank">https://janakiev.com/blog/gps-points-distance-python/</a><br>
<br>
<br>
*or if it's doable, i guess not, since the scope of gdal is reading/writing geospatial formats<br>
<br>
-a<br>
<br>
- <br>
Atle Frenvik Sveen<br>
<a href="mailto:atle@frenviksveen.net" target="_blank">atle@frenviksveen.net</a><br>
45278689<br>
<a href="http://atlefren.net" rel="noreferrer" target="_blank">atlefren.net</a><br>
<br>
On Fri, Jun 14, 2019, at 17:32, Nicolas Cadieux wrote:<br>
> Hi,<br>
> <br>
> I am trying to get the length of a line in python. (Not just the <br>
> straight length between the first and last nodes). Using geopandas, <br>
> (therefore the Shapely lib) I am getting the euclidien distance even <br>
> though the dataframe holdings the line geometries has a CRS (WGS84, <br>
> zone UTM 18 S). Obviously, the WGS84 Ellipsoid is not taken into <br>
> account. <br>
> <br>
> Can I do this with gdal/ogr? <br>
> <br>
> Thanks for the help<br>
> Nicolas<br>
> _______________________________________________<br>
> gdal-dev mailing list<br>
> <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
> <a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a></blockquote></div>
</div></blockquote></div>_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a></blockquote></div>
</div></blockquote></body></html>