Βαθμονόμηση στερεοκάμερας και χρήση της για την δημιουργία 3D φωτογραμμετρικού σαρωτή

1. Βαθμονόμηση στερεοκάμεραςκαι χρήση της για την δημιουργία3D φωτογραμμετρικού σαρωτή Πρόκος Αντώνης Δρ. Αγρονόμος και Τοπογράφος Μηχανικός

2. Βαθμονόμηση μηχανής Εύρεση Εσωτερικού προσανατολισμού: Σταθερά μηχανής (c) Συντεταγμένες πρωτεύοντος σημείου (xo, yo) Συντελεστές ακτινικής διαστροφής (k1, k2) …

3. Από εικόνες σκακιέρας

4. Fauccal Douskoset al. 2009

5. Από σημεία φυγής Γραμματικόπουλος, 2007

6. Αυτόματες διαδικασίες Μονοεικονικές

7. Στερεοκάμερα

8. Στερεοκάμερα Απαραίτητη η εύρεση: • εσωτερικού προσανατολισμού για τις δύο μηχανές ξεχωριστά • πλήρους σχετικού προσανατολισμού μεταξύ των μηχανών

9. Περιγραφή Βαθμονόμησης • Δύο μηχανές πακτωμένες σε ράβδο. (Σταθερός σχετικός προσανατολισμός.) • Εικόνες από επίπεδο πεδίο ελέγχου μορφής κανάβου «σκακιέρας». • Αυτόματη αναγνώριση σημείων. • Αυτόματη αναγνώριση αρχής συστήματος συντεταγμένων. (κόκκινο τετράγωνο) • Εύρεση αρχικών τιμών των αγνώστων. • Συνόρθωση παρατηρήσεων.

10. Βαθμονόμηση • Δεδομένα: • Ζευγάρια συγχρονισμένων εικόνων. (οι εικόνες δείχνουν ένα πεδίο ελέγχου με ταξινομημένα φωτοσταθερά σε κάνναβο) • Διάσταση κανάβου • Ζητούμενα: • Ξεχωριστός εσωτερικός προσανατολισμός για τις δύο μηχανές. • Πλήρης σχετικός προσανατολισμός. • Επίλυση με μέθοδο Δέσμης

11. Μέθοδος Δέσμης • Άγνωστοι: • xoyo c k1 k2(για την αριστερή μηχανή) • xoyo c k1 k2(για την δεξιά μηχανή) • Χο ΥοΖο ω φ κ (αριστερής μηχανής γιακάθε θέση) • Bx By Bzδω δφδκ • Παρατηρήσεις: • x – εικονοσυντεταγμένη Συνθήκη • y – εικονοσυντεταγμένη Συγγραμμικότητας

12. Εύρεση σημείων • Τελεστής Harris (Operator) • Συλλέγει σημεία ενδιαφέροντος με κριτήριο οι κλίσεις γύρω από το σημείο είναι μεγάλες.

13. Εύρεση Συστήματος Συντεταγμένων • Πρέπει να βρεθεί αυτόματα το κόκκινο τετράγωνο! • Γίνεται χωρισμός των καναλιών της εικόνας: • image_red • image_green • image_blue • Ανάδειξη κόκκινων της εικόνας: image_red - image_green Y X O Z

14. Εύρεση Συστήματος Συντεταγμένων • Τιμή Κατωφλίου (threshold) • Αύξηση εικόνας (dilate image) • Laplace • 4 σημεία. Καθορισμός συστήματος ΟΧΥΖ

15. Εύρεση όλων των σημείων • Υπάρχει καλή εκτίμηση του βήματος του κανάβου. • Κίνηση κατά τον y- άξονα • Κίνηση κατά τον x- άξονα • Όλα τα σημεία έχουν: • Κωδικό • x,y εικονοσυντεταγμένες. • Γεωδαιτικές συντεταγμένες X,Y,Z (όλα τα Ζ=0 Επίπεδο Πεδίο Ελέγχου) (Υπάρχει γνώση της διάστασης του κανάβου)

16. Διαδικασία Συνόρθωσης

17. Διαδικασία Συνόρθωσης • Η συνθήκη της συγγραμμικότητας δεν είναι γραμμική. • Γραμμικοποίηση κατά Taylor • Εύρεση αρχικών τιμών από συντελεστές αναγωγής • Εκτέλεση ΜΕΤ • Διόρθωση τιμών αγνώστων Επανάληψη

18. Εικόνες Βαθμονόμησης

19. Αποτέλεσμα Βαθμονόμησης Interior c1 = 949.67 pix ± 0.17 xo1 = -18.59 pix ± 0.39 yo1 = 0.19 pix ± 0.26 k11 = -1.396e-007 ± 2.078e-009 k21 = -6.805e-014 ± 1.414e-014 c2 = 962.13 pix ± 0.15 xo2 = -66.98 pix ± 0.40 yo2 = -6.12 pix ± 0.27 k12 = -1.620e-007 ± 1.506e-009 k22 = -4.440e-014 ± 7.713e-015 Relative Bx = 37.311 cm ± 0.007 By = 1.960 cm ± 0.004 Bz = -13.860 cm ± 0.018 om12 = -7.411 deg ± 0.027 ph12 = 41.290 deg ± 0.032 ka12 = 13.982 deg ± 0.014 18 image pairs 6335 points σo= 0.19pix maxresidual = 0.64pix

20. Περιγραφή επιπολικών εικόνων • Γνωστός εσωτερικός και σχετικός προσανατολισμός -> δημιουργία εικόνων παράλληλες στην βάση. • Οι εικόνες αυτές έχουν τις ιδιότητες της «κανονικής περίπτωσης στερεοζεύγους». • Προκύπτουν εφαρμόζοντας μία στροφή των εικόνων γύρω από το προβολικό κέντρο (επανασύσταση) • Η αναζήτηση ομολόγων σημείων γίνεται σε μία διάσταση αντί για δύο.

21. Ο b εP p Ο p eP εP επιπολικές ευθείες επιπολικό επίπεδο P Επιπολική Γεωμετρία

22. Δημιουργία επιπολικών εικόνων • Οι εικόνεςέχουν πίνακες στροφής R1(= Ι), R2(= R12) και μετασχηματίζονται σε αυστηρά κατακόρυφες μέσω των στροφών: R1T(= Ι) και R2T • Οι αυστηρά κατακόρυφες εικόνες ανάγονται κατόπιν στο σύστημα της βάσης μέσω της κοινής στροφής: RB = RΩRΦRΚ όπου: • Έχουμε M1 = RB R1T και M2 = RB R2T • Οπότε για κάθε κανονικοποιημένη εικόνα:

23. Δημιουργία επιπολικών εικόνων • Γίνεται να δημιουργηθούν εικόνες απαλλαγμένες από σφάλματα ακτινικής διαστροφής. • Με τις αντίστροφες σχέσεις υπολογίζονται σημεία της αρχικής εικόνας στα οποία θα γίνει παρεμβολή τόνου. Διαδικασία επανασύστασης.

24. Επιπολική γεωμετρία ->Εύρεση ομολογιών σκακιέρας • Η ανάποδη διαδικασία μπορεί να βοηθήσει στην εύρεση των ομολογιών μεταξύ των κόμβων της σκακιέρας • Αφαίρεση του κόκκινου τετραγώνου

25. Αυτόματη εύρεση ομολογιών Καλησπεράκης, 2010

26. Αυτόματη εύρεση ομολογιών Καλησπεράκης, 2010

27. Αυτόματη εύρεση ομολογιών Καλησπεράκης, 2010

28. Επιπολική γεωμετρία ->Εύρεση ομολογιών σκακιέρας

29. Επιπολική γεωμετρία ->Εύρεση ομολογιών σκακιέρας

30. Εύρεση ομολογιών σκακιέρας

31. Χρήση βαθμονομημένης στερεοκάμερας για δημιουργία 3D σαρωτή • Το αντικείμενο προς σάρωση παρατηρείται από δύο ακίνητες βαθμονομημένες κάμερες με σταθερή σχετική θέση. • Με ένα χειροκίνητο επίπεδο laser σαρώνεται το αντικείμενο. • Συγχρονισμένα καρέ από τις 2 κάμερες. Έχουνένα προφίλ που δημιουργείται από την τομή του αντικειμένου με το laser.

32. Περιγραφή συστήματος Εκτροπή γραμμής λόγω αναγλύφου. Παραμένει απροσδιοριστία σημείων.

33. Περιγραφή συστήματος • Αφαίρεση υποβάθρου • Επιπολικη γεωμετρία • 2D ► 1D ανίχνευση ομολόγων • Εμπροσθοτομία

34. Περιγραφή συστήματος • Μέσωεπιπολικής γεωμετρίας η ανίχνευση ομολόγου γίνεται μόνο πάνω στην επιπολική ευθεία. • Αναζήτηση της θέσης στην επιπολική ευθεία που παρουσιάζεται η μέγιστη τιμή φωτεινότητας.

35. Εμπροσθοτομία • Στα ομόλογα σημεία που έχουν τελικά εξαχθεί, γίνεται φωτογραμμετρική εμπροσθοτομία. • Προκύπτει ένα νέφος σημείων ΧΥΖ.

36. Λήψη εικόνων • Πρέπει να απομονωθεί το υπόβαθρο • Έχοντας μερικά καρέ (frames) είναι δυνατόν να εξαχθεί το υπόβαθρο από τον χρονική μεσαία τιμή των εικόνων αυτών. • Από κάθε εικόνα αφαιρείται το υπόβαθρο και μένει μόνο η γραμμή του laser. • Γίνεται επιπολική επανασύσταση.

37. Μείωση θορύβου • Φίλτρα Εξομάλυνσης • Μεσαίας τιμής 3x3 παράθυρο • Gaussian 3x3 παράθυρο σ = 1

38. 2D ► 1D Ανίχνευση • Έχοντας επιπολικές εικόνες απαλλαγμένες από θόρυβο, η ανίχνευση ομολόγου γίνεται μόνο πάνω στην επιπολική ευθεία. • Γίνεται αναζήτηση της θέσης στην επιπολική ευθεία που παρουσιάζεται η μέγιστη τιμή.

39. b a c d Εύρεση Μεγίστου • Ένας τρόπος είναι να προσαρμοστεί μία συνάρτηση στις τιμές του γκρίζου και να βρεθεί το μέγιστο της καμπύλης αυτής. • Κάθε σημείο έχει y συντεταγμένη από την επιπολική γραμμή (κοινό για τις δύο εικόνες) και χ συντεταγμένη από την θέση της μέγιστης τιμής της καμπύλης (ένα χ για κάθε εικόνα). • Επαναλαμβάνεται η διαδικασία για την επόμενη επιπολική γραμμή και για τις δύο εικόνες.

40. Εμπροσθοτομία • Για κάθε ζευγάρι στις επιπολικές εικόνες προκύπτουν: x1 x2 y • Οι εικόνες αυτές έχουν εξωτερικούς: Χο Υο Ζο ω φ κ Αριστερή 0 00000 Δεξιά Β 0 0000 (Β2 = Βχ2 + By2 + Bz2) • Οπότε οι συντεταγμένες των σημείων προκύπτουν από τους τύπους:

41. Νέφος σημείων • Για κάθε επιπολική γραμμή προκύπτει ένα σημείο του 3D χώρου. • Η διαδικασία επαναλαμβάνεται για κάθε επιπολική γραμμή του ζεύγους. • Έπειτα γίνεται ξανά η ίδια σειρά υπολογισμών για το επόμενο ζεύγος, κ.ο.κ. • Τελικά προκύπτει ένα νέφος σημείων με 3D συντεταγμένες.

42. Συνένωση Νεφών Ποια R T ώστε Σmin(d)

43. Συνένωση Νεφών

45. Παρεμβολή τόνου • Γνωρίζοντας τις ΧΥΖ συντεταγμένες των σημείων γίνεται να βρεθεί και τόνος RGB από την εικόνα του υποβάθρου εφαρμόζοντας την συνθήκη της συγγραμμικότητας.