Skip to main content

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 \space 1 \space = \space (Band5 \space / \space Band6) \space * \space 100\)

\(Band \space 2 \space = \space (Band5 \space / \space Band1) \space *\space 100\)

\(Band \space 3 \space = \space (Band3 \space / \space Band4) \space * \space (Band5 \space / \space Band4) \space * \space 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})

Parameter Explanation Data 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 Type Explanation

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)