[Talk-cz] pLPIS a WFS - jak na to?

Čtvrtek Květen 29 06:57:01 UTC 2014

Použij buď geotools (java) nebo gdal + knihovnu proj (C++) je to jedno, důležitý
je výsledek.

Velice ve zkratce: v ČR používáme souř. systém S-JTSK (systém jednotné sítě
trigonometrické katastrální), které používá tzv. Křovákovo zobrazení, což je
kužel, s vrcholem někde u Finska a ta čísla jsou metry.

V OSM jsou v systému WGS84 - ta čísla jsou ve stupních.

A zobrazení OSM (dlaždice) je v tzv. Sperhical Mercator.

Mezi těmito systémy žijeme.

------------------------------

Křovákovo zobrazení je zvláštní tím, že data jsou "ve třetím kvadratnu" jakoby
karteziánského systému - viz 

https://www.google.com/search?q=s-jtsk&client=ubuntu-browser&es_sm=93&source=lnms&tbm=isch&sa=X&ei=PdeGU7G7M6m47Ab0q4GoAQ&ved=0CAgQ_AUoAQ&biw=1559&bih=815#facrc=_&imgdii=_&imgrc=22_Rt53VkcYHyM%253A%3B-CkXk2l6bWfqpM%3Bhttp%253A%252F%252Fimg.geocaching.com%252Fcache%252F31939f54-f3d0-4850-b5e6-bfae88e5bf52.jpg%3Bhttp%253A%252F%252Fwww.geocaching.com%252Fseek%252Fcache_details.aspx%253Fguid%253D43feb9fe-3b72-4085-b788-6c242ba6c2ed%3B475%3B360

Kvůli kdysi starému ArcView se původně pozitivní hodnoty prostě vynásobily -1 a
začalo to fungovat.

------------------------------

EPSG je všeobecně uznávaná databáze souř. systémů člověk se nemusí vypisovat s
parametry, stačí číslo.

"Negativní" křovák (tzv. křovák-gis) tam ale dlouho nebyl. Vedlo to k různým
řešením - používal se kód pro "normálního křováka", nebo se převzalo číslo z
ESRI ArcGIS. Proto se člověk může setkat s 2065, 102067. V posledních verzích
EPSG už je všeobecně platné a jediné doporučeníhodné číslo 5514.

EPSG pro WGS84 je 4326

EPSG pro zobrazení dlaždic OSM (spherical mercator, používá to hlavně google) je
3857.

------------------------------

A teď jak to dát dohromady:

LPIS nenabízí data v EPSG:4326, pouze v křovákovi-gis a ještě se špatným EPSG
kódem (102067 ale to je jedno). Je to na MapServeru, přihodit další souř.
systémy je otázka konfigurace (doslova 1 řádek), ale z firmy která to pásla mi
nikdo neodpověděl (sitewell).

Dobrá zpráva je, že souřadnice jsou negativní. 

Je potřeba tyhle data převézt na WGS84, aby byla použitelná pro JOSM.

Konzolově bych to dělal přes ogr2ogr (která využívá knihovnu proj4), ale my to
asi opravdu chceme přes geotools (aby se to dalo zaintegrovat do JOSM).

Jestli Ti to převádí správně zjistíš podle čísel na http://freegis.fsv.cvut.cz/gwiki/S-JTSK

Nezajímá nás chyba v řádu centrimetrů, ale nesmí to lítat o desítky metrů.

Ještě mě napadlo možné řešení: nemá někdo server, přes který bychom to vařili do
EPSG:4326? Prostě nahulána nastavit mapserver jako WFS klient pro LPIS a WFS
server pro OSM? Tohle bych uměl, pokud by mě někdo pustil k serveru (ssh) nebo
bychom se mohli i domluvit, že bych nastavil mapfile a někomu ho poslal.

Jachym

On Thu, May 29, 2014 at 06:58:45AM +0200, Marián Kyral wrote:
>    ---------- Původní zpráva ----------
>    Od: hanoj <ehanoj na gmail.com>
>    Komu: OpenStreetMap Czech Republic <talk-cz na openstreetmap.org>
>    Datum: 28. 5. 2014 23:18:03
>    Předmět: Re: [Talk-cz] pLPIS a WFS - jak na to?
> 
>      > Odpovím si sám ;-)
>      *** to mě těší, že si tu wiki o JTSK nakonec někdo přečte
> 
>    Já ji četl několikrát. A pořád mi to připomíná Woodoo. A nedošlo mi, že
>    EPSG:102067 = S-JSTK. Naivně jsem předpokládal, že když v kódu nahradím
>    jedno EPSG za druhé, tak mi to buď spočítá dobře, nebo vyhodí nějakou
>    chybu. Ani jedno nenastalo.
> 
>      >>
>      http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4&BBOX=-904539,-1227290,-731680,-935232&SRSNAME=EPSG:102067
>      > Trochu jsem si s tím hrál. Našel jsem nějaký příklad, jak pomocí
>      > geotools převádět mezi projekcemi. Nějak to funguje, ale nevím, jestli
>      > dobře.
>      *** netrap se s geografickými spouřadnicemi. Vytahej data WFS alis GML
>      v krásně planárním systému S-JTSK a konečný výsledek pak převeď do
>      WGS84 krze GDAL/ogr2ogr.
> 
>    Huh, další Woodoo.
> 
>    Plánuji podporu pro PointInfo a Tracer. Tedy jednorázový dotaz z JOSM.
>    Hromadný import "zatím?" není v plánu.
> 
>      Ještě by nás mělo před importem zajímat, co zdroj obsahuje a co ne,
>      jaká je jeho cca přesnost, jaká je licence těchto dat, a udělat o tom
>      zápis sem:
>      http://wiki.openstreetmap.org/wiki/WikiProject_Czech_Republic/freemap
> 
>    No už by to konečně chtělo, když to ještě nikdo za mne neudělal ;-)
> 
>    Marián
> 
>      ha
>      hanoj
> 
>      _______________________________________________
>      Talk-cz mailing list
>      Talk-cz na openstreetmap.org
>      https://lists.openstreetmap.org/listinfo/talk-cz

> _______________________________________________
> Talk-cz mailing list
> Talk-cz na openstreetmap.org
> https://lists.openstreetmap.org/listinfo/talk-cz

-- 
Jachym Cepicky
URL: http://les-ejk.cz
e-mail: jachym.cepicky at gmail com
PGP: http://les-ejk.cz/pgp/JachymCepicky.pgp
@jachymc
------------- další část ---------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 181 bytes
Desc: Digital signature
URL: <https://lists.openstreetmap.org/pipermail/talk-cz/attachments/20140529/9f00fde5/attachment.sig>