<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div>Dear all,</div>
<div><br>
</div>
I apologize if this problem is known or if I don't use the command cs2cs correctly. I have experimented a unit conversion problem when transforming coordinates from WGS84 geographic coordinate system to NTF geographic (Paris meridian) coordinate system using
cs2cs.<br>
<br>
Using cs2cs:<br>
------------<br>
version 7.1.0<br>
<div>- $ echo 43.97942110 5.77338693 0.0 | cs2cs EPSG:4326 IGNF:NTFP -f %.8f</div>
<div>Result: 3.43669113 43.97940959 0.00000000<br>
</div>
The right result (IGN code) is: 3.81854570 48.86601066<br>
Here the angle unit conversion is missing. Indeed:<br>
echo 3.43669113 43.97940959 | cct -d 10 +proj=unitconvert +xy_in=deg +xy_out=grad
<br>
Result: 3.8185457000 48.8660106556 inf inf<br>
<br>
Proposed proj string from projinfo:<br>
-----------------------------------<br>
If launched, the proj string proposed by the following command "projinfo -s EPSG:4326 -t IGNF:NTFP" gives the right result in the right unit.<br>
echo 43.97942110 5.77338693 0.0 | cct -d 8 +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=WGS84 +step +proj=helmert +x=168 +y=60 +z=-320 +step +inv +proj=cart +ellps=clrk80ign
+step +proj=pop +v_3 +step +proj=longlat +ellps=clrk80ign +pm=paris +step +proj=unitconvert +xy_in=rad +xy_out=grad<br>
Result: 3.81854570 48.86601066 0.00000000 inf<br>
<div><br>
</div>
<div>I hope that it helps. Do I use the cs2cs command in the right way? <br>
</div>
<div><br>
</div>
<div>Thank you,<br>
</div>
<div>Best regards,</div>
<div><br>
</div>
<div>- Xavier<br>
</div>
<br>
PS: using cs2cs in version 6.3.2, I have also faced coordinate order problem. But I guess it is a known result<br>
The cs2cs command does not recognize the coordinate order:<br>
- $ echo 43.97942110 5.77338693 0.0 | cs2cs +init=EPSG:4326 +to +init=IGNF:NTFP -f %.8f<br>
Result: 41.64152849 5.77099039 0.00000000<br>
This is the wrong result. When putting the longitude first:<br>
- $ echo 5.77338693 43.97942110 0.0 | cs2cs +init=EPSG:4326 +to +init=IGNF:NTFP -f %.8f<br>
I got 3.43669113 43.97940959 0.00000000<br>
The same unit conversion problem as in new version 7.3.1</div>
</body>
</html>