Re: [Talk-cz] pLPIS a WFS - jak na to?
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 >Komu: OpenStreetMap Czech Republic >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@openstreetmap.org > https://lists.openstreetmap.org/listinfo/talk-cz > ___ > Talk-cz mailing list > Talk-cz@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 signature.asc Description: Digital signature ___ T
Re: [Talk-cz] Nový S-JTSK grid, kdo pomůže?
> obávám se že ČÚZK nemůže - ani jeden z těchto systémů nemá pod svou správou > a schvalování algoritmů se týká pouze ETRS (možná bude někdo jiný vědět > víc). Ale přítel Google mi hned jako první odkaz dal > http://transformace.webst.fd.cvut.cz/Iframe/WGS_ETRS_iframe.htm ... *** ze je 7/14 prvkova je asi pochopitelne, ale tech klicu jsem na googlu nasel pro ruzne epochy a realizace desitky, např [1]. Takze by mi velmi pomohlo ukazat prstem, kterou verzi pouziva CUZK. nebo alespon popsat to slovne např.: ETRS89-ETRF2000 R8 epocha 2000.0 a WGS84-G1150 epocha 2000.0: [1] http://www.defensie.nl/binaries/defence/documents/circulars/2013/08/23/transformation-parameters-between-itrs-wgs84-and-etrs89/trafoparsitrs_v2013en.xls ha hanoj ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] pLPIS a WFS - jak na to?
-- Původní zpráva -- Od: hanoj Komu: OpenStreetMap Czech Republic 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@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz";___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Nový S-JTSK grid, kdo pomůže?
> ano, mapové služby (jak na geoportal.cuzk.cz, tak na services.cuzk.cz - WMS > i WFS nebo GML) obsahují oficiálně schválené algoritmy pro transformaci *** tak zda se, ze jadro pudla pri prevodu gridem[1] je v tom, ze tento grid dela transf EPSG:5514 <-> EPSG:4258 (a pribuzne), kdezto my v OSM pouzivame EPSG:4326 (a pribuzne) a chybi nam popsat vztah EPSG:4258 a EPSG:4326 tj. ETRS89 a WGS84. Muze CUZK tento transf postup ETRS89 <-> WGS84 (predpokladam ze je to 7 prvkovy klic) pro aktualni verzi ETRF2000 publikovat? priklad viz: http://pastebin.com/rDgjErEt ha hanoj [1] http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Grid ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Nový S-JTSK grid, kdo pomůže?
> Nejprve převeďme tvé výsledky na počítačovější formu. Jestli počítám správně, > tak tvé výsledky pro transformaci bez gridu jsou: > ... > 2.520413388 > to je trochu moc, ne? 2.5 metru? *** myslím, že jsi opsal špatně jednu číslici jedné souřadnice > A teď moje výsledky: > > 0.747197172 metru *** souhlas http://pastebin.com/eZ7Km2Mx > Dejme tomu, že je posun při transformaci s gridem správným směrem a že je > tedy výsledek přesnější než bez gridu. Když si dáme v JOSM přes sebe > průhlednou KM > a k tomu moje vrstvy budov, proč se tedy lépe kryjí růžové budovy (t.j. bez > gridu) než budovy modré (t.j. s gridem) ???Růžové se lépe kryjí i s tím, co je > v OSM. Znamená to, že máme vše v OSM špatně? A hlavně - co s tím dál? *** Nejprve je treba vedet co je spravne a proc. Proto jsem vznesl oficiální dotaz na CUZK ohledne transformaci. ha hanoj ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
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 >> 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. 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 ha hanoj ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] pLPIS a WFS - jak na to?
Odpovím si sám ;-)
-- Původní zpráva --
Od: Marián Kyral
Komu: talk-cz@openstreetmap.org
Datum: 24. 5. 2014 23:54:56
Předmět: Re: [Talk-cz] pLPIS a WFS - jak na to?
"Dne 16.5.2014 01:27, hanoj napsal(a):
>> Tak jsem měl teď přes oběd trochu času a koukal jsem na WFS.
>>
>> Základní URL:
>>
>> Prvních 200 položek:
>>
>> No a s tím mám jen dva (dost podstatné ;-) ) problémy
> *** a což zkusit project ID?
>
> http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&;
REQUEST=GetFeature&SRSNAME=EPSG:102067&TYPENAME=LPIS_FB4&featureID=LPIS_FB
4.8780373
To právě musím nejprve někde zjistit ne? Mám akorát bod, na který bylo
kliknuto.
> nebo BBOX?
>
> 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.
Například:
proN49.660144848959035, E18.29907177958237 (EPSG:4326)
dostanu 471210.46759085247, 1121224.5497619768 (EPSG:102067)
V tom URL nahoře jsou ale všechny souřadnice BBOXU záporné. Tady ne. Je
to dobře?
Když se třeba kouknu, co za request posílá LPIS, tak tam to je taky
záporné. Navíc je vidět, že vypočítané souřadnice by cca souhlasily.
(neklikl jsem na stejné místo).
http://eagri.cz/public/app/lpisext/lpis/mapserver/cgi-bin/tiskove_hladiny.
fcgi?LPIS_FB4_STAV=4&service=WMS&version=1.1.0&request=GetFeatureInfo&layers
=LPIS_KP4%2CLPIS_FB4&query_layers=LPIS_KP4%2CLPIS_FB4&bbox=-471131.845377%2C
-1121174.062338%2C-470688.051172%2C-1120897.837487&srs=EPSG%3A102067&feature
_count=1&x=324&y=414.612988281&height=783&width=1258&info_format=
application%2Fvnd.ogc.gml
"
Dle očekávání to není dobře. Nevěděl jsem, co přesně to EPSG:102067 znamená.
A ono to je ESRI 102067 alias EPSG:5514 alias S-JTSK / Krovak East North [1]
Takže pro převod lze použít cs2cs [2]
No a ten pro N49.660144848959035, E18.29907177958237 vrátí:
$ echo "18.29907177958237 49.660144848959035" |cs2cs -f "%.2f" +proj=longlat
+datum=WGS84 +to +proj=krovak +ellps=bessel +nadgrids=czech
-471084.50 -1121102.96 0.00
Což vypadá slibně.
Následně jsem pomocí opačné transformace převedl data z LPIS requestu výše
do WGS 84 a následně pomocí UtilsPlugin2 [3] (LatLon tool) vytvořil v JOSM
polygon, který se téměř kryje s wms vrstvou LPIS [4] (odchylka cca 10cm) (y)
No a teď už jen vyřešit ten zatrolený BBOX. Ale když tak nad tím uvažuji,
tak bych vlastně měl mít všechno v JOSM. Jak BBOX, tak height/width tak i x/
y. Jen si to vzít.
No vypadá to nadějně :-D
Marián
[1] http://geoportal.cuzk.cz/%28S%28vfttsk45dozkiyjz5itt15ew%29%29/Default.
aspx?mode=TextMeta&side=sit.trans&text=souradsystemy
[2] http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Grid
[3] http://josm.openstreetmap.de/wiki/Help/Plugin/UtilsPlugin2
[4] wms:http://eagri.cz/public/app/wms/plpis.fcgi?FORMAT=image/png&VERSION=
1.1.1&SERVICE=WMS&REQUEST=GetMap&LAYERS=LPIS_FB4_KOD&STYLES=&SRS={proj}&
WIDTH={width}&HEIGHT={height}&BBOX={bbox}&TRANSPARENT=true
___
Talk-cz mailing list
Talk-cz@openstreetmap.org
https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Nový S-JTSK grid, kdo pomůže?
Zdravím všechny, ano, mapové služby (jak na geoportal.cuzk.cz, tak na services.cuzk.cz - WMS i WFS nebo GML) obsahují oficiálně schválené algoritmy pro transformaci (http://www.cuzk.cz/Zememerictvi/Geodeticke-zaklady-na-uzemi-CR/GNSS/ Programy-pouzitelne-pro-data-ziskana-pomoci-GN-%281%29.aspx). J. Veselý -- Původní zpráva -- Od: Martin Kokes Komu: talk-cz@openstreetmap.org Datum: 28. 5. 2014 18:04:45 Předmět: Re: [Talk-cz] Nový S-JTSK grid, kdo pomůže? "Kdybych nebyl zrovna na dovolené, určitě bych přišel poznat české OSMaře a ČÚZKy na vlastní oko... Tak třeba byste si o tom mohli nějak popovídat sami, cílem by mělo být dosáhnout přesnosti oficiální transformační služby http://geoportal.cuzk.cz/Default.aspx?head_tab=sekce-01-gp&mode=TextMeta&; text=wcts&menu=19 v prostředí PROJ4 a GDAL. Té IMHO ve vztahu k využitelnosti dat v OSM dosahuje jen WFS ČÚZKu - služba INSPIRE (CP, případně AU, AD). Předpokládám tedy, že tam Petr Souček má zabudovanou jejich oficiální transformaci, když nabízí data i ve WGS84 http://services.cuzk.cz/wfs/inspire-cp-wfs.asp?service=WFS&request= GetCapabilities ... MK "Martin Landa" píše v diskusním příspěvku news:ca+ei1od97auwc_8b4+nyo2x+tyornqb8jm8qb9qtxc_+mgn...@mail.gmail.com... Dne 28. května 2014 17:26 Martin Kokes napsal(a): > Netuším to je spíš otázka na Ondřeje Chlupa - chtěl to asi prostě spočítat > v > PGSQL v R :-). Proto jsem ostatně chtěl "urychlit" výsledek a prostě vzít > DATAZ body a jít nějakou jinou cestou (Jan Ježek-way) nebo způsob > naznačený > v http://lists.maptools.org/pipermail/proj/2013-January/006539.html > > Ale sám to nezvládnu, já jsem spíš IT allrounder, tohle je na mne už vyšší > GDAL dívčí. tema na diskuzi na Geoinformatics [1] ? ;-) Martin [1] http://geoinformatics.fsv.cvut.cz/gwiki/Geoinformatics_FCE_CTU_2014 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz";___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Nominatim a české adresy
Dne 28.5.2014 16:09, Matěj Cepl napsal(a): On 2014-05-27, 21:13 GMT, Petr Souček wrote: A dále je členěna na několik prvků: - na 57 městských částí http://vdp.cuzk.cz/vdp/ruian/mestskecasti/vyhledej?ob.kod=554782&search=Vyhledat - na 22 správních obvodů Prahy http://vdp.cuzk.cz/vdp/ruian/spravniobvody/vyhledej?search=Vyhledat - na 10 městských obvodů Prahy http://vdp.cuzk.cz/vdp/ruian/mestskeobvody/vyhledej?search=Vyhledat - na 112 katastrálních území http://vdp.cuzk.cz/vdp/ruian/katastralniuzemi/vyhledej?ob.kod=554782&search=Vyhledat (v rámci katastrálních území jsou v Praze číslovány stavební objekty, takže se někde můžete setkat s tím, že v Praze je také 112 částí obce) No právě a myslím, že by asi bylo dobré rozhodnout co z toho bude obec pro účely OSM. Teď tam máme Prahu jako jednu obec což mi připadne (z důvodů, které jsou snad jasné z tohoto threadu) špatně. Nechtělo by to nahradit toto členění něčím jiným? Čím? Já bych hlasoval pro Prahu 1-10 (nebo 1-22), ale chtěl bych se o tom spíše pobavit (a kdyby někdo navrhnul jak změnit OSM data tak, aby Nominatim chápal obce jako co se dohodneme). Probuh jen to ne ... uvedom si laskave, ze to neni mapa pro prazaky ... a naprosto nikoho jinyho tohle cleneni nezajima (a to vcetne zbytku republiky) nebo vubec netusi ze existuje. Ty kdyz budes hledat adresu v Berline/Vidni/... tak taky nechces resit, v jaky mestsky casti to je. Napises "Hitler strase, Berlin" ;D. Navic pokud vim se to tu resilo uz nekolikrat, a doslo se k tomu, ze je to i technicky nerealizovatelny, protoze ty hranice nikdo realne nezna, a pokud existuji, tak vedou mimo existujici system katastralnich uzemi. Tusim ze jsou dokonce "definovany" vyctem ulic/adres ... ktera kam spada => neexistuje hranicni cara. Hezký den, Matěj ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Nominatim a české adresy
Dne 27.5.2014 17:53, Matěj Cepl napsal(a): On 2014-05-26, 18:19 GMT, Petr Morávek [Xificurk] wrote: Osobně jsem se (zjevně mylně) domníval, že číslo městského obvodu do adresy "nepatří". Dokonce i formulář pro ověření adresy [1] obsahuje jen políčko "Obec", nikoliv "Městský obvod", takže je potřeba zadat jen "Praha" bez čísla, v opačném případě sice vyhledávání zafunguje, ale objeví se upozornění: Z mé temné právnické minulosti (ale už je to fakt docela dávno) si mlhavě pamatuji, že zákon o hl. m. Praze (nebo zákon o obcích?) praví, že obec je v Praze městský obvod (tj. Praha 1-22? chtělo by to sehnat aktuální znění) a ne Praha jako celek. Takže pokud chtějí "Obec" tak je správná odpověď "Praha 3", nikoli "Praha". Praha jako celek je kraj (a taky nemá starostu ale primátora). Matěj Sem si nevsim ze by na vjezdech do praglu byly znacky "PRAHA 10" ... Navic i pokud se pokusis zjistit, kudy ze vlastne vedou hranice mestskych casti, tak zjistis, ze to nevedi ani urednici (staci trochu zagooglit a najdes desitky pripadu kdy pri reseni ruznych veci kolem nemovitosti posilaj urednici lidi od certa k dablu). ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Nový S-JTSK grid, kdo pomůže?
Kdybych nebyl zrovna na dovolené, určitě bych přišel poznat české OSMaře a ČÚZKy na vlastní oko... Tak třeba byste si o tom mohli nějak popovídat sami, cílem by mělo být dosáhnout přesnosti oficiální transformační služby http://geoportal.cuzk.cz/Default.aspx?head_tab=sekce-01-gp&mode=TextMeta&text=wcts&menu=19 v prostředí PROJ4 a GDAL. Té IMHO ve vztahu k využitelnosti dat v OSM dosahuje jen WFS ČÚZKu - služba INSPIRE (CP, případně AU, AD). Předpokládám tedy, že tam Petr Souček má zabudovanou jejich oficiální transformaci, když nabízí data i ve WGS84 http://services.cuzk.cz/wfs/inspire-cp-wfs.asp?service=WFS&request=GetCapabilities ... MK "Martin Landa" píše v diskusním příspěvku news:ca+ei1od97auwc_8b4+nyo2x+tyornqb8jm8qb9qtxc_+mgn...@mail.gmail.com... Dne 28. května 2014 17:26 Martin Kokes napsal(a): Netuším to je spíš otázka na Ondřeje Chlupa - chtěl to asi prostě spočítat v PGSQL v R :-). Proto jsem ostatně chtěl "urychlit" výsledek a prostě vzít DATAZ body a jít nějakou jinou cestou (Jan Ježek-way) nebo způsob naznačený v http://lists.maptools.org/pipermail/proj/2013-January/006539.html Ale sám to nezvládnu, já jsem spíš IT allrounder, tohle je na mne už vyšší GDAL dívčí. tema na diskuzi na Geoinformatics [1] ? ;-) Martin [1] http://geoinformatics.fsv.cvut.cz/gwiki/Geoinformatics_FCE_CTU_2014 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Nový S-JTSK grid, kdo pomůže?
Dne 28. května 2014 17:26 Martin Kokes napsal(a): > Netuším to je spíš otázka na Ondřeje Chlupa - chtěl to asi prostě spočítat v > PGSQL v R :-). Proto jsem ostatně chtěl "urychlit" výsledek a prostě vzít > DATAZ body a jít nějakou jinou cestou (Jan Ježek-way) nebo způsob naznačený > v http://lists.maptools.org/pipermail/proj/2013-January/006539.html > > Ale sám to nezvládnu, já jsem spíš IT allrounder, tohle je na mne už vyšší > GDAL dívčí. tema na diskuzi na Geoinformatics [1] ? ;-) Martin [1] http://geoinformatics.fsv.cvut.cz/gwiki/Geoinformatics_FCE_CTU_2014 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Nový S-JTSK grid, kdo pomůže?
Netuším to je spíš otázka na Ondřeje Chlupa - chtěl to asi prostě spočítat v PGSQL v R :-). Proto jsem ostatně chtěl "urychlit" výsledek a prostě vzít DATAZ body a jít nějakou jinou cestou (Jan Ježek-way) nebo způsob naznačený v http://lists.maptools.org/pipermail/proj/2013-January/006539.html Ale sám to nezvládnu, já jsem spíš IT allrounder, tohle je na mne už vyšší GDAL dívčí. MK "Martin Landa" píše v diskusním příspěvku news:CA+Ei1Of6gPBYytQFHR1Ux7BWAYiBvoAPS5H=ww22t4brwo0...@mail.gmail.com... Dne 28. května 2014 15:54 Martin Kokes napsal(a): Ozval se mi p. Ondřej Chlup, co nový grid počítá, můj poslední mail se mu někde ztratil v korespondenci, konkrétně píše: "V současné době ještě není výsledný grid k dispozici. Zatím jsem ve stavu, že se podařilo vypočítat grid pro celou ČR, ale zatím bohužel vrací nepřesné výsledky. Metoda popsaná v mé diplomové práci je funkční, ale problém nastává u využití knihovny R pro práci s maticemi. Pro omezené území, na kterém byla metoda testována v rámci DP běží výpočet bez problému. Pokud se však výpočet provede pro body z celé ČR, tak jsem narazil na limit knihovny R, jelikož není možné alokovat vektor větší než 3, nebo 4 GB. Proto jsem musel přistoupit k rozdělení výpočtu na několik menších částí, které budou následně spojeny v jeden grid. Jsem v kontaktu s panem Ježkem, kterého informuji o aktuálních výsledcích. Jakmile tedy bude nový grid hotový, tak budu pana Ježka informovat a následně bude dán k dispozici veřejnosti." mozna nemistna otazka, ale proc proboha R-ko? Martin ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Nominatim a české adresy
Dne 28.5.2014 16:09, Matěj Cepl napsal(a): > On 2014-05-27, 21:13 GMT, Petr Souček wrote: >> A dále je členěna na několik prvků: >> - na 57 městských částí >> >> http://vdp.cuzk.cz/vdp/ruian/mestskecasti/vyhledej?ob.kod=554782&search=Vyhledat >> - na 22 správních obvodů Prahy >> http://vdp.cuzk.cz/vdp/ruian/spravniobvody/vyhledej?search=Vyhledat >> - na 10 městských obvodů Prahy >> http://vdp.cuzk.cz/vdp/ruian/mestskeobvody/vyhledej?search=Vyhledat >> - na 112 katastrálních území >> >> http://vdp.cuzk.cz/vdp/ruian/katastralniuzemi/vyhledej?ob.kod=554782&search=Vyhledat >> >> (v rámci katastrálních území jsou v Praze číslovány stavební >> objekty, takže se někde můžete setkat s tím, že v Praze je >> také 112 částí obce) > > No právě a myslím, že by asi bylo dobré rozhodnout co z toho > bude obec pro účely OSM. Teď tam máme Prahu jako jednu obec což > mi připadne (z důvodů, které jsou snad jasné z tohoto threadu) > špatně. Nechtělo by to nahradit toto členění něčím jiným? Čím? > > Já bych hlasoval pro Prahu 1-10 (nebo 1-22), ale chtěl bych se > o tom spíše pobavit (a kdyby někdo navrhnul jak změnit OSM data > tak, aby Nominatim chápal obce jako co se dohodneme). > > Hezký den, > > Matěj Tohle raději ne. Vzhledem k adresám mi přijde jako nejlepší naimportovat městské obvody Prahy jako relace hranic s admin_level=9. Jen pro úplnost uvádím, že admin_level=8 je u nás všude hranice obcí a admin_level=10 hranice katastrálních území (které jsou v Praze zároveň de facto částmi obce). Správní obvody Prahy, městské části Prahy a ostatních statutárních měst bych naopak zatím do OSM asi netahal. Petr Morávek ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Nominatim a české adresy
On 2014-05-27, 21:13 GMT, Petr Souček wrote: > A dále je členěna na několik prvků: > - na 57 městských částí > > http://vdp.cuzk.cz/vdp/ruian/mestskecasti/vyhledej?ob.kod=554782&search=Vyhledat > - na 22 správních obvodů Prahy > http://vdp.cuzk.cz/vdp/ruian/spravniobvody/vyhledej?search=Vyhledat > - na 10 městských obvodů Prahy > http://vdp.cuzk.cz/vdp/ruian/mestskeobvody/vyhledej?search=Vyhledat > - na 112 katastrálních území > > http://vdp.cuzk.cz/vdp/ruian/katastralniuzemi/vyhledej?ob.kod=554782&search=Vyhledat > > (v rámci katastrálních území jsou v Praze číslovány stavební > objekty, takže se někde můžete setkat s tím, že v Praze je > také 112 částí obce) No právě a myslím, že by asi bylo dobré rozhodnout co z toho bude obec pro účely OSM. Teď tam máme Prahu jako jednu obec což mi připadne (z důvodů, které jsou snad jasné z tohoto threadu) špatně. Nechtělo by to nahradit toto členění něčím jiným? Čím? Já bych hlasoval pro Prahu 1-10 (nebo 1-22), ale chtěl bych se o tom spíše pobavit (a kdyby někdo navrhnul jak změnit OSM data tak, aby Nominatim chápal obce jako co se dohodneme). Hezký den, Matěj ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Nový S-JTSK grid, kdo pomůže?
Dne 28. května 2014 15:54 Martin Kokes napsal(a): > Ozval se mi p. Ondřej Chlup, co nový grid počítá, můj poslední mail se mu > někde ztratil v korespondenci, konkrétně píše: > > "V současné době ještě není výsledný grid k dispozici. Zatím jsem ve stavu, > že se podařilo vypočítat grid pro celou ČR, ale zatím bohužel vrací nepřesné > výsledky. Metoda popsaná v mé diplomové práci je funkční, ale problém > nastává u využití knihovny R pro práci s maticemi. Pro omezené území, na > kterém byla metoda testována v rámci DP běží výpočet bez problému. Pokud se > však výpočet provede pro body z celé ČR, tak jsem narazil na limit knihovny > R, jelikož není možné alokovat vektor větší než 3, nebo 4 GB. Proto jsem > musel přistoupit k rozdělení výpočtu na několik menších částí, které budou > následně spojeny v jeden grid. Jsem v kontaktu s panem Ježkem, kterého > informuji o aktuálních výsledcích. Jakmile tedy bude nový grid hotový, tak > budu pana Ježka informovat a následně bude dán k dispozici veřejnosti." mozna nemistna otazka, ale proc proboha R-ko? Martin ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Nový S-JTSK grid, kdo pomůže?
Ozval se mi p. Ondřej Chlup, co nový grid počítá, můj poslední mail se mu někde ztratil v korespondenci, konkrétně píše: "V současné době ještě není výsledný grid k dispozici. Zatím jsem ve stavu, že se podařilo vypočítat grid pro celou ČR, ale zatím bohužel vrací nepřesné výsledky. Metoda popsaná v mé diplomové práci je funkční, ale problém nastává u využití knihovny R pro práci s maticemi. Pro omezené území, na kterém byla metoda testována v rámci DP běží výpočet bez problému. Pokud se však výpočet provede pro body z celé ČR, tak jsem narazil na limit knihovny R, jelikož není možné alokovat vektor větší než 3, nebo 4 GB. Proto jsem musel přistoupit k rozdělení výpočtu na několik menších částí, které budou následně spojeny v jeden grid. Jsem v kontaktu s panem Ježkem, kterého informuji o aktuálních výsledcích. Jakmile tedy bude nový grid hotový, tak budu pana Ježka informovat a následně bude dán k dispozici veřejnosti." MK "Martin Landa" píše v diskusním příspěvku news:CA+Ei1OdYW5Cuyx2jjgzFTu4=03TxOB0JAwO2QMgGQ0G5q=j...@mail.gmail.com... Ahoj, Dne 17. května 2014 14:06 Jachym Cepicky napsal(a): Martine, ty máš nějakou zkušenost s formátem pro nat2bin? Já našel jenom tohle Honza Jezek mi psal, ze ma studenta, ktery pracuje na novem gridu, doufam, ze se dozvime vic na Geoinformatics [1]... Martin [1] http://geoinformatics.fsv.cvut.cz/gwiki/Geoinformatics_FCE_CTU_2014 ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Nový S-JTSK grid, kdo pomůže?
Ohledně té chyby http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Grid, díval jsem se na ten XLS a vzal bod s největší odchylkou 930102291 http://dataz.cuzk.cz/hledej.php?a=1 > podle čísla bodu > číslo TL = 3010, číslo bodu = 229 tak vyjede zrušený bod, bohužel při "hromadnějším" exportu do texťáku podle mapových listů je informace o zrušeném bodu ztracena, proto v tom mém txt na Google Drive tato informace není. MK "Martin Landa" píše v diskusním příspěvku news:CA+Ei1Ofzp2Bm_1MbKtQrgtZrfKp=DfWo2yXEVk1FVSfx6=a...@mail.gmail.com... Dne 19. května 2014 16:26 hanoj napsal(a): wget http://gis.templ.net/grid_jezek2008/jezek_czech08.llb sudo cp jezek_czech08.llb /usr/share/proj toto je jediny funkcni link s gridem pokud vim, poznamenal jsem to na wiki [1]. Martin [1] http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Grid#Transformace_S-JTSK_-_WGS84_.28ETRS.29_pomoc.C3.AD_gridu -- Martin Landa * http://geo.fsv.cvut.cz/gwiki/Landa ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
