[FreeGeoCZ] v.proj a JTSK

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 mailing list
[email protected]
http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz