Διαχωρισμός των δομημένου περιοχών από τη γυμνή γη σε μεσογειακές πόλεις χρησιμοποιώντας εικόνες 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)