topografoi.com



Author Topic: Ειδικά κεφάλαια Τηλεπισκόπησης - Άσκηση 4  (Read 13201 times)

an8rwpas

  • Posts: 98
  • Gender: Male
                                                     Γεωμετρική διόρθωση των εικόνων 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
« Last Edit: 09 Δεκ 2009, 17:37 by an8rwpas »

Holly

  • Posts: 12
  • Gender: Female
Αυτή είναι άσκηση στα ειδικά κεφάλαια τηλεπισκόπησης? ?????? Μάλλον απέδωσε η περσινή γκρίνια...  :P
« Last Edit: 09 Δεκ 2009, 22:28 by Holly »

an8rwpas

  • Posts: 98
  • Gender: Male
Το http://rapidfire.sci.gsfc.nasa.gov/  από όπου μπορούμε να κατεβάσουμε δεδομένα MODIS, δεν δουλεύει πια, αλλά μας στέλνει κατευθείαν στο http://lance.nasa.gov .

Η συγκεκριμένη ιστοσελίδα, επιτρέπει να κατεβάσουμε μόνο πολύ πρόσφατα δεδομένα  MODIS  (near real-time (NRT) data) και απευθύνεται κυρίως σε αυτούς που ενδιαφέρονται να συγκεντρώσουν δεδομένα για πολύ πρόσφατα γεγονότα.

Για ολόκληρο το αρχείο MODIS, πάμε στην ιστοσελίδα : http://ladsweb.nascom.nasa.gov/ .
« Last Edit: 01 Ιουλ 2011, 19:39 by an8rwpas »

chiossif

  • Posts: 334
Καλό μήνα :-) και ΚΑΛΟ ΚΑΛΟΚΑΙΡΙ :-) :-)  :-)   :-)    :-)      :-)        :-)            :-)

Πράγματι... κι εχθές που το χρειάστηκα το βρήκα αλλαγμένο και με τσάντισε...
Απομένει να το συνηθίσω...

Άλλωστε ως γνωστόν : "You cannot teach an old dog new tricks" και δεν θέλω να αρχίσω να γαυγίζω ;-)

υγ.
Θέλω να τονίσω εδώ την συνέπεια στο θέμα του μαθητή και φίλου an8rwpas -όνομα και πράμα- ...
« Last Edit: 02 Ιουλ 2011, 09:07 by chiossif »

an8rwpas

  • Posts: 98
  • Gender: Male

υγ.
Θέλω να τονίσω εδώ την συνέπεια στο θέμα του μαθητή και φίλου an8rwpas -όνομα και πράμα- ...


Ευχαριστώ δάσκαλε!

 

Copyright © topografoi.com