Skip to Content.
Sympa Menu

freegeocz - [FreeGeoCZ] v.proj a JTSK

freegeocz AT fsv.cvut.cz

Předmět: Svobodná geoinformační infrastruktura

List archive

[FreeGeoCZ] v.proj a JTSK


Chronological Podle vláken 
  • From: Rastislav Vojtko <vojtko AT fns.uniba.sk>
  • To: freegeocz AT fsv.cvut.cz
  • Subject: [FreeGeoCZ] v.proj a JTSK
  • Date: Thu, 28 Jun 2012 14:18:15 +0200
  • List-archive: <http://mailman.fsv.cvut.cz/pipermail/freegeocz>
  • List-id: Svobodná geoinformační infrastruktura <freegeocz.fsv.cvut.cz>

Ahojte,

mam jeden problem, ktory by som potreboval vyriesit cim skor. Mam subor bodovych udajov, ktore su v latitude/longitude, ktory je definovany:
GRASS 6.4.2:~ > g.region -p
projection: 3 (Latitude-Longitude)
zone: 0
datum: ** unknown (default: WGS84) **
ellipsoid: wgs84
north: 49:08:02N
south: 47:37:42N
west: 16:42:48E
east: 19:16:08E
nsres: 0:00:01
ewres: 0:00:01
rows: 5420
cols: 9200
cells: 49864000
GRASS 6.4.2:~ > g.proj -p
-PROJ_INFO-------------------------------------------------
name : Latitude-Longitude
proj : ll
ellps : wgs84
-PROJ_UNITS------------------------------------------------
unit : degree
units : degrees
meters : 1.0

Tieto udaje potrebujem previest do SJTSK, ktory je definovany nasledovne:
GRASS 6.4.2:~ > g.region -p
projection: 99 (Krovak)
zone: 0
datum: hermannskogel
ellipsoid: bessel
north: -1226641.89403
south: -1230792.75052
west: -369713.492349
east: -364565.345475
nsres: 7.5060696
ewres: 7.50458728
rows: 553
cols: 686
cells: 379358
GRASS 6.4.2:~ > g.proj -p
-PROJ_INFO-------------------------------------------------
name : Krovak
datum : hermannskogel
towgs84 : 485.021,169.465,483.839,7.786342,4.397554,4.102655,0
proj : krovak
ellps : bessel
a : 6377397.1550000003
es : 0.0066743722
f : 299.1528128000
-PROJ_UNITS------------------------------------------------
unit : meter
units : meters
meters : 1.0

Po pouziti prikazu v.proj mi udaje ulozi o cca 130 m JZ od realnych bodov, kde by mali padnut. Ked vsak tieto bodove data prefuknem cez cs2cs padnu presne tam kde maju. Napr.:
echo "19.8569480516017 48.7236205860973" | cs2cs +proj=latlong +datum=WGS84 \
+to +proj=krovak +ellps=bessel +no_defs \
+towgs84=485.021,169.465,483.839,7.786342,4.397554,4.102655,0

Samozrejme, ze je to alternativa, ale strasna makacka to robit priamo v terminali pomocou cs2cs. Keby to fungovalo korektne v GRASSe bolo by to best. Viete mi poradit, kde je nieco zhnite? Kde robim chybu? Okrem toho neviem si predstavit, ako by som robil ine zlozitejsie objekty ako su body.

Pekny den,

Rasto




  • [FreeGeoCZ] v.proj a JTSK, Rastislav Vojtko, 06/28/2012

Archivace běží na MHonArc 2.6.19+.

Top of Page