Καλή ώρα σε όλες και όλους :-)
Οι χρήστες των παραθύρων του Βασιλάκη του Πύλη όταν τρέχουν προγράμματα python από την gdal θα πρέπει να ΜΗΝ βάζουν το .py στο τέλος παρότι έτσι δείχνει στο
παράδειγμα. Εμείς σε αυτή την άσκηση θέλουμε μία εντολή σαν αυτή:
gdal_calc -A ndvi_020_23july2015.tif -B ndvi_020_28april2013.tif --type=Int16 --outfile=ndvi_2015-2013_bin.tif --calc="(1.*A-1.*B)"
Εδώ έχουμε το πρόγραμμα <gdal_calc>, ως Α κατώφλι στον δείκτη βλάστησης της 1ης κομμένης και ως Β της δεύτερης, τύπο αποτελέσματος προσημασμένο ακέραιο <Int16>, αποτέλεσμα στην εικόνα ndvi_2015-2013_bin.tif και τύπο υπολογισμού την αφαίρεση (1.*A-1.*B).
Προσοχή εδώ στον πολλαπλασιασμό με το 1. (ένα τελεία) ώστε η πράξη να μετατρέψει σε πραγματικό το όρισμα Α και έτσι τα αποτελέσματα να είναι -1 0 ή +1.
Επίσης μην λησμονήσετε στην απεικόνιση στο QGIS να κάνετε τις ρυθμίσεις τις οποίες δείξαμε στο μάθημα ώστε να απεικονιστεί σωστά (
όσες και όσοι απουσιάζατε ζητείστε σημειώσεις από τους συμφοιτητές σας).
Όσες και όσοι είναι χρήστες του Μάκη τότε εκτελούν το πρόγραμμα κατά τα γνωστά (find κτλ) αλλά στα ορίσματα ισχύουν οι ίδιες ρυθμίσεις.
Για ότι χρειαστεί
ΕΔΩ είμαστε :-)
Λέφτερα,
Ch Iossif