Sultan

Summary

Calculates Sultan's formula from a six-band 8-bit raster object and returns a three-band 8-bit raster object.

Discussion

The Sultan's process takes a six-band 8-bit image and uses Sultan's formula to produce a three-band 8-bit image. The resulting image highlights rock formations called ophiolites on coastlines. This formula was designed based on the TM or ETM bands of a Landsat 5 or 7 scene. The equations applied to create each output band are as follows:

Band 1 = (Band5 / Band6) x 100
Band 2 = (Band5 / Band1) x 100
Band 3 = (Band3 / Band4) x (Band5 / Band4) x 100

For information about other multiband raster indexes, see the Band Arithmetic raster function.

The referenced raster dataset for the raster object is temporary. To make it permanent, you can call the raster object's save method.

Syntax

Sultan (raster, {band1_id}, {band3_id}, {band4_id}, {band5_id}, {band6_id})
ParameterExplanationData Type
raster

The input raster.

Raster
band1_id

The band index of Band 1.

The band ID index uses one-based indexing.

(The default value is 1)

Integer
band3_id

The band index of Band 3.

The band ID index uses one-based indexing.

(The default value is 3)

Integer
band4_id

The band index of Band 4.

The band ID index uses one-based indexing.

(The default value is 4)

Integer
band5_id

The band index of Band 5.

The band ID index uses one-based indexing.

(The default value is 5)

Integer
band6_id

The band index of Band 6.

The band ID index uses one-based indexing.

(The default value is 6)

Integer
Return Value
Data TypeExplanation
Raster

The output raster with Sultan's formula applied.

Code sample

Sultan example

Calculates Sultan's formula for a Landsat ETM image.

import arcpy

Sultan_raster = arcpy.sa.Sultan("LandsatETM.tif",1,3,4,5,6)

Related topics