Preview only show first 10 pages with watermark. For full document please download

Folie 1

   EMBED


Share

Transcript

9 Remote Sensing 9.1 Physical Basics 9.2 Recording Techniques 9.3 Image Processing 9.4 Thematic Classification 9.5 Summary http://saturn.unibe.ch/.../Fotogrammetrie-Bildflug.pdf 704 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 9 Remote Sensing • A geographic information system (GIS) is a computer hardware and software system designed to – Collect – Manage – Analyze – Display geographically referenced data (geospatial; spatial) • It is a specialized information system consisting of a (spatial) database and a (special) database system Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 705 9 Remote Sensing • Recording on site – Terrestrial survey techniques • Global navigation satellite systems (e.g. GPS) • Very long baseline interferometer (VLBI) • Theodolite: measuring both horizontal and vertical angles optically • Total station: electronic theodolite (transit) integrated with an electronic distance meter www1.tudarmstadt.de http://de.wikipedia.org/ www.photolib.noaa.gov http://tu-dresden.de/ Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 706 9 Remote Sensing – Hydrographic survey • Sounding – Thematic survey • Map digitization • Survey by different sources – Statistics – Ministerial data – Technical literature http://tu-dresden.de/die_ tu_dresden/…/papers/fuhrland.pdf • Aerial survey and survey by remote sensing Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 707 9 Remote Sensing • Remote sensing is the acquisition of information of an object or phenomenon, by the use of device(s) that are not in physical or intimate contact with the object → indirect observation technique – That uses the electromagnetic radiation which is emitted by the observed object – That carries receiving devices on http://www.etsu.edu/cas/geosciences/ aircraft or spacecraft – That serves for the observation of the surface of the earth including all objects thereon, the oceans or the atmosphere Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 708 9 Remote Sensing • Photogrammetry – Greek: photo - grammetry ≈ image-measurement – Acquisition and analysis of images to determine the properties, form and position of arbitrary objects – Remote sensing is the acquisition of physical properties of objects whereas photogrammetry is the reconstruction of their geometric form based on this data http://www.gisdevelopment.net/…/mm063d_155.htm Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig www.maps.google.de 709 9 Remote Sensing • System Characteristics – Recording techniques • Radiometric resolution • Geometric resolution – Platform • • • • Kind of platform Altitude Orbit Period – Mission www.atmos.albany.edu/deas/ atmclasses/atm335/history.pdf www.irs.uni-stuttgart.de • Temporal coverage • Spatial coverage Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 710 9.1 Physical Basics • Electromagnetic waves as information carrier – Straight propagation with the speed of light – Speed of light = wavelength x frequency – Longer wavelength, lesser energy → more difficult to sense electrical field E λ: wavelength distance magnetic field M c http://www.fe-lexikon.info/images/ ElektromagnetischeWelle.jpg ν: frequency speed of light number of cycles that passes a certain point per second Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 711 9.1 Physical Basics • Electromagnetic spectrum – The electromagnetic spectrum is the range of all possible frequencies of electromagnetic radiation http://de.wikipedia.org/wiki/Bild:Elektro-magnetisches_Spektrum.JPG Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 712 9.1 Physical Basics • Behavior of electromagnetic waves at interfaces Transmission Scattering Reflection Absorption Emission Transmission + Reflection + Absorption = 1 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 713 9.1 Physical Basics sensor solar radiation sky radiation received signal atmospheric absorption and scattering reflection at the surface absorption and reflection in the water (suspended particles) scattered light scattering at the surface water depth reflection at the ground Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig [Al07] 714 9.1 Physical Basics – The albedo (lat. albedo = „whiteness“), reflectivity • The extent to which an object diffusely reflects light from the sun Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 715 9.1 Physical Basics – Albedo depends on wavelength • There is a strong difference between visual and infrared albedos of natural materials [Al07] Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 716 9.1 Physical Basics • The sun is the most important source of electromagnetic radiation • With the exception of objects at absolute zero, all objects emit electromagnetic radiation – The higher the temperature, the shorter the wavelength of maximum emission www.eduspace.esa.int/eduspac e/.../images/03.jpg 717 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 9.1 Physical Basics • Blackbody – Hypothetical source of energy that behaves in an idealized manner – It absorbs all incident radiation, none is reflected – It emits energy with perfect efficiency – Its effectiveness as a radiator of energy varies only as temperature varies http://mynasadata.larc.nasa.gov/images/BB_illustration2.jpg Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 718 9.1 Physical Basics • Emissivity – The ratio between the emitance of a given object and that of blackbody at the same temperature – Useful measure of the effectiveness of objects as radiators – Kirchhoff„s law: At thermal equilibrium, the emissivity of a body (or surface) equals its absorptivity surface emissivity (8-14 μm) blackbody 1 water, depending on pollution 0,973-0,979 water with oil film 0,96-0,979 snow 0,99 grass, 0,92-0,97 dense, short Sands, depending on water moisture Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 0,88-0,985 719 9.1 Physical Basics • Atmospheric window – Ultraviolet 0.01 - 0.4 μm • Reflected solar radiation • Because of atmospheric absorption it can only be used on aircrafts flying at low altitude • Main application: oil contamination detection in water (it is possible to identify the ship which has lost the oil!) www.geographie.ruhr-uni-bochum.de/agklima/vorlesung/strahlung/spektrum-atmosphaere.jpg Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 720 9.1 Physical Basics – Visible light 0.4 - 0.7 μm • Reflected solar radiation • Atmospheric influences particularly on blue and green light • Several applications, e.g. land use mapping – Near infrared 0.7 - 3 μm • Reflected solar radiation • Nearly no atmospheric influences http://altmed.creighton.edu/ • Main application: Classification of vegetation, forest health survey (healthy green plants strongly reflect near infrared radiation), classification of water (expanses of water seem dark as they absorb all) Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 721 9.1 Physical Basics – Far infrared (thermal energy) 3 - 1000 μm (usually : 8 - 14 μm) • Radiation emitted by the earth • Nearly no atmospheric influences (but clouds are impermeable, CO2 as well: greenhouse effect is measurable!) • Applicable day and night • Measurements beneath the surface to some extent (pipelines and leaks...) http://www.qualitas1998.net/paul/ • Applications for which the temperature and its change are important, e.g. sea temperature, thermal properties of stone, tectonics Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 722 9.1 Physical Basics • Emitted radiation • Nearly no atmospheric influences (capable to measure through clouds) • Measurements beneath the surface to some extent • Complex signal difficult to interpret • Low ground resolution (weak signal) • Disadvantageous SNR (Signal-to-Noise Ratio) → noisy images • Main applications: Meteorology (temperature profiles of the atmosphere) und oceanography (ice observation) Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig http://www.icefloe.net/hly0503/ – Passive microwaves 1 - 300 mm 723 9.1 Physical Basics – Active microwaves (radar) 1 - 300 mm http://www.wetteronline.de/ • Reflected, transmitted microwave radiation • Nearly no atmospheric influences (except reaction on water drops) • Applicable day and night • Measurements beneath the surface to some extent • Polarization effects • Higher ground resolution as passive microwaves • Complex signal • Doppler effect allows detection of moving objects (military applications), sea pollution Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 724 9.1 Physical Basics • Orbits – Altitude, orbital period, – Apogee/perigee • Greatest/least distance from the earth v orbital speed R Earth‘s radius= 6 370 km g0 gravitational acceleration on the Earth‘s surface = 9,81 m/s2 r radius of the satellite orbit – Inclination • Angular distance of the orbital plane from the equator www.satellitentracking.de/txt/ 04_satellitenbahnen.html Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 725 9.1 Physical Basics – Low Earth Orbit (LEO) • Heights between 200 and 600 km • Manned space stations: low inclination and heights above 400 km • Satellites with biological or material experiments and astronomical satellites • Spy satellites 90° inclination , perigee 200-250 km, apogee 600-900km – Medium Earth Orbits (MEO) • All orbits above 1000 km up to 36000 km • Navigation satellite systems (GPS, Glonass) • Small communication satellites http://www.tobedetermined.org/ Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 726 9.1 Physical Basics – Geosynchronous/geostationary Orbit (GSO) – Sun Synchronous Orbit (SSO) or Polar Earth Orbit (PEO) http://cimss.ssec.wisc.edu/sage/ • Orbit height approximately 35786 km, 0° inclination • Period is equal to the Earth's rotational period → It maintains the same position relative to the Earth's surface • Television satellites, weather satellites • Orbit height between 700 and 1000 km, inclination approximately 90° • Orbit ascends or descends over any given point of the Earth's surface at the same local mean solar time so the surface illumination angle will be nearly the same every time • Earth observation satellites Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 727 9.2 Recording Techniques • Passive systems: photography, scanner (optic, mechanical, optoelectronic) • Active systems: radar sensors active systems T/R reflected artificial radiation passive systems R R thermal radiation Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig reflected solar radiation 728 9.2 Photographic Systems • • • • Passive technique VIS and NIR (400-1000 nm) Analog storage medium Common types of films – Black and white/panchromatic: • Highest geometric resolution – Infrared • Unusual representation • Contrastier • Distinction between coniferous and deciduous forests • Surfaces of water easier to identify Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig [Al07] 729 9.2 Photographic Systems – Color/chromatic: • Worse geometric resolution as black and white, better thematic interpretability – Color infrared films: • The blue-sensitive layer is replaced by an emulsion sensitive to a portion of the near infrared region • Good thematic interpretability (vegetation). Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig [Al07] 730 9.2 Photographic Systems • Example: aerial photo of Braunschweig – Altitude approximately 1600 m – Ground resolution 10 cm – Color reversal film – Central projection – 21. April 2005 www.braunschweig.de/.../luftbilder.html Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 731 9.2 Photographic Systems • Example: Cosmos with KVR 1000 Camera – Russian spy satellite – Polar, sun-synchronous – Altitude 200km – Ground resolution 2m – Black and white film – Durability 45 days http://www.spotimage.fr/web/en/186-kvr-1000.php Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 732 9.2 Photographic Systems • Disadvantages – Difficult radiometric calibration – Low spectral bandwidth – Analog data • Advantages – Relatively cheap – High resolution – „Spontaneous“ recording of areas http://saturn.unibe.ch/.../Fotogrammetrie-Bildflug.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 733 9.2 Whisk Broom Scanner • Opto-mechanical scanner • A rotating 45 degree scan mirror continuously scans the Earth beneath the platform perpendicular to the direction of flight • The system collects data one pixel at a time sequentially • A scan line (mirror rotation) is equivalent to the image swath • The forward motion of the platform used to acquire a scene with sequential scan lines Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 734 9.2 Whisk Broom Scanner sensor platform aperture angle altitude instantaneous field of view IFOV: pixel flight direction scan direction a: geometric resolution > ground segment s: swath width http://www.uni-potsdam.de/.../febasis/febasis06_04-1206.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 735 9.2 Whisk Broom Scanner • Radiation imaging – Mirror rotates around an axis parallel to the flight direction – The radiation is split into its various wavelengths and focused onto detectors – Stored on magnetic tape (HDDT, CCT), remote data transmission beam splitter optical system telescope rotating mirror motor dispersion prism electronics amplifier, converter beam splitter interference grid radiation photodetectors streamer magnetic tape HDDT, CCT http://www.uni-potsdam.de/.../febasis/febasis06_04-1206.pdf 736 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 9.2 Whisk Broom Scanner • Advantages – Precise spectral and radiometric measurements – Wide total field of view – Digital data, remote data transmission • Disadvantages – – – – Relatively short dwell-time S-bend Panoramic distortion Low SNR → limited radiometric resolution http://landsat.gsfc.nasa.gov/images/archive/c0005.html Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 737 9.2 Whisk Broom Scanner • Landsat – American satellite series • • • • • • • Landsat 1: 1972-1978 Landsat 2: 1975-1981 Landsat 3: 1978-1983 Landsat 4: 1982-1993 Landsat 5: since 1984 Landsat 6: 1993 failure Landsat 7: since1999 4, 5 1-3 6, 7 http://de.wikipedia.org/wiki/Landsat Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 738 9.2 Whisk Broom Scanner – Orbit • Near polar, sun synchronous • Altitude: 907-913 km (Landsat 1-3), 705 km (Landsat 4-7 ) • Inclination: 99.2° (Landsat 1-3), 98.2° (Landsat 4-7) • Orbital period: approximately 100 minutes → 14 circulations per day • Provide complete coverage of the Earth every 18 (Landsat 1-3) respectively 16 days (Landsat 4-7) ground trace for Landsat1-3 for one day Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig [Al07] 739 9.2 Whisk Broom Scanner LANDSAT 4,5 (1-3) sensor Multispectral Scanner (MSS) pixel size 79 x 79 m² spectral 1 (4) 0,50 - 0,60 µm, channels green 2 (5) 0,60 - 0,70 µm, red 3 (6) 0,70 - 0,80 µm, near infrared 4 (7) 0,80 - 1,10 µm, near infrared thermal channel panchromatic channel LANDSAT 4,5 Thematic Mapper (TM) 30 x 30 m² 1 0,45 - 0,52 µm, bluegreen 2 0,52 - 0,60 µm, green 3 0,63 - 0,69 µm, red 4 0,76 - 0,90 µm, near infrared 5 1,55 - 1,73 µm, mid infrared 7 2,08 - 2,35 µm , mid infrared 6 10,4 - 12,5 µm (120 x 120 m²) LANDSAT 7 Enhanced Thematic Mapper Plus (ETM+) 30 x 30 m² 1 0,45 - 0,52 µm, bluegreen 2 0,52 - 0,60 µm, green 3 0,63 - 0,69 µm, red 4 0,76 - 0,90 µm, near infrared 5 1,55 - 1,73 µm, mid infrared 7 2,08 - 2,35 µm , mid infrared 6 10,4 - 12,5 µm (60 x 60 m²) 8 0,52 - 0,90 µm (15 x 15 m²) Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 740 9.2 Whisk Broom Scanner – Typical combination of channels green 0,5-0,6 μm false colour composite red 0,6-0,7 μm infrared 0,8-0,9 μm true colour composite http://www.uni-potsdam.de/.../febasis/febasis06_04-1206.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 741 9.2 Whisk Broom Scanner http://landsat.gsfc.nasa.gov/images/lg_jpg/f0012_77-89-06.jpg 742 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 9.2 Push Broom Scanner • Optoelectronical scanner • Employs a linear array of solid semiconductive elements to acquire one entire line of spectral data simultaneously • Scan lines perpendicular to the direction of flight • Forward motion of the platform to acquire a sequence of imaged lines to map a scene • CCDs (charge coupled device) to serialize parallel analog signals Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 743 9.2 Push Broom Scanner sensor platform : aperture angle altitude scan direction flight direction a: geometric resolution > ground segment s: swath width http://www.uni-potsdam.de/.../febasis/febasis06_04-1206.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 744 9.2 Push Broom Scanner • Radiation imaging – Tilted mirror, sometimes fixed sometimes tiltable – CCD image sensors in the image plane of the lens: line scan camera – Data storage in parallel memory chips, remote data transmission CCD sensors focal distance sample mirror optical system radiation lens aperture angle http://www.uni-potsdam.de/.../febasis/febasis06_04-1206.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 745 9.2 Push Broom Scanner • Spot (Systeme Probatoire d'Observation de la 1-3 Terre) – French satellite series • • • • • Spot-1: 1986-1990 Spot-2: since 1990 Spot-3: 1993-1997 Spot-4: since 1998 Spot-5: since 2002 4 http://www.uni-potsdam.de/... /febasis/febasis06_04-1206.pdf 5 – Two identical parallel sensors that can be operated independently of one another http://www.fe-lexikon.info/images/Spot5.jpg 746 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 9.2 Push Broom Scanner – Pivoting of the sensors can be employed for stereoscopy and also for a higher repeat circle – Sensors are operated from the ground stations angled view http://www.terraengine.com/Dgroundstation.cfm nadirlooking http://www.uni-potsdam.de/.../febasis/febasis06_04-1206.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 747 9.2 Push Broom Scanner – Orbit • • • • Sun synchronous Altitude: 822 km Inclination 98,7° Orbital period 101,4 min → approximately 14 circulations per day http://spot5.cnes.fr/.../35.htm geometric resolution SPOT 1-3 HRV (Instrument Haute Résolution Visible) 20 m (XS), 10 m (PN) SPOT 4 HRVIR (High Resolution Visible and Infrared) 20 m (XS), 10 m (P) radiometric resolution 0,5-0,9 μm: 3 VIS, 1 NIR 0,5-1,75 μm: 3 VIS, 1 NIR, 1 MIR sensor Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig SPOT 5 HRG (High Resolution Geometric) 10 m (VIS, NIR), 2,5/5 m (PAN), 20 m (MIR) 0,45-1,75 μm: 2 VIS, 2 NIR, 1 MIR 748 9.2 Push Broom Scanner Spot-1 HRV XS-Modus Spot-1 HRV P-Modus Detroit(USA), false colour composite, resolution 30 m San Diego(USA), panchromatic, resolution 20 m http://www.uni-potsdam.de/.../febasis/febasis06_04-1206.pdf 749 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 9.2 Push Broom Scanner Spot-5 HRG XS-Modus: stereo Dead sea (Jordan), panchromatic, 11/2002 resolution 2,5 m http://www.uni-potsdam.de/.../febasis/febasis06_04-1206.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 750 9.2 Radar • Radio Detection And Ranging • Principle: – Transmitting radar pulses (microwaves) and recording the reflected radiation → active – The transit time and the strength of the reflected signal is measured [LKC08] Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 751 9.2 Radar • Nadir: – The local vertical direction pointing in the direction of the force of gravity at that location • Range: – Line of sight • Azimuth: – Direction of flight http://ladamer.org/Feut/studium/fe1/FE1-06-Radar.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 752 9.2 Radar • Recording parameters – Polarization • Direction of the electric field which is perpendicular to the direction of propagation in the transmitted radar signal (H = horizontal,V = vertical) → 4 possibilities: HH,VV, HV,VH – Depression angle θd – Pulse length – Wavelength is divided into bands K-band X-band C-band L-band P-band 0,7-1 cm 2,4-4,5 cm 4,5-7,5 cm 15-30 cm 77-136 cm http://ladamer.org/.../FE1-06-Radar.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 753 9.2 Radar • Azimuth resolution AR depends on beam width(β) and the ground range distance (GR) β [LKC08] B GR1 R1 A B GR2 R2 A → Azimuth resolution is better in the near range Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 754 9.2 Radar • Ground range resolution (GRR) depends on the pulse length (τ) and the depression angle(θ) – Distinction between A and B only possible if the pulse passed A completely before reaching B [LKC08] → Better ground range resolution in the far range Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 755 9.2 Radar • In order to improve the resolution – Ground range • Decrease pulse length – Azimuth • Decrease wavelength • Increase antenna length • The azimuth resolution is unacceptably coarse for systems operating at satellite altitudes http://ladamer.org/Feut/studium/fe1/FE1-06-Radar.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 756 9.2 Radar • Synthetic aperture radar (SAR) – Scene is illuminated over an interval of time → history of reflections – The further an object the longer the time it is illuminated – As changes in frequency are systematic separate components of the reflected signal can be assigned to their correct position http://ladamer.org/Feut/studium/fe1/FE1-06-Radar.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 757 9.2 Radar • Doppler-effect • Physical antenna as small as possible • Azimuth resolution independent of GR and λ synthetic aperture – Approaching → increase in frequency – Receding → decrease in frequency radar pulse with frequency v2 object frequency v2 v1 – v2 > 0 v3 – v2 < 0 http://ladamer.org/Feut/studium/fe1/FE1-06-Radar.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 758 9.2 Radar • Comparison of the resolution between systems with real (a) and synthetic (b) aperture http://ladamer.org/Feut/studium/fe1/FE1-06-Radar.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 759 9.2 Radar • Interactions between radar signals and materials very complex as it depends on: – Wavelength – Incidence angle – Electrical properties – Moisture – Surface property http://www.meteo.physik.uni-muenchen.de/.../fe_boden_micro.html Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 760 9.2 Radar • Penetration depths of microwaves – Increases with decreasing wavelength – Decreases with increasing conductivity, which is also influenced by moisture vegetation – Is higher for smoother surfaces dry alluvium glacier [Al07] Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 761 9.2 Radar • Problem-oriented quantitative analysis of radar images is difficult as it relies mostly on hardly comprehensible interdependencies C-Band L-Band P-Band http://www.ccrs.nrcan.gc.ca/resource/tutor/gsarcd/pdf/bas_intro_e.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 762 9.2 Radar • ERS (European Remote Sensing Satellite) – ERS-1: 1991-1999 – ERS-2: since 1995 – Envisat: since 2002 – Orbit: • • • • • Sun synchronous 800 km altitude 98,5° inclination Orbital period 100 min Repeat circle 35 days http://www.esa.int/esaEO/ GGGWBR8RVDC_index_0.html http://www.raumfahrer.net/raumfahrt/envisat/ablauf.shtml 763 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 9.2 Radar – Instruments • SAR two modes of operation: image mode and wave mode in combination with the wind scatterometer (WS) • WS, active microwaves to measure ocean surface wind speed and direction • RA (Radar Altimeter); active: Ku-Band (13.8 GHz) measures variations in the satellite‟s height above sea level and ice with an accuracy of a few centimetres http://ceos.cnes.fr:8100/.../ers/earonc00.htm Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 764 9.2 Radar • GOME (Global Ozone Monitoring Experiment) Spectrometer (UV and VIS) provides information on ozone • ATSR (Along-Track Scanning Radiometer) an Imaging Infrared Radiometer (IRR: 4 channels, temperature) and a passive Microwave Sounder (MWS: 2 channels providing measurements of the total water content of the atmosphere within a 20 km footprint) • PRARE (Precise Range and Range Rate Equipment), all-weather microwave ranging system designed to provide measurements used for highly precise orbit determination and geodetic applications, such as movements of the Earth‟s crust • LRR (Laser-Retroreflector) passive optical device(IR) used as a target by ground-based laser ranging stations to determine the precise altitude Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 765 9.2 Radar • ATSR image of Crete http://earth.esa.int/earthimages/ Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 766 9.2 Radar • SAR image of Vorpommern – Three images acquired in September 1991 were overlaid each in one of the primary colors – Considerable changes of surface structure and moisture due to farming Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig http://earth.esa.int/earthimages/ 767 9.2 Radar • SAR image of the coast of Norway – In situations with little wind many different features appear on the ocean surface • Linear elements: current shear • Black areas: very light winds • Linear features and internal waves: currents alternated by the bottom topography, in shallow sea Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig http://earth.esa.int/earthimages/ 768 9.2 Recording Techniques • Comparison of the wavelengths used by different satellites http://www.fe-lexikon.info/images/sp_sat.gif Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 769 9.2 LIDAR • Light Detection and Ranging • Active sensor • Laser beams (UV, VIS near IR) to measure – Distance – Speed – Chemical composition and concentrations • Often imprecisely called "laser-radar" • LASER (Light Amplification by Stimulated Emission of Radiation) – Device that emits an intense narrow low-divergence beam of a specific wavelength [SX08] Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 770 9.2 LIDAR • Airborne Laserscanning – The distance between the sensor and the surface to be measured is determined from the runtime of a light pulse – By deflection of the laser beam and the forward movement of the aircraft a wide strip is scanned elliptical scanning fibre scanner Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig swiveling mirror 771 9.2 LIDAR – Parameters • • • • • • Sampling rate Scan angle Scan frequency Altitude Aircraft speed Distance between the trajectories – Recorded data • • • • Position Orientation of the aircraft Angle of every emitted beam Measured distance http://www.fht-stuttgart.de/fbv/fbvweb/veranstaltungen/GISDay/Rueckblick/gis_day2004_guelch.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 772 9.2 LIDAR – One laser beam might be reflected at different heights, e.g. in presence of vegetation: • Primary return: originate from the first objects a lidar pulse encounters, often the upper surface of a vegetation canopy • Well suited to create a digital object model (DOM) http://www.fht-stuttgart.de/.../gis_day2004_guelch.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 773 9.2 LIDAR • Secondary returns: lower vegetation layers and the ground surface • Last return provides data for a digital terrain model (DTM) if the vegetation is not too dense emitted signal pulse strength first echo last echo time signal strength scrup signal strength discrete echo determination full waveform digitisation terrain time time http://publik.tuwien.ac.at/files/PubDat_166922.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 774 9.2 LIDAR – Coordinates of the reflection points: • Calculated from the position and orientation of the sensor (by GPS and INS), the deflection angle of the beam and the distance between sensor and reflection point – Result: 3D point set along the trajectory http://www.photo.verm.tu-muenchen.de/.../EFE03_Kap23.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 775 9.2 LIDAR – Advantages • Uniform, dense acquisition of points • Acquisition of height information for DOM (with vegetation), as well as for DTM (without vegetation) • Accuracy in height between 50 and 15 cm in position1m • Fast area-wide acquisition • Active measuring method, nearly independent of illumination http://www.fht-stuttgart.de/.../gis_day2004_guelch.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 776 9.2 LIDAR – Disadvantages • Arbitrary points, no structure elements (prominent terrain points, borders) • Only single points, interpolation necessary • Relatively noisy http://www.fht-stuttgart.de/fbv/fbvweb/veranstaltungen/GIS-Day/Rueckblick/gis_day2004_guelch.pdf 777 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 9.3 Image Processing • Comparison between topographic maps and remotely sensed images Properties Remotely sensed image Topographic map Mapping not true to scale, Mapping true to scale, only minor image scales are only approximations, changes due to generalization additional errors if terrain is uneven Mapping not positional accurate, influenced by sensor alignment, grade, earth curvature, etc. Mapping positional accurate, only minor changes due to generalization No parallel projection Orthogonal parallel projection of the earth‘ s surface on the map reference plane Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 778 9.3 Image Processing Content Remotely sensed image Topographic map Communicating information in images Content defined causally by physical-chemical processes High information density, but irrelevant data included Unlimited diversity of forms Snap shot, contains transient data Information coded by graphic symbols content scale independent, no selection Up to date , short production time content scale dependent, reduction of information by generalization Not up to date, long production time, problem of revision Content defined conventionally, stipulated map symbols, explained in a legend Low information density, but all topographically relevant Limited number of map symbols Contains only topographically stable data Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 779 9.3 Image Processing Readability and interpretation Remotely sensed image Topographic map Varying image quality Uniform map quality No readability objects have to be interpreted Objects are directly readable as they are represented by clearly defined symbols Ambiguous, as interpretation depends on the interpreter Unambiguous independent of the user Real 3d impression possible, if third dimension by stereoscopy captured No real 3d impression, third dimension may only be coded by symbols Interpretation scale dependent, resolution determines if objects can be recognized Readability scale independent, granted by generalization Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 780 9.3 Image Processing http://homepage.univie.ac.at/thomas.engleder/lba_fe/lba_fe_18112004.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 781 9.3 Image Processing • Geometric errors, distortions – Inaccurate position and form of objects – Causes • Recording techniques and system • Relief • Platform (instability, motion) • Radiometric errors – Faulty pixel values – Causes • Atmospheric interference • Topographical effects • Technical defects (sensors, data transfer) Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig http://www.fas.org/irp/ 782 9.3 Geometric Errors • Goals of geometric corrections – Represent objects in uniform scale and true geometry (system correction) – Register overlapped images of a scene from different dates and views (image to image registration) – Register the image to real world map coordinates (image to map registration) • The planimetrically corrected image is called orthophoto aerial photo, uncorrected corrected → orthophoto [Al07] Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 783 9.3 Geometric Errors in Photographic Systems • Relief displacement – Points above the chosen reference plane are moved radially away from the center – Points below the chosen reference plane are moved radially towards the center – Radial displacement is reference plane larger near the border invisible space invisible space – Displacement diminishes at the center side view http://homepage.univie.ac.at/.../lba_fe_28102004.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 784 9.3 Geometric Errors in Photographic Systems [SX08] Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 785 9.3 Geometric Errors in Photographic Systems • Varying scale – Mapping scale changes with variations in terrain – The scale of objects closer to the camera is larger than that of objects being further away – The mapping of a rectangle that covers a terrace is not a rectangle Map: Aerial photo: constant scale varying scale higher terrace lower http://homepage.univie.ac.at/thomas.engleder/lba_fe/lba_fe_28102004.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 786 9.3 Geometric Errors in Scanners • Capturing a scene (image) takes a certain time • During the recording time the earth rotates eastward, so that the starting point of the last scan line is further west than that of the first line • The displacement depends on the relative speed of the satellite and the earth rotation and also on the size of the image • Example (Landsat 7): earth rotation → satellite – 33,8°S (Sidney) – Image size: 185 km → Offset: 10,82 km (ca 6%) motion↓ pixel http://ladamer.org/Feut/pdf/Kursbegleitung/dbv_vl/dbv_vl_kapitel3.pdf 787 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 9.3 Geometric Errors in Scanners • Whiskbroom scanner flight direction ↑ – The distance between sensor and terrain increases towards the edges – Size of scanning spots increases towards the edges [Al07] scan direction http://ladamer.org/Feut/pdf/Kursbegleitung/ dbv_vl/dbv_vl_kapitel3.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 788 9.3 Geometric Errors in Scanners – If the angular speed is constant, the image seems to be increasingly compressed towards the edges – More elevated surfaces are perpendicular moved away from the flight direction http://homepage.univie.ac.at/.../lba_fe_28102004.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig [Al07] 789 9.3 Geometric Errors in Radar Systems • Image geometry depends on the depression angle and the terrain • Oblique perspective (i.e. side-looking) leads to relief displacement – The type and degree of relief displacement in the radar image is a function of the angle at which the radar beam hits the ground, i.e. it depends upon the local slope of the ground http://www.ccrs.nrcan.gc.ca/.../bas_intro_e.pdf 790 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 9.3 Geometric Errors in Radar Systems • Foreshortening – Compression of those features in the scene which are tilted toward the radar – Foreshortening effects are reduced with increasing incident angles – Maximum when a steep slope is orthogonal to the radar beam http://www.ccrs.nrcan.gc.ca/.../bas_intro_e.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 791 9.3 Geometric Errors in Radar Systems • Radar shadow – Areas not illuminated by the radar – Caused by either concave or convex relief features if the slope on the opposite side of the antenna is larger than the depression angle – Typical in high relief terrain – Occur in the downrange direction – Most prominent with large incidence angle illumination http://www.ccrs.nrcan.gc.ca/resource/tutor/gsarcd/pdf/bas_intro_e.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 792 9.3 Geometric Errors in Radar Systems • Layover – Occurs when the reflected energy from the upper portion of a feature is received before the return from its lower – The top of the feature will be displaced, or “laid over” relative to its base http://www.ccrs.nrcan.gc.ca/.../bas_intro_e.pdf 793 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 9.3 Geometric Errors • Instability of the platform (aircraft) yawing rolling pitching change of altitude change of flight speed http://wdc.dlr.de/data_products/SURFACE/LCC/diplomarbeit_u_gessner_2005.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig [Al07] 794 9.3 Geometric Corrections • Model-based correction algorithm – Develop a model for a given recording techniques and platform that considers all its inherent causes for distortions – Parameterize the model to fit the actual conditions under which the image was taken – Suitable if the kind and cause of the distortion is known, as earth rotation, satellite orbit or positional parameters of the platform Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 795 9.3 Geometric Corrections • Mathematical function to map the positions of pixels on the coordinates of the same points in a map – Independent of the sensor platform, commonly used – Uses ground control points i.e. features visible on the image with known ground coordinates – Assign to each pixel a new position in the reference grid – Involves the following steps: e = f (c,r) r c n I. Choice of a suitable function (mapping) II. Coordinate transformation III. Resampling (interpolation) n = f (c,r) e corrected image Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig raw image 796 9.3 Image Processing • Radiometric corrections – Dark pixel subtraction • Assumption: the minimum value of every channel is 0 → for each channel the smallest measured value is subtracted from every value as it has to be an atmospheric influence, very simplifying – Radiance to reflectance conversion • Correction of values by known reflection values for certain surface properties – Atmospheric modeling • Develop a complex model for the transfer of EM energy under the atmospheric conditions (e.g. vapor content, ozone, temperature, etc.) to the time the image was taken – Determining missing pixels or rows by interpolation Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 797 9.3 Image Enhancement • Emphasizing structures – High pass filter 0 -1 0 -1 5 -1 0 -1 0 • Noise reduction (smoothing) – Low pass filter 19 19 19 19 19 19 19 19 19 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig [Al07] 798 9.3 Image Enhancement • Contrast enhancement – Alters each pixel value in the old image to produce a new set of values that exploits the full range of values – E.g. linear stretching • Chose a new minimum and maximum value • Intermediate values are scaled proportionally g‘(x,y) = g(x,y)⋅ c1 + c2 with c1 =255/[max(g(x,y)) – min(g(x,y))], c2 = -min(g(x,y)) 255 0 0 http://ivvgeo.uni-muenster.de/Vorlesung/FE_Script/3_2.html Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 255 799 9.4 Thematic Classification • Assignment of objects, features, or areas to classes based on their appearance on the imagery • Distinction between 3 levels of confidence – Detection: determination of the presence or absence of a feature – Recognition: object can be assigned an identity in a general class or category – Identification: object or feature can be assigned to a very specific class Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 800 8.4 Thematic Classification • Eight elements of image interpretation – Image tone • Lightness or darkness of a region within an image • Refers ultimately to the brightness of an area of ground as portrayed by the film • Influenced by vignetting, i.e. the image becomes darker near the edges [Ca07] Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 801 8.4 Thematic Classification – Image texture • Apparent roughness or smoothness of an image region • Caused by the pattern of highlighted and shadowed areas created when an irregular surface is illuminated from an oblique angle [Ca07] Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 802 9.4 Thematic Classification – Shadow • May reveal characteristics of its size or shape that would not be obvious from the overhead view alone • Important clue in the interpretation of individual objects [Ca07] Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 803 9.4 Thematic Classification – Pattern • Arrangement of individual objects into distinctive recurring forms • Usually follows from a functional relationship between the individual features that compose the pattern [Ca07] Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 804 9.4 Thematic Classification – Association • Specifies the occurrence of certain objects or features, without a strict spatial arrangement • Identification of a class implies that objects of another class are likely to be found nearby – Site • Refers to topographic position • E.g. sewage treatment facilities are positioned at low topographic sites near streams or rivers Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 805 9.4 Thematic Classification – Shape • Obvious clue to the identity of objects • Often shape alone might be sufficient to provide clear identification – Size • Relative size of an object in relation to other objects on the image provides the interpreter with an intuitive notion of its scale and resolution • Can be measured, permit derivation of quantitative information Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 806 9.4 Thematic Classification • Classification key – Provide a pictorial, exemplary representation of the examined areas or objects spruce pine douglas fir silver fir beech oak http://homepage.univie.ac.at/thomas.engleder/lba_fe/lba_fe_02122004.pdf 807 Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 9.4 Thematic Classification http://homepage.univie.ac.at/thomas.engleder/index_20072008.html Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 808 9.4 Thematic Classification • Multispectral classification – Ideally every class is defined by a typical multispectral signature, caused by a statistical distribution of the pixels of each class → Examination of the pixels of a multispectral image by mathematical algorithms • With regard to their homogeneity • Spatial distribution – Two types of classifiers • Unsupervised, autonomous • Supervised, interactive Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 809 9.4 Thematic Classification – After parameterization the multispectral feature space may be divided into • Primary feature spaces (reflectance, temperature etc.) • Linear transformed feature spaces http://homepage.univie.ac.at/thomas.engleder/lba_fe/lba_fe_25112004.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 810 9.4 Thematic Classification • Unsupervised classification – Assignment of pixels to spectral classes without prior knowledge of the existence or names of these classes – Cluster-algorithms to define spectral classes – Collateral information is used to define thematic classes a posteriori, e.g.: • Terrain surveys • Spectral measurements • Maps – Particularly suited to determine spectral properties of relevant thematic classes Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 811 9.4 Thematic Classification • Supervised classification – Analytical method to extract quantitative information – Assumption: every class in the feature space can be described by a probability distribution • Distribution assigns to every pixel the probability that it Band 5 belongs to the class in whose Band area it is located 7 • Usually Gaussian distribution • Number of variables = number of channels water control limits vegetation soil http://ladamer.org/Feut/pdf/Kursbegleitung/dbv_vl/dbv_vl_kapitel8.pdf Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 812 9.6 Summary • Physical basics – Electromagnetic radiation – Orbits • Recording techniques – Photographic systems (Aerial camera, Cosmos) – Whiskbroom scanner(Landsat) – Pushbroom scanner (SPOT) – Radar (ERS) – LIDAR (Airborne Laserscanning) Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 813 9.6 Summary • Image processing – Comparison between remotely sensed images and topographic maps – Causes of geometric errors – Image enhancement • Thematic classification – Visual interpretation – Quantitative image analysis Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 814 9.6 Summary collect GIS display analyse physics manage remote sensing objects classification recording techniques image enhancements / corrections Spatial Databases and GIS – Karl Neumann, Sarah Tauscher– Ifis – TU Braunschweig 815