🚀 go-pugleaf

RetroBBS NetNews Server

Inspired by RockSolid Light RIP Retro Guy

Article View: pl.comp.gis
Article #12019

Re: Formuła transformacji z GPS do układu 65

#12019
From: Mateusz Loskot
Date: Sat, 27 Jun 2009 23:28
97 lines
5134 bytes
Juser.80 wrote:
> Mateusz Loskot pisze:
>>> We wskazanych przez ciebie linkach jest kilkadziesiąt plików i
>>> powiem szczerze, że nawet nie chce mi się ich przeglądać.
>> 
>> I własnie taki miałem cel. Chciałem Ci pokazać, że to nie jest taka
>>  prosta liniowa sprawa, dwa wzory, jeden algorytm i pstryk.
> 
> No bo przecież nie jest tak prosto zaimplementować bazę danych i jej 
> transformacje w jakimś języku, choć w accessie czy excelu jest to 
> dziecinnie proste.

Rozumiem, że mówisz o dziedzinnej prostocie mając na myśli również
np. normalizację owej bazy danych.

>>> Poważnie, doceniam je i dziękuję, jednak nie chciałbym go więcej 
>>> niepotrzebnie marnować.
>> 
>> Jak już napiszesz owe transformacje w Excel'u, to możesz ją
>> sprawdzić programem cs2cs z pakietu PROJ.4:
> 
> Skoro tak, to mam nadzieje że okażesz się również konsekwentny. Więc
> zaczynajmy: skąd wziąć wiarygodne cs2cs, a najlepiej cały PROJ.4

Coś mi się nie zgadza. My tu o inżynierii wyższej,
a jest kłopot z kojarzeniem faktów i dotarciem do strony projektu
PROJ.4, na przykład przez Google (wskazałem, że cs2cs jest w pakiecie
PROJ.4) i odnalezieniu na niej informacji jak pobrać pakiet z wersją
binarną tych programów.

Tu jest wszystko: http://trac.osgeo.org/proj/

>> Na przykład tak, współrzędne testowe z jednego z plików od Prof.
>> Kadaja GRS-80: 16.45663625    54.4967809944444 1965/3: 3465384.82
>> 6100829.38 -36.81
> 
> Co to jest -36.81 (to jest reszta lub coś?) i skąd się wzięło?

Jest to współrzędna wysokości (Z) dla danego punktu odniesiona do elipsoidy.

>> $ cs2cs +init=epsg:4326 +to +init=epsg:2173 16.45663625
>> 54.4967809944444 ...
> 
> Każdą współrzędna podaje z klawiatury, czy jest to jakoś 
> zautomatyzowane? Czy jest, a nie czy może być.

Można podać plik tekstowy z zapisanymi współrzędnymi
które zostaną przetransformowane i wynik transformacji zostanie wysłany
na standardowe wyjście, które można zkierować do pliku:

cs2cs +init=epsg:4326 +to +init=epsg:2173 input.txt > output.txt

Polecam lekturę podręcznika, bez tego "nierabotajet":

http://trac.osgeo.org/proj/wiki/man_cs2cs

>> i na wyjściu powinieneś dostać w/w współrzędne w 1965/3. W ten
>> sposób możesz sprawdzić swoje procedury.
> 
> Chętnie sprawdzę, jak je kiedyś dostanę, bo z tego co widzę, to każdy
>  coś tam słyszał, że to jest trudne, ale jest tam coś (programik) i 
> potrafi to liczyć i podobno wyjaśnione to jest w wytycznych
> technicznych G-1.10, choć nikomu chyba jeszcze to się nie wyjaśniło,
> wiec chyba jeszcze nikt nie wyszedł poza swoje konfabulacje :-o

Nie ma co zrzędzić, ale "zrób to sam" jak mawiał klasyk.
Jest szansa, że będziesz pierwszy. Kibicuję.

> I żebyśmy się dobrze zrozumieli, to nie jest moja branża i ja tu
> jestem stroną pytającą, wiec póki co, konfabulacje mnie nie dotyczą.

Byłbym ostrożny więc, bo na wstępie "zjechałeś" propozycję
użycia programu, że to lameriada jakich mało i prawdziwi spece
gryzą temat ołówkiem na kartce papieru. To jest konfabulacja, jeśli
chodzi o branżę.

> PS. doraźnie cs2cs wydaje się pomocny, jednak czy wiesz w ogóle, jak
> on liczy, czy tylko wiesz, że liczy?

Ło matko! Ale Cię męczy :-)

Znam algorytm transformacji dla układów polskich, które to opisał Prof.
Kadaj. Ale jeśli pytasz o to czy mam na pamieć wykute wszystkie wzory
oraz tablice parametrów i poprawek, to nie, nie mam.
I uwierz, aktualny autor tej biblioteki - Frank Warmerdam,
z którym pracowałem przez 3 ostatnie lata też nie ma.
Od tego jest dokumentacja, opis odwzorowań, elipsoid, i innych elementów.

Całość tematu jest tak obszerna, że sam się "wystraszyłeś" patrząc na to
jak obszerna jest implementacja:

http://trac.osgeo.org/proj/browser/trunk/proj/src

Pozdrawiam
-- 
Mateusz Loskot, http://mateusz.loskot.net
Charter Member of OSGeo, http://osgeo.org
Technical Details
Message-ID: <h269v3$8pi$1@inews.gazeta.pl>
Path: mbox-mo1.pugleaf.net!news-archive.icm.edu.pl!news.gazeta.pl!not-for-mail
References: <h1pqko$k6r$1@news.onet.pl> <h1vivn$ptg$1@inews.gazeta.pl> <h20dt8$500$1@news.onet.pl> <h20lvj$4jt$1@inews.gazeta.pl> <h21p1d$7n8$1@news.onet.pl> <h238g5$cuo$1@inews.gazeta.pl> <h23dil$3r6$1@news.onet.pl> <h23h88$l16$1@inews.gazeta.pl> <h2518d$bf0$1@news.onet.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 <see.my@signature.net>
Newsgroups: pl.comp.gis
Subject: Re: Formuła transformacji z GPS do
 układu 65
Date: Sun, 28 Jun 2009 00:28:34 +0100
Organization: "Portal Gazeta.pl -> http://www.gazeta.pl"
Lines: 97
Message-ID: <h269v3$8pi$1@inews.gazeta.pl>
References: <h1pqko$k6r$1@news.onet.pl> <h1vivn$ptg$1@inews.gazeta.pl> <h20dt8$500$1@news.onet.pl> <h20lvj$4jt$1@inews.gazeta.pl> <h21p1d$7n8$1@news.onet.pl> <h238g5$cuo$1@inews.gazeta.pl> <h23dil$3r6$1@news.onet.pl> <h23h88$l16$1@inews.gazeta.pl> <h2518d$bf0$1@news.onet.pl>
NNTP-Posting-Host: 5ad618f9.bb.sky.com
Mime-Version: 1.0
X-Trace: inews.gazeta.pl 1246145315 9010 90.214.24.249 (27 Jun 2009 23:28:35 GMT)
X-Complaints-To: usenet@agora.pl
NNTP-Posting-Date: Sat, 27 Jun 2009 23:28:35 +0000 (UTC)
X-User: mloskot
In-Reply-To: <h2518d$bf0$1@news.onet.pl>
User-Agent: Thunderbird 2.0.0.22 (X11/20090608)
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit