Γρήγορη ταξινόμηση τηλεπισκοπικών εικόνων και παραγωγή έξυπνων θεματικών χαρτών

Από RemoteSensing Wiki

(Διαφορές μεταξύ αναθεωρήσεων)
Μετάβαση σε: πλοήγηση, αναζήτηση
(Νέα σελίδα με ''''Αξιολόγηση της ακρίβειας των μεθόδων εξαγωγής ακτογραμμών στην τηλεπισκόπηση: μια έξυπνη δ...')
 
(5 ενδιάμεσες αναθεωρήσεις δεν εμφανίζονται.)
Γραμμή 1: Γραμμή 1:
-
'''Αξιολόγηση της ακρίβειας των μεθόδων εξαγωγής ακτογραμμών στην τηλεπισκόπηση: μια έξυπνη διαδικασία για εικόνες sentinel-2'''
+
'''Γρήγορη ταξινόμηση τηλεπισκοπικών εικόνων και παραγωγή έξυπνων θεματικών χαρτών'''
-
'''Πρωτότυπος τίτλος:''' ''Accuracy evaluation of coastline extraction methods in remote sensing: a smart procedure for sentinel-2 images''
+
'''Πρωτότυπος τίτλος:''' ''Remotely sensed image fast classification and smart thematic map production''
-
'''Συγγραφείς: ''' E.Alcaras , P. P. Amoroso , F. G. Figliomeni  , C. Parente , G. Prezioso
+
'''Συγγραφείς: ''' Emanuele Alcaras , Pier Paolo Amoroso, Claudio Parente , Giuseppina Prezioso
-
'''Δημοσιεύθηκε: ''' ''The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Volume XLVIII-4/W3-2022''
+
'''Δημοσιεύθηκε: ''' ''23 Δεκεμβρίου 2021 - The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Volume XLVI-4/W5-2021, 43–50, 2021''
-
'''Σύνδεσμος πρωτότυπου κειμένου:''' [https://doi.org/10.5194/isprs-archives-XLVIII-4-W3-2022-13-2022]
+
'''Σύνδεσμος πρωτότυπου κειμένου:''' [https://doi.org/10.5194/isprs-archives-XLVI-4-W5-2021-43-2021]
-
'''Λέξεις-Κλειδιά: ''' Μη επιβλεπόμενη ταξινόμηση, K-Means, Sentinel-2, εξαγωγή ακτογραμμής, GIS, DEM.
+
'''Λέξεις-Κλειδιά: ''' Thematic Map, Landsat 8 OLI, GIS, Image Classification, Smart Applications, Vegetation Index.
-
'''1. ΕΙΣΑΓΩΓΗ'''[[Εικόνα:RSpaper2.1_SIOMOS.png | thumb| right|'''Εικόνα 1.''' Η περιοχή μελέτης πάνω, εντοπισμός της περιοχής μελέτης σε ισορθογώνια  προβολή και γεωγραφικές συντεταγμένες WGS84 (EPSG: 4326)- κάτω, απεικόνιση σε σύνθεση RGB των εικόνων Sentinel-2 σε συντεταγμένες επιπέδου UTM/WGS 84 εκφρασμένες σε μέτρα (EPSG: 32632).]]
+
'''1. ΕΙΣΑΓΩΓΗ'''[[Εικόνα:RSpaper3.1_image_SIOMOS.png | thumb| right|'''Εικόνα 1.''' Γεωεντοπισμός της περιοχής μελέτης στο βόρειο τμήμα της Σικελίας (κόκκινο ορθογώνιο), σε ισόπλευρη προβολή και γεωγραφικές συντεταγμένες WGS84.]][[Εικόνα:RSpaper3.2_image_SIOMOS.png | thumb| right|'''Εικόνα 2.''' Φυσικό έγχρωμο σύνθετο RGB των εικόνων Landsat 8 OLI, σε συντεταγμένες επιπέδου της προβολής UTM-WGS84.]][[Εικόνα:RSpaper3.1_Table_SIOMOS.png | thumb| right|'''Πίνακας 1.''' Κύρια χαρακτηριστικά των εικόνων Landsat 8 OLI.]][[Εικόνα:RSpaper3.3_image_SIOMOS.png | thumb| right|'''Εικόνα 3.''' Θεματικοί χάρτες με 3 κατηγορίες (νερό, έδαφος, βλάστηση), όπως προέκυψαν από πάνω προς τα κάτω με εφαρμογή αντίστοιχα των NDVI KM, RN KM, RGBN MD και 8MS ML.]][[Εικόνα:RSpaper3.1_Table_SIOMOS.png | thumb| right|'''Πίνακας 2.''' Δείκτες θεματικής ακρίβειας που εφαρμόστηκαν στις μεθόδους ταξινόμησης.]][[Εικόνα:RSpaper3.4_image_SIOMOS.png | thumb| right|'''Εικόνα 4.''' Απεικόνιση του θεματικού χάρτη από τον NDVI KM σε μορφή web σε επικάλυψη με το διανυσματικό αρχείο της Ιταλικής Περιφέρειας.]]
-
Ως ακτογραμμή ορίζεται το όριο μεταξύ ξηράς και θάλασσας και έχει θεμελιώδη σημασία σε πολλούς τομείς, όπως η χαρτογραφία, η παρακολούθηση της παράκτιας γεωμορφολογίας και η θαλάσσια ναυσιπλοΐα. Επειδή αλλάζει συνεχώς σχήμα και θέση λόγω των δυναμικών περιβαλλοντικών συνθηκών, απαιτεί προσεκτική παρακολούθηση η οποία όταν γίνεται με επιτόπιους ελέγχους, είναι χρονοβόρα, δαπανηρή αλλά και δύσκολη πολλές φορές λόγω της μορφολογίας. Η τηλεπισκόπηση με τις διάφορες τεχνικές που έχει αναπτύξει αποτελεί μια έγκυρη εναλλακτική λύση με εξαιρετικά αποτελέσματα, τα οποία βέβαια εξαρτώνται από πολλούς παράγοντες, όπως την εφαρμοζόμενη μέθοδο, τα χαρακτηριστικά του αισθητήρα, την τυπολογία της ακτογραμμής. Για να διαπιστωθεί η ακρίβεια οποιασδήποτε τεχνικής, σε σχέση με αυτούς τους παράγοντες, είναι απαραίτητη  μια ακτογραμμή αναφοράς ως μέτρο σύγκρισης.  Συνήθως, χρησιμοποιούμε τη διανυσματοποιημένη ακτογραμμή που παράγεται χειρωνακτικά από την οπτική φωτοερμηνεία μιας παγχρωματικής εικόνας ή ενός έγχρωμου σύνθετου RGB.
+
Οι εφαρμογές που διατίθενται για smartphone, καθώς και το λογισμικό για συσκευές GNSS/GIS, επιτρέπουν την εύκολη χαρτογράφηση του εντοπισμού και του σχήματος μιας περιοχής με την απόκτηση των συντεταγμένων των κορυφών του περιγράμματός της. Η επιλογή αυτή είναι χρήσιμη για την ταξινόμηση μέσω τηλεπισκόπησης, υποστηρίζοντας την ανίχνευση αντιπροσωπευτικών δειγματοληπτικών περιοχών ενός γνωστού τύπου κάλυψης που θα χρησιμοποιηθούν για την εκπαίδευση αλγορίθμων ή για τη δοκιμή των αποτελεσμάτων ταξινόμησης.
-
Το παρόν άρθρο έχει ως στόχο να καταδείξει, ότι η ακτογραμμή που μπορεί να εξαχθεί πολύ γρήγορα με την χρήση ενός ψηφιακού μοντέλου υψομέτρου, που παράγεται με LiDAR (Lg-DEM) αποτελεί μια αξιόπιστη εναλλακτική προσέγγιση της ακτογραμμής αναφοράς, για τη σύγκριση διαφορετικών μεθόδων στην ανίχνευση της ακτογραμμής σε πολυφασματικές εικόνες Sentinel-2.
+
Η παρούσα εργασία αποσκοπεί στην ανάλυση της δυνατότητας παραγωγής έξυπνων χαρτών από την ταξινόμηση τηλεπισκοπικών εικόνων με γρήγορο τρόπο. Η προσοχή επικεντρώνεται σε διαφορετικές μεθόδους που συγκρίνονται για τον εντοπισμό γρήγορης και ακριβούς διαδικασίας για την παραγωγή ενημερωμένων και αξιόπιστων χαρτών. Οι πολυφασματικές εικόνες Landsat 8 OLI υποβάλλονται σε διάφορους αλγορίθμους ταξινόμησης για τη διάκριση του νερού, του γυμνού εδάφους και της βλάστησης. Ο χάρτης που προκύπτει εισάγεται κατάλληλα σε μια ευρύτερη βάση δεδομένων με στόχο την αναπαράσταση της κατάστασης σε ένα χωροχρονικά εξελικτικό σενάριο και είναι κατάλληλος ώστε να αποτυπωθεί οποιαδήποτε μεταβολή προκαλείται, όπως π.χ. εντοπισμός καμένων περιοχών, καθορισμός περιοχών με βλάστηση για τουριστικούς σκοπούς κ.λπ.
-
Η παραγωγή των ακτογραμμής γίνεται από χάρτες που προκύπτουν από τη χρήση κάποιων δεικτών, στους οποίους εφαρμόζεται η μέθοδος της μη επιβλεπόμενης ταξινόμησης με τον αλγόριθμο  k-Means  ώστε να ταξινομηθεί κάθε προϊόν δείκτη που προκύπτει, σε δύο ομάδες, νερό και χωρίς νερό (δηλαδή νερό και ξηρά). Η αυτόματη διανυσματοποίηση επιτρέπει την ανίχνευση της ακτογραμμής ως διαχωρισμό μεταξύ ξηράς και θάλασσας. Για τον έλεγχο του επιπέδου της ακρίβειας των αποτελεσμάτων, χρησιμοποιήθηκαν δύο διαφορετικές ακτογραμμές αναφοράς, σε διανυσματική μορφή:
+
-
• Η ακτογραμμή που λαμβάνεται από το DEM που παράγεται με δεδομένα 3D LiDAR και με χωρική ανάλυση 2 m x 2 m
+
'''2. ΠΕΡΙΟΧΗ ΜΕΛΕΤΗΣ ΚΑΙ ΔΕΔΟΜΕΝΑ'''
-
• Η ακτογραμμή που εξάγεται χειροκίνητα και βασίζεται στη φωτοερμηνεία RGB.
+
Η περιοχή μελέτης καλύπτει μια επιφάνεια 5.400 Km2 και βρίσκεται στο βόρειο τμήμα της Σικελίας (Ιταλία), όπως φαίνεται στις εικόνες (1) και (2). Η εξεταζόμενη περιοχή περιλαμβάνει πολλά από τα κύρια χαρακτηριστικά του νησιού όσον αφορά την κάλυψη γης, όπως πυκνή αλλά και διάσπαρτη βλάστηση, γυμνό έδαφος καλυμμένο με λάβα, μικρές δασώδεις εκτάσεις, καλλιεργούμενες εκτάσεις με καλή κατάσταση υγείας αλλά και σε γήρανση ,άγονες περιοχές, λίμνες νερού φυσικής και τεχνητής προέλευσης.Οι εικόνες Landsat 8 OLI περιλαμβάνουν 8 πολυφασματικές ζώνες (MS) και μία παγχρωματική ζώνη (PAN), όπως αναφέρεται στον πίνακα (1).
-
'''2. ΠΕΡΙΟΧΗ ΜΕΛΕΤΗΣ ΚΑΙ ΔΕΔΟΜΕΝΑ'''
+
'''3. ΜΕΘΟΔΟΙ'''
-
Τα πειράματα διεξάγονται σε εικόνες Sentinel-2A, που αποκτήθηκαν στις 6 Αυγούστου 2019 και παρουσιάζουν επίπεδο επεξεργασίας 2Α, που αφορούν ένα τμήμα της μητροπολιτικής πόλης της Νάπολης (Ιταλία) με μεγάλη ποικιλία στη μορφολογία της ακτογραμμής (εικόνα 1).Το DEM που χρησιμοποιήθηκε προέρχεται από δεδομένα LiDAR, που αποκτήθηκαν το 2008, έχει κατακόρυφη ακρίβεια μικρότερη από ± 15 cm, η οριζόντια ακρίβεια είναι ± 30 cm ενώ έχει και αρνητικές τιμές που αντιπροσωπεύουν το βυθό κάτω από την ακτή (εικόνα 2). Η ακτογραμμή εξάγεται ως η γραμμή περιγράμματος στο μηδέν. Οι εικόνες του Sentinel-2 παράγονται από 2 δίδυμους δορυφόρους, τους Sentinel-2A και Sentinel-2B σε 13 φασματικές υποπεριοχές της ΗΜΑ και τα κύρια χαρακτηριστικά τους παρουσιάζονται στον πίνακα 1.
+
Ακολουθούν τα βήματα για την παραγωγή ενός έξυπνου χάρτη:
-
[[Εικόνα:RSpaper2.2_SIOMOS.png | thumb| left|'''Εικόνα 2.''' Τρισδιάστατο φωτορεαλιστικό μοντέλο της περιοχής μελέτης που προέκυψε από την εφαρμογή της RGB σύνθεσης πραγματικών χρωμάτων του Sentinel-2 στο Lg-DEM]][[Εικόνα:RSpaper2.1_Table_SIOMOS.png | thumb| center|'''Πίνακας 1.''' Κύρια χαρακτηριστικά των εικόνων Sentinel-2.]]
+
 +
'''3.1 Μετατροπή των τιμών BV σε ανακλαστικότητα στην κορυφή της ατμόσφαιρας (TOA) και ατμοσφαιρική διόρθωση εφαρμόζοντας το μοντέλο DOS'''
 +
'''3.2 Ταξινόμηση των εικόνων'''
 +
•''3.2.1 Μη Επιβλεπόμενη Ταξινόμηση''
 +
Με εφαρμογή της μεθόδου K-means που είναι μια αριθμητική, μη εποπτευόμενη, μη ντετερμινιστική, επαναληπτική μέθοδος.
 +
 +
•''3.2.2 Επιβλεπόμενη Ταξινόμηση''
 +
 +
Στην οποία οι περιοχές εκπαίδευσης εντοπίζονται μέσω οπτικής ανάλυσης στη σύνθεση RGB, καθώς και μέσω της συμβουλευτικής εικόνας υψηλότερης ανάλυσης που αποκτήθηκε την ίδια περίοδο με εφαρμογή των μεθόδων της ελάχιστης απόστασης MD και της μέγιστης πιθανοφάνειας ΜL.
 +
 +
•''3.2.3 Δείκτης κανονικοποιημένης διαφοράς βλάστησης NDVI''
 +
 +
Με την εφαρμογή του NDVI (Normalized Difference Vegetation Index) παράγεται μια συνθετική εικόνα που χρειάζεται την εφαρμογή κριτηρίου ταξινόμησης για να κατηγοροποιηθούν στη συνέχεια εικονοστοιχεία νερού, γυμνού εδάφους και βλάστησης.
 +
 +
'''3.3 Αξιολόγηση της θεματικής ακρίβειας'''
 +
 +
Για να διαπιστωθεί η θεματική ακρίβεια των αποτελεσμάτων που προκύπτουν, επιλέγονται περιοχές δοκιμής, διαφορετικές από τις περιοχές εκπαίδευσης που χρησιμοποιούνται στην επιβλεπόμενη ταξινόμηση, και χρησιμοποιείται ο πίνακας σύγχυσης (Confusion Matrix) ως ποσοτική μέθοδος για την αξιολόγηση της ακρίβειας ταξινόμησης εικόνων, με τις τιμές των παραδοσιακών δεικτών που ονομάζονται Producer Accuracy (PA), User Accuracy (UA) και Overall Accuracy (OA)
 +
 +
'''3.4 Παραγωγή έξυπνων χαρτών για διαδικτυακή εφαρμογή'''
 +
 +
Οι εργασίες για τη δημοσίευση δεδομένων ΓΣΠ στο διαδίκτυο και την προσβασιμότητά τους σε άλλους χρήστες ονομάζονται συνήθως διαδικτυακή χαρτογράφηση. Τα προϊόντα που προκύπτουν είναι επίσης γνωστά ως "έξυπνοι χάρτες", δεδομένου ότι είναι κατάλληλα για έξυπνες εφαρμογές σε tablet και smartphone.
 +
Στην παρούσα μελέτη μετατρέπουμε τους παραγόμενους θεματικούς χάρτες με το πρόσθετο QGIS2Web σε κατάλληλα αρχεία μορφής για να τους καταστήσουμε κατάλληλους για δημοσίευση στο διαδίκτυο. Με αυτόν τον τρόπο μπορούμε να τους χρησιμοποιήσουμε και για τον έλεγχο του επιτευχθέντος επιπέδου θεματικής ακρίβειας στην περίπτωση επιτόπιας παρατήρησης.
 +
 +
'''4. ΑΠΟΤΕΛΕΣΜΑΤΑ ΚΑΙ ΣΥΖΗΤΗΣΗ'''
 +
 +
Στην παρούσα μελέτη εφαρμόζονται οι μέθοδοι ταξινόμησης Minimum Distance, Maximum Likelihood και K-means σε όλες τις εικόνες Landsat 8 OLI ή σε μια επιλογή από αυτές καθώς και στον NDVI- ειδικότερα χρησιμοποιείται ο ακόλουθος συνδυασμός:
 +
 +
- NDVI με ελάχιστη απόσταση (NDVI MD),
 +
 +
- NDVI με μέγιστη πιθανοφάνεια (NDVI ML),
 +
 +
- NDVI με K-Means (NDVI KM),
 +
 +
- RED & NIR με ελάχιστη απόσταση (RN MD),
 +
 +
- RED & NIR με μέγιστη πιθανοφάνεια (RN ML),
 +
 +
- RED & NIR με K-Means (RN KM),
 +
 +
- RGB & NIR με ελάχιστη απόσταση (RGBN MD),
 +
 +
- RGB & NIR με μέγιστη πιθανοφάνεια (RGBN ML),
 +
 +
- RGB & NIR με K-Means (RGBN KM),
 +
 +
- Και οι 8 ζώνες MS με ελάχιστη απόσταση (8MS MD),
 +
 +
- Και οι 8 ζώνες MS με μέγιστη πιθανοφάνεια (8MS ML),
 +
 +
- Και οι 8 ζώνες MS με K-Means (8MS KM).
 +
 +
Οι έξοδοι ταξινομούνται σε τρεις κατηγορίες: νερό, γυμνό έδαφος και βλάστηση. Μια οπτική σύγκριση ορισμένων αποτελεσμάτων ταξινόμησης, συγκεκριμένα εκείνων που λαμβάνονται από τα NDVI KM, RN KM, RGBN MD και 8MS ML, παρουσιάζεται στην εικόνα (3).
 +
Η μέθοδος NDVI KM είναι η πιο αποδοτική μέθοδος, παρέχοντας την υψηλότερη τιμή OA (0,9962) και τις βέλτιστες τιμές PA και UA, αφού είναι ίσες με 1, εκτός από την περίπτωση της PA της βλάστησης (0,9886) και της UA του εδάφους (0,9887).
 +
Ο θεματικός χάρτης που προκύπτει από τον NDVI KM σε μορφή ιστού, ο οποίος παράγεται με τη χρήση του πρόσθετου QGIS2Web, παρουσιάζεται στην εικόνα (4):
 +
 +
'''5. ΣΥΜΠΕΡΑΣΜΑΤΑ'''
 +
 +
Η μελέτη αυτή αποδεικνύει ότι οι θεματικοί διαδικτυακοί χάρτες που επικεντρώνονται στη διάκριση της γήινης επιφάνειας σε νερό, γυμνό έδαφος και βλάστηση μπορούν εύκολα να παραχθούν με αυτόματο τρόπο χρησιμοποιώντας λίγες ζώνες MS των δορυφορικών εικόνων, συγκεκριμένα μόνο δύο ζώνες της εικόνας Landsat 8 OLI, οι οποίες είναι η κόκκινη και η NIR που χρησιμοποιούνται για την παραγωγή του NDVI.
 +
Επίσης αποδεικνύεται ότι η χρήση του k-means, η οποία δεν απαιτεί την ανθρώπινη επίβλεψη, επιτρέπει τη διάκριση των αποτελεσμάτων του NDVI και τη γρήγορη απόδοση ετικετών στις τρεις ονομαστικές κλάσεις διατηρώντας πολύ υψηλό επίπεδο θεματικής ακρίβειας. Επομένως η μη επιβλεπόμενη ταξινόμηση είναι χρήσιμη για τη διάκριση λίγων, μη περίπλοκων και ευρείας κλίμακας κλάσεων κάλυψης γης, όπως το νερό, το γυμνό έδαφος και η βλάστηση και επιπλέον μπορεί να μειώσει τα σφάλματα του χειριστή κατά την αναγνώριση των περιοχών εκπαίδευσης. Τέλος η δημοσίευση στο διαδίκτυο του παραγόμενου θεματικού χάρτη, μια δύσκολη εργασία που συνήθως απαιτεί ειδικές δεξιότητες διαφορετικές από εκείνες που χαρακτηρίζουν τον παραδοσιακό χαρτογράφο, είναι δυνατή με γρήγορο τρόπο με τη χρήση ειδικών εργαλείων όπως το πρόσθετο QGIS2Web.
-
'''3. ΜΕΘΟΔΟΙ'''
 
-
'''3.1 Υπολογισμός του δείκτη'''
 
  [[category:Χαρτογραφία]]
  [[category:Χαρτογραφία]]

Παρούσα αναθεώρηση της 09:13, 12 Φεβρουαρίου 2023

Γρήγορη ταξινόμηση τηλεπισκοπικών εικόνων και παραγωγή έξυπνων θεματικών χαρτών

Πρωτότυπος τίτλος: Remotely sensed image fast classification and smart thematic map production

Συγγραφείς: Emanuele Alcaras , Pier Paolo Amoroso, Claudio Parente , Giuseppina Prezioso

Δημοσιεύθηκε: 23 Δεκεμβρίου 2021 - The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Volume XLVI-4/W5-2021, 43–50, 2021

Σύνδεσμος πρωτότυπου κειμένου: [1]

Λέξεις-Κλειδιά: Thematic Map, Landsat 8 OLI, GIS, Image Classification, Smart Applications, Vegetation Index.


1. ΕΙΣΑΓΩΓΗ
Εικόνα 1. Γεωεντοπισμός της περιοχής μελέτης στο βόρειο τμήμα της Σικελίας (κόκκινο ορθογώνιο), σε ισόπλευρη προβολή και γεωγραφικές συντεταγμένες WGS84.
Εικόνα 2. Φυσικό έγχρωμο σύνθετο RGB των εικόνων Landsat 8 OLI, σε συντεταγμένες επιπέδου της προβολής UTM-WGS84.
Πίνακας 1. Κύρια χαρακτηριστικά των εικόνων Landsat 8 OLI.
Εικόνα 3. Θεματικοί χάρτες με 3 κατηγορίες (νερό, έδαφος, βλάστηση), όπως προέκυψαν από πάνω προς τα κάτω με εφαρμογή αντίστοιχα των NDVI KM, RN KM, RGBN MD και 8MS ML.
Πίνακας 2. Δείκτες θεματικής ακρίβειας που εφαρμόστηκαν στις μεθόδους ταξινόμησης.
Εικόνα 4. Απεικόνιση του θεματικού χάρτη από τον NDVI KM σε μορφή web σε επικάλυψη με το διανυσματικό αρχείο της Ιταλικής Περιφέρειας.

Οι εφαρμογές που διατίθενται για smartphone, καθώς και το λογισμικό για συσκευές GNSS/GIS, επιτρέπουν την εύκολη χαρτογράφηση του εντοπισμού και του σχήματος μιας περιοχής με την απόκτηση των συντεταγμένων των κορυφών του περιγράμματός της. Η επιλογή αυτή είναι χρήσιμη για την ταξινόμηση μέσω τηλεπισκόπησης, υποστηρίζοντας την ανίχνευση αντιπροσωπευτικών δειγματοληπτικών περιοχών ενός γνωστού τύπου κάλυψης που θα χρησιμοποιηθούν για την εκπαίδευση αλγορίθμων ή για τη δοκιμή των αποτελεσμάτων ταξινόμησης. Η παρούσα εργασία αποσκοπεί στην ανάλυση της δυνατότητας παραγωγής έξυπνων χαρτών από την ταξινόμηση τηλεπισκοπικών εικόνων με γρήγορο τρόπο. Η προσοχή επικεντρώνεται σε διαφορετικές μεθόδους που συγκρίνονται για τον εντοπισμό γρήγορης και ακριβούς διαδικασίας για την παραγωγή ενημερωμένων και αξιόπιστων χαρτών. Οι πολυφασματικές εικόνες Landsat 8 OLI υποβάλλονται σε διάφορους αλγορίθμους ταξινόμησης για τη διάκριση του νερού, του γυμνού εδάφους και της βλάστησης. Ο χάρτης που προκύπτει εισάγεται κατάλληλα σε μια ευρύτερη βάση δεδομένων με στόχο την αναπαράσταση της κατάστασης σε ένα χωροχρονικά εξελικτικό σενάριο και είναι κατάλληλος ώστε να αποτυπωθεί οποιαδήποτε μεταβολή προκαλείται, όπως π.χ. εντοπισμός καμένων περιοχών, καθορισμός περιοχών με βλάστηση για τουριστικούς σκοπούς κ.λπ.

2. ΠΕΡΙΟΧΗ ΜΕΛΕΤΗΣ ΚΑΙ ΔΕΔΟΜΕΝΑ

Η περιοχή μελέτης καλύπτει μια επιφάνεια 5.400 Km2 και βρίσκεται στο βόρειο τμήμα της Σικελίας (Ιταλία), όπως φαίνεται στις εικόνες (1) και (2). Η εξεταζόμενη περιοχή περιλαμβάνει πολλά από τα κύρια χαρακτηριστικά του νησιού όσον αφορά την κάλυψη γης, όπως πυκνή αλλά και διάσπαρτη βλάστηση, γυμνό έδαφος καλυμμένο με λάβα, μικρές δασώδεις εκτάσεις, καλλιεργούμενες εκτάσεις με καλή κατάσταση υγείας αλλά και σε γήρανση ,άγονες περιοχές, λίμνες νερού φυσικής και τεχνητής προέλευσης.Οι εικόνες Landsat 8 OLI περιλαμβάνουν 8 πολυφασματικές ζώνες (MS) και μία παγχρωματική ζώνη (PAN), όπως αναφέρεται στον πίνακα (1).

3. ΜΕΘΟΔΟΙ

Ακολουθούν τα βήματα για την παραγωγή ενός έξυπνου χάρτη:

3.1 Μετατροπή των τιμών BV σε ανακλαστικότητα στην κορυφή της ατμόσφαιρας (TOA) και ατμοσφαιρική διόρθωση εφαρμόζοντας το μοντέλο DOS

3.2 Ταξινόμηση των εικόνων

•3.2.1 Μη Επιβλεπόμενη Ταξινόμηση

Με εφαρμογή της μεθόδου K-means που είναι μια αριθμητική, μη εποπτευόμενη, μη ντετερμινιστική, επαναληπτική μέθοδος.

•3.2.2 Επιβλεπόμενη Ταξινόμηση

Στην οποία οι περιοχές εκπαίδευσης εντοπίζονται μέσω οπτικής ανάλυσης στη σύνθεση RGB, καθώς και μέσω της συμβουλευτικής εικόνας υψηλότερης ανάλυσης που αποκτήθηκε την ίδια περίοδο με εφαρμογή των μεθόδων της ελάχιστης απόστασης MD και της μέγιστης πιθανοφάνειας ΜL.

•3.2.3 Δείκτης κανονικοποιημένης διαφοράς βλάστησης NDVI

Με την εφαρμογή του NDVI (Normalized Difference Vegetation Index) παράγεται μια συνθετική εικόνα που χρειάζεται την εφαρμογή κριτηρίου ταξινόμησης για να κατηγοροποιηθούν στη συνέχεια εικονοστοιχεία νερού, γυμνού εδάφους και βλάστησης.

3.3 Αξιολόγηση της θεματικής ακρίβειας

Για να διαπιστωθεί η θεματική ακρίβεια των αποτελεσμάτων που προκύπτουν, επιλέγονται περιοχές δοκιμής, διαφορετικές από τις περιοχές εκπαίδευσης που χρησιμοποιούνται στην επιβλεπόμενη ταξινόμηση, και χρησιμοποιείται ο πίνακας σύγχυσης (Confusion Matrix) ως ποσοτική μέθοδος για την αξιολόγηση της ακρίβειας ταξινόμησης εικόνων, με τις τιμές των παραδοσιακών δεικτών που ονομάζονται Producer Accuracy (PA), User Accuracy (UA) και Overall Accuracy (OA)

3.4 Παραγωγή έξυπνων χαρτών για διαδικτυακή εφαρμογή

Οι εργασίες για τη δημοσίευση δεδομένων ΓΣΠ στο διαδίκτυο και την προσβασιμότητά τους σε άλλους χρήστες ονομάζονται συνήθως διαδικτυακή χαρτογράφηση. Τα προϊόντα που προκύπτουν είναι επίσης γνωστά ως "έξυπνοι χάρτες", δεδομένου ότι είναι κατάλληλα για έξυπνες εφαρμογές σε tablet και smartphone. Στην παρούσα μελέτη μετατρέπουμε τους παραγόμενους θεματικούς χάρτες με το πρόσθετο QGIS2Web σε κατάλληλα αρχεία μορφής για να τους καταστήσουμε κατάλληλους για δημοσίευση στο διαδίκτυο. Με αυτόν τον τρόπο μπορούμε να τους χρησιμοποιήσουμε και για τον έλεγχο του επιτευχθέντος επιπέδου θεματικής ακρίβειας στην περίπτωση επιτόπιας παρατήρησης.

4. ΑΠΟΤΕΛΕΣΜΑΤΑ ΚΑΙ ΣΥΖΗΤΗΣΗ

Στην παρούσα μελέτη εφαρμόζονται οι μέθοδοι ταξινόμησης Minimum Distance, Maximum Likelihood και K-means σε όλες τις εικόνες Landsat 8 OLI ή σε μια επιλογή από αυτές καθώς και στον NDVI- ειδικότερα χρησιμοποιείται ο ακόλουθος συνδυασμός:

- NDVI με ελάχιστη απόσταση (NDVI MD),

- NDVI με μέγιστη πιθανοφάνεια (NDVI ML),

- NDVI με K-Means (NDVI KM),

- RED & NIR με ελάχιστη απόσταση (RN MD),

- RED & NIR με μέγιστη πιθανοφάνεια (RN ML),

- RED & NIR με K-Means (RN KM),

- RGB & NIR με ελάχιστη απόσταση (RGBN MD),

- RGB & NIR με μέγιστη πιθανοφάνεια (RGBN ML),

- RGB & NIR με K-Means (RGBN KM),

- Και οι 8 ζώνες MS με ελάχιστη απόσταση (8MS MD),

- Και οι 8 ζώνες MS με μέγιστη πιθανοφάνεια (8MS ML),

- Και οι 8 ζώνες MS με K-Means (8MS KM).

Οι έξοδοι ταξινομούνται σε τρεις κατηγορίες: νερό, γυμνό έδαφος και βλάστηση. Μια οπτική σύγκριση ορισμένων αποτελεσμάτων ταξινόμησης, συγκεκριμένα εκείνων που λαμβάνονται από τα NDVI KM, RN KM, RGBN MD και 8MS ML, παρουσιάζεται στην εικόνα (3). Η μέθοδος NDVI KM είναι η πιο αποδοτική μέθοδος, παρέχοντας την υψηλότερη τιμή OA (0,9962) και τις βέλτιστες τιμές PA και UA, αφού είναι ίσες με 1, εκτός από την περίπτωση της PA της βλάστησης (0,9886) και της UA του εδάφους (0,9887). Ο θεματικός χάρτης που προκύπτει από τον NDVI KM σε μορφή ιστού, ο οποίος παράγεται με τη χρήση του πρόσθετου QGIS2Web, παρουσιάζεται στην εικόνα (4):

5. ΣΥΜΠΕΡΑΣΜΑΤΑ

Η μελέτη αυτή αποδεικνύει ότι οι θεματικοί διαδικτυακοί χάρτες που επικεντρώνονται στη διάκριση της γήινης επιφάνειας σε νερό, γυμνό έδαφος και βλάστηση μπορούν εύκολα να παραχθούν με αυτόματο τρόπο χρησιμοποιώντας λίγες ζώνες MS των δορυφορικών εικόνων, συγκεκριμένα μόνο δύο ζώνες της εικόνας Landsat 8 OLI, οι οποίες είναι η κόκκινη και η NIR που χρησιμοποιούνται για την παραγωγή του NDVI. Επίσης αποδεικνύεται ότι η χρήση του k-means, η οποία δεν απαιτεί την ανθρώπινη επίβλεψη, επιτρέπει τη διάκριση των αποτελεσμάτων του NDVI και τη γρήγορη απόδοση ετικετών στις τρεις ονομαστικές κλάσεις διατηρώντας πολύ υψηλό επίπεδο θεματικής ακρίβειας. Επομένως η μη επιβλεπόμενη ταξινόμηση είναι χρήσιμη για τη διάκριση λίγων, μη περίπλοκων και ευρείας κλίμακας κλάσεων κάλυψης γης, όπως το νερό, το γυμνό έδαφος και η βλάστηση και επιπλέον μπορεί να μειώσει τα σφάλματα του χειριστή κατά την αναγνώριση των περιοχών εκπαίδευσης. Τέλος η δημοσίευση στο διαδίκτυο του παραγόμενου θεματικού χάρτη, μια δύσκολη εργασία που συνήθως απαιτεί ειδικές δεξιότητες διαφορετικές από εκείνες που χαρακτηρίζουν τον παραδοσιακό χαρτογράφο, είναι δυνατή με γρήγορο τρόπο με τη χρήση ειδικών εργαλείων όπως το πρόσθετο QGIS2Web.

Ανακτήθηκε από το "http://147.102.106.44/rs/wiki/index.php/%CE%93%CF%81%CE%AE%CE%B3%CE%BF%CF%81%CE%B7_%CF%84%CE%B1%CE%BE%CE%B9%CE%BD%CF%8C%CE%BC%CE%B7%CF%83%CE%B7_%CF%84%CE%B7%CE%BB%CE%B5%CF%80%CE%B9%CF%83%CE%BA%CE%BF%CF%80%CE%B9%CE%BA%CF%8E%CE%BD_%CE%B5%CE%B9%CE%BA%CF%8C%CE%BD%CF%89%CE%BD_%CE%BA%CE%B1%CE%B9_%CF%80%CE%B1%CF%81%CE%B1%CE%B3%CF%89%CE%B3%CE%AE_%CE%AD%CE%BE%CF%85%CF%80%CE%BD%CF%89%CE%BD_%CE%B8%CE%B5%CE%BC%CE%B1%CF%84%CE%B9%CE%BA%CF%8E%CE%BD_%CF%87%CE%B1%CF%81%CF%84%CF%8E%CE%BD".