Odpovím si sám ;-)
---------- Původní zpráva ---------- Od: Marián Kyral <mky...@email.cz> 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.6333312988281&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
