[gdal-dev] Problem with Shapefile and CRS 2056/21781

[Stefan Keller]

Given a Shapefile "swiss_admin_boundaries.shp" in EPSG:2056 (or 21781,
to be downloaded for free at [1]) I have problems transforming it to
any other SRS/CRS.

When doing this:
> ogr2ogr -f "GPKG" swiss_admin_boundaries1.gpkg swiss_admin_boundaries.shp -t_srs EPSG:3857

... swiss_admin_boundaries2.gpkg have a shift of about 150m.

But when I convert Shapefile first e.g. to GeoPackage, and transform then:
> ogr2ogr -f "GPKG" tmp.gpkg swiss_admin_boundaries.shp
> ogr2ogr -f "GPKG" swiss_admin_boundaries2.gpkg tmp.gpkg -t_srs EPSG:3857
> del tmp.gpkg

... swiss_admin_boundaries2.gpkg is correct!

What could be the reason?

:Stefan

[1] https://opendata.swiss/en/dataset/swissboundaries3d-gemeindegrenzen