Article View: pl.comp.gis
Article #12344Re: QGIS i błąd reprojekcji
From: mt
Date: Fri, 12 Nov 2010 12:07
Date: Fri, 12 Nov 2010 12:07
35 lines
2717 bytes
2717 bytes
mt pisze: > Mateusz Loskot pisze: > >> Ja się nie spotkałem, ale nie wykluczone, że jest błąd. >> Proponuję wrzucić ticket do bazy błędów: >> >> http://trac.osgeo.org/qgis > > Dzięki, zgłoszone. > Aczkolwiek doczytałem jeszcze na forum QGIS, że nie tylko ja mam > problemy z reprojekcją rastrów z układu 1965 do 1992 czy 2000. Pewnym > rozwiązaniem jest reprojekcja pośrednia przez np. WGS84, wtedy błąd > przeliczenia wynosi około 5 - 6 metrów w obu osiach, ale już nie ponad > 100 :). Inną sprawą jest co gdalwarp robi z tym plikiem pośrednim w > WGS84. Co prawda współrzędne zaszyte w geotiffie są mniej więcej ok, ale > raster jest jakby mocno rozciągnięty poziomo. Podejrzałem nagłówek tiffa > i z pliku o wymiarach ok. 12000x8000 pikseli po reprojekcji do WGS84 > powstaje plik o wymiarach ok. 14000x5000 pikseli, stąd ten efekt > rozciągnięcia/spłaszczenia, a to chyba też nie powinno mieć miejsca? To jeszcze odpowiem poniekąd sobie. Problem zdaje się jest w nowych wydaniach bibliotek GDAL. W starszej wersji FWTools (ja mam 1.3.5) reprojekcja z 1965/3 do 1992 przebiega bezbłędnie. Używając QGIS 1.5.0 czy FWTools w wydaniu 2.4.7 ten sam raster jest przesunięty o ponad 130 metrów! Nie śledzę dokładnie historii implementacji różnych wersji GDAL, więc nie wiem kiedy jakiś ewentualny błąd mógł się pojawić. Jeszcze mam pytanie, da się, mając zainstalowanego QGIS i FWTools w jakiejś starszej wersji, podmienić pliki i które ewentualnie, aby wtyczka do QGIS korzystała ze starszych wersji bibliotek GDAL? Sama zamiana pliku gdalwarp.exe nie daje porządanego rezultatu i raster po reprojekcji wciąż jest przesunięty. -- pozdrawiam marcin
Technical Details
Message-ID:
<ibjami$so5$1@usenet.news.interia.pl>Path:
mbox-mo1.pugleaf.net!news-archive.icm.edu.pl!news.gazeta.pl!newsfeed.pionier.net.pl!news.nask.pl!news.nask.org.pl!news.interia.pl!not-for-mailReferences:
<ibhm6d$8nb$1@usenet.news.interia.pl> <ibhoib$82d$1@inews.gazeta.pl> <ibi1q2$rb0$1@usenet.news.interia.pl>Show Full Headers
Complete Article Headers
Path: mbox-mo1.pugleaf.net!news-archive.icm.edu.pl!news.gazeta.pl!newsfeed.pionier.net.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: Fri, 12 Nov 2010 13:07:44 +0100 Organization: INTERIA.PL S.A. Lines: 35 Message-ID: <ibjami$so5$1@usenet.news.interia.pl> References: <ibhm6d$8nb$1@usenet.news.interia.pl> <ibhoib$82d$1@inews.gazeta.pl> <ibi1q2$rb0$1@usenet.news.interia.pl> NNTP-Posting-Host: host-94-231-62-53.ip.jarsat.pl Mime-Version: 1.0 X-Trace: usenet.news.interia.pl 1289563666 29445 94.231.62.53 (12 Nov 2010 12:07:46 GMT) X-Complaints-To: usenet@firma.interia.pl NNTP-Posting-Date: Fri, 12 Nov 2010 12:07:46 +0000 (UTC) User-Agent: Thunderbird 2.0.0.24 (Windows/20100228) In-Reply-To: <ibi1q2$rb0$1@usenet.news.interia.pl> Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit