Article View: pl.comp.gis
Article #12357Re: QGIS i błąd reprojekcji
From: mt
Date: Thu, 18 Nov 2010 00:45
Date: Thu, 18 Nov 2010 00:45
36 lines
2600 bytes
2600 bytes
Maciej Sieczka pisze: > 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ę. No niestety nie załatwia, chyba że ja coś źle robię, ale o tym poniżej. Generalnie plik gcs.override.csv występuje w QGIS w trzech miejscach: C:\QGIS\apps\gdal-16\share\gdal\ C:\QGIS\share\epsg_csv\ C:\QGIS\share\gdal\ W pierwszym i trzecim przypadku dodałem do wspomnianego pliku: 4284,Pulkovo 1942,6284,Pulkovo 1942,6284,9122,7024,8901,1,0,6422,9607,23.9,-141.3,-80.9,0,-0.37,-0.85,-0.12 W drugim przypadku ten zapis był już standardowo w pliku gcs.override jednak po wszystkich zmianach qgis nadal podaje towgs84=33.4,-146.6,-76.3,-0.359,-0.053,0.844,-0.84 +units=m +no_defs Jak domniemam to jest przyczyna błędnego przeliczania? Ponadto, w panelu Układ współrzędnych użytkownika zdefiniowałem na nowo 65/3 z poprawionym parametrem towgs84, pełne parametry poniżej: +proj=sterea +lat_0=53.58333333333334 +lon_0=17.00833333333333 +k=0.9998 +x_0=3501000 +y_0=5999000 +ellps=krass +towgs84=23.9,-141.3,-80.9,0,-0.37,-0.85,-0.12 +units=m +no_defs No i dalej jest ten sam ponad 100-metrowy błąd po reprojekcji, więc sam już nie wiem gdzie leży przyczyna. Ale dzięki za poświęcony czas! -- pozdrawiam marcin
Technical Details
Message-ID:
<ic1svn$5rp$1@usenet.news.interia.pl>Path:
mbox-mo1.pugleaf.net!news-archive.icm.edu.pl!news.icm.edu.pl!news.nask.pl!news.nask.org.pl!news.interia.pl!not-for-mailReferences:
<ibhm6d$8nb$1@usenet.news.interia.pl> <4CE43D09.6060208@sieczka.org>Show Full Headers
Complete Article Headers
Path: mbox-mo1.pugleaf.net!news-archive.icm.edu.pl!news.icm.edu.pl!news.nask.pl!news.nask.org.pl!news.interia.pl!not-for-mail From: mt <telepski@interia.pl> Newsgroups: pl.comp.gis Subject: Re: QGIS i błąd reprojekcji Date: Thu, 18 Nov 2010 01:45:45 +0100 Organization: INTERIA.PL S.A. Lines: 37 Message-ID: <ic1svn$5rp$1@usenet.news.interia.pl> References: <ibhm6d$8nb$1@usenet.news.interia.pl> <4CE43D09.6060208@sieczka.org> NNTP-Posting-Host: host-94-231-62-53.ip.jarsat.pl Mime-Version: 1.0 X-Trace: usenet.news.interia.pl 1290041143 6009 94.231.62.53 (18 Nov 2010 00:45:43 GMT) X-Complaints-To: usenet@firma.interia.pl NNTP-Posting-Date: Thu, 18 Nov 2010 00:45:43 +0000 (UTC) User-Agent: Thunderbird 2.0.0.24 (Windows/20100228) In-Reply-To: <4CE43D09.6060208@sieczka.org> Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit