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

[[category:Ψηφιακή Τηλεπισκόπηση]]

<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>
Πηγή : 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>
<b>Περίληψη</b><br><br>
Η χαρτογράφηση χρήσεων γης και ιδιαίτερα ο εντοπισμός και η αναγνώριση των αμπελώνων αποτελεί ιδιαίτερα χρήσιμο εργαλείο σε περιοχές που κυριαρχούν τέτοιου είδους καλλιέργειες, όπως για παράδειγμα στη Γαλλία. Σε αυτό το πλαίσιο, η διαθεσιμότητα ενός αυτόματου εργαλείου για την ανίχνευση των αμπελώνων θα ήταν πολύ χρήσιμο. Λόγω των επαναλαμβανόμενων (περιοδικών) σχημάτων αυτών των καλλιεργιών, η ανάλυση συχνότητας φαίνεται να είναι ένα πολύτιμο εργαλείο για την ανίχνευσή τους στις εναέριες εικόνες. Αναπτύχθηκε μία επαναλαμβανόμενη διαδικασία που χρησιμοποιεί τον αλγόριθμο μετασχηματισμού κατά Fourier [http://en.wikipedia.org/wiki/Fourier_transform] ώστε να ικανοποιήσει αυτή την ανάγκη. Αυτό οδηγεί στον ευκρινή χωρισμό των αμπελώνων κάνοντας διακριτά τα όριά τους (της πολυγωνικής μορφής τους) και έχοντας ακριβή εκτίμηση του πλάτους και του προσανατολισμού των σειρών των αμπελιών. Τελικά πραγματοποιείται η ανίχνευση, περισσότερων από 84% των επιφανειών αυτών, από δεδομένα πολύ υψηλών χωρικών συχνοτήτων.<br><br>

<b>Η περιοχή μελέτης</b><br><br>
Τα πειράματα πραγματοποιήθηκαν σε μία έκταση 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"Α]]

<br><br>
<b>Μεθοδολογία</b><br><br>
Η παρούσα μελέτη είναι βασισμένη στις τεχνικές επεξεργασίας εικόνας που εξετάζουν αποκλειστικά τη χωρική δομή των αμπελώνων στις εναέριες εικόνες και όχι τις ραδιομετρικές ιδιότητές τους. Οι τεχνικές αυτές χρησιμοποιούν εικόνες τόνων του γκρι, με την απαίτηση (προϋπόθεση) ικανοποιητικής χρωματικής αντίθεσης μεταξύ των σειρών των αμπελώνων. Επιπλέον παρουσιάζεται ένα βήμα τυποποίησης της διαδικασίας, ώστε τα αποτελέσματα να καταστούν ανεξάρτητα (όσο αυτό είναι δυνατόν) από το αρχικό επίπεδο χρωματικής αντίθεσης. Ως αποτέλεσμα, η προτεινόμενη μέθοδος μπορεί να εφαρμοστεί σε διαφορετικούς τύπους πηγών εικόνας. Ανάλογα με τις φασματικές διαθέσιμες ζώνες σε ένα δεδομένο πλαίσιο εφαρμογής, οι τόνοι του γκρι, μπορούν να χρησιμοποιηθούν για παράδειγμα σε έναν υπολογισμό NDVI [http://en.wikipedia.org/wiki/Normalized_Difference_Vegetation_Index](ομαλοποιημένος δείκτης διαφοράς βλάστησης, συνδυάζοντας τις κόκκινες και εγγύς υπέρυθρες  περιοχές του ηλεκτρομαγνητικού φάσματος) ή στον υπολογισμό της φωτεινότητας ων εικονοστοιχείων για μια ιδιαίτερη υποπεριοχή μιας εικόνας. 
Η βασική ιδέα είναι να απομονωθεί κάθε μεμονωμένη εμφάνιση μιας υποπεριοχής της εικόνας (θα την ονομάζω πλοκή από δω και πέρα), με την επιλογή των αντίστοιχων συχνοτήτων στο φάσμα 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_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]πηγή]]

Για την επικύρωση των αποτελεσμάτων, έγιναν επίγειοι έλεγχοι και συλλέχθηκαν πληροφορίες όπως η κάλυψη του εδάφους, η σύσταση του χώματος, οι διαστάσεις των αμπελώνων και του τρόπου χωροθέτησής τους μέσα στο κάθε αμπέλι κ.α. Οι 160 συστάδες αμπελώνων έχουν εισαχθεί σε μια βάση δεδομένων GIS. Οι επικυρωμένες αυτές συστάδες θα θεωρηθούν ως πλοκές στη διαδικασία επικύρωσης του αμπελιού (με τον τρόπο που αναφέρθηκε παραπάνω). Ο κύριος στόχος της αυτόματης ανίχνευσης αμπελιών, είναι ο προσδιορισμός των ορίων τους. Με τη διαδικασία που αναφέρθηκε, έχει εκτελεστεί μια αυτοματοποιημένη επικύρωση συστάδων οι οποίες έχουν τα χαρακτηριστικά των αμπελώνων, συγκρίνοντας τα αποτελέσματα της αυτοματοποιημένης διαδικασίας και με φωτοερμηνεία, όπου ανάλογα γίνεται η επιβεβαίωσή τους (σε αμπελώνες ή όχι).  Έχουν καθοριστεί κάποιες κατηγορίες στο τελικό αποτέλεσμα, όπως καλή ή μέτρια κατάτμηση, ανάλογα το ποσοστό εμφάνισης (ανίχνευσης) της πλοκής (με τον τρόπο που αναφέρθηκε παραπάνω) στην τελική εικόνα, μετά την εφαρμογή του φίλτρου.<br><br>
<b>Αποτελέσματα</b><br><br>
Οι τρεις πιο αντιπροσωπευτικοί τύποι αποτελεσμάτων είναι οι ορθά ανιχνευμένες συστάδες (πλοκές) με ποσοστό 64%,  οι μερικώς εντοπισμένες με ποσοστό 15,8% και οι μη εντοπισμένες συστάδες, με ποσοστό 11,4%.  Οι μερικώς εντοπισμένες, αντιστοιχούν στην ομαδοποίηση από τις γειτονικές συστάδες που έχουν τον ίδιο προσανατολισμό σειρών και πλάτους και διαχωρίζονται μόνο από ένα στενό δρόμο ή τάφρο. Μερικοί από αυτούς δεν είναι χωρικά διαχωρισμένοι και διαφέρουν μόνο στη σύσταση χώματος της επιφάνειας μεταξύ των σειρών η από μερικά χαρακτηριστικά δυσδιάκριτα στις εναέριες εικόνες, όπως η ηλικία ή το ύψος των αμπελώνων. Η αναλογία έλλειψης ανίχνευσης είναι σχετικά μικρή (περίπου 11,4%) και αφορά κυρίως τις μικρές συστάδες. Πράγματι 9 στις 13 μη ανιχνεύσιμες συστάδες είναι μικρότερες από 2 στρέμματα. Συνεπώς, λιγότερο από 5% της συνολικής επιφάνειας των αμπελώνων δεν ανιχνεύεται ορθά με το φίλτρο Gabor. Γενικά η χρωματική αντίθεση των εικόνων που χρησιμοποιούνται αλλά και η εσωτερική δομή των αμπελώνων αναλογικά με τις διαστάσεις των αμπελιών και του τρόπου που είναι χωροθετημένα, βελτιώνει ή όχι, τα αποτελέσματα αυτόματης ανίχνευσής τους με την παραπάνω μεθοδολογία, προσθέτοντας σφάλματα στη μέθοδο.
[[Εικόνα: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>
<b>Συμπεράσματα </b><br><br>
Η προτεινόμενη αυτοματοποιημένη προσέγγιση, έχει αποδείξει την αποδοτικότητά της στην ανίχνευση και το χαρακτηρισμό των αμπελώνων. Βασικό της πλεονέκτημα, αποτελεί το γεγονός ότι τα αποτελέσματα είναι άμεσα διαθέσιμα σε μια μορφή πολυγωνική (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]πηγή]]