<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
Ah yes, I missed that proj was used in the example.
<div class="">The problem is the same, just located in proj.c instead.</div>
<div class=""><br class="">
</div>
<div class="">I have pushed a fix here: <a href="https://github.com/OSGeo/proj.4/pull/1162" class="">https://github.com/OSGeo/proj.4/pull/1162</a></div>
<div class=""><br class="">
</div>
<div class="">While proj is mainly meant to be used to project geodetic</div>
<div class="">coordinates it doesn’t hurt to treat other input correctly</div>
<div class="">when it is possible. It is a bit closer to doing that now. I</div>
<div class="">am sure there are some exceptions still that I haven’t </div>
<div class="">found (the code in proj.c structured quite poorly)</div>
<div class=""><br class="">
</div>
<div class="">/Kristian<br class="">
<div class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On 26 Oct 2018, at 17:37, Even Rouault <<a href="mailto:even.rouault@spatialys.com" class="">even.rouault@spatialys.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="">On vendredi 26 octobre 2018 14:40:17 CEST Kristian Evers wrote:<br class="">
<blockquote type="cite" class="">If that is the case it should just be a matter of changing the<br class="">
type of the output units internally in the operation. If you<br class="">
submit a new ticket on the github page with a short description<br class="">
I'll take a look at it within the next couple of days.<br class="">
</blockquote>
<br class="">
Things are fine. This is just that using 'proj' for such operation is <br class="">
inapproprate. You should use the new 'cct' binary.<br class="">
<br class="">
'proj' is just for plain-old projection methods that take long, lat as input, <br class="">
whereas here input coordinates are in the projected space<br class="">
<br class="">
echo 1000 0 0 | cct +proj=helmert +convention=coordinate_frame +theta=-162000<br class="">
<br class="">
outputs<br class="">
707.1068 707.1068 0.0000 inf<br class="">
<br class="">
(162000 arc-second = 45 deg)<br class="">
<br class="">
(due to a current oddity of cct, you need to specify 3D coordinates, or add "-<br class="">
z 0" in the parameters )<br class="">
<br class="">
-- <br class="">
Spatialys - Geospatial professional services<br class="">
<a href="http://www.spatialys.com" class="">http://www.spatialys.com</a><br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</body>
</html>