🚀 go-pugleaf

RetroBBS NetNews Server

Inspired by RockSolid Light RIP Retro Guy

Article View: pl.comp.gis
Article #12345

Re: QGIS i błąd reprojekcji

#12345
From: Mateusz Loskot
Date: Fri, 12 Nov 2010 12:16
55 lines
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-mail
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>

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