Module 06
Visualizing change (with point clouds)
Yesterday in the Great Smokies
Webcam
archive
GEO 409: Advanced GIS
Pauer Cartography Award
Submissions are open!
Cash prize
Traditional cartography, 3D mapping, and web maps
Let me know if you want to submit.
Videos
Modules 05 and 06
have been posted.
POSTPONED
Sites
dtabeyrathne
garnieri
katyeberry
dandelaneyuky
crkr226
brianemetz
skparker1
LucasRuckdeschel
sgsa228
lsisler19
sbti226
CSWALLiN
mhwh225
Review
First goal
How much of x is in (or near) y?
Highly accurate and precise data now.
Essential GIS.
Second goal
Where x
n
is in y, z exists.
x
n
is any number of variables.
z is not directly observed, but predicted.
Spatial modeling.
Current goal
Where x
t1
!= x
t2
, e occurred.
x
t
is same variable at different times.
!= is not equal to – they are different.
e is an event that caused the difference.
Change
There is nothing permanent except change.
– Heraclitus, 500 BCE
We will be known forever by the tracks we leave.
– Dakota Tribe
Land use change
NLCD captured every ~five years (2001-2019)
Raster cells are 30m x 30ma and align across years.
Class codes (land cover types) are the same across years.
c
t1
!= c
t2
if true, change occurred.
In-class task
Open the explore-change.ipynb Notebook.
Find how much land was developed between 2001 and 2019
within 20 miles of your hometown.
Mass terrain modification
Dam building.
Mountain Top Removal (MTR)
.
Surface mining creates overburden
which becomes valley fill.
Explore
Google Maps
Elevation change
Lidar derived DEMs from 2010 to current.(>5ft resolution)
USGS DEMs from the early 1990s. (10m resolution)
change = DEM
t1
- DEM
t2
If it is not close to 0, change occurred.
In-class task
Find a surface mine site
and show the change in elevation between the DEMs
within two miles of the center of the site.
Above-ground change
Compare lidar Phase 1 to Phase 2
via their DSMs.
change = DSM
t1
- DSM
t2
You'll need two point clouds to do this.
Visualization
Colorize Phase 1 point cloud by amount of removal
and Phase 2 by amount of addition.
Show both simultaneously.
🤯
Colorize
Point cloud can be colored by any RGB image.
Export raster and force rendering to RGB.
Sunrise spring equinox
Create hillshade from DSM
enabling shadows.
[Hillshade function](https://pro.arcgis.com/en/pro-app/latest/tool-reference/spatial-analyst/hillshade.htm) ```python # Clip streams to the area of interest sunshine = arcpy.sa.Hillshade('UKy2010_above', 90, 10, "SHADOWS") # 90 degrees is the azimuth, 10 degrees is the altitude ```
In-class task
Color Phase 1 point cloud by amount of removal
and Phase 2 by amount of addition.
Make a layout showing change.
STOP