This method assumes that you have already input two sets of contacts in the field area into a GIS as polylines. This method should work on any GIS, including R, but some steps are specific to Arc. In this example, Contact A and Contact B are a set of contacts for the same related lithological area. They are being compared using the shortest distance between contacts.
- Select only Contact A
- Generate a Euclidian Distance raster for Contact A. Be thoughtful about your cell size (I have used 0.5 m) and extent of the new raster.
- If using Arc, enter Editor mode and create a point layer. Ensure that point layer and the contact layer are in editing mode.
- Select only Contact B, merge the contact lines (perhaps in a new feature class).
- Generate an interval, your choice, of points that coincide with the merged Contact B. In ArcMap 9.3., you can use the "Divide" tool under the Editor menu (... maybe). In ArcMap 10,you can use the "Construct Points" tool under the Editor menu. Save and stop editing.
- In ArcMap, use Spatial Analyst's Sample tool. Use the Euclidian distance raster as your input raster and your line coinciding points as your point features ("Input Location Raster or Point Features"). Define your output table's location. The default resampling technique is fine. See above figure.
- Press OK and you will now have a table with the shortest distance from the points along Contact B to the lines of Contact A
- Now import the distance data into your favorite statistics software and see how your data compares (or just use Arc's built in statistics tools).
|Output table with sampled distances between Contact A and Contact B.|