[Gdal-dev] Transforming shapefiles from euref89 to ed50

[Frank Warmerdam]

Jakob Simon-Gaarde/Trapsoft Inc wrote:
> 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]],

> 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]],
> What am I doing wrong - is it the wkt-files? I got those from two
> shapefiles with the projections I need.

Jacob, 

I believe the problem is that you are indicating a prime meridian
of "Greenwich" at 9 degrees east of greenwich in the euref89 .wkt file.
I suspect the unnamed commerical software is ignoring the prime meridian
while OGR is trying to take it into account. 

Best regards,
-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent