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})

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.ia.Sultan("LandsatETM.tif",1,3,4,5,6)