Figure .  Button used to access the import/export action window.
Once the import/export button is selected, the Import/
Export dialog box will pop up (fig. 4 ).
To import an image, first choose the type from the 
pulldown menu “Type:”. If files are purchased from EDC and 
are Landsat 7 images, then the appropriate choice should be 
Landsat 7 Fast-L7A EROS. If images are taken from other 
sources, then other options may be appropriate. For example, 
many images from NASA sources are provided in an open-
source hierarchical-data format (HDF), which has various 
options on the Imagine import/export dialog. If images are 
from the LEDAPS site, they are in HDF format, which is a 
public-format heirarchical data type used by NASA. Select 
“HDF Scientific Data Format (Direct Read)” as the data type.
Once the appropriate input file type has been chosen, 
fill in the “Input file:” box, either by manually typing the full 
path name, or navigate with the mouse using the open folder 
button just to the right of the blank input-file field. Clicking 
the “Open folder” button to the right of each filename brings 
up a file finder window (fig. 5 ) that allows quick access of 
The “recent” and “goto” buttons on the right are 
particularly useful, as they allow quick reference to recently 
accessed files and folders, respectively.
Setting the default data directories: If the same directory 
folder is to be used often, it is advantageous to set that 
directory as the default directory in Imagine. To do this, 
select “Session/Preferences” from the Imagine main console 
to bring up the “Preference Editor” window. Select “User 
Interface & Session” from the “Category” list,  and then 
manually type in the desired path name into the Default Data 
Directory and Default Output Directory fields. Choose “User 
Save” to register these preferences for this user for this and 
all future sessions.
When the input filename is chosen, Imagine will place 
a default output filename in the “Output File” field of the 
Import/Export window. Unless specifically overridden, this 
file will be in the default directory folder. Since the default 
directory is not typically the desired folder for the output 
image, it will often be necessary to navigate in the output 
file field to the desired directory, and then retype the name of 
the output file. A workaround for this process is to navigate 
in the Output File field to the desired directory before the 
input file is chosen, and to simply type in a dummy filename 
in the output file field, and thereafter to select the input file 
as desired. This will ensure that the default name chosen 
by Imagine for the output file will be placed in the desired 
output directory.
Figure .  Import/export action window.
Figure .  File finder window.
Images from a wide range of sources can be imported into Imagine through the Imagine Import/Export 
action window. This action window is accessed through the second box on the main console (fig. 3).
SOP 2.  Preprocessing Landsat Imagery    
Pdf text searchable - search text inside PDF file in, ASP.NET, MVC, Ajax, WinForms, WPF
Learn how to search text in PDF document and obtain text content and location information
how to search a pdf document for text; find text in pdf files
Pdf text searchable - VB.NET PDF Text Search Library: search text inside PDF file in, ASP.NET, MVC, Ajax, WinForms, WPF
Learn How to Search Text in PDF Document and Obtain Text Content and Location Information in VB.NET application
find and replace text in pdf; how to select text in pdf image
Figure 6.  Open-file folder icon, circled.
IV. Reprojecting and Resampling 
It is critical that all images and vector coverages be 
in the same projection, spheroid, and datum before any 
radiometric or change-detection processing is begun! 
Although Imagine will tolerate different projections and 
spheroids when processing, it is best not to allow Imagine to 
reproject “on-the-fly” without knowing explicitly that images 
are aligned as they should be. Only by viewing spatial data 
and ensuring that spatial relations are as expected can robust 
results be achieved.
The geometric reference image is the base of all 
geometric processing, and thus it is important that it have the 
appropriate map properties before any further processing is 
done. Therefore, after the reference image has been imported 
to the Imagine format, it frequently must be reprojected and 
The NCCN standard spheroid/datum pair is GRS80/
NAD83, which should be used for all spatial data. Although 
many existing Digital Elevation Models and some vector 
coverages are in the Clarke 1866/NAD27 Spheroid/Datum 
combination, those are based on an older model of the Earth, 
a poorer datum-measurement network, and thus should be 
avoided as baselines for future monitoring. The conversion 
between NAD27 and NAD83, although not perfect, can be 
achieved with accuracy more than sufficient for Landsat-based 
All Landsat images are to be resampled to the 25 m pixel 
size. Although the native resolution of the imagery is 28.5 m, 
experience has shown that the 25 m size is far easier to work 
with, both for processing and for eventual field validation.
Therefore, it is necessary as a first step to check the 
condition of the reference image and reproject it if necessary.
1.  Start Imagine if that has not been done already.
a.  If Imagine asks whether a Classic viewer or a 
geospatial light table is desired, select the classic 
viewer. It may be desirable to choose not to be asked 
this question again.
2.  In Viewer #1, load the geometric-reference image.
a.  Click on the folder icon directly under File on the 
viewer menu (fig. 6 ).
b.  Navigate to the folder with the reference image and 
open it.
i.  Note: This assumes that the reference image has 
been imported into Imagine’s format. If not, go 
back to section III above.
ii.  The image should come up in the viewer.
(1)  Because of the tilt of the Landsat image 
relative to the earth, the upper left portion 
of the image area likely is black; because 
the viewer defaults to viewing this area, 
the initial visible area may be black. To see 
more of the image, navigate around with the 
sidebar and bottom bars.
(2)  Alternatively, right-click in the viewing 
space of the viewer to see the viewer options 
menu, and select “Fit image to window.” 
This option will resize the entire image to 
fit within the viewing space of the viewer. It 
may take several seconds for the resizing to 
6    Protocol for Landsat-Based Monitoring of Landscaped Dynamics at North Coast and Cascades Network Parks
C# PDF Convert to Text SDK: Convert PDF to txt files in
NET project. Powerful .NET control for batch converting PDF to editable & searchable text formats in C# class. Free evaluation library
convert pdf to searchable text; pdf searchable text
C# Create PDF Library SDK to convert PDF from other file formats
The PDF document file created by RasterEdge C# PDF document creator library is searchable and can be fully populated with editable text and graphics
how to select all text in pdf; search text in pdf image
Figure .  ImageInfo button, circled.
3.  Note the geometric properties of this image.
a.  Select the ImageInfo icon (a white box with a lower-
case “i” in it (fig. 7 ).
b.  An ImageInfo window will pop up.
i.  On the front panel of this window, labeled as the 
“General” tab, the geometric information about 
the reference image is shown.
c.  At the bottom of the ImageInfo box is the Projection 
i.  Ensure that the image is in Universal Transverse 
Mercator (UTM), in the correct zone (zone 10 
for all of the parks of the NCCN), and that the 
spheroid and datum are as desired (NAD83, 
unless compelling reasons exist to use a 
different datum).
(1)  See Bolstad (2002) for more information on 
datums and spheroids.
ii.  Ensure that the pixel size is 25.0 m.
(1)  If the reference image has been acquired 
recently, it was recommended in SOP 1 
that the image be purchased at 28.5 m. 
Therefore, it is expected that the pixel 
size will be 28.5 m, which will need to be 
resampled to 25.0 m.
d.  If either the map-projection information or the 
pixel size is not as desired, then reprojection will be 
4.  Reprojecting and resampling images in Imagine
a.  If the reference-geometric image is in a projection 
or a spheroid/datum other than desired, it is 
straightforward to reproject the image in Imagine. 
Reprojection also can be achieved in Arc programs, 
but details of the Imagine-based approach are 
b.  Steps:
i.  Ensure that the image to be reprojected is open 
in a viewer.
(1)  See step 2 above if not.
ii.  From within that viewer, select “Raster/
Geometric Correction”
(1)  An action window named “Set Geometric 
Model” will popup, with a list of options.
iii.  Scroll down and select “Reproject” from the list.
(1)  A “Reproject Model Properties” window 
will pop up, as well as a “Geo Correction 
Tools” window.
(2)  In the Reproject Model Properties window:
(a)  Change grid sampling X and Y to 30 
(i)  Note that the number 30 is not 
the size of pixels. This number 
refers to the density of a grid that 
Imagine will construct to solve 
the reprojection equation between 
different projections. The number 
30 in X and Y means that a 30 
by 30 grid (i.e. with 900 points) 
will be used to project between 
(b)  Select the “Projection” tab.
(c)  Click on “Add/Change Projection.”
(i)  A “Projection Chooser” window 
will come up.
(d)  Select the desired projection into which 
the image is to be projected.
(e)  Click on OK button.
(3)  In the Geo Correction Tools window, select 
the “Resample” button, which is a box with 
four smaller multicolored boxes inside it, 
tilted on its side (fig. 8).
(a)  The “Resample” dialog will pop up.
(4)  For the Output File, navigate to the desired 
output folder before typing in the name of 
the resampled output file.
(a)  Follow naming conventions in SOP 5 
Data Management when reprojecting a 
SOP 2.  Preprocessing Landsat Imagery    
VB.NET Image: Robust OCR Recognition SDK for VB.NET, .NET Image
for VB.NET provides users fast and accurate image recognition function, which converts scanned images into searchable text formats, such as PDF, PDF/A, WORD
select text in pdf file; how to search text in pdf document
VB.NET PDF Convert to Text SDK: Convert PDF to txt files in
keeping original layout. VB.NET control for batch converting PDF to editable & searchable text formats. Support .NET WinForms, ASP
how to make a pdf document text searchable; cannot select text in pdf file
Figure .  Geo Correction Tools window, with the 
resample icon circled.
(5)  For the resample method, select “Cubic 
(a)  Note: Nearest-neighbor (NN) 
resampling is often recommended 
to maintain the radiometric integrity 
of the imagery, since it leaves the 
spectral response within a pixel 
unchanged. However, experience in 
monitoring at the NCCN parks has 
shown that the geometric distortions 
introduced by the NN resampling 
cause artifacts in later change 
detection. Therefore, it is important 
that the cubic-convolution approach 
be used. Cubic convolution (see 
Richards (1993) for a discussion) 
retains the geometric properties of the 
image much better, especially when 
pixel size is changing as in this case, 
but does introduce some spectral 
anomalies at the edges of contrasting 
pixels that must be considered in later 
(6)  For the output cell sizes, manually type in 
25 for both X and Y fields.
(7)  Select “Ignore Zero in Stats.”
(a)  This will make later display of the 
image work better. Without selecting 
this option, the large area of no-data 
around the margins of the imagery 
is included in the statistics used by 
Imagine to display the image. Images 
will appear “washed-out.”
(8)  Click on OK button.
(a)  A progress window will appear and 
notify when the job is complete.
(9)  Close all windows except the Viewer 
#1, and do NOT save the changes or the 
geometric model used for the resampling.
(10)  The resultant 25 m image will be clipped 
to the study area in section V below.
5.  Reprojecting the Digital Elevation Mode (DEM)
a.  A DEM is necessary for a variety of steps in the 
protocol, and must be in the same projection and 
datum as the geometric-reference image. Because 
DEMs typically exist in the NAD27 datum at 30 
m, it is often necessary to reproject these into the 
desired NAD83 and 25 m condition.
i.  NOTE: If the DEM is in NAD27 and is a grid 
coverage, reproject it in an ESRI Arc software 
product using the NADCON option for 
reprojection. This method will be more precise 
than using Imagine’s approach in steps b. and c. 
b.  Open the DEM in Viewer #1.
i.  Properties of the DEM:
(1)  A digital-elevation model for the areas around 
the parks has been included by the authors, 
and thus this step will not be necessary unless 
new DEMs become available, or if the study 
area around the parks is changed. These steps 
are provided for that case.
(2)  If a digital-elevation model is acquired 
from another source, or if an in-house DEM 
from the park is used, it must have metadata 
describing the map-projection and datum 
information, pixel size, and units of the 
vertical dimension (meters or feet).
ii.  Open the DEM in a viewer.
c.  Reproject the DEM to NAD83 and 25 m.
i.  Follow steps in section 4 above in this section, but 
use the DEM as the base rather than the geometric 
reference image.
d.  If reprojection is done in an Arc software product, 
it will be necessary to import the reprojected DEM 
GRID into Imagine. See section III above and follow 
the instructions for importing images into imagine. 
The GRID format is listed as one of the import format 
i.  After importing, ensure that the DEM has 
imported successfully by viewing it in a viewer 
(see step 2 in this section above).
    Protocol for Landsat-Based Monitoring of Landscaped Dynamics at North Coast and Cascades Network Parks
C# HTML5 Viewer: Load, View, Convert, Annotate and Edit Word
C# users can convert Convert Microsoft Office Word to searchable PDF online, create multi empowered to add annotations to Word, such as add text annotations to
search pdf files for text programmatically; select text in pdf reader
Online Convert PDF to Text file. Best free online PDF txt
PDF document conversion SDK provides reliable and effective .NET solution for Visual C# developers to convert PDF document to editable & searchable text file.
pdf text searchable; find text in pdf image
V. Defining Study Area, Clipping 
Images and DEMs to the Study Area
This protocol focuses on monitoring in and around 
the parks, and as such does not require that processing be 
conducted on an entire Landsat image. Those images are 
approximately 180 by 180 km, and are located in a manner 
completely irrelevant to the parks. Therefore, before any 
processing is done, it is necessary to define a study area in 
which all image-processing and change-detection activities 
will take place.
The end goal of study area delineation is to have two 
geographic coverages representing the same study area 
boundary: an Imagine-format Area of Interest (AOI) and an 
ESRI Arc-format vector coverage. Both can be created in 
Imagine, or, if the study area has been defined elsewhere and 
is currently in a vector format, Imagine can be used to make 
an AOI of the vector coverage.
The delineation of the study area is an important 
step in the long-term success of this protocol, as the study 
area defined here will be used for reference when creating 
summary metrics of change over time. Therefore, this step 
should involve the input of as many interested parties at the 
parks as possible. As of this writing, there was no agreed-upon 
rules by which such a designation could be made generic, 
but there are several components to consider. First, of course, 
the study area should include the entire park. There may be 
interest in creating separate study areas in some parks (such 
as the coastal region at OLYM), but this decision should be 
approached carefully as it requires more processing. Second, 
the study area may extend beyond a park by either an arbitrary 
distance or by some distance defined by ecological or 
economic factors. It may be difficult to predict these factors at 
the current time, which argues for a conservatively large study 
area that could allow future contraction. It will be significantly 
easier to make the study area smaller in the future than it 
would be to make it larger.
If the desired study area already exists as a vector layer, 
the following steps are required:
1.  Start Imagine and load the geometric reference image in a 
viewer (see section IV, steps 1 and 2, if this has not been 
done already). Ensure that this is the geometric reference 
image with 25 m pixels.
2.  In the viewer, choose “File/Open/Vector Layer” and 
navigate to the directory where the new study area vector 
layer is housed.
a.  Open the vector layer, and ensure that it matches with 
the reference image as expected.
b.  Use the ImageInfo (section IV, step 3 above) to 
ensure that the vector coverage is in an appropriate 
projection and datum for the geometric-reference 
i.  Assuming that the geometric reference will be 
UTM and NAD83, respectively, then the vector 
coverage must be in UTM, NAD83.
(1)  NOTE: If a vector coverage is in the Clarke 
1866/NAD 27 projection/datum space, it 
must be reprojected. It is recommended 
that this projection process occur in 
ESRI Arc software, using the NADCON 
transformation approach. Arc/Info 
documentation on “Projections” provides 
more information. However, if this approach 
is not possible, the steps for reprojecting a 
vector coverage in Imagine are given below.
(2)  If the coverage is not in the appropriate 
datum, it will be necessary to reproject it 
into the correct datum. This can be done in 
various ESRI Arc products and in Imagine. 
The steps for doing this in Imagine are 
detailed below. GIS professionals also 
may prefer to do this in Arc workstation or 
ii.  Reprojecting a vector coverage in Imagine.
(1)  In a blank viewer, open the vector coverage.
(2)  Click on the ImageInfo icon (section IV, step 
3 above).
(a)  The window associated with the 
“General” tab will pop up.
(3)  Select “Edit/Reproject the coverage.”
(a)  A “reproject coverage coordinates” 
window will pop up.
(4)  Name the output vector in such a manner 
that the distinction in datums can be 
inferred. Follow data-management 
considerations in SOP 5 Data Management.
(5)  Click on “Define new projection.”
(a)  The “Projection chooser” window will 
pop up.
(6)  Use the pulldown menu for the spheroid 
name field to select NAD83.
(7)  Use the pulldown menu for “Datum name” 
to pick the NAD83 datum from the list.
(8)  Ensure that zone is 10 and “North or South” 
is set to north.
(9)  Click OK, then click OK again in the 
“reproject coverage coordinates” window.
(a)  Save the coverage.
SOP 2.  Preprocessing Landsat Imagery    
VB.NET Create PDF from Text to convert txt files to PDF in
Best VB.NET adobe text to PDF converter library for Visual Studio .NET project. Batch convert editable & searchable PDF document from TXT formats in VB.NET
pdf find and replace text; search pdf for text
C# Create PDF from Text to convert txt files to PDF in, ASP
Visual Studio .NET project. .NET control for batch converting text formats to editable & searchable PDF document. Free .NET library for
how to select text in a pdf; search text in pdf using java
Figure .  Vector viewing properties for test window.
Figure 10.  Vector-attributes window.
3.  With the vector coverage and the geometric-reference 
image open in the same viewer, select “Vector/Viewing 
a.  The “Properties for test” window for the vector 
coverage will pop up (fig. 9).
b.  When vector layers are first shown in Imagine, the 
default is to display them as arcs, not polygons. To 
manipulate polygons, as will be done below, the 
polygons must be displayed.
c.  First, uncheck the box for “Arcs” by clicking it.
d.  Then check the box to the left of “Polygon.”
i.  To change the border and fill properties of the 
polygon, select the box to the right of the colored 
field to the right of the polygon (in cyan in 
fig. 9).
ii.  Keep this window open for now.
4.  Then from the viewer with the vector coverage in it, select 
a.  The “Attributes” window will pop up.
i.  Each polygon will have a row in the table 
of attributes. Select all the rows that include 
polygons that are desired to defined the study 
area. Typically, this should be a small number of 
polygons – often only one – because the study 
area should be a large, internally connected 
(1)  To select a row, click on the row number 
in the left-hand side of the window (the 
“Record” column in fig. 10).
ii.  To select multiple rows, click and drag down all 
iii.  Because polygons are displayed, selecting these 
rows will change all of the polygons on the 
viewer to the selection color, which defaults to 
5.  With these polygons selected, select “AOI/Copy selection 
to AOI” from the viewer.
a.  Note: If the coverage and the image are not in the 
same projection, this step will fail. It is not sufficient 
to simply remove the image from the viewer to avoid 
this failure, because the AOI must be associated with 
the correct datum of the image for all later steps. 
Therefore, if the coverage is not in the correct datum, 
it must be either reprojected or the datum must be 
overruled, as noted in step 2 above.
b.  This step creates a new AOI in the viewer, anchored 
to the geometric reference image.
6.  Save the AOI.
a.  In the viewer, select “File/Save/AOI Layer as.”
b.  Navigate to the appropriate directory and save the 
AOI layer. Follow the SOP 5 on naming and file 
location conventions.
7.  This AOI is now the base AOI for all subsequent steps 
that involve clipping of imagery to the study-area 
boundary, including clipping the DEM.
 If the study area does not exist in vector form and must be 
created, follow steps 8 through 11.
0    Protocol for Landsat-Based Monitoring of Landscaped Dynamics at North Coast and Cascades Network Parks
VB.NET Create PDF from Word Library to convert docx, doc to PDF in
Export all Word text and image content into high quality Professional .NET PDF batch conversion control. Easy to create searchable and scanned PDF files from
pdf search and replace text; search a pdf file for text
VB.NET Create PDF from PowerPoint Library to convert pptx, ppt to
PDF, VB.NET convert PDF to text, VB.NET multiple pages PowerPoint to fillable and editable PDF documents. Easy to create searchable and scanned PDF files from
convert pdf to word searchable text; convert a scanned pdf to searchable text
8.  In a viewer, open the geometric-reference image.
a.  See section IV, steps 1–2 if this is not done already.
9.  Draw the study area as an AOI.
a.  Select “AOI/Tools” from the viewer, and then select 
the polygon drawing button (fig. 11 ).
b.  When the cursor is moved across the viewer, it will 
take the form of a crosshair. Clicking with the left 
mouse button will create a vertex in the developing 
polygon. Define the vertices that will define the 
polygon desired to encompass the study area. 
Double-click when the last desired vertex has been 
placed to close the polygon-creating process.
i.  NOTE: The study area for ecological monitoring 
may well extend beyond the boundaries of the 
parks proper.
10.  Save the AOI as the study area AOI.
a.  In the viewer, select “File/Save/AOI layer as.”
b.  Navigate to the appropriate directory and save the 
AOI layer. Follow the SOP 5 on naming and file-
location conventions.
c.  This AOI is now the base AOI for all subsequent 
steps that involve clipping of imagery to the study 
area boundary, including clipping the DEM.
Figure 11.  AOI tools window, 
with the polygon drawing button 
Figure 12.  Arrange layers window.
11.  Copy this AOI to a vector coverage.
a.  NOTE: This step is designed to create a vector 
coverage of the study area for use in displaying data 
in ESRI software. It is not strictly necessary, as it 
will not be used for analysis in any other part of this 
b.  It is first necessary to create a blank vector coverage.
i.  In the viewer, ensure that the geometric-
reference image is displayed.
ii.  Select “File/New/Vector layer.”
(1)  A “Create a New Vector Layer” window will 
pop up.
iii.  Navigate to the appropriate directory and type in 
the name of the coverage in the Filename field.
(1)  Follow the conventions in SOP 5 Data 
Management for placement and naming of 
the coverage.
iv.  Select “single precision” when given the choice 
between single and double precision.
c.  Once this blank vector layer is open, the viewer 
should show “Vector” in the list of options at the top 
of the viewer.
d.  In the AOI tools window, click on the selection arrow 
icon (upper left icon in fig. 12 ).
i.  Click in the viewer on the space of the AOI to 
select it.
(1)  A white box will appear around the outside 
perimeter of the study area AOI. This 
indicates that the AOI is selected.
SOP 2.  Preprocessing Landsat Imagery    1
ii.  Select “Vector/Copy selection to vector.”
(1)  Little will seem to have happened. Do not be 
iii.  Select “File/Save/Top layer.”
e.  Confirm that the vector coverage is as expected by 
removing the study area AOI from the viewer.
i.  Select “View/Arrange Layers.”
(1)  The “Arrange layers” window will pop up.
(2)  Right-click on the top layer, the AOI layer. 
Select “Delete layer” from the list. Then 
click on the “Apply” button in the lower left. 
This will remove the AOI from the viewer. 
Do not save the AOI layer, since this has 
been done already in step 10 above.
(3)  After the AOI is removed, the outline of 
the study area drawn on the screen should 
f.  Build the vector layer.
i.  In the main imagine console, select the Vector 
layer utility.
(1)  Note: If the Vector layer utility is not 
included in the Imagine installation being 
used, this layer can be built in ArcGIS or 
Arc workstation. See the documentation for 
those software packages for assistance on 
building polygon coverages.
ii.  Select “Build Vector Layer Topology”
(1)  Select the vector coverage just defined.
(2)  Click OK.
Once the study area AOI has been defined, the 
geometric-reference image and the DEM must be clipped 
to that area.
12.  From the main Imagine icon panel (fig. 1  above), click on 
the “Intepreter” button.
a.  The “Image Interpreter” window will pop up.
b.  Click on the “Utilities” button from this window.
i.  The “Utilities” window will pop up.
ii.  Select “Subset.”
(1)  The “Subset” window will pop up (fig. 13 ).
iii.  In the “Input File” field, navigate to the directory 
with the geometric-reference image and select it.
(1)  Ensure that the selected image is the one 
with 25 m pixels and in the NAD83 datum.
iv.  In the “Output File” field, navigate to the 
appropriate directory for the geometric-reference 
Figure 1.  Subset window.
(1)  Follow the conventions in SOP 5 for 
locating and naming the reference image.
v.  Then click on the “AOI” button at the bottom of 
the Subset window (fig. 13).
(1)  In the “Choose AOI” window, select “AOI 
(2)  Use the folder button on the right of the AOI 
file field to navigate to the study area AOI 
created in step 6 or 10 above in this section.
vi.  Then click on OK button in the Subset window 
to start the clipping process.
13.  Conduct the same process on the DEM (step 12), except 
use the DEM file as the input file in the Subset window 
instead of the geometric-reference image.
At the end of the steps in this section, there should be the 
a.   A study area AOI that defines the bounds of the area 
where all subsequent work will be done.
b.   A vector coverage of the same area.
c.   A geometric-reference image clipped to the study area.
d.   A DEM clipped to the study area.
2    Protocol for Landsat-Based Monitoring of Landscaped Dynamics at North Coast and Cascades Network Parks
VI. Geometrically Registering Images 
to the Reference Image
Images must be aligned so that the information in a 
pixel at a given coordinate in an image corresponds to the 
same point on the ground over multiple dates of imagery. 
To do this, all images must be geometrically aligned to 
the reference image. Because that image is assumed to be 
geometrically accurate and registered to a map projection with 
known coordinates, geometrically registering other images 
to that reference will result in those new images also being 
registered to map coordinates. This allows not only direct 
comparison between images, but also allows connection with 
other spatially registered datasets, such as GIS layers with 
information on roads, streams, and trails.
Registration of one image to the reference image is a 
multistep process (Richards, 1993). First, tie points are in 
each image. A tie point is a location that can be identified 
unambiguously in both the reference image and the dependent 
image, here known as an “Input” image. After a sufficiently 
large number of tie points is found linking the coordinates 
in the input image to those in the reference image, a 
mathematical transformation in the form of a linear equation 
can be solved that projects a coordinate in the input image into 
the map system of the reference image. Because of the need 
to terrain-correct imagery, the process of solving this equation 
also must take into account the viewing geometry of the sensor 
and the topography of the land being imaged. Finally, once 
the equation for transforming coordinates is solved, all of 
the coordinates of the input image are transformed into their 
coordinates of the reference image, and resampled to place 
them on a regular grid in the reference-image map system.
The first step in this process, finding tie points between 
images, can be done manually or in an automated fashion. 
Manual location of tie points is time consuming, subject to 
the interpretation of the individual interpreter and thus not 
repeatable, and typically results in tie-point populations that 
are small. Therefore, an automated approach is desirable, and 
is provided as a module in IDL Virtual Machine (VM) called 
Locating Image Tie Points Using ITPFind
ITPFind uses image correlation to locate coordinates 
of image tie points (ITPs) that best match the pixel patterns 
of small windows in the reference and the input images. It is 
described in detail in Kennedy and Cohen (2003).
The program is supplied as a module to be run in the 
IDL-VM. IDL-VM must be installed on a computer before 
this module can be used (see
index.asp, as noted above, for downloading instructions). Once 
IDL-VM has been installed, the program can be run by simply 
clicking and dragging the ITPF module onto the IDL-VM 
icon. Before this is done, several preparatory files must be 
setup. The program is named “ITPFind.sav.”
ITPFind requires text files as inputs. These files describe 
the location and properties of the images to be registered and 
the parameters used by the program for choosing windows of 
correlation, grid spacing of tie points, and thresholds for tie-
point acceptance. The authors have provided examples of each 
of these files as part of this protocol, along with the ITPFind 
VM module. These files must be set up before the program 
can be used to register images.
1.  The first file needed is named “squeet_master.txt”. This 
file provides the pathnames for all subsequent text files 
needed by the program. This file can be located anywhere, 
but it is best to place it in a directory relevant for the 
files of interest so its location can be found easily. When 
the ITPFind module is run, it will query the user for 
the location of this file. The file contains the path and 
filenames of the two other pointer files required for the 
program to begin functioning.
a.  One pointer file contains the list of available 
parameter files (its contents are described in (2.) 
below). Identify the name of this pointer file using the 
following syntax: <path/file_here_without_brackets>.
2.  The other pointer file contains the list of available image 
files (its contents are described in (3.) below). Identify 
the name of this pointer file using the following syntax: 
b.  File of available parameter files — The filename for 
this file is specified in squeet_master.txt. This file 
contains the list of filenames of candidate parameter 
files. The list of candidate-parameter files starts at the 
beginning of the file, one file per line separated by 
carriage returns, with no leading text in the file or on 
any line. The contents of the parameter files will be 
outlined below.
3.  File of available image files — The filename for this file 
is specified in squeet_master.txt. This file contains the 
list of available images and their vital characteristics. For 
each image file, you must include critical information 
for the program to run. The following is an example. 
SOP 2.  Preprocessing Landsat Imagery    
Each line will be described below the example. Note that 
each line ends with a carriage return. Also, there must 
be a colon and a space between the colon (:) and the first 
character of the associated information (i.e. Rotation: 0, 
not Rotation:0 or Rotation 0). Text that does not exactly 
match the text before the colons will be ignored and will 
cause errors.
Filename: images/exampleInpImage.img
File codename: inputimage
Start point:(x,y): 388, 663
Pixel center to pixel center distance:(x,y): 25, 25
Rotation: -6
Layer to use: 4
Ignore: 0
Mask below: 10
Mask above: 250
Description of Each Line:
Type in the full pathname of the image without quotes. 
Currently, the image formats that are accepted are ERDAS 
Imagine format, or a flat-binary file denoted with a “.bsq” 
suffix. If the flat-binary option is used, however, there MUST 
be an accompanying file that describes the geometric qualities 
of the image. See “Supported Image Types” below.
Currently, the software is designed to accept either 
ERDAS Imagine 8.3–8.7 images or flat-binary files. Imagine 
types supported are: Unsigned 8 bit, Unsigned or Signed 16 
bit, and IEEE Float.
Flat-binary files must be stored in band sequential 
format, and the filename must end with “.bsq”. Also, 
such files must have an accompanying file named the 
same as the image file, except with “.hdr” swapped in 
the place of “.BSQ” (i.e. image.bsq and image.hdr). 
The .hdr file must contain at least the following entries 
in the following order:
Text on other lines will be ignored.
BANDS: the number of bands in the image 
ROWS: the number of rows (lines) in the image 
COLS: the number of columns (pixels) in each line of 
the image 
DATATYPE: use 3 for 8-bit unsigned, 5 for 16-bit 
unsigned, and 9 for 32-bit IEEE floating point.
<***>_COORDINATE: All coordinates can be 
written as decimal values or integer values. These 
must be in the coordinate system units of the image. 
If that system is unknown, assume that each pixel is 
offset from the next by a value of 1.0. It is critical, 
however, that the coordinates be referenced as if they 
were a UTM-type map coordinate system for the 
northern hemisphere, where Y-value increases from 
bottom to top. File-coordinate system coordinates 
typically work in reverse (Y increases from top to 
bottom), and this should be avoided.
PIXEL_WIDTH, _HEIGHT: These should be 
reported in the units of the map system for the 
image, and should be consistent with the UL_, LR_
COORDINATES given above. For example, if the 
coordinates given above are in the UTM system with 
meters as the basis, then pixel dimensions also must 
be in meters (i.e. 28.5 for raw LANDSAT TM data, 
or 1000.0 for resampled 1km AVHRR data). Again, if 
the coordinate system of the image is not known, use 
1.0 for both X and Y pixel sizes.
File Codename:
This codename will be incorporated in the output ITP 
files. It allows later identification of which images were used 
to find the ITPs. Use some short, unique identifier string 
without quotes. Because the input and the reference images 
are stored in this file and there is no distinction in how they 
are stored, it is sometimes useful to include as part of the 
codename a reminder about which image is to be used as a 
reference and which as an input image.
Start Point:(x,y):
The map coordinates of the ‘start point’ of this image. 
It is simply the approximate location of a point in this image 
that will match with the ‘Start point’ of its complementary 
reference or input image. This start point must be in the 
coordinate system of this image: i.e. in the coordinates of the 
input image for the input image, and in the coordinates of the 
reference image for the reference image. Note that the method 
for location of these coordinates is described in Running the 
ITPFind Module section below.
    Protocol for Landsat-Based Monitoring of Landscaped Dynamics at North Coast and Cascades Network Parks
Documents you may be interested
Documents you may be interested