I am using a shape of Austria from gadm.org which I want to overlay with a netCDF file in R. As the netCDF has projected coordinates in EPSG:31287, I have reprojected the shape using FWtools and the .proj4-string from spatialreference.org.
Unfortunately the shape does not fit the raster data: The best fit is in the west, and the error grows in eastern direction (shape has a bigger extent in eastern direction than the raster). The north/south extent is correct.
Any ideas what could be wrong?
I already tried to play around with
lon_0 values, but as far I am informed correctly, this causes only a general shift.