topic Re: Batch Field Extraction from a Point Feature Class: in ArcGIS Spatial Analyst Questions

Batch Field Extraction from a Point Feature Class:

MarkEnglish — Wed, 20 Aug 2014 11:42:54 GMT

Hello: I have a point dataset which has 20+ fields indicating different storm surge elevations. For instance, one field indicates storm 1, and contains 200 water surface elevation records. The second field indicates storm 2, with also 200 water surface records, and so on. The returns for each storm are present on the same record/point, just different values. I want to generate 20+ point feature classes, one for each storm by extracting each field into its own dataset. Following this process, I'd like to run a batch process for each storm point file that generates a raster for each storm scenario. I know how to do this manually, but it would take a lot of time considering the number of steps involved. Any help would be greatly appreciated.

Thank you!

-Mark

Re: Batch Field Extraction from a Point Feature Class:

JakeSkinner — Wed, 20 Aug 2014 12:11:25 GMT

Hi Mark,

You wouldn't have to create a separate feature class for each field. You could iterate through each field and create a raster using the IDW tool. Below is an example that iterates through a point feature class called 'pt' and creates a File Geodatabase Raster Dataset using the IDW tool.

import arcpy
arcpy.CheckOutExtension("spatial")

from arcpy import env
from arcpy.sa import *

env.workspace = r"C:\temp\python\test.gdb"

fc = "pt"

for field in arcpy.ListFields(fc):
    if field.name not in ('OBJECTID', "SHAPE"):
        outIDW = Idw("pt", field.name, "2", "2", "VARIABLE 12")
        outIDW.save("pt_" + field.name)

Re: Batch Field Extraction from a Point Feature Class:

JakeSkinner — Mon, 25 Aug 2014 14:54:57 GMT

You can specify the cell size (i.e. 10 Meters) by updating the 3rd parameter in the Idw function. Ex:

outIDW = Idw("pt", field.name, "10", "2", "VARIABLE 12")

You can see a description of each parameter from the below link:

ArcGIS Help (10.2, 10.2.1, and 10.2.2)

Re: Batch Field Extraction from a Point Feature Class:

MarkEnglish — Mon, 25 Aug 2014 16:18:47 GMT

Thanks Jake!

It worked, although I am unearthing additional issues within the data. For instance, 'no data' is represented as -999. Is there a method to convert all -999 in the point storm surge dataset to NULL?

M-

Re: Batch Field Extraction from a Point Feature Class:

JakeSkinner — Mon, 25 Aug 2014 17:24:45 GMT

Yes, you can use the Set Null tool to do this.

Re: Batch Field Extraction from a Point Feature Class:

MarkEnglish — Mon, 25 Aug 2014 17:50:47 GMT

Yes, but that would have to be performed on all of the raster outputs. Is there a way to locate all -999 records within a single point dataset and change them to NULL? Or batch each field at a time?