25 odpovědí na toto téma

[jeleniste]

Mám dotaz, jakou máte zkušenost s přesností zalícování různých ortofotomap do systému souřadnic? Porovnával jsem několik snímků z různých serverů a narazil na nezanedbatelné rozdíly. V metrickém systému souřadnic (S-JTSK) jsou to až desítky metrů (určení souřadnic unikátního objektu, např křižovatka), jaká je vaše zkušenost? Při převodu na WGS se ještě přičte chyba vzniklá při přepočtu... Netušíte co to způsobuje? Je to lajdáctví při referencování snímků, důsledek zkreslení při focení "pod úhlem", lepení placatejch snímků na kulatou zemegulu? Dík

[Tom Fuk]

Možná je to tím, že ty mapy nejsou primárně určené k přesnému odečtu souřadnic S-JTSK. Tento souřadnicový systém tam mohl být napasován dodatečně a třeba jen pro hrubší orientaci. Napiš, mezi kterými mapami jsi zjistil ten rozdíl v desítkách metrů.

[sobikovi]

Nevím čím to je, ale mám docela dojem, že např. na mapy.cz je ortofoto mapa dost přesná a ukazuje možná odchylku jen pár metrů. Důležitý ale je asi to, co píše Tom Fuk, že záleží na tom, k čemu jsou ortofoto mapy určené... Jestliže mám na mapě rozlišovací schopnost na např. 10m, tak musím počítat i s takovou chybou..

[Mion]

Porovnával jsem mapy na altas.cz a na mapy.cz. Je zvlastni, ze mapovy podklad na Atlasu vcetne fotomap pri priblizeni nelicuje presne na sever - tzn. ze smer severojizni neni presne svisly. Vsimnete si toho, kdyz si zobrazite v priblizeni stejne misto v prohlizeci a prepinate mezi pohledy. Kraje mapy jsou pak posunuty ve svislem smeru o desitky metru (dle zvoleneho meritka). Mozna to ma takhle Atlas schvalne - rozhodne je dost velky problem je napr. dostatecne presne zkalibrovat v OziE

[jeleniste]

Mysis od Gepro (katastrálka + ortofoto) pro obecni úřady, jtsk+wgs režim, ochranařske mapy na CENII a ortofotomapa na IZGARD, mapa primarně nemůže být ve wgs, na izgard je nakreslena siť souřadnic, osobně věřim víc vojákům, než gepru, krom toho na tom mysisu se místama nekreje ortofoto s katastrálkou.. Desitky metru, to zase nebyly, maximalně jedna desitka, ale i to mi připadá dost. Když k tomu přidám chybu gps, už jsem např na druhém břehu potoka a když bych chtěl vyznačovat něco v mapě, je to už dost slabý.

[jeleniste]

Nageoreferencovat v Arcu, to problem neni, ale je treba mit neco jako pevnej zaklad, ledaže bych si stoupl na křižovatku, zaměřil souřadnice a pak to na ně referencoval

[javla]

mapy.cz mají mapu mírně natočeno. Zabývám jsem se stahováním a sestavováním map z mapy.cz a googlemaps.com pro osobní potřebu při geocachingu. Pomocí vhodných scriptů jsem si stáhl i gps souřadnice středů každého čtverce 256x256 bodu a vypadá to takto: Google má v řadách a sloucích vždy longtitude či latitude shodnou - tedy z pořadového čísla mapky ve směru osy y či y vždy dostanu stejnou long/lat souřadnici. Na mapy.cz když vezmu čtveřici k sobě pasujících map (čtverec 2x2 mapy po 256x256 bodech), tak souřadnice středu vypadají například takto: <cal x="2176" y="2304" lon="49.286536667" lat="13.112213333"/> <cal x="2176" y="2305"lon="49.295742778" lat="13.111861667"/> <cal x="2177" y="2304"lon="49.286765833" lat="13.126286944"/> <cal x="2177" y="2305"lon="49.295971944" lat="13.125937778"/> Je to prostě mírně natočený. Při kalibraci 4 body (po sestavení do větší plochy) je to samozřejmě OK - mapové softwary v PDA si s tím poradí.

[jeleniste]

na izgardu to stáhnes v PNG i se souborem georef.

[alkes]

Mion: A není to právě kvůli kalibraci v S-JTSK? Mám pocit, že ten směr severu na atlasu sedí právě s se severem S-JTSK, odpovídal by tomu i fakt, že podklady jsou z CUZK. Nechce se mi ted ale vymyslet zpusob jak to overit na pocitaci ve skole, hledat stranky bez odkazu a tak... jeleniste: posli prosim odkazy na konkretni mapy edit: tak koukam, ze nez jsem dopsal, uz tu jsou

[jeleniste]

S-jtsk nema sever, mas vrchol kuzelovy cepicky referencni plochy nekde tusim nad baltem, je to tzv krovakova projekce, takovej kornout, kterej protina elipsoid ve dvou soustrednejch siskach, vrchol neni na severnim polu. odkazy sou uz ve foru digi atlas cr

[jt-k]

Tak jsem si kliknul na komín našeho domečku (okapová deska komínu je 60x60) postupně na ortomapy amapy.atlas.cz, mapy.1188.cz, mapy.cz a izgard.cenia.cz. Pak jsem vystrčil GPS ze střešního okna- 11 satelitů, nasázel údaje do OziExploreru a musím říct, že rozdíly jsou naprosto minimální od 60cm do 2,20 m. No a to jsem nelezl na komín ten je od okna asi tak dva metry. Amapy.cz a mapy.1188.cz koukám používají stejné ortosnímky. Spíš mě zaráží jednoduchost zjišťování dat o občananech pri spojení mapy.1188.cz a nahlizenidokn.cuzk.cz/. Např. křestní jméno sousedky, pokud je spoluvlastníkem pozemku máte za minutku. Už čekám, že někde kliknu a budu vědět číslo své kreditní karty ... , že tam nic moc není, to už vím !

[jeleniste]

napis mi souradnice svyho komina a ja ti reknu tvy rodny cislo chces??

[jt-k]

jeleniste napsal/a:
napis mi souradnice svyho komina a ja ti reknu tvy rodny cislo chces??

Jó, je to pravda Google ho vyplivne podle jména za 5 vteřin... :@

[Tom Fuk]

jeleniste napsal/a:
...IZGARD, mapa primarně nemůže být ve wgs, na izgard je nakreslena siť souřadnic, osobně věřim víc vojákům, než gepru ...
...krom toho na tom mysisu se místama nekreje ortofoto s katastrálkou...

Já osobně bych věřil spíš Gepru, protože to je geodetická firma. Vojáci měli mapy v S-42 a S-JTSK se tam museli nějak napasovat. Vzhledem k tomu, že ty dva systémy nejsou kompatibilní, napasování může být jen přibližné.
Podobný problém je v tom že se ortofoto snímky nekryjí s katastrálkama (1:2880) - je to jiný typ zobrazení.

[jeleniste]

na www.uhul.cz sou taky ortofoto podklady, už černobily, lepsi na tisteni

[jeleniste]

http://www.omks.cz/c...?typ=html&id=80

[alkes]

jeleniste napsal/a:
S-jtsk nema sever, mas vrchol kuzelovy cepicky referencni plochy nekde tusim nad baltem, je to tzv krovakova projekce, takovej kornout, kterej protina elipsoid ve dvou soustrednejch siskach, vrchol neni na severnim polu.
odkazy sou uz ve foru digi atlas cr

Do principu úplně nevidím, o Křovákově zobrazení vím v podstatě jen že existuje a že je použito v S-JTSK, takže jsem se asi špatně vyjádřil.
Důsledkem zřejmě je, že mapy v S-JTSK (přinejmenším ty používané na ČÚZK) jsou natočené rovnoběžně se souřadnými osami a zeměpisný sever potom leží pod azimutem 8° od vertikální osy mapy (aspoň tedy na zeměpisné délce přibližně procházející Prahou... zkoušel jsem to u Mácháče, v Praze a u Dačic). Názorně je to vidět tady
A přesně tak tomu je i u map na atlasu... teď už ověřeno.
Mimochodem ta čísla v adresovém řádku jsou práve S-JTSK souřadnice středu mapy. (při větších zooomech, při oddálení např. na celou Evropu už je to S42)
Kalibraci by to problém dělat nemělo, jen bude pracnější o přepočty, i ty ale jdou zautomatizovat.
edit:Jéje, už ani úhly odčítat neumim, když mi to neukazuje digitálně Takže ne 3-4 stupně, ale přibližně 8. :|

[FHonza]

Dobře je popsán princip Křovákova zobrazení zde:
http://grass.fsv.cvu...ndex.php/S-JTSK

[rsc_cz]

souřadnicový systém JTSK je v zobrazovací rovině (na mapě) pravoúhlý (WGS84 tvoří lichoběžníky) a pouze souřadnicová soustava v bodě [0,0] (leží za Zakarpatskou Ukrajinou (mimo území první republiky)) odpovídá směru místního poledníku. Čím je člověk od tohoto bodu dál, tím je rozdíl směru těchto os větší a u Prahy to dělá asi 10° (říká se tomu meridiánová konvergence).

[supermudla]

javla napsal/a:
Na mapy.cz když vezmu čtveřici k sobě pasujících map (čtverec 2x2 mapy po 256x256 bodech), tak souřadnice středu vypadají například takto:

<cal x="2176" y="2304" lon="49.286536667" lat="13.112213333"/>
<cal x="2176" y="2305"lon="49.295742778" lat="13.111861667"/>
<cal x="2177" y="2304"lon="49.286765833" lat="13.126286944"/>
<cal x="2177" y="2305"lon="49.295971944" lat="13.125937778"/>

Je to prostě mírně natočený. Při kalibraci 4 body (po sestavení do větší plochy) je to samozřejmě OK - mapové softwary v PDA si s tím poradí.

Je to nejen natočený, ale i nelineární. Já dělám něco podobnýho, a pro rozumnou přesnost jsem to musel proložit polynomem třetího stupně, tohle dává už docela použitelný výsledky.

def degrees_to_tiles(x, y):
    x -= 48.205
    y -= 15.623
    return int(
        135881231.97 -28829.5454549*x +2377560.57798*y
        -234.288249202*x*x -46283.8135727*x*y -13.5277601485*y*y
        +3.62894668673*x*x*x -357.288355627*x*x*y -18.2996927239*x*y*y -17.2772303351*y*y*y), int(
        129256743.95 +3556760.37208*x +19264.5205545*y
        +312.885409175*x*x -36.5956770677*x*y +15472.6199945*y*y
        -0.501076953863*x*x*x -30.4036871386*x*x*y -79.3150032588*x*y*y +4.21831932515*y*y*y)