My last blog article described how to create a slope map from a digital elevation model in Surfer. Moving forward on that topic, I found this blog article written for ArcMap on creating aspect-slope maps, which was improved upon for QGIS. This single map combines both the compass direction of slopes (aspect) and the steepness of the slopes (in degrees) and uses both color and saturation to display the combined results. Slopes facing different directions use different colors, and the brightness of that color shows the steepness of that slope (the brighter the color, the steeper the slope). I thought this was a really interesting map type and it made me wonder how this could be done in Surfer.
Coincidentally, at that time, a user asked me this exact question! The user wanted to come up with a way to see the very small slope variations in the soft sediments they have on the surface, using both aspect and slope. Looking at the slopes and aspect together may reveal small variations that otherwise could be overlooked.
Before we delve into the steps on how to create this type of map, I want to go over a little background on what the colors represent.
To create an aspect-slope map in Surfer, all you need to start with is a grid, DEM, or DTM file. Using the method described in the QGIS article, the idea is to create a grid of slopes (in %) and create grid of aspect directions. Reclassify the slope grid into bins from 0 to 8 in steps of 2, and reclassify the aspect grids into bins from 10 to 80 in steps of 10.
|Original Slope% Z Value||New Z Value|
|>=0.0 and <5.0||0|
|>=5.0 and <15.0||2|
|>=15.0 and <30.0||4|
|>=30.0 and <45.0||6|
|Original Aspect Z Value||New Z Value|
|>=0 and <22.5||10|
|>=22.5 and <67.5||20|
|>=67.5 and <112.5||30|
|>=112.5 and < 157.5||40|
|>=157.5 and <202.5||50|
|>=202.5 and <247.5||60|
|>=247.5 and <292.5||70|
|>=292.5 and <337.5||80|
|>=337.5 and <360.5||10|
Reclassify the slope (%) and aspect grid files so that the slope values are even numbers from 0 to 8, and the aspect values are ten values from 10 to 80.
Combine the reclassified grids by adding them together to create a single aspect-slope grid. For the combined grid, the values can range anywhere from 10 to 88, which is the minimum of the slope plus the minimum of aspect (10+0) and the maximum of the slope plus the maximum of the aspect (80+8). The first digit in the number in the ten’s place is the aspect orientation, and the second digit in the one’s place is the slope. For example, a value of 24 in the combined grid indicate a slope in the direction between 22.5° and 67.5° azimuth (since the first digit is 2), and the slope would be between 15% and 30% (since the second digit is 4).
We will then create a map and color it based on this combined value. The actual color is based on the aspect value (the first digit in the ten’s place) and the brightness of that color based on the slope value (the second digit in the one’s place). Since any value that has a 0 in the one’s place (e.g. 10, 20, 30, etc.) is relatively flat (a slope between 0-5°), we can assign it a flat gray color that we can make completely transparent. For the other colors, we can use colors based on a wheel like the image below. For example, a value of 24 in the combined grid would be assigned a medium green color.
This color wheel was selected from the ArcMap article. The color around the wheel is based on the aspect (slope direction), and is identified by the first digit in the 10’s place of the Z values in the combined grid. The brightness of that color is based on the slope itself, identified by the second digit in the 1’s place of the Z values in the combined grid.
Performing the steps manually is not very difficult, but they do take some time. I decided to shorten this workflow significantly by writing a script, compatible with Surfer 13. The manual process of walking through the steps only takes about 5 minutes, but the script takes only seconds and is wonderfully easy. I’ll provide instructions for both running the script and manually walking through the steps.
Steps to Create an Aspect-Slope Map using a Script
To run the script to create the aspect-slope map in seconds, follow these steps:
That’s it! Surfer is opened and the map is created with the appropriate values and colors. Even with a large grid file, this takes only seconds on my computer.
Steps to Create an Aspect-Slope Map Manually
If you want to work through this manually, the equivalent steps to perform in Surfer are as follows:
Capturing both slope and aspect in a single map is an excellent way to identify large and small scale trends in a map. You have tremendous flexibility in Surfer. You can adjust the classifications if you wish, to emphasize certain slope or aspect ranges by adjusting the Grid Math functions, or you can choose to make the flat areas solid gray instead of transparent by setting the Opacity in the Colormap dialog to 100%. There are many options you can choose from so that this map shows exactly what you need it to.