Lecture 9. Raster Data Analysis. Tomislav Sapic GIS Technologist Faculty of Natural Resources Management Lakehead University

Transcription

1 Lecture 9 Raster Data Analysis Tomislav Sapic GIS Technologist Faculty of Natural Resources Management Lakehead University

2 Raster Data Model The GIS raster data model represents datasets in which square cells of the same size hold numeric values representing features or phenomena within a geographic area. Raster datasets extents are of a square or a rectangle shape in which the cells form horizontal (east-west) rows and vertical (north-south) columns. The structure of the raster and its cells allow for the cell positions to be calculated by knowing the position (x, y coordinates) of one of the raster corners. Raster cells can hold only numeric values, either integer or decimal. Raster s data type and data depth are set for the entire raster cannot be changed in totality or partially. Tables attached to integer rasters can hold text attributes as well. Because of a complete and continuous coverage of a geographic area raster data model lends itself to topological analyses.

3 GRID Selected Raster File Formats in ArcGIS Proprietary ESRI raster file format Composite files that are structured on a principle of a workspace and have part of their content stored in their own folder and part of it in the INFO folder. Can contain either integers decimal numbers (floating points) at 32-bit. GRID folder, containing parts (cell values, etc.) of the grid file. INFO folder, containing parts (statistics, etc.) of all GRIDs in the same workspace (folder)

4 Selected Raster File Formats in ArcGIS TIFF Tagged Image File Format Supports 1-bit unsigned and 4-, 8-, 16-, 32-bit signed/unsigned integer, and 32-bit floating point. 8- and 16-bit unsigned integer files allow for multiple bands. IMG ERDAS IMAGINE remote sensing software proprietary file format. Supports 1-, 2-, and 4-bit unsigned integer, 8-, 16-, 32-bit signed/unsigned integer, 32- and 64-bit decimal numbers. Multiple bands supported at all data depths.

5 ArcGIS raster dataset properties of the same dataset in three different formats. GRID TIFF IMG

6 Raster Data Analysis Analysis and manipulation of rasters can be divided along different lines, depending on which aspect of the raster use is considered: o Per-cell versus multiple cells analysis. o Single raster versus multiple rasters analysis. o Raster only analysis versus an analysis with a combination of a raster and a vector dataset. o Mathematical raster operations versus logical operations. o Generic versus applied algorithms.

7 Cell Statistics (per-cell, multiple rasters, generic math. algorithm) Available operations: minimum, maximum mean, median, range, sum, standard deviation, majority, etc. Block Statistics (multiple cells, single raster, generic math. algorithm) Zonal Statistics (multiple cells, multiple rasters, generic math. algorithm) Learn more about the functions by searching for them at

8 Reclassify (per-cell, single raster, value sorting) Con (per-cell, single or multiple rasters, logical algorithm) Con (in_conditional_raster, in_true_raster_or_constant, {in_false_raster_or_constant}, {where_clause}) Learn more about the functions by searching for them at 0.3/tools/spatial-analyst-toolbox/con-.htm l

9 Raster Calculator (per-cell, single or multiple rasters, generic math. algorithm) Learn more about the functions by searching for them at 0.3/tools/spatial-analyst-toolbox/rastercalculator.htm

10 Filter (multiple cells, single raster, generic math. algorithm) low-pass filter smoothes out data; equivalent to the Mean option of Focal Statistics Learn more about the functions by searching for them at 0.3/tools/spatial-analyst-toolbox/howfilter-works.htm l

11 Filter (cont d) The filter high-pass filter accentuates edges in data. The data The calculation for the processing cell (number 8) Raster example:

12 Extract by Points (per-cell, raster and vector combination) Extracts cell values to the points that fall within the cells. IDW (Inverse Distance Weighted) (per-cell interpolation of points, raster and vector combination) IDW function works on a principle that the interpolated variable decreases in influence with distance from its sampled location. Learn more about the functions by searching for them at 0.3/tools/3d-analyst-toolbox/how-idwworks.htm The Power parameter is used to increase the influence.

13 Different Interpolations Methods Topo to Raster Should be used to interpolate terrain if contours, rivers and lakes datasets are available. Ensures that lakes are flat and that rivers have a continuous downward slope in the direction of the flow. Spline Creates smooth surfaces, however, Continues interpolated trends outside of the range of observed values. IDW (Inverse Distance Weighted) Applies the IDW technique, the assumption of which is that the interpolated values is more close to the observed value, the more closer it is to it geographically, and vice versa. Should be only used on dense and evenly distributed observed (sample) values otherwise, tends to create visible rings of interpolated values. Radial Basis Function (Completely Regularized Spline version) Can be found in the Geostatistical Analyst extension in ArcGIS Creates smooth surfaces but that stay within the observed values range. Kriging Geostatistical method others above are deterministic interpolation methods. As such, it can produce a prediction surface and provide some measure of the certainty or accuracy of the predictions.

14 Fill (multiple cells, single raster, generic math. algorithm) Fills out sinks in the DEM raster, commonly created due to the resolution of the data or rounding of elevations to the nearest integer value. A sink is a cell surrounded by cells with higher values. Should be used to rectify a DEM raster and create a new one ready for hydrology tools such as Flow Direction. ArcGIS 10.1 Help File Sink Filled sink

15 Cost Path (multiple cells, multiple rasters, generic math. algorithm) Calculates the least-cost path from a source to a destination. The calculation is based on the values (costs) assigned to cells between the source(s) and the destination. The size of the cell and the side or diagonal position with respect to the direction of the cost calculation are also included in the calculation. There can be more than one source cell. The costs can be additionally weighted. Learn more about the functions by searching for them at p/10.0/help/index.html

16 Cost Distance one of the inputs for Cost Path is Cost Distance Learn more about the functions by searching for them at p/10.0/help/index.html Diagonal movement is increased by the proportional increase in distance versus the side movement, based on the Pythagoras theorem: = sq.rt. (1 + 1)

17 Cost Distance The final Cost Distance raster is calculated by spreading, through iterations, with distance costs, cell by cell, between the source (destination) cells and across the cost raster. Final Cost Distance raster Learn more about the functions by searching for them at ktop/10.0/help/index.html

18 Flow Direction (multiple cells, single raster, applied to hydrology) Working on a terrain (elevation) raster, creates a new raster with cell values indicating flow direction from each cell to its steepest downslope neighbour. Directions are indicated through codes based on a binary increment. The drop in the diagonal direction is divided by This affects the calculated steepness of a slope since it depends on height and distance (rise and run) and therefore the final direction (e.g. check the direction for the cell with the value of 47, in row 5, column 3). Flow directions Learn more about the functions by searching for them at p/10.0/help/index.html Outside of the main flow directions, additional codes are used to indicate directions where there are multiple cell neighbours with the same slope steepness.

19 Flow Accumulation (multiple cells, single raster, applied to hydrology) Based on the Flow Direction raster, a Flow Accumulation raster is created in which cell values indicate the number of cells that flow in that particular cell. Flow accumulation values can then be used to identify, e.g., streams (areas with high accumulation). Flow directions Learn more about the functions by searching for them at p/10.0/help/index.html

20 Viewshed (multiple cells, multiple rasters or a combination of a raster and a vector file, applied to viewscaping) Based on the terrain (elevation, surface) raster and observer points (raster cells), a new raster is created in which cell values indicate whether particular cells can be seen from the observer points or not. Learn more about the functions by searching for them at p/10.0/help/index.html The output values indicate whether the particular cell is visible from any of the observer points or not (Value = 0) and if it is, from how many (Value > 0) not possible to tell, though, from which observer points.