Orthorectification Module
Orthorectification Module
ENVIOrthorectification
Module
Orthorectification Module
Files Used in This Tutorial
Background
The ENVI Classic Orthorectification Module Workflow
Installing the ENVI Classic Orthorectification Module
Selecting Input Images and DEMs
Loading the Images and DEM into ENVI Classic
Loading the Images and DEM into the ENVI Classic Orthorectification Wizard
Working with Ground Control Points
Restoring GCPs
Adding New GCPs
Evaluating Residual Errors
Ordering Images and Defining Cutlines
Working with Tie Points
Defining Output Parameters
Optional Step
2
2
2
3
4
5
5
7
9
9
11
12
15
17
20
20
Page 1 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
Orthorectification Module
This tutorial demonstrates how to create an orthorectified mosaic from two QuickBird multispectral
images and a DEM of Phoenix, Arizona, USA. You will learn how to properly open the images and
DEM in ENVI Classic, restore a ground control point (GCP) file, add new GCPs, add tie points, and
define the parameters of the orthorectified product.
You must have an ENVI Classic Orthorectification Module license to run this tutorial.
Description
QuickBird Level-1 multispectral imagery for
Phoenix, AZ from 11 July 2005
QuickBird Level-1 multispectral imagery for
Phoenix, AZ from 09 October 2005
DEM subset of Phoenix, AZ, in GeoTIFF format
Ground control points (GCPs) for the 05JUL*
QuickBird image (in the 005606990010_01_P008_
MUL directory)
Note: It is highly recommended that you save a copy of the entire rigorous_ortho directory to
your local drive, to speed up performance and to simplify editing the GCP file (see "Restoring GCPs"
on page 9).
QuickBird files are courtesy of DigitalGlobe and may not be reproduced without explicit permission
from DigitalGlobe.
Background
An orthorectified image (or orthophoto) is one where each pixel represents a true ground location and all
geometric, terrain, and sensor distortions have been removed to within a specified accuracy. Scale is
constant throughout the orthophoto, regardless of elevation, thus providing accurate measurements of
distance and direction. Geospatial professionals can easily combine orthophotos with other spatial data
in a geographic information system (GIS) for city planning, resource management, and other related
fields.
ENVI Classic's existing orthorectification tools (accessed through Map > Orthorectification in the
ENVI Classic main menu bar) allow you to orthorectify images using rational polynomial coefficients
(RPCs), elevation and geoid information, and optional ground control points (GCPs). However, RPCs
and elevation information do not provide enough details to build a rigorous model representing the path of
light rays from a ground object to the sensor.
Page 2 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
The ENVI Classic Orthorectification Module allows you to build highly accurate orthorectified images
by rigorously modeling the object-to-image transformation. The details of this transformation are mostly
transparent to the user, which means you can quickly create orthorectified images without defining
detailed model parameters.
With the ENVI Classic Orthorectification Layout Manager, you can easily see the spatial coverage of
images, DEMs, GCPs, and tie points, along with residual error vectors for each GCP. You can adjust
your GCPs and tie points to improve the overall residual error for the orthorectified output. ENVI
Classic uses the latest innovative research methods to perform color balancing and mosaicking in the
output product.
The ENVI Classic Orthorectification Module is a joint collaboration between Exelis VIS and
Spacemetric AB of Stockholm, Sweden. Spacemetric designed the underlying block-adjustment model,
which provides a precision orthorectification solution for various sensors. For more technical information
on Spacemetrics orthorectification models, see their website at http://www.spacemetric.se/technology.
Page 3 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
Page 4 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
Page 5 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
Page 6 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
2. Use the Ctrl key to select both GeoTIFF images (05OCT* and 05JUL*). Click OK. The images
appear in the Input Images list in the Wizard.
3. Click the
button next to the Input DEMs section. The Select Input Files dialog appears.
4. Select Band 1 under phoenix_DEM_subset.tif, and click OK. The DEM filename
appears in the Input DEMs list in the Wizard.
5. Click Next in the Wizard. After the images and DEMs load, their boundaries are shown in the
Layout Manager. Depending on your system, this may take some time.
The DEM (shown with a red dashed box) surrounds the two QuickBird images. Use the DEM
checkbox in the Layout Manager to control the display of the DEM boundary in the Draw Area.
Notice also that the two QuickBird images overlap. The image boundaries are shown with solid
green boxes. Use the Image checkbox in the Layout Manager to control the display of the image
boundaries in the Draw Area.
Page 7 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
The yellow dashed box represents the default output area for the orthorectified product, which is
currently selected to be the image boundaries. Use the Output Area checkbox in the Layout
Manager to control the display of the output area in the Draw Area.
The names of the input images and DEM appear, respectively, under the Images and DEMs tabs
of the Layout Manager:
Page 8 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
Restoring GCPs
For this tutorial, you will restore a set of GCPs that lie within the 05JUL* QuickBird image boundary
(the left image as shown in the Layout Manager). DigitalGlobe, who provided the QuickBird images,
Page 9 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
also produced a set of GCPs for the entire Phoenix area, in Microsoft Excel format. However, the ENVI
Classic Orthorectification Module requires GCPs to be in a Rigorous Orthorectification GCP format.
For more information, see the ENVI Classic Help.
The GCP file you are about to restore (phoenixGCPs.pts) was created by identifying those GCPs
that fell within the 05JUL* QuickBird image boundary, then entering them one-by-one into the ENVI
Classic Orthorectification Wizard using map coordinates provided by DigitalGlobe. The GCPs were
already in a Geographic Lat/Lon projection with WGS-84 datum (the same projection as the QuickBird
images and DEM). The GCP file was then saved to ENVI Classic format, using the Save GCPs button
in the Wizard.
The first few lines in the GCP file (phoenixGCPs.pts) are the header lines, and they are preceded
by semicolons. They list the input images you selected earlier. They currently only list the filenames, but
the Rigorous Orthorectification GCP file format requires a full path to the image filenames. You will
need to edit phoenixGCPs.pts to include the full path to the 05JUL* QuickBird image filename on
your computer. As mentioned in the introduction, it is highly recommended that you save a copy of the
entire rigorous_ortho directory to your local drive, so that you can save your changes to the GCP
file.
1. Navigate to rigorous_ortho and open the file phoenixGCPs.pts in a text editor.
2. Edit the third line that begins with FileName1 to include the full path to the 05JUL* image.
Following is an example of editing the line in Windows, assuming you copied the entire
rigorous_ortho directory to your C: drive. The bold text indicates the part that you add:
; FileName1=C:\rigorous_ortho\005606990010_01_P008_MUL\05JUL11182931M1BS-005606990010_01_P008.TIF
Following is an example of editing this line in Linux, assuming you copied the entire
rigorous_ortho directory to a directory named usr1 on your local machine:
; FileName1=/usr1/rigorous_ortho/005606990010_01_P008_MUL/05JUL11182931M1BS-005606990010_01_P008.TIF
Note: You dont need to edit the FileName0 line, since the GCP locations in
phoenixGCPs.pts only pertain to the 05JUL* image (FileName1).
3. Save the file phoenixGCPs.pts.
4. When working with multiple input images, you need to specify which image corresponds to the
GCPs you are about to restore or add. In the Wizard, select the 05JUL* image filename from the
Active Image drop-down list:
Page 10 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
Page 11 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
9. Enter the following values in their respective fields in the Wizard, then click the Add GCP button
.
Lat: 33.39088287
Lon: -111.89436526
Image X: 1279
Image Y: 1631
Elevation: 333.136
10. Enter the following values in their respective fields in the Wizard, then click the Add GCP button
.
Lat: 33.41914515
Lon: -111.90030564
Image X: 1058
Image Y: 385
Elevation: 329.81
The table in the GCP Selection panel should look as follows. The new GCP locations that you
added are shown in rows 15 and 16. ENVI Classic internally converts the geographic lat/lon
coordinates into eastings and northings for this table. The locations of all 16 GCPs are shown in
the Layout Manager.
Page 12 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
1. In the Controls tab of the Layout Manager, the RMSE (root mean square error) field lists the
overall model error.
2. Enable the Residual Image option in the Layout Manager. A contour map is added to the Draw
Area, showing the residual error throughout the area based on your GCP data. Purple-to-blue
areas indicate the lowest residual errors, and red-to-white areas indicate the highest errors.
3. Right-click inside the Draw Area, and select Zoom In until you can see the contour labels more
closely. Click-and-drag the left mouse button around the Draw Area to pan around. From looking
at this contour view, you can see that ID #3 (in the far lower-left corner of the contour map) has a
relatively high residual error. IDs #1 and #14 have a relatively low error.
4. Select the GCP tab in the Layout Manager, and scroll to the right until you see the Mag field. The
contour view of residual error is based on these relative magnitude values. You must have four or
more GCPs for this column to appear.
5. Select Row #2, which corresponds to GCP ID #3. A diamond marker appears over the GCP. This
is the GCP associated with the large error you saw in Step 2. If you scroll to the Mag field for
this point, you will see a value of 206.2582. This value is very high, compared to the Mag values
of other GCPs.
6. Select the Controls tab in the Layout Manager, then enable the Residual Vectors option. Use the
arrow button (to the right of the slider) to increase the Residual Exaggeration value to 10. Cyancolored arrows originate from each GCP location, showing their error magnitudes and directions.
The left side of the scene has the highest errors and, as expected, the GCP corresponding to ID #3
has the largest residual error vector.
In the next step, you will remove this GCP and see if the overall model improves.
7. In the Wizard, click Row #3.
Page 13 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
8. Click the Remove GCP button . The entire model is recomputed, and the updated residual errors
are shown in the Layout Manager. What is the overall RMSE value now?
9. Click the GCP tab in the Layout Manager, and evaluate the updated Mag values for each GCP.
Did the overall model improve after deleting that GCP?
10. You can experiment with removing GCPs or editing their values to improve the overall model
error. For this tutorial, you can assume the results are reasonable and proceed to selecting tie
points. Click Next in the Wizard to proceed to the Tie Point Selection panel.
Page 14 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
Cutlines are polygons that are used in combination with image ordering to define areas that you want to
keep in the final output. A cutline is essentially an inclusion polygon, similar to a cookie cutter applied to
one or more images. The following example shows a cutline for a single image:
Following is an example of applying a cutline to the first-order (top-most) image when you have multiple
images.
Page 15 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
For this tutorial, you will skip the process of ordering images and defining cutlines, by clicking Next.
ENVI Classic will use the image frame as a cutline by default. The image frame is the area of usable
image data, excluding any background pixels.
Page 16 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
To save some time, the pixel coordinates for each tie point were already determined for you. This
process involved identifying an object in the 05JUL* image (such as an intersection or building
corner), zooming in, and using ENVI Classics Cursor Location/Value tool to record the
fractional image coordinates for that pixel. The same object was then identified in the 05OCT*
image, and the fractional image coordinates for the corresponding pixel were recorded.
2. Enter the following pixel coordinates that represent the first tie point, then click the Add Tiepoint
button
Image 1
Image X: 6737, Image Y: 511.75
Image 2
Image X: 363.75, Image Y: 480
3. Enter the following pixel coordinates that represent the second tie point, then click the Add
Tiepoint button.
Page 17 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
Image 1
Image X: 6860.25, Image Y: 1997.75
Image 2
Image X: 489.75, Image Y: 1928
4. Enter the following pixel coordinates that represent the third (and final) tie point, then click the
Add Tiepoint button.
Image 1
Image X: 6835.5, Image Y: 3174
Image 2
Image X: 471.75, Image Y: 3073.5
The tie points are shown in the Layout Manager with magenta-colored X symbols, and their
values are added to the Wizard and the Tie Points tab in the Layout Manager. These dialogs
should look similar to the following.
Page 18 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
5. Before proceeding to the next step in the workflow, you should save your project. From the
Wizard menu bar, select File > Save Project. The Select Project File dialog appears.
6. Enter a filename and location for the project state, and click Save. ENVI Classic saves the
project state in XML format. If you need to restore your project later, select Map > Rigorous
Orthorectification from the ENVI Classic main menu bar, then select File > Restore Project
from the Wizard menu bar to select the XML file for your project. ENVI Classic automatically
opens all input files, DEMs, GCPs, and tie points, so that you dont need to open them prior to
restoring the project.
7. Click Next in the Wizard to proceed to the Image Order & Cutline Selection panel.
Page 19 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
By default, the 05OCT* image (the eastern image) was at the top of the Image Order list back in
the Image Order & Cutline Selection step. So in the area of overlap between the two images, the
05OCT* image is ordered above the 05JUL* image. You can see this in the above figure.
6. Select File > Exit from the Wizard menu bar. You will be prompted to save the state of the
orthorectificaton process in case you want to resume work later.
Optional Step
If you have some extra time, you may want to experiment with skipping the steps of adding GCPs and tie
points altogether, then re-running the orthorectification process. Compare the RMSE value of the overall
model with no GCPs or tie points.
Page 20 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
GCPs and tie points should have a quality that is consistently better than the initial accuracy of the
sensor model; this varies with data type and vendor. For example, mid-latitude QuickBird scenes of a
relatively flat area may provide an accurate orthorectification without the use of GCPs or tie points.
Copyright Notice:
ENVI Classic is a registered trademark of Exelis Inc.
QUAC and FLAASH are registered trademarks of Spectral Sciences, Inc.
Page 21 of 22
2014 Exelis Visual Information Solutions, Inc. All Rights Reserved. This information is not subject to the controls
of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However,
this information may be restricted from transfer to various embargoed countries under U.S. laws and regulations.