A script for solving elevation networks in geodesy with Least Squares Method, but not for profesional use as it may contain bugs
Input --> Ortho Heights in .txt
Output --> Corrected Ortho Heights in .txt and info about the Network Quality
Requires Numpy and FPDF librarys to run
- Solves every 1 Dimentional Elevation Netowrk wthich consists of Ortho Heights
- Outputs a .txt with the Corrected Ortho Heights
- Generates A PDF file with the report of the Network Quality which contains every points height with its error (from covariance matrix), Sigma 0 value and Residuals of every height
Input Data should be ortho height of the points, and its distances between points. See demo_data.txt for reference. Every input data should look like this
r8-r9 2.534
r8-r14 5.722
r7-r9 1.416
r7-r14 4.649
r9-r10 1.567
r9-r11 2.647
r9-r12 5.573
r9-r14 3.188
Static Point is defined inside the code in line 167 and 168, its defined by tha value and name
After height value, data can be blank like above or set a value for distance between 2 points
If you like to suggest anything on this scirpt please contact me in nikolas619065@gmail.com
Το παραπάνω πρόγραμμα έχει γραφτεί σε Python 3.9 και μπορεί να χρησιμοποιηθεί για την συνόρθωση ενός υψομετρικού δικτιού. Δεν είναι για επαγγελματική χρήση καθώς μπορεί να περιέχει bugs τα οποία δεν έχουν ακόμα λυθεί.
Τα δεδομένα πρέπει να είναι σε αρχείο .txt με την μορφή όπως φαίνεται στο αρχείο demo_data.txt οπως το παραδειγμα απο κατω (ονοματα σημειων, υψομετρικη διαφορα, αποσταη μεταξυ σημειων--> αυτο μπορει να παραλειφθει)
π.χ r1-r2 2.938 1.2
Το πρόγραμμα λειτουργεί με μονάδες μετρήσεις τα μετρά. Το πρόγραμμα κάνει την συνόρθωση με την μέθοδο των ελάχιστων εξωτερικών δεσμεύσεων και για αυτό το λόγο μέσα στο πρόγραμμα πρέπει να δηλωθεί πιο σημείο είναι το σταθερό και το υψόμετρο του όπως και την τοποθεσία του αρχείου που είναι τα δεδομένα. H αβεβαιότητά μεταξύ των μετρήσεων ειναι βάση απόστασης μεταξύ κορυφών το οποίο πρέπει να δηλωθεί στα δεδομένα μετά το υψόμετρο σε km. Το πρόγραμμα από default θα έχει πίνακα βαρών (P) μοναδιαίο το οποιο μπορει να αλλαξει απο το org_data βαζοντας να περνει τιμες το s_dist.
Τελικα το προγραμμα επιστρεφει ενα αρχειο output.txt με μονο τα υψομετρα μεσα, ενα αρχειο report.txt με διαφορα δεδομενα της συνορθωσης και τελικα ενα .pdf αρχειο οπου εχει μεσα αναλυτικα πληροφοριες για την συνορθωση οπως το σο οπου εχει μοναδες μετρα και τα υπολοιπα του πινακα δλ για την ερευση καποιου προβληματικου δεδομενου.
Απαραιτητα Libraries ειναι τα NumPy και FPDF