Map Algebra and Beyond: Advanced topics and applications to Nexrad Xingong Li University of Kansas 5 November 2009

1. Map Algebra and Beyond:Advanced topics and applications to NexradXingong LiUniversity of Kansas5 November 2009

2. Major Extensions to Map Algebra • Scott (1999) extended the original 2D MA operations into three dimensional raster datasets (volumetric MA) • Solid earth • Atmosphere • Ocean • Li and Hodgson (2004) and Wang and Pullar (2005) developed MA operations for vector fields (cell values are vectors rather than scalars) • Aspect, surface normal • flow and wind fields • Mennis et al., 2005 developed cubic MA for spatio-temporal datasets where the third dimension is time • Spatio-temporal time series • French and Li (in press) proposed MA operations for the vector data model

3. Map Algebra for Vector Fields • Types of fields • Scalar fields—each cell stores a scalar value • Normal, ordinal, interval, and ratio • Vector fields—each cell stores a vector • 2D, 3D, Multi dimensional • Map algebra operations on vector fields

4. Aspect1 Aspect2 Mean Aspect How to calculate the mean aspect? MeanAspect = (Aspect1 + Aspect2)/2? What is the mean aspect within landuse (or elevation) zones?

5. Calculate Mean Aspect • What’s the mean aspect of 2 and 358 ? • (2 + 358) = 180 ? • Aspects are unit vectors • How to calculate mean aspects? • Vector algebra

6. y =330 =30 x Mean Aspect by Unit Vector

7. y A B  x Angular Between Two Vectors

8. Terrain Hillshade

9. Friction and Movement Direction The cost distance operation in ArcGIS assumes that friction is independent of movement direction (cost per unit distance)

10. Friction and Movement Direction

11. missing polygons Map Algebra for the Vector Data Model • No counterpart in the vector data model • Have to convert vector data into raster to use map algebra operations • Various problems during the conversion • Impose an arbitrary analysis resolution

12. Local Spatial Scope • A cell in the raster data model • A feature in the vector data model • Two types of vector layers (focus and value layer) • Each feature on the focus layer defines a local spatial scope of an operation • Value layer stores the features to which the features on the focus layer will be spatially compared • Focus and value layer can be the same

13. Local Scope

14. Focal Spatial Scope • Neighborhoods for points, lines, and polygons • Neighborhoods are not necessary polygons • Neighborhoods can be defined based topological relationships among features • Generic neighborhood could also be defined

15. Neighborhoodsfor points

16. Neighborhoods for lines

17. Neighborhoods for polygons

18. Zonal Scope • A collection of features with the same values for a given field • May become a local scope if each feature has a unique value in the field

19. Value Feature Selection and Adjustment • The value features and their attributes associated with a focus feature may partially overlap with the focus feature • Four selection/adjustment options • No adjustment on geometry and attribute • Only on geometry • On geometry and attribute (over value feature) • On geometry and attribute (over spatial scope)

20. Select Value Features • Value features are selected based on the dimensionally extended 9-intersection model (DE9IM) developed by Egenhofer and Herring (1991) and Clementini et al. (1993) • The ‘within’ relationship (“T*F**F***”) • Geometric types which can have the ‘within’ relationship

21. Attribute Adjustment OVER_VALUE_FEATURE OVER_LNZ

22. Operations

23. A Possible Syntax NewLayer = FocusLayer.Operation (Scope, ValueLayer, Attribute, Adjustment, Normalization)

24. An Implementation

25. Examples (a) NewLayer = Siren.Sum (Radial (0, 0, 0, X), CensusBlock, POP, OVER_VALUE_FEATURE). (b) NewLayer = SirenZone. Sum (Zonal(ID), CensusBlock, POP, OVER_VALUE_FEATURE).

26. Examples (b) NewLayer = Subwatersheds. Majority (Local(), RadarCells, PRECIP, ON_GEOMETRY, Area) (c) NewLayer = Subwatersheds. Sum (Local(), RadarCells, PRECIP, OVER_LNZ)

27. Comparison to Raster Map Algebra • Vector MA does not impose any arbitrary resolutions but simply maintain the original resolution of the data through its operations • Raster MA has difficulty handling the neighborhoods which are defined for individual features or are based on the topological relationships between features • The vector cartographic modeling is more appropriate for characterizing discrete features and the relationships among the features

28. Spatiotemporal Map Algebra Cubic local functions Cubic Focal functions Mennis, J., Viger, R., and Tomlin, D. 2005, “Cubic map algebra functions for spatiotemporal analysis”. Cartography and Geographic Information Science, 32(1): 17- 32.

29. Cubic Zonal Operations vary only in time vary only in space vary both in space and time

30. Antecedent Precipitation and Water Quality Explore the relationship between water sample quality and antecedent rainfall (precipitation occurred before water samples were taken)

31. Water Samples in Space and Time 1049 water samples were collected from 89 locations at different times (from 1992 to 1999)

32. Defining Spatiotemporal Zones 1227a 1224b Zone = flow length + antecedent time

33. Total Amount of Phosphorous vs. Antecedent Precipitations

34. From Spatio-temporal precipitation data to precipitation events (storms) • The Eulerian view focuses on the change of state in space • While a sequence of changes in space may portray the movement of an entity across the space, there is no explicit representation of those entities. • no structured data object representing "a storm“ • no explicit representation of behaviors that storms can exhibit. • The Lagrangian view offers an alternative perspective that focuses on movement and uses an object-based approach

35. Study Area and Data • The study domain is the ABRFC (Stage III and P1 NEXRAD products, 4 km spatial resolution, hourly in time) • The precipitation data span a period of 11 years from 10/01/1995 to 09/30/2006

36. NEXRAD (Next generation Radar) • About 150 stations covering the entire U.S. • Provides hourly precipitation estimate by combining radar, satellite, and rain gauge data • Spatial resolution is about 4 km

37. NEXRAD Data • Precipitation data are broken down into 13 separate geographical regions • Each region covers a NWS-designated river basin (River Forecast Center) • Temporal coverage of the dataset varies in each river basin • Data can be downloaded from the NOAA website or from individual RFCs

38. Time Series Data Animation

39. Storm (Event) Extraction • A storm (event) is defined as a contiguous precipitation object in space and time • a set of connected precipitation cells delineated from stacked hourly NEXRAD precipitation layers. • The algorithm is based on the component labeling algorithm in digital image processing • Controlled by 3 parameters • the minimum hourly precipitation (MHP) in a cell • the minimum time span (MTS) of a storm • the definition of spatial and temporal connectivity t-1 t t+1

40. Projected on x-y plane Projected on to the x-time plane Projected on to the y-time plane A Storm Example

41. Storm Tracking and Representation • A directed graph is used to represent a storm • Nodes are precipitation-weighted centroids of spatially contiguous areas receiving rainfall in each hour • Directed edges indicate spatial and temporal linkage (split or merge) among the rainfall areas during the life span of the event

42. Data Processing and Software Tools

43. Warm Season Storm Spatio-temporal Characteristics • Warm season: April to September • 04/01/96—09/30/2006 • 519,562 storms

44. Temporal Characteristics (Annual)

45. Temporal Characteristics (diurnal)

46. Spatial Characteristics Total number of storms that occurred during the 11 year period

47. Spatial Characteristics Total amount of storm precipitation in mm during the 11 year period

48. Spatial Characteristics • Precipitation-weighted centroids of the events were calculated and used to represent the events as points in space and time in storm density analysis • The number of events per km2 of the 11-year period • The amount of precipitation per km2 of the 11-year period

49. Storm Movement • Precipitation-weighted mean storm movement vector is calculated for each storm from the directed graph • All the data from 10/01/1995 to 09/30/2006 Length represents movement speed. Start point is precipitation-weighted centroid.

50. Storm Movement Directional distribution of storms (left) and storm precipitation (right)