<div dir="ltr">The proj_factors(...) function, which a.o. returns the meridian convergence in the PJ_FACTORS typedef, may be useful in this work.</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Den fre. 25. jun. 2021 kl. 12.26 skrev Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</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>
    <p>Michael,</p>
    <p>A general solution would have more appeal, although I suspect it
      will be a bit more computation intensive as you'll probably have
      to transform to invoke PROJ multiple times per grid point to
      compute the direction to the North</p>
    <p>Even<br>
    </p>
    <div>Le 25/06/2021 à 06:29, Michael Sumner a
      écrit :<br>
    </div>
    <blockquote type="cite">
      <div dir="ltr">Hello, we're using `gdaldem aspect` on polar grid
        products and we need to convert the grid-oriented directions to
        compass oriented ones, location-specific directions, where
        "compass north" is dependent on longitude. 
        <div><br>
        </div>
        <div>It's easy but inefficient for us to do this after running
          gdal dem for our polar grids at latitude -90, and a more
          general solution would require a general approach to "where
          north is" for any location in a projection. </div>
        <div><br>
        </div>
        <div>I'll be able to modify GDAL to do this as an option to
          gdaldem for a simple polar grid (e.g. EPSG:3031), and that
          will make our use of it much more efficient. Is there any
          interest in a more general approach? </div>
        <div><br>
        </div>
        <div>Best, Mike</div>
        <div><br>
        </div>
        <div><br>
          <div><br>
          </div>
          -- <br>
          <div dir="ltr">Michael Sumner<br>
            Software and Database Engineer<br>
            Australian Antarctic Division<br>
            Hobart, Australia<br>
            e-mail: <a href="mailto:mdsumner@gmail.com" target="_blank">mdsumner@gmail.com</a></div>
        </div>
      </div>
      <br>
      <fieldset></fieldset>
      <pre>_______________________________________________
gdal-dev mailing list
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</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>

_______________________________________________<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>
</blockquote></div>