Μετατροπή συντεταγμένων από HATT σε ΕΓΣΑ

[stavrosLinux]

Υπάρχουν on-line εφαρμογές για την μετατροπή, κάποιες freeware.
π.χ. http://www.meleth.gr/Topo.html

http://www.geographer.gr/ci/index.php/hatt/

Μήπως γνωρίζει καποιος αν υπάρχει εφαρμογή για linux?

Επίσης υπάρχει το coord_gr http://dasologoi.gr/xrhsima-eggrafa/viewdownload/16-genika-eggrafa-ergaleia-kai-arxeia/127-coordgr.html. Σαν άδεια αναφέρει Creative Commons (το έχω κατεβάσει και θέλω να το δοκιμάσω σε wine).

σχετικά με τα μαθηματικά που χρειάζονται: http://www.geoapikonisis.gr/index.php?option=com_content&view=article&id=160&Itemid=115&lang=el

Κώδικας: Επιλογή όλων

Μετασχηματισμοί συντεταγμένων

Για να μετατραπούν οι συνεταγμένες που αναφέρονται σε ένα προβολικό σύστημα σε ένα άλλο, πρέπει να ληφθούν υπόψην:

    η αλλαγή του datum
    η αλλαγή του δικτύου
    η αλλαγή της προβολής

Για τις μετατροπές αυτές χρησιμοποιούνται μαθηματικοί τύποι, οι οποίοι στις περισσότερες περιπτώσεις είναι πολύπλοκοι και καθιστούν την χρήση υπολογιστικού μέσου απαραίτητη.

Μετατροπή από το προβολικό σύστημα HATT στο ΕΓΣΑ'87

Με την υιοθέτηση του νέου γεωδαιτικού συστήματος αναφοράς ΕΓΣΑ'87 και με την προοπτική έναρξης του προγράμματος του Εθνικού Κτηματολογίου, έγινε αναγκαία η ύπαρξη ενός εύχρηστου εργαλείου μετατροπής των συντεταγμένων από το παλαιό προβολικό σύστημα ΗΑΤΤ. Για να αποθευχθεί η χρήση πολύπλοκων μαθηματικών τύπων που απαιτούνται λαμβάνοντας υπόψην και τις τρείς αλλαγές (του δικτύου, του datum και της προβολής), αποφασίστηκε από το Οργανισμό Κτηματολογίου και Χαρτογραφήσεων Ελλάδας (Ο.Κ.Χ.Ε) να γίνει εργασία προσδιορισμού Συντελεστών Μετατροπής Συντεταγμένων από το σύστημα ΗΑΤΤ στο ΕΓΣΑ'87. Η εργασία αυτή πραγματοποιήθηκε από την Γεωγραφική Υπηρεσία Στρατού (μετά από αίτηση και χρηματοδότηση του ΟΚΧΕ). Οι σχέσεις μετατροπής των συντεταγμένων (που καθορίστηκαν σε συνεργασία του ΟΚΧΕ, της ΓΥΣ και του εργαστηρίου Ανώτερης Γεωδαισίας του Εθνικού Μετσοβίου Πολυτεχνείου) είναι πολυώνυμα 2ου βαθμού της μορφής:

X = Ao + A1 x + A2 y + A3 x2 + A4 y2 + A5 x y

Y = Bo + B1 x + B2 y + B3 x2 + B4 y2 + B5 x y

όπου:

    (X,Y) οι ζητούμενες συντεταγμένες στο προβολικό σύστημα ΕΓΣΑ'87.
    (x,y) οι γνωστές συντεταγμένες προς μετατροπή στο προβολικό σύστημα ΗΑΤΤ.
    Αo...A5, Bo...B5 είναι οι πολυωνυμικοί συντελεστές μετατροπής οι οποίοι, γενικά, είναι υπολογισμένοι ανά φύλλο χάρτη κλίμακας 1:100.000 και σε μερικές περιπτώσεις για αυξημένη ακρίβεια, είναι υπολογισμένοι ανά φύλλο χάρτη κλίμακας 1:50.000 (και σε ειδικές περιπτώσεις και για τμήματα αυτών). Η Γεωγραφική Υπηρεσία Στρατού έχει εκδόσει έντυπο με τους παραπάνω Συντελεστές Μετατροπής.


[anagno]

Για χαρτογραφικές μετατροπές υπάρχει και το proj. Τώρα αν υποστηρίζει την μετατροπή από HATT σε ΕΓΣΑ δεν ξέρω αλλά τί χρειάζεσαι κιόλας αυτή την μετατροπή ??? Τί θες να κάνεις ??? Ως φοιτητής στην Σ.Α.Τ.Μ. δεν έχω χρησιμοποιήσει HATT ποτέ ...
Πάντως να έχεις υπόψιν πως όλες αυτές οι μετατροπές εισάγουν σημαντικά σφάλματα και γι` αυτό καλύτερα θα ήταν να πεις τί στοιχεία έχεις και πως θέλεις να τα χρησιμοποιήσεις ... Αν θυμάμαι καλά οι μετατροπές αυτές είναι κατάλληλες μόνο χαρτογραφικούς σκοπούς και όχι για μετρήσεις ακριβείας.
Πάντως αν δεν βγάλεις άκρη γράψου και στο φόρουμ των τοπογράφων και φαντάζομαι ότι θα πάρεις κάποια απάντηση

[simosx]

Δες στο http://simos.info/blog/archives/1102 για λογισμικό σε Linux.
Στην παραπάνω δημοσίευση περιγράφω πως θα μπορούσε κάποιος να φτιάξει ένα πρόγραμμα σε Python (π.χ. pyGTK) όπου να μπορεί να γίνει εύκολα η μετατροπή δεδομένων.

[stavrosLinux]

Για χαρτογραφικές μετατροπές υπάρχει και το proj. Τώρα αν υποστηρίζει την μετατροπή από HATT σε ΕΓΣΑ δεν ξέρω αλλά τί χρειάζεσαι κιόλας αυτή την μετατροπή ??? Τί θες να κάνεις ??? Ως φοιτητής στην Σ.Α.Τ.Μ. δεν έχω χρησιμοποιήσει HATT ποτέ ...
Πάντως να έχεις υπόψιν πως όλες αυτές οι μετατροπές εισάγουν σημαντικά σφάλματα και γι` αυτό καλύτερα θα ήταν να πεις τί στοιχεία έχεις και πως θέλεις να τα χρησιμοποιήσεις ... Αν θυμάμαι καλά οι μετατροπές αυτές είναι κατάλληλες μόνο χαρτογραφικούς σκοπούς και όχι για μετρήσεις ακριβείας.
Πάντως αν δεν βγάλεις άκρη γράψου και στο φόρουμ των τοπογράφων και φαντάζομαι ότι θα πάρεις κάποια απάντηση

Πρέπει να γνωρίζεις ότι υπάρχουν πολλά παλιά δημόσια (εν χρήση) τοπογραφικά διαγράμματα σε HATT. Έτσι η ανάγκη για τη μετατροπή υπάρχει ακόμη. Η ακρίβεια είναι σχετικό πράγμα.

Αν είναι γραμμένος στο topografoi.com θα μπορούσες ανοίξεις ένα θέμα με παραπομπή στο παρόν.

@simox δεν βλέπω αναφορά για το σύστημα HATT, θα ρίξω μια ματιά αν αναφέρει κάτι.

[anagno]

Αν είναι γραμμένος στο topografoi.com θα μπορούσες ανοίξεις ένα θέμα με παραπομπή στο παρόν.

Άνοιξα το θέμα . Άμα προκύψει κάτι θα ενημερώσω ...

Πρέπει να γνωρίζεις ότι υπάρχουν πολλά παλιά δημόσια (εν χρήση) τοπογραφικά διαγράμματα σε HATT. Έτσι η ανάγκη για τη μετατροπή υπάρχει ακόμη. Η ακρίβεια είναι σχετικό πράγμα.

Το γνωρίζω ότι υπάρχουν αρκετά παλιά τοπογραφικά διαγράμματα αλλά συνήθως ο τρόπος αντιμετώπισης τους είναι διαφορετικός για την επίτευξη ακριβειών που να συνάδουν σε τοπογραφικά διαγράμματα.
Οι σχέσεις (που αναφέρεις πιο πάνω) αυτές εν γένει δεν ενδείκνυνται για γεωδαιτικούς σκοπούς γιατί οι αποκλίσεις μπορεί να είναι σημαντικές και ίσως να ξεπερνούν και το 1m. Βέβαια είναι ενδεχόμενο σε κάποια περιοχή να έχει γίνει καλή προσαρμογή και οι αποκλίσεις να είναι μικρές πιθανά και λίγα εκατοστά. Οπότε στην πράξη, πρέπει να ελέγχει κανείς τοπικά, στην συγκεκριμένη περιοχή που πρόκειται να χρησιμοποιήσει τον μετασχηματισμό αυτό, αν τον ικανοποιεί από άποψη ακρίβειας, Ο έλεγχος αυτός μπορεί να γίνει με εφαρμογή των σχέσεων για τις συντεταγμένες των τριγωνομετρικών της περιοχής με βάση της επίσημες τιμές της ΓΥΣ και στα δύο συστήματα. (Μετατροπές συστημάτων Αναφοράς και Προβολικών συστημάτων - Υψομετρικά Συστήματα Αναφοράς, Π. Μήλας, Χ. Μπιλλήρης, Εκδόσεις Ε.Μ.Π., Αθήνα 2006 ,σελ. 17)

Γι` αυτό σε ρώτησα για ποιο σκοπό χρειάζεσαι τις μετατροπές αυτές, για την αποφυγή τυχόν λαθών ...

[pmav99]

Για τα bindings του proj σε python ένα minimum working example είναι αυτό.

Κώδικας: Επιλογή όλων

from pyproj import Proj
import numpy as np

to_WGS84 = Proj(init='epsg:4326')
to_EGSA87 = Proj(init='epsg:2100')

x_old = np.array((23.72826, 37.98414))
y_old = np.array((37.98414, 37.98414))

x_new, y_new = to_EGSA87(x_old, y_old)

for x, y in zip(x_new, y_new):
    print x, y


[simosx]

@simox δεν βλέπω αναφορά για το σύστημα HATT, θα ρίξω μια ματιά αν αναφέρει κάτι.

Η αναφορά για αυτό το HATT από την αγγλική Wikipedia είναι ελλιπής, http://en.wikipedia.org/wiki/Hellenic_G ... ystem_1987
Θα ήταν εξαιρετικό αν το πρόγραμμα proj (πακέτο proj-bin) υποστήριζε και HATT ώστε να είναι εφικτή η μετατροπή μέσω των Python bindings.
Δες στο http://trac.osgeo.org/proj/ για το πως/αν μπορείς να ρυθμίσεις το πακέτο proj ώστε να υποστηρίξει και το ελληνικό grid.

[medigeek]

(το έχω κατεβάσει και θέλω να το δοκιμάσω σε wine).

http://www.geographer.gr/gis/69--coordgr-ubuntu.html
Επίσης βρήκα αυτό: http://itia.ntua.gr/~soulman/icoordstrans/ -- http://itia.ntua.gr/~soulman/icoordstrans/ccode.zip

[medigeek]

Επίσης βρήκα αυτό: http://itia.ntua.gr/~soulman/icoordstrans/ -- http://itia.ntua.gr/~soulman/icoordstrans/ccode.zip

Θέλει compile με -fPIC, οπόταν το Makefile:

Κώδικας: Επιλογή όλων

all: icoordstrans test

icoordstrans: main.o main_types.o corefuncs.o
   g++ -fPIC -shared -O -o libicoordstrans.so main.o main_types.o corefuncs.o -lm

test: test.o
   g++ -fPIC -L . -O -o test test.o -licoordstrans

main.o: main.cpp
   g++ -fPIC -c main.cpp

main_types.o: main_types.cpp
   g++ -fPIC -c main_types.cpp

corefuncs.o: corefuncs.cpp
   g++ -fPIC -c corefuncs.cpp

test.o: test.c
   g++ -fPIC -c test.c

clean:
   rm -f corefuncs.o main_types.o main.o icoordstrans.so test.o test


Θα το κοιτάξω την επόμενη βδομάδα, δεν δουλεύει για κάποιο λόγο:
$ ./test
./test: error while loading shared libraries: libicoordstrans.so: cannot open shared object file: No such file or directory

[simosx]

Θα το κοιτάξω την επόμενη βδομάδα, δεν δουλεύει για κάποιο λόγο:
$ ./test
./test: error while loading shared libraries: libicoordstrans.so: cannot open shared object file: No such file or directory

Η βιβλιοθήκη libicoordstrans.so πρέπει να πάει στο /usr/local/lib και να τρέξεις μετά sudo ldconfig (ή επανεκκίνηση).