[GRASS-de] v.proj

Hartmut Wziontek hwz at gmx.de
Sa Mai 15 21:01:36 CEST 2004


hallo,

hat jemand erfahrungen mit v.proj?

beim versuch, gk-koordinaten (dhdn, 5. streifen) in 
laenge/breite (wgs84) zu konvertieren, erhalte ich 
ziemlichen unsinn, obwohl die transformationsparameter,
die v.proj aus den projektdaten ermittelt meiner meinung 
nach korrekt sind: 

**************************************************************
GRASS:~/grassdata > v.proj input=ps-20040515-t2 location=30e_baskup mapset=gps out=z7

Input Projection Parameters: +proj=tmerc +lat_0=0.0000000000 +lon_0=15.0000000000 +k_0=1.0000000000 +x_0=0.0000000000 +y_0=5500000.0000000000 +a=6377397.155 +es=0.006674372231802145 +towgs84=612.4,77.0,440.2,-0.054,0.057,-2.797,2.55

Output Projection Parameters: +proj=latlong +a=6378137 +es=0.006694379990141316 +towgs84=0.000,0.000,0.000

Creating dig file...
north:      48S
south:      50S
east:       15E
west:       13E
**************************************************************

die punkte liegen jedoch auf der noerdlichen hemisphaere, im bereich 51.1N/14.7E :(

interessant ist, dass ich mit 

m.proj2 inproj="proj=tmerc,datum=potsdam, \
                lat_0=0.0,lon_0=15.0,k=1.0,  \
                x_0=5500000.0" \
        outproj="proj=ll,datum=wgs84" \
        inproj="proj=tmerc,datum=potsdam, \
                lat_0=0.0,lon_0=15.0,k=1.0,  \
                x_0=5500000.0" \

mit punktdaten (ascii-xy) sinnvolle ergebnisse erhalte. da 
sowohl in m.proj als auch v.proj auf der PROJ.4-lib basieren, 
kann ich mir das alles nciht recht erklaeren. 
was mache ich fuer einen fehler bei der benutzung von v.proj?
funktioniert v.proj ueberhaupt korrekt?

viele gruesse

hartmut