Σύγκριση εκδόσεων

Εδώ βλέπετε τις διαφορές μεταξύ της επιλεγμένης έκδοσης και της τρέχουσας έκδοσης της σελίδας.

--- geo-processing [2010/05/24 19:03]
Sairin_Lote
+++ geo-processing [2020/11/26 17:40]
176.58.195.227 [Πειγραφή του Shapefile]
@@ Γραμμή 4: / Γραμμή 4: @@
 Η βιβλιοθήκη είναι γραμμένη σε C, αλλά έχει μεταφερθεί και σε Python.
-======Εγκατάσταση σε Ubuntu 10.04======
+=====Εγκατάσταση σε Ubuntu 10.04=====
@@ Γραμμή 15: / Γραμμή 16: @@
-=====GDAL module=====
+======GDAL module======
 Όπως αναφέρθηκε, το module διαχειρίζεται raster αρχεία, δηλαδή με απλά λόγια είναι χρήσιμο για την Ψηφιακή Τηλεπισκόπηση! Τα είδη των αρχείων που υποστηρίζονται περιλαμβάνουν όλα τα γνωστά αρχεία εικόνα (TIFF, Erdas, ERmapper κ.λπ.) και πολλά άλλα! Μια πλήρη λίστα με αυτά υπάρχει [[http://www.gdal.org/formats_list.html | εδώ!]]
@@ Γραμμή 25: / Γραμμή 26: @@
-====Άνοιγμα αρχείου====
+=====Άνοιγμα αρχείου=====
 Προχωράμε αμέσως στο άνοιγμα ενός αρχείου. Έστω ότι έχουμε μια εικόνα .ers (ERMapper format) στο '/home/user/Desktop' με ονομα 'image1.ers'. Αυτό γίνεται σε 3 στάδια:
@@ Γραμμή 46: / Γραμμή 47: @@
 Ορίζουμε στην μεταβλητή image την διαδικασία ανοίγματος της εικόνας. Η εικόνα ανοίγει χρησιμοποιώντας την μέθοδο Open() της κλάσσης gdal, η οποία δέχεται 2 arguments: το path της εικόνας και την κατάσταση ανοίγματος.
-====Ανάκτηση Βασικών Δεδομένων====
+=====Προσπέλαση Βασικών Δεδομένων=====
 Μπορούμε να πάρουμε κάποιες γενικές πληροφορίες για το αρχείο ως εξής:
@@ Γραμμή 75: / Γραμμή 77: @@
     print 'Strofi kata X kai kata Y:', image.GetGeoTransform()[2],',',image.GetGeoTransform()[4]
     print 'Diastaseis pixel kata X kai Y:', image.GetGeoTransform()[1],',',math.fabs(image.GetGeoTransform()[5])
+======OGR module======
+To OGR είναι το αντίστοιχο module διαχείρισης vector δεδομένων.
+=====Άνοιγμα αρχείου=====
+Το άνοιγμα αρχείων γίνεται σχεδόν όπως στην GDAL. Για παραδείγμα εδώ θα χρησιμοποιήσουμε ένα αρχείο Shapefile (.shp). Ο λόγος είναι ότι το shapefile είναι ένα από τα πιο διαδεδομένα αρχεία vector δεδομένων. Έστω λοιπόν ότι έχουμε ένα αρχείο test.shp. Αφήνω προς το παρόν την γεωμετρία του (point, line, polygon) φλου, επειδή θα αναφερθούμε σε εντολές που ισχύουν για όλα τα παραπάνω είδη. Έτσι έχουμε:
+    from osgeo import ogr
+    driver=ogr.GetDriverByName('ESRI Shapefile')
+    datasource=driver.Open('/home/user/Desktop/test.shp', 0)
+    layer=datasource.GetLayer()
+=> Line 1 : καλώ το ogr module
+=> Line 2 : στην περίπτωση του ogr πρέπει να καλέσω συγκεκριμένο driver. O driver που καλώ είναι ο 'ESRI Shapefile' (προσοχή στα κεφαλαία/μικρά). Τον καλώ μέσω της μεθόδου GetDriverByName και τον αποθηκεύω στην μεταβλητή driver.
+=> Line 3 : ανοίγω το αρχείο όπως και στην περίπτωση της GDAL. Η μέθοδος Open() δέχεται εκτός από το path του αρχείου και μια τιμή, η οποία έιναι 0 ή 1. Το 0 χρησιμοποιείται για read-only άνοιγμα, ενώ το 1 σε περίπτωση που θέλουμε να 'πειράξουμε' το αρχείο.
+=> Line 4 : Η βασική πληροφορία του Shapefile βρίσκεται στο Layer. Συνήθως οι βασικές εργασίες GIS γίνονται στο επίπεδο Layer ή στο επίπεδο Feature.
+=====Περιγραφή του Shapefile=====
+Το shapefile είναι ένα διανυσματικό αρχείο γεωγραφικής πληροφορίας. Το κάθε αρχείο απεικονίζει αποκλειστικά ένα είδος πληροφορίας (σημειακή, γραμμική ή επιφανειακή). Παρόλο που συχνά το shapefile παρουσιάζεται ως αρχείο, στην πραγματικότητα είναι μια ομάδα αρχείων. Από αυτήν μερικά αρχεία είναι απαραίτητα και κάποια άλλα είναι προεραιτικά.
+Τα απαραίτητα αρχεία για τον ορισμό του shapefile είναι τα .shp, .shx και .dbf. Τα δυο πρώτα σχετίζονται με την γεωγραφική πληροφορία του αρχείου ενώ το τελευταίο με την πληροφορία την οποία απεικονίζει. Για παράδειγμα, ένα σημειακό shapefile που απεικονίζει σημεία ενδιαφέροντος αποθηκεύει το είδος του αρχείου (σημειακό) και τις συν/νες των σημείων στα .shp και .shx. Αν τώρα για κάθε σημείο υπάρχει ένας μοναδικός κωδικός, ένας κωδικός για τον δήμο που ανήκει και μια περιγραφή σχετικά με το τι σημείο είναι (οργανισμός, μουσείο, κινηματογράφος κ.λπ.) αυτά αποθηκεύονται στο .dbf, το οποίο ουσιαστικά είναι ένα αρχείο-πίνακας με τόσες σειρές όσα τα σημεία και τόσες στήλες όσα τα περιγραφικά χαρακτηριστικά τους (σε αυτή την περίπτωση 3).
+Για την περίπτωση των γραμμικών και επιφανειακών στοιχείων στο .shp και .shx αποθηκεύονται όλες οι συν/νες της polyline που ορίζει την κάθε γραμμή.
+Από τα προεραιτικά αρχεία ιδιαίτερο ενδιαφέρον παρουσιάζει το αρχείο .prj (από το projection), το οποίο περιέχει την πληροφορία σχετικά με το σύστημα συντεταγμένων του αρχείου.
+=====Το ogrinfo=====
+Πριν προχωρήσουμε παρακάτω καλό είναι να αναφερθεί το πρόγραμμα ogrinfo, το οποίο δίνει βασικές πληροφορίες για ένα shapefile. Από ένα τερματικό (Applications -> Accessories -> Terminal) με την εντολή:
+    ogrinfo -so /path/to/the/shapefile.shp shapefile
+ Για παράδειγμα το αρχείο roads.shp με την εντολή
+   ogrinfo -so /home/user/roads.shp roads
+επιστρέφει:
+    INFO: Open of `/home/vagvaf/Data/greece.shp/roads.shp'
+       using driver `ESRI Shapefile' successful.
+    Layer name: roads
+    Geometry: Line String
+    Feature Count: 88925
+    Extent: (19.632112, 34.929228) - (27.353774, 41.754906)
+    Layer SRS WKT:
+    GEOGCS["GCS_WGS_1984",
+         DATUM["WGS_1984",
+              SPHEROID["WGS_1984",6378137,298.257223563]],
+         PRIMEM["Greenwich",0],
+         UNIT["Degree",0.017453292519943295]]
+    osm_id: Real (11.0)
+    name: String (48.0)
+    ref: String (16.0)
+    type: String (16.0)
+    oneway: Integer (1.0)
+    maxspeed: Integer (3.0)

TopoWiki

Εργαλεία Χρήστη

Εργαλεία ιστότοπου

Σύγκριση εκδόσεων

Εργαλεία ιστοσελίδας