🚀 go-pugleaf

W dniu 11.11.2010 22:11, mt pisze:

> Za pomocą wtyczki GDAL Tools do QGISa (wersja 1.2.0 opartego na GDAL
> w wersji 1.6.3) dokonałem hurtowej reprojekcji rastrów z 1965/3 na
> 1992 no i po otwarciu z pozostałymi danymi wektorowymi zdziwienie,
> gdyż nastąpił spory rozjazd między wektorami a rastrami w nowym
> odwzorowaniu, zacząłem sprawdzać dokładniej i się okazało, że po
> reprojekcji błąd przeliczenia współrzędnych wynosi około 130 metrów w
> osi równoleżnikowej i około 30 metrów w osi południkowej. Spotkał sie
> ktoś z Was już z czymś takim?

Przyczyna błedu jest (była) w GDALu - w 1.7.x wcięło definicję towgs84
dla polskich układów opartych na elipsoidzie Krasowskiego. Wcześniej,
jeszcze przed 1.6.3, GDAL przestał stosować towgs84=0,0,0 dla układów
opartych na GRS80. Skutki uboczne koniecznych zmian w realizacji układów
współrzędnych w GDALu.

Popraw (albo stwórz, jeśli nie istnieje) plik
sciezka_do_gdala/share/gdal/gcs.override.csv zgodnie z
http://trac.osgeo.org/gdal/changeset/19810. Daj znać czy to załatwia sprawę.

W starszych QGISach warto też przejrzeć czy wpisy dla nastepujących
EPSG: 3120 2172 2173 2174 2175 3333 3334 3335 3329 3330 3331 3332 3328
4179 mają jawne, poprawne towgs84. Jeśli nie, należy zdefiniować własne
parametry dla tych układów - jako kopie istniejących w QGISie plus
odpowiednia definicja towgs84, poprzez Opcje > Układ współrzędnych
użytkownika.

Niedawno wydany GDAL 1.7.3 zawiera już poprawkę. Oparte na nim
oprogramowanie (miejmy nadzieję, że równiez nadchodzący QGIS 1.6)
powinny być pozbawione problemu raz na zawsze.

BTW: w GRASSie, który używa własnej kopii gdalowych *.csv z definicjami
układów, poprawkę naniosłem jeszcze przed wydaniem 6.4.0
(http://trac.osgeo.org/grass/changeset/42666).

Maciek

-- 
Maciej Sieczka
http://www.sieczka.org