<DIV>Hi all,</DIV>
<DIV>I have one trouble about projection. My testing area is around UTM zone-47 and zone-48. The official road shapefile is translated both zones into zone-47 (Indian 1975). My road is in zone-48 (WGS84). So, I'd like to translate my utm-48 road data into zone-47. I used:</DIV>
<DIV> </DIV>
<DIV>$proj48=ms_newProjectionObj("proj=utm,zone=48,ellps=WGS84,datum=WGS84,units=m,no_defs");</DIV>
<DIV>$proj47=ms_newProjectionObj("proj=utm,zone=47,a=6377276.345,b=6356075.413,towgs84=209,818,290,0,0,0,0,units=m,no_defs");</DIV>
<DIV> </DIV>
<DIV>$line=ms_newLineObj();</DIV>
<DIV>$pt=ms_newPointObj();</DIV>
<DIV>$pt->setXY(195397,1604329);</DIV>
<DIV>$pt->project($proj48,$proj47);</DIV>
<DIV>$line->add($pt);</DIV>
<DIV>
<DIV>$pt->setXY(195385,1604354);</DIV>
<DIV>$pt->project($proj48,$proj47);</DIV>
<DIV>$line->add($pt);</DIV>
<DIV>.............</DIV>
<DIV>$shp->add($line);</DIV>
<DIV> </DIV>
<DIV>My road should overlay onto official road. But it wasn't on the right place, it shifted to the left of official road about 350m. If I use data by not translating in map file as:</DIV>
<DIV> </DIV>
<DIV>.............</DIV>
<DIV># ------------- Official road ----------</DIV>
<DIV>LAYER</DIV>
<DIV> NAME "road47"</DIV>
<DIV> DATA "road/road47</DIV>
<DIV> TYPE line</DIV>
<DIV> PROJECTION</DIV>
<DIV> "proj=utm"</DIV>
<DIV> "zone=47"</DIV>
<DIV> "a=6377276.345"</DIV>
<DIV> "b=6356075.413"</DIV>
<DIV> "towgs84=209,818,290,0,0,0,0"</DIV>
<DIV> "units=m"</DIV>
<DIV> "no_defs"</DIV>
<DIV> END</DIV>
<DIV> CLASS</DIV>
<DIV> COLOR 0 200 0</DIV>
<DIV> END</DIV>
<DIV>END</DIV>
<DIV>
<DIV># ------------- my road ----------</DIV>
<DIV>LAYER</DIV>
<DIV> NAME "road48"</DIV>
<DIV> DATA "road/road48</DIV>
<DIV> TYPE line</DIV>
<DIV> PROJECTION</DIV>
<DIV>
<DIV> "proj=utm"</DIV>
<DIV> "zone=48"</DIV>
<DIV> "ellps=WGS84"</DIV>
<DIV> "datum=WGS84"</DIV>
<DIV> "units=m"</DIV>
<DIV> "no_defs"</DIV> END</DIV>
<DIV> CLASS</DIV>
<DIV> COLOR 0 200 0</DIV>
<DIV> END</DIV>
<DIV>END</DIV>
<DIV>.......</DIV>
<DIV> </DIV>
<DIV>It works as it should be. What I've done wrong about the script above. Could I use this way to re-projection on the fly?</DIV>
<DIV> </DIV>
<DIV>Thank you.</DIV>
<DIV> </DIV>
<DIV>Krung</DIV>
<DIV> </DIV></DIV></DIV><p>
<hr size=1> <a href="http://us.rd.yahoo.com/evt=34442/*http://www.yahoo.com/r/hs">Start your day with Yahoo! - make it your home page </a>