pLPIS a WFS - jak na to?

25 zpráv
Zpět na přehled

pLPIS a WFS - jak na to?

25 zpráv MJHVPM 6 účastníků 29 min čtení
  1. Marián Kyral mkyral na email.cz #mdabb61
    Tak jsem měl teď přes oběd trochu času a koukal jsem na WFS. Základní URL: http://eagri.cz/public/app/wms/plpis_wfs.fcgi?SERVICE=WFS&REQUEST=GetCapabilities Prvních 200 položek: http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4 No a s tím mám jen dva (dost podstatné ;-) ) problémy: 1) Souřadnice: V odpovědi vidím EPSG:102067 a nepodařilo se mi to nijak změnit. Na SRS=EPSG:4326 nebo DefaultSRS=EPSG:4326 to nijak nereaguje. Na srsName=EPSG:4326 zase vrátí: "Invalid GetFeature Request: SRSNAME value should be valid for all the TYPENAMES". Znamená to tedy, že si to musím přepočítat sám? Jak? Je na to nějaká knihovna pro javu? 2) Jak mám ten dotaz poskládat, abych zadal bod a server mi vrátil objekt , který je kolem toho bodu? Musí tam být BBOX? Koukal jsem, že by to mělo jít nějak přes filtr, ale zase. moc moudrý z toho nejsem. Díky, Marián
    Ahoj, s tím Tracerem2 mě to už také napadlo, ale z nějakého důvodu to nefunguje. Co se týče licence, tak se to už tady někdy řešilo. Podle příručky pro ovládání http://eagri.cz/public/web/file/2091/_120420_PZ_2012_03_UG_pLPIS.pdf se jedná o "veřejný seznam z evidence využití půdy" - str. 4. A pokud jde o stažení dat pře email, tak na str. 17 je následující text "Export je volně přístupný z prostředí internetu. Aplikace exportují veřejná data, právo použít tento export má každý uživatel aplikace." Zdraví Pavel Kwiecien Na trasování tohohle wms můžeš zkusit Tracer2 (potažmo původní Tracer - TracerServer). Ty umí z principu tracovat jakoukoli bitmapu. Osobně bych se tomu rád vyhnul. Ale dle všeho nemusí být jediný výstup jako obrázek. WMS podporuje i xml. I WFS zdá se funguje a vypadá nadějně. http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4 Umíš udělat nějaký dotaz, kde stačí zadat bod v 4326 a odpověď dostanu taky tak?? Zkoumám dokumentaci WFS, ale stále z toho moc nechápu :-( Myslím, že pak už zůstane jediná otázka, licence. Marián Dne 17.2.2014 18:30, Pavel Kwiecien napsal: Ahoj, asi by nebyl problém udělat tracer wms mapy -- nastavení pro josm wms:http://eagri.cz/public/app/wms/plpis.fcgi?FORMAT=image/jpeg&VERSION=1.1.1&SERVICE=WMS&REQUEST=GetMap&LAYERS=LPIS_FB4_KOD&STYLES=&SRS={proj}&WIDTH={width}&HEIGHT={height}&BBOX={bbox} V nějakém manuálu také píšou, že lze data stáhnout přes WFS, to jsem ale nikdy nezkoušel, takže nevím, jestli to opravdu jde. Rozhodně nemá smysl ručně stahovat asi 12000 katastrální území. Zdraví Pavel Kwiecien A někdy v budoucnu - trasování souvislých ploch (pole jedním kliknutím) Ahoj, pokud by šlo o pole, tak tady bych nepoužil RUIAN, ale znovu bych oprášil myšlenku importovat pLPIS (veřejný registr půdy).
    No problém. Dej mi server, který mi pro zadaný bod vrátí polygon a jeho vlastnosti a zbytek už v pluginu dodělám. Bohužel co jsem koukal, je potřeba si o data žádat ručně a to po jednotlivých katastrálních územích :-( Nebo je i jiná možnost, jak se k těm datům dostat? Ač data v RUIANu nejsou nejpřesnější, mají tu výhodu, že jsou jednoduše dostupná. Bylo by zajímavé udělat nějaké srovnání RUIAN a pLPIS. Marián Links: ------
  2. Vojtěch Kalčík vojta na kalcik.cz #m064574
    Pro přepočet mezi souřadnými systémy existuje tato javovská knihovna http://www.jhlabs.com/java/maps/proj/index.html neumí ale úplně všechno jak Proj.4, tak nevím, jestli tam projde definice S-JTSK. GeoTools http://geotools.org/ by to mělo asi umět taky. Na stránkách Proj.4 je seznam alternativnívh knihoven: http://trac.osgeo.org/proj/.
  3. Marián Kyral mkyral na email.cz #m1aec82
    Pr(ed c(asem jsem si s tím trochu hrál. Bohuz(el z WMS nedostanu vektor, ale jen a pouze bitmapu. Takz(e jediná cesta je pr(es tracování bitmapy, tedy me(lo by stac(it správne( nastavit tracer. Nicméne(, proc( to sem pís(i. V konferenci josm-dev se pr(ed c(asem objevil dotaz na pr(ístup k bitmapovým datu*m za úc(elem trasování. Viz: https://lists.openstreetmap.org/pipermail/josm-dev/2014-April/006940.html Pokud to podar(í dokonc(it, tak bude stac(it vybrat si WMS/TMS vrstvu a pak uz( jen klikat ;-) Tak drz(me palce. Marián
  4. hanoj ehanoj na gmail.com #mb418e4
    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_FB4.8780373 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 ha hanoj
  5. Jachym Cepicky jachym.cepicky na gmail.com #m2a79c8
    Omlouvám se, že takhle vpadám do diskuze - moc jsem nečel celé vlákno, ale nepoužívejte WMS, když je dostupnej primární datový zdroj - WFS. to je jako dělat screenshot z monitoru počítače foťákem ručně Pokud server na eagri nepodporuje epsg:4326 tak nejsprávnější přístup je zvednout telefon a požádat je, aby tam epsg:4326 přihodili - je to jeden řádek konfigurace v mapserveru (na kterym to běží). na transformace v javě fungujou pokud já vím geotools - třeba by pomohl nějak rychle Michal Kepka (cced)? Já to "umím" v C světě (python, shell script, ...) WFS může vrátit prostě "všechny prvky" nebo pouze vyfiltrované (těžko říct, co je na straně serveru podporováno) Každopádně standard WFS 1.1.0: www.opengeospatial.org/standards/wfs příklady pro getfeature request jsou na straně 98 (a dalších) Ten request níže se mi bohužel nepodařilo rozjet, na serveru jim něco hnije (Server exception 500) zkusím se zeptat, kdo má to wfs na starosti - metadata té služby nejsou vyplněná (achjo). Můžu se zeptat, kde jste ten link na to WFS vyštrachali? Ještě jednou doporučení: radši chvíli (pár dní) počkejte. Tracovat data z WMS je špatně, pokud je dostupný primární zdroj. Včera na konferenci CAGI padala strašná spousta dotazů na OSM, bylo cítit jak se ten projekt líbí, že má podporu v odborné veřejnosti. Ale jeden dotaz byl pořád ve vzduchu: jaká je přesnost těch dat. A já bych byl opravdu rád, kdybych mohl říkat, že tam jde došlo k importu, jsou ta data tak přesná, jako primární zdroj J
  6. Marián Kyral mkyral na email.cz #m75c4bd
    Omlouvám se, že takhle vpadám do diskuze - moc jsem nečel celé vlákno, ale nepoužívejte WMS, když je dostupnej primární datový zdroj - WFS. to je jako dělat screenshot z monitoru počítače foťákem ručně
    V pohodě. Já zatím s ničím nezačal. O dané oblasti nic moc nevím a zatím jsem se v ní nedokázal pořádně zorientovat.
    Pokud server na eagri nepodporuje epsg:4326 tak nejsprávnější přístup je zvednout telefon a požádat je, aby tam epsg:4326 přihodili - je to jeden řádek konfigurace v mapserveru (na kterym to běží).
    Jednak nevím, že to obnáší jen změnu konfigurace a taky nejsem moc komunikativní typ (ve vztahu k neznámým lidem/institucím). Tak trochu jsem počítal s pomocí, že mi třeba někdo pomůže sestavit vhodný dotaz a já už pak dodělám podporu do Pointinfa a Traceru. Bohužel odezva minimální (nějak to zapadlo).
    na transformace v javě fungujou pokud já vím geotools - třeba by pomohl nějak rychle Michal Kepka (cced)? Já to "umím" v C světě (python, shell script, ...)
    No já v javě taky neuměl, než jsem si začal v lednu hrát s Tracerem :-D Geotools zkusím, koukám, že se už v JOSM používají.
    WFS může vrátit prostě "všechny prvky" nebo pouze vyfiltrované (těžko říct, co je na straně serveru podporováno)
    Líbilo by se mi, kdybych poslal dotaz se souřadnicemi bodu a server by mi vrátil data o pozemku, jehož jsou tyto souřadnice součástí.
    Každopádně standard WFS 1.1.0: www.opengeospatial.org/standards/wfs příklady pro getfeature request jsou na straně 98 (a dalších)
    Mrknu, nějakou dokumentaci jsem zkoumal, ale bylo toho strašně moc a nedokázal jsem najít, co vlastně potřebuji.
    Ten request níže se mi bohužel nepodařilo rozjet, na serveru jim něco hnije (Server exception 500)
    Je třeba to zkoušet. Taky mi to nejprve nefungovalo, ale napodruhé už to jelo. Možná nějaké výpadky.
    zkusím se zeptat, kdo má to wfs na starosti - metadata té služby nejsou vyplněná (achjo).
    Díky
    Můžu se zeptat, kde jste ten link na to WFS vyštrachali?
    Už si přesně nepamatuji. Pravděpodobně jsem se nějak doklikal na tuto stránku: http://eagri.cz/public/web/mze/farmar/LPIS/novinky/nove-wms-a-wfs-sluzby.html Původní návrh na to využít data LPIS padl tady v konferenci. Myslím, že o to hodně stál Pavel Kwiecien.
    Ještě jednou doporučení: radši chvíli (pár dní) počkejte. Tracovat data z WMS je špatně, pokud je dostupný primární zdroj.
    No on ten obecný tracer hned tak nebude (pokud vůbec někdy bude).
    Včera na konferenci CAGI padala strašná spousta dotazů na OSM, bylo cítit jak se ten projekt líbí, že má podporu v odborné veřejnosti. Ale jeden dotaz byl pořád ve vzduchu: jaká je přesnost těch dat. A já bych byl opravdu rád, kdybych mohl říkat, že tam jde došlo k importu, jsou ta data tak přesná, jako primární zdroj
    No, ono s tou přesností to je ošidné. OSM nikdy nebude tak přesné jako přesně zaměřené KM. Ani ten export není stoprocentní, pokud to nějaký automat nebude striktně kontrolovat a opravovat nepřesnosti a omyly vznikající při manuální editaci. Pak ovšem musíme danému zdroji věřit na sto procent. Takhle se třeba pracuje s hranicemi. Ale obecně v OSM importům moc nevěří a po předchozích špatných zkušenostech je snaha importy omezit na nejmenší možnou míru a podmínky pro importy jsou docela drsné. Hlavně diskuze na osm-imports ;-) No a u ostatních dat hodně záleží na přesnosti podkladů (bing, km, gps). Tam to může docela hodně lítat. Klidně i o desítky centimetrů. Někdy i pár metrů - viz třeba ulice, které po nahrání budov z RUIAN procházejí skrz ty budovy. Docela by mně zajímalo, k čemu by chtěla odborná veřejnost data z OSM využít, že je zajímá přesnost ;-) Marián
  7. Jachym Cepicky jachym.cepicky na gmail.com #m2483e2
    Dorazíte na konferenci? Mohli bychom zkusit někoho zlomit aby přišel říct svůj pohled na OSM, někoho z čuzaku zbytek stále zjištuju J
  8. Petr Souček soucekp na atlas.cz #m044e44
    Dobrý den, s účastí zástupců ČÚZK můžete na 99% počítat. Prezentace a jejich obsah domluvíme podle aktuálního zájmu. Petr Souček
  9. Marián Kyral mkyral na email.cz #ma8cbdf
    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_FB4.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.6333312988281&height=783&width=1258&info_format=application%2Fvnd.ogc.gml Ale, I když toto pominu a dám souřadnice záporné, tak tápu, jak správně definovat ten BBOX. Jediné co udělám, že v URL výše vyměním souřadnice. Vezmu ten bod, napoprvé odečtu 100, napodruhé přičtu a dám to do URL.¨ http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4&BBOX=-471110,-1121124,-471310,-1121324&SRSNAME=EPSG:102067 Místo očekávané odpovědi ale dostanu xml, kde je klíč gml:null s hodnotou "missing": http://pastebin.com/zFXfmPSe A očekávané odpovědi se nedočkám ani když použiji souřadnice z URL od LPIS. Zase klíč s hodnotou "missing". http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4&BBOX=-471131.845377%2C-1121174.062338%2C-470688.051172%2C-1120897.837487&SRSNAME=EPSG:102067 Co mi uniká? Marián
  10. Marián Kyral mkyral na email.cz #m751db0
    Odpovím si sám ;-)
  11. hanoj ehanoj na gmail.com #m87ce83
    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
  12. Marián Kyral mkyral na email.cz #m8b2fef
    "> 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
  13. Jachym Cepicky jachym.cepicky na gmail.com #mce99b1
    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
  14. Martin Kokes shr3k na typo3-hosting.com #me1980b
    Shodou okolností jsem bývalým zaměstnancem Sitewellu a na základech LPISu jsem také před 12 lety dělal a Sitewell si u mne dodnes outsourcuje nějaké služby. Není šance, že by to kdokoliv ze Sitewellu opravil, to spíš tlač na MZE. Sitewell byl jen subdodavatel pro Telefonicu a kvůli právním sporům ohledně provozovaných licencí byla veškerá aktiva týkající se této věci převedena na společnost Tescium viz zápis v OR. MK "Jachym Cepicky" píše v diskusním příspěvku news:20140529065701.GD23593 na doctor... "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)."
  15. Marián Kyral mkyral na email.cz #mff196a
    "Použij buď geotools (java) nebo gdal + knihovnu proj (C++) je to jedno, důležitý je výsledek. Geotools jsem právě zkoušel, ale zadal jsem tam EPSG:102067 a to nějak nebylo ono ;-) Takže myslíš, že když tam zadám EPSG:5514 tak to bude fungovat rovnou bez jakýchkoli dalšíck kejklí? "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). No snad se ještě ozvou. Naděje umírá poslední. 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ů. Nebo vezmu jako referenční souřadnice, které jsem tady testoval. O víkendu vyzkouším. Zvláštní, že v JOSM ještě není nějaká třída, která by geotools zapouzdřila (nebo je a já ji přehlédl). Ale každý plugin si to dělá po svém. Včera jsem se snažil z JOSM API vytáhnout BBOX aktuální screeny. Zatím se nezadařilo. Bádání pokračuje ;-) 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. Pokud by Petr Vejsada neměl nic proti, tak bychom mohli použít poloha.net. Nevím jak to je teď s vytížeností serveru a jaké má podmínky (datový limit). Pokud by bylo potřeba, rád přispěji na provoz. Případně se inspirovat u kolegů na Slovensku [1]. Nemůžu si pomoct, ale ty jejich mapy [2] se mi moc líbí. Marián [1] http://oz.freemap.sk/ [2] http://www.freemap.sk/ Jachym
    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
  16. Marián Kyral mkyral na email.cz #m7785b0
    Ahoj, Netušíte někdo, jaký je rozdíl mezi LPIS_FB4 a LPIS_FB4_BBOX? S LPIS_FB4 mi to vždy vrátí "missing". Teď jsem objevil LPIS_FB4_BBOX a to zase končí na krásné ORA chybě níže :-( <FeatureType> <Name>LPIS_FB4</Name> <Title>UcinneFB-obrys</Title> <Abstract>Půdní bloky - účinné (max. 200 záznamů)</Abstract> <DefaultSRS>EPSG:102067</DefaultSRS> <OutputFormats> <Format>text/xml; subtype=gml/3.1.1</Format> </OutputFormats> <ows:WGS84BoundingBox dimensions="2"> <ows:LowerCorner>11.6287033481963 47.8921430739317</ows:LowerCorner> <ows:UpperCorner>19.2039743450146 51.5189483888239</ows:UpperCorner> </ows:WGS84BoundingBox> </FeatureType> <FeatureType> <Name>LPIS_FB4_BBOX</Name> <Title>UcinneFB-obrys</Title> <Abstract>Půdní bloky - účinné (max. 200 záznamů)</Abstract> <DefaultSRS>EPSG:102067</DefaultSRS> <OutputFormats> <Format>text/xml; subtype=gml/3.1.1</Format> </OutputFormats> <ows:WGS84BoundingBox dimensions="2"> <ows:LowerCorner>11.6287033481963 47.8921430739317</ows:LowerCorner> <ows:UpperCorner>19.2039743450146 51.5189483888239</ows:UpperCorner> </ows:WGS84BoundingBox> </FeatureType> $ BBOX="&bbox=-470952.15,-1121084.54,-471309.44,-1120790.40" $ wget -O - "http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4_BBOX&${BBOX}&SRSNAME=EPSG:102067" --2014-05-30 08:49:09--  http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4_BBOX&&bbox=-470952.15,-1121084.54,-471309.44,-1120790.40&SRSNAME=EPSG:102067 Překládám eagri.cz? 94.199.42.131 Navazuje se spojení s eagri.cz|94.199.42.131|:80? spojeno. HTTP požadavek odeslán, program čeká na odpověď? 200 OK Délka: neudáno [text/xml] Ukládám do: ?STDOUT? [<=> ] 0           --.-K/s              <?xml version="1.0" encoding="UTF-8"?> <ows:ExceptionReport xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ows="http://www.opengis.net/ows" version="1.1.0" language="en-US" xsi:schemaLocation="http://www.opengis.net/ows http://schemas.opengis.net/ows/1.0.0/owsExceptionReport.xsd"> <ows:Exception exceptionCode="mapserv" locator="NoApplicableCode"> <ows:ExceptionText>msWFSGetFeature(): WFS server error. ms_error-&gt; code not found msOracleSpatialLayerWhichShapes(): OracleSpatial error. Error: ORA-29902: při vykonávání rutiny ODCIIndexStart() se vyskytla chyba ORA-13208: vnitřní chyba při vyhodnocování operátoru [window SRID does not match layer SRID] ORA-06512: na "MDSYS.SDO_INDEX_METHOD_10I", line 333 . Query statement: SELECT rownum, UZIVATEL, ADRESA_UZIVATELE, ID_UZ, ICO, ID_FB, CTVEREC, ZKODFB, STAV, PLATNYOD, PLATNYDO, VYMERAM, KULTZAKNAZEV, REZIMEZ, ZPUSOBILOST, PRISLUSZA, VYSKA, SVAZITOST, KULTURAOD, MAPOVYLIST, MAPOVYLIST5, LHUTAOD, STAV_KOD, VYMZCHU, VYMZCHU1, VYMMKRAS, VYMBKARP, VYMMCHU, VYMLFA_HA, VYMLFA_HB, VYMLFA_OA, VYMLFA_OB, VYMLFA_S, VYMLFA_SX, VYMLFA_E, PHOVYMER1, KODCILE, VYMNAT_2000, VYMNAT_PTACI, VYMNAT_VYZNLOK, ZAUDIT0, GEOMETRY FROM VIEW_MG_LPIS_RFB_WFS WHERE  rownum &lt; 51  AND SDO_ FILTER( GEOMETRY, MDSYS.SDO_GEOMETRY(2003, :srid, NULL,MDSYS.SDO_ELEM_INFO_ ARRAY(1,1003,3),:ordinates ),'querytype=window') = 'TRUE' . Check your data statement.</ows:ExceptionText> </ows:Exception> </ows:ExceptionReport> [ <=> ] 1 478       --.-K/s   za 0s 2014-05-30 08:49:10 (53,4 MB/s) ? zapsáno na standardní výstup [1478] Marián
  17. Jachym Cepicky jachym.cepicky na gmail.com #m94f7f8
    Zkuil jsem nakopnout tu proxy, https://gist.github.com/jachym/8c33b98a794679bc5de3 Ale tedy že by to nějak zvlášť fungovalo, to nefunguje :-( myslím, že jim tam něco hnije, ale jistě to říct si netroufám ještě se tomu zkusím pověnovat j
  18. Marián Kyral mkyral na email.cz #mfd55b3
    Souřadnice mi teď těžkou hlavu nedělají. Když už jsem si ujasnil, co je co, tak se to nějak udělá. Spíše potřebuji nějak rozjet ten dotaz přes bbox. Co jsem našel, tak by měl bohatě stačit nějaký miniaturní bbox. Ale prostě to nefunguje. Buď dělám něco blbě, nebo je problém u nich. Marián
  19. Marián Kyral mkyral na email.cz #m7dc723
    Takz(e jedna dobrá zpráva. Odepsali mi z helpdesku mze, z(e WFS opravili. Dále jsem se dove(de(l, z(e pro BBOX dotazy se musí pouz(ívat LPIS_FB4_BBOX. Vrstva LPIS_FB4 je zase urc(ena na dotazy dle ID pu*dního bloku. Takz(e jsou potr(eba dva dotazy: Jeden pro získání geometrie ID pu*dního bloku. http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4_BBOX&bbox=-471162.327712,-1120637.409667,-471161.327712,-1120636.409667&SRSNAME=EPSG:102067 A druhý pro dopln(ující detaily: http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4&&featureID=LPIS_FB4.9630656&SRSNAME=EPSG:102067 Takz(e ted( uz( jen trochu c(asu, vyr(es(it pr(evod sour(adnic a ove(r(it licenci (nechce se toho ne(kdo ujmout, at( tam nepís(i jen já? :-D ) a mu*z(u to zaintegrovat do Traceru. Marián
  20. Jachym Cepicky jachym.cepicky na gmail.com #mcfa7b1
    S převodem souřadnic můžu pomoct, ale asi to vyhoří na platformě (funguju na python/c, ne java) popis pro jst je popsán tady: http://docs.geotools.org/stable/userguide/library/api/jts.html#transform Nevím, jestli máte ve vaší verzi správně zadefinovanýho Křováka - EPSG:102067 nebo EPSG:5514 (podle verze geotools) Na http://epsg.io/102067 a http://epsg.io/5514 se dají sehnat definice pro GeoServer (který by měl být s geotools kompatibilní. Zejména 5514 vypadá nadějně: 5514=PROJCS["S-JTSK / Krovak East North",GEOGCS["S-JTSK",DATUM["System_Jednotne_Trigonometricke_Site_Katastralni",SPHEROID["Bessel 1841",6377397.155,299.1528128,AUTHORITY["EPSG","7004"]],TOWGS84[589,76,480,0,0,0,0],AUTHORITY["EPSG","6156"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4156"]],PROJECTION["Krovak"],PARAMETER["latitude_of_center",49.5],PARAMETER["longitude_of_center",24.83333333333333],PARAMETER["azimuth",30.28813972222222],PARAMETER["pseudo_standard_parallel_1",78.5],PARAMETER["scale_factor",0.9999],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["X",EAST],AXIS["Y",NORTH],AUTHORITY["EPSG","5514"]] Doporučuju použít tam jiný parametr TOWGS84, a to [570.8,85.7,462.8,4.998,1.587,5.261,3.56] převod bude přesnější ------------ Co to znamená "ověřit licenci" ? Zavolat na ministerstvo? J
    Takže jedna dobrá zpráva. Odepsali mi z helpdesku mze, že WFS opravili. Dále jsem se dověděl, že pro BBOX dotazy se musí používat LPIS_FB4_BBOX. Vrstva LPIS_FB4 je zase určena na dotazy dle ID půdního bloku. Takže jsou potřeba dva dotazy: Jeden pro získání geometrie ID půdního bloku. [1]http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4_BBOX&bbox=-471162.327712,-1120637.409667,-471161.327712,-1120636.409667&SRSNAME=EPSG:102067 A druhý pro doplňující detaily: [2]http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4&&featureID=LPIS_FB4.9630656&SRSNAME=EPSG:102067 Takže teď už jen trochu času, vyřešit převod souřadnic a ověřit licenci (nechce se toho někdo ujmout, ať tam nepíši jen já? :-D ) a můžu to zaintegrovat do Traceru. Marián Souřadnice mi teď těžkou hlavu nedělají. Když už jsem si ujasnil, co je co, tak se to nějak udělá. Spíše potřebuji nějak rozjet ten dotaz přes bbox. Co jsem našel, tak by měl bohatě stačit nějaký miniaturní bbox. Ale prostě to nefunguje. Buď dělám něco blbě, nebo je problém u nich. Marián Zkuil jsem nakopnout tu proxy, [5]https://gist.github.com/jachym/8c33b98a794679bc5de3 Ale tedy že by to nějak zvlášť fungovalo, to nefunguje :-( myslím, že jim tam něco hnije, ale jistě to říct si netroufám ještě se tomu zkusím pověnovat j
    Ahoj, Netušíte někdo, jaký je rozdíl mezi LPIS_FB4 a LPIS_FB4_BBOX? S LPIS_FB4 mi to vždy vrátí "missing". Teď jsem objevil LPIS_FB4_BBOX a to zase končí na krásné ORA chybě níže :-( <FeatureType>   <Name>LPIS_FB4</Name>   <Title>UcinneFB-obrys</Title>   <Abstract>Půdní bloky - účinné (max. 200 záznamů)</Abstract> >   <DefaultSRS>EPSG:102067</DefaultSRS>   <OutputFormats>     <Format>text/xml; subtype=gml/3.1.1</Format>   </OutputFormats>   <ows:WGS84BoundingBox dimensions="2">     <ows:LowerCorner>11.6287033481963
    47.8921430739317</ows:LowerCorner>
        <ows:UpperCorner>19.2039743450146
    51.5189483888239</ows:UpperCorner>
      </ows:WGS84BoundingBox> </FeatureType> <FeatureType>   <Name>LPIS_FB4_BBOX</Name>   <Title>UcinneFB-obrys</Title>   <Abstract>Půdní bloky - účinné (max. 200 záznamů)</Abstract> >   <DefaultSRS>EPSG:102067</DefaultSRS>   <OutputFormats>     <Format>text/xml; subtype=gml/3.1.1</Format>   </OutputFormats>   <ows:WGS84BoundingBox dimensions="2">     <ows:LowerCorner>11.6287033481963
    47.8921430739317</ows:LowerCorner>
        <ows:UpperCorner>19.2039743450146
    51.5189483888239</ows:UpperCorner>
      </ows:WGS84BoundingBox> </FeatureType> $ BBOX="&bbox=-470952.15,-1121084.54,-471309.44,-1120790.40" > $ wget -O -
    [6]"http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4_BBOX&${BBOX}&SRSNAME=EPSG:102067" > --2014-05-30 08:49:09-- [7]http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4_BBOX&&bbox=-470952.15,-1121084.54,-471309.44,-1120790.40&SRSNAME=EPSG:102067 > Překládám eagri.cz? 94.199.42.131
    Navazuje se spojení s eagri.cz|94.199.42.131|:80? spojeno. HTTP požadavek odeslán, program čeká na odpověď? 200 OK Délka: neudáno [text/xml] Ukládám do: ?STDOUT?
    [<=>                                                                                                                                                                   &! nbsp;  ;
    ] 0           --.-K/s              <?xml version="1.0"
    encoding="UTF-8"?>
    <ows:ExceptionReport
    xmlns:xsi=[8]"http://www.w3.org/2001/XMLSchema-instance"
    xmlns:ows=[9]"http://www.opengis.net/ows" version="1.1.0"
    language="en-US"
    xsi:schemaLocation="[10]http://www.opengis.net/ows [11]http://schemas.opengis.net/ows/1.0.0/owsExceptionReport.xsd"> >   <ows:Exception exceptionCode="mapserv" locator="NoApplicableCode"> >     <ows:ExceptionText>msWFSGetFeature(): WFS server error. ms_error-&gt;code not found msOracleSpatialLayerWhichShapes(): OracleSpatial error. Error: ORA-29902: při vykonávání rutiny ODCIIndexStart() se vyskytla chyba ORA-13208: vnitřní chyba při vyhodnocování operátoru [window SRID does not match layer SRID] ORA-06512: na "MDSYS.SDO_INDEX_METHOD_10I", line 333  . Query statement: SELECT rownum, UZIVATEL, ADRESA_UZIVATELE, ID_UZ, ICO, ID_FB, CTVEREC, ZKODFB, STAV, PLATNYOD, PLATNYDO, VYMERAM,
    KULTZAKNAZEV,
    REZIMEZ, ZPUSOBILOST, PRISLUSZA, VYSKA, SVAZITOST, KULTURAOD, MAPOVYLIST, MAPOVYLIST5, LHUTAOD, STAV_KOD, VYMZCHU, VYMZCHU1, VYMMKRAS, VYMBKARP, VYMMCHU, VYMLFA_HA, VYMLFA_HB, VYMLFA_OA, VYMLFA_OB, VYMLFA_S, VYMLFA_SX, VYMLFA_E, PHOVYMER1, KODCILE, VYMNAT_2000, VYMNAT_PTACI,
    VYMNAT_VYZNLOK,
    ZAUDIT0, GEOMETRY FROM VIEW_MG_LPIS_RFB_WFS WHERE  rownum &lt; 51 AND SDO_FILTER( GEOMETRY, MDSYS.SDO_GEOMETRY(2003, :srid, NULL,MDSYS.SDO_ELEM_INFO_ARRAY(1,1003,3),:ordinates
    ),'querytype=window')
    = 'TRUE' . Check your data statement.</ows:ExceptionText>   </ows:Exception> </ows:ExceptionReport>     [
    <=>                                                                                                                                                                   &n! bsp;
    ] 1 478       --.-K/s   za 0s 2014-05-30 08:49:10 (53,4 MB/s) ? zapsáno na standardní výstup [1478] Marián _______________________________________________
    -- Jachym Cepicky URL: http://les-ejk.cz e-mail: jachym.cepicky at gmail com PGP: http://les-ejk.cz/pgp/JachymCepicky.pgp @jachymc
  21. Marián Kyral mkyral na email.cz #mc4767e
    S převodem souřadnic můžu pomoct, ale asi to vyhoří na platformě (funguju na python/c, ne java)
    No já umím s pythonem zhruba stejně jako s javou. Něco jsem v tom napsal, ale ani jeden jazyk nepoužívám denně. Ale to není problém, většina jazyků si je podobná. Všude jsou proměnné a cykly ;-)
    popis pro jst je popsán tady: http://docs.geotools.org/stable/userguide/library/api/jts.html#transform Nevím, jestli máte ve vaší verzi správně zadefinovanýho Křováka - EPSG:102067 nebo EPSG:5514 (podle verze geotools) Na http://epsg.io/102067 a http://epsg.io/5514 se dají sehnat definice pro GeoServer (který by měl být s geotools kompatibilní. Zejména 5514 vypadá nadějně: 5514=PROJCS["S-JTSK / Krovak East North",GEOGCS["S-JTSK",DATUM["System_Jednotne_Trigonometricke_Site_Katastralni",SPHEROID["Bessel 1841",6377397.155,299.1528128,AUTHORITY["EPSG","7004"]],TOWGS84[589,76,480,0,0,0,0],AUTHORITY["EPSG","6156"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4156"]],PROJECTION["Krovak"],PARAMETER["latitude_of_center",49.5],PARAMETER["longitude_of_center",24.83333333333333],PARAMETER["azimuth",30.28813972222222],PARAMETER["pseudo_standard_parallel_1",78.5],PARAMETER["scale_factor",0.9999],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["X",EAST],AXIS["Y",NORTH],AUTHORITY["EPSG","5514"]] Doporučuju použít tam jiný parametr TOWGS84, a to [570.8,85.7,462.8,4.998,1.587,5.261,3.56] převod bude přesnější
    Díky za popis. Pokusím se to pochopit. Ale možná bych ještě počkal, třeba to nebude potřeba. Nevím jak to je s tou žádostí o přidání WGS84 transformace (jestli teda nějaká žádost byla odeslána). Na tu žádost ohledně nefunkční WFS jsem dostal první odpověď až po cca třech týdnech s tím, že to bylo vyřešeno. Takže jistá možnost tu je. Jen neznám podrobnosti.
    ------------ Co to znamená "ověřit licenci" ? Zavolat na ministerstvo?
    Jestli tam někoho znáš ;-) Pokud se shodneme, že toto:
    ... Podle příručky pro ovládání http://eagri.cz/public/web/file/2091/_120420_PZ_2012_03_UG_pLPIS.pdf se jedná o "veřejný seznam z evidence využití půdy" - str. 4. A pokud jde o stažení dat pře email, tak na str. 17 je následující text "Export je volně přístupný z prostředí internetu. Aplikace exportují veřejná data, právo použít tento export má každý uživatel aplikace."... http://permalink.gmane.org/gmane.comp.gis.openstreetmap.region.cz/9331
    se vztahuje i na WFS a není to v rozporu s licencí OSM, tak mi to stačí. Ovšem v příručce k WMS a WFS službám [1] se mluví pouze o prohlížení:
    Ministerstvo zemědělství poskytuje bezplatný přístup ke svým vybraným prostorovým datům přes WMS a WFS služby. Data jsou poskytována pro území celé České republiky a jsou dostupná zdarma. Uživatelé mají takto přístup k datům MZE a mohou si je zobrazit přes URL adresu, či v GIS prohlížeči (např. Quantum GIS).
    [1] http://eagri.cz/public/web/mze/farmar/LPIS/uzivatelske-prirucky/wms-sluzby/prirucka-pro-pouziti-wms-a-wfs-sluzeb.html Marián
  22. Jachym Cepicky jachym.cepicky na gmail.com #m19e837
    Za mě je ta licence OK - public domain data J
    S převodem souřadnic můžu pomoct, ale asi to vyhoří na platformě (funguju na python/c, ne java) No já umím s pythonem zhruba stejně jako s javou. Něco jsem v tom napsal, ale ani jeden jazyk nepoužívám denně. Ale to není problém, většina jazyků si je podobná. Všude jsou proměnné a cykly ;-) popis pro jst je popsán tady: [1]http://docs.geotools.org/stable/userguide/library/api/jts.html#transform Nevím, jestli máte ve vaší verzi správně zadefinovanýho Křováka - EPSG:102067 nebo EPSG:5514 (podle verze geotools) Na [2]http://epsg.io/102067 a [3]http://epsg.io/5514 se dají sehnat definice pro GeoServer (který by měl být s geotools kompatibilní. Zejména 5514 vypadá nadějně: 5514=PROJCS["S-JTSK / Krovak East North",GEOGCS["S-JTSK",DATUM["System_Jednotne_Trigonometricke_Site_Katastralni",SPHEROID["Bessel 1841",6377397.155,299.1528128,AUTHORITY["EPSG","7004"]],TOWGS84[589,76,480,0,0,0,0],AUTHORITY["EPSG","6156"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4156"]],PROJECTION["Krovak"],PARAMETER["latitude_of_center",49.5],PARAMETER["longitude_of_center",24.83333333333333],PARAMETER["azimuth",30.28813972222222],PARAMETER["pseudo_standard_parallel_1",78.5],PARAMETER["scale_factor",0.9999],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["X",EAST],AXIS["Y",NORTH],AUTHORITY["EPSG","5514"]] Doporučuju použít tam jiný parametr TOWGS84, a to [570.8,85.7,462.8,4.998,1.587,5.261,3.56] převod bude přesnější Díky za popis. Pokusím se to pochopit. Ale možná bych ještě počkal, třeba to nebude potřeba. Nevím jak to je s tou žádostí o přidání WGS84 transformace (jestli teda nějaká žádost byla odeslána). Na tu žádost ohledně nefunkční WFS jsem dostal první odpověď až po cca třech týdnech s tím, že to bylo vyřešeno. Takže jistá možnost tu je. Jen neznám podrobnosti. ------------ Co to znamená "ověřit licenci" ? Zavolat na ministerstvo? Jestli tam někoho znáš ;-) Pokud se shodneme, že toto: ... Podle příručky pro ovládání [4]http://eagri.cz/public/web/file/2091/_120420_PZ_2012_03_UG_pLPIS.pdf se jedná o "veřejný seznam z evidence využití půdy" - str. 4. A pokud jde o stažení dat pře email, tak na str. 17 je následující text "Export je volně přístupný z prostředí internetu. Aplikace exportují veřejná data, právo použít tento export má každý uživatel aplikace."... [5]http://permalink.gmane.org/gmane.comp.gis.openstreetmap.region.cz/9331 se vztahuje i na WFS a není to v rozporu s licencí OSM, tak mi to stačí. Ovšem v příručce k WMS a WFS službám [1] se mluví pouze o prohlížení: Ministerstvo zemědělství poskytuje bezplatný přístup ke svým vybraným prostorovým datům přes WMS a WFS služby. Data jsou poskytována pro území celé České republiky a jsou dostupná zdarma. Uživatelé mají takto přístup k datům MZE a mohou si je zobrazit přes URL adresu, či v GIS prohlížeči (např. Quantum GIS). [1] [6]http://eagri.cz/public/web/mze/farmar/LPIS/uzivatelske-prirucky/wms-sluzby/prirucka-pro-pouziti-wms-a-wfs-sluzeb.html Marián J Takže jedna dobrá zpráva. Odepsali mi z helpdesku mze, že WFS opravili. Dále jsem se dověděl, že pro BBOX dotazy se musí používat LPIS_FB4_BBOX. Vrstva LPIS_FB4 je zase určena na dotazy dle ID půdního bloku. Takže jsou potřeba dva dotazy: Jeden pro získání geometrie ID půdního bloku. [1][7]http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4_BBOX&bbox=-471162.327712,-1120637.409667,-471161.327712,-1120636.409667&SRSNAME=EPSG:102067 A druhý pro doplňující detaily: [2][8]http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4&&featureID=LPIS_FB4.9630656&SRSNAME=EPSG:102067 Takže teď už jen trochu času, vyřešit převod souřadnic a ověřit licenci (nechce se toho někdo ujmout, ať tam nepíši jen já? :-D ) a můžu to zaintegrovat do Traceru. Marián Souřadnice mi teď těžkou hlavu nedělají. Když už jsem si ujasnil, co je co, tak se to nějak udělá. Spíše potřebuji nějak rozjet ten dotaz přes bbox. Co jsem našel, tak by měl bohatě stačit nějaký miniaturní bbox. Ale prostě to nefunguje. Buď dělám něco blbě, nebo je problém u nich. Marián Zkuil jsem nakopnout tu proxy, [5][11]https://gist.github.com/jachym/8c33b98a794679bc5de3 Ale tedy že by to nějak zvlášť fungovalo, to nefunguje :-( myslím, že jim tam něco hnije, ale jistě to říct si netroufám ještě se tomu zkusím pověnovat j
    Ahoj, Netušíte někdo, jaký je rozdíl mezi LPIS_FB4 a LPIS_FB4_BBOX? S LPIS_FB4 mi to vždy vrátí "missing". Teď jsem objevil LPIS_FB4_BBOX a to zase končí na krásné ORA chybě níže :-( <FeatureType>   <Name>LPIS_FB4</Name>   <Title>UcinneFB-obrys</Title>   <Abstract>Půdní bloky - účinné (max. 200 záznamů)</Abstract> >   <DefaultSRS>EPSG:102067</DefaultSRS>   <OutputFormats>     <Format>text/xml; subtype=gml/3.1.1</Format>   </OutputFormats>   <ows:WGS84BoundingBox dimensions="2">     <ows:LowerCorner>11.6287033481963
    47.8921430739317</ows:LowerCorner>
        <ows:UpperCorner>19.2039743450146
    51.5189483888239</ows:UpperCorner>
      </ows:WGS84BoundingBox> </FeatureType> <FeatureType>   <Name>LPIS_FB4_BBOX</Name>   <Title>UcinneFB-obrys</Title>   <Abstract>Půdní bloky - účinné (max. 200 záznamů)</Abstract> >   <DefaultSRS>EPSG:102067</DefaultSRS>   <OutputFormats>     <Format>text/xml; subtype=gml/3.1.1</Format>   </OutputFormats>   <ows:WGS84BoundingBox dimensions="2">     <ows:LowerCorner>11.6287033481963
    47.8921430739317</ows:LowerCorner>
        <ows:UpperCorner>19.2039743450146
    51.5189483888239</ows:UpperCorner>
      </ows:WGS84BoundingBox> </FeatureType> $ BBOX="&bbox=-470952.15,-1121084.54,-471309.44,-1120790.40" > $ wget -O -
    [6][12]"http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4_BBOX&${BBOX}&SRSNAME=EPSG:102067" > --2014-05-30 08:49:09-- [7][13]http://eagri.cz/public/app/wms/plpis_wfs.fcgi?VERSION=1.1.0&SERVICE=WFS&REQUEST=GetFeature&TYPENAME=LPIS_FB4_BBOX&&bbox=-470952.15,-1121084.54,-471309.44,-1120790.40&SRSNAME=EPSG:102067 > Překládám eagri.cz? 94.199.42.131
    Navazuje se spojení s eagri.cz|94.199.42.131|:80? spojeno. HTTP požadavek odeslán, program čeká na odpověď? 200 OK Délka: neudáno [text/xml] Ukládám do: ?STDOUT?
    [<=>                                                                                                                                                                   &! nbsp;  ;
    ] 0           --.-K/s              <?xml version="1.0"
    encoding="UTF-8"?>
    <ows:ExceptionReport
    xmlns:xsi=[8][14]"http://www.w3.org/2001/XMLSchema-instance"
    xmlns:ows=[9][15]"http://www.opengis.net/ows" version="1.1.0" language="en-US" xsi:schemaLocation="[10][16]http://www.opengis.net/ows [11][17]http://schemas.opengis.net/ows/1.0.0/owsExceptionReport.xsd"> >   <ows:Exception exceptionCode="mapserv" locator="NoApplicableCode"> >     <ows:ExceptionText>msWFSGetFeature(): WFS server error. ms_error-&gt;code not found msOracleSpatialLayerWhichShapes(): OracleSpatial error. Error: ORA-29902: při vykonávání rutiny ODCIIndexStart() se vyskytla chyba ORA-13208: vnitřní chyba při vyhodnocování operátoru [window SRID does not match layer SRID] ORA-06512: na "MDSYS.SDO_INDEX_METHOD_10I", line 333  . Query statement: SELECT rownum, UZIVATEL, ADRESA_UZIVATELE, ID_UZ, ICO, ID_FB, CTVEREC, ZKODFB, STAV, PLATNYOD, PLATNYDO, VYMERAM,
    KULTZAKNAZEV,
    REZIMEZ, ZPUSOBILOST, PRISLUSZA, VYSKA, SVAZITOST, KULTURAOD, MAPOVYLIST, MAPOVYLIST5, LHUTAOD, STAV_KOD, VYMZCHU, VYMZCHU1, VYMMKRAS, VYMBKARP, VYMMCHU, VYMLFA_HA, VYMLFA_HB, VYMLFA_OA, VYMLFA_OB, VYMLFA_S, VYMLFA_SX, VYMLFA_E, PHOVYMER1, KODCILE, VYMNAT_2000, VYMNAT_PTACI,
    VYMNAT_VYZNLOK,
    ZAUDIT0, GEOMETRY FROM VIEW_MG_LPIS_RFB_WFS WHERE  rownum &lt; 51 AND SDO_FILTER( GEOMETRY, MDSYS.SDO_GEOMETRY(2003, :srid, NULL,MDSYS.SDO_ELEM_INFO_ARRAY(1,1003,3),:ordinates
    ),'querytype=window')
    = 'TRUE' . Check your data statement.</ows:ExceptionText>   </ows:Exception> </ows:ExceptionReport>     [
    <=>                                                                                                                                                                   &n! bsp;
    ] 1 478       --.-K/s   za 0s 2014-05-30 08:49:10 (53,4 MB/s) ? zapsáno na standardní výstup [1478] Marián _______________________________________________
    -- Jachym Cepicky URL: http://les-ejk.cz e-mail: jachym.cepicky at gmail com PGP: http://les-ejk.cz/pgp/JachymCepicky.pgp @jachymc
  23. Marián Kyral mkyral na email.cz #mdb4232
    Za mě je ta licence OK - public domain data J
    Tak to je fajn. Díky.
  24. Marián Kyral mkyral na email.cz #m2a8edd
    Za me( je ta licence OK - public domain data J
    Tak to je fajn. Díky.
    Pr(idal jsem do zdroju*: http://wiki.openstreetmap.org/wiki/WikiProject_Czech_Republic/freemap#pLPIS_-_ve.C5.99ejn.C3.BD_registr_p.C5.AFdy Prosím o kontrolu a doplne(ní. Pr(i té pr(ílez(itosti jsem jes(te( doplnil Mapbox Satellite, upravil Bing (ani jeden odkaz na pokrytí uz( nefungoval) a pr(esunul Yahoo a Adresní body - C(ÚZK + MVC(R do Nefunkc(ních/zastaralých zdroj dat. A kdyz( uz( jsme u toho: C(SÚ UIR-ZSJ se jes(te( stále pouz(ívá, nebo byl taky nahrazen RUIANem? http://wiki.openstreetmap.org/wiki/WikiProject_Czech_Republic/freemap#.C4.8CS.C3.9A_UIR-ZSJ Marián
  25. Marián Kyral mkyral na email.cz #m4c6aa8
    "S převodem souřadnic můžu pomoct, ale asi to vyhoří na platformě (funguju na python/c, ne java) popis pro jst je popsán tady: http://docs.geotools.org/stable/userguide/library/api/jts.html#transform Nevím, jestli máte ve vaší verzi správně zadefinovanýho Křováka - EPSG: 102067 nebo EPSG:5514 (podle verze geotools) Na http://epsg.io/102067 a http://epsg.io/5514 se dají sehnat definice pro GeoServer (který by měl být s geotools kompatibilní. Zejména 5514 vypadá nadějně: 5514=PROJCS["S-JTSK / Krovak East North",GEOGCS["S-JTSK",DATUM["System_ Jednotne_Trigonometricke_Site_Katastralni",SPHEROID["Bessel 1841", 6377397.155,299.1528128,AUTHORITY["EPSG","7004"]],TOWGS84[589,76,480,0,0,0, 0],AUTHORITY["EPSG","6156"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]], UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY[ "EPSG","4156"]],PROJECTION["Krovak"],PARAMETER["latitude_of_center",49.5], PARAMETER["longitude_of_center",24.83333333333333],PARAMETER["azimuth", 30.28813972222222],PARAMETER["pseudo_standard_parallel_1",78.5],PARAMETER[ "scale_factor",0.9999],PARAMETER["false_easting",0],PARAMETER["false_ northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["X",EAST],AXIS[ "Y",NORTH],AUTHORITY["EPSG","5514"]] Doporučuju použít tam jiný parametr TOWGS84, a to [570.8,85.7,462.8,4.998,1.587,5.261,3.56] převod bude přesnější Tak jsem to s pomocí google poskládal dohromady a ono to funguje. Díky. Radost mi jen kazí zjištění, že ve WFS datech chybí dost podstatné pole: "kultura". Takže mám sice souřadnice polygonu, ale už nevím, jestli to je louka, orná půda nebo zalesněná plocha :-( Napsal jsem na helpdesk mze, jestli by to pole nemohli přidat. Doufám, že to nebude problém. Marián
Napsat odpověď e-mailem… Odpovědět

Otevře váš e-mailový klient. Odpovědi pak sledujte zde na webu.