Geospatial
Visualization
School of Architecture
Dr. Timothy Norris : tnorris@miami.edu
Christopher Mader : cmader@med.miami.edu
UM GIS Day
Wednesday November 13th 1:30-5 p.m.
Today's Goals
  1. Vector Analysis and Layout Export (review)
  2. Field Presentations
  3. Overlays and Interpolation Overview and Exercise
  4. Little River Critique and Review
So Many Tools ....
!CLICK ME!
Vector Analysis
  • Vector Operations (one layer)
    Buffer and Dissolve
  • Vector Operations (two layers)
    Spatial Join
  • Attribute Statistics & Aggregation
One Layer Vector Operations
  • Buffer
  • Dissolve
  • Transformations (vector to raster)
  • Convex Hull
  • Centroids
  • Transformations (polygon to line or vice versa)
  • Voroni Polygons (fill plane from points)
One Layer Vector Operations - Buffer
  • ring buffer (one radius)
  • multiple-ring buffer (several radii)
More Buffers
  1. variable width
  2. multiple ring
  3. doughnut
  4. setback
  5. non-dissolved
  6. dissolved
One Layer Vector Operations - Dissolve
  • combines features based on attribute
  • useful for data aggregation
  • useful for cartography
    • multiple line segments with same label (roads)
    • remove duplicate polygon labels for area
Dissolve Exercise - Create Simple Basemap

Dissolve Miami-Dade Streets layer

  • Create project folder and start new project
  • From the Box folder download Streets.zip and Water_P.zip; unzip into the project data folder
  • Add both layers to your project and put streets on top of water
  • Symbolize streets as unique values based on CLASS (0-1 freeways, 2-4 main roads, 5-9 streets)
  • Label streets as with classes for CLASS (use SQL to select certain streets and label)
  • Use the Dissolve Tool to dissolve streets on PRE_DIR, ST_NAME, ST_TYPE, SUF_DIR, CLASS
  • Import symbology from Streets to new dissolve layer
  • Explore differences ..
Dissolve Exercise - Export to aix file and open in Illustrator

Create Layout and Refine

  • Create layout zoomed in to Miami River
  • Set position and scale of map by numbers
  • Add scale and north arrow
  • Symbolize the Water_P layer to have a thin blue stroke, zoom in to check
  • Export layout as aix
  • With the Maps for Adobe Creative Cloud plugin installed and the correct version of Illustrator ...
  • Open the aix file in Illustrator
  • Explore the layers
Vector Data Overlays
  • Data Structures: Connection and Transformation
    • Overlay tools
      • Intersect
      • Union
      • Subtract
      • Symmetric Difference
      • Identity
      • Cover
    • Clipping
      - Sub-setting data with polygons
      - Raster or vector
image credit: Overlay. wiki.GIS.com
Intersection of Two Vector Layers
  • Only geometries and attributes that overlap are preserved
    - Attributes from both layers copied to new feature
    - Output geometry is lowest dimension of input geometries
  • Ratio (ArcGIS) or Overlap (QGIS) analysis
    - Attributes calculated as percentage of original area for each new polygon
  • Two rectangles in one layer and one circle in second layer
    - Three features in two layers to begin
    - Two features in one layer as a result
  • Similar to logical AND
Intersect (Analysis). ArcGIS Pro Help.Calculating percentage of overlap of two layers in QGIS 3. GIS StackExchange
Union of Two Vector Layers
  • All geometries and attributes preserved
    - FID's are labelled according to layer names
    - All features are of original size or smaller (never larger)
  • Ratio (ArcGIS) or Overlap (QGIS) analysis
    - Attributes calculated as percentage of original area for each new polygon
  • Two rectangles in one layer and one circle in second layer
    - Three features in two layers to begin
    - Five features in one layer as a result
  • Similar to logical OR
Union (Analysis). ArcGIS Pro Help.
Overlay. wiki.GIS.com
Erase (Subtract) of Two Vector Layers
  • Only geometries and attributes from input preserved
    - All features are of original size or smaller (never larger)
  • Ratio (ArcGIS) or Overlap (QGIS) analysis
    - Attributes calculated as percentage of original area for each new polygon
  • Four polygons one layer and one square in second layer
    - Five features in two layers to begin
    - Four features in one layer as a result (with multipart)
  • Similar to logical AND NOT
Erase (Analysis). ArcGIS Pro Help.
Less common Operations on Two Vector Layers
  • Symmetric Difference
    - All geometries and attributes that do not co-exist are preserved
    - Similar to logical XOR
    - (A union B) subtract (A intersect B)
  • Identity
    - Extent of input with geometry and attributes from both layers merged
    - (A subtract B) union (A intersect B)
  • Cover or Update
    - Similar to union in extent, but only preserves input geometry and attributes if feature not in update, otherwise copies update attributes to new geometry
    - A union (B subtract A)


Overlay. wiki.GIS.com
Identity (Analysis). ArcGIS Pro Help.
Update (Analysis). ArGIS Pro Help
Symmetrical Difference. ArcMap Help
Clip Vector with Vector (crop)
  • Only geometries and attributes from input preserved
    - All features are of original size or smaller (never larger)
  • Both input and clip features can be points lines or polygons
    - Be careful of unexpected results
  • Four polygons and two points in one layer and one circle in clip layer
    - Seven features in two layers to begin
    - Four features in one layer as a result (with multi geometries)
Clipping line layer based on polygons with QGIS? GIS StackExchange.
Clip (Analysis). ArcGIS Pro Help.
Clip Raster with Vector
  • Must use one vector feature to clip
  • Make sure projections are the same
    - It will function across projections, but makes computer work harder
  • Does not have to be rectangle
    - Hence not "crop"
Clip Raster (Data Management). ArcGIS Pro Help.
Clipping raster with vector boundaries using QGIS. GIS StackExchange.
Interpolation

Moving data from vector to raster structures

  • (review) Raster to Vector : Extract or Zonal Statistics
    • Extracts data from continuous surface to discrete features
  • Vector (point/line) to Raster: Interpolation
    • Interpolates one or more specific data values
    • Fills in the gaps between discrete features on a continuous surface
  • Vector to Raster (point/line/polygon): Rasterize
    • Converts vector geometries to raster image
    • Leaves gaps in raster where there is no discrete feature
Interpolation Exercise

Estimating PM2 concentrations for the continental United States

  1. Creating x/y data from a csv (review)
  2. Inverse Distance Weighting Interpolation
  3. Spline Interpolation
  4. Kriging
  5. Clipping interpolation with boundary polygon
  6. Clipping raster with boundary polygon
  7. Raster algebra (review)
Interpolation Exercise

Download from indicated source or Box (epa.zip in data folder)

  1. EPA AQS PM2.5 Annually Summary data 2022
    annual_conc_by_monitor_2022.csv (create x/y points from lat/lon - NAD83!!)
  2. USA nation outline TIGER file (geography division, census bureau)
    cb_2018_us_nation_5m.shp
  3. USA major cities (ESRI living atlas, can be data as a service)

Load into ArcGIS Pro and style as needed

Inverse Distance Weighting (IDW) Interpolation

In tool box under Spatial Analyst

  1. Take a moment to review the Interpolation toolbox
  2. Open IDW tool
    • Input point features: EPA data
    • Z_value: Arithmetic Mean (what is Z-Value?)
    • Output Cell Size: 0.183071288 (what is this?)
    • Search radius: variable
    • Number of Points: 12

Spline Interpolation

In tool box under Spatial Analyst

  1. Open Spline tool
    • Input point features: EPA data
    • Z_value: Arithmetic Mean
    • Output Cell Size: 0.183071288
    • Spline Type: regularized
    • Weight: 0.1
    • Number of Points: 12

Kriging Interpolation
  1. Open Kriging tool
    • Input point features: EPA data
    • Z_value: Arithmetic Mean
    • Kriging Method: Ordinary
    • Sevivariogram Model: Spherical
    • Lag size: 0.183071288
    • Number of Points: 12
    • Output Cell Size: 0.183071288
    • Search Radius: Variable
    • Number of Points: 12

Compare Different Interpolations

Clip all three rasters to extend of continental USA

  1. edit cb_2018_us_nation_5m.shp and delete Alaska, Hawaii, and Puerto Rico (review)
  2. Open clip tool (analysis tools)
    • Input Raster: an interpolated raster layer
    • Output extent: cb_2018_us_nation_5m
    • [x] Use Input Features for Clipping Geometry
    • When complete, copy symbology from source to new layer
      hint: right click new layer, select symbology and use hamburger menu to "Import from Layer"

Compare Different Interpolations

Use Raster Calculator to compare values (review)

  1. Open Raster Calculator (spatial analyst tools)
  2. Map algebra expression: kriging - spline
  3. Symbolize with diverging color ramp
  4. Interpret ....