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

Odpovedet emailem