Exercise: Use Satellite Data To Map A Land Cover Type
Exercise: Use Satellite Data To Map A Land Cover Type
Exercise
Use satellite data to map a land
cover type
Section 2 Exercise 2
August 17, 2021
Imagery in Action
Time to complete
40 minutes
Introduction
Imagery layers in ArcGIS Pro can be used from online and file-based sources for analysis and
displayed at the same time in the same map. Depending on the source of the imagery layer,
the functionality and visualization options can be different. For online layers, the provider of
the imagery layer controls many aspects of the imagery layer, which affects the user's ability to
download, visualize, or use it in analysis. ArcGIS Living Atlas of the World imagery layers are
provided by Esri for users to support their studies and form the foundation for shared
knowledge about the earth. ArcGIS Online also provides additional information that other
organizations have shared to inform and further additional studies. Your organization's data
combined with some of these additional layers can create new insights and allow for better-
informed decision-making. Through ArcGIS Image for ArcGIS Online, your organization can
provide imagery layers that are specific for your organization and can be used in imagery
analysis.
For more information about ArcGIS Image for ArcGIS Online, see the ArcGIS Image for
ArcGIS Online product page (https://links.esri.com/ArcGISImageOnline | www.esri.com/en-us/
arcgis/products/arcgis-online-imagery/overview).
Exercise scenario
This exercise is a continuation of the previous exercise in Section 2. After reviewing your
vegetation mask based on Landsat imagery, you have decided to refine the vegetation layer
with higher-resolution imagery to ensure that the calculations are as good as they can be. The
Sentinel-2 satellite data offers higher-resolution data with similar multispectral characteristics.
A Sentinel-2 raster has been acquired, and your organization has published a multispectral
imagery layer for you to use for your vegetation area analysis. In this exercise, you will use an
ArcGIS Online imagery layer to refine your vegetation mask in ArcGIS Pro to provide more
accurate measurements of vegetation area.
b In the Catalog pane, click the Project tab, expand the Maps section, and then right-click
the current map and choose Rename.
f Click the My Organization button , and then in the search field, type sentinel and press
Enter.
g Right-click the Paradise Valley Multispectral imagery layer , point to Add To New, and
then choose Map.
The Paradise Valley Multispectral imagery layer is based on a Sentinel-2 raster dataset. When
the imagery layer was created in ArcGIS Image for ArcGIS Online, the multispectral 10
template was used to create the imagery layer, which controls which bands are added to the
imagery layer and the default visualization. Each configuration option during the publishing
affects the bands that are used for the imagery layer, which can affect the spatial resolution of
the imagery layer.
Note: ArcGIS Image for ArcGIS Online will be covered in more detail later in the MOOC. The
configuration options are similar to the options available for mosaic datasets in ArcGIS Pro.
h In the Contents pane, open the Paradise Valley Multispectral imagery layer properties.
_______________________________________________________________________________
_______________________________________________________________________________
When choosing a band combination for your analysis, consider the spatial resolution of the
raster bands based on the desired feature. For vegetation analysis, you will use the Near
Infrared and Red raster bands, which have a spatial resolution of 10 meters. So, the band
combination for the imagery layer that you will use is considered to be finer and have a higher
spatial resolution than the Landsat imagery layer used before.
The dynamic Sentinel imagery that you have added to the map has been created for the
vegetation analysis using only the high-resolution bands (Red, Blue, Green, Near Infrared). To
maintain the higher spatial resolution, only the raster bands that have the higher spatial
resolution have been added. If another study requires the shortwave-infrared layers, then a
new dynamic imagery layer will need to be created with those bands. However, this imagery
layer will have a spatial resolution of 20 meters to reflect the addition of the new bands.
a In the Catalog pane, click the ArcGIS Online button , and then in the Search ArcGIS
Online field, type potential water owner:esritrainingsvc and press Enter.
b Right-click the Potential Area For Water Policy feature layer and choose Add To
Current Map.
c Click the My Content button , and then locate and add your NDVI Mask For Water
Policy imagery layer to the map.
Note: This exercise uses the imagery layer that was created in a previous exercise. If you did
not complete the Section 1 exercise, then you can add the NDVI Mask For Water Policy
Solution layer provided by clicking the ArcGIS Online button and searching for ndvi mask
owner:esritrainingsvc.
The additional layers overlay the Sentinel imagery layer.
For this imagery layer, you can modify the band combination in the Symbology pane or with
the Band Combination button in the Image Service Layer options.
e In the Contents pane, open the Paradise Valley Multispectral layer properties.
f From the Processing Templates tab, click the Processing Template down arrow to see the
options.
3. What are the processing templates listed for the imagery layer?
_______________________________________________________________________________
The imagery layer from ArcGIS Online is a different type of layer than the layer provided by
ArcGIS Living Atlas. Due to the size and scale of coverage for the Sentinel-2 Views layer, the
imagery layer cannot be manipulated as a normal raster dataset added to ArcGIS Pro would.
However, the imagery layer from ArcGIS Image for ArcGIS Online is a different type of layer
and can be manipulated like a normal raster dataset. When the imagery layers are created
with ArcGIS Image for ArcGIS Online, they can be used for geoprocessing operations and
raster functions.
h In the Contents pane, turn off the visibility of the NDVI Mask layer.
i Select the Paradise Valley Multispectral imagery layer, and then click the Appearance tab.
j In the Rendering group, click Band Combination and choose Color Infrared.
The imagery layer was created with the Sentinel-2 raster type, allowing ArcGIS Pro to use the
metadata for the imagery layer to create band combinations based on the band metadata.
The visualization of the imagery layer uses the Color Infrared band combination but is very
bright. To fix the brightness, you will use the Dynamic Range Adjustment button.
Dynamic range adjustment will use only statistics calculated from the pixels displayed on your
screen for the display. For darker visualizations, this will brighten the display, and for brighter
visualizations, it will reduce the brightness.
k In the Contents pane, ensure that the Paradise Valley Multispectral layer is selected.
The color-infrared visualization with DRA allows vegetation to be more visible, and the
contrast within the imagery layer is more vivid.
Processing Benefits
method
Geoprocessing
tools • Creates new data on disk
• Easily view and edit geoprocessing history
• Create and use custom geoprocessing tools
• Use batch geoprocessing
• String together processes for complex modeling
• Save sets of geoprocessing tools in your project
• Can be run in the Enterprise environment to use distributed
processing
Raster functions
• No new dataset created
• Fast and efficient; on the fly
• Easily view and edit raster function history
• Create and use custom raster functions
• Chain together processes for complex modeling
• Save sets of raster functions in your project
• Generate processing templates for image services
• Can be run in the Enterprise environment to use distributed
processing
a In the Contents pane, ensure that the Paradise Valley Multispectral imagery layer is still
selected.
b On the Imagery tab, in the Tools group, click Indices and choose NDVI.
The raster bands referred to in the NDVI window are based on the metadata for the imagery
layer. For this imagery layer, there are four bands in the imagery layer based on the Sentinel-2
raster, which had 11 bands.
d Click OK.
_______________________________________________________________________________
The values reported in the layer indicate the measure of the Vegetation Index similar to the
calculation in the Landsat Viewer app.
_______________________________________________________________________________
_______________________________________________________________________________
Note: The output specifies the calculation used in the raster function.
This information indicates that both the Near Infrared and Red bands are used in the
calculation as defined by the NDVI method. The Band Arithmetic function is capable of using
known calculations, such as NDVI and custom calculations input by users.
For more information about the Band Arithmetic function, see ArcGIS Pro Help: Band
Arithmetic function.
g Click Cancel to close the Band Arithmetic Properties dialog box, and then close the raster
function template editor.
a On the Imagery tab, in the Analysis group, click Raster Functions to open the Raster
Functions pane.
The raster layer has been clipped to the extent of the study area.
7. In the Contents pane, what are the values reported for the clipped NDVI layer?
_______________________________________________________________________________
The values reported in the Contents pane are based on the values in the layer. For the larger
NDVI layer, those values are in the visible area in the map. Now that the spatial extent of the
layer has changed, the values in the Contents pane will change. The new raster layer based
on the study area contains a different subset of values when compared to the entire dataset.
b Click the Greater Than function, and then in the Greater Than Properties, specify the
following parameters:
You are specifying a value that will be used to generate a constant raster where each cell
holds that value. The Greater Than raster function will compare the clipped NDVI raster to
that value and in the output, it will assign a value of 1 to each cell with a value greater than
0.45 and a value of 0 to each cell with a value less than 0.45.
Note: The value of 0.45 was previously established as the threshold value to evaluate the
NDVI results in the Landsat Explorer app.
The output is visualized with the Stretch renderer. The values in the output are either 0 or 1 to
indicate which areas are above the threshold of 0.45. To create a mask that compares to the
previous vegetation mask, you will change the symbology of the raster layer.
e On the Appearance tab, in the Rendering group, click Symbology to open the Symbology
pane.
f In the Symbology pane, change the Primary Symbology from Stretch to Classify.
The classified renderer will list two classes. One class will have a value of 1 and the other a
value of 0. These classes indicate values that exceed the threshold and values that do not.
h In the lower section of the Symbology pane, on the Classes tab, click the color swatch for
the 0 value and choose No Color.
The vegetation mask layer will display as green and the areas that are not part of the mask will
not be displayed.
a In the Contents pane, turn off the visibility of the Clip_NDVI_Paradise Valley layer and turn
on the visibility of your NDVI Mask For Water Policy layer.
You can see the blue of the NDVI mask around some of the edges of the new mask.
Due to the differences in the spatial resolution between the Landsat imagery and the Sentinel
imagery, the new vegetation mask appears to cover a smaller area than the previous mask. To
refine the estimate of the area, you will need to survey the total area covered by the new
mask. To see the total area covered, you will need to generate a raster attribute table to see
the total area covered by the Greater Than raster layer. The output of the raster function does
not have a raster attribute table, so you will export the mask to create a raster attribute table.
e In the Export Raster pane, for Output Raster Dataset, click the Browse button .
h From the Catalog pane, click the Project tab and then under Folders, right-click Data and
choose Refresh to see the NDVI_RefinedMask.tif file.
Now that you have created the new raster dataset, you are ready to create the raster attribute
table to measure the area. Before processing, it is good practice to remove the raster to
prevent the table from locking and not updating.
k On the Analysis tab, in the Geoprocessing group, click Tools to open the
Geoprocessing pane.
l In the Geoprocessing pane, search for build raster and choose Build Raster Attribute
Table.
m In the Build Raster Attribute Table tool, for the Input Raster, choose
NDVI_RefinedMask.tif.
Note: If you do not see this file in the drop-down list, click the Browse button , browse to
C:\EsriMOOC\Data, select NDVI_RefinedMask.tif, and click OK.
n Click Run.
Note: The symbology colors in your map may differ from the preceding graphic.
The raster layer symbology changes due to the presence of the raster attribute table. The
raster layer will be visualized using the Unique Values renderer.
o In the Symbology pane, change the colors of the values as you did before, making the 0
value display as No Color and the 1 value display as green again.
p In the Contents pane, turn off the visibility of the Greater Than_Clip_NDVI layer.
_______________________________________________________________________________
The value reported in the raster attribute table indicates the number of cells that have a value
of 1. To calculate the area, the cell size multiplied by the count indicates the total area.
Note: Remember that the cell size is 10 meters, which indicates that the area covered by one
cell is 10 m x 10 m or 100 m.
9. What is the refined calculated vegetation area based on the raster attribute table?
_______________________________________________________________________________
_______________________________________________________________________________
In the Landsat-based area calculation, the vegetation mask estimated an area of 4.8 sq.
kilometers for comparison.
s In the Contents pane, turn off the visibility of the Paradise Valley Multispectral imagery
layer and drag your NDVI Mask For Water Policy layer above the new
NDVI_RefinedMask.tif layer.
In the map, you can see the additional areas visualized in green polygons in the new
vegetation mask. These new areas found explain how the estimated vegetation area has
changed.
Because the new vegetation mask is based on higher-resolution imagery, you were able to
identify smaller vegetation areas not included in the Landsat vegetation mask. When
presenting this data to your manager, you now have a good estimate of the overall total
vegetation area and may decide to further refine the vegetation area calculation to remove
some of the new smaller areas discovered in favor of larger areas. Either way, you have
created information for decision makers to consider when making any water policy changes.
3. What are the processing templates listed for the imagery layer?
There are no processing templates listed in the drop-down list.
4. In the Contents pane, what are the values reported in the NDVI_Paradise Valley
Multispectral layer?
The values are 1 to -0.4229.
5. In the raster function template editor, which raster function is used to create the
visualization?
The Band Arithmetic function is used to calculate the NDVI.
7. In the Contents pane, what are the values reported for the clipped NDVI layer?
The values reported now are 0.747 to -0.235.
9. What is the refined calculated vegetation area based on the raster attribute table?
75,000 x 100 sq. meters (cell size 10 x 10) = 7,500,000 sq. meters or 7.5 sq
kilometers. (Note that your calculation will vary slightly from this example.)