Γεωμετρική διόρθωση των εικόνων MODIS με τη χρήση του MODIS Swath Reprojection Tool. Τα πρωτογενή δεδομένα που θα πάρουμε από την ιστοσελίδα του MODIS δεν μπορούμε να τα χρησιμοποιήσουμε αν δεν υποστούν πρώτα μια γεωμετρική διόρθωση. Η διόρθωση αυτή γίνεται με τη χρήση του εργαλείου MODIS Swath Reprojection Tool και περιγράφεται παρακάτω:
Τα αρχεία που κατεβάζουμε από την επίσημη ιστοσελίδα είναι όλα της μορφής .hdf .Συνολικά χρειαζόμαστε δύο αρχεία, ένα που θα περιέχει την εικόνα και ένα που θα περιέχει τις απαραίτητες πληροφορίες για την διόρθωση.
Παραδείγματος χάριν αν έχουμε δεδομένα από τον δορυφόρο Terra με μέγεθος εδαφοψηφίδας 500m την 204 ημέρα του 2009, οι ονομασίες των αρχείων θα έχουν την εξής μορφή:
Αρχείο εικόνων: MOD02HKM.A2009204.0925.005.2009204230137.hdf
Αρχείο γεωαναφοράς: MOD03.A2009204.0925.005.2009204204216.hdf
Από την ονομασία του κάθε αρχείου, μπορούμε να εξάγουμε τις παρακάτω πληροφορίες:
MOD : Δεδομένα από τον δορυφόρο Terra
MYD : Δεδομένα από τον δορυφόρο Aqua
02: Αρχείο εικόνων
03: Αρχείο που περιλαμβάνει τις απαραίτητες πληροφορίες για την διόρθωση.
1KM(1 km) : Μέγεθος εδαφοψηφίδας ίσο με 1km
HKM(Half km): Μέγεθος εδαφοψηφίδας ίσο με 500m
QKM (quarter km): Μέγεθος εδαφοψηφίδας ίσο με 250m
Γραφικό περιβάλλον της εφαρμογής 
Πατώντας το κουμπί
open Input File εισάγουμε το αρχείο .hdf που περιέχει τις εικόνες που κατεβάσαμε.
Δεύτερο βήμα,
Το πρόγραμμα από μόνο του, επιλέγει όλα τα κανάλια που περιέχει το συγκεκριμένο αρχείο, εμείς χρησιμοποιώντας τα βελάκια που φαίνονται στην παραπάνω προεπισκόπηση μπορούμε να επιλέξουμε μόνο τα κανάλια που χρειαζόμαστε.
Παραδείγματος χάριν, σε αρχείο .hdf που περιέχει εικόνες με μέγεθος εδαφοψηφίδας ίσο με 500m, οι αντιστοιχίες των καναλιών με αυτά του Landsat 5 είναι:
α) eik_EV_250_Aggr500_RefSB_b0 ->RED
β) eik_EV_500_RefSB_b1 ->GREEN
γ) eik_EV_500_RefSB_b0 ->BLUE
δ) eik_EV_250_Aggr500_RefSB_b1 ->INFRARED
Τρίτο βήμα,Στη συνέχεια πρέπει να εισάγουμε το αρχείο που δίνει τα απαραίτητα στοιχεία στο πρόγραμμα για να κάνει την γεωαναφορά. Αυτό το κάνουμε πατώντας το κουμπί
Open Geolocation file . Το αρχείο αυτό έχει στην ονομασία του τον κωδικό του δορυφόρου(MOD ή MYD) και συνοδεύεται πάντα από τον αριθμό 03.
Δηλαδή, MOD03 : το αρχείο γεωαναφοράς για δεδομένα από τον δορυφόρο Terra
και MYD03: το αρχείο γεωαναφοράς για δεδομένα από τον δορυφόρο Aqua.
Εισάγοντας το συγκεκριμένο αρχείο αυτόματα συμπληρώνονται και τα πεδία
UL corner και
ΒR corner. Αλλάζοντας τα φ και λ που εμφανίζονται σε αυτά τα πεδία και εισάγοντας τις γεωγραφικές συντεταγμένες των άκρων της περιοχή που μας ενδιαφέρει, μπορούμε ως τελικό προϊόν να πάρουμε απευθείας το συγκεκριμένο τμήμα της εικόνας.
*Τα φ και λ μπορούμε εύκολα να τα προσδιορίσουμε από το Google earth.Τέταρτο βήμα,Πατώντας στο κουμπί
Specify destination folder επιλέγουμε σε ποιο φάκελο θα αποθηκευτούν τα κανάλια της διορθωμένης εικόνας, παράλληλα πρέπει να εισάγουμε και την ονομασία που θέλουμε να δώσουμε π.χ. eik.tif
Κατά της συμπλήρωση του πεδίου η ονομασία που θέλουμε να δώσουμε στο αρχείο πρέπει να συνοδεύεται πάντα από την κατάληξη που υποδηλώνει σε τι format θα είναι το τελικό προϊόν. Για την συγκεκριμένη άσκηση επιλέγουμε το format GEOTIFF , οπότε συμπληρώνουμε την ονομασία του αρχείου μας με την κατάληξη .tif
**Παρατηρούμαι ότι στο πεδίο output file type επιλέγεται αυτόματα ο συγκεκριμένος τύπος αρχείου.Πέμπτο βήμα,Στο πεδίο
resampling type επιλέγουμε με ποια μέθοδο θέλουμε το πρόγραμμα να αναδομήσει το αρχείο μας. Συνήθως προτιμούμε το
Nearest Neighbor.
Έκτο βήμα,Στη συνέχεια ορίζουμε στο πρόγραμμα σε ποια προβολή θέλουμε την διορθωμένη εικόνα μας , από το κουμπί
Output Projection parameters.
Στην συγκεκριμένη άσκηση επιλέγουμε την
UTM ή την
Mercator και ανοίγοντας το
Edit Projection Parameters χρησιμοποιούμε ως Ε.Ε.Π. το
WGS84 και ορίζουμε την τιμή zone ίση με 34 (μόνο όταν εργαζόμαστε στην Ελλάδα)
**Εάν θέλουμε να δουλέψουμε στο ΕΓΣΑ 87, τότε επιλέγουμε την προβολή
Transverse Mercator και πατώντας το κουμπί
Edit Projection Parameters εισάγουμε τις παρακάτω παραμέτρους:
SMajor =0 SMinor=0 factor=0.9996
CentMer=24 OriginLat=0 FE=500000
fn=0 Ellipsoid= GRS1980Έβδομο βήμα,Στο πεδίο
output pixel size επιλέγουμε το μέγεθος της εδαφοψηφίδας που θα έχει η τελική μας εικόνα. Συνηθέστερα η τιμή αυτή θα ταυτίζεται με το μέγεθος pixel που θα έχει και το αρχείο .hdf που κατεβάσαμε.
Εφόσον ακολουθήσουμε τα παραπάνω βήματα, πατάμε το κουμπί
run και το πρόγραμμα σε ένα νέο παράθυρο ξεκινάει την γεωμετρική διόρθωση.
Όταν ολοκληρωθεί η διαδικασία, πηγαίνουμε στο φάκελο που ορίσαμε για την αποθήκευση των διορθωμένων δεδομένων. Εκεί θα δούμε για κάθε κανάλι που επιλέξαμε μια εικόνα της μορφής GEOTIFF.
Για να ενώσουμε τα κανάλια που δημιουργήθηκαν σε ένα ενιαίο αρχείο, θα χρησιμοποιήσουμε το πρόγραμμα ΕRmapper.
Η ιστοσελίδα του προγράμματος MODIS είναι:
http://modis.gsfc.nasa.gov/Ενώ για να κατεβάσουμε δεδομένα πάμε στην παρακάτω ιστοσελίδα:
http://rapidfire.sci.gsfc.nasa.gov/Το MODIS Swath reprojection tool το κατεβάζουμε από εδώ:
https://lpdaac.usgs.gov/lpdaac/tools/modis_reprojection_tool_swath
