jak dodać warstwy z geoportalu

[borys]

To wszystko jasne. Problem leży w bazie układów współrzędnych i bibliotece PROJ.4, dokonującej reprojekcji. Ten poprawiony układ 92 (2180) różni się od standardowego ("GDAL compatible", 42180) dodatkowym parametrem +towgs84, opisującym elipsoidę, a konkretnie parametry jej transformacji do elipsoidy wgs84. Dla elipsoidy etrs89 w uproszczeniu przyjmuje się 0,0,0. Jego brak uniemożliwia przeniesienie współrzędnych pomiędzy elipsoidami i stąd przesunięcie podczas transformacji między układami opartymi na różnych elipsoidach. Jako docelowy układ transformacji wystarczy więc wybrać ten poprawiony 92, a nie kompatybilny z gdalem. Gorzej jest ze źródłowym, bo niepolskich układów nie poprawiałem i stąd 31469 nie posiada parametru +towgs84. Wyjściem jest podanie go podczas transformacji nie w postaci kodu epsg, tylko całej tej deklaracji '+proj=tmerc +lat_0=0 +lon_0=15 +k=1 +x_0=5500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs', wzbogaconej o parametr +towgs84.

Niestety nie umiem Ci powiedzieć, jak powinien wyglądać dla elipsoidy Bessela. Pogooglaj pod hasłem: +towgs84 31469 i popróbuj wartości, które ludzie podają (bo widzę, że podają różne).