Article View: pl.comp.gis
Article #12345Re: QGIS i błąd reprojekcji
From: Mateusz Loskot
Date: Fri, 12 Nov 2010 12:16
Date: Fri, 12 Nov 2010 12:16
55 lines
3444 bytes
3444 bytes
"mt" <telepski@interia.pl> wrote in message news:ibjami$so5$1@usenet.news.interia.pl... > 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. GDAL używa PROJ.4 - http://trac.osgeo.org/proj Jeśli błąd leży w samej transformacji, to jest to problem PROJ.4. Możesz spróbować wyizolować przykładowe punkty z rastra i sprawdzić ich transformację za pomocą narzędzi proj: http://trac.osgeo.org/proj/wiki/man_proj http://trac.osgeo.org/proj/wiki/man_cs2cs > 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. Teoretycznie jest to możliwe, ale jedynie w obrębie kompatybilnych wersji GDAL. Na przykład, jeśli QGIS używa GDAL 1.7.0, a w FWTools masz GDAL 1.6.0 to podmiana biblioteki GDAL w instalacji QGIS nie zadziała. Ale jeśli to podmiana GDAL 1.7.0 na np. 1.7.2, to powinno zadziałać. Idzie o kompatybilność binarną bibliotek. Podmiana polega na skopiowaniu gdal.dll. Pozdrawiam, -- Mateusz Loskot, http://mateusz.loskot.net pl.comp.lang.c FAQ: http://pl.cpp.wikia.com/wiki/FAQ C++ FAQ: http://parashift.com/c++-faq-lite
Technical Details
Message-ID:
<ibjb3m$mr7$1@inews.gazeta.pl>Path:
mbox-mo1.pugleaf.net!news-archive.icm.edu.pl!news.gazeta.pl!not-for-mailReferences:
<ibhm6d$8nb$1@usenet.news.interia.pl> <ibhoib$82d$1@inews.gazeta.pl> <ibi1q2$rb0$1@usenet.news.interia.pl> <ibjami$so5$1@usenet.news.interia.pl>Show Full Headers
Complete Article Headers
Path: mbox-mo1.pugleaf.net!news-archive.icm.edu.pl!news.gazeta.pl!not-for-mail From: "Mateusz Loskot" <mateusz@loskot.net> Newsgroups: pl.comp.gis Subject: Re: QGIS i błąd reprojekcji Date: Fri, 12 Nov 2010 12:16:40 -0000 Organization: "Portal Gazeta.pl -> http://www.gazeta.pl" Lines: 56 Message-ID: <ibjb3m$mr7$1@inews.gazeta.pl> References: <ibhm6d$8nb$1@usenet.news.interia.pl> <ibhoib$82d$1@inews.gazeta.pl> <ibi1q2$rb0$1@usenet.news.interia.pl> <ibjami$so5$1@usenet.news.interia.pl> NNTP-Posting-Host: 82.110.183.162 Mime-Version: 1.0 X-Trace: inews.gazeta.pl 1289564086 23399 82.110.183.162 (12 Nov 2010 12:14:46 GMT) X-Complaints-To: usenet@agora.pl NNTP-Posting-Date: Fri, 12 Nov 2010 12:14:46 +0000 (UTC) In-Reply-To: <ibjami$so5$1@usenet.news.interia.pl> X-MimeOLE: Produced By Microsoft MimeOLE V6.0.6002.18263 X-Priority: 3 X-Newsreader: Microsoft Windows Mail 6.0.6002.18197 X-User: mloskot X-MSMail-Priority: Normal Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit