Koordinatentransformation (Geodäsie/Vermessung)

Mark Johnson ⌂ @, Berlin, Thursday, 22.04.2010, 05:28 (vor 5118 Tagen) @ Surkel

So, ich habe mit Hilfe eines Kollegen das kleine Umrechnungstool fertiggestellt. Es funktioniert bis auf einen halben Meter genau.

Ich wäre darin intressiert diese Tool zu bekommen.

Gestern habe mit eint tool names 'cs2cs' gearbeitet:
[http://trac.osgeo.org/proj/]

cs2cs -E -f "%.15f" -r -s -v +init=epsg:3068 +to +init=epsg:4326 input.txt

# ---- From Coordinate System ----
#Cassini
# Cyl, Sph&Ell
# +init=epsg:3068 +proj=cass +lat_0=52.41864827777778
# +lon_0=13.62720366666667 +x_0=40000 +y_0=10000 +ellps=bessel
# +datum=potsdam +units=m +no_defs +towgs84=606.0,23.0,413.0
#--- following specified but NOT used
# +ellps=bessel
# ---- To Coordinate System ----
#Lat/long (Geodetic alias)
# 
# +init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
# +towgs84=0,0,0
#--- following specified but NOT used
# +ellps=WGS84
10000 40000 52.417257037334295 13.625433714571523 -8.342795641189696
19760.743 20281.021 52.504612931236132 13.335030335112300 -8.212803853439254 
21292.525 25254.895 52.518536505762093 13.408211248735720 -8.342374544203722
21047.0017702 23181.9915716 52.516269670837112 13.377686876448641 -8.297927901715866

Diese Ergebnisse kommen die von transdat sehe nah.
In der source-code kann man ansatzweise sehen wie dies berechnet wird.
(PJ_cass.c aus der projekt)

Die Parameter basieren auf die Angabe die man unter 'http://spatialreference.org/ref/epsg/3068/' finden kann.
Weitere Information über Cassini-Soldner findet man :

http://www.remotesensing.org/geotiff/proj_list/cassini_soldner.html

Mark Johnson
Berlin
mj10777@googlemail.com


gesamter Thread:

 RSS-Feed dieser Diskussion