WEBVTT

00:00:00.000 --> 00:00:48.600


00:00:48.600 --> 00:00:52.300
Lesson 10b3 clipping LAS data and

00:00:52.300 --> 00:00:55.200
creating derivative products in ArcGIS Pro.

00:00:55.200 --> 00:00:57.533
Lessons 10b1 and 10b2 introduced

00:00:57.533 --> 00:01:00.299
the basics of interacting with LiDAR

00:01:00.300 --> 00:01:02.633
point cloud data in ArcGIS Pro.

00:01:02.633 --> 00:01:04.899
This lesson will cover creating and editing

00:01:04.900 --> 00:01:07.000
a shapefile as an area of interest,

00:01:07.000 --> 00:01:10.066
creating the LAS data set and generating

00:01:10.066 --> 00:01:12.232
derivative data in ArcGIS Pro.

00:01:12.233 --> 00:01:14.599
If necessary, please review Lesson 10b1

00:01:14.600 --> 00:01:17.666
for guidance with importing LAS files, filtering

00:01:17.666 --> 00:01:20.932
and visualizing LiDAR point clouds in ArcGIS Pro.

00:01:20.933 --> 00:01:22.499


00:01:22.500 --> 00:01:25.700
Any use of trade, product or firm names is for descriptive

00:01:25.700 --> 00:01:27.900
purposes only and does not imply

00:01:27.900 --> 00:01:30.066
endorsement by the US government.

00:01:30.066 --> 00:01:31.332


00:01:31.333 --> 00:01:32.899
By the end of this lesson,

00:01:32.900 --> 00:01:34.966
you will be able to create and utilize

00:01:34.966 --> 00:01:36.932
a shapefile area of interest to subset

00:01:36.933 --> 00:01:39.566
LAS files, create a LAS data set

00:01:39.566 --> 00:01:41.899
with multiple LAS files,

00:01:41.900 --> 00:01:44.333
generate and export raster services and

00:01:44.333 --> 00:01:47.099
contours derived from the LiDAR point cloud,

00:01:47.100 --> 00:01:50.100
and visualize the data in a 3D scene.

00:01:50.100 --> 00:01:51.700
In this tutorial we are

00:01:51.700 --> 00:01:52.966
using compressed LAS data,

00:01:52.966 --> 00:01:55.199
but the national map also offers

00:01:55.200 --> 00:01:57.166
compressed LAZ data which can

00:01:57.166 --> 00:01:59.432
be converted to LAS files for

00:01:59.433 --> 00:02:00.899
compatibility with ArcGIS Pro.

00:02:00.900 --> 00:02:03.533
Refer to Lesson 10e2 on using

00:02:03.533 --> 00:02:06.366
LASzip to decompress LiDAR data.

00:02:06.366 --> 00:02:07.966
If you're interested in learning

00:02:07.966 --> 00:02:09.966
more about using USGS LiDAR data

00:02:09.966 --> 00:02:11.899
in other software packages.

00:02:11.900 --> 00:02:14.500
Additional videos show how to use

00:02:14.500 --> 00:02:17.900
LiDAR data in Global Mapper and LP360.

00:02:17.900 --> 00:02:20.233
For this lesson we will be using the

00:02:20.233 --> 00:02:22.766
same tiles of LAS data from lesson 10b1.

00:02:22.766 --> 00:02:26.366
If you have not done so already, please download the

00:02:26.366 --> 00:02:31.832
“USGS_LPC_CO_SoPlatteRiver_Data_for_Lessons.zip”

00:02:31.833 --> 00:02:34.499
file from our Rocky website at.

00:02:34.500 --> 00:02:40.300
https://rockyweb.usgs.gov/Training_Data

00:02:40.300 --> 00:02:43.066
After you've downloaded the LiDAR data,

00:02:43.066 --> 00:02:44.899
extract the zip file into a folder

00:02:44.900 --> 00:02:47.400
on your local computer to use during this lesson.

00:02:47.400 --> 00:02:49.666
The national map has a download client

00:02:49.666 --> 00:02:51.466
where you can find USGS products,

00:02:51.466 --> 00:02:55.432
including elevation data such as LiDAR point clouds at

00:02:55.433 --> 00:02:59.199
https://apps.nationalmap.gov/download

00:02:59.200 --> 00:03:01.833
If you are interested in learning more about downloading

00:03:01.833 --> 00:03:04.099
products in the National Map, be sure to check out our

00:03:04.100 --> 00:03:05.833
training videos located at

00:03:05.833 --> 00:03:10.266
www.usgs.gov/NGPvideos.

00:03:10.266 --> 00:03:12.632


00:03:12.633 --> 00:03:16.233
This lesson uses ArcGIS Pro version 2.5.

00:03:16.233 --> 00:03:19.833


00:03:19.833 --> 00:03:22.633
To begin, open ArcGIS Pro on your computer.

00:03:22.633 --> 00:03:25.366


00:03:25.366 --> 00:03:27.299
Launch ArcGIS Pro.

00:03:27.300 --> 00:03:32.900


00:03:32.900 --> 00:03:35.100
Use your enterprise or personal

00:03:35.100 --> 00:03:36.900
login information if prompted

00:03:36.900 --> 00:03:39.366
by the ArcGIS Pro login prompt.

00:03:39.366 --> 00:03:43.699
Open a new map. Name the project "ArcGIS Pro

00:03:43.700 --> 00:03:47.166
LiDAR Training Lesson 10b3".

00:03:47.166 --> 00:03:49.199
Use the same folder we downloaded

00:03:49.200 --> 00:03:51.466
the LiDAR data to in the previous

00:03:51.466 --> 00:03:53.799
Lesson 10b1, then click Okay.

00:03:53.800 --> 00:04:02.133


00:04:02.133 --> 00:04:04.533
Depending on your ArcGIS Pro version,

00:04:04.533 --> 00:04:06.499
the base map may vary.

00:04:06.500 --> 00:04:08.700
This exercise specifically uses

00:04:08.700 --> 00:04:11.133
ArcGIS Pro version 2.5. However,

00:04:11.133 --> 00:04:13.566
any version 2.1 or newer should

00:04:13.566 --> 00:04:15.432
suffice for this tutorial.

00:04:15.433 --> 00:04:17.833
ArcGIS Pro versions older than

00:04:17.833 --> 00:04:20.399
version 2.5 May default to different

00:04:20.400 --> 00:04:22.533
basemaps other than the topographic basemap.

00:04:22.533 --> 00:04:24.366
To change the basemap,

00:04:24.366 --> 00:04:26.932
navigate to the map tab on the layer ribbon

00:04:26.933 --> 00:04:29.333
under the layer section of the map tab,

00:04:29.333 --> 00:04:31.233
click the basemap drop down arrow

00:04:31.233 --> 00:04:34.966
and select the topographic basemap option.

00:04:34.966 --> 00:04:36.999
Under the Map tab, on the layer ribbon,

00:04:37.000 --> 00:04:39.133
select the Add Data button.

00:04:39.133 --> 00:04:42.133
In the add data window, browse the location of the

00:04:42.133 --> 00:04:44.799
LiDAR point cloud LAS files that were downloaded

00:04:44.800 --> 00:04:47.833
from the FTP site. Highlight the six LAS files

00:04:47.833 --> 00:04:50.133
and load them into your map.

00:04:50.133 --> 00:04:51.933


00:04:51.933 --> 00:04:54.299
Under the View Tab, Open the catalog pane.

00:04:54.300 --> 00:04:57.933


00:04:57.933 --> 00:05:00.999
Expand the folders dropdown and right click

00:05:01.000 --> 00:05:04.033
on the project folder and select new.

00:05:04.033 --> 00:05:06.566
Click on shapefile.

00:05:06.566 --> 00:05:08.399
In the geoprocessing pane that

00:05:08.400 --> 00:05:10.266
appears except the default folder

00:05:10.266 --> 00:05:12.032
location under feature class name

00:05:12.033 --> 00:05:15.333
called a Shapefile Lesson10b3_AOI

00:05:15.333 --> 00:05:18.199
for our area of interest.

00:05:18.200 --> 00:05:20.433
Change the geometry type

00:05:20.433 --> 00:05:22.499
to Polygon if necessary.

00:05:22.500 --> 00:05:24.666
Change the coordinate system to

00:05:24.666 --> 00:05:29.232
NAD 1983 UTM zone 13 N.

00:05:29.233 --> 00:05:31.566
You can do this by searching in the

00:05:31.566 --> 00:05:33.066
coordinate system popup window or

00:05:33.066 --> 00:05:34.766
you can select an existing layer

00:05:34.766 --> 00:05:36.866
in the dropdown menu which should

00:05:36.866 --> 00:05:38.632
have the same coordinate system.

00:05:38.633 --> 00:05:41.766
Accept all other defaults.

00:05:41.766 --> 00:05:42.832
Click run.

00:05:42.833 --> 00:05:45.033


00:05:45.033 --> 00:05:46.533
The empty shapefile will be

00:05:46.533 --> 00:05:48.633
added to the table of contents,

00:05:48.633 --> 00:05:50.699
but we still need to draw the

00:05:50.700 --> 00:05:52.600
extent and shape of the Polygon.

00:05:52.600 --> 00:05:55.266
Click on the edit tab and under

00:05:55.266 --> 00:05:57.532
the features ribbon select create.

00:05:57.533 --> 00:05:59.233
The create features pain will

00:05:59.233 --> 00:06:00.633
open and under templates,

00:06:00.633 --> 00:06:05.033
select 10b3 AOI from the list.

00:06:05.033 --> 00:06:07.966
Select the Polygon icon to draw

00:06:07.966 --> 00:06:09.932
the shapefile using vertices

00:06:09.933 --> 00:06:10.033


00:06:10.033 --> 00:06:11.833
to create the Polygon.

00:06:11.833 --> 00:06:13.866
We want our area of interest to be

00:06:13.866 --> 00:06:15.299
everything east of the foothills.

00:06:15.300 --> 00:06:17.666
Using the topographic map and the contour lines

00:06:17.666 --> 00:06:19.566
left click at the center of the top left,

00:06:19.566 --> 00:06:21.032
LAS tile beginning at

00:06:21.033 --> 00:06:22.633
the base of the foothills.

00:06:22.633 --> 00:06:24.666


00:06:24.666 --> 00:06:26.199
Left click again at the bottom

00:06:26.200 --> 00:06:29.866
of the middle tile and trace the tile extends to include the

00:06:29.866 --> 00:06:32.166
LAS tiles east of the foothills.

00:06:32.166 --> 00:06:35.132
Double Click to finish drawing the Polygon.

00:06:35.133 --> 00:06:37.466
To save your edits, under the edit tab,

00:06:37.466 --> 00:06:39.999
select the Save icon in the

00:06:40.000 --> 00:06:42.166
Manage Edit section. Click yes.

00:06:42.166 --> 00:06:45.766
In the dialog window to save all edits.

00:06:45.766 --> 00:06:47.399
In the table of contents pane,

00:06:47.400 --> 00:06:49.566
select the default symbology rectangle

00:06:49.566 --> 00:06:52.899
under Lesson10b3_AOI layer.

00:06:52.900 --> 00:06:55.333
This will bring up the format Polygon

00:06:55.333 --> 00:06:57.666
symbol window in the symbology pane,

00:06:57.666 --> 00:06:59.899
select the black outline,

00:06:59.900 --> 00:07:01.600
2 point thickness.

00:07:01.600 --> 00:07:04.433
Now is a good time to save your project.

00:07:04.433 --> 00:07:07.299
Now that the AOI is clearly defined,

00:07:07.300 --> 00:07:09.600
we can create a last data set

00:07:09.600 --> 00:07:12.133
constrained by the extent of our AOI.

00:07:12.133 --> 00:07:14.966


00:07:14.966 --> 00:07:18.099
A LAS, data set is designed to use LiDAR

00:07:18.100 --> 00:07:21.633
data and .las or .zlas file types.

00:07:21.633 --> 00:07:26.866
It supports any LAS file version 1.0 through 1.4.

00:07:26.866 --> 00:07:30.966
It is important to note that the LAS format supports the

00:07:30.966 --> 00:07:34.799
classification codes for each point according to the

00:07:34.800 --> 00:07:37.466
American Society for Photogrammetry and Remote

00:07:37.466 --> 00:07:41.299
Sensing (ASPRS) and ArcGIS applies the class

00:07:41.300 --> 00:07:44.733
code specified for last version 1.4.

00:07:44.733 --> 00:07:47.666
A LAS, data set has the capability to

00:07:47.666 --> 00:07:50.766
store and reference many LAS files,

00:07:50.766 --> 00:07:53.532
compiling them to one LAS data set

00:07:53.533 --> 00:07:56.899
similar to a mosaic of raster surfaces.

00:07:56.900 --> 00:07:59.533
The LAS data set refers to the LAS

00:07:59.533 --> 00:08:02.499
file or files in their native binary

00:08:02.500 --> 00:08:05.200
format of stored airborne LiDAR data.

00:08:05.200 --> 00:08:06.733
By using a LAS data set,

00:08:06.733 --> 00:08:08.633
users can more efficiently conduct

00:08:08.633 --> 00:08:11.433
QA/QC on LiDAR data and quickly

00:08:11.433 --> 00:08:14.166
visualize many LAS files at once as

00:08:14.166 --> 00:08:16.532
either a tin or point cloud in both

00:08:16.533 --> 00:08:19.233
2D and 3D views in Esri software.

00:08:19.233 --> 00:08:21.566
The LAS data sets calculate

00:08:21.566 --> 00:08:23.599
statistics for referenced LAS files

00:08:23.600 --> 00:08:26.366
in order for the user to understand

00:08:26.366 --> 00:08:29.999
and analyze the LiDAR data. By using a LAS data set,

00:08:30.000 --> 00:08:32.300
we can create raster surfaces

00:08:32.300 --> 00:08:35.100
and derived products to better understand the

00:08:35.100 --> 00:08:38.100
topographic layers from the LiDAR point cloud.

00:08:38.100 --> 00:08:41.166
For more information on what a last data set is,

00:08:41.166 --> 00:08:44.232
or last data set capabilities and properties,

00:08:44.233 --> 00:08:47.633
visit the ArcGIS Pro Help pages listed here.

00:08:47.633 --> 00:08:50.966


00:08:50.966 --> 00:08:54.232
Open the tools to your processing pane in the analysis tab.

00:08:54.233 --> 00:08:56.599


00:08:56.600 --> 00:08:59.300
And search for create LAS data sets.

00:08:59.300 --> 00:09:05.133


00:09:05.133 --> 00:09:08.133
For the input files, browse to the six LAS files

00:09:08.133 --> 00:09:08.233


00:09:08.233 --> 00:09:10.333
downloaded from the Rocky website,

00:09:10.333 --> 00:09:12.433
select all six LAS files.

00:09:12.433 --> 00:09:15.399


00:09:15.400 --> 00:09:17.866
Name the output LAS data set

00:09:17.866 --> 00:09:23.166
"Lesson_10b3_AOI.lasd"

00:09:23.166 --> 00:09:26.966


00:09:26.966 --> 00:09:29.466
Surface constraints can be used to limit

00:09:29.466 --> 00:09:31.866
the LAS data set to a given extent.

00:09:31.866 --> 00:09:33.899
Using our AOI shapefile,

00:09:33.900 --> 00:09:36.766
we will clip the LAS data set to the AOI's extent.

00:09:36.766 --> 00:09:42.932
Under input features select Lesson_10b3_AOI

00:09:42.933 --> 00:09:45.066
from the drop down menu.

00:09:45.066 --> 00:09:46.299
There is no height field,

00:09:46.300 --> 00:09:48.900
so leave it blank.

00:09:48.900 --> 00:09:51.100
For the type, we will use a soft clip

00:09:51.100 --> 00:09:53.633
to define the boundary of the data set.

00:09:53.633 --> 00:10:02.766
Select NAD 1983 UTM Zone 13 N/NAVD88.

00:10:02.766 --> 00:10:04.432
Check the compute statistics if

00:10:04.433 --> 00:10:06.366
not already marked then click run.

00:10:06.366 --> 00:10:09.632


00:10:09.633 --> 00:10:12.133
ArcGIS Pro allows various surface

00:10:12.133 --> 00:10:14.166
constraints depending upon the

00:10:14.166 --> 00:10:16.432
feature class or shapefile type.

00:10:16.433 --> 00:10:18.433
Surface constraints use the features

00:10:18.433 --> 00:10:21.066
geometry to capture or define a surface

00:10:21.066 --> 00:10:23.332
when it is displayed in a TIN format.

00:10:23.333 --> 00:10:26.399
Esri defines options with hard or soft

00:10:26.400 --> 00:10:28.700
designation which refer to whether

00:10:28.700 --> 00:10:30.900
the surface edges represent distinct

00:10:30.900 --> 00:10:33.733
breaks in slope or a gradual change.

00:10:33.733 --> 00:10:35.799
Anchor points can be used for point

00:10:35.800 --> 00:10:38.200
data and hard or soft line constraints

00:10:38.200 --> 00:10:40.833
are typically used for brake lines to

00:10:40.833 --> 00:10:43.066
delineate waterbodies or linear features.

00:10:43.066 --> 00:10:45.399
Hard clips and soft clips are used to

00:10:45.400 --> 00:10:47.833
delineate a boundary of the LAS data set.

00:10:47.833 --> 00:10:49.933
Surface constraints are applied only when

00:10:49.933 --> 00:10:52.699
the LAS data set is drawn as a surface.

00:10:52.700 --> 00:10:55.133
If the LAS data set is drawn using points,

00:10:55.133 --> 00:10:58.099
it will draw all LAS files being

00:10:58.100 --> 00:11:00.633
referenced in the LAS data set.

00:11:00.633 --> 00:11:02.633
For more information on storing LAS data

00:11:02.633 --> 00:11:05.133
sets or LAS files and surface constraints,

00:11:05.133 --> 00:11:07.066
visit the Esri help page shown here.

00:11:07.066 --> 00:11:11.132


00:11:11.133 --> 00:11:12.966
In order to better view the new

00:11:12.966 --> 00:11:15.099
LAS data sets, we will convert

00:11:15.100 --> 00:11:18.166
the map to a local 3D scene.

00:11:18.166 --> 00:11:20.032
In the View tab,

00:11:20.033 --> 00:11:22.766
select convert and click local scene.

00:11:22.766 --> 00:11:25.499
This local scene shows 2D layers draped

00:11:25.500 --> 00:11:28.700
over a 3D scene and draws 3D layers

00:11:28.700 --> 00:11:31.933
that can be extruded or exaggerated to

00:11:31.933 --> 00:11:35.533
show perspective or represent elevation data.

00:11:35.533 --> 00:11:36.899
After you've converted the

00:11:36.900 --> 00:11:38.933
map layout to a local scene,

00:11:38.933 --> 00:11:40.599
save your project by clicking

00:11:40.600 --> 00:11:43.333
on the save icon at the top left

00:11:43.333 --> 00:11:45.366
corner of the Arcgis Pro window.

00:11:45.366 --> 00:11:47.766
In the 3D rendering of the LAS data set,

00:11:47.766 --> 00:11:50.199
make sure the LAS data set layer is

00:11:50.200 --> 00:11:52.100
highlighted and click on the appearance

00:11:52.100 --> 00:11:54.733
tab under the LAS data set Layer Ribbon,

00:11:54.733 --> 00:11:56.866
click on the Symbology icon and in

00:11:56.866 --> 00:11:58.866
the symbology pane click on the Yellow

00:11:58.866 --> 00:12:00.999
Pyramid to draw the layer using a

00:12:01.000 --> 00:12:02.900
surface and choose elevation from the

00:12:02.900 --> 00:12:06.933
dropdown to visualize the data set in the AOI.

00:12:06.933 --> 00:12:10.033
Make sure that the draw using box is checked.

00:12:10.033 --> 00:12:12.233
The LAS data set will only show the

00:12:12.233 --> 00:12:14.366
clipped extent when it is drawn as a

00:12:14.366 --> 00:12:16.632
layer surface rather than a point cloud.

00:12:16.633 --> 00:12:21.066


00:12:21.066 --> 00:12:23.366
You can see the noise and artifacts

00:12:23.366 --> 00:12:25.399
represented in the LAS data set

00:12:25.400 --> 00:12:27.033
because all classes are turned

00:12:27.033 --> 00:12:28.966
on using the elevation field.

00:12:28.966 --> 00:12:31.332
We will change the LAS points

00:12:31.333 --> 00:12:33.533
filter to only ground points.

00:12:33.533 --> 00:12:39.799
Right click on the Lesson10b3_AOI.lasd from the table of

00:12:39.800 --> 00:12:41.966
contents and select properties.

00:12:41.966 --> 00:12:43.799
In the layer properties window,

00:12:43.800 --> 00:12:47.100
select the LAS filter tab. As you can see,

00:12:47.100 --> 00:12:49.000
all the classification codes are turned

00:12:49.000 --> 00:12:51.200
on and currently drawn in the map view.

00:12:51.200 --> 00:12:54.233
Uncheck the all box under the classification

00:12:54.233 --> 00:12:57.133
code to turn all the classes off.

00:12:57.133 --> 00:12:59.133
Then check only the two

00:12:59.133 --> 00:13:01.133
ground class and click Okay.

00:13:01.133 --> 00:13:02.466
Now that we have filtered

00:13:02.466 --> 00:13:03.566
only the ground points,

00:13:03.566 --> 00:13:05.699
we can create a digital elevation

00:13:05.700 --> 00:13:08.500
model or DEM in the analysis tab,

00:13:08.500 --> 00:13:11.233
open the tools geoprocessing window.

00:13:11.233 --> 00:13:14.299
Search for LAS data set to raster.

00:13:14.300 --> 00:13:16.466
Select the LAS data set

00:13:16.466 --> 00:13:18.199
to raster conversion tool.

00:13:18.200 --> 00:13:21.400
The LAS data set conversion to raster

00:13:21.400 --> 00:13:23.700
tool uses interpolation methods to

00:13:23.700 --> 00:13:26.700
generate a smooth output of raster cells.

00:13:26.700 --> 00:13:29.200
The two interpolation types for this

00:13:29.200 --> 00:13:32.000
tool are binning or triangulation.

00:13:32.000 --> 00:13:33.900
Binning is typically faster and

00:13:33.900 --> 00:13:35.766
assigns values based on points

00:13:35.766 --> 00:13:35.866


00:13:35.866 --> 00:13:37.699
within the extent of the cell.

00:13:37.700 --> 00:13:40.000
Binning does not fill voids or data gaps,

00:13:40.000 --> 00:13:42.733
but allows options for filling these gaps.

00:13:42.733 --> 00:13:45.333
Binning gives various cell assignment types,

00:13:45.333 --> 00:13:47.233
such as average.

00:13:47.233 --> 00:13:49.333
Nearest neighbor inverse distance

00:13:49.333 --> 00:13:51.433
weighting and other statistics to

00:13:51.433 --> 00:13:53.866
define the value of the cell based

00:13:53.866 --> 00:13:55.899
on the points been to that cell.

00:13:55.900 --> 00:13:57.800
Triangulation uses interpolation of

00:13:57.800 --> 00:14:00.633
the points to create a triangulated

00:14:00.633 --> 00:14:00.733


00:14:00.733 --> 00:14:02.666
surface of the last data set.

00:14:02.666 --> 00:14:04.566
The interpolation method can be

00:14:04.566 --> 00:14:06.899
either natural neighbor or linear and

00:14:06.900 --> 00:14:09.833
point thinning can be used if needed.

00:14:09.833 --> 00:14:11.733
In the geoprocessing window,

00:14:11.733 --> 00:14:17.433
select the Lesson10b3_AOI.lasd file for the input.

00:14:17.433 --> 00:14:19.666


00:14:19.666 --> 00:14:25.132
Name the output raster DEM_10b3_AOI.

00:14:25.133 --> 00:14:28.266


00:14:28.266 --> 00:14:31.199
We will use triangulation with natural

00:14:31.200 --> 00:14:33.900
neighbor interpolation and no thinning.

00:14:33.900 --> 00:14:36.700
Make sure the data type is floating point.

00:14:36.700 --> 00:14:39.433
And that cell size is selected from

00:14:39.433 --> 00:14:41.733
the sampling type dropdown menu.

00:14:41.733 --> 00:14:43.833
This last data set references QL2

00:14:43.833 --> 00:14:46.699
level last files giving a nominal

00:14:46.700 --> 00:14:49.600
point spacing of .7 meters or better.

00:14:49.600 --> 00:14:51.933
It's a good practice to never use a

00:14:51.933 --> 00:14:53.866
resolution finer than the point spacing.

00:14:53.866 --> 00:14:56.632
We will use a resolution of two meters for

00:14:56.633 --> 00:14:58.999
the DEM in the sampling value section,

00:14:59.000 --> 00:15:02.333
right two then click run.

00:15:02.333 --> 00:15:05.533
This may take a minute or two.

00:15:05.533 --> 00:15:08.033
Once the tool has finished running, save your work.

00:15:08.033 --> 00:15:11.833


00:15:11.833 --> 00:15:14.099
We will now create more derivative

00:15:14.100 --> 00:15:16.966
products from the DEM that was generated.

00:15:16.966 --> 00:15:21.099
Select the DEM_10b3_AOI layer

00:15:21.100 --> 00:15:23.433
from the table of contents.

00:15:23.433 --> 00:15:26.299
Click the analysis ribbon and click

00:15:26.300 --> 00:15:28.200
the raster functions dropdown.

00:15:28.200 --> 00:15:29.933
Select raster functions from the

00:15:29.933 --> 00:15:32.133
dropdown and then expand the surface

00:15:32.133 --> 00:15:34.099
dropdown and select hillshade.

00:15:34.100 --> 00:15:42.100


00:15:42.100 --> 00:15:46.800
Choose the DEM_10b3_AOI layer

00:15:46.800 --> 00:15:46.900


00:15:46.900 --> 00:15:49.300
for the raster input section.

00:15:49.300 --> 00:15:52.500
For this tutorial, we will keep the

00:15:52.500 --> 00:15:55.533
traditional hillshade at 315 degrees

00:15:55.533 --> 00:15:59.699
azimuth and 45 degree altitude.

00:15:59.700 --> 00:16:03.333
Accept the remaining defaults and click create new layer.

00:16:03.333 --> 00:16:06.999


00:16:07.000 --> 00:16:09.466
Reorder the DEM layer in the table of

00:16:09.466 --> 00:16:11.999
contents to be on top of the hill shade.

00:16:12.000 --> 00:16:14.366


00:16:14.366 --> 00:16:16.932
Under appearance, make the DEM 50%

00:16:16.933 --> 00:16:19.099
transparent and choose a color scheme from

00:16:19.100 --> 00:16:21.800
the symbology tab that fits your liking.

00:16:21.800 --> 00:16:27.800


00:16:27.800 --> 00:16:29.533
Now we can visualize the final

00:16:29.533 --> 00:16:31.466
DEM more clearly in the 3D view,

00:16:31.466 --> 00:16:34.332
on top of the hillshade layer.

00:16:34.333 --> 00:16:36.566
Next we will create a digital surface model.

00:16:36.566 --> 00:16:45.966


00:16:45.966 --> 00:16:47.532
A digital surface model

00:16:47.533 --> 00:16:50.266
or DSM is similar to a DEM.

00:16:50.266 --> 00:16:51.966
However, it shows the maximum

00:16:51.966 --> 00:16:54.132
heights of objects above bare ground,

00:16:54.133 --> 00:16:57.899
such as tops of trees and buildings.

00:16:57.900 --> 00:17:00.300
Right click the LAS data set in the

00:17:00.300 --> 00:17:02.966
table of contents and in the layer

00:17:02.966 --> 00:17:05.266
properties under last filter check all,

00:17:05.266 --> 00:17:08.432
then uncheck classes 7 low noise,

00:17:08.433 --> 00:17:12.633
11 roads and 18 high noise.

00:17:12.633 --> 00:17:14.299
These classes contain some

00:17:14.300 --> 00:17:16.433
misclassifications in the high noise

00:17:16.433 --> 00:17:18.599
class and will show artifacts and

00:17:18.600 --> 00:17:21.100
unwanted noise points in our DSM results.

00:17:21.100 --> 00:17:27.533


00:17:27.533 --> 00:17:29.099
Once the points are filtered,

00:17:29.100 --> 00:17:31.166
open the geoprocessing pane and

00:17:31.166 --> 00:17:33.832
search last data set to raster.

00:17:33.833 --> 00:17:35.566
In this case, the tool is still

00:17:35.566 --> 00:17:37.532
open from the last time it was run.

00:17:37.533 --> 00:17:41.699
Select the tutorial 10b3_AOI.lasd

00:17:41.700 --> 00:17:44.300
file from the input dropdown.

00:17:44.300 --> 00:17:50.666
Name the output raster DSM_10b3_AOI.

00:17:50.666 --> 00:17:52.932
Select elevation for the value

00:17:52.933 --> 00:17:55.566
field and choose binning maximum and

00:17:55.566 --> 00:17:57.799
linear for the interpolation fields.

00:17:57.800 --> 00:18:01.733


00:18:01.733 --> 00:18:03.933
Make sure the sampling type is cell

00:18:03.933 --> 00:18:06.399
size and type 2 for the sampling

00:18:06.400 --> 00:18:08.100
value or resolution. Click run.

00:18:08.100 --> 00:18:14.266


00:18:14.266 --> 00:18:15.732
We'll create a hillshade

00:18:15.733 --> 00:18:17.966
for the DSM using the same

00:18:17.966 --> 00:18:20.166
methods for the DEM hillshade.

00:18:20.166 --> 00:18:22.532
Once you have your DSM raster created,

00:18:22.533 --> 00:18:23.866
highlight the layer.

00:18:23.866 --> 00:18:25.666
Click the analysis ribbon

00:18:25.666 --> 00:18:27.766
and choose raster functions.

00:18:27.766 --> 00:18:31.599


00:18:31.600 --> 00:18:33.133
Expand the surface dropdown

00:18:33.133 --> 00:18:34.266
and choose hillshade.

00:18:34.266 --> 00:18:38.766


00:18:38.766 --> 00:18:44.266
Select DSM_10b3_AOI from the input dropdown

00:18:44.266 --> 00:18:49.299
and choose create new layer.

00:18:49.300 --> 00:18:49.900


00:18:49.900 --> 00:18:52.200
You may edit the symbology and transparency

00:18:52.200 --> 00:18:54.633
to view the DSM over the hillshade.

00:18:54.633 --> 00:18:58.433


00:18:58.433 --> 00:19:01.033
You can also set your DSM to

00:19:01.033 --> 00:19:02.766
your preferred color scheme.

00:19:02.766 --> 00:19:04.699
By comparing the hillshade layers

00:19:04.700 --> 00:19:07.266
created from the DEM and the DSM,

00:19:07.266 --> 00:19:09.132
the surface features become very evident.

00:19:09.133 --> 00:19:13.266


00:19:13.266 --> 00:19:15.366
Turn off the DEM in the DSM so

00:19:15.366 --> 00:19:17.099
that only the hillshade layers

00:19:17.100 --> 00:19:19.400
are visible in the map window.

00:19:19.400 --> 00:19:21.666
Order the hillshade layers so

00:19:21.666 --> 00:19:24.432
that the DSM hillshade layer is

00:19:24.433 --> 00:19:26.533
above the DM hillshade layer.

00:19:26.533 --> 00:19:28.899
The swipe tool only works in a 2D map,

00:19:28.900 --> 00:19:31.566
so we will need to convert the

00:19:31.566 --> 00:19:33.999
3D scene back into a 2D map.

00:19:34.000 --> 00:19:36.100
Click the View tab and select

00:19:36.100 --> 00:19:37.533
the convert dropdown arrow.

00:19:37.533 --> 00:19:38.699
Choose to map.

00:19:38.700 --> 00:19:41.433
Once the 2D map has been generated,

00:19:41.433 --> 00:19:43.499
make sure the hillshade

00:19:43.500 --> 00:19:46.866
DSM_10b3 layer is on top of

00:19:46.866 --> 00:19:49.932
the hillshade DEM_10b3 layer.

00:19:49.933 --> 00:19:52.133
Zoom into an area that seems densely

00:19:52.133 --> 00:19:53.999
covered in trees and buildings,

00:19:54.000 --> 00:19:56.600
with the hillshade DSM_10b3 layer

00:19:56.600 --> 00:19:58.700
selected in the table of contents.

00:19:58.700 --> 00:20:00.333
Select the Appearance tab

00:20:00.333 --> 00:20:02.766
and click on the swipe tool.

00:20:02.766 --> 00:20:04.299
Click and drag the swipe tool

00:20:04.300 --> 00:20:06.666
up or down or side to side to

00:20:06.666 --> 00:20:08.232
visualize the differences in the

00:20:08.233 --> 00:20:10.533
DSM and DEM hillshade layers.

00:20:10.533 --> 00:20:13.766


00:20:13.766 --> 00:20:15.466
When finished exploring the differences,

00:20:15.466 --> 00:20:17.499
close the 2D map window and

00:20:17.500 --> 00:20:19.433
return to the 3D scene.

00:20:19.433 --> 00:20:21.333
It's important to note that the

00:20:21.333 --> 00:20:23.333
hillshade layers were not saved as

00:20:23.333 --> 00:20:25.399
individual rasters to the project folder.

00:20:25.400 --> 00:20:28.333
Rather, they are currently referencing the

00:20:28.333 --> 00:20:31.433
source layers as either the DSM or DEM.

00:20:31.433 --> 00:20:34.066
We will need to export these layers to

00:20:34.066 --> 00:20:36.966
save them as individual hillshade rasters.

00:20:36.966 --> 00:20:40.899
Right click the hillshade DSM_10b3 raster.

00:20:40.900 --> 00:20:44.633
Select data and then export raster.

00:20:44.633 --> 00:20:47.199
Except the default name and location,

00:20:47.200 --> 00:20:48.800
the cell size should be the same

00:20:48.800 --> 00:20:49.900
as the source raster.

00:20:49.900 --> 00:20:52.133
2 meters by two meters.

00:20:52.133 --> 00:20:53.499
Accept the remaining

00:20:53.500 --> 00:20:55.300
defaults and click export.

00:20:55.300 --> 00:21:01.400


00:21:01.400 --> 00:21:03.533
Repeat the same process for

00:21:03.533 --> 00:21:05.233
the hillshade DEM_10b3 layer.

00:21:05.233 --> 00:21:11.433


00:21:11.433 --> 00:21:13.266
Now is a good time to save your project.

00:21:13.266 --> 00:21:15.632


00:21:15.633 --> 00:21:17.499
Now that we have created elevation

00:21:17.500 --> 00:21:19.500
models for both the bare Earth

00:21:19.500 --> 00:21:21.266
ground and the surface features,

00:21:21.266 --> 00:21:23.432
we can use map algebra to determine

00:21:23.433 --> 00:21:25.099
the actual heights of various

00:21:25.100 --> 00:21:27.166
features across the area of interest.

00:21:27.166 --> 00:21:29.366
A digital height model or DHM

00:21:29.366 --> 00:21:31.766
can be created by subtracting the

00:21:31.766 --> 00:21:34.832
bare earth raster DEM from the DSM,

00:21:34.833 --> 00:21:36.199
showing the true heights of

00:21:36.200 --> 00:21:37.566
objects relative to the ground.

00:21:37.566 --> 00:21:40.866


00:21:40.866 --> 00:21:43.199
Open the analysis tab and click

00:21:43.200 --> 00:21:45.200
on the tools Geoprocessing icon.

00:21:45.200 --> 00:21:49.133
Search for raster calculator and select the

00:21:49.133 --> 00:21:51.966
raster calculator spatial analyst tool.

00:21:51.966 --> 00:21:59.299
Double click DSM_10b3_AOI. The minus sign.

00:21:59.300 --> 00:22:03.800
And DEM_10b3_AOI.

00:22:03.800 --> 00:22:07.100
Or you can type the following expression

00:22:07.100 --> 00:22:20.833
"DSM_10b3_AOI" - "DEM_10b3_AOI"

00:22:20.833 --> 00:22:23.533
in the map algebra expression window.

00:22:23.533 --> 00:22:29.666
Name the output raster DHM_10b3_AOI. Click run.

00:22:29.666 --> 00:22:33.032


00:22:33.033 --> 00:22:36.599
Turn off all layers except the DHM_10B3 raster.

00:22:36.600 --> 00:22:38.500
Changed the symbology of the

00:22:38.500 --> 00:22:40.800
DHM_10b3 raster layer to

00:22:40.800 --> 00:22:43.300
the Inferno color scheme in the appearance tab.

00:22:43.300 --> 00:22:51.233


00:22:51.233 --> 00:22:54.033
Under the Map tab, select the Explorer icon.

00:22:54.033 --> 00:22:56.066


00:22:56.066 --> 00:22:58.166
Navigate to the circular buildings

00:22:58.166 --> 00:23:00.866
in the north edge of the AOI.

00:23:00.866 --> 00:23:02.966
By clicking on the brightest part of

00:23:02.966 --> 00:23:05.332
the building a popup will show up

00:23:05.333 --> 00:23:07.099
containing the pixels height value.

00:23:07.100 --> 00:23:09.266
You can determine heights of any object

00:23:09.266 --> 00:23:11.732
by clicking on the pixel in the DHM_10b3

00:23:11.733 --> 00:23:15.533
layer and reading the attributes in the popup window.

00:23:15.533 --> 00:23:19.266
This building is approximately 40 meters tall.

00:23:19.266 --> 00:23:21.466
ArcGIS Pro has numerous pre

00:23:21.466 --> 00:23:23.032
selected surface derivatives that

00:23:23.033 --> 00:23:25.133
are saved under the data tab in

00:23:25.133 --> 00:23:27.066
the last data set layer ribbon.

00:23:27.066 --> 00:23:30.399
Deselect the DHM and turn on the

00:23:30.400 --> 00:23:35.333
Lesson_10b3_AOI.lasd layer in the table of contents.

00:23:35.333 --> 00:23:37.566
Then select the data tab from

00:23:37.566 --> 00:23:40.166
the last data set layer ribbon.

00:23:40.166 --> 00:23:42.132
Click the surface derivatives

00:23:42.133 --> 00:23:44.633
dropdown arrow and choose contours.

00:23:44.633 --> 00:23:46.366
Before running the contour tool,

00:23:46.366 --> 00:23:48.166
we must filter out all LiDAR

00:23:48.166 --> 00:23:50.032
points that are not the ground

00:23:50.033 --> 00:23:54.633
with the Lesson_10b3_AOI.lasd layer

00:23:54.633 --> 00:23:56.866
selected in the table of contents.

00:23:56.866 --> 00:23:58.432
Select the appearance tab from

00:23:58.433 --> 00:24:00.633
the LAS data set Layer ribbon.

00:24:00.633 --> 00:24:02.933
Click the LAS filter icon and filter

00:24:02.933 --> 00:24:05.433
only the two ground classification.

00:24:05.433 --> 00:24:08.666


00:24:08.666 --> 00:24:10.499
Go back to the geoprocessing

00:24:10.500 --> 00:24:12.766
window and name the output feature

00:24:12.766 --> 00:24:15.399
class Contours_10b3.

00:24:15.400 --> 00:24:19.033


00:24:19.033 --> 00:24:22.633
Enter 10 for the contour interval.

00:24:22.633 --> 00:24:24.566
The map units are in meters,

00:24:24.566 --> 00:24:26.632
so we will have a detailed

00:24:26.633 --> 00:24:28.566
contour feature class for every

00:24:28.566 --> 00:24:30.766
10 meters of change in elevation.

00:24:30.766 --> 00:24:33.699
Enter 50 for the index interval.

00:24:33.700 --> 00:24:35.966
This will show an index contour

00:24:35.966 --> 00:24:38.066
for every 50 meters. Click run.

00:24:38.066 --> 00:24:39.799


00:24:39.800 --> 00:24:42.000
The Contour feature class we created

00:24:42.000 --> 00:24:44.533
is stored in the projects geodatabase.

00:24:44.533 --> 00:24:47.099
We will alter the symbology to view the

00:24:47.100 --> 00:24:49.200
contour lines and the index contours.

00:24:49.200 --> 00:24:53.333
Turn off the Lesson_10b3_AOI.lasd layer.

00:24:53.333 --> 00:24:56.299
Right click the Contour_10B3 feature class in

00:24:56.300 --> 00:25:00.800
the table of contents and select symbology. Under the

00:25:00.800 --> 00:25:03.700
primary symbology header. Choose unique values from

00:25:03.700 --> 00:25:07.533
the drop down. Change field one to index contours.

00:25:07.533 --> 00:25:09.299
The index contour function in

00:25:09.300 --> 00:25:11.533
the contour tool created a binary

00:25:11.533 --> 00:25:13.599
output for each contour feature.

00:25:13.600 --> 00:25:15.833
One means that it was an index

00:25:15.833 --> 00:25:17.966
contour and 0 means it's not.

00:25:17.966 --> 00:25:21.766
Click on the line for the 0 attribute.

00:25:21.766 --> 00:25:22.432
If necessary,

00:25:22.433 --> 00:25:24.466
click the properties tab to adjust

00:25:24.466 --> 00:25:26.399
the line color and linewidth.

00:25:26.400 --> 00:25:28.700
Choose a dark color for the one

00:25:28.700 --> 00:25:30.400
index contour and a lighter

00:25:30.400 --> 00:25:32.400
shade for the 0 contour lines.

00:25:32.400 --> 00:25:35.033
For this tutorial we will use the

00:25:35.033 --> 00:25:37.533
medium coral light with the line width

00:25:37.533 --> 00:25:40.133
of one for the zero value contour.

00:25:40.133 --> 00:25:42.799
And Tuscan red with a line width of two

00:25:42.800 --> 00:25:45.800
for the one value index contour types.

00:25:45.800 --> 00:25:50.833


00:25:50.833 --> 00:25:52.966
Now we can clearly distinguish the index

00:25:52.966 --> 00:25:54.966
contours from the other contour lines.

00:25:54.966 --> 00:26:00.999


00:26:01.000 --> 00:26:03.866
Save your project in close ArcGIS Pro.

00:26:03.866 --> 00:26:09.699


00:26:09.700 --> 00:26:12.700
Congratulations, you finished Lesson 10b3

00:26:12.700 --> 00:26:15.633
clipping LAS data and creating

00:26:15.633 --> 00:26:18.166
derivative products in ArcGIS Pro.

00:26:18.166 --> 00:26:19.599
In this lesson, we discussed

00:26:19.600 --> 00:26:21.400
how to create and utilize a

00:26:21.400 --> 00:26:23.000
shapefile area of interest subset.

00:26:23.000 --> 00:26:25.233
LAS files, create a LAS data set

00:26:25.233 --> 00:26:27.533
with multiple LAS files, generate and

00:26:27.533 --> 00:26:29.599
export raster services and contours

00:26:29.600 --> 00:26:31.733
derived from the LiDAR point cloud,

00:26:31.733 --> 00:26:34.199
and visualize the data in a 3D scene.

00:26:34.200 --> 00:26:36.066
If you are interested in learning

00:26:36.066 --> 00:26:37.932
more about using LiDAR data,

00:26:37.933 --> 00:26:40.299
please see the additional USGS

00:26:40.300 --> 00:26:42.933
LiDAR training videos in ArcGIS Pro,

00:26:42.933 --> 00:26:45.399
LP360, and Global Mapper at

00:26:45.400 --> 00:26:49.800
www.usgs.gov/NGPvideos