1 Dálkový průzkum země v termální části spektra Princip termálního snímání Metoda pasivní nepřímá Fyzikální podstata Charakteristika • Termální oblast spektra zahrnuje vlnové délky od 3 µµµµm (atmosférická okna v intervalech 3-5 µµµµm a 8-12 µµµµm) • Intenzita elektromagnetického záření v těchto vlnových délkách je malá • V těchto vlnových délkách převažuje vlastní vyzařování objektů nad odraženým slunečním zářením • Dlouhovlnné vyzařování je vnějším projevem energie tělesa • Termální část spektra umožňuje získat především informace o teplotních vlastnostech objektů Charakteristika • Termální snímky lze vytvářet ve dne i v noci (denní a noční snímky) • Termální snímky se vyznačují menší prostorovou rozlišovací schopností • Z termálních snímků lze získat charakteristiky území, které jsou klasickými metodami zjistitelné pouze bodovým měřením. Vzhled termálních obrazových záznamů Termální snímek LANDSAT Thematic Mapper 6, 120 m pixel 2 Princip práce termálního skeneru • Princip mechanooptického snímání • Čidlo – fotonové detektory ochlazované na nízkou teplotu • Přístroje umožňují měřit radiační teplotu s přesností 0,1o C. Radiační (jasová teplota) • Tradičním měřením teploty (teploměrem je měřena tzv. kinetická teplota (termodynamická – meteorologie) - je interním projevem energie molekul • Stefan-Boltzmannův zákon: • Objekty vyzařují energii jako funkci své teploty • Vyzařování je vnějším projevem energie tělesa, která je v DPZ měřena • Změřená vyzářená energie vyjadřuje tzv. radiační (jasovou) teplotu tělesa 4 TM σ= Radiace z reálných povrchů • Reálný povrch emituje pouze část záření, které by emitovalo absolutně černé těleso o stejné teplotě • Tato míra vyzařovací schopnosti každého reálného tělesa se nazývá koeficient vyzařování - emisivita (εεεε) ( ) ( ) ε = M T M T R A MR(T) - intenzita vyzařování reálného objektu o teplotě T MA(T) - intenzita vyzařování absolutně černého tělesa o teplotě T Emisivita základních druhů povrchů • Většina reálných povrchů se chová jako tzv. selektivní zářič • Jejich emisivita nabývá hodnot 0 až 1 a mění se s vlnovou délkou • Určité povrchy vyzařují jako černé těleso v určitých vlnových délkách – (voda 6 až 14 µµµµm). • U takových povrchů není problém měřit teplotu metodami DPZ Emisivita základních druhů povrchů • Hodnoty emisivity základních druhů povrchů na souši se výrazně mění v čase i v prostoru • Silný vliv na měření tepelné energie metodou DPZ mají meteorologické podmínky • Pro přesná měření teploty je nutné provádět tzv. atmosférické korekce • V řadě případů lze z termálních snímků získat pouze kvalitativní informaci Interakce tepelného záření na zemském povrchu zákon zachování energie M M M MI A R T= + + MI - intenzita záření dopadajícího na povrch tělesa MA - intenzita pohlceného záření MR - intenzita odráženého záření MT - intenzita propouštěného (vedeného) záření I T I R I A I I M M M M M M M M ++= 3 Interakce tepelného záření na zemském povrchu ( )α λ = M M A I ( )ρ λ = M M R I Charakteristiky, které popisují schopnost povrchu absorbovat, odrážet a vést elektromagnetické záření: α(λ) - pohltivost ρ(λ) - odrazivost τ(λ) - propustnost ( )τ λ = M M T I ( ) ( ) ( )α λ ρ λ τ λ+ + = 1 Z Kirkhofova zákona plyne, že emisivita objektu je rovna jeho pohltivosti na dané vlnové délce. Položku propustnosti τ(λ) lze považovat v DPZ za rovnu nule Emisivita reálných povrchů • Čím nižší je odrazivost tělesa, tím vyšší je jeho emisivita a naopak. • Emisivita vodních ploch v termální části spektra je blízká 1 • Emisivita povrchů, které výrazně odrážejí záření je nízká (kovové povrchy) ( ) ( )ε λ ρ λ+ = 1 Emisivita reálných povrchů Stefan - Boltzmannův zákon pro reálná tělesa M T= εσ 4 T Trad kin= ε 1 4 Vztah mezi radiační a reálnou teplotou povrchu Radiační teplota povrchů je v důsledku emisivity nižší než teplota kinetická (reálná). • Kalibrace založená na současných pozemních měřeních teploty vzduchu • Kalibrace založená na snímání referenčních ploch o známé teplotě Kalibrace termálních obrazových záznamů Zdroje termálních snímků • Termální snímky poskytují všechny meteorologické družicové systémy (METEOSAT, NOAA) a to i v několika pásmech • Snímky jsou vytvářeny i na snímcích družic pro výzkum přírodních zdrojů (LANDSAT, TERRA) • Termální snímání se často provádí i z letadel (termovize) • Velká většina termálních obrazových záznamů poskytuje pouze kvalitativní informaci Termální snímky z družice NOAA, skener AVHRR, 1100 m pixel Pásmo 3B 3,6-3,9 µµµµm 4 Termální snímky z družice NOAA, skener AVHRR, 1100 m pixel Pásmo 4 10,3- 11,3 µµµµm Termální snímky z družice NOAA, skener AVHRR, 1100 m pixel Pásmo 5 11,5- 12,5 µµµµm Vzhled a vlastnosti termálních snímků • Termální snímky jsou většinou zobrazovány ve stupních šedi. • Obvykle světlé tóny na snímku reprezentují teplé povrchy. • Na meteorologických snímcích je škála tónů opačná. • Často se teplotní rozdíly vyjadřují v barevné stupnici (tzv. pseudobarevný režim) kanál 4 s barevně zvýrazněným teplotním intervalem 200-240 K Zdroj: ČHMÚ Vzhled a vlastnosti termálních snímků • Protože obrazy zachycují dlouhovlnnou část spektra mimo viditelné vlnové délky, mají často velmi neočekávaný vzhled daný teplotními rozdíly povrchů na snímku. • V důsledku rozdílných teplotních vlastností povrchů mají denní a noční snímky velmi odlišný vzhled. • Termální obrazové záznamy často nezobrazují ani hlavní terénní tvary • Mají málo orientačních bodů, často se proto používají v kombinaci s fotografií, snímkem z optické části spektra či mapou Vzhled a vlastnosti termálních snímků • Vzhled termálních obrazových záznamů je ovlivňován teplotními vlastnostmi povrchů (vodivost, tepelná kapacita, tepelná setrvačnost). • Vzhled snímků je ovlivňován dále též topografií, orientací snímaných povrchů. • Vegetační kryt vytváří velmi komplexní povrch s množstvím ke Slunci orientovaných nebo naopak zastíněných částí. • Voda výrazně mění teplotní vlastnosti povrchů. Termální snímky mohou sloužit k mapování půdní vlhkosti či obsahu vody ve vegetačním krytu. Vzhled a vlastnosti termálních snímků • Krátkovlnné sluneční záření bude ovlivňovat vzhled termálních snímků pouze v oblasti 3-5 µµµµm a nebude mít vliv na delší vlnové délky (8-14 µµµµm). • Denní snímky v první uvedené části spektra (3-5 µµµµm) budou obsahovat tzv. termální stíny v částech, kam nedopadá přímé sluneční záření. • Množství přijatého záření bude modifikovat také orientace svahů. • Tyto efekty jsou minimalizovány na obrazových záznamech pořízených těsně před východem Slunce. 5 Denní chod radiační teploty Termální snímky denní a noční Termální snímky denní a noční Termální snímky denní Terén bývá světlejší než voda, voda bývá velmi tmavá. Tmavšími tóny bývá znázorněn i lesní porost. Teplotní rozdíl mezi lesem a okolní volnou krajinou může být až 5o C. V denních hodinách je vedle vlastního dlouhovlnného vyzařování objektů registrováno i vyzařování okolních objektů nebo zpětné záření atmosféry, bývají denní snímky charakterizovány značnou tónovou proměnlivostí. Z tohoto důvodu se také často provádí termální snímání po západu Slunce, kdy jsou ještě patrné teplotní kontrasty různých povrchů. Termální snímky noční Na nočních snímcích bývají vodní plochy naopak teplejší než terén, který na rozdíl od vody teplo poměrně rychle ztrácí. Při teplotní inverzi, kdy se v údolích hromadí studený vzduch, bývají nižší polohy znázorněny tmavšími tóny než teplejší vyšší polohy. Teplotní mapování Při přesné kalibraci termálních obrazových záznamů jich lze využít k sestavování map teplotních charakteristik studovaného území. Cílem je nalézt vztah mezi hodnotou obrazového prvku a teplotou objektu: DN A B T= + ⋅ ⋅ε 4 DN - hodnota naměřená termálním radiometrem A, B - kalibrační konstanty radiometru ε - emisivita T - kinetická teplota 6 Teplotní charakteristiky povrchů • Jevy a procesy bodové (např. lokalizace termálních pramenů) • Jevy a procesy plošné (např. lesní požáry). Na termálních snímcích lze sledovat: Příklady aplikací • mapování povrchové teploty půdy • odhady radiační bilance území • studium teplotního znečištění • určování druhů hornin • lokalizace geologických zlomů • mapování půdních druhů a vodního obsahu půd • lokalizace netěsností zavlažovacích kanálů • zjišťování teplotních charakteristik vulkánů • studium evapotranspirace • lokalizace teplých a studených pramenů a gejzírů • studium cirkulace v nádržích, jezerech i mořích • lokalizace lesních požárů a požárů uhelných slojí Monitorování lesních požárů • Hodnocení rizika vzniku požárů • Detekce ohnisek a rozsahu požárů • Mapování následků požárů Monitorování požárů Peloponéský poloostrov Družice RESURS MSU-SK, Pixel 600 m Monitorování požárů Fire Atlas 7 Fire M3 – algoritmus detekce požárů ze snímků (NOAA-14 AVHRR) Kalibrace, radiometrická a geometrická korekce snímků (T3) > 315 K NO Yes Fire pixel Fire clear pixel Li et al., 1998 CCRS WWW 1. Vyhledávání „podezřelých“ pixelů pomocí 3. pásma AVHRR T3 - T4 > 14 K T4 ≥≥≥≥ 260 K R2 < 0.22 Eliminates warm background Eliminates clouds Land cover ≠≠≠≠ cropland or grassland Eliminates highly reflective clouds & surface T4 - T5 < 4.1 K and T3 - T4 ≥≥≥≥ 19 K One of neighbour pixels successful at first 6 tests TRUE FIRE PIXEL Eliminates cropland and grassland false fires Eliminates thin clouds with warm background Eliminates single pixel fire Li et al., 1998 CCRS WWW 2. Algoritmus detekce ploch postižených požárem Využitelnost SST: • Studium kolísání klimatu • Studium jevu ENSO • Vstup do GCM (Global Circulation Models) Mapování SST (Sea Surface Temperature) z družicových měření Výpočet SST ze snímků družice NOAA AVHRR a0,a1,a2 – konstanty odvozené empiricky z regresního vztahu T3 – radiační (jasová) teplota 3. pásma AVHRR (3,55 – 3,93 mikrometrů) T4 – radiační (jasová) teplota 4. pásma AVHRR (10,3 – 11,3 mikrometrů) SST = a0 + a1T3 + a2T4 Mapování extrémních projevů konvekce Teplota horní vrstvy bouřkové oblačnosti (Cb) 8 Mapování rozsahu lávových proudů Havaj, série nočních snímků z družice ASTER, plocha 9 x 12 km Mapování rozsahu lávových proudů Etna, termální snímek z družice ASTER, plocha 24 x 30 km Studium vulkánů Sopka Popocatepetl (Mexiko) Družice TERRA, skener ASTER Mapování oblaků sopečného popela Sopka Cleveland (Aleutské ostrovy) snímky z družice MODIS 9 Mapování teplotních poměrů New York, 16.9.2001 New York, 18.9.2001 Mapování tepelného ostrova města Variabilita povrchových teplot (LST) v oblasti města Brna a okolí vypočtená z termálního snímku pořízeného družicí Landsat 5 dne 15. června 2006 a histogram dokumentující rozdělení povrchových teplot na katastrálním území Brna Vztah mezi teplotou aktivních povrchů (LST) a charakterm zástavy Lokalizace poruch teplovodů Termovize – letecké snímání http://www.geog.ucsb.edu/~jeff/115a/remote_sensing/thermal/thermal irinfo.html http://www.realclimate.org/index.php?p=43 http://en.wikipedia.org/wiki/Urban_heat_island http://employees.oneonta.edu/baumanpr/geosat2/Urban_Heat_Island/ Urban_Heat_Island.htm http://adaptation.nrcan.gc.ca/perspective/health_3_e.php Mapování povrchové teploty oceánu 10 Teplotní poměry golfského proudu Družice TERRA, skener MODIS Teplotní poměry jsou mapovány na základě odhadu radiační teploty z termálních snímků v oblasti 11 a 12 mikrometrů.