SDMI/AlaskaMapped Imagery, DEM, DRG, and Charts Web Services with Desktop GIS Applications Tom Heinrichs AmericaView Fall Technical Meeting Madison, Wisconsin October 11, 2010 www.gina.alaska.edu www.alaskamapped.org

[email protected]

Alaska North Slope Oil Industrial Area Landsat 5: June 15, 2005 1

Landsat 5: 2004 vs. 2005 Alaska acquisitions

2004: 25 Scenes

2005: 3777 Scenes 2005-2009: >18,500 Scenes

NOAA-NESDIS captures; UAF-GINA transports; USGSEROS processes and archives; GINA interfaces with users

Anaktuvuk Fire Scar – June 14, 2008 – Landsat 5 (bands 432)

4

Goals for talk • Describe current use of web services in AmericaView member states • Discuss open standards vs. proprietary standards implementations

5

(An incomplete) Listing of States with Active Web Services Implementation • • • • • •

Alaska Maryland Michigan Texas Wisconsin Please tell me if you belong on this list

6

Alaska Statewide Digital Mapping Initiative www.alaskamapped.org The Alaska Statewide Digital Mapping Initiative’s primary goals are to acquire new and better maps for Alaska and to make existing map products more easily available. The SDMI seeks to make ongoing improvements to Alaska maps on a broad, statewide scale. Alaska is the only state in the nation lacking current, accurate, high-resolution maps.

http://browse.alaskamapped.org

8

More Alaska Data….

9

Statewide Digital Mapping Initiative – Web Services • What is an OGC Web Mapping Service? – Open standards service that dynamically produces a spatially referenced image for use in mapping applications – http://wms.alaskamapped.org//cgibin/bdl.cgi?VERSION=1.1.1&REQUEST=GetMap&SRS=EPSG:102006&BB OX=137980.289057582,1070139.9759959,146394.055885116,1076807.489 33093&WIDTH=640&HEIGHT=500&LAYERS=bdl_high_res&STYLES=&EX CEPTIONS=application/vnd.ogc.se_xml&FORMAT=image/png&BGCOLOR= 0xFFFFFF&TRANSPARENT=TRUE – – – – –

HTTP and XML Projection Geographic extent Resolution Image Format

– All of the above is handled on the server 10

Data Gallery www.alaskamapped.org/data/gallery

13

The bdl_low_res layer contains the Blue Marble Next Generation July dataset. It is intended to be viewed at scales greater than 250m per pixel 14

The bdl_mid_res layer contains a Landsat mosaic generated by the USGS from data acquired from 1988 to 1992 and a simulated true colour Landsat mosaic provide by Earthstar Geographics LLC. It is intended to be viewed at scales greater than 15m per pixel. Imagery copyright Earthstar Geographics (http://www.terracolor.net). 15

The bdl_high_res layer contains data from a variety of sources and is intended to be viewed at scales from 25m per pixel to 1m per pixel. Images courtesy of Digital Globe, USDA-NRCS, USDA-FS, USGS, BLM, NPS, FAA, Alaska DNR, DMVA, DCED, DOT&PF 16

NOAA Charts The Charts service is meant to provide a unified view of the RNC data, without the collar information. The service is designed to only display charts that are meaningful at the requested scale. 17

DRGs This group provides four layers, “drg_24k”, “drg_25k”, “drg_63k”, “drg_250k”, which contain the USGS topographic maps for the scales 1:24,000, 1:25,000, 1:63,360, and 1:250,000 respectively. 18

Large-scale Image Processing • Enhancement of statewide topo maps – 3167 quad maps – 63,000 megapixels – 1.7 million sq-km

• Processed at ARSC on NSF-supported PACMAN system – Original USGS scans – Sharpened lines and text – Hillshaded – Done on a large-scale; hands-off processing; in ARSC-managed supercomputing queue

www.gina.alaska.edu www.alaskamapped.org

Canadian DRGs 1:250,000 scale topographic covering Canada. For more information, please see Natural Resources Canada Mapping Services . 20

Landsat Pan This layer group provides a single banded dataset derived from Landsat 7’s pan-chromatic band. It was generated from the Landsat 7 datasets provided by NASA’s OnEarth. It is broken up into two layers to increase the rendering 21 speed.

Shaded Relief NED This group has one layer which provides a shaped relief dataset derived from the National Elevation Dataset (NED). 22

Realtime MODIS and AVHRR Realtime optical and infrared satellite data captured in Fairbanks and Barrow by UAF-GINA and NOAA-NESDIS FCDAS. 23

2 arc second NED Standard National Elevation Dataset (NED) for Alaska.

24

ASTER GDEM – 1 arc second

25

SRTM 1 arc second Shuttle Radar Topography Mission

26

NED -- 1/3 arc second Airborne IfSAR For NPR-A.

27

Chugach SPOT DEM This layer is named “Chugach Spot DEM” and contains a 20 meter resolution DEM covering the greater Chugach National Forest area. The vertical datum is “altitude over GEM 96”. Source: U.S. Forest Service 28

LIDAR – various patches around the state Kenai Peninsula piece shown here

29

1A

1B

1C

1D

Figure 1. Samples of the eMODIS Alaska product suite. (1A). RGB image created with Bands 1, 4, 3 with 500-m data acquired in August 2006. (1B). Grey-scaled 250-m NDVI calculated with data from the same time period. (1C) depicts the per-pixel day of acquisition during the August interval. (1D) is an example of per-pixel quality taken from a June 2003 interval. Credit: eMODIS ALASKA, Calli B. Jenkerson & Gail L. Schmidt, USGS, ASPRS 2009 Annual Conference , March 9-13, 2009 30

kml.gina.alaska.edu

31

ESRI ArcGIS Layer Files (.lyr)

32

Infrastructure • We have built a fast-performing infrastructure to make all this work well. • Our WMS, WCS,and WFS services are powered by a combination of Mapserver, GDAL/OGR, and Proj4 . We do extensive data preparation work using scripts written in Ruby and keep our inventory information in a PostgreSQL database with the PostGIS spatial extensions. The vast majority of our data being served is stored as tiled, deflate compressed, GeoTiffs. • We served data to tiled base mapping tool kits using a custom tile server application. It is written in Ruby and uses Mongrel in combination with Rack as the web server engine. Image processing is performed with the Ruby language bindings to ImageMagick.

33

EPSG Code Datum

Projection

EPSG:102006 EPSG:3338 EPSG:2964 EPSG:102009 EPSG:3031 EPSG:4326 EPSG:900913

NAD 83 NAD 83 NAD 27 NAD 83 WGS 84 WGS84 Spherical

Alaska Albers Alaska Albers Alaska Albers North America Lambert Conformal Conic Polar Stereographic Geographic Google Mercator

EPSG:102631 EPSG:102632 EPSG:102633 EPSG:102634 EPSG:102635 EPSG:102636 EPSG:102637 EPSG:102638 EPSG:102639 EPSG:102640 EPSG:26701 EPSG:26702 EPSG:26703 EPSG:26704 EPSG:26705 EPSG:26706 EPSG:26707 EPSG:26708 EPSG:26709 EPSG:26901 EPSG:26902 EPSG:26903 EPSG:26904 EPSG:26905 EPSG:26906 EPSG:26907 EPSG:26908 EPSG:26909 EPGS:32606 EPSG:26936 EPSG:26937 EPSG:26938 EPSG:26939

NAD 83 NAD 83 NAD 83 NAD 83 NAD 83 NAD 83 NAD 83 NAD 83 NAD 83 NAD 83 NAD 27 NAD 27 NAD 27 NAD 27 NAD 27 NAD 27 NAD 27 NAD 27 NAD 27 NAD 83 NAD 83 NAD 83 NAD 83 NAD 83 NAD 83 NAD 83 NAD 83 NAD 83 WGS 84 NAD 83 NAD 83 NAD 83 NAD 83

NAD 1983 StatePlane Alaska 1 FIPS 5001 Feet NAD 1983 StatePlane Alaska 2 FIPS 5002 Feet NAD 1983 StatePlane Alaska 3 FIPS 5003 Feet NAD 1983 StatePlane Alaska 4 FIPS 5004 Feet NAD 1983 StatePlane Alaska 5 FIPS 5005 Feet NAD 1983 StatePlane Alaska 6 FIPS 5006 Feet NAD 1983 StatePlane Alaska 7 FIPS 5007 Feet NAD 1983 StatePlane Alaska 8 FIPS 5008 Feet NAD 1983 StatePlane Alaska 9 FIPS 5009 Feet NAD 1983 StatePlane Alaska 10 FIPS 5010 Feet UTM zone 1 UTM zone 2 UTM zone 3 UTM zone 4 UTM zone 5 UTM zone 6 UTM zone 7 UTM zone 8 UTM zone 9 UTM zone 1 UTM zone 2 UTM zone 3 UTM zone 4 UTM zone 5 UTM zone 6 UTM zone 7 UTM zone 8 UTM zone 9 WGS 84 / UTM zone 6N NAD83 / Alaska zone 6 NAD83 / Alaska zone 7 NAD83 / Alaska zone 8 NAD83 / Alaska zone 9

34

Projections currently offered through GINA web services. We can add more. Just ask.

Send us the EPSG code or your .prj file. Note: In 9.3 desktop client will allow you request data in a projection not supported by the server. However, the reprojection is not handled optimally by the client. Q: Which polar projections are you using?

GINA Tile Services Web mapping toolkits • Google Maps • Google Earth • Open Layers • ArcGIS API for JS • Virtual Earth • Yahoo Maps • Bing • and many others…. 35

36

Sample ESRI Flex Application Uses GINA tiles

37

Sample ESRI Flex Application Uses GINA tiles

38

Google Maps API

• GINA/SDMI layers encapsulated in a one-line include 39

Google Maps API

• GINA/SDMI layers encapsulated in a one-line include 40

Google Maps API

PWS Herring Portal: http://www.pwsherringportal.org/herringportal/ 41

by Axiom

Another API

NRCS Web Soil Survey: http://websoilsurvey.nrcs.usda.gov/ 42

Java API

DNR AlaskaMapper: http://mapper.landrecords.info/ 43

CIMSS – WisconsinView Met, Wx Radar, Imagery mashup based on WMS

44

SSEC Smart Phone Capability – Prototypes

46

47

48

Open vs Proprietary Standards Sample ESRI Flex Application – Uses ESRI REST API

A n OGC WMS call: http://wms.alaskamapped.org//cgibin/bdl.cgi?VERSION=1.1.1&REQUEST=GetMap&SRS=EPSG:102006&BBOX=137980.289057582 ,1070139.9759959,146394.055885116,1076807.48933093&WIDTH=640&HEIGHT=500&LAYERS= bdl_high_res&STYLES=&EXCEPTIONS=application/vnd.ogc.se_xml&FORMAT=image/png&BGCO LOR=0xFFFFFF&TRANSPARENT=TRUE 49

50

(An incomplete) Listing of States with Active Web Services Implementation • • • • • •

Alaska Maryland Michigan Texas Wisconsin Please tell me if you belong on this list

51