🚀 go-pugleaf

RetroBBS NetNews Server

Inspired by RockSolid Light RIP Retro Guy

Article View: pl.comp.gis
Article #12357

Re: QGIS i błąd reprojekcji

#12357
From: mt
Date: Thu, 18 Nov 2010 00:45
36 lines
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-mail
References: <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