[gdal-dev] projection of GDAL and ArcMap

[agnes123]

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.