Φωτοερμηνεία Τηλεπισκόπηση Εαρινό 2020 - Άσκηση 2η

[chiossif]

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

καλησπέρα σας :-)

Σήμερα 19/2 ολοκληρώθηκε η 2η Άσκηση του τρέχοντος εξαμήνου στην ΦΤ. Έτσι, για να βοηθήσουμε την διαδικασία επίλυσης αποριών, ανοίγουμε αυτό το θέμα στο οποίο μπορείτε να υποβάλετε τις ερωτήσεις σας και να διαβάζετε τις ερωτήσεις και τις απαντήσεις των άλλων.

Μην ξεχνάτε ΠΡΙΝ ρωτήσετε να διαβάζετε κατά σειρά: τις οδηγίες για το μάθημα και τις ασκήσεις του, τις σημειώσεις σας και ότι έχει ήδη ερωτηθεί εδώ.

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

Λέφτερα,
Ch Iossif

[chiossif]

Καλησπέρα σε όλες και όλους :-)

Γράφω εδώ τα βήματα για την εκτέλεση του κοψίματος της εικόνας της άσκησης (συνοπτικά) ώστε να έχουμε μια εικόνα των πιθανών προβληματικών σημείων.

0. Έχουμε κατεβάσει το αρχείο LC08_L1TP_165051_20170614_20170628_01_T1.tar.gz και το αποσυμπιέζουμε στον φάκελο:
C:\kapou ston disko\me kena kai ελληνικά\LC08_L1TP_165051_20170614_20170628_01_T1 όπου μέσα έχει αρχεία TIF κτλ (προηγούμενη άσκηση)

1. Φτιάχνουμε στον υπολογιστή μας (όχι στην επιφάνεια εργασίας ή αλλού) τον φάκελο C:\FT\LC08_L1TP_165051_20170614 (FT= φωτοερμηνεία τηλεπισκόπηση) χωρίς κενά ή ελληνικά. Εκεί μέσα αντιγράφουμε τα περιεχόμενα του LC08_L1TP_165051_20170614_20170628_01_T1

2. Επίσης φτιάχνουμε άδειο τον φάκελο C:\FT\ASK2 (ASK2 από το άσκηση 2)

3. Εκτελούμε το qgis και ανοίγουμε ως επίπεδο ράστερ ένα αρχείο TIF από τον φάκελο C:\FT\LC08_L1TP_165051_20170614 έστω το ..._B3.TIF
 
4. Σχηματίζουμε την εντολή gdalbuildvrt από το Build Virtual Ratser μέσα στο Miscellaneous υπομενού του Raster για την all1.vrt η οποία στο παράδειγμά μας πρέπει να είναι έτσι:

Κώδικας:
gdalbuildvrt
-separate
-srcnodata 0
C:\FT\ASK2\all1.vrt
C:\FT\LC08_L1TP_165051_20170614\LC08_L1TP_165051_20170614_20170628_01_T1_B2.TIF
C:\FT\LC08_L1TP_165051_20170614\LC08_L1TP_165051_20170614_20170628_01_T1_B3.TIF
C:\FT\LC08_L1TP_165051_20170614\LC08_L1TP_165051_20170614_20170628_01_T1_B4.TIF
C:\FT\LC08_L1TP_165051_20170614\LC08_L1TP_165051_20170614_20170628_01_T1_B5.TIF
C:\FT\LC08_L1TP_165051_20170614\LC08_L1TP_165051_20170614_20170628_01_T1_B6.TIF
C:\FT\LC08_L1TP_165051_20170614\LC08_L1TP_165051_20170614_20170628_01_T1_B7.TIF

(έχω βάλει τις παραμέτρους μία σε κάθε γραμμή)

5. Κόβουμε την εικόνα μας με το Raster-> Extraction-> Clipper για να φτιάξουμε την kommeni1.tif με την τελική εντολή να είναι στο παράδειγμά μας αυτή:

Κώδικας:
gdal_translate
-projwin 430000.000 1490000.000 490000.0000 1430000.000
-of GTiff
C:\FT\ASK2\all1.vrt
C:\FT\ASK2\kommeni1.tif

( οι αριθμοί εδώ ΔΕΝ είναι σωστοί αλλά κι εσείς πρέπει να έχετε διαστάσεις ~60km όπως εδώ. Το GTiff βγήκε διότι είχα την ρύθμιση GeoTIFF κάτω δεξιά στο παράθυρο του Clipper όταν έφτιαχνα output select την kommeni.tif)

ή

5. Κόβουμε την εικόνα μας με δεξί κλικ Export Save Raster Layer as... και ρυθμίζουμε τις διαστάσεις της τελικής εικόνας ~60χλμ και το όνομα του αρχείου εξόδου

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

Καλή μελέτη :-)

Λέφτερα,
Ch Iossif

[andronis]

Γεια σας και από μένα !
Οι χρήστες Windows λειτουργικού, ας προσέξουν στις εντολές στο Command Prompt (cmd) να χρησιμοποιούν slash αντί για backslash (το slash είναι αυτό: / )

[albi]

Χαιρετώ την κοινότητα. Θα ήθελα να κάνω μια έρωτηση σχετικά με την 2η άσκηση. Τις πληροφορίες του 1ου ερωτήματος που χρειάζονται λίγο παραπάνω ψάξιμο από που θα τις βρούμε?

[andronis]

Το 1ο ερώτημα της 2ης άσκησης, έχει εξηγηθεί αναλυτικά στο μάθημα. Γενικά οι απαντήσεις προκύπτουν από το λογισμικό QGIS, όταν θα έχετε φορτώσει την κομμένη εικόνα. Και προφανώς σε συνδυασμό με την εικόνα αναφοράς, από την οποία κόψατε την περιοχή ενδιαφέροντος.