ISIS 3 Application Documentation
|Printer Friendly View | TOC | Home|
Create slope, aspect, or percent slope data cube.
This program creates a new cube containing slope (radians, degree,s or percent) or aspect (radians, degrees) data from an input cube. This program computes the slope and aspect of an input cube and outputs the results. An individual slope, aspec, or percent slope cube can be rendered using this program.
Data requirements to run slpmap:
A B C D E F G H IEach output pixel is computed as follows:
SLOPE EQUATION [dz/dx] = ((C + 2F + I) - (A + 2D + G)) / (8 * X_PIXEL_RESOLUTION) [dz/dy] = ((G + 2H + I) - (A + 2B + C)) / (8 * Y_PIXEL_RESOLUTION) slope = ATAN ( SQRT ( [dz/dx]^2 + [dz/dy]^2 ) ) percentslope = slope / 90The slope equation above assumes the pixels are not square: hence X_PIXEL_RESOLUTION and Y_PIXEL_RESOLUTION. That is, the x distance across the pixel is not equal to the y distance form top to bottom of the pixel.
ASPECT EQUATION [dz/dx] = ((C + 2F + I) - (A + 2D + G)) / 8 [dz/dy] = ((G + 2H + I) - (A + 2B + C)) / 8 aspect = 90 - ATAN2 ([dz/dy], -[dz/dx]) if (aspect < 0) then aspect = aspect + 360
Slope: Slope is typically between 0 and 90 degrees, were 0 is flat and 90 is vertical.Slope may be output in radians, degrees, or percent slope.
Aspect: The aspect represents the direction or 0 to 360 degrees of the slope in pixel space (see figure below). From the center pixel (E), 0 degrees is straight towards B (generally north), 45 is towards C, 90 is towards F, 135 is towards I, 180 is towards H, and so on.
A B C D E F G H I
Scale Construction: By default, the program will attempt to remove any scaling differences by using the map projection information (PIXRES=AUTOMATIC). This computation is done at every pixel so the correct x-to-y ratio is computed. This is important, for example, for global maps where the x/y ratio deviates with distance from the latitude and/or longitude of true scale in a map projection. This only works for radius DTMs.
What if the input cube Z units are in elevation values instead of radius values?
The default setting PIXRES=AUTOMATIC will not work. You will need to set PIXRES=FILE (recommended for a map projected file) or PIXRES=USER. Note: unlike PIXRES=AUTOMATIC, these later two methods do not correct for distortions int he map projection. They will return good results for smaller regious when the map projection is defined to minimixe distortions (e.g. LROC NAC or HIRISE sereo DTMs).
What if the input cube is not map projected?
If the image lacks a map projection,you must provide the pixel resolution via PIXRES=USER (for this option the pixels are assumed to be square). You must also provide a single value using the RESOLUTION parameter that will be applied to all pixels in the image and in both directions.
What if the xy units are not the same as the z units?
The program assumes the xy units are the same as the z (pixel) units. By default the program assumes the units are the same, but allows you to scale the z units to the xy units using the CONVERSIONFACTOR parameter if the PIXRES=USER option is selected.
Use this parameter to select the input cube.The algorithm is applied to all bands of a multiband cube.
Use this parameter to define the filename of the resultant slope, aspect, or percent slope cube.
This specifies the type of output pixels that will be created. The output file will contain either slope, aspect, or, percent slope.
Units of the output pixels
This determines how the resolution of the input pixels will be defined.
This parameter is a multiplicative factor for the resolution. It is used to convert the spatial units (pixels) to the correct distance units. Since the slope algorithm depends on cancelling out the units, the conversion factor from the z vertical (input pixel) units to the xy horizontal (spatial) units needs to be considered. For example, if the xy units are in meters and the z units are in kilometers the unit conversion factor should be 1000 as we are converting z kilometers to xy meters. The default presumes the units are the same.
Defines both the X and Y pixel resolution to be used across the entire image. Note that if the pixel height varies significantly from the width, the slope will not be computed correctly.