Ψ.Τ. 2018 - Άσκηση 8

[chiossif]

Κυρίες, Δεσποσύνες και Κύριοι,

καλή σας ώρα :-)

Με την δημοσίευση αυτή εγκαινιάζουμε μια νέα προσέγγιση για εμάς στη ΣΑΤΜ στον τρόπο επίλυσης αποριών: το φόρουμ. Ήδη υπάρχουν φοιτητικά φόρα αλλά εδώ θα είναι στοχευμένo στην νέα άσκηση 8 του μαθήματος και θα συμμετέχουμε κι εμείς οι διδάσκοντες. Όποια απορία έχετε, με έναν σχετικά ανώνυμο τρόπο, μπορείτε να την εκφράσετε εδώ πατώντας Απάντηση. Σας παρακαλώ ΠΡΙΝ ρωτήσετε διαβάστε τις οδηγίες και τις δημοσιεύσεις που έχουν προηγηθεί. Έτσι ίσως βρείτε απάντηση και ΔΕΝ χρειαστεί να ξαναρωτήσετε κάτι.

Οι οδηγίες εκτέλεσης της άσκησης είναι στα ακόλουθα βήματα:

με τον διαχειριστή αρχείων του υπολογιστή σας:

-> φτιάξτε έναν φάκελο στον υπολογιστή σας με το όνομα ask8 στον χώρο C;\. Το αποτέλεσμα πρέπει να είναι ακριβώς έτσι: C:\ask8
-> μέσα σε αυτό τον φάκελο βάλτε το αρχείο zip με τα προγράμματα και αποσυμπιέστε τα ώστε να τα βλέπετε άμεσα σε μορφή RasterAlgebra.py και όχι σε κάποιο υποφάκελλο
-> μέσα  σε αυτό τον φάκελο βάλτε την κομμένη (clipped) εικόνα σας σε μορφή tif (και όχι vrt ή άλλη) και μετονομάστε την σε kommeni.tif

τρέξτε τώρα το QGIS και κάντε τα ακόλουθα:

-> Επιλέξτε Plugins->Python console

-> στην διεπαφή >>> πληκτρολογείστε τις ακόλουθες εντολές:

--> import os
--> os.chdir('C:\ask8')

ώστε να μπορούμε να εκτελέσουμε τα προγράμματα και να διαβάσουμε την εικόνα μας. Αυτές οι εντολές είναι απαραίτητες κάθε φορά που ανοίγουμε το QGIS.

Για να εκτελέσουμε ένα πρόγραμμα πχ το RasterAlgebra.py θα πρέπει να γράψουμε στην διεπαφή >>> την εντολή:

---> execfile('RasterAlgebra.py')     # ( ΠΡΟΣΟΧΗ ! αν έχετε QGIS 3 η εντολή είναι exec(open('RasterAlgebra.py').read()) )

Ενδεικτικά οι εντολές είναι :

---> execfile('RasterAlgebra.py')
---> execfile('RasterFiltering.py')
---> execfile('RasterKMeans.py')
---> execfile('RasterKMeans_multiband.py')

Η τελευταία αργεί λίγο μερικά λεπτά οπότε υπομονή. Αν ο υπολογιστής γράψει ότι Δεν Αποκρίνεται αυτό είναι φυσιολογικό όσο τρέχει το πρόγραμμα.
Μετά την ολοκλήρωσή τους τα νέα αρχεία είναι έτοιμα και ...

... η συνέχεια εντός του Β1 ;-)

Καλή και γόνιμη μελέτη :-)

Λέφτερα,
Ch Iossif

υγ. Θέλω να ευχαριστήσω το topografoi.com και τον διαχειριστή του, εξέχων μέλος της κοινότητάς μας, για την φιλοξενία :-)

[chiossif]

Γεια και χαρά σε όλες και όλους :-)

Κι όμως εμφανίστηκε περίπτωση όπου η εικόνα εισόδου στον κώδικα των KMeans περιείχε τιμές inf και NaN (πρόκειται για ΕΙΔΙΚΕΣ ΜΗ ΑΡΙΘΜΗΤΙΚΕΣ τιμές ειδικές στην Πάιθον για το άπειρο και τον μη-αριθμό).

Το σφάλμα εμφανίζεται κάπως έτσι:
...
ValueError: array must not contain infs or NaNs

Γιαυτό χρειάζεται η προσθήκη της εντολής:

flatraster = numpy.nan_to_num(flatraster)

αμέσως μετά την αλλαγή σχήματος του πίνακα raster σε flatraster έτσι:

...
# flatten image to get line of values
flatraster=raster.reshape(rows*cols)
flatraster = numpy.nan_to_num(flatraster) # "Mask" NaN and inf numbers
flatraster[flatraster<0] = 0 # Mask NON vegetation
...

Όλα αυτά στο αρχείο RasterKMeans.py (και μην ξεχάσετε να το σώσετε αλλαγμένο έτσι; :-) )

Θυμίζω για κάθε πρόβλημα ΒΕΒΑΙΩΘΕΊΤΕ ότι ακολουθείτε ΠΙΣΤΆ τις οδηγίες. ΑΝ αλλάζει έστω και κάτι το οποίο σας φαίνεται ασήμαντο πιθανό να δημιουργήσει πρόβλημα...

Καλή δουλειά :-)

Λέφτερα,
Ch Iossif

[chiossif]

Γεια και χαρά σε όλες και όλους :-)

Όσοι έχετε χρησιμοποιήσει το snap για την 7η άσκηση ίσως ΔΕΝ μπορείτε να χρησιμοποιήσετε πίνακα σύγχυσης στο QGIS ή αλλού. Έτσι έφτιαξα ένα πολύ εύκολο προγραμματάκι για να το κάνετε. Απλά ονομάστε class1.tif την μία και class2.tif την άλλη και τρέξτε το κατά τα γνωστά. Το αποτέλεσμα το γράφει στην οθόνη και με αντιγραφή επικόλληση το πάτε όπου θέλετε :-)

Καλή δουλειά :-)

Λέφτερα,
Ch Iossif

[siachospan]

Καλησπέρα. Η εκπόνηση της άσκησης 8 πρέπει να γίνει και για τις δύο εικόνες της περιοχής μας ή μόνο για την πρώτη;

[Sotiris]

Καλησπέρα. Η εκπόνηση της άσκησης 8 πρέπει να γίνει και για τις δύο εικόνες της περιοχής μας ή μόνο για την πρώτη;

Αποφασίστηκε στην αίθουσα να γίνει για τη μία εικόνα