Διαχωρισμός των δομημένου περιοχών από τη γυμνή γη σε μεσογειακές πόλεις χρησιμοποιώντας εικόνες Sentinel-2A.

Από RemoteSensing Wiki

Add Your Content Here

Πρωτότυπος τίτλος: Separating Built-Up Areas from Bare Land in Mediterranean Cities Using Sentinel-2A Imagery

Σχήμα 1. Χάρτης θέσης της κύριας περιοχής μελέτης και των περιοχών δοκιμής (Χάρτης χώρας από ESRI©, Καλιφόρνια, ΗΠΑ, πιο προσεκτική ματιά από το φυσικό σύνθετο χρώμα του Sentinel-2A).

Συγγραφέας: Sajid Pareeth, Poolad Karimi, Mojtaba Shafiei, Charlotte De Fraiture

Σχήμα 2. Φασματικές καμπύλες ανάκλασης διαφορετικών τύπων κάλυψης/χρήσης γης (LCU), σύμφωνα με τις ζώνες εικόνας Sentinel-2A (κορυφαίες τιμές ανάκλασης ατμόσφαιρας που που προέρχονται από δορυφορική εικόνα ανάλυσης 12 bit)

Σχήμα 2. Φασματικές καμπύλες ανάκλασης διαφορετικών τύπων κάλυψης/χρήσης γης (LCU), σύμφωνα με τις ζώνες εικόνας Sentinel-2A (κορυφαίες τιμές ανάκλασης ατμόσφαιρας που που προέρχονται από δορυφορική εικόνα ανάλυσης 12 bit)

Σχήμα 3. Εικόνες του (α) Sentinel-2A RGB; β) ομαλοποιημένος δείκτης δημιουργίας διαφορών (NDBI)· γ) δείκτης δημιουργίας (BUI)· δ) δείκτης εξόρυξης δομημένης έκτασης (BAEI)· ε) νέος δείκτης δημιουργίας (NBI)· στ) δείκτης βλάστησης (VIBI)· ζ) δείκτης που βασίζεται σε ευρετήριο (IBI)· η) αστικός δείκτης (UI)· και i) γυμνό δείκτη εδάφους (BSI).

Σχήμα 4. Οι εικόνες του (α) του Δείκτη Ομαλοποιημένης Διαφοράς (NDTI) και(β) ενός υποσυνόλου της NDTI αποτελούν παράδειγμα της αντίθεσης μεταξύ γυμνών εδαφών και δομημένων περιοχών.

Πίνακας 1. Δείγμα κατάρτισης και κατανομή σημείων αξιολόγησης της ακρίβειας για τις περιφέρειες μελέτης.

Σχήμα 5. Αποτελέσματα ταξινόμησης:α) NDBI· β)BUI· γ)BAEI· δ) NBI· ε)VIBI· στ)IBI· ζ) περιβάλλον εργασίας χρήστη· η)BSI· i)σύνολα πολλαπλών ευρετηρίων που βασίζονται σε NDTI· και(ι)πρωτότυπες φασματικές ζώνες.

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

Εισαγωγή Σε όλη την ιστορία, η αύξηση της πυκνότητας του πληθυσμού και η επέκταση των αστικών περιοχών, ιδίως στις μητροπολιτικές πόλεις, έχουν αλλάξει τη μορφή της επιφάνειας της Γης. Ο ρυθμός αύξησης της κάλυψης/χρήσης γης (LCU) έχει αυξηθεί τις τελευταίες δεκαετίες. Η αύξηση του πληθυσμού οδηγεί σε αύξηση της κατανάλωσης νερού και ενέργειας και προκαλεί αλλαγές στην επιφάνεια της γης, οι οποίες έχουν ως αποτέλεσμα την περιφερειακή έως παγκόσμια κλιματική αλλαγή και την υποβάθμιση του περιβάλλοντος. Οι αλλαγές στην LCU που έχουν συμβεί λόγω της αστικοποίησης, της αποψίλωσης των δασών, των φυσικών καταστροφών και των έντονων γεωργικών πρακτικών, έχουν επηρεάσει σε μεγάλο βαθμό τα κλιματικά χαρακτηριστικά σε περιφερειακή και παγκόσμια κλίμακα. Οι αλλαγές αυτές οδηγούν γενικά σε αυξήσεις των θερμοκρασιών κοντά στην επιφάνεια και σχηματισμό των νήσων θερμότητας, οι οποίες προκαλούν άλλα κλιματικά φαινόμενα. Η διαθεσιμότητα δωρεάν παγκόσμιων και ιστορικών δορυφορικών εικόνων παρέχει μια πολύτιμη ευκαιρία για τη χαρτογράφηση και την παρακολούθηση της LCU. Ο κύριος στόχος αυτού του ερευνητικού έργου ήταν να προτείνει μια προσέγγιση ταξινόμησης SVM που βασίζεται σε πολλαπλούς δείκτες, για τη χαρτογράφηση επτά διαφορετικών κατηγοριών LCU, σε σύνθετες αστικές περιοχές. Η έρευνα επικεντρώθηκε στον διαχωρισμό των κατοικημένων περιοχών και των γυμνών εδαφών, εκτός από την παροχή ενός ακριβούς και αξιόπιστου χάρτη LCU, σε τρεις πυκνά αστικοποιημένες μητροπολιτικές πόλεις της Τουρκίας, χρησιμοποιώντας τις εικόνες Sentinel-2A.

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

Η Κωνσταντινούπολη είναι η πιο πυκνοκατοικημένη και η μεγαλύτερη πόλη της Τουρκίας, που βρίσκεται σε γεωγραφικό πλάτος 41°00′44,06" Β και γεωγραφικό μήκος 28°58′33,66" Ε, στο Βόρειο Ημισφαίριο, ενώνοντας τις δύο ηπείρους της Ασίας και της Ευρώπης. Η Κωνσταντινούπολη είναι επίσης μία από τις μεγαλύτερες μητροπολιτικές πόλεις της Ευρώπης, καλύπτοντας περίπου 5.500 χιλιόμετρα2, με πληθυσμό άνω των 15 εκατομμυρίων το 2017, που αντιστοιχεί στο 18% της χώρας . Η σημαντική αύξηση του πληθυσμού που σημειώθηκε λόγω της βιομηχανικής ανάπτυξης και της απρογραμμάτιστης αστικοποίησης, κατά το δεύτερο μισό του εικοστού αιώνα, είχε ως αποτέλεσμα κρίσιμους μετασχηματισμούς της δομής και της μορφολογίας της πόλης.

Μεθοδολογία

Προεπεξεργασία δορυφορικών εικόνων και εξαγωγή φασματικών υπογραφών

Οι δορυφορικές εικόνες που χρησιμοποιήθηκαν σε αυτή την έρευνα αποκτήθηκαν σε συνθήκες καθαρού ουρανού, με ελάχιστες ατμοσφαιρικές διαταραχές. Για κάθε περιοχή χρησιμοποιήθηκε μία μόνο εικόνα. Ως εκ τούτου, δεν ήταν αναγκαίο ένα βήμα προεπεξεργασίας της ατμοσφαιρικής διόρθωσης. Η χωρική ανάλυση των ζωνών εικόνας Sentinel-2 διέφερε μέσω των τμημάτων μήκους κύματος. Έτσι, υπήρχε ανάγκη για ομοιόμορφη χωρική ανάλυση για αναλύσεις όπως η δημιουργία φασματικού προφίλ βάσει σημείου, η δημιουργία φασματικού δείκτη και η ταξινόμηση πολυφασματικής εικόνας. Oι ζώνες ανάλυσης 20 m και 60 m των εικόνων Sentinel-2A ανανεώθηκαν στα 10 m, χρησιμοποιώντας την πλησιέστερη μέθοδο γείτονα, για να διατηρηθεί η ακεραιότητα της χωρικής ανάλυσης. Στο επόμενο βήμα, εξήχθησαν φασματικά προφίλ διαφόρων τύπων αντικειμένων γης, συμπεριλαμβανομένων των δασών broadleaf, των φυλλοβόλας δασών, των γεωργικών εκτάσεων, της αστικής πράσινης κάλυψης, της δομημένης, βιομηχανικής περιοχής, της αραιής κατοικημένης περιοχής, του θαλασσινού νερού, των υδάτων της λίμνης, της άσφαλτου και της γυμνής γης, για να εξετάσουν και να συγκρίνουν την ικανότητα διαχωρισμού των ζωνών εικόνας Sentinel-2A(Σχήμα 2 ). Κατά τη διάρκεια της διαδικασίας ταξινόμησης, τα δείγματα από πλατύφυλλα και φυλλοβόλα δάση τοποθετήθηκαν στην κατηγορία των δασών, δείγματα από αστική πράσινη κάλυψη και γεωργικές εκτάσεις ανατέθηκαν στην κατηγορία βλάστησης και δείγματα από λίμνες και θάλασσες ανατέθηκαν στην κατηγορία ύδατος, για να ληφθούν οι επτά τύποι κάλυψης γης που χρησιμοποιήθηκαν στην παρούσα έρευνα. Όπως δείχνει το σχήμα 2, η γυμνή γη είχε παρόμοια αντανάκλαση με τις κατοικημένες περιοχές και ήταν δύσκολο να προσδιοριστούν αυτές οι δύο κατηγορίες με τη χρήση ενός ενιαίου δείκτη. Ήταν απλό να προσδιοριστούν τα υδάτινα σώματα από άλλους τύπους κάλυψης γης, λόγω της μοναδικής φασματικής υπογραφής τους. Η σταδιακή μείωση της αντανάκλασης από τη ζώνη 1 στη ζώνη 12 ήταν ειδική για τα υδατικά συστήματα. Μια σημαντική αύξηση της αντανάκλασης στις ζώνες των κόκκινων άκρων (B5, B6, B7) και nir ζώνες (B8, B8a), σε σύγκριση με την κόκκινη ζώνη (B4) ήταν ειδική για την κάλυψη βλάστησης και θα μπορούσε να χρησιμοποιηθεί για την ανίχνευση φυτικών περιοχών. Επιπλέον, οι αναλύσεις της καμπύλης αντανάκλασης απέδειξαν ότι τα B1, B9 και B10 (εγγενής ανάλυση 60 m) δεν μπορούσαν να χρησιμοποιηθούν για τον διαχωρισμό των κατηγοριών κάλυψης γης. Οι παρατηρήσεις αυτές θα μπορούσαν να εξηγηθούν από τα χαρακτηριστικά αυτών των ζωνών. Β1 (παράκτιο αεροζόλ), επηρεασμένο έντονα από την ατμόσφαιρα και από τα Β9 και Β10, τα οποία ήταν υδρατμοί και κίρρος, δεν παρείχαν φασματικές πληροφορίες σχετικά με την επιφάνεια της Γης. Έτσι, αυτές οι ζώνες αφαιρέθηκαν από τα δεδομένα και πραγματοποιήθηκαν περαιτέρω αναλύσεις με τις υπόλοιπες 10 ζώνες. Η φασματική αξιολόγηση έδειξε ότι η κύρια πρόκληση ήταν ο διαχωρισμός της γυμνής γης και των κατοικημένων περιοχών, που ήταν ο κύριος στόχος αυτής της έρευνας.

Δημιουργία εικόνων πολλαπλών δεικτών

Για το στοιχείο κάλυψης βλάστησης, εξετάστηκε ο δείκτης ομαλοποιημένης βλάστησης με βάση τις κόκκινες άκρες (NDVIre) και δύο γνωστοί δείκτες βλάστησης, ο προσαρμοσμένος στο έδαφος δείκτης βλάστησης (SAVI) , και ο NDVI. O Hansen et al. , είχε αξιολογήσει για πρώτη φορά το NDVIre αναλύοντας δεδομένα υπερφασματικής αντανάκλασης. Η Delgado et al.Mκαι η Frampton et al είχαν δοκιμάσει για πρώτη φορά τη δυνατότητα εφαρμογής της στα σχεδιασμένα τμήματα μήκους κύματος Sentinel-2, πριν από την έναρξη της δορυφορικής αποστολής, χρησιμοποιώντας δεδομένα από διάφορες επιτόπιες εκστρατείες του ΗΣΥ σε γεωργικές εγκαταστάσεις. Τα αποτελέσματα και των δύο μελετών έδειξαν ότι η εφαρμογή αυτού του δείκτη στο κόκκινο B4 του Sentinel-2 (665 nm) και στις νέες ζώνες B5 (705 nm) του κόκκινου άκρου, παρείχε υψηλές συσχετίσεις κατά την εκτίμηση του δείκτη περιοχής φύλλων και της περιεκτικότητας σε χλωροφύλλη. Pu et al και Zhu et al. εισήγαγε τη λειτουργική χρήση αυτού του δείκτη για τις δορυφορικές εικόνες Worldview-2. Απ' όσο γνωρίζουμε, αυτή η έρευνα είναι μια νέα αξιολόγηση του NDVIre σχετικά με τις επιχειρησιακές εικόνες Sentinel-2A. Τέλος, για το συστατικό του υδατικού σώματος, εξετάστηκε ο ομαλοποιημένος δείκτης νερού διαφοράς (NDWI) και ο τροποποιημένος ομαλοποιημένος δείκτης νερού διαφοράς (MNDWI). Οι τύποι που σχετίζονται με τους παραπάνω φασματολογικούς δείκτες ήταν:

NDTI=((SWIR 1−SWIR 2))/((SWIR 1+SWIR 2)) (1)

NDVIre=((RedEdge 1−Red))/((RedEdge 1+Red)) (2)

NDVI=((NIR−Red))/((NIR+Red)) (3)

SAVI=((NIR−Red))/((NIR+Red+0.5))×1.5 (4)

NDWI=((Green−NIR))/((Green+NIR)) (5)

MNDWI=((Green−SWIR 1))/((Green+SWIR 1)) (6)

Πειραματική σύγκριση των υφιστάμενων ενσωματωμένων δεικτών

Το κύριο ενδιαφέρον για την ταξινόμηση της LCU στις αστικές περιοχές είναι ο διαχωρισμός των γυμνών εκτάσεων και των κατοικημένων περιοχών, λόγω των παρόμοιων φασματικών χαρακτηριστικών τους. Η εξόρυξη της γυμνής γης είναι ένα δύσκολο έργο, λόγω της πολυπλοκότητας των συστατικών του εδάφους και των φασμάτων του εδάφους. Όπως δήλωσε ο Ben-Dor et al το χημικό συστατικό επηρεάζει άμεσα τη φασματική υπογραφή γυμνών εδαφών, οι οποίες μπορεί να είναι ισχυρές ή αδύναμες. Επιπλέον, πολλές από αυτές τις φασματικές υπογραφές αλληλοεπικαλύπτονται, γεγονός που καθιστά δύσκολο τον προσδιορισμό της κάλυψης του εδάφους. Κατά συνέπεια, τα φασματικά χαρακτηριστικά της κάλυψης του εδάφους, με διαφορετικά συστατικά και περιεκτικότητα σε νερό, μπορεί να διαφέρουν μεταξύ διαφορετικών περιβαλλόντων και εποχών, γεγονός που καθιστά δύσκολη τη διαφοροποίηση της γυμνής γης.

Ομαλοποιημένος δείκτης διαφοράς άροσης (NDTI) Για να προταθεί ένας δείκτης που μπορεί να αναδείξει τις κατοικημένες περιοχές και να τις διαχωρίσει από τη γυμνή γη, τα φασματικά προφίλ αυτών των δύο χαρακτηριστικών γης αναλύθηκαν σε διάφορα σημεία του δείγματος. Έδειξε ότι η διαφορά αντανάκλασης μεταξύ των ζωνών SWIR (ζώνες 11 και 12) ήταν υψηλότερη για τα εικονοστοιχεία που επιλέχθηκαν από τη γυμνή γη από ό, τι για τα εικονοστοιχεία που επιλέχθηκαν από τις κατοικημένες περιοχές. Αυτό υποδείκνυε την πιθανή αποτελεσματικότητα αυτών των δύο ζωνών SWIR για τη διαφοροποίηση της δομημένης περιοχής από τη γυμνή γη. Διερευνήθηκε η δυνατότητα εφαρμογής του NDTI στις ζώνες SWIR των εικόνων Sentinel-2A για την κατοικημένη περιοχή και την εξόρυξη γυμνών εκτάσεων. Ο δείκτης αυτός προτάθηκε για πρώτη φορά από τους van Deventer et al. για πρακτικές εδάφους, διαχείριση άροσης και χαρτογράφηση καταλοίπων καλλιεργειών και εφαρμόστηκε με επιτυχία από τους Daughtry et al. και Eskandari et al. για γεωργικές πρακτικές και διαχείριση του εδάφους. Απ' όσο γνωρίζουμε, είναι η πρώτη φορά που η NDTI χρησιμοποιείται ως συστατικό στοιχείο για τη διάκριση και τον διαχωρισμό των Τα στοιχεία της NDTI που παρέχονται στο σχήμα 4 και οι υφιστάμενοι ενσωματωμένοι δείκτες που παρέχονται στο σχήμα 3 δείχνουν ότι η NDTI μπορεί να αναδείξει τις αστικές περιοχές και αυξάνει την αντίθεση μεταξύ της γυμνής γης (κόκκινος κύκλος) και της δομημένης περιοχής (μπλε κύκλος). Αυτή η οπτική επιθεώρηση δείχνει την πιθανή αποτελεσματικότητα της NDTI, σε σύγκριση με τους υπάρχοντες ενσωματωμένους δείκτες.κατοικημένων περιοχών και της γυμνής γης.

Ομαλοποιημένος δείκτης διαφοράς άροσης (NDTI)

Για να προταθεί ένας δείκτης που μπορεί να αναδείξει τις κατοικημένες περιοχές και να τις διαχωρίσει από τη γυμνή γη, τα φασματικά προφίλ αυτών των δύο χαρακτηριστικών γης αναλύθηκαν σε διάφορα σημεία του δείγματος. Έδειξε ότι η διαφορά αντανάκλασης μεταξύ των ζωνών SWIR (ζώνες 11 και 12) ήταν υψηλότερη για τα εικονοστοιχεία που επιλέχθηκαν από τη γυμνή γη από ό, τι για τα εικονοστοιχεία που επιλέχθηκαν από τις κατοικημένες περιοχές. Αυτό υποδείκνυε την πιθανή αποτελεσματικότητα αυτών των δύο ζωνών SWIR για τη διαφοροποίηση της δομημένης περιοχής από τη γυμνή γη. Διερευνήθηκε η δυνατότητα εφαρμογής του NDTI στις ζώνες SWIR των εικόνων Sentinel-2A για την κατοικημένη περιοχή και την εξόρυξη γυμνών εκτάσεων. Ο δείκτης αυτός προτάθηκε για πρώτη φορά από τους van Deventer et al. για πρακτικές εδάφους, διαχείριση άροσης και χαρτογράφηση καταλοίπων καλλιεργειών και εφαρμόστηκε με επιτυχία από τους Daughtry et al. και Eskandari et al. για γεωργικές πρακτικές και διαχείριση του εδάφους. Απ' όσο γνωρίζουμε, είναι η πρώτη φορά που η NDTI χρησιμοποιείται ως συστατικό στοιχείο για τη διάκριση και τον διαχωρισμό των Τα στοιχεία της NDTI που παρέχονται στο σχήμα 4 και οι υφιστάμενοι ενσωματωμένοι δείκτες που παρέχονται στο σχήμα 3 δείχνουν ότι η NDTI μπορεί να αναδείξει τις αστικές περιοχές και αυξάνει την αντίθεση μεταξύ της γυμνής γης (κόκκινος κύκλος) και της δομημένης περιοχής (μπλε κύκλος). Αυτή η οπτική επιθεώρηση δείχνει την πιθανή αποτελεσματικότητα της NDTI, σε σύγκριση με τους υπάρχοντες ενσωματωμένους δείκτες.κατοικημένων περιοχών και της γυμνής γης.

Αξιολόγηση ταξινόμησης και ακρίβειας Για την αξιολόγηση της ακρίβειας των αποτελεσμάτων ταξινόμησης, οι συνολικές ακρίβειες με τις μετρήσεις του χρήστη, του παραγωγού και της συνολικής ακρίβειας και της στατικής Kappa προήλθαν από τον πίνακα σύγχυσης. Η αξιολόγηση της ακρίβειας των αποτελεσμάτων ταξινόμησης πραγματοποιήθηκε με στρωματοποιημένα τυχαία σημεία και η αρχική εικόνα Sentinel-2A και οι εικόνες © Google Earth ως δεδομένα αναφοράς. Σχεδιάστηκε μια τυχαία κατανομή σημείων, ανάλογα με το δυναμικό ετερογένειας και την κάλυψη των τάξεων. Η κατανομή του δείγματος εκπαίδευσης και ο αριθμός των βαθμών που χρησιμοποιούνται στην αξιολόγηση ακρίβειας παρέχονται για κάθε κατηγορία και περιφέρεια του πίνακα 1.