Ανίχνευση αμπελώνων με φίλτρα Gabor, από εναέριες φωτογραφίες, χρησιμοποιώντας μια μη-επιβλεπόμενη προσέγγιση.

Από RemoteSensing Wiki

(Διαφορές μεταξύ αναθεωρήσεων)
Μετάβαση σε: πλοήγηση, αναζήτηση
 
(8 ενδιάμεσες αναθεωρήσεις δεν εμφανίζονται.)
Γραμμή 1: Γραμμή 1:
-
[[category:Νέα δοκιμαστική υποκατηγορία]]
+
[[category:Ψηφιακή Τηλεπισκόπηση]]
-
<b> Ανίχνευση αμπελώνων με φίλτρα Gabor, από εναέριες φωτογραφίες, χρησιμοποιώντας μια μη-επιβλεπόμενη προσέγγιση.</b><br>
+
<b> Ανίχνευση αμπελώνων με φίλτρα Gabor [http://en.wikipedia.org/wiki/Gabor_filter], από εναέριες φωτογραφίες, χρησιμοποιώντας μια μη-επιβλεπόμενη προσέγγιση.</b><br>
Πρωτότυπος τίτλος :  A non-supervised approach using Gabor filters for vine-plot detection in aerial images<br>
Πρωτότυπος τίτλος :  A non-supervised approach using Gabor filters for vine-plot detection in aerial images<br>
-
Πηγή : Gilles Rabatel, Carole Delenne, Michel Deshayes, Computers and Electronics in Agriculture 62 (2008) pp 159–168 [http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6T5M-4RWB0PX-1-N&_cdi=5006&_user=83473&_pii=S0168169907002487&_origin=search&_zone=rslt_list_item&_coverDate=07%2F31%2F2008&_sk=999379997&wchp=dGLzVlz-zSkWA&md5=779b60008c2984f2efd5a060d62e0c11&ie=/sdarticle.pdf]
+
Πηγή : Gilles Rabatel, Carole Delenne, Michel Deshayes, Computers and Electronics in Agriculture, volume 62, issue 2, July,2008, pp 159–168 [http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6T5M-4RWB0PX-1-N&_cdi=5006&_user=83473&_pii=S0168169907002487&_origin=search&_zone=rslt_list_item&_coverDate=07%2F31%2F2008&_sk=999379997&wchp=dGLzVlz-zSkWA&md5=779b60008c2984f2efd5a060d62e0c11&ie=/sdarticle.pdf]
<br><br>
<br><br>
<b>Περίληψη</b><br><br>
<b>Περίληψη</b><br><br>
-
Η χαρτογράφηση χρήσεων γης και ιδιαίτερα ο εντοπισμός και η αναγνώριση των αμπελώνων αποτελεί χρήσιμο εργαλείο ιδιαίτερα σε περιοχές που κυριαρχούν τέτοιου είδους καλλιέργειες όπως η Γαλλία. Σ αυτό το πλαίσιο, η διαθεσιμότητα ενός αυτόματου εργαλείου για την ανίχνευση των αμπελώνων θα ήταν πολύ χρήσιμο. Λόγω των περιοδικών σχεδίων αυτών των καλλιεργιών, η ανάλυση συχνότητας φαίνεται να είναι ένα πολύτιμο εργαλείο για την ανίχνευσή τους στις εναέριες εικόνες. Αναπτύχθηκε μία επαναλαμβανόμενη διαδικασία που χρησιμοποιεί τον αλγόριθμο μετασχηματισμού κατά Fourier ώστε να ικανοποιήσει αυτή την ανάγκη. Αυτό οδηγεί στον χωρισμό των αμπελώνων με τα όρια της πολυγωνικής μορφής και τον χαρακτηρισμό της ακριβής εκτίμησης τους πλάτους interrow και τον προσανατολισμό των σειρών. Τελικά πραγματοποιείται η ανίχνευση περισσότερων από 84% των επιφανειών αυτών από δεδομένα πολύ υψηλών χωρικών συχνοτήτων.<br><br>
+
Η χαρτογράφηση χρήσεων γης και ιδιαίτερα ο εντοπισμός και η αναγνώριση των αμπελώνων αποτελεί ιδιαίτερα χρήσιμο εργαλείο σε περιοχές που κυριαρχούν τέτοιου είδους καλλιέργειες, όπως για παράδειγμα στη Γαλλία. Σε αυτό το πλαίσιο, η διαθεσιμότητα ενός αυτόματου εργαλείου για την ανίχνευση των αμπελώνων θα ήταν πολύ χρήσιμο. Λόγω των επαναλαμβανόμενων (περιοδικών) σχημάτων αυτών των καλλιεργιών, η ανάλυση συχνότητας φαίνεται να είναι ένα πολύτιμο εργαλείο για την ανίχνευσή τους στις εναέριες εικόνες. Αναπτύχθηκε μία επαναλαμβανόμενη διαδικασία που χρησιμοποιεί τον αλγόριθμο μετασχηματισμού κατά Fourier [http://en.wikipedia.org/wiki/Fourier_transform] ώστε να ικανοποιήσει αυτή την ανάγκη. Αυτό οδηγεί στον ευκρινή χωρισμό των αμπελώνων κάνοντας διακριτά τα όριά τους (της πολυγωνικής μορφής τους) και έχοντας ακριβή εκτίμηση του πλάτους και του προσανατολισμού των σειρών των αμπελιών. Τελικά πραγματοποιείται η ανίχνευση, περισσότερων από 84% των επιφανειών αυτών, από δεδομένα πολύ υψηλών χωρικών συχνοτήτων.<br><br>
<b>Η περιοχή μελέτης</b><br><br>
<b>Η περιοχή μελέτης</b><br><br>
-
Τα πειράματα πραγματοποιήθηκαν σε μία έκταση 200 εκταρίων που αποτελεί υποσύνολο του υδροκρίτη La Peyne και βρίσκεται στην περιοχή Languedoc–Roussillon της Γαλλίας[http://fr.wikipedia.org/wiki/Languedoc-Roussillon]. Αυτή η ζώνη αποτελεί αντιπροσωπευτική ζώνη της γαλλικής μεσογειακής παράκτιας πεδιάδας όσον αφορά τις γεωργικές πρακτικές και τη διαχείριση των αμπελώνων. Παρά τη γενική μείωση στις μέρες μας, η καλλιέργεια των αμπελώνων είναι ακόμη κυρίαρχη στην περιοχή μελέτης με ποσοστό 70% έναντι άλλου είδους καλλιεργειών. Τα δεδομένα συλλέχθηκαν κατά τη διάρκεια της πρώτης εβδομάδας του Ιουλίου το 2005 όπου η άνθιση τους ήταν ορατή από τις αεροφωτπγραφίες καθώς και η περιοδικότητα των σχεδίων τους. Χρησιμοποιήθηκε μία ψηφιακή κάμερα SONY dsc-P1509) και αποκτήθηκαν εικόνες στα φυσικά χρώματα (κόκκινο, πράσινο, μπλέ). Οι εικόνες αυτές διορθώθηκαν γεωμετρικά από το Γεωγραφικό Γαλλικό Ινστιτούτο ώστε να προκύψουν ορθοφωτογραφίες  χωρικής ανάλυσης 50cm.  
+
Τα πειράματα πραγματοποιήθηκαν σε μία έκταση 200 εκταρίων (2000 στρεμμάτων) η οποία αποτελεί υποσύνολο του υδροκρίτη La Peyne και βρίσκεται στην περιοχή Languedoc – Roussillon της Γαλλίας [http://fr.wikipedia.org/ wiki/Languedoc - Roussillon]. Αυτή η ζώνη αποτελεί αντιπροσωπευτική ζώνη της γαλλικής μεσογειακής παράκτιας πεδιάδας όσον αφορά τις γεωργικές πρακτικές και τη διαχείριση των αμπελώνων. Παρά τη γενική μείωση στις μέρες μας, η καλλιέργεια των αμπελώνων είναι ακόμη κυρίαρχη στην περιοχή μελέτης με ποσοστό 70% έναντι άλλου είδους καλλιεργειών. Τα δεδομένα συλλέχθηκαν κατά τη διάρκεια της πρώτης εβδομάδας του Ιουλίου το 2005 όπου η άνθιση τους ήταν ορατή από τις αεροφωτογραφίες καθώς και η περιοδικότητα του προτύπου. Χρησιμοποιήθηκε μία ψηφιακή κάμερα (SONY dsc-P1509) και αποκτήθηκαν εικόνες στα φυσικά χρώματα (κόκκινο, πράσινο, μπλέ -RGB). Οι εικόνες αυτές διορθώθηκαν γεωμετρικά και λόγω αναγλύφου, από το Γεωγραφικό Γαλλικό Ινστιτούτο [http://en.wikipedia.org/wiki/Institut_G%C3%A9ographique_National], ώστε να προκύψουν ορθοφωτογραφίες  χωρικής ανάλυσης 50cm.
[[Εικόνα:garg_gabor_1_02.jpg|right|thumb||Εικόνα 1 : La Peine, 11240 Escueillens-et-Saint-Just-de-Bélengard, Γαλλία, πηγή : Google Earth, La Peine, 11240 Escueillens-et-Saint-Just-de-Bélengard, Γαλλία,  43° 5'49.10"Β,    2° 1'54.30"Α]]
[[Εικόνα:garg_gabor_1_02.jpg|right|thumb||Εικόνα 1 : La Peine, 11240 Escueillens-et-Saint-Just-de-Bélengard, Γαλλία, πηγή : Google Earth, La Peine, 11240 Escueillens-et-Saint-Just-de-Bélengard, Γαλλία,  43° 5'49.10"Β,    2° 1'54.30"Α]]
<br><br>
<br><br>
<b>Μεθοδολογία</b><br><br>
<b>Μεθοδολογία</b><br><br>
-
Η παρούσα μελέτη είναι βασισμένη στις τεχνικές επεξεργασίας εικόνας που εξετάζουν αποκλειστικά τη χωρική δομή των αμπελώνων στις εναέριες εικόνες και όχι τις ραδιομετρικές ιδιότητές τους. Οι τεχνικές αυτές χρησιμοποιούν εικόνες γκρίζων τόνων με την προϋπόθεση ικανοποιητικής αντίθεσης μεταξύ των σειρών αμπελώνων και των interrows for the discontinue structure to be visible. Επιπλέον παρουσιάζεται ένα βήμα κανονικοποίησης  ώστε τα αποτελέσματα να καταστούν ανεξάρτητα από το αρχικό επίπεδο αντίθεσης. Επομένως η προτεινόμενη μέθοδος μπορεί να εφαρμοστεί με διάφορους τύπους πηγών εικόνας. Ανάλογα με τις φασματικές διαθέσιμες ζώνες σε ένα δεδομένο πλαίσιο εφαρμογής, οι τλονοι του γκρι μπορούν να είναι για παράδειγμα ένας υπολογισμός NDVI [http://en.wikipedia.org/wiki/Normalized_Difference_Vegetation_Index]( ομαλοποιημένος δείκτης διαφοράς βλάστησης, συνδυάζοντας τις κόκκινες και near-infraded ζώνες), η φωτεινότητα υπολογίζεται από μια RGB εικόνα, η μια ιδιαίτερη ζώνη μιας εικόνας.  
+
Η παρούσα μελέτη είναι βασισμένη στις τεχνικές επεξεργασίας εικόνας που εξετάζουν αποκλειστικά τη χωρική δομή των αμπελώνων στις εναέριες εικόνες και όχι τις ραδιομετρικές ιδιότητές τους. Οι τεχνικές αυτές χρησιμοποιούν εικόνες τόνων του γκρι, με την απαίτηση (προϋπόθεση) ικανοποιητικής χρωματικής αντίθεσης μεταξύ των σειρών των αμπελώνων. Επιπλέον παρουσιάζεται ένα βήμα τυποποίησης της διαδικασίας, ώστε τα αποτελέσματα να καταστούν ανεξάρτητα (όσο αυτό είναι δυνατόν) από το αρχικό επίπεδο χρωματικής αντίθεσης. Ως αποτέλεσμα, η προτεινόμενη μέθοδος μπορεί να εφαρμοστεί σε διαφορετικούς τύπους πηγών εικόνας. Ανάλογα με τις φασματικές διαθέσιμες ζώνες σε ένα δεδομένο πλαίσιο εφαρμογής, οι τόνοι του γκρι, μπορούν να χρησιμοποιηθούν για παράδειγμα σε έναν υπολογισμό NDVI [http://en.wikipedia.org/wiki/ Normalized_Difference_Vegetation_Index] (ομαλοποιημένος δείκτης διαφοράς βλάστησης, συνδυάζοντας τις κόκκινες και εγγύς υπέρυθρες περιοχές του ηλεκτρομαγνητικού φάσματος) ή στον υπολογισμό της φωτεινότητας των pixels για μια ιδιαίτερη υποπεριοχή μιας εικόνας.  
-
Η βασική ιδέα είναι να απομονωθεί κάθε μεμονωμένη πλοκή με την επιλογή των αντίστοιχων συχνοτήτων στο φάσμα Fourier [http://en.wikipedia.org/wiki/Fourier_analysis], χρησιμοποιώντας ένα συγκεκριμένο φίλτρο Gabor [http://en.wikipedia.org/wiki/Gabor_filter]. Τα φίλτρα Gabor έχουν χρησιμοποιηθεί ευρέως και ως θεμελιώδης λειτουργία αποσύνθεσης κυματιδίων και για την κατάτμηση υφής. Η διαδικασία φιλτραρίσματος του Gabor, πρέπει να εφαρμοστεί σε περιορισμένο μέγεθος εικόνας (500x500 pixels) και για υπολογιστικούς λόγους αλλά και για να παρθούν τα εκμεταλλεύσιμα φάσματα Fourier. Επειδή η περιοχή μελέτης είναι πολύ μεγάλη, αρχικά γίνεται ένας χωρισμός της εικόνας. Για να ανιχνευθεί η κάθε πλοκή αμπελώνων το φίλτρο πρέπει να εφαρμοστεί διαδοχικά στις κύριες αιχμές εύρους στο φάσμα Fourier με τη ρύθμιση της κεντρικής συχνοτητάς του. Μετά την εφαρμογή του φίλτρου Gabor, ένα threshold εφαρμόζεται στην παραγόμενη εικόνα. Kατά συνέπεια λαμβάνουμε μια δυαδική εικόνα, στην οποία κάθε αντικείμενο είναι πιθανό να είναι αμπελώνας με τα χαρακτηριστικά του προσανατολισμού και του πλάτους interrow του τρέχοντος φίλτρου. Εντούτοις, η υποψήφια πλοκή μπορεί να μην περιλαμβάνει πλήρως τα ακριβή χαρακτηριστικά αλλά παραπλήσια. Βασικός παράγοντας τέτοιας αβεβαιότητας είναι πιθανό να είναι η αρχική διαδικασία χωρισμού των εικόνων. Σ αυτή τη περίπτωση, κάθε υποψήφια πλοκή υποβάλλεται σε επαναληπτική διαδικασία φιλτραρίσματος μέσω ενός νέου φίλτρου Gabor.
+
Η βασική ιδέα είναι να απομονωθεί κάθε μεμονωμένη εμφάνιση μιας υποπεριοχής της εικόνας (θα την ονομάζω αγροτεμάχιο από δω και πέρα), με την επιλογή των αντίστοιχων συχνοτήτων στο φάσμα Fourier [http://en.wikipedia.org/wiki/Fourier_analysis], χρησιμοποιώντας ένα συγκεκριμένο φίλτρο Gabor [http://en.wikipedia.org/wiki/ Gabor_filter]. Τα φίλτρα Gabor έχουν χρησιμοποιηθεί ευρέως και ως θεμελιώδης λειτουργία αποσύνθεσης κυματιδίων και για την κατάτμηση υφής. Έχει βρεθεί ότι η διαδικασία φιλτραρίσματος μέσω του φίλτρου Gabor, πρέπει να εφαρμοστεί σε περιορισμένο μέγεθος εικόνας (έως 500x500 pixels) και για υπολογιστικούς λόγους αλλά και για να παρθούν τα εκμεταλλεύσιμα φάσματα Fourier. Επειδή η περιοχή μελέτης είναι πολύ μεγάλη, αρχικά γίνεται ένας χωρισμός της εικόνας σε επιμέρους περιοχές. Για να ανιχνευθεί το κάθε αγροτεμάχιο αμπελώνων, το φίλτρο πρέπει να εφαρμοστεί διαδοχικά στις κύριες συνιστώσες διεύθυνσης της χρωματικής αντίθεσης του αμπελώνα, με τη ρύθμιση της κεντρικής συχνοτητάς του φίλτρου Fourier. Μετά την εφαρμογή του φίλτρου Gabor, ένα φασματικό κατώφλι περιορισμού συχνοτήτων, εφαρμόζεται στην παραγόμενη εικόνα. Ως αποτέλεσμα, λαμβάνουμε μια δυαδική εικόνα [http://wapedia.mobi/en/Color_depth], στην οποία κάθε αντικείμενο το οποίο εμφανίζεται, είναι πιθανό να είναι αμπελώνας με τα χαρακτηριστικά του προσανατολισμού και του πλάτους που του προσέδωσε το φίλτρο. Βασικός παράγοντας της πιθανότητας μιας τέτοιας αβεβαιότητας, αποτελεί η αρχική διαδικασία χωρισμού της εικόνας σε περιοχές. Σε αυτή τη περίπτωση, κάθε υποψήφιο αγροτεμάχιο (όπως αναφέρεται παραπάνω) υποβάλλεται σε επαναληπτική διαδικασία φιλτραρίσματος μέσω ενός νέου φίλτρου Gabor.
 +
 
[[Εικόνα:garg_gabor_2_02.jpg|right|thumb||Εικόνα 2 :Φάσμα Fourier εφαρμοσμένο στις πλοκές αμπελώνων, [http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6T5M-4RWB0PX-1-N&_cdi=5006&_user=83473&_pii=S0168169907002487&_origin=search&_zone=rslt_list_item&_coverDate=07%2F31%2F2008&_sk=999379997&wchp=dGLzVlz-zSkWA&md5=779b60008c2984f2efd5a060d62e0c11&ie=/sdarticle.pdf]πηγή]]
[[Εικόνα:garg_gabor_2_02.jpg|right|thumb||Εικόνα 2 :Φάσμα Fourier εφαρμοσμένο στις πλοκές αμπελώνων, [http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6T5M-4RWB0PX-1-N&_cdi=5006&_user=83473&_pii=S0168169907002487&_origin=search&_zone=rslt_list_item&_coverDate=07%2F31%2F2008&_sk=999379997&wchp=dGLzVlz-zSkWA&md5=779b60008c2984f2efd5a060d62e0c11&ie=/sdarticle.pdf]πηγή]]
[[Εικόνα:garg_gabor_3_02.jpg|right|thumb||Εικόνα 3 :Εφαρμογή φίλτρου Gabor στις αντίστοιχες πλοκές αμπελώνων [http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6T5M-4RWB0PX-1-N&_cdi=5006&_user=83473&_pii=S0168169907002487&_origin=search&_zone=rslt_list_item&_coverDate=07%2F31%2F2008&_sk=999379997&wchp=dGLzVlz-zSkWA&md5=779b60008c2984f2efd5a060d62e0c11&ie=/sdarticle.pdf]πηγή]]
[[Εικόνα:garg_gabor_3_02.jpg|right|thumb||Εικόνα 3 :Εφαρμογή φίλτρου Gabor στις αντίστοιχες πλοκές αμπελώνων [http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6T5M-4RWB0PX-1-N&_cdi=5006&_user=83473&_pii=S0168169907002487&_origin=search&_zone=rslt_list_item&_coverDate=07%2F31%2F2008&_sk=999379997&wchp=dGLzVlz-zSkWA&md5=779b60008c2984f2efd5a060d62e0c11&ie=/sdarticle.pdf]πηγή]]
-
Για την επικύρωση των αποτελεσμάτων, συλλέχθηκαν πληροφορίες επίγειας έρευνας ταυτόχρονα με την απόκτηση της εικόνας, όπως η χρήση εδάφους, σύσταση χώματος επιφάνειας, πλάτος σειρών και interrows κ.α. Οι 160 πλοκές αμπελώνων έχουν μετατραπεί πλέον σε μια βάση δεδομένων GIS. Οι επικυρωμένες αυτές πλοκές θα θεωρηθούν ως πλοκές στη διαδικασία επικύρωσης.Ο κύριος σκοπός της αυτόματης ανίχνευσης πλοκών είναι ο προσδιορισμός των ορίων τους. επομένως έχει εκτελεστεί μια αυτοματοποιημένη επικύρωση βάσης πλοκών, συγκρίνοντας τα αυτόματα με το χέρι, σύμφωνα με το επικαλυπτόμενο ποσοστό τους. διαφορετικές περιπτώσεις έχουν καθοριστεί όπως καλή κατάτμηση, κατάτμηση, κατάτμηση κατά της οποία μία πλοκή εμφανίζεται σαν αρκετές ανιχνευμένες πλοκές, μόνο ένα μέρος της πλοκής έχει ανιχνευθεί, πλοκές δεν έχουν καθόλου ανιχνευθεί κ.α.<br><br>
+
Για την επικύρωση των αποτελεσμάτων, έγιναν επίγειοι έλεγχοι και συλλέχθηκαν πληροφορίες όπως η κάλυψη του εδάφους, η σύσταση του χώματος, οι διαστάσεις των αμπελώνων και του τρόπου χωροθέτησής τους μέσα στο κάθε αμπέλι κ.α. Οι 160 συστάδες αμπελώνων έχουν εισαχθεί σε μια βάση δεδομένων GIS. Οι επικυρωμένες αυτές συστάδες θα θεωρηθούν ως αγροτεμάχια στη διαδικασία επικύρωσης του αμπελιού (με τον τρόπο που αναφέρθηκε παραπάνω). Ο κύριος στόχος της αυτόματης ανίχνευσης αμπελιών, είναι ο προσδιορισμός των ορίων τους. Με τη διαδικασία που αναφέρθηκε, έχει εκτελεστεί μια αυτοματοποιημένη επικύρωση συστάδων οι οποίες έχουν τα χαρακτηριστικά των αμπελώνων, συγκρίνοντας τα αποτελέσματα της αυτοματοποιημένης διαδικασίας και με φωτοερμηνεία, όπου ανάλογα γίνεται η επιβεβαίωσή τους (σε αμπελώνες ή όχι). Έχουν καθοριστεί κάποιες κατηγορίες στο τελικό αποτέλεσμα, όπως καλή ή μέτρια κατάτμηση, ανάλογα το ποσοστό εμφάνισης (ανίχνευσης) του αγροτεμαχίου (με τον τρόπο που αναφέρθηκε παραπάνω) στην τελική εικόνα, μετά την εφαρμογή του φίλτρου.<br><br>
<b>Αποτελέσματα</b><br><br>
<b>Αποτελέσματα</b><br><br>
-
Οι τρεις πιο αντιπροσωπευτικοί τύποι αποτελεσμάτων είναι η σωστά ανιχνευμένες πλοκές με ποσοστό 64%,  οι under-segmented πλοκές με ποσοστό 15,8% και οι μη εντοπισμένες πλοκές  με ποσοστό 11,4%.  Οι under-segmented πλοκές αντιστοιχούν στην ομαδοποίηση από τις γειτονικές πλοκές που έχουν τον ίδιο προσανατολισμό σειρών and πλάτος interrows και χωρίζονται μόνο από ένα στενό δρόμο ή τάφρο. Μερικοί από αυτούς δεν είναι χωρικά χωρισμένοι και διαφέρουν μόνο στη σύσταση χώματος της επιφάνειας μεταξύ των σειρών η από μερικά χαρακτηριστικά δυσδιάκριτα στις εναέριες εικόνες όπως η ηλικία ή το ύψος. Η αναλογία έλλειψης ανίχνευσης είανι σχετικά μικρή 11,4% και αφορά κυρίως τις μικρές πλοκές. Πράγματι 9 στις 13 μη ανιχνευμένες πλοκές είναι μικρότερες από 0,2 εκτάρια. Συνεπώς, λιγότερο από 5% της συνολικής περιοχής αμπελώνων δεν ανιχνεύεται. Μέσα στις τέσσερις άλλες μη-ανιχνευμένες πλοκές, one has most of its vinetrees
+
Οι τρεις πιο αντιπροσωπευτικοί τύποι αποτελεσμάτων είναι οι ορθά ανιχνευμένες συστάδες (αγροτεμάχια) με ποσοστό 64%,  οι μερικώς εντοπισμένες με ποσοστό 15,8% και οι μη εντοπισμένες συστάδες, με ποσοστό 11,4%.  Οι μερικώς εντοπισμένες, αντιστοιχούν στην ομαδοποίηση από τις γειτονικές συστάδες που έχουν τον ίδιο προσανατολισμό σειρών και πλάτους και διαχωρίζονται μόνο από ένα στενό δρόμο ή τάφρο. Μερικοί από αυτούς δεν είναι χωρικά διαχωρισμένοι και διαφέρουν μόνο στη σύσταση χώματος της επιφάνειας μεταξύ των σειρών η από μερικά χαρακτηριστικά δυσδιάκριτα στις εναέριες εικόνες, όπως η ηλικία ή το ύψος των αμπελώνων. Η αναλογία έλλειψης ανίχνευσης είναι σχετικά μικρή (περίπου 11,4%) και αφορά κυρίως τις μικρές συστάδες. Πράγματι 9 στις 13 μη ανιχνεύσιμες συστάδες είναι μικρότερες από 2 στρέμματα. Συνεπώς, λιγότερο από 5% της συνολικής επιφάνειας των αμπελώνων δεν ανιχνεύεται ορθά με το φίλτρο Gabor. Γενικά η χρωματική αντίθεση των εικόνων που χρησιμοποιούνται αλλά και η εσωτερική δομή των αμπελώνων αναλογικά με τις διαστάσεις των αμπελιών και του τρόπου που είναι χωροθετημένα, βελτιώνει ή όχι, τα αποτελέσματα αυτόματης ανίχνευσής τους με την παραπάνω μεθοδολογία, προσθέτοντας σφάλματα στη μέθοδο.
-
missing, και η άλλη έχει ένα μικρό πλάτος interrow (1.7 μ), κοντά στο όριο ανίχνευσης σύμφωνα με την ανάλυση εικόνας, ενώ τα δύο άλλα έχουν πολύ κακή αντίθεση στην εικόνα.
+
[[Εικόνα:garg_gabor_4_02.jpg|right|thumb||Εικόνα 4 :Φωτοερμηνευτική και αυτόματη κατάτμηση εικόνας [http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6T5M-4RWB0PX-1-N&_cdi=5006&_user=83473&_pii=S0168169907002487&_origin=search&_zone=rslt_list_item&_coverDate=07%2F31%2F2008&_sk=999379997&wchp=dGLzVlz-zSkWA&md5=779b60008c2984f2efd5a060d62e0c11&ie=/sdarticle.pdf]πηγή]]
-
τέσσερεις πλοκές έχουν λανθασμένα segmented. Μία από αυτές είναι λόγω υψηλής διαμήκους μετάβασης επιπέδου γκρι ενός δρόμου, ο οποίος παράγει τις αιχμές εύρους στη σειρά αναζήτησης. Οι άλλες κατηγοριοποιούνται σε μη-καλλιεργημένες πλοκές, αλλά έχουν οργωθεί  πρόσφατα και επομένως παρουσιάζουν μια παράλληλη δομή.
+
-
[[Εικόνα:garg_gabor_4_02.jpg|right|thumb||Εικόνα 4 :Manual and automatic segmentation [http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6T5M-4RWB0PX-1-N&_cdi=5006&_user=83473&_pii=S0168169907002487&_origin=search&_zone=rslt_list_item&_coverDate=07%2F31%2F2008&_sk=999379997&wchp=dGLzVlz-zSkWA&md5=779b60008c2984f2efd5a060d62e0c11&ie=/sdarticle.pdf]πηγή]]
+
<br><br>
<br><br>
<b>Συμπεράσματα </b><br><br>
<b>Συμπεράσματα </b><br><br>
-
Η προτεινόμενη επαναλαμβανόμενη προσέγγιση έχει αποδείξει την αποδοτικότητά της για την ανίχνευση και το χαρακτηρισμό των πλοκών των αμπελώνων από πολλές απόψεις. Η κύρια πρωτοτυπία της μεθόδου αυτής είναι ότι τα αποτελέσματα είναι άμεσα διαθέσιμα σε μια μορφή πολυγωνική. Ένα άλλο σημαντικό πλεονέκτημα είναι ότι τα σωστά αποτελέσματα επιτυγχάνονται με το κόκκινο κανάλι στις διαθέσιμες εικόνες με τα φυσικά χρώματα. Επιπλέον, απ τη στιγμή που η κατάλληλη χωρική ανάλυση συνδέεται με την τοπική περιοδικότητα των σχεδίων, ο πιο χονδροειδής θα μπορούσε να χρησιμοποιηθεί  σε πολλές περιοχές αμπελοκαλλιέργειας, ιδιαίτερα στους ξηρούς (όπως στην Ισπανία), όπου τα interrows πλάτη είναι εώς 3m. Κατόπιν, οι δορυφορικές εικόνες όπως Ikonos & Quickbird θα μπορούσαν να χρησιμοποιηθούν. Το κύριο όριο αυτής της μεθόδου είναι ότι έχει αναπτυχθεί για να εφαρμοστεί στα γραμμικά σχέδια σειρών. Οι επακόλουθες τροποποιήσεις θα απαιτούνταν για να χρησιμοποιήσουν μια διαδικασία φιλτραρίσματος του Gabor, στις πλοκές αμπελώνων επίπεδων γραμμών όπως αντιμετωπίζεται σε ορισμένες περιοχές αμπελώνων με σημαντικές κλίσεις (Πορτογαλία, Ιταλία) [[Εικόνα:garg_gabor_5_02.jpg|right|thumb||Εικόνα 5 :Ποσοστό αποτελεσμάτων ανά κατηγορία [http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6T5M-4RWB0PX-1-N&_cdi=5006&_user=83473&_pii=S0168169907002487&_origin=search&_zone=rslt_list_item&_coverDate=07%2F31%2F2008&_sk=999379997&wchp=dGLzVlz-zSkWA&md5=779b60008c2984f2efd5a060d62e0c11&ie=/sdarticle.pdf]πηγή]]
+
Η προτεινόμενη αυτοματοποιημένη προσέγγιση, έχει αποδείξει την αποδοτικότητά της στην ανίχνευση και το χαρακτηρισμό των αμπελώνων. Βασικό της πλεονέκτημα, αποτελεί το γεγονός ότι τα αποτελέσματα είναι άμεσα διαθέσιμα σε μια μορφή πολυγωνική (vector αρχεία). Ένα άλλο πρακτικό πλεονέκτημα είναι ότι τα ορθά αποτελέσματα επιτυγχάνονται στην περιοχή του ηλεκτρομαγνητικού φάσματος, όπου βρίσκεται το κόκκινο κανάλι (περίπου στα 0.65μm) στις διαθέσιμες εικόνες με τα φυσικά χρώματα. Επιπλέον, απ τη στιγμή που η κατάλληλη χωρική ανάλυση συνδέεται με την περιοδικότητα των σχημάτων των αμπελώνων, το πιο συχνά επαναλαμβανόμενο πρότυπο που δημιουργείται, θα μπορούσε να χρησιμοποιηθεί  σε πολλές διαφορετικές περιοχές αμπελοκαλλιέργειας, οι οποίες παρουσιάζουν ανάλογες εμφανίσεις (όπως για παράδειγμα στην Ισπανία). Προφανώς αντί για αεροφωτογραφίες, θα μπορούσαν να χρησιμοποιηθούν οι δορυφορικές απεικονίσεις Ikonos [http://www.satimagingcorp.com/gallery-ikonos.html] και Quickbird [http://www.digitalglobe.com/index.php/85/QuickBird]. Ένα όριο αυτής της μεθόδου, όπως εφαρμόστηκε, είναι ότι έχει αναπτυχθεί για να εφαρμοστεί σε γραμμικά πρότυπα αμπελώνων. Οι επακόλουθες τροποποιήσεις θα ήταν απαραίτητες για να χρησιμοποιήσουν, σε διαδικασία φιλτραρίσματος μέσω του φίλτρου Gabor, στις συστάδες των αμπελώνων όπως εμφανίζονται σε ορισμένες περιοχές με σημαντικές κλίσεις στο έδαφος (Πορτογαλία, Ιταλία). [[Εικόνα:garg_gabor_5_02.jpg|right|thumb||Εικόνα 5 :Ποσοστό αποτελεσμάτων ανά κατηγορία [http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6T5M-4RWB0PX-1-N&_cdi=5006&_user=83473&_pii=S0168169907002487&_origin=search&_zone=rslt_list_item&_coverDate=07%2F31%2F2008&_sk=999379997&wchp=dGLzVlz-zSkWA&md5=779b60008c2984f2efd5a060d62e0c11&ie=/sdarticle.pdf]πηγή]]
-
 
+
-
 
+
-
<br><br><br>
+

Παρούσα αναθεώρηση της 09:06, 16 Μαρτίου 2011


Ανίχνευση αμπελώνων με φίλτρα Gabor [5], από εναέριες φωτογραφίες, χρησιμοποιώντας μια μη-επιβλεπόμενη προσέγγιση.
Πρωτότυπος τίτλος : A non-supervised approach using Gabor filters for vine-plot detection in aerial images
Πηγή : Gilles Rabatel, Carole Delenne, Michel Deshayes, Computers and Electronics in Agriculture, volume 62, issue 2, July,2008, pp 159–168 [6]

Περίληψη

Η χαρτογράφηση χρήσεων γης και ιδιαίτερα ο εντοπισμός και η αναγνώριση των αμπελώνων αποτελεί ιδιαίτερα χρήσιμο εργαλείο σε περιοχές που κυριαρχούν τέτοιου είδους καλλιέργειες, όπως για παράδειγμα στη Γαλλία. Σε αυτό το πλαίσιο, η διαθεσιμότητα ενός αυτόματου εργαλείου για την ανίχνευση των αμπελώνων θα ήταν πολύ χρήσιμο. Λόγω των επαναλαμβανόμενων (περιοδικών) σχημάτων αυτών των καλλιεργιών, η ανάλυση συχνότητας φαίνεται να είναι ένα πολύτιμο εργαλείο για την ανίχνευσή τους στις εναέριες εικόνες. Αναπτύχθηκε μία επαναλαμβανόμενη διαδικασία που χρησιμοποιεί τον αλγόριθμο μετασχηματισμού κατά Fourier [7] ώστε να ικανοποιήσει αυτή την ανάγκη. Αυτό οδηγεί στον ευκρινή χωρισμό των αμπελώνων κάνοντας διακριτά τα όριά τους (της πολυγωνικής μορφής τους) και έχοντας ακριβή εκτίμηση του πλάτους και του προσανατολισμού των σειρών των αμπελιών. Τελικά πραγματοποιείται η ανίχνευση, περισσότερων από 84% των επιφανειών αυτών, από δεδομένα πολύ υψηλών χωρικών συχνοτήτων.

Η περιοχή μελέτης

Τα πειράματα πραγματοποιήθηκαν σε μία έκταση 200 εκταρίων (2000 στρεμμάτων) η οποία αποτελεί υποσύνολο του υδροκρίτη La Peyne και βρίσκεται στην περιοχή Languedoc – Roussillon της Γαλλίας wiki/Languedoc - Roussillon. Αυτή η ζώνη αποτελεί αντιπροσωπευτική ζώνη της γαλλικής μεσογειακής παράκτιας πεδιάδας όσον αφορά τις γεωργικές πρακτικές και τη διαχείριση των αμπελώνων. Παρά τη γενική μείωση στις μέρες μας, η καλλιέργεια των αμπελώνων είναι ακόμη κυρίαρχη στην περιοχή μελέτης με ποσοστό 70% έναντι άλλου είδους καλλιεργειών. Τα δεδομένα συλλέχθηκαν κατά τη διάρκεια της πρώτης εβδομάδας του Ιουλίου το 2005 όπου η άνθιση τους ήταν ορατή από τις αεροφωτογραφίες καθώς και η περιοδικότητα του προτύπου. Χρησιμοποιήθηκε μία ψηφιακή κάμερα (SONY dsc-P1509) και αποκτήθηκαν εικόνες στα φυσικά χρώματα (κόκκινο, πράσινο, μπλέ -RGB). Οι εικόνες αυτές διορθώθηκαν γεωμετρικά και λόγω αναγλύφου, από το Γεωγραφικό Γαλλικό Ινστιτούτο [8], ώστε να προκύψουν ορθοφωτογραφίες χωρικής ανάλυσης 50cm.

Εικόνα 1 : La Peine, 11240 Escueillens-et-Saint-Just-de-Bélengard, Γαλλία, πηγή : Google Earth, La Peine, 11240 Escueillens-et-Saint-Just-de-Bélengard, Γαλλία, 43° 5'49.10"Β, 2° 1'54.30"Α



Μεθοδολογία

Η παρούσα μελέτη είναι βασισμένη στις τεχνικές επεξεργασίας εικόνας που εξετάζουν αποκλειστικά τη χωρική δομή των αμπελώνων στις εναέριες εικόνες και όχι τις ραδιομετρικές ιδιότητές τους. Οι τεχνικές αυτές χρησιμοποιούν εικόνες τόνων του γκρι, με την απαίτηση (προϋπόθεση) ικανοποιητικής χρωματικής αντίθεσης μεταξύ των σειρών των αμπελώνων. Επιπλέον παρουσιάζεται ένα βήμα τυποποίησης της διαδικασίας, ώστε τα αποτελέσματα να καταστούν ανεξάρτητα (όσο αυτό είναι δυνατόν) από το αρχικό επίπεδο χρωματικής αντίθεσης. Ως αποτέλεσμα, η προτεινόμενη μέθοδος μπορεί να εφαρμοστεί σε διαφορετικούς τύπους πηγών εικόνας. Ανάλογα με τις φασματικές διαθέσιμες ζώνες σε ένα δεδομένο πλαίσιο εφαρμογής, οι τόνοι του γκρι, μπορούν να χρησιμοποιηθούν για παράδειγμα σε έναν υπολογισμό NDVI Normalized_Difference_Vegetation_Index (ομαλοποιημένος δείκτης διαφοράς βλάστησης, συνδυάζοντας τις κόκκινες και εγγύς υπέρυθρες περιοχές του ηλεκτρομαγνητικού φάσματος) ή στον υπολογισμό της φωτεινότητας των pixels για μια ιδιαίτερη υποπεριοχή μιας εικόνας. Η βασική ιδέα είναι να απομονωθεί κάθε μεμονωμένη εμφάνιση μιας υποπεριοχής της εικόνας (θα την ονομάζω αγροτεμάχιο από δω και πέρα), με την επιλογή των αντίστοιχων συχνοτήτων στο φάσμα Fourier [9], χρησιμοποιώντας ένα συγκεκριμένο φίλτρο Gabor Gabor_filter. Τα φίλτρα Gabor έχουν χρησιμοποιηθεί ευρέως και ως θεμελιώδης λειτουργία αποσύνθεσης κυματιδίων και για την κατάτμηση υφής. Έχει βρεθεί ότι η διαδικασία φιλτραρίσματος μέσω του φίλτρου Gabor, πρέπει να εφαρμοστεί σε περιορισμένο μέγεθος εικόνας (έως 500x500 pixels) και για υπολογιστικούς λόγους αλλά και για να παρθούν τα εκμεταλλεύσιμα φάσματα Fourier. Επειδή η περιοχή μελέτης είναι πολύ μεγάλη, αρχικά γίνεται ένας χωρισμός της εικόνας σε επιμέρους περιοχές. Για να ανιχνευθεί το κάθε αγροτεμάχιο αμπελώνων, το φίλτρο πρέπει να εφαρμοστεί διαδοχικά στις κύριες συνιστώσες διεύθυνσης της χρωματικής αντίθεσης του αμπελώνα, με τη ρύθμιση της κεντρικής συχνοτητάς του φίλτρου Fourier. Μετά την εφαρμογή του φίλτρου Gabor, ένα φασματικό κατώφλι περιορισμού συχνοτήτων, εφαρμόζεται στην παραγόμενη εικόνα. Ως αποτέλεσμα, λαμβάνουμε μια δυαδική εικόνα [10], στην οποία κάθε αντικείμενο το οποίο εμφανίζεται, είναι πιθανό να είναι αμπελώνας με τα χαρακτηριστικά του προσανατολισμού και του πλάτους που του προσέδωσε το φίλτρο. Βασικός παράγοντας της πιθανότητας μιας τέτοιας αβεβαιότητας, αποτελεί η αρχική διαδικασία χωρισμού της εικόνας σε περιοχές. Σε αυτή τη περίπτωση, κάθε υποψήφιο αγροτεμάχιο (όπως αναφέρεται παραπάνω) υποβάλλεται σε επαναληπτική διαδικασία φιλτραρίσματος μέσω ενός νέου φίλτρου Gabor.

Εικόνα 2 :Φάσμα Fourier εφαρμοσμένο στις πλοκές αμπελώνων, [1]πηγή
Εικόνα 3 :Εφαρμογή φίλτρου Gabor στις αντίστοιχες πλοκές αμπελώνων [2]πηγή

Για την επικύρωση των αποτελεσμάτων, έγιναν επίγειοι έλεγχοι και συλλέχθηκαν πληροφορίες όπως η κάλυψη του εδάφους, η σύσταση του χώματος, οι διαστάσεις των αμπελώνων και του τρόπου χωροθέτησής τους μέσα στο κάθε αμπέλι κ.α. Οι 160 συστάδες αμπελώνων έχουν εισαχθεί σε μια βάση δεδομένων GIS. Οι επικυρωμένες αυτές συστάδες θα θεωρηθούν ως αγροτεμάχια στη διαδικασία επικύρωσης του αμπελιού (με τον τρόπο που αναφέρθηκε παραπάνω). Ο κύριος στόχος της αυτόματης ανίχνευσης αμπελιών, είναι ο προσδιορισμός των ορίων τους. Με τη διαδικασία που αναφέρθηκε, έχει εκτελεστεί μια αυτοματοποιημένη επικύρωση συστάδων οι οποίες έχουν τα χαρακτηριστικά των αμπελώνων, συγκρίνοντας τα αποτελέσματα της αυτοματοποιημένης διαδικασίας και με φωτοερμηνεία, όπου ανάλογα γίνεται η επιβεβαίωσή τους (σε αμπελώνες ή όχι). Έχουν καθοριστεί κάποιες κατηγορίες στο τελικό αποτέλεσμα, όπως καλή ή μέτρια κατάτμηση, ανάλογα το ποσοστό εμφάνισης (ανίχνευσης) του αγροτεμαχίου (με τον τρόπο που αναφέρθηκε παραπάνω) στην τελική εικόνα, μετά την εφαρμογή του φίλτρου.

Αποτελέσματα

Οι τρεις πιο αντιπροσωπευτικοί τύποι αποτελεσμάτων είναι οι ορθά ανιχνευμένες συστάδες (αγροτεμάχια) με ποσοστό 64%, οι μερικώς εντοπισμένες με ποσοστό 15,8% και οι μη εντοπισμένες συστάδες, με ποσοστό 11,4%. Οι μερικώς εντοπισμένες, αντιστοιχούν στην ομαδοποίηση από τις γειτονικές συστάδες που έχουν τον ίδιο προσανατολισμό σειρών και πλάτους και διαχωρίζονται μόνο από ένα στενό δρόμο ή τάφρο. Μερικοί από αυτούς δεν είναι χωρικά διαχωρισμένοι και διαφέρουν μόνο στη σύσταση χώματος της επιφάνειας μεταξύ των σειρών η από μερικά χαρακτηριστικά δυσδιάκριτα στις εναέριες εικόνες, όπως η ηλικία ή το ύψος των αμπελώνων. Η αναλογία έλλειψης ανίχνευσης είναι σχετικά μικρή (περίπου 11,4%) και αφορά κυρίως τις μικρές συστάδες. Πράγματι 9 στις 13 μη ανιχνεύσιμες συστάδες είναι μικρότερες από 2 στρέμματα. Συνεπώς, λιγότερο από 5% της συνολικής επιφάνειας των αμπελώνων δεν ανιχνεύεται ορθά με το φίλτρο Gabor. Γενικά η χρωματική αντίθεση των εικόνων που χρησιμοποιούνται αλλά και η εσωτερική δομή των αμπελώνων αναλογικά με τις διαστάσεις των αμπελιών και του τρόπου που είναι χωροθετημένα, βελτιώνει ή όχι, τα αποτελέσματα αυτόματης ανίχνευσής τους με την παραπάνω μεθοδολογία, προσθέτοντας σφάλματα στη μέθοδο.

Εικόνα 4 :Φωτοερμηνευτική και αυτόματη κατάτμηση εικόνας [3]πηγή



Συμπεράσματα

Η προτεινόμενη αυτοματοποιημένη προσέγγιση, έχει αποδείξει την αποδοτικότητά της στην ανίχνευση και το χαρακτηρισμό των αμπελώνων. Βασικό της πλεονέκτημα, αποτελεί το γεγονός ότι τα αποτελέσματα είναι άμεσα διαθέσιμα σε μια μορφή πολυγωνική (vector αρχεία). Ένα άλλο πρακτικό πλεονέκτημα είναι ότι τα ορθά αποτελέσματα επιτυγχάνονται στην περιοχή του ηλεκτρομαγνητικού φάσματος, όπου βρίσκεται το κόκκινο κανάλι (περίπου στα 0.65μm) στις διαθέσιμες εικόνες με τα φυσικά χρώματα. Επιπλέον, απ τη στιγμή που η κατάλληλη χωρική ανάλυση συνδέεται με την περιοδικότητα των σχημάτων των αμπελώνων, το πιο συχνά επαναλαμβανόμενο πρότυπο που δημιουργείται, θα μπορούσε να χρησιμοποιηθεί σε πολλές διαφορετικές περιοχές αμπελοκαλλιέργειας, οι οποίες παρουσιάζουν ανάλογες εμφανίσεις (όπως για παράδειγμα στην Ισπανία). Προφανώς αντί για αεροφωτογραφίες, θα μπορούσαν να χρησιμοποιηθούν οι δορυφορικές απεικονίσεις Ikonos [11] και Quickbird [12]. Ένα όριο αυτής της μεθόδου, όπως εφαρμόστηκε, είναι ότι έχει αναπτυχθεί για να εφαρμοστεί σε γραμμικά πρότυπα αμπελώνων. Οι επακόλουθες τροποποιήσεις θα ήταν απαραίτητες για να χρησιμοποιήσουν, σε διαδικασία φιλτραρίσματος μέσω του φίλτρου Gabor, στις συστάδες των αμπελώνων όπως εμφανίζονται σε ορισμένες περιοχές με σημαντικές κλίσεις στο έδαφος (Πορτογαλία, Ιταλία).
Εικόνα 5 :Ποσοστό αποτελεσμάτων ανά κατηγορία [4]πηγή
Προσωπικά εργαλεία