LAB 5. Assigning Coordinates to Images

Raster images in several formats, including JPEG and TIFF, can be georeferenced to coordinate systems if certain parameters are known. This allows co-registration of the images to vector data. For this lab, there will be two maps, Lab5a and Lab5b.

Lab 5a

The easiest method is to use the Georeferencing tool

Georeferencing

If the tool is not visible, please activate it by going to the “Tools” pull-down menu and look under both “Customize” and “Extensions” to check the activation boxes.

Step 1. Use ArcGIS to open Basemap.jpg (download it first - use save target as using the right mouse button).  Make sure its extension is .jpg, not.jpeg. DO NOT BUILD PYRAMIDS.

Step 2. Click on “Georeference” and unselect “Auto Adjust”

Step 3. Add ground control points at the center of the tic marks on the map using the tool that looks like two Xs connected by a diagonal line. Click on this tool and then left-click once on a tic mark. Then move you cursor off of the tic (you will notice a line it creates) and then right-click and choose “Input X and Y”. You can then type the coordinates in the pop-up box.

Step 4. After you have three ground control points completed you need to go back to “Georeference” and select “Auto Adjust”. The map will now likely disappear. To get it back, click the “Full Extent” button.

Step 5. Go back to “Georeferenceing” and choose rectify and follow the instructions.

Step 6. Close your current map document and open a new one. You now need to find the rectifed file using ArcCatalog. Right-click on the image name and choose “Properties”. Now scroll down to “Spatial Reference” and define the projection as Utah State Plane, North Zone, NAD 27.

Check your map accuracy by putting the cursor on each of the tic marks and comparing the coordinates given by ArcGIS to those written on the map.

<>Finally, make a layout using your raster map, place a state plane grid (find in Data Frame Properties) around it and include a title, north arrow, and scale bar, and mention the spatial accuracy. Also, put a second neat line around the entire map and ancillary graphics using the graphics rectangle tool.

Lab 5b

For the second part of Lab 5
download this data and unzip it. This creates a new folder called Georef containing your data.

Open two ArcGIS sessions add the geol_tr.img to one session and the sections Shapefile to the other. DO NOT BUILD PYRAMIDS!!!

You will find the coordinates of section corners and use those to rectify the image.

Use Select by Attribute (Selection pull-down menu) queries to find sections of interest. Start with this query:

"TNUM" = 30 AND "RNUM" = 12 AND "SECTION" = 10

Zoom in on this section and click on the lower-right (southeast) corner with the Identify Tool. The UTM coordinates will appear next to the Location annotation in the pop-up window. You can cut and paste these, but make sure to remove any commas.

Next click on Georeferencing (if this tool is not present open the Tools pull-down menu, click on Extensions and make sure it is checked and then click on customize (Tools once again) and check Georeferencing) and uncheck Auto Adjust. Then find the southeast corner of the correlative section on the geologic map image and left-click on it with the Add Control Points tool. Then move the tool slightly off of the point and right-click to open the input box for the coordinate. Either type or paste the coordinate in.

Repeat this process for four control points and then check Auto Adjust. You will then need to click the Full Extent button to see your map again. Open the ground control table and check your RMS error. If it is high check for errors.

You can then go back to the Georeferencing pull-down menu and select Rectify. Save your rectified image and then open ArcCatalog. Right click on the rectified image and select Properties. Then scroll down to Spatial Reference and edit for UTM Zone 12, NAD 83.

You will need to open a new Map Document and add the rectified data. If you add the rectified data to the same Map Document you did the georeferencing in, it will not show the proper coordinates.

Make a new map with the rectified image overlain with the sections Shapefile. Make a map in Layout View and turn it in.

Below are instructions for making a World file – this is not part of the assignment, just old information that may come in handy sometime.

Goals:

1. Determine the pixel (grid cell) size both in X and Y.

2. Determine the XY coordinate of the upper-left image (raster) pixel.

3. Build a World file -a text file with coordinate and pixel size information.

4. Use ArcCatalog to generate spatial reference

Proceed as with this assignment as follows:

Note: the map coordinates are State Plane, Utah North.

Step 1. Use ArcView to open Basemap.jpg (download it first, of course - use save target as using the right mouse button)  Make sure its extension is .jpg, not .jpeg.

Step 2. Open the above file in Arcmap.

Raster image coordinates are determined from the upper-left hand corner pixel. Note that a tick mark with a state plane easting and northing is near the upper-left corner. There are also tics in the upper-right and lower-left corners.

Step 3. Determine pixel size. Take the screen coordinates of the upper-left, upper-right, and lower-left tic marks. Do this by placing the cursor at the center of each and reading the coordinate at the bottom-right of your ArcView window. Then subtract the upper-left X screen value from the upper-right screen X value to determine the total number of pixels between the tic marks. Now subtract the upper-left map X coordinate (state plane feet) from the upper-right X coordinate. Divide this value by the total number of pixels you derived previously to get the pixel size, in feet, along the X axis. Now repeat this process for the Y direction (Upper-left & lower-left tics).

Step 4 - Determine the map coordinates for the extreme upper-left image pixel (not the tic!).

The logic here is similar to that for finding the pixel size. You determine the number of pixels from the upper-left tic to the top and left image margins. For the X coordinate, multiply the number of pixels between the tic and the left margin by the X axis pixel size in feet, and then subtract this value from the X map coordinate of the upper-left tic (easting values are decreasing to the west)

Now repeat this process for the Y coordinate. Determine the number of pixels between the upper-left tic and the top image margin. Multiply this by the Y axis pixel size in feet. Now add this value to the Y map coordinate value (values are increasing to the north).

These XY coordinates are used for the upper-left image pixel.

Step 5. Make a world file. Now a text file, called a world file, must be created. You can do this with any text editor such as Notepad or Wordpad. The extension of the text file must end in the first and last letters of the image type file extension and end in a w. For instance, a world file named test, for a JPEG (.jpg) file, would have a .jgw extension: test.jgw.

To build this file we use the following information:

1. Spatial resolution (pixel size) for X
2. A rotation term (zero for this)
3. A rotation term (zero for this)
4. Spatial resolution (pixel size) for Y
(must be negative)
5. The upper-left X coordinate in state plane
6. the upper-left Y coordinate in state plane

The completed file will look something like this (note: these are example coordinates and will not work for your assignment and the notes I have included are not added to the file):

5.00000000 (note: 5 meters, 5 feet, 5 kilometers, etc)
0.00000000 (note: rarely used, usually zero)
0.00000000 (note: rarely used, usually zero)
-5.01200000 (note: 5 meters, 5 feet, 5 kilometers, etc)
14569093.00000000 (note: map coordinate, i.e. UTM)
342900.00000000 (note: map coordinate, i.e. UTM)

For "basemap.jpg" the world file would be named "basemap.jgw"

Add the world file to the same folder as the JPEG image.  Open a new blank document in ArcMap ("File" pull-down menu) and reload the basemap.jpg as a new layer. This will allow ArcView to read the world file.

Now go to Data Frame Properties (under the View pull-down menu). Click on the Coordinate Systems tab. Use the "Predefined" folder to find State Plane NAD 27. Then choose Utah North. If you look into the folder containing the image, you will notice that ArcGIS added a .aux file containing the coordinate and datum information.

To check your accuracy, place the cursor on each of the tics and record their values.