[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