<html><body>Odpovím si sám ;-)<br><br><p>---------- Původní zpráva ----------<br>Od: Marián Kyral <mkyral@email.cz><br>Komu: talk-cz@openstreetmap.org<br>Datum: 24. 5. 2014 23:54:56<br>Předmět: Re: [Talk-cz] pLPIS a WFS - jak na to?</p><br><blockquote>Dne 16.5.2014 01:27, hanoj napsal(a):<br>>> Tak jsem měl teď přes oběd trochu času a koukal jsem na WFS.<br>>><br>>> Základní URL: ....<br>>><br>>> Prvních 200 položek: ....<br>>><br>>> No a s tím mám jen dva (dost podstatné ;-) ) problémy<br>> *** a což zkusit project ID?<br>><br>> 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<br><br>To právě musím nejprve někde zjistit ne? Mám akorát bod, na který bylo<br>kliknuto.<br><br>> nebo BBOX?<br>><br>> 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<br>><br><br>Trochu jsem si s tím hrál. Našel jsem nějaký příklad, jak pomocí<br>geotools převádět mezi projekcemi. Nějak to funguje, ale nevím, jestli<br>dobře.<br><br>Například:<br>proN49.660144848959035, E18.29907177958237 (EPSG:4326)<br>dostanu 471210.46759085247, 1121224.5497619768 (EPSG:102067)<br><br>V tom URL nahoře jsou ale všechny souřadnice BBOXU záporné. Tady ne. Je<br>to dobře?<br>Když se třeba kouknu, co za request posílá LPIS, tak tam to je taky<br>záporné. Navíc je vidět, že vypočítané souřadnice by cca souhlasily.<br>(neklikl jsem na stejné místo).<br><br>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<br></blockquote><p><br></p><p>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 <b>102067 alias EPSG:</b>5514 alias S-JTSK / Krovak East North [1]<br></p><p><br></p><p>Takže pro převod lze použít cs2cs [2]<br></p><p><br></p><p>No a ten pro N49.660144848959035, E18.29907177958237 vrátí:</p><p><br></p><p>$ echo "18.29907177958237 49.660144848959035" |cs2cs -f "%.2f" +proj=longlat +datum=WGS84 +to +proj=krovak +ellps=bessel +nadgrids=czech<br>-471084.50 -1121102.96 0.00<br></p><p><br></p><p>Což vypadá slibně.</p><p><br></p><p>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)<br></p><p><br></p><p><br></p>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.<br><p>No vypadá to nadějně :-D<br></p><br>Marián<br><p><br></p><p>[1]
http://geoportal.cuzk.cz/%28S%28vfttsk45dozkiyjz5itt15ew%29%29/Default.aspx?mode=TextMeta&side=sit.trans&text=souradsystemy<br></p><p>[2] http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Grid</p><p>[3] http://josm.openstreetmap.de/wiki/Help/Plugin/UtilsPlugin2<br></p><p>[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</p><p><br></p><br></body></html>