TPI Tool Bathymetric High

This tool maps bathymetric high features from bathymetry data using a Topographic Position Index (TPI) [Weiss, 2001] based method.

Positive TPI usually indicates bathymetric high location.

The followings are the key steps of this tool:

  1. Calculate TPI from the input bathymetry raster using the TPI Circle Radius parameter

  2. Calculate the TPI threshold using equation (1), where c is the TPI STD Scale parameter, mean_TPI and STD_TPI are the mean and standard deviation statistics of the TPI raster

    (1)\[TPI\_threshold = mean\_TPI + c * STD\_TPI\]

  3. Select locations that have TPI values greater than the TPI threshold

  4. Convert the selected areas into polygons

  5. Remove the polygons with areas smaller than the Area Threshold parameter to obtain the final set of bathymetric high features as output

The TPI radius should be large enough to capture the largest bathymetric high features in the dataset. For example, for a 5m resolution bathymetry raster, a radius of 50 cells should be used to capture any bathymetric high features that is smaller than 500 m in length. Users should also experiment the TPI STD Scale and the Area Threshold parameters to obtain an optimal output solution.

_images/TPI.png 
 1from arcpy import env
 2from arcpy.sa import *
 3arcpy.CheckOutExtension("Spatial")
 4
 5# import the python toolbox
 6arcpy.ImportToolbox("C:/semi_automation_tools/User_Guide/Tools/BathymetricHigh.pyt")
 7
 8env.workspace = 'C:/semi_automation_tools/testSampleCode/Gifford.gdb'
 9env.overwriteOutput = True
10
11# specify input and output parameters of the tool
12inBathy = 'gifford_bathy'
13outTPI = 'gifford_tpi10'
14outFeat = 'tpi10_1std_300000m2_BH'
15areaT = '300000 SquareMeters'
16tpiRadius = 10
17tpiSTD = 1.0
18tempWorkspace = 'C:/Users/u56061/Documents/ArcGIS/Projects/UserGuide/UserGuide.gdb'
19
20# execute the tool
21arcpy.BathymetricHigh.TPITool(inBathy,outTPI,outFeat,areaT,tpiRadius,tpiSTD,tempWorkspace)