<div dir="auto">This may help getting the factors: <a href="https://pyproj4.github.io/pyproj/stable/api/proj.html#pyproj.proj.Proj.get_factors">https://pyproj4.github.io/pyproj/stable/api/proj.html#pyproj.proj.Proj.get_factors</a></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, May 15, 2020, 5:33 PM Charles Karney <<a href="mailto:charles@karney.com" rel="noreferrer noreferrer" target="_blank">charles@karney.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I believe PROJ does have some way of computing the scale of a<br>
projection; but it uses numerical differencing.  GeographicLib computes<br>
the scale of TM from an analytic expression, so it's probably more<br>
accurate.  The command I used to get the scale (using GeographicLib) is:<br>
<br>
     echo 0.5 |<br>
     GeodSolve -I 31.4617 -83.4849 31.4704 -83.4781 -F -p 10 |<br>
     cut -f1,2 -d' ' | GeoConvert -c -p 10<br>
=><br>
     -1.2959245331415 1.000285683174366<br>
<br>
The GeodSolve expression gets the midpoint of a geodesic.  GeoConvert -c<br>
returns the meridian convergence and scale for the UTM projection.<br>
<br>
The command line utility geod uses the same code as pyproj.Geod.  So<br>
unless you have extra versions of geod lying around they should return<br>
the same results.<br>
<br>
<br>
On 5/15/20 5:00 PM, Dima Kogan wrote:<br>
> Thank you very much. That's VERY useful info. How did you get that scale<br>
> factor? Is there a proj tool and/or function call?<br>
> <br>
> I have just tried to run the "geod" cmd on another box with another proj<br>
> install, and it matches up over there. Which is probably enough for me<br>
> to go figure this out. We DO expect the "geod" cmdline and the<br>
> pyproj.Geod results to match 100%, right?<br>
> <br>
> Thanks again.<br>
> _______________________________________________<br>
> PROJ mailing list<br>
> <a href="mailto:PROJ@lists.osgeo.org" rel="noreferrer noreferrer noreferrer" target="_blank">PROJ@lists.osgeo.org</a><br>
> <a href="https://lists.osgeo.org/mailman/listinfo/proj" rel="noreferrer noreferrer noreferrer noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/proj</a><br>
> <br>
_______________________________________________<br>
PROJ mailing list<br>
<a href="mailto:PROJ@lists.osgeo.org" rel="noreferrer noreferrer noreferrer" target="_blank">PROJ@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/proj" rel="noreferrer noreferrer noreferrer noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/proj</a><br>
</blockquote></div>