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 |
Data Type | Explanation |
Raster | The output raster with Sultan's formula applied. |
Code sample
Calculates Sultan's formula for a Landsat ETM image.
import arcpy
Sultan_raster = arcpy.sa.Sultan("LandsatETM.tif",1,3,4,5,6)