Hi!<br>
<br>
I have been trying to reproject a point from one projection to an other
with phpmapscript and its projection object (with the PROJ.4 that is
enhanced with phpmapscript). I am using phpmapscript_48.dll (i supose
it refers to version 4.8...). <br>
<br>
The original projection of the point that I am trying to reproject is a
greek projection system called "TM3". It has the following attributes:<br>
<ul>
<li>projection = <i><font color="navy" face="Arial" size="2"><b>T</b>ranvserse <b>M</b>ercator 3<sup>o</sup></font></i></li><li>ellpisoid = bessel</li><li>lon0 = <font color="navy" face="Arial" size="2">23<sup>ï</sup>42' 58''.815
</font></li><li><font color="navy" face="Arial" size="2">lat0 = 34</font></li><li><font color="navy" face="Arial" size="2">towgs84 = DX=456.39, DY=372.62, DZ=496,82, DangleX=0, DangleY=0, DangleZ=0.000363, ppm=0.0000159</font>
</li><li><font color="navy" face="Arial" size="2">False E = 200000</font></li><li><font color="navy" face="Arial" size="2">False N = 0</font></li><li><font color="navy" face="Arial" size="2">Ellipse big axis = </font><font color="navy" face="Arial" size="2">
6377397.155m</font></li><li><font color="navy" face="Arial" size="2">1/f = </font><font color="navy" face="Arial" size="2">
1/299.1528128</font></li><li><font color="navy" face="Arial" size="2">k0 = </font><font color="navy" face="Arial" size="2">
0.9999</font></li><li><font color="navy" face="Arial" size="2">units = m</font></li>
</ul>
The second projection is the universal WGS84 (LAT/LONG)<br>
<br>
The point's initial values are:<br>
X = 134216<br>
Y = 731046<br>
The correct tranformation should have give us: <br>
lon=22.562799<br>
lat=40.3511194<br>
<br>
Though I keep taking something like:<br>
lon=22.9206930<br>
lat=40.5838561<br>
The diference seems to be small with naked eye, though it hundreds of meters...<br>
<br>
I have tryed (believe me) with a number of different ways to export the
right lat/long though the result came again and again the same. I
provide the code that I used:<br>
[CODE]<br>
$projInObj = ms_newprojectionobj("proj
<div id="mb_0">=tmerc,ellps=bessel,towgs84=456.39,372.62,496,82,0,0d00m00363s,0.0000159,".<br>
"lat_0=34d0m0s,lon_0=23d42m58815s,x_0=200000,units=m,k=0.9999,a=6378137,rf=298.25722101");<br>
<br>
$projOutObj=ms_newprojectionobj("proj=latlong,ellps=WGS84,DATUM=WGS84");<br>
<br>
$myshape = ms_newpointobj();<br>
$myshape->setXY(134216, 731046); <br>
<br>
$myshape->project($projOutObj, $projInObj);<br>
<br>
print_r($myshape);<br>
[/CODE]<br>
<br>
Can you please help me understand why I keep getting bad results? Is it
something that I am doing wrong? Does it have to do with a PROJ problem?<br>
<br>
Thanks for the time</div>