1 / 31

A Data System for Visualizing 4-D Atmospheric CO 2 Models and Data

A Data System for Visualizing 4-D Atmospheric CO 2 Models and Data. Tyler A. Erickson, Ph.D. Research Scientist Adjunct Assistant Professor of Civil & Environmental Engineering Michigan Technological University, Michigan, USA 22 October 2009. Collaborators.

sileas
Télécharger la présentation

A Data System for Visualizing 4-D Atmospheric CO 2 Models and Data

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. A Data System for Visualizing 4-D Atmospheric CO2 Models and Data Tyler A. Erickson, Ph.D. Research Scientist Adjunct Assistant Professor of Civil & Environmental Engineering Michigan Technological University, Michigan, USA 22 October 2009

  2. Collaborators Prof. Anna M. MichalakUniversity of MichiganAnn Arbor, Michigan, USA Carbon Cycle Science Researcher Prof. John C. LinUniversity of WaterlooWaterloo, Ontario, Canada STILT Atmospheric Transport Model Creator

  3. Time is Important Geospatial representation of present conditions is fine,but predicting future conditions is really useful and interesting...

  4. Temporal Change is Everywhere Land cover change Climate Change Disease Spread Environmental Change Economic Change Photo credit: John McColgan of the Bureau of Land Management, Alaska Fire Service

  5. Spatial-Temporal Data radio collar(a.k.a. caribou bling) Location: North Slope of Alaska, USASource: The National Academies Location: Central CanadaSource: NASA

  6. Problem #1: Too Few Data! Collecting data with in-situ sensors is expensive Even with dense meshesof sensors, processes are severely undersampled Soddie Meteorological TowerLocation: Niwot Ridge LTER, Colorado, USASource: Tyler Erickson

  7. Solution #1: Model it! “Give me a modeland data with which run it,and I shall estimate allthe properties of the worldin both space and time.” - Archimedes of Syracuse(severely paraphrased)

  8. Spatial-Temporal Modeling DATA(x,t) = MODEL(x,t) + ERROR(x,t) Models should rigorously represent reality, if possible Errors near each other are often similar(i.e. geostatistics)

  9. Problem #2: Too Much Data! Data volumes are overwhelming! How do you go about exploringthe data in space and time? Soddie Meterological TowerLocation: Niwot Ridge LTER, Colorado, USASource: Tyler Erickson

  10. Solution #2: FOSS4G !!!

  11. CO2 Monitoring Source: http://commons.wikimedia.org/wiki/File:Mauna_Loa_Carbon_Dioxide.png

  12. Carbon Balance Source: NASA

  13. Atmospheric Carbon Monitoring NOAATall Tower Adapted from work by: K. Mueller, University of Michigan

  14. By the Numbers... A Typical Particle Simulation: 500*24 per day, simulated particles 10*24 hours of simulation per particle 6 positions per hour per particle 30 days, total dataset length TOTAL: ~500 million records per measurement tower

  15. How does this data get reported? Source: Lin, J.C. et al., 2003. JGR (Atmospheres) Figure 7

  16. Another Approach IDEA: Provide model results in user-friendly, standard data formats OGC KML 2.2 standard

  17. What's under the Hood? Virtual Globe ClientApplication KML GeospatialServer SciPy/NumPy libkml pylibkml DataStorage

  18. GeoDjango + = “A world-class geographic web framework” • A geospatial extension of a web framework designedfor publishing on-line newspapers (Django)created by a law student (Justin Bronn) • Python-based • Leverages GEOS, GDAL, Proj.4, and PostGIS • Templates allow for output in various geospatial formats Image Sources: Wikipedia and Wikimedia Commons

  19. pylibkml pylibkml is a Python wrapper for the libkml C++ library Allows for easy programmatic creation of valid KML documents http://code.google.com/p/pylibkml/ from pylibkml import Kml placemark = Kml().create_placemark({ 'name' : 'FOSS4G 2009', 'description' : 'In Sydney', 'timestamp' : {'when': '8/22/2009'}, 'point' : Kml().create_point({ 'extrude' : True, 'altitudemode' : 'relativetoground', 'coordinates' : Kml().create_coordinates( 151.1998,-33.8761), }) })

  20. Why Google Earth? Free (for some uses) Open Source Easy to Use Interface Rich Reference Imagery Wide User Base Runs on Linux Full KML Implementation Talks to External Servers Source: http://www.flickr.com/photos/gillpoole/

  21. Google's KML in Research Competition Earlier this year Google hosted a contest onusing KML to communicate scientific research Judging criteria: Usability Educational value Visual/interactive appeal Efficiency Attribution KML output from this data system was selected as one of 5 professional winners

  22. Tall Tower Measurements NOAATall Tower Adapted from work by: K. Mueller, University of Michigan

  23. Sensitivity to Surface Flux NOAATall Tower Adapted from work by: K. Mueller, University of Michigan

  24. Simulated Particle Tracks 29

  25. What Else Goes Under The Hood? ASCII Grid netCDF GeoTiff ESRIShapefile KML CONSUME (Fuel Consumption & Emissions Model) SciPy/NumPy libkml WKT Raster

  26. Modeling Wildfire Emissions

  27. WKT Raster Adds raster-vector spatial analysis to PostGIS http://trac.osgeo.org/postgis/wiki/WKTRaster 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 2 2 2 2 2 2 0 0 1 2 2 2 2 0 0 1 1 0 0 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 0 0 1 1 1 1 1 0 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 1 1 1 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 0 1 1 2 2 2 2 2 0 1 1 1 0 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 0 0 0 0 1 1 1 0 0 1 1 1 1 1 1 1 1 2 2 2 2 2 2 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 Example:Intersection(geometry,raster) → raster and = b a Source: WKTSpecifications1.0.ppt (Pierre Racine)

  28. Next Steps Atmospheric Carbon Application: move from prototypeto a web-accessible tool for researchers Allow users to upload data from atmospheric transport model runs Improve KML styling Create visualizations of additional high-dimensional datasets

  29. Questions? Tyler A. Erickson, Ph.D. Email: tyler.erickson@mtu.edu Web: http://people.mtri.org/tyler+erickson Twitter: tylerickson Code: http://bitbucket.org/tylere/geodjango-stilt/

More Related