Dear All, I'm trying to project a shapefile from Transverse Mercator projection with Krasovsky 1940 datum to Lambert Conformal Conic with D_Sphere_ARC_INFO datum by ArcMap and GDAL.
I project the shapefile by importing it to the ArcMap which is with the target projection in advance. There is a geographic transformation warning box comes out and i choose (1) close the transformation box (2)add a new transformation parameters of geocentric method with 0 x,y,z axis translation. When compare to the result of GDAL with the command: target_prj='+proj=lcc +a=6370997 +b=6370997 +towgs84=0,0,0,0,0,0,0 +lat_0=28.5 +lat_1=15 +lat_2=40 +lon_0=114 +x_0=0 +y_0=0 +units=m ' org_prj=' +proj=tmerc +ellps=krass +lat_0=0 +lon_0=111 +k=1 +x_0=500000 +y_0=0.0 +units=m +no_defs ' ogr2ogr -skipfailures -s_srs org_prj -f ESRI Shapefile -t_srs target_prj a.shp b.shp I found out that the result is idential to (2) one. But what i want is to generate the result of (1) by using GDAL. The result of (2) and (1) has vertical difference around 20km. What should i need to do in order to get the result of (1) by GDAL?Is there something wrong for my projection definition? It would be grateful if someone can give me suggestion on how to do this! Regards, Agnes -- View this message in context: http://n2.nabble.com/projection-of-GDAL-and-ArcMap-tp4721380p4721380.html Sent from the GDAL - Dev mailing list archive at Nabble.com. _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev