Δουλεύοντας με Δεδομένα Εδάφους

Τα εδαφικά ή υψομετρικά δεδομένα είναι χρήσιμα για πολλές GIS Αναλύσεις και χρησιμοποιούνται συχνά σε χάρτες.Το QGIS έχει καλές built-in δυνατότητες επεξεργασίας εδάφους. Σε αυτό το tutorial, θα δουλέψουμε μέσα από τα βήματα που δημιουργούν διάφορα προϊόντα από υψομετρικά δεδομένα όπως οι ισοϋψείς καμπύλες, σκίαση αναγλύφου κ.λπ.

Επισκόπηση εργασίας

Στόχος είναι να δημιουργήσετε χάρτη με ισοϋψείς καμπύλες και σκίαση αναγλύφου για την περιοχή γύρω από το όρος Έβερεστ.

Άλλες δεξιότητες που θα μάθετε

  • Αναζητώντας και κατεβάζοντας δωρεάν διαθέσιμα δεδομένα εδάφους.

  • Εξάγοντας ένα διανυσματικό επίπεδο ως KML και προβάλλοντας το στο Google Earth.

Πάρτε τα δεδομένα

Θα εργαστούμε με το σύνολο δεδομένων GMTED2010 από το USGS.Μπορείτε να κατεβάσετε αυτά τα δεδομένα από την ιστοσελίδα USGS EarthexplorerGMTED (Global Multi-resolution Terrain Elevation Data) είναι ένα σύνολο δεδομένων εδάφους για όλο τον κόσμο που είναι νεότερη έκδοση του συνόλου δεδομένων GTOPO30.

Εδώ είναι το πως θα αναζητήσετε και θα κατεβάσετε τα σχετικά δεδομένα από το USGS Earthexplorer.

  1. Πηγαίνετε στο USGS Earthexplorer . Στην καρτέλα Search Criteria, αναζητήστε για την περιοχή με όνομα Mt. Everest.Κάντε κλικ στο αποτέλεσμα για να επιλέξετε την περιοχή.

../_images/1144.png
  1. Στην καρτέλα Data Sets, ανοίξτε την ομάδα Digital Elevation και τσεκάρετε GMTED2010.

../_images/2100.png
  1. Τώρα μπορείτε να μεταβείτε στην καρτέλα : guilabel: Results και δείτε το μέρος του συνόλου δεδομένων που τέμνουν τα κριτήρια της αναζήτησής σας. Κάντε κλικ στο κουμπί Download Options.Τώρα θα πρέπει να συνδεθείτε στην ιστοσελίδα. Μπορείτε να δημιουργήσετε έναν δωρεάν λογαριασμό εάν δεν έχετε έναν.

../_images/350.png
  1. Επιλέξτε την επιλογή 30 ARC SEC και κάντε κλικ στο Select Download Option.

../_images/437.png

Θα έχετε τώρα ένα αρχείο με το όνομα GMTED2010N10E060_300.zip. Τα υψομετρικά δεδομένα διανέμονται σε διάφορες μορφές πλέγματος, όπως ASC, BIL, GeoTiff κλπ. Το QGIS υποστηρίζει μια μεγάλη ποικιλία μορφών πλεγμάτων <http://www.gdal.org/formats_list.html>` _ μέσω της βιβλιοθήκης GDAL. Τα δεδομένα GMTED λαμβάνονται ως αρχεία GeoTiff που περιέχονται σε αυτό το αρχείο zip.

Για ευκολία, μπορείτε να μεταφορτώσετε απευθείας ένα αντίγραφο και από τις δύο δέσμες στοιχείων από τις παρακάτω συνδέσεις.

GMTED2010N10E060_300.zip

Πηγή Δεδομένων : [GMTED2010]

Διαδικασία

  1. Ανοίξτε Layer ‣ Add Raster Layer και αναζητήστε το αρχείο zip που έχετε κατεβάσει.

../_images/538.png
  1. Υπάρχουν πολλά διαφορετικά αρχεία που δημιουργούνται από διαφορετικούς αλγορίθμους. Για αυτό το tutorial, θα χρησιμοποιήσουμε το αρχείο με το όνομα 10n060e_20101117_gmted_mea300.tif.

../_images/636.png
  1. Θα δείτε τα δεδομένα εδάφους που παρέχονται στο QGIS καμβά. Κάθε pixel στο πλέγμα εδάφους αντιπροσωπεύει το μέσο υψόμετρο σε μέτρα σε αυτή την τοποθεσία. Τα σκούρα pixels αντιπροσωπεύουν περιοχές με χαμηλό υψόμετρο και τα ανοιχτόχρωμα pixels αντιπροσωπεύουν περιοχές με μεγάλο υψόμετρο.

../_images/735.png
  1. Ας βρούμε την περιοχή που μας ενδιαφέρει. Από Wikipedia <http://en.wikipedia.org/wiki/Mount_Everest> _, βρίσκουμε ότι οι συντεταγμένες για την περιοχή ενδιαφέροντος - Mt. Everest - βρίσκεται στις συντεταγμένες 27.9881 ° N, 86.9253 ° Ε Σημειώστε ότι το QGIS χρησιμοποιεί τις συντεταγμένες σε μορφή (Χ, Υ), έτσι θα πρέπει να χρησιμοποιήσετε τις συντεταγμένες ως (γεωγραφικό μήκος και γεωγραφικό πλάτος). Κάντε επικόλληση το 86.9253,27.9881 στο κάτω μέρος του παραθύρου του QGIS όπου λέει : guilabel:` Coordinate` και πατήστε Enter. Το παράθυρο θα επικεντρωθεί σε αυτές τις συντεταγμένες. Για μεγέθυνση, Πληκτρολογήστε 1: 1000000 στο πεδίο : guilabel:` Scale` και πατήστε Enter. Θα δείτε το παράθυρο να έχει κάνει zoom στην περιοχή γύρω από τα Ιμαλάια.

../_images/834.png
  1. Τώρα θα κόψουμε το πλέγμα στην περιοχή που μας ενδιαφέρει. Επιλέξτε το εργαλείο Clipper από το Raster ‣ Extraction ‣ Clipper.

Note

Το μενού Raster στο QGIS προέρχεται από ένα κεντρικό πρόσθετο που ονομάζεται GdalTools.Εάν δεν μπορείτε να δείτε το μενού : guilabel: Raster, ενεργοποιήστε το πρόσθετο : guilabel:` GdalTools` από Plugins ‣ Manage and install plugins ‣ Installed. Δείτε το : doc: using_plugins για περισσότερες λεπτομέρειες.

../_images/932.png
  1. Στο παράθυρο Clipper,ονομάστε το εξαγόμενο αρχείο ως everest_gmted30.tif. Επιλέξτε το Clipping mode ως Extent.

../_images/1033.png
  1. Αφήστε το παράθυρο Clipper ανοιχτό και μεταβείτε στο κύριο παράθυρο του QGIS. Κρατήστε πατημένο το αριστερό πλήκτρο του ποντικιού σας και σχεδιάστε ένα ορθογώνιο που καλύπτει ολόκληρο τον καμβά.

../_images/1145.png
  1. Τώρα πίσω στο παράθυρο Clipper θα δείτε τις συντεταγμένες να έχουν συμπληρωθεί αυτόματα από την επιλογή σας.Σιγουρευτείτε ότι έχετε τσεκάρει την επιλογή Load into canvas when finished και κάντε κλικ στο OK.

../_images/1233.png
  1. Μόλις ολοκληρωθεί η διαδικασία,θα δείτε ένα νέο επίπεδο να έχει φορτώσει στο QGIS.Αυτό το επίπεδο καλύπτει μόνο την περιοχή γύρω από το όρος Έβερεστ.Τώρα είμαστε έτοιμοι να δημιουργήσουμε τις ισοϋψείς καμπύλες.Επιλέξτε το εργαλείο ισοϋψών καμπυλών από το Raster ‣ Extraction ‣ Contour.

../_images/1330.png
  1. Στο παράθυρο διαλόγου Contour ,επιλέξτε everest_gmted30 ως Input file.Ονομάστε το Output file for contour lines ως everest_countours.shp.Θα δημιουργήσουμε ισοϋψείς καμπύλες για διάστημα 100μ., γι’αυτό βάλτε 100 στο Interval between contour lines. Επίσης,τσεκάρετε την επιλογή Attribute name έτσι ώστε η τιμή του υψομέτρου να καταγραφεί ως χαρακτηριστικό κάθε ισοϋψούς καμπύλης.Κάντε κλικ στο OK.

../_images/1429.png
  1. Μόλις ολοκληρωθεί η διαδικασία,θα είναι ισοϋψείς καμπύλες να έχουν φορτωθεί στον καμβά.Κάθε καμπύλη στο επίπεδο αναπαριστά ένα συγκεκριμένο υψόμετρο.Όλα τα σημεία κατά μήκος μιας ισοϋψούς καμπύλης του πλέγματος που βρίσκεται από κάτω θα βρίσκεται στο ίδιο υψόμετρο.Όσο πιο κοντά είναι οι καμπύλες, τόσο πιο απότομη είναι η κλίση. Ας μελετήσουμε τις ισοϋψείς καμπύλες λίγο περισσότερο. Κάντε δεξί κλικ στο επίπεδο ισοϋψών καμπυλών και επιλέξτε Open Attribute Table.

../_images/1527.png
  1. Θα δείτε ότι κάθε χαρακτηριστικό της καμπύλης έχει μια ιδιότητα που ονομάζεται ELEV. Αυτό είναι το ύψος σε μέτρα που αντιπροσωπεύει κάθε καμπύλη.Κάντε κλικ μερικές φορές στο πάνω μέρος της στήλης για να ταξινομήσετε τις τιμές σε φθίνουσα σειρά.Εδώ θα βρείτε την καμπύλη που αντιπροσωπεύει το μεγαλύτερο υψόμετρο στα δεδομένα μας,δηλαδή το όρος Έβερεστ.

../_images/1626.png
  1. Επιλέξτε την πρώτη γραμμή και κάντε κλικ στο κουμπί Zoom to selection.

../_images/1724.png
  1. Μεταβείτε στο κύριο παράθυρο του QGIS .Θα δείτε ότι η επιλεγμένη ισοϋψής καμπύλη θα έχει κίτρινη επισήμανση.Αυτή είναι η περιοχή με το μεγαλύτερο υψόμετρο στο σύνολο των δεδομένων μας.

../_images/1824.png
  1. Τώρα, ας δημιουργήσουμε ένα χάρτη σκίασης αναγλύφου από το πλέγμα. Επιλέξτε Raster ‣ Analysis ‣ DEM (Terrain Models).

../_images/1922.png
  1. Στο παράθυρο διαλόγου DEM (Terrain Models), επιλέξτε everest_gmted30 ως Input file. Ονομάστε το Output file ως everest_hillshade.tif.Επιλέξτε το Hillshade ως Mode. Αφήστε όλες τις υπόλοιπες επιλογές όπως είναι.Σιγουρευτείτε ότι είναι τσεκαρισμένη η επιλογή Load into canvas when finished και κάντε κλικ στο OK.

../_images/2019.png
  1. Μόλις ολοκληρωθεί η διαδικασία,θα δείτε ακόμη ένα πλέγμα να έχει φορτώσει στον καμβά του QGIS .Από τη στιγμή που ίσως έχετε εστιάσει στην περιοχή του όρους Έβερεστ,κάντε δεξί-κλικ στο επίπεδο everest_hillshade και επιλέξτε Zoom to Layer Extent.

../_images/2123.png
  1. Τώρα θα δείτε την πλήρη έκταση του πλέγματος της σκίασης αναγλύφου.

../_images/2222.png
  1. Μπορείτε επίσης να δείτε επίπεδο ισοϋψών καμπυλών και να επιβεβαιώσετε την ανάλυσή σας, εξάγοντας το επίπεδο ισοϋψών καμπυλών ως KML και να το ανοίξετε στο Google Earth. Κάντε δεξί κλικ στο επίπεδο ισοϋψών καμπυλών, επιλέξτε: guilabel: Save as...

../_images/2318.png
  1. Επιλέξτε το Keyhole Markup Language [KML] ως το Format. Ονομάστε την έξοδο ως contours.kml και κάντε κλικ στο OK.

../_images/2416.png
  1. Αναζητήστε το εξαγόμενο αρχείο στο δίσκο σας και κάντε διπλό-κλικ πάνω του και ανοίξτε το με Google Earth.

../_images/2515.png
comments powered by Disqus

This work is licensed under a Creative Commons Attribution 4.0 International License