<!DOCTYPE html>
<html>
  <head>
    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
  </head>
  <body>
    <br>
    <blockquote type="cite"
cite="mid:AS5PR09MB68978280BB1A4FABD412B666FDE62@AS5PR09MB6897.eurprd09.prod.outlook.com">
      <div class="WordSection1">
        <p class="MsoNormal"><span style="mso-fareast-language:EN-US"
            lang="EN-US">
            Now the same with ogrinfo. I simplified the SQL by  saving
            the test linesting into a file in EPSG:3067.<br>
            <br>
            ogrinfo -dialect sqlite -sql "select ST_Length(geometry)
            from length" length.json<br>
            208808.794960845<br>
            <br>
            <o:p></o:p></span></p>
        <p class="MsoNormal"><span style="mso-fareast-language:EN-US"
            lang="EN-US">Good so far. But with “use_ellipsoid” gives
            very unexpected result:<o:p></o:p></span></p>
        <p class="MsoNormal"><span style="mso-fareast-language:EN-US"
            lang="EN-US"><br>
            ogrinfo -dialect sqlite -sql "select ST_Length(geometry,1)
            from length" length.json<br>
            <br>
            417750.100376289<br>
            <br>
            <o:p></o:p></span></p>
        <p class="MsoNormal"><span style="mso-fareast-language:EN-US"
            lang="EN-US">Same wrong result also with this:<br>
            <br>
            ogrinfo -dialect sqlite -sql "select
            ST_Length(ST_Transform(geometry,4326),1) from length"
            length.json<o:p></o:p></span></p>
        <p class="MsoNormal"><span style="mso-fareast-language:EN-US"
            lang="EN-US"><o:p> </o:p></span></p>
        <p class="MsoNormal"><span style="mso-fareast-language:EN-US"
            lang="EN-US">length.json is this:<br>
          </span></p>
      </div>
    </blockquote>
    <p>Jukka,</p>
    <p>Thanks for the feedback. ST_Length(geometry, use_ellipsoid=1) is
      indeed seriously broken in GDAL 3.10 as it only properly works on
      closed linestrings / polygons</p>
    <p>Fix in <a class="moz-txt-link-freetext" href="https://github.com/OSGeo/gdal/pull/11700">https://github.com/OSGeo/gdal/pull/11700</a><br>
    </p>
    <p>With it,</p>
    <p>$ ogrinfo -dialect sqlite -sql "select ST_Length(geometry, 1)
      from test" test.json -al -q<br>
      <br>
      Layer name: SELECT<br>
      OGRFeature(SELECT):0<br>
        ST_Length(geometry, 1) (Real) = 208875.050188144</p>
    <p>Even</p>
    <br>
    <pre class="moz-signature" cols="72">-- 
<a class="moz-txt-link-freetext" href="http://www.spatialys.com">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
  </body>
</html>