Advertisement
Guest User

Untitled

a guest
Oct 24th, 2014
155
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.49 KB | None | 0 0
  1. pj_svy21 = pj_init_plus("+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")
  2.  
  3. pj_latlong = pj_init_plus("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
  4.  
  5. double test1 = 37621.665;
  6. double test2 = 12056.915;
  7.  
  8. projResult = pj_transform(pj_svy21, pj_latlong, 1, 1, &test1, &test2, NULL );
  9.  
  10. result: latlong
  11. 1.1253115030392806 103.9197679941208889
  12.  
  13. Qgis result:
  14. 1.35650722603 103.690059679 (check on map more accurate)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement