[Fossgis-talk] Landsat/Aster

Sa Aug 14 16:16:51 CEST 2010

Hallo Markus, hallo Heiko
Tut mir leid dass ich mich jetzt erst in die Diskussion einschalte, aber diese Mails bekomme ich immer nur auf meinem Handy und da ist man schreibfaul.

Wie bereits erwaehnt sind Landsat Daten schon einige zeit frei vom USGS verfuegbar. Da ihr neues Landsat Product Generation System (LPGS) nun auch orthorektifiziert und die Daten als L1T (Level 1 Terrain corrected) generiert, werden die Daten meist nur noch in diesem Format via GLOVIS ausgeliefert. (Also kein weiterer Bedarf zur Orthorektifizierung) Der Standard resampling Algorithmus fuer das LPGS ist  MTF (Modulated Transfer Function resampling) welcher allen anderen Methoden (Nearest neighbour, Cubic Convolution, Biliniar interpolation) vorzuziehen ist. 

Trotz der Orthorektifizierung kann es sein dass Bilder nicht genau auf einnander liegen. Der automatische GCP Algorithmus in LPGS beruht auf einer chip Datenbank welche nicht global detailliert genug sein mag. Mein System (SARMES South African Resource management and Expert System frueher IRMES Infoterra Resource Management and Expert System welches leider momentan noch auf PCI basiert) nutzt hierzu einen anderen Ansatz der etwas besser ist. Also wenn du nun die L1T Bilder vom USGS weiter koregistrieren moechtest rate Ich Dir einfach ein polynomisches Modell, oder falls sehr genaue und genuegend GCP's vorhanden ein Spline with Tension Modell zu nutzen und die bilder dann auf 10m mit Nearest Neighbour zu resamplen. Landsat Bilder sind leider immer, schon beim L1G format resampled, da wir es mit der "wisk broom scanner" Architektur zu tun haben. Du kommst also ums resampling nicht herrum und die Kombination von MTF und NN nach 10m verspricht den geringsten Datenverlust.

Solltest Du dir nicht sicher sein mit  welchem Datenprodukt Du es zutun hast, schaue in die den Daten anliegende Datei *MTL.txt und dem Feld  PRODUCT_TYPE = "L1T". Es gibt noch RCC, L0Ra, L0Rp und L1G wobei ner normale Nutzer nur was mit den L1 Bildern anfangen kann. Zu weiteren Fragen in Bezug auf Landsat Daten empfehle ich http://landsathandbook.gsfc.nasa.gov/handbook.html

Was die Verfuegbarkeit der Bilder anbetrifft. 2007, 2008, 2009 war Suedafrika "cochair und chair" von CEOS. Dort haben wir die "Datademocracy Initiatieve" angeregt welches zur Folge hatte dass Brazielien und China sofort Bilder der CBERS Satellieten frei zur Verfuegung gestellt haben. Da sich die USA nicht lumpen lassen wollte, haben sie hierrauf auch ihr Archiev von Landsat und EO1 frei zur Verfuegung gestellt. Der Haken an der Sache ist jedoch nur dass das USGS immer noch Empfangs Gebuehren von anderen Empfangsstationen fuer Landsat verlangt, somit sind diese nicht besonders motiviert ihre Archieve kostenfrei zu oeffnen. Da der Speicher auf Landsat 5 schon viele Jahre ausser Funktion ist, hat das USGS nur wenige Landsat 5 Daten von anderen Laendern in Ihrem Archiev welche sie im Rahmen eines Projektes von den Empfangstationen erhalten haben. Ich versuche schon eine geraume Zeit lang das USGS zu einem Archievtausch fuer Bilder des Suedlichen Afrikas durchzufueren und uns dann auch die Empfangsgebueren zu ersparen. ein gewisses Interresse ist da, mal sehn wie weit wir kommen. Eine erhebliche Maenge an freien Daten niederiger Aufloesung gibt es auch ueber Geonetcast und ESA DDS. Institutionen mit einem kleinen Budget koennen sich fuer 1000 US $ so einen Empfang einrichten. 

Wenn Ihr mal auf unseren Katalog welcher jedoch noch in der Entwicklung steckt http://catalogue.sac.co.za/, koennt Ihr eventuell kostenfrei (ohne garantie einer Aufnahme) fuer kleine Studiengebiete unseren Sumbandilasat mit 6m Aufloesung programmieren. CBERS und SAC-C Daten sollten hierrueber auch bald frei zur Verfuegung stehen.

Markus, zu deinem Anliegen in Bezug auf ein Sensormodell. Diese liegen in den OSS Packeten, OTB (Orfeo Tool Box) und OSSIM fuer Zeilenscanner bei. Ich habe momentan eine Gruppe von Informatik Studenten der Uni Pretoria welche einen "plugin" fuer das automatische finden von GCP's und eventuelle orthorektifizieren von Bildern fuer QGIS implementieren. Abgabe ist im Oktober. Dir steht die Kode dann auch fuer GRASS zur verfuegung. Es gibt ein Abkommen zwischen Suedafrika, Brazilien und Indien um gemeinsahm FOSS Fernerkundungssoftware fuer die automatische radio- geometrische Korrektur und Klassifizierung von VHR Satellietenbildern zu entwickeln.

Naechste Woche geben Tim Sutton und Horst Duester einen FOSSGIS Kursus bei uns. Ich aergere mich schon eine geraume Zeit ueber das Begehren nach mehr ARCGIS lizenzen meiner Kollegen.

Mit freundlichen Gruessen

Wolfgang Lück
CSIR
Satellite Application Centre
Earth Observation Service Centre
Technology Manager
Tel +27 12 334 5142
Fax +27 12 334 5001
Email:wluck at csir.co.za

 ----- Original Message ------
 From:Markus Neteler
 Sent:Thursday, August 12, 2010 09:10
 To:Heiko Jacobs jacobs at cousin.de; 
 Cc: FOSSGIS-talk fossgis-talk-liste at fossgis.de; 
 Subject:Re: [Fossgis-talk] Landsat/Aster

Moin,

2010/8/12 Heiko Jacobs :

> Moin

>

> Um noch mal auf die alte Mail zurückzukommen, weil ich gerade

> langsam anfange, die vorhandenen Szenen langsam mal runterzuladen ...

>

> Zitat von Markus Neteler :

>>>

>>> Landsat gibt's ja kostenlos bspw. beim USGS Earth Explorer.

>>

>> Ja, und auch hier:

>> http://ftp.glcf.umd.edu/index.shtml

>

> Da scheinen aber eher ganz wenige rumzuliegen und schon etwas älter ...

> Also kaum was aus 2009/2010 in Ba-Wü (So 194-196 25-27)

> Oder kann ich nicht gucken?

Hier sind weitere:

http://glovis.usgs.gov/

>>> Eine TM 5-Szene habe ich schon mal testweise runtergeladen.

>>> Offenbar sind die dort kartenorientiert abgelegt, also schon gedreht.

>>

>> Ein Teil ist Norden-ist-oben, ein Teil Flugrichtung-ist-oben. Haengt

>> von der Karte ab...

>

> Bisher habe ich nur kartenorientierte gesehen, die auch schon

> angeblich auf UTM WGS84 projiziert sind mit Passpunkten ...

> Bloß wie genau ... Hmmm ...

Bei GLOVIS koennten auch andere sein (habe es nicht kontrolliert).

>>> Gibt's die kostenlos auch irgendwo am path orientiert, also ungedreht?

>>> Wäre mir vermutlich eigentlich lieber ...

>>

>> Du kannst mit "gdalwarp" von Flugrichtung-ist-oben nach Norden-ist-oben

>> rotieren.

>

> Extra umdrehen will ich ja nicht, denn:

>

>>> Jeder Bearbeitungsschritt mit resampling verfälscht die Daten ein

>>> klein wenig, daher würde ich das Drehen gerne mit einer präzisen

>>> Entzerrung zusammen legen ...

>>

>> Da hilft vermutlich nur im Archiv zu stoebern.

>> Ich glaube aber nicht, dass der schon gemachte Schritt dramatisch

>> ist bei orthorektifizierten Szenen. Bei nicht-flachem Gelaende wird

>> die Orthorektifizierung wesentlich wichtiger sein als eventuelle

>> Artefakte beim Rotieren.

>

> DGM hätte ich ja ...

> Und ich würde ja am liebsten alle geometrischen und radiometrischen

> Korrekturen in einem Schritt vom unverfälschten Original aus

> machen ...

Wenn Du einen freien Orthorektifikationsalgorithmus fuer Zeilenscanner

hast, waere ich interessiert (ggfls um ihn in GRASS einzubauen).

> Aber geschenkten Gäulen ... *seufz* ;-)

>

> MfG Heiko Jacobs

Gruesse

Markus

-- 

http://gis.ambiente.fmach.it/

_______________________________________________

Fossgis-talk-liste mailing list

Fossgis-talk-liste at fossgis.de

https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste

-------------- nächster Teil --------------
Ein Dateianhang mit HTML-Daten wurde abgetrennt...
URL: <https://lists.fossgis.de/pipermail/fossgis-talk-liste/attachments/20100814/2ab5f327/attachment.htm>