Χρήση εικόνων Landsat 8 για την εκτίμηση του βάθους και της ποιότητας των υδάτων στη "δεξαμενή" νερού El Guajaro στην Κολομβία
Από RemoteSensing Wiki
(Νέα σελίδα με ' category:Αξιολόγηση υδρολογικών αποθεμάτων / δυναμικού <p align="justify"> '''Πρωτότυπος τίτλος:'''Use of LAND...') |
|||
| Γραμμή 38: | Γραμμή 38: | ||
'''6.1 Προεπεξεργασία και Επεξεργασία δορυφορικών εικόνων:''' <br> | '''6.1 Προεπεξεργασία και Επεξεργασία δορυφορικών εικόνων:''' <br> | ||
Η εικόνες Landsat 8 υπέστησαν ραδιομετρική και ατμοσφαιρική διόρθωση. Κατά την ραδιομετρική διόρθωση η σχετικές τιμές των pixel ή οι ψηφιακές τους τιμές μετασχηματίζονται σε μετρήσεις ακτινοβολίας ανά μονάδα μήκους κύματος του φωτός ή σε ανακλαστικότητα. Κατά την ατμοσφαιρική διόρθωση περιορίζεται η επίδραση των ατμοσφαιρικών παραγόντων και οι ραδιομετρικές τιμές μετασχηματίζονται σε ακτινοβολία ή ανακλαστικότητα εδάφους. Με τον τρόπο αυτό οι μετρήσεις των βιοφυσικών παραμέτρων είναι δυνατόν να εκτιμηθούν και να συγκριθούν στο χώρο και στον χρόνο (Chavez, 1988). Στην παρούσα έρευνα η ατμοσφαιρική διόρθωση εφαρμόστηκε στα κανάλια του ορατού φάσματος στα κανάλια NIR και SWIR1, SWIR2 μέσω της μεθόδου αφαίρεσης του μέλανος σώματος (Chavez, 1988). Με βάση τη μέθοδο αυτή θεωρείται ότι μερικά από τα πιο σκούρα αντικείμενα έχουν τιμές ανακλαστικότητας κοντά στο 0, αλλά λόγω της διασποράς και της ατμοσφαιρικής απορρόφησης, τιμές ανακλαστικότητας διάφορες του 0 κατατάσσονται στα pixel των σκούρων αντικειμένων και οι τιμές αυτές πρέπει να αφαιρεθούν απο τα διάφορα φασματικά κανάλια της εικόνας. Για την επεξεργασία χρησιμοποιήθηκε το λογισμικό ENVI και η έκδοση 5.2 . | Η εικόνες Landsat 8 υπέστησαν ραδιομετρική και ατμοσφαιρική διόρθωση. Κατά την ραδιομετρική διόρθωση η σχετικές τιμές των pixel ή οι ψηφιακές τους τιμές μετασχηματίζονται σε μετρήσεις ακτινοβολίας ανά μονάδα μήκους κύματος του φωτός ή σε ανακλαστικότητα. Κατά την ατμοσφαιρική διόρθωση περιορίζεται η επίδραση των ατμοσφαιρικών παραγόντων και οι ραδιομετρικές τιμές μετασχηματίζονται σε ακτινοβολία ή ανακλαστικότητα εδάφους. Με τον τρόπο αυτό οι μετρήσεις των βιοφυσικών παραμέτρων είναι δυνατόν να εκτιμηθούν και να συγκριθούν στο χώρο και στον χρόνο (Chavez, 1988). Στην παρούσα έρευνα η ατμοσφαιρική διόρθωση εφαρμόστηκε στα κανάλια του ορατού φάσματος στα κανάλια NIR και SWIR1, SWIR2 μέσω της μεθόδου αφαίρεσης του μέλανος σώματος (Chavez, 1988). Με βάση τη μέθοδο αυτή θεωρείται ότι μερικά από τα πιο σκούρα αντικείμενα έχουν τιμές ανακλαστικότητας κοντά στο 0, αλλά λόγω της διασποράς και της ατμοσφαιρικής απορρόφησης, τιμές ανακλαστικότητας διάφορες του 0 κατατάσσονται στα pixel των σκούρων αντικειμένων και οι τιμές αυτές πρέπει να αφαιρεθούν απο τα διάφορα φασματικά κανάλια της εικόνας. Για την επεξεργασία χρησιμοποιήθηκε το λογισμικό ENVI και η έκδοση 5.2 . | ||
| - | Οι υδάτινες περιοχές οριοθετήθηκαν με τη χρήση του δείκτη Normalized Differential Water Index (NDWI) (McFeeters, 1996; Khattab and Merkel, 2013). Ο δείκτης ορίζεται ως εξής:<b> | + | Οι υδάτινες περιοχές οριοθετήθηκαν με τη χρήση του δείκτη Normalized Differential Water Index (NDWI) (McFeeters, 1996; Khattab and Merkel, 2013). Ο δείκτης ορίζεται ως εξής:<br> |
| + | </p> | ||
| + | <p align="center"> | ||
| + | <b> NDWI=(b3-NIR)/(b3+NIR) </b> | ||
</p> | </p> | ||
| - | |||
<p align="justify"> | <p align="justify"> | ||
| - | + | ||
, όπου b3 το πράσινο κανάλι και NIR το κανάλι του εγγύος υπέρυθρου (b5). | , όπου b3 το πράσινο κανάλι και NIR το κανάλι του εγγύος υπέρυθρου (b5). | ||
Για την επεξεργασία χρησιμοποιήθηκε το λογισμικό ENVI. Δημιουργήθηκε μάσκα στην οποία στα εικονοστοιχεία που αντιστοιχούσαν στην επιφάνεια του νερού αποδόθηκε η τιμή 1, ενώ στα υπόλοιπα(λ.χ. Βλάστηση, έδαφος) η τιμή 0. Στη συνέχεια η μάσκα αυτή πολλαπλασιάστηκε με τις εικόνες Landsat 8 που έχουν διορθωθεί ραδιομετρικά και ατμοσφαιρικά. <br> | Για την επεξεργασία χρησιμοποιήθηκε το λογισμικό ENVI. Δημιουργήθηκε μάσκα στην οποία στα εικονοστοιχεία που αντιστοιχούσαν στην επιφάνεια του νερού αποδόθηκε η τιμή 1, ενώ στα υπόλοιπα(λ.χ. Βλάστηση, έδαφος) η τιμή 0. Στη συνέχεια η μάσκα αυτή πολλαπλασιάστηκε με τις εικόνες Landsat 8 που έχουν διορθωθεί ραδιομετρικά και ατμοσφαιρικά. <br> | ||
Αναθεώρηση της 09:19, 20 Ιανουαρίου 2021
Πρωτότυπος τίτλος:Use of LANDSAT 8 images for depth and water quality assessment of El Guajaro reservoir, Colombia
Συγγραφείς:Luis Carlos González-Márquez, Franklin M. Torres-Bejarano, Ana Carolina Torregroza-Espinosa, Ivette Renée Hansen-Rodríguez, Hugo B. Rodríguez-Gallegos
Citation:Luis Carlos González-Márquez, Franklin M. Torres-Bejarano, Ana Carolina Torregroza-Espinosa, Ivette Renée Hansen-Rodríguez, Hugo B. Rodríguez-Gallegos, Use of LANDSAT 8 images for depth and water quality assessment of El Guájaro reservoir, Colombia, Journal of South American Earth Sciences, Volume 82, 2018, Pages 231-238, ISSN 0895-9811, [1].
Πηγή:ScienceDirect
Λέξεις κλειδιά:Water quality; Remote sensing; Statistical models; Multispectral analysis
Αντικείμενο εγαρμογής: Αξιλόγηση ποιότητας και χαρακτηριστικών υδάτων.
1. Σκοπός και χρήση της εφαρμογής:
Η παρούσα έρευνα στοχεύει στην ανάπτυξη μεθοδολογίας για την εκτίμηση του βάθους και την ποιότητα των υδάτων συνδυάζοντας δορυφορικές εικόνες, επιτόπιες μετρήσεις και στατιστικά μοντέλα.
2. Εισαγωγή:
Το κλίμα, οι ορογενετικές διαδικασίες, καθώς και η βιο-γεωγραφική θέση της Κολομβίας συνέβαλαν στην αφθονία υδάτων στην χώρα(Minambiente, 2010; Castellanos, 2015). Η εθνική πολιτική για την ολοκληρωμένη διαχείριση των υδάτων επισημαίνει την έλλειψη περιοδικής και συστηματικής μελέτης των υδάτων. Ταυτόχρονα, είναι προφανές πως η ποιότητα των υδάτων μεταβάλλεται, επηρεάζοντας τόσο τη διαθεσιμότητα, όσο και την καταλληλότητα για ορισμένες χρήσεις, με διαφορετική ένταση ανά περιοχή και υδάτινο σώμα. Συνεπώς, αναδεικνύεται η αναγκαιότητα για πρόσβαση σε επαρκούς κάλυψης, συνεκτικότητας και ανάλυσης πληροφορίες για την πιο ολοκληρωμένη χωροχρονική μελέτη των υδάτων (Minambiente, 2010). Συνολικά, η πρόσβαση σε πληροφορίες σε σχέση με την ποιότητα των υδάτων στην Κολομβία και ειδικά στην περιοχή της Καραϊβικής χαρακτηρίζεται από δυσκολίες, ενώ ταυτόχρονα οι διαθέσιμες πληροφορίες είναι περιορισμένες και κατακερματισμένες.
3. Ο ρόλος της τηλεπισκόπησης:
Παραδοσιακά η μελέτη των υδάτινων σωμάτων βασίζεται στην επιτόπια δειγματοληψία και την ανάλυση των δειγμάτων σε εργαστηριακό περιβάλλον, πρακτική η οποία είναι χρονοβόρα και κοστοβόρα και δεν ενδείκνυται για την παρατήρηση των χωροχρονικών μεταβολών σε μεγάλες περιοχές (He et al., 2008). Ακόμη, δεν διευκολύνεται η τακτική παρατήρηση των υδάτων, η οποία είναι ιδιαίτερα σημαντική για την ολοκληρωμένη διαχείριση τους (Wang and Ma, 2001). Δεδομένου ότι η ορατότητα του νερού επηρεάζεται από παράγοντες που σχετίζονται με την ποιότητα του (Pavelsky and Smith, 2009), αναπτύχθηκαν αλγόριθμοι για την παρακολούθηση παραμέτρων της ποιότητας των υδάτων με εργαλεία τηλεπισκόπησης (He et al., 2008; El Saadi et al., 2014) Η τηλεπισκόπηση επιτρέπει την παρακολούθηση εκτεταμένων περιοχών, χωρίς τα υψηλά κόστη της επιτόπιας έρευνας (Hadjimitsis et al., 2010). Εμπόδιο στην μελέτη της ποιότητας των υδάτων με εργαλεία τηλεπισκόπησης μπορεί να σταθεί η επίδραση των ατμοσφαιρικών παραγόντων, αν και έχουν αναπτυχθεί μεθοδολογίες για την απαλοιφή της . Σε γενικές γραμμές, μια σειρά από έρευνες έχουν διεξαχθεί με χρήση εργαλείων τηλεπισκόπησης για την μελέτη παραμέτρων της ποιότητας των υδάτων, όπως η συγκέντρωση χλωροφύλλης, η διαπερατότητα, η θολότητα κ.ο.κ. (Ritchie et al., 2003; Hellweger et al., 2004; Dona et al., 2015; Harvey et al., 2015), αλλά και για τον προσδιορισμό θέσεων υδατοκαλλιεργειών με χρήση πολυφασματικών ή υπερφασματικών εικόνων (Rajitha et al., 2007; Alexandridis et al., 2008; Abd-Elrahman et al., 2011).
4. Περιοχή μελέτης:
Ως περιοχή μελέτης επιλέγεται η "δεξαμενή" El Guajaro (Εικόνα 1) στην Κολομβία. Στην αρχική της κατάσταση, η λίμνη χαρακτηρίζονταν από χωρητικότητα ίση με 400,000,000 m3 και έκταση ίση με 16,000 ha και μέσο βάθος 5 μέτρων. Σήμερα, η έκταση της λίμνης υπολογίζεται σε 11,647 ha, η περίμετρο της σε 114.28 km και πραγματικό όγκο ίσο με 240,000,000 m3 (IDEAM, 2015). Επίσης, η λίμνη διαθέτει δύο διαύλους επικοινωνίας με το κανάλι του φράγματος για τον έλεγχο του επιπέδου της λίμνης (Torres Bejarano et al., 2015). Τα νερά της "δεξαμενή"ς αξιοποιούνται ποικιλοτρόπως: για ανθρώπινη κατανάλωση, για άρδευση γεωργικών προϊόντων και καλλιέργειες ψαριών κ.ο.κ.
5. Δεδομένα
5.1 Δεδομένα πεδίου:
Πραγματοποιήθηκαν δύο συμπληρωματικές δειγματοληψίες στην περιοχή, πριν την περίοδο των βροχοπτώσεων, στις 28/02/2015 (Δείγμα 1) και στις 16/03/2015 (Δείγμα 2) και οι ημερομηνίες επιλέχθηκαν έτσι ώστε να συμπίπτουν με την ημερομηνία λήψης των εικόνων Landsat 8. Τα δείγματα, προέρχονται από όλη την έκταση της λίμνης για να είναι αντιπροσωπευτικά, ενώ εξαιρέθηκαν οι περιοχές βλάστησης στην επιφάνεια της λίμνης (Εικόνα 1). Μετρήθηκαν οι θερμοκρασία, η ηλεκτρική αγωγιμότητα (ΕΟ), το διαλυμένο οξυγόνο (DO), το pH, η θολότητα και το βάθος.
5.2 Δορυφορικά δεδομένα:
Χρησιμοποιήθηκαν δύο εικόνες Landsat 8 (path, 9; row, 53) με ημερομηνίες λήψης 28/02/2015 και 16/03/2015. Τα δεδομένα αντλήθηκαν από την ιστοσελίδα [2]. Η νεφοκάλυψη και για τις δύο εικόνες ανέρχεται σε 3%. Ο τύπος των εικόνων είναι level 1 GeoTIFF. Στον δορυφόρο είναι προσαρμοσμένοι δύο αισθητήρες ο OLI και ο ΤIRS. Οι εικόνες διαθέτουν 11 κανάλια διαφορετικής χωρικής διακριτικής ικανότητας. Tα κανάλια (bands) b1 (coast/aerosol), b2 (blue), b3 (green), b4 (red), b5 (near infrared, NIR), b6 and b7 (short wave infrared bands SWIR1, and SWIR2), όπως και το κανάλι b9 (cirrus) έχουν χωρική διακριτική ικανότητα της τάξης των 30m , το κανάλι b8 (panchromatic) 15 m και τα κανάλια b10 και b11 (TIR-1 and TIR-2) 100 m, η οποία επαναπροσδιορίζεται στα 30m.
6. Μέθοδος:
Η μέθοδος που χρησιμοποιήθηκε στην παρούσα έρευνα παρουσιάζεται σχηματικά στην Εικόνα 2.
6.1 Προεπεξεργασία και Επεξεργασία δορυφορικών εικόνων:
Η εικόνες Landsat 8 υπέστησαν ραδιομετρική και ατμοσφαιρική διόρθωση. Κατά την ραδιομετρική διόρθωση η σχετικές τιμές των pixel ή οι ψηφιακές τους τιμές μετασχηματίζονται σε μετρήσεις ακτινοβολίας ανά μονάδα μήκους κύματος του φωτός ή σε ανακλαστικότητα. Κατά την ατμοσφαιρική διόρθωση περιορίζεται η επίδραση των ατμοσφαιρικών παραγόντων και οι ραδιομετρικές τιμές μετασχηματίζονται σε ακτινοβολία ή ανακλαστικότητα εδάφους. Με τον τρόπο αυτό οι μετρήσεις των βιοφυσικών παραμέτρων είναι δυνατόν να εκτιμηθούν και να συγκριθούν στο χώρο και στον χρόνο (Chavez, 1988). Στην παρούσα έρευνα η ατμοσφαιρική διόρθωση εφαρμόστηκε στα κανάλια του ορατού φάσματος στα κανάλια NIR και SWIR1, SWIR2 μέσω της μεθόδου αφαίρεσης του μέλανος σώματος (Chavez, 1988). Με βάση τη μέθοδο αυτή θεωρείται ότι μερικά από τα πιο σκούρα αντικείμενα έχουν τιμές ανακλαστικότητας κοντά στο 0, αλλά λόγω της διασποράς και της ατμοσφαιρικής απορρόφησης, τιμές ανακλαστικότητας διάφορες του 0 κατατάσσονται στα pixel των σκούρων αντικειμένων και οι τιμές αυτές πρέπει να αφαιρεθούν απο τα διάφορα φασματικά κανάλια της εικόνας. Για την επεξεργασία χρησιμοποιήθηκε το λογισμικό ENVI και η έκδοση 5.2 .
Οι υδάτινες περιοχές οριοθετήθηκαν με τη χρήση του δείκτη Normalized Differential Water Index (NDWI) (McFeeters, 1996; Khattab and Merkel, 2013). Ο δείκτης ορίζεται ως εξής:
NDWI=(b3-NIR)/(b3+NIR)
, όπου b3 το πράσινο κανάλι και NIR το κανάλι του εγγύος υπέρυθρου (b5).
Για την επεξεργασία χρησιμοποιήθηκε το λογισμικό ENVI. Δημιουργήθηκε μάσκα στην οποία στα εικονοστοιχεία που αντιστοιχούσαν στην επιφάνεια του νερού αποδόθηκε η τιμή 1, ενώ στα υπόλοιπα(λ.χ. Βλάστηση, έδαφος) η τιμή 0. Στη συνέχεια η μάσκα αυτή πολλαπλασιάστηκε με τις εικόνες Landsat 8 που έχουν διορθωθεί ραδιομετρικά και ατμοσφαιρικά.
6.2 Στατιστική ανάλυση:
Για τον προσδιορισμό των στατιστικών μοντέλων, τα αποτελέσματα της δεύτερη δειγματοληψίας χρησιμοποιήθηκαν ως εξαρτημένες μεταβλητές και τα δεδομένα για την ανακλαστικότητα των διάφορων συνδυασμών φασματικών καναλιών των επεξεργασμένων Landsat 8 εικόνων ως ανεξάρτητες μεταβλητές. Η κανονικότητα των δεδομένων ελέγχθηκε με τον έλεγχο e Shapiro-Wilk test, χρησιμοποιώντας το λογισμικό Real Statistics Resource Pack (Zaiontz, 2015), ενώ όσα δεδομένα δεν ακολουθούσαν κανονική κατανομή κανονικοποιήθηκαν . Για την παραγωγή και την επαλήθευση των μοντέλων χρησιμοποιήθηκαν τα δεδομένα των δύο δειγματοληψιών. Τα μοντέλα παράχθηκαν με βηματική παλινδρόμηση και η επαλήθευση πραγματοποιήθηκε με απλή γραμμική παλινδρόμηση σε περιβάλλον Matlab (version 2015).
6.3 Χαρτογράφηση παραμέτρων:
Τέλος, τα μοντέλα αναπαράχθηκαν με το Band Math tool of ENVI και οι τιμές ανακλαστικότητας αντιστοιχήθηκαν σε τιμές ΕC, Turbidity (θολότητα), pH, DO και depth (βάθος). Η χαρτογράφηση των παραμέτρων πραγματοποιήθηκε στο λογισμικό ArcMap 10.2.
7. Αποτελέσματα:
