<div dir="ltr"><div>Hi Kristian,</div><div><br></div><div>thank you very much for your very clear reply, I was confused by some old posts where it was possible to force cs2cs using a different operation than the best estimate.</div><div><br></div><div>All the best<br></div><div><br></div><div>Stefano<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mer 5 feb 2020 alle ore 14:14 Kristian Evers <<a href="mailto:kreve@sdfe.dk">kreve@sdfe.dk</a>> ha scritto:<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 lang="EN-US">
<div class="gmail-m_1503937248355893685WordSection1">
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)" lang="DA">Stefano,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)" lang="DA"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)">cs2cs uses the ”best guestimate” operation between source and destination CRS.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)">If you want something else you need to specify that yourself. projinfo can help you<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)">find the correct transformation:<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)">$ projinfo -k operation -o PROJ EPSG:1660<u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)">PROJ string:<u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)">+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=push +v_3 +step +proj=cart
 +ellps=intl +step +proj=helmert +x=-104.1 +y=-49.1 +z=-9.9 +rx=0.971 +ry=-2.917 +rz=0.714 +s=-11.68 +convention=position_vector +step +inv +proj=cart +ellps=WGS84 +step +proj=pop +v_3 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)">Which you can then use with cct:<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)">$ echo 41.7447031 15.0618939 | cct -z0 -t0 -d8 +proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad
 +step +proj=push +v_3 +step +proj=cart +ellps=intl +step +proj=helmert +x=-104.1 +y=-49.1 +z=-9.9 +rx=0.971 +ry=-2.917 +rz=0.714 +s=-11.68 +convention=position_vector +step +inv +proj=cart +ellps=WGS84 +step +proj=pop +v_3 +step +proj=unitconvert +xy_in=rad
 +xy_out=deg +step +proj=axisswap +order=2,1<u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)">  41.74535412    15.06180323    0.00000000        0.0000<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)">That should do the trick :-)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)">Please don’t use the +towgs84 parameter – it is deprecated and only supported for legacy reasons.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)">/Kristian<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri",sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><b><span style="font-size:11pt;font-family:"Calibri",sans-serif">From:</span></b><span style="font-size:11pt;font-family:"Calibri",sans-serif"> PROJ <<a href="mailto:proj-bounces@lists.osgeo.org" target="_blank">proj-bounces@lists.osgeo.org</a>>
<b>On Behalf Of </b>Stefano Iacovella<br>
<b>Sent:</b> 5. februar 2020 13:57<br>
<b>To:</b> <a href="mailto:proj@lists.osgeo.org" target="_blank">proj@lists.osgeo.org</a><br>
<b>Subject:</b> [PROJ] cs2cs and +towgs84<u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">Hello,<br>
<br>
I already read the cs2cs usage page at [1] but I am still a little bit confused about the +towgs84 option so I am looking for your help.<br>
<br>
I am trying to convert from EPSG:4265 CRS to EPSG:4326 CRS<br>
<br>
First of all I used the following syntax<br>
<br>
$> cs2cs -d 8 epsg:4265 epsg:4326<br>
41.7447031 15.0618939<br>
41.74535412     15.06180323 0.00000000<br>
<br>
this way the coordinate are transformed but I have no control on which operation is really used, e.g. I don't know which are the towgs84's parameters<br>
<br>
Really I want to use the EPSG:1660 operation, so I tried to add the +towgs84 option but I got an error:<br>
<br>
$> cs2cs -d 8 epsg:4265 epsg:4326 +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68<br>
proj_create: crs not found<br>
Rel. 6.3.0, January 1st, 2020<br>
<cs2cs>:<br>
cannot instantiate source coordinate system<br>
program abnormally terminated<br>
<br>
Which is the correct syntax to be used?<br>
<br>
Thank you in advance<br>
<br>
Stefano<br>
<br>
[1] <a href="https://proj.org/apps/cs2cs.html" target="_blank">https://proj.org/apps/cs2cs.html</a><u></u><u></u></p>
</div>
</div>
</div>

</blockquote></div>