ArcGIS Tips – DEM Edition

So you find yourself in a position where you have to make DEMs, but, like me, it’s been a while since you had to make one and are now sweatting as you peruse through the seemingly infinite menus and submenus of ArcGIS.

Displaying DEM data (bathymetry, elevation, georeferenced data points, etc.) on ArcGIS from a text file.

You have exported your georeferenced (yeah?) elevation (yeah?) data from that super simple program (ha) and now you’re at a loss for what to do with it in GIS. Here’s what you do:

To start with, DEM data generally comes in the ‘.txt’ format which is delimited (separated) using commas, and without spaces between the values (unless it is delimited using spaces). At the top you have a row of text which serve as column headers, which allow you to select that column of values within Arc.

  1. Import data using the ‘Add data’ function.
  2. Right click the imported file, select ‘Display XY Data’.
  3. If you have the appropriate headers in the text file you should be able to input your X, Y and Z fields (Easting, Northing, Depth etc.).
  4. Right click the new layer, select ‘Export Data’.
  5. This should be exported as a ‘.shp’ file.
  6. If this doesn’t work, select ‘Display XY Data’ again with the ‘Events’ file output from the previous steps.
  7. Follow steps 3 and 4, this should now work.
  8. After this, you can convert feature to raster.

If you want to use the bathymetry later, you can export to a TIFF (georeferenced, often called a GEOTIFF) using the export function (right click the raster layer you just made and export, select TIFF as the output filetype if not already selected).

Making a super professional hill-shaded bathymetry

  1. Complete the ‘Displaying DEM data’ section, click and drag the TIFF into ArcGIS.
  2. Select the ‘Jet’ colour scheme (or which ever colour scheme you desire) for the bathymetry
  3. Make a hillshade of the bathymetry (3D Analyst > Raster Surface > Hillshade).
  4. You can leave the ‘Azimuth’ and ‘Elevation’ values as they are, you can alter the ‘Z’ factor to increase the vertical exaggeration of the hillshade. I usually set that to 3.
  5. Right click the hillshade, select properties.
  6. In the ‘Display’ tab, you can change the transparency of the hillshade. I like to use 30% but you can play around with it if you wish.

Displaying XY data by category

Sometimes you may have the XY locations of different types of variable (i.e. waypoint, sampling location, mooring location, etc.). If you add a column into the header row of the text file named ‘Types’, for example, you can add the ‘waypoint, sampling location, decent pubs, etc.’ labels to each row and this will be the categories which you sort the XY points into within Arc.

  1. You can import all of the XY coordinates into Arc and convert to shape files using the above method.
  2. Right click the layer of XY coordinates and select the ‘Symbology’ tab, in there you can select the ‘Categories > Unique Values’, or ‘Unique Values, many fields’ if you have different variables. Using this you can add in each of the categories you put into the text file.
  3. Be sure to choose an appropriate symbol for each category by double clicking each of the ‘point’ symbols and scrolling through the symbols until you find the desired symbol.

 

Converting elevation to depth

You just made a sweet bathymetry output from Caris or Fledermaus or whatever and oh no(!), the values imply elevation (+ive) rather than depth (-ive). Here’s a sneaky tip for getting around that without having to export your data again.

 

  1. Import your elevation data using the above method.
  2. Convert to depth using the raster calculator (Spatial analyst > Raster Algebra (?) > Raster Calculator).
  3. Double click, or drag, your elevation layer so it appears in the box below the layer list.
  4. Using the calculator function to the right, multiply the elevation layer by -1.
  5. This method can be used to convert depth to elevation also.

Sneaky sneaky!

 

Taking repeat elevation (or depth) profiles from the same point.

This assumes that you have a repeat elevation or bathymetry survey ongoing.

 

  1. First you must make a line shape file (See Previous Post on how to make a shape file) and position the line in the spot you want, usually across or along the axis of a channel.
  2. Stack Profile > Right click > Batch process.
  3. It will output everything into a single table, but you can export this to an excel or text file to use in Matlab or something.
  4. I will make an additional post explaining how to code for that when I get around to figuring it out

 

Georeferencing an elevation/bathymetry

 So you received a new bathymetric survey, you go to import it into Arc and you discover that the survey has an obscure, unique georeferencing system that you can’t get a hold of. Or you have an aerial photograph that you want to georeference. What ever your issue, here’s how to fix it.

 

  1. There is a tool called ‘Georeference’.
  2. I’ve not used this tool much but I’m just letting you know it exists.

 

General tips:

If you can’t find a function, click on the search bar and type in a function name. It could be there! E.g. I needed to flip a bathymetry so I typed in ‘flip’ and both ‘flip’ and ‘mirror’ functions came up in the results. They weren’t useful for what I wanted but sometimes, life gives you digital lemons.

 

Always make sure the coordinate systems (e.g. WGS84) match up. Sometimes you have to convert the coordinate system, which is a bit of a headache, but there are websites online that can do this for you (e.g. http://www.earthpoint.us/Convert.aspx).

Basic ArcGIS for the amnesia-stricken student

As I’ve been moving through my Ph.D, the use of GIS systems has become increasingly important to process the data I have been collecting or have been given. Unfortunately for myself, it has been a while since I had to use GIS tools to their full effect with my undergraduate degree teaching me very briefly how to create a polygon in ArcGIS in my first and second year. That was 4 years ago. Since then I have barely used ArcGIS, and thus have forgotten the basics. In my frustration I turned to Google, that all-knowing tool of wisdom, to help me through my state of amnesia but was left feeling more frustrated than when I started. The lack of illustrated guides became apparent, even on Youtube, as I sifted through the seemingly endless buttons and sub-menus that ArcGIS had to offer.

So, to save time for others who may have indeed forgotten how to create a polygon, I have put together a short guide on how to not only create a polygon, but also how to cut or clip the data within that polygon. Ideally at this point you will know how to import your basemap and/or generate a difference map (Search for Map Algebra for generating difference maps, you literally subtract one map from the other). If not, let me know, and I can direct you to an appropriate guide or even write another illustrated guide myself.

How to draw polygons in ArcMap

  • In Arc Catalog, select an appropriate folder in which to work. Make sure the basemap or difference map and the polygon are in the same folder.
  • Right click on your chosen work folder and select “New > Shapefile > Polygon” and give it a sensible name (Below).

NewShapeFile

CreateNewShapeFile

 

 

  • Open the Editor toolbar if it is not already open and click “Start Editing”. On the right-hand side of the Editor toolbar is a “Create Features” button (Below).

CreateFeatures

  • When you click it, it opens a window that allows you to select the feature you wish to draw onto your map (Below).

CreatingPolygon

 

 

 

 

 

 

 

 

 

 

  • Then it is simply a case of tracing the outer edge of the feature you wish to retain. In this case, a turbidity current event on a difference map. Once you are done, double click to finish the polygon. The edge of the polygon will turn blue to indicate that this was successful.

 

 

FinishPolygon

 

 

 

 

 

 

 

 

 

  • Now to clip the data so that you’re left with the turbidity current event, and not the whole difference map. In Arc Toolbox you need to find Datamanagement Tools > Raster > Raster Processing > Clip (Below).

 

DataClip

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

  • In this menu you will have to select the input raster (the difference map or base map that you’re using) and the output extent, which is the polygon you just created. Check the “Use Input Features for Clipping Geometry” and “Maintain Clipping Extent” boxes and click “Ok”.

ClipMenu

 

 

 

 

 

 

 

 

 

  • Once this is done, you will be left with the clipped difference/base map which you can export as a TIFF for use in Matlab, or carry on working on in ArcGIS. Good job!

 

 

 

CompletedClip

If you have any questions at all, please do not hesitate to leave a comment or conact me.