GRASS GIS / Poznámky pro družicová data Landsat: Porovnání verzí

← Porovnání se starší verzí Porovnání s novější verzí →

Verze z 21. 9. 2010, 06:26

Tato stránka obsahuje poznámky pro stažení dat LandSat Thematic Mapper a následné zpracování těchto dat v prostřední GRASS GIS.

Časová řada

80. léta (co nejstarší)
přelom 80. a 90. let
kolem roku 2000
nejnovější

Postup pro manuální stažení dat

Data jsou dostupná na adrese: http://eros.usgs.gov (bližší informace zde)
- Data mohou být stažena také z adres: http://glovis.usgs.gov nebo http://edcsns17.cr.usgs.gov/EarthExplorer

Návod pro http://glovis.usgs.gov:

Import dat

Data je nejprve nutno dekomprimovat

tar xvzf LT51910252009213MOR00.tar.gz

Archiv obsahuje několik souborů ve formátu geotiff (jeden soubor na kanál), seznam identických bodů (GCP, ground control points) a metadata v textové podobě (MTL).

ls -w1
L5191025_02520090801_B10.TIF
L5191025_02520090801_B20.TIF
L5191025_02520090801_B30.TIF
L5191025_02520090801_B40.TIF
L5191025_02520090801_B50.TIF
L5191025_02520090801_B60.TIF
L5191025_02520090801_B70.TIF
L5191025_02520090801_GCP.txt
L5191025_02520090801_MTL.txt
README.GTF

Nástroj knihovny GDAL gdalinfo poslouží pro kontrolu metadat, např.

gdalinfo L5191025_02520090801_B10.TIF

Driver: GTiff/GeoTIFF
Files: L5191025_02520090801_B10.TIF
Size is 8441, 7441
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",15],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["unknown",1],
    AUTHORITY["EPSG","32633"]]
Origin = (398699.999999999941792,5679000.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Metadata:
  AREA_OR_POINT=Point
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  398700.000, 5679000.000) ( 13d32'54.25"E, 51d15'12.06"N)
Lower Left  (  398700.000, 5455770.000) ( 13d36'29.29"E, 49d14'46.73"N)
Upper Right (  651930.000, 5679000.000) ( 17d10'35.91"E, 51d14'31.60"N)
Lower Right (  651930.000, 5455770.000) ( 17d 5'13.62"E, 49d14'9.03"N)
Center      (  525315.000, 5567385.000) ( 15d21'18.49"E, 50d15'29.05"N)
Band 1 Block=8441x1 Type=Byte, ColorInterp=Gray

Vytvoření lokace

Hromadný import dat

Data lze naimportovat pomocí modulu r.in.gdal (File -> Import raster data -> Import raster data), ukázka z příkazové řadky (přepínač -e nastaví po importu region podle vstupní dat).

r.in.gdal -e input=L5191025_02520090801_B10.TIF output=b10 title="kanal 1"

Základní metadata vypíše r.info.

r.info map=b10

 +----------------------------------------------------------------------------+
 | Layer:    b10                            Date: Thu Sep 16 17:48:23 2010    |
 | Mapset:   PERMANENT                      Login of Creator: martin          |
 | Location: tm                                                               |
 | DataBase: /home/martin/grassdata                                           |
 | Title:    kanal 1 ( b10 )                                                  |
 | Timestamp: none                                                            |
 |----------------------------------------------------------------------------|
 |                                                                            |
 |   Type of Map:  raster               Number of Categories: 0               |
 |   Data Type:    CELL                                                       |
 |   Rows:         7441                                                       |
 |   Columns:      8441                                                       |
 |   Total Cells:  62809481                                                   |
 |        Projection: UTM (zone 33)                                           |
 |            N:    5679000    S:    5455770   Res:    30                     |
 |            E:     651930    W:     398700   Res:    30                     |
 |   Range of data:    min = 0  max = 255                                     |
 |                                                                            |
 |   Data Description:                                                        |
 |    generated by r.in.gdal                                                  |
 |                                                                            |
 |   Comments:                                                                |
 |    r.in.gdal -e input="L5191025_02520090801_B10.TIF" output="b10"          |
 |                                                                            |
 +----------------------------------------------------------------------------+

Data lze hromadně naimportovat pomocí jednoduchého skriptu. Následuje ukázka pro Bash a Python.

Bash

Ukázka skriptu pro Bash.

#!/bin/bash

idx=1
for file in *.TIF; do
    oname=`echo $file | cut -d'_' -f3 | cut -d'.' -f1`
    g.message "Importuji $file -> $oname..."
    r.in.gdal -e input=$file output=$oname title="kanal $idx"
    idx=$(($idx+1))
done

Python

Ukázka skriptu pro Python.

#!/usr/bin/python

import os
import glob

import grass.script as grass

idx = 1
for file in glob.glob1('.', '*.TIF'):
    oname = os.path.splitext(file)[0].split('_')[-1]
    grass.message('Importuji %s -> %s...' % (file, oname))
    grass.run_command('r.in.gdal',
                      flags = 'e',
                      input = file,
                      output = oname,
                      title = "kanal %d" % idx)
    idx += 1

wxGUI

Dále lze data naimportovat pomocí wxGUI.

File -> Import raster data -> Bulk import of raster data

Vizualizace dat

Příklad vizualizace kompozitního snímku (kanály 453).

Seznam měst

Brno                                                
Břeclav                                            
Česká Lípa                                       
Česká Třebová                                    
České Budějovice                               
Český Těšín                                    
Děčín                                               
Frýdek-Místek                         
Cheb                                              
Chomutov                                       
Chrudim                                         
Jablonec                                        
Jablonec nad Nisou

* Pantuckova
* Karochova
* Kuzelka
* Hron

Jihlava                                            
Jirkov                                             
Kadaň                                             
Karviná                                           
Klatovy                                    
Kutná Hora                                          
Liberec                                                
Litvínov                                            
Náchod                                    
Ostrava
Plzeň
Přerov
Příbram
Rožnov pod Radhoštěm
Sokolov
Teplice

* Jarosova
* Linhartova
* Zavadil
* Kopecky

Třebíč
Ústí nad Labem
Ústí nad Orlicí
Valašské Meziříčí
Vsetín
Zlín
Znojmo

* Kratochvil
* Sedlackova
* Svobodova
* Duchnova

Žďár nad Sázavou

Externí odkazy