<div dir="ltr"><div><div><span style="font-family:arial,helvetica,sans-serif"><font>Dear all<br><br></font></span></div><span style="font-family:arial,helvetica,sans-serif"><font>I have got a small question that I can not find the confirmed answer on the web. I am now using pj_transform to transform coordinates from svy21 to latlon using the following code snippet:<br><br>pj_svy21 = pj_init_plus(&quot;+proj=tmerc +lat_0=1.366666666666667 +lon_0=103.8333333333333 +k=1 +x_0=28001.642 +y_0=38744.572 +ellps=WGS84 +units=m +no_defs&quot;)<br>pj_latlong = pj_init_plus(&quot;+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs&quot;)<br><br>double svy_east = 12056.915;<br>double svy_north = 37621.665;<br>    <br>double lat_radian = svy_east;<br>double lon_radian = svy_north;<br>projResult = pj_transform(pj_svy21, pj_latlong, 1, 0, &amp;lat_radian, &amp;lon_radian, NULL );<br><br></font></span></div><span style="font-family:arial,helvetica,sans-serif"><font>But the result (lat_radian, lon_radian) I got is reversed, i.e. the result given by <br></font></span><pre><span style="font-family:arial,helvetica,sans-serif"><font>pj_transform<span class="">(</span>pj_svy21, pj_latlong<span class="">,</span> <span class="">1</span><span class="">,</span> 0<span class=""></span><span class="">,</span> <span class="">&amp;</span>x<span class="">,</span> <span class="">&amp;</span>y<span class="">,</span> <span class="">NULL</span> <span class="">)<br><br></span></font></span></pre><pre><span style="font-family:arial,helvetica,sans-serif"><font><span class="">is x corresponds to longitude and y corresponds to latitude. It is not the other way around as I did.<br><br></span></font></span></pre><pre><span style="font-family:arial,helvetica,sans-serif"><font><span class="">Is this the way it worked or I just used this function wrongly? Thanks.<br><br></span></font></span></pre><pre><span style="font-family:arial,helvetica,sans-serif"><font><span class="">Best regards,</span></font></span><span class=""><br>Zhang Qun<br></span></pre></div>