cee 795 water resources modeling and gis learning objectives: demonstrate the concepts of spatial...

53
CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information Use raster and vector representations of spatial fields Perform raster calculations in hydrology Perform raster based watershed delineation from digital elevation models Handouts: Assignments: Exercise #3 Lecture 4: Spatial Fields and DEM Processing (some material from Dr. David Maidment, University of Texas and Dr. David Tarboton, Utah State University) February 6, 2006

Upload: roy-sullivan

Post on 13-Dec-2015

220 views

Category:

Documents


3 download

TRANSCRIPT

Page 1: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

CEE 795Water Resources Modeling and GIS

Learning Objectives:• Demonstrate the concepts of spatial fields as a way to represent geographical

information• Use raster and vector representations of spatial fields• Perform raster calculations in hydrology• Perform raster based watershed delineation from digital elevation models

Handouts: Assignments:Exercise #3

Lecture 4: Spatial Fields and DEM Processing

(some material from Dr. David Maidment, University of Texas and Dr. David Tarboton, Utah State University)

February 6, 2006

Page 2: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Vector and Raster Representation of Spatial Fields

Vector Raster

Page 3: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Numerical representation of a spatial surface (field)

Grid

TIN Contour and flowline

Page 4: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Six approximate representations of a field used in GIS

Regularly spaced sample points Irregularly spaced sample points Rectangular Cells

Irregularly shaped polygons Triangulated Irregular Network (TIN) Polylines/Contours

from Longley, P. A., M. F. Goodchild, D. J. Maguire and D. W. Rind, (2001), Geographic Information Systems and Science, Wiley, 454 p.

Page 5: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

A grid defines geographic space as a matrix of identically-sized square cells. Each cell holds a

numeric value that measures a geographic attribute (like elevation) for that unit of space.

Page 6: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

The grid data structure

• Grid size is defined by extent, spacing and no data value information– Number of rows, number of column– Cell sizes (X and Y) – Top, left , bottom and right coordinates

• Grid values – Real (floating decimal point)– Integer (may have associated attribute table)

Page 7: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Definition of a Grid

Numberof

rows

Number of Columns(X,Y)

Cell size

NODATA cell

Page 8: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Points as Cells

Page 9: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Line as a Sequence of Cells

Page 10: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Polygon as a Zone of Cells

Page 11: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

NODATA Cells

Page 12: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Cell Networks

Page 13: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Grid Zones

Page 14: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Floating Point Grids

Continuous data surfaces using floating point or decimal numbers

Page 15: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Value attribute table for categorical (integer) grid data

Attributes of grid zones

Page 16: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Raster Sampling

from Michael F. Goodchild. (1997) Rasters, NCGIA Core Curriculum in GIScience, http://www.ncgia.ucsb.edu/giscc/units/u055/u055.html, posted October 23, 1997

Page 17: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Raster Generalization

Central point ruleLargest share rule

Page 18: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Raster Calculator

Cell by cell evaluation of mathematical functions

Page 19: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Example

Precipitation-

Losses (Evaporation,

Infiltration)=

Runoff5 22 3

2 43 3

7 65 6

-

=

Page 20: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Runoff generation processesInfiltration excess overland flowaka Horton overland flow

Partial area infiltration excess overland flow

Saturation excess overland flow

PP

P

qrqs

qo

PP

P

qo

f

PP

P

qo

f

f

Page 21: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Runoff generation at a point depends on

• Rainfall intensity or amount

• Antecedent conditions

• Soils and vegetation

• Depth to water table (topography)

• Time scale of interest

These vary spatially which suggests a spatial geographic approach to runoff estimation

Page 22: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Modeling infiltration excessEmpirical, e.g. SCS Curve Number

method

Runoff from SCS Curve Number

0

1

2

3

4

5

6

7

8

9

10

0 2 4 6 8 10 12

Precipitation (in)

Dir

ect

Ru

no

ff (

in)

CN=10090 80

70

60

50

40

3020

S8.0P

)S2.0P(R

2

10CN

1000S

Page 23: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Cell based discharge mapping flow accumulation of generated runoff

Radar Precipitation grid

Soil and land use grid

Runoff grid from raster calculator operations implementing runoff generation formula’s

Accumulation of runoff within watersheds

Page 24: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Raster calculation – some subtleties

Analysis extent

+

=

Analysis cell size

Analysis mask

Resampling or interpolation (and reprojection) of inputs to target extent, cell size, and projection within region defined by analysis mask

Page 25: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Spatial Snowmelt Raster Calculation ExampleThe grids below depict initial snow depth and average temperature over a day for an area.

40 50 55

42 47 43

42 44 41

100 m

100

m

(a) Initial snow depth (cm)

4 6

2 4

150 m

150

m

(b) Temperature (oC)

One way to calculate decrease in snow depth due to melt is to use a temperature index model that uses the formula

TmDD oldnew

Here Dold and Dnew give the snow depth at the beginning and end of a time step, T gives the temperature and m is a melt factor. Assume melt factor m = 0.5 cm/OC/day. Calculate the snow depth at the end of the day.

Page 26: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Snow Depth and Temperature

4

2 4

640 50 55

4347

414442

42

Initial Snow Depth (cm) Temperature (º C)

100 m

100

m

150 m

150

m

Page 27: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

New depth calculation using Raster Calculator

[snow100m] - 0.5 * [temp150m]

Page 28: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

The Result

38 52

41 39

• Outputs are on 150 m grid.

• How were values obtained ?

Page 29: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Nearest Neighbor Resampling with Cellsize Maximum of Inputs

40 50 55

4347

414442

42

100

m

4

2 4

6150

m

40-0.5*4 = 38

55-0.5*6 = 5238 52

41 39

42-0.5*2 = 41

41-0.5*4 = 39

Page 30: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Scale issues in interpretation of measurements and modeling results

The scale triplet

From: Blöschl, G., (1996), Scale and Scaling in Hydrology, Habilitationsschrift, Weiner Mitteilungen Wasser Abwasser Gewasser, Wien, 346 p.

a) Extent b) Spacing c) Support

Page 31: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

From: Blöschl, G., (1996), Scale and Scaling in Hydrology, Habilitationsschrift, Weiner Mitteilungen Wasser Abwasser Gewasser, Wien, 346 p.

Page 32: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Spatial analyst options for controlling the scale of the output

Extent Spacing & Support

Page 33: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Raster Calculator “Evaluation” of temp150

4 6

2 4

6

2

4

44

4 6

Nearest neighbor to the E and S has been resampled to obtain a 100 m temperature grid.

2 4

Page 34: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Raster calculation with options set to 100 m grid

• Outputs are on 100 m grid as desired.

• How were these values obtained ?

38 52

41 39

47

41

42

4145

[snow100m] - 0.5 * [temp150m]

Page 35: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

100 m cell size raster calculation

40 50 55

4347

414442

42

100

m15

0 m

40-0.5*4 = 38

42-0.5*2 = 4138 52

41 39

43-0.5*4 = 41

41-0.5*4 = 39

47

41 45 41

42

50-0.5*6 = 47

55-0.5*6 = 52

47-0.5*4 = 45

42-0.5*2 = 41

44-0.5*4 = 42

Nearest neighbor values resampled to 100 m grid used in raster calculation

4

6

2 4

6

2

4

44

4

6

2 4

Page 36: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

What did we learn?

• Spatial analyst automatically uses nearest neighbor resampling

• The scale (extent and cell size) can be set under options

• What if we want to use some other form of interpolation?

Page 37: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

InterpolationEstimate values between known values.

A set of spatial analyst functions that predict values for a surface from a limited number of sample points creating a continuous raster.

Apparent improvement in resolution may not be justified

Page 38: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Interpolation methods

• Nearest neighbor• Inverse distance

weight• Bilinear

interpolation• Kriging (best linear

unbiased estimator)• Spline

ii

zr

1z

)dyc)(bxa(z

iizwz

ii eei yxcz

Page 39: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Nearest Neighbor “Thiessen” Polygon Interpolation Spline Interpolation

Page 40: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Spatial Surfaces used in Hydrology

Elevation Surface — the ground surface elevation at each point

Page 41: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

3-D detail of the Tongue river at the WY/Mont border from LIDAR.

Roberto GutierrezUniversity of Texas at Austin

Page 42: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Topographic Slope

• Defined or represented by one of the following– Surface derivative z (dz/dx, dz/dy)

– Vector with x and y components (Sx, Sy)

– Vector with magnitude (slope) and direction (aspect) (S, )

Page 43: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Standard Slope Function

a b c

d e f

g h i

cingx_mesh_spa * 8

i) 2f (c - g) 2d (a

dx

dz

acing y_mesh_sp* 8

i) 2h (g - c) 2b (a

dy

dz

22

dy

dz

dx

dz

run

rise

dy/dz

dx/dzatan

Page 44: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Aspect – the steepest downslope direction

dx

dz

dy

dz

dy/dz

dx/dzatan

Page 45: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Example30

80 74 63

69 67 56

60 52 48

a b c

d e f

g h i229.0

30*8

)2456*263()6069*280(

cingx_mesh_spa * 8

i) 2f (c - g) 2d (a

dx

dz

329.030*8

)6374*280()4852*260(

acing y_mesh_sp* 8

c) 2b (a -i) 2h (g

dy

dz

o8.21)401.0(atan

o8.34329.0

229.0atanAspect

o

o

2.145

180

145.2o

401.0

329.0229.0Slope 22

Page 46: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

80 74 63

69 67 56

60 52 48

80 74 63

69 67 56

60 52 48

30

45.0230

4867

50.0

30

5267

Slope:

Hydrologic Slope - Direction of Steepest Descent

30

ArcHydro Page 70

Page 47: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

32

16

8

64

4

128

1

2

Eight Direction Pour Point Model

ESRI Direction encoding

ArcHydro Page 69

Page 48: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

?

Limitation due to 8 grid directions.

Page 49: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Length on Meridians and Parallels

0 N

30 N

Re

Re

RR

A

BC

(Lat, Long) = (, )

Length on a Meridian:AB = Re (same for all latitudes)

Length on a Parallel:CD = R Re Cos(varies with latitude)

D

Page 50: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Example: What is the length of a 1º increment along on a meridian and on a parallel at 30N, 90W?Radius of the earth = 6370 km.

Solution: • A 1º angle has first to be converted to radians radians = 180 º, so 1º = /180 = 3.1416/180 = 0.0175 radians

• For the meridian, L = Re km

• For the parallel, L = Re CosCoskm• Parallels converge as poles are approached

Page 51: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Summary Concepts

• Grid (raster) data structures represent surfaces as an array of grid cells

• Raster calculation involves algebraic like operations on grids

• Interpolation and Generalization is an inherent part of the raster data representation

Page 52: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Summary Concepts (2)

• The elevation surface represented by a grid digital elevation model is used to derive surfaces representing other hydrologic variables of interest such as– Slope– Drainage area (more details in later classes)– Watersheds and channel networks (more details

in later classes)

Page 53: CEE 795 Water Resources Modeling and GIS Learning Objectives: Demonstrate the concepts of spatial fields as a way to represent geographical information

Summary Concepts (3)

• The eight direction pour point model approximates the surface flow using eight discrete grid directions.