<div>Hello,</div><div>&nbsp;</div><div>&nbsp;</div><div>It is probably wise to give a warning against the new selftest function in the PROJ projection modules.</div><div>&nbsp;</div><div>If the projection code is wrong, the selftest results are in absolute sense completely wrong,</div><div>but the selftest could verify the computer code as perfectly correct.</div><div>&nbsp;</div><div>Example: PJ_mod_ster.c, as found in the repository&nbsp;on <a href="https://github.com/OSGeo/proj.4/tree/master/src">https://github.com/OSGeo/proj.4/tree/master/src</a></div><div>&nbsp;</div><div>What is the problem?</div><div>Please refer to earlier PROJ list messages:</div><div>My message from Feb 17, 2008; 9:31am; subject: Small errors in modified stereographic</div><div>And a&nbsp;repeat from Frank Willett, Sep 10 2011, 10:36:42 EST, subject: Projection gs48 error?<br></div><div>The Lee code has an incorrect coefficient and the GS48 code has an incorrect phi0.</div><div>Assuming Snyder is correct (Map projections a working manual).</div><div>&nbsp;</div><div>For Lee, the first selftest produces:</div><div>-25564478.9526050538, 154490848.8286255</div><div>but it should&nbsp;be:</div><div>-25448425.343461025, 153791681.09640324</div><div>&nbsp;</div><div>For GS48, the first selftest produces:</div><div>-12110635.970867658000, 11668127.145744404000</div><div>but it should be:</div><div>-1923908.4465293458, 355874.65894447843</div><div>&nbsp;</div><div>No, I cannot file a ticket, I cannot produce a patch, I cannot give a working example of command line input and output from a compiled PROJ version.</div><div>&nbsp;</div><div>Please do whatever you want with the above information, but at the moment the PJ_mod_ster.c code is still in error.</div><div>&nbsp;</div><div>&nbsp;</div><div>Yours,</div><div>&nbsp;</div><div>Oscar van Vlijmen</div><div>&nbsp;</div><div>&nbsp;</div><div>&nbsp;</div><div>&nbsp;</div><div>&nbsp;</div>