GeorgeOuzou

@dimamp: Γενικώς ναι, 3-μεταθέσεις, με την διαφορά ότι όταν στον μετασχηματισμό συμμετέχει σύστημα αναφοράς του Παλαιού Ελληνικού Datum ή του HTRS07 γινεται ένας ενδιάμεσος μετασχηματισμός σε ΕΓΣΑ87/ΤΜ87 χρησιμοποιώντας τις ακριβέστερες μεθόδους (ΟΚΧΕ πολυώνυμα, μετασχηματισμός Hepos με διορθωτικά grid).

Δηλαδή π.χ. για τον μετασχηματισμό ED50/UTM34 --> ΤΜ3 Δυτ.Ζώνη το πρόγραμμα αυτόματα κάνει τα εξής βήματα:

• ED50/UTM34 --> ΕΓΣΑ87/ΤΜ87 (μέσω 3-παραμετρικού)
• ΕΓΣΑ87/ΤΜ87 --> Hatt (μέσω πολυωνύμων ΟΚΧΕ)
• Hatt --> TM3

Βέβαια σε αυτή την περίπτωση δεν κερδίζουμε ή χάνουμε κάτι, ο 3-παραμετρικός μετασχηματισμός είναι ούτως ή άλλως ακρίβειας 5 μέτρων.

Όσον αφορά τα λ, φ στο ίδιο datum εφαρμόζει τις αντίστροφες προβολικές σχέσεις μόνο. Δηλαδή o "μετασχηματισμός" ΕΓΣΑ87/ΤΜ87 --> ΕΓΣΑ87(λ,φ) είναι ουσιαστικά μετατροπή συντεταγμένων με τις αντίστροφες προβολικές σχέσεις.

@dpapdpap:  Δυστυχώς δεν μπορώ να προσδιορίσω σε τι μπορεί να οφείλεται αυτή η διαφορά. Με έναν πρώτο έλεγχο στο εξής site: http://cs2cs.mygeodata.eu/ βλέπω ότι βγάζω τα ίδια αποτελέσματα βάζοντας τα εξής δεδομένα:

• Input proj4 text: +proj=tmerc +lat_0=0 +lon_0=24 +k=0.9996 +x_0=500000 +y_0=0 +ellps=GRS80 +towgs84=-199.723,74.030,246.018 +units=m
• Output proj4 text: +proj=longlat +ellps=WGS84 +datum=WGS84
• Input coordinate pairs: 412860.350997563,4511507.66899286 

@dpapdpap: Για να μπορέσω να ελέγξω το πρόβλημα, είναι εύκολο να μου δώσετε εκτός από τις output, τις input συντεταγμένες σας απο ΕΓΣΑ87/ΤΜ87 σε WGS84? Επίσης για τους μεμονομένους μετασχηματισμούς WGS84 --> ΕΓΣΑ87(λ,φ) και  ΕΓΣΑ87(λ,φ) --> ΕΓΣΑ87/ΤΜ87 θα ήταν εφικτό να ελέγξουμε αν το πρόβλημα βρίσκεται στον 3-παραμετρικό μετασχηματισμό ή στην προβολή.

Ευχαριστώ πολύ για τα καλά σας λόγια.

@dpapdpap: Για οποιοδήποτε μετασχηματισμό από ή προς Hatt (αλλά και γενικότερα το Παλαιό Ελληνικό Datum) χρησιμοποιώ έναν ενδιάμεσο μετασχηματισμό σε ΕΓΣΑ87 μέσω των πολυωνυμικών συντελεστών του ΟΚΧΕ. Στην περίπτωση του μετασχηματισμού Hatt-->WGS84 γίνεται πρώτα ένα ενδιάμεσο βήμα Hatt --> ΕΓΣΑ87/ΤΜ87 και μέτα ένα επιπλέον βήμα ΕΓΣΑ87/ΤΜ87 --> WGS84 με ΔΧ=-199.723, ΔΥ=74.030, ΔZ=246.018 (Γεωμετρική Γεωδαισία του Α. Φωτίου)

Θα ήθελα να σας παρουσιάσω μία νέα, web εφαρμογή μετασχηματισμών.

Παρέχεται δωρεάν και μπορείτε να την επισκεφθείτε στην ιστοσελίδα http://georgeouzou.pythonanywhere.com

Ευελπιστώ να βοηθήσει τους συναδέλφους μηχανικούς. Η εφαρμογή θα συνεχίσει να αναπτύσσεται και είμαι ανοιχτός σε προτεινόμενες βελτιώσεις.