Top Banner
Using Open Source Programming with Satellite Imagery Claire Porter Polar Geospatial Center University of Minnesota
29
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Porter

Using Open Source Programming with Satellite Imagery

Claire Porter

Polar Geospatial Center

University of Minnesota

Page 2: Porter

1. Electro-optical high-resolution commercial satellite imagery

2. PGC’s open source approach: GDAL and Python

3. Applications to other datasets

Page 3: Porter

high-res imagery

• High-resolution Commercial Satellite Imagery in ice-covered regions

– 0.5m – 4m resolution

– 1 to 8 electro-optical bands

– Passive sensors

– High temporal resolution

Page 4: Porter

high-res imagery

Page 5: Porter

Fuel Cache

Imagery ©2011 Digital Globe, Inc.

overland traverses

Presenter
Presentation Notes
Getting supplies to pole by ground is more time consuming than by air but much more cost effective.
Page 6: Porter

Imagery ©2012 Digital Globe, Inc.

field camp monitoring

Presenter
Presentation Notes
This camp is just 10s of miles up stream of the Pine Island Glacier Rift. If we can land a ski equipped LC-130 without a season to prepare a skiway NSF can save hundreds of thousands of dollars and get scientists working on site a year earlier.
Page 7: Porter

icebreaker operations

1000 m Imagery ©2012 Digital Globe, Inc.

Icebreaker

McMurdo Station

Presenter
Presentation Notes
Each day the icebreaker stands down saves NSF tens of thousands of dollars.
Page 8: Porter

landing site reconnaissance

Page 9: Porter

landing site reconnaissance

Page 10: Porter

100 m

geologic mapping

Imagery ©2012 Digital Globe, Inc.

Presenter
Presentation Notes
Geologic Mapping is a common application for commercial imagery at the poles.
Page 11: Porter

emperor penguin census

Imagery ©2012 Digital Globe, Inc.

Page 12: Porter

weddell seal counts

Page 13: Porter
Page 14: Porter
Page 15: Porter

elevation from stereoscopic images

Page 16: Porter

feature tracking

200m

Page 17: Porter

imagery challenges

• Challenges of the imagery

– Organization

• Over 500,000 images in polar regions

• Cataloging and space requirements

– Format

• Basic spatial location – not terrain corrected

• Uncorrected DN values

• 16 bit

• NITF format

Page 18: Porter

open source tools

• GDAL and Python

• Open source and freely available

• Large user base, well established and well maintained

– GDAL

• Libraries and standalone utilities

• Built with C++

• Becoming incorporated into other GIS packages

• Both raster and vector libraries

• Very good coordinate system support

• Bindings to Python, Ruby, Java, .NET, Perl, etc.

– Python

• GDAL Python bindings include Numpy/Scipy

• Benefits of GDAL with intuitive Python language

Page 19: Porter

open source tools

• What you can do:

– Organization

• Query rasters and vectors for geometry and metadata

• Write info to a spatial database or vector format

• Rename and archive imagery

– Manipulation

• Orthorectify, project, and manipulate imagery

• Stack bands, pan-sharpen, mosaic, or subset multiple images

• Run the package on a linux cluster for batch processing

• Easily schedule cron jobs for repeat updates

Page 20: Porter

imagery catalog

• Example 1: PGC imagery catalog

Page 21: Porter

imagery processor

Input Image DEM (optional)

gdalwarp : orthorectify with RPCs

gdal_translate and gdal_merge or VRT: rescale and stretch

RPCs (if external)

VRT

Reflectance Info

gdal_translate: write to output format

Output Image

gdaladdo: compute pyramids (if geoTiff)

• Example 2: PGC image orthorectification script

Page 22: Porter

radiometric correction

Page 23: Porter

radiometric correction

Page 24: Porter

Imagery Copyright Geoeye, Inc.

radiometric correction

Page 25: Porter

mosaics

Page 26: Porter

GDAL script by Dan Stahlke, GINA Imagery Copyright 2011 Digital Globe, Inc.

pan-sharpening

Page 27: Porter

GDAL script by Dan Stahlke, GINA Imagery Copyright 2011 Digital Globe, Inc.

pan-sharpening

Page 28: Porter

other applications

• Mosaics, Reprojection, and Subsetting: Rapidice.org

Page 29: Porter