_GEO_CH1903_WGS84 ( y ; x )

Convert Swiss coordinate data to normal WGS84 coordinates you find with GPS devices.

Average rating: 3.7 (41 votes) Log in to vote

Ralph Nusser - Show more from this author
Sogetes Computer-Services
http://www.sogetes.com

Sample input:

_GEO_CH1903_WGS84 ( 746320 ; 254651 )

Sample output:

47.42639 , 9.37806

Function definition: (Copy & paste into FileMaker's Edit Custom Function window)

/*

Syntax:
_GEO_CH1903_WGS84 ( y ; x )

Author / Autor:
Ralph Nusser, 13. Mai 2006

FileMaker example download / FileMaker Beispiel herunterladen:
http://www.filemaker-magazin.de/fmmforum/60198.html

Example Input / Beispiel Eingabe: cf_CH1903_WGS84 (745689;253816)
Example Output / Beispiel Ausgabe: [WGS84-Koordinaten (λ, φ)] ( 47.41903, 9.36944)

Source / Quellen:
http://www.swisstopo.ch/pub/down/basics/geo/system/ch1903_wgs84_de.pdf
http://www.swisstopo.ch/pub/down/basics/geo/system/swiss_projection_de.pdf
http://de.wikipedia.org/wiki/CH1903

Explanation / Erklärung gemäss swisstopo.ch-Aufsatz:
Näherungsformeln für die direkte Umrechnung von: Schweizer Projektionskoordinaten (y [Rechtswert], x [Hochwert]) ellipsoidischen WGS84-Koordinaten (λ, φ)

Genauigkeit:
im 0.1"-Bereich
Diese Näherungen sind für die ganze Schweiz besser als 0.12" in der Länge, 0.08" in der Breite und 0.5 Meter in der Höhe.
Es handelt sich dabei um eine Herleitung von U. Marti vom Mai 1999, basierend auf den Formeln aus [Bolliger, 1967]

Erklärung zur Annäherungsformel:

1. Die Projektionskoordinaten y (Rechtswert) und x (Hochwert) sind ins zivile System (Bern = 0 / 0) und in die Einheit [1000 km] umzuwandeln:

     y' = (y – 600000 m)/1000000
     x' = (x – 200000 m)/1000000

2. Länge und Breite in der Einheit [10000"] berechnen:

     λ λ' = 2.6779094
     + 4.728982 * y'
     + 0.791484 * y' * x'
     + 0.1306 * y' * x'2
      - 0.0436 * y'3

     φ φ' = 16.9023892
     + 3.238272 * x'
     - 0.270978 * y' 2
     - 0.002528 * x' 2
     - 0.0447 * y' 2 * x'
     - 0.0140 * x' 3

     h [m] = h' + 49.55
     - 12.60 * y'
     - 22.64 * x'

3. Umrechnen der Länge und Breite in die Einheit [°]

     λ = λ' * 100 / 36
     φ = φ' *100 / 36

4. Zahlenbeispiel

     gegeben: y = 700 000 m x = 100 000 m h' = 600 m
     => y' = 0.1 x' = -0.1
     => λ' = 3.14297976 φ' = 16.57588564 h = 650.55 m
     => λ = 8° 43' 49.80" φ = 46° 02' 38.86"
     aus NAVREF: λ = 8° 43' 49.79" φ = 46° 02' 38.87" h = 650.60 m [bei swisstopo.ch können Werte online umgerechnet werden, um die Genauigkeit zu prüfen]



Bemerkung zu den Höhen:
In diesen Formeln wird davon ausgegangen, dass mit ellipsoidischen Höhen gearbeitet wird, wie sie z.B. mit GPS-Messungen erhalten werden. Wird mit 'Höhen über Meer' gearbeitet, so sind die Höhen im Meterbereich in beiden Systemen gleich. Sie müssen also in diesem Fall nicht umgerechnet werden.


*/


Let ( [ 

// VARIABLENDEKLARATION

y' = (y - 600000) / 1000000;
x' = (x - 200000) / 1000000;

λ' = 2.6779094  + 4.728982 * y' + .791484 * y' * x' + .1306 * y' * x' ^ 2  - .0436 * y' ^ 3;
φ' = 16.9023892  + 3.238272 * x' - .270978 * y' ^ 2  - .002528 * x' ^ 2  - .0447 * y' ^ 2 * x' - .014 * x' ^ 3;

Lng1 = λ' * 100 / 36;
Lat1 = φ' * 100 / 36

]; // End Let

Case ( 
IsEmpty ( y );"";
IsEmpty ( x );"";

Lat1 & " , " & Lng1

) // End Case
) // End Let

see Wikipedia article:
http://de.wikipedia.org/wiki/CH1903
Download the example file at the German "FileMaker Magazin": http://www.filemaker-magazin.de/fmmforum/60198.html

Comments

Note: these functions are not guaranteed or supported by BrianDunning.com. Please contact the individual developer with any questions or problems.

_GEO_CH1903_WGS84 ( y ; x )

Comments

Support this website.

Support this website.