<div dir="ltr">Many thanks, that makes sense, and it solves my problem.<div><br></div><div>---Andrew</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, May 27, 2024 at 11:42 AM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</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"><u></u>

  
    
  
  <div>
    <p>Andrew,</p>
    <p>The error you get is related to the fact that you can't formally
      transform coordinates between 2 vertical CRS, without specifying
      the "horizontal" / geographic CRS into which you enter the long,
      lat.</p>
    <p>Here you, very likely, should combine them with EPSG:4312 "MGI"
      which is the CRS into which the at_bev_GV_Hoehengrid_V1.tif grid
      is registered to</p>
    <p>So:<br>
    </p>
    <p>$ echo 46.642970 12.792060  2513.98 |  PROJ_NETWORK=ON bin/cs2cs
      -d 8 EPSG:4312+9274  EPSG:4312+5778<br>
      46.64297000    12.79206000 2514.15224077<br>
    </p>
    <p>Note that I've specified a latitude, longitude order as EPSG:4312
      is defined with axis in that order</p>
    <p>Even<br>
    </p>
    <div>Le 27/05/2024 à 20:27, Andrew Kirmse
      via PROJ a écrit :<br>
    </div>
    <blockquote type="cite">
      
      <div dir="ltr">I have some Lidar-derived lat/lng/elevation points
        in Austria in <a href="https://epsg.io/9274" target="_blank">EPSG:9274</a> that I want to convert to
        <a href="https://epsg.io/5778" target="_blank">EPSG:5778</a>. 
        I gather this is done by applying the offsets from a particular
        model called the "Höhengrid", which I found at <a href="https://cdn.proj.org/" target="_blank">https://cdn.proj.org/</a>. 
        These points are locations and heights of mountains, and I want
        to convert them to EPSG:5778 because that's what's used on the
        BEV Austrian topographic maps.
        <div><br>
        </div>
        <div>I get this error when I try to do the conversion:</div>
        <div><br>
        </div>
        <div>echo 46.642970 12.792060 2513.98 | cs2cs +init=epsg:9274
          +to +init=epsg:5778<br>
          DeprecationWarning: PROJ_LIB environment variable is
          deprecated, and will be removed in a future release. You are
          encouraged to set PROJ_DATA instead.<br>
          proj_crs_get_geodetic_crs: CRS has no geodetic CRS<br>
          Rel. 9.3.1, December 1st, 2023<br>
          <cs2cs>:<br>
          cannot initialize transformation<br>
          cause: Unknown error (code 4096)<br>
          program abnormally terminated<br>
        </div>
        <div><br>
        </div>
        <div>What does "CRS has no geodetic CRS" mean, and how do I fix
          it?</div>
        <div><br>
        </div>
        <div>(Also, if anyone knows more about why these two different
          EPSGs exist, I am curious.  One is for maps and one is for
          "scientific purposes" but I don't know what that means
          in practice.)</div>
        <div><br>
        </div>
        <div>---Andrew</div>
      </div>
      <br>
      <fieldset></fieldset>
      <pre>_______________________________________________
PROJ mailing list
<a href="mailto:PROJ@lists.osgeo.org" target="_blank">PROJ@lists.osgeo.org</a>
<a href="https://lists.osgeo.org/mailman/listinfo/proj" target="_blank">https://lists.osgeo.org/mailman/listinfo/proj</a>
</pre>
    </blockquote>
    <pre cols="72">-- 
<a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
  </div>

</blockquote></div>