[Gdal-dev] Transforming shapefiles from euref89 to ed50
Jakob Simon-Gaarde/Trapsoft Inc
Jakobs at trapezesoftware.com
Mon Apr 19 09:19:10 EDT 2004
Hello.
I have a shapefile in the projection "UTM euref89 - Zone 32 north" which I
need to transform to the projection "UTM ed50 - Zone 32 north".
I have tried the following:
ogr2ogr -s_srs euref89_32n.wkt -t_srs ed50_32n.wkt out.shp stops_ll.shp
where euref89_32n.wkt contains:
PROJCS["ETRF_1989_UTM_Zone_32N",
GEOGCS["EUREF89_32N",
DATUM["ETRF_1989",
SPHEROID["WGS_1984",6378137.0,298.257223563]],
PRIMEM["Greenwich",9.0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["False_Easting",500000.0],
PARAMETER["False_Northing",0.0],
PARAMETER["Central_Meridian",9.0],
PARAMETER["Scale_Factor",0.9996],
PARAMETER["Latitude_Of_Origin",0.0],
UNIT["Meter",1.0]]
and ed5089_32n.wkt contains:
PROJCS["ED_1950_UTM_Zone_32N",
GEOGCS["GCS_European_1950",
DATUM["European_1950",
SPHEROID["International_1924",6378388,297]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["False_Easting",500000],
PARAMETER["False_Northing",0],
PARAMETER["Central_Meridian",9],
PARAMETER["Scale_Factor",0.9996],
PARAMETER["Latitude_Of_Origin",0],
UNIT["Meter",1]]
ogr results in following transformation:
(523999.60582 ; 6092835.36626) -> (1099082.776 ; 6133192.723)
(527527.56990 ; 6094198.96889) -> (1102420.779 ; 6135007.083)
(550925.36281; 6085171.60262) -> (1126899.377 6129037.630)
I then tried using a commercial transformationtool, and it did the job
right:
(523999.60582 ; 6092835.36626) -> (524080.088 6093041.491)
(527527.56990 ; 6094198.96889) -> (527608.073 ; 6094405.086)
(550925.36281 ; 6085171.60262) -> (551005.934 ; 6085377.520)
What am I doing wrong - is it the wkt-files? I got those from two
shapefiles with the projections I need.
Please help :-(
Med venlig hilsen / Best regards
Jakob Simon-Gaarde
Trapeze Software Europe A/S
Tlf. 87 44 16 00
More information about the Gdal-dev
mailing list