DIGITÁLNÍ ZPRACOVÁNÍ MATERIÁLŮ DPZ Úvodní a organizační informaceLukáš Herman Požadavky Účast na cvičeních Vypracování zadaných úkolů Prezentace semestrálního úkolu Zápočtový test – praktické úlohy Účast na cvičeních Povinná = bude kontrolována 1 absence s „rozumnou“ omluvou 2 a více absencí omluveny úředně (studijní oddělení, IS) Úkoly Protokoly (3) zadání postupně v průběhu semestru 2 týdny na vypracování možnost jedné opravy pozdní odevzdání – ztráta bodů Poster (ve skupinách 4 až 5 členech) v závěru semestru možnost získat bonusové body pravděpodobně týden na vypracování Bodování PROTOKOLY 10 b. za jeden, celkem 30 b. sankce za pozdní odevzdání – 1 bod za den ! ZÁVĚREČNÝ TEST 20 b. HODNOCENÍ – zápočet min. 25 b. F 0 - 24 b. E 25 - 29 b. D 30 - 34 b. C 35 - 39 b. B 40 - 44 b. A 45 - 50 b. A, B – pozitivně ovlivní výsledek zkoušky C – neutrální vliv na výsledek zkoušky D, E – negativně ovlivní výsledek zkoušky F – konečné F Prezentace Ve dvojicích (případně jednotlivě) připravit stručnou prezentaci zaměřenou na téma vybraného článku Pokyny: Prezentace délka 5 až 10 min Text 1 až 2 strany A4 Zaměřit se na zpracování obrazu (konkrétní metody, algoritmy, přednosti, omezení) Průběh výuky 1 17. a 20.9. úvod, pracovní prostředí sw. Geomatica, radiometrické a atm. chyby 2 24. a 27.9 statistika (histogramy) základních druhů povrchu, rozlišovací schopnosti snímku 3 1. a 4.10. geometrická transformace 4 8. a 11. 10 geometrická transformace 5 15. a 18. 10. úprava histogramu (zvýraznění, prahování, řezy, filtrace) 6 22. a 25. 10. spektrální zvýraznění, PCA 7 29. 10. a 1. 11. neřízená klasifikace 8 5. a 8. 11. neřízená klasifikace 9 12. a 15. 11 řízená klasifikace 10 19. a 22. 11 řízená klasifikace 11 26. a 29. 11. řízená klasifikace, pokročilé metody 12 3. a 6. 12. modelování, spektrální indexy, temporální analýzy 13 10. a 13. 12. opakování 14 17. a 20. 12. zápočtový test DIGITÁLNÍ ZPRACOVÁNÍ MATERIÁLŮ DPZ Pracovní prostředí programu Geomatica Opakování z DPZ Jaké znáte vlastnosti DPZ snímků? Jak se liší snímky z družic od letecké fotografie? Jaké znáte formáty DPZ snímků? Jaké znáte postupy (metody) zpracování snímku DPZ? Pojmy: odrazivost, spektrální chování, DN hodnota Geomatica Toolbar Moduly: Focus OrthoEngine Modeler EASI FLY! Dostupná trial verze (15 dní zdarma, plná funkčnost): http://www.pcigeomatics.com Geomatica - Focus Focus - zoom, pohyb po snímku Menu View Map View Mode Area View Mode Zoom Window Named Regions Focus - Menu Layer Add Layer Manager Histograms NumericValues Focus - Menu Layer Scatter Plot Spectral Plot Profile Table/Graph Zobrazení Barevná kompozice RGB x stupně šedi x nepravé barvy (PCT) Maps Properties - RGB Layer Properties Files View As Greyscale AS Pseudocolor Menu Layer - RGB Mapping Souřadnicové systémy S-JTSK S-42 WGS-84 (UTM) WGS-84 (souř. na elipsoidu) - 595812,719 -1153740,492 3585092,114 5294787,686 545839,132 5478846,320 14°28‘41,851“ E 49°56‘37,386“ N S-JTSK_Krovak_East_North Pulkovo_1942_GK_Zone_3 WGS_1984_UTM_Zone_33N GCS_WGS_1984 Označnív ArcGIS Struktura souboru PIX Soubor PIX Spektrální pásma Bitové masky Vektorové vrstvy Lookup Tables Pseudocolor Tables DIGITÁLNÍ ZPRACOVÁNÍ MATERIÁLŮ DPZ Radiometrické chyby, atmosférické korekce Radiometrické chyby Bitové chyby Chybějící řádek Páskování obrazu (stripping) Bitové chyby Analýza histogramu (odlehlé hodnoty) DN hodnoty výrazně odlišné od svého okolí Náhodné rozmístění, často nutné hledat manuálně Odstranění přepsáním chybné hodnoty v okně NumericValues Chybějící řádek Odstranění průměrováním DN hodnot v řádcích nad a pod chybějícím Nebo pomocí regresní závislosti Bodový graf z TM 3 proti TM 2 regrese - Kopie 2v*42c TM 3 = -53,5911+1,6625*x; 0,95 Int.spol. 50 55 60 65 70 75 80 85 90 TM 2 30 40 50 60 70 80 90 100 TM3 76 77 80 80 84 83 78 73 80 75 79 80 79 78 77 77 79 76 77 75 73 78 76 79 74 74 75 67 82 80 82 80 78 71 64 81 79 79 80 71 59 47 86 82 86 79 71 58 57 89 82 89 82 60 40 39 88 85 81 67 60 57 57 87 92 79 58 41 39 37 75 70 63 56 57 56 56 90 73 63 43 40 38 39 75 70 63 56 57 56 56 71 63 51 40 41 40 40 Lin. regrese pásem s vysokou korelací (zde TM2 a TM3) Šedé buňky jsou vypočtené hodnoty Páskování Snáze patrné na homogenních plochách (voda), při využití vhodného zvýraznění (Infrequency) K odstranění se využívá úpravy histogramu http://glovis.usgs.gov Atmosférické korekce Na snímky z jaké části spektra má největší vliv atmosféra? Metody korekce: Metoda nejtmavšího pixelu Regresní analýza Modelování atmosférických podmínek Metoda nejtmavšího pixelu Analýza histogramu TM4 (blízké infračervené) Vodní plochy by měly mít DN hodnotu 0 příspěvek atmosféry Odečíst (od všech hodnot) pomocí Raster Calculatoru Regresní analýza Průsečík přímky lineární regrese a osy x = velikost příspěvku atmosféry Ne vždy ji lze použít, je závislá na povaze vstupních dat Napovědět může hodnota R