Python-arcpy : How make feature list from gdb and value list from field and use them in gaLayer function?

Question

I'm not an expert in python, however I have to write a script to do geostatistical analyses with the GACreateLayer function from Arcgis.

I have several data (>2000) and I want to use list of features and list of field values to perform cokriging.

I have written a script, the first part return list of features in my gdb and the list of values in the field that I have chosen (it works and return field values and feature names) and the second part do cokring (it works for one feature and its field)

However when I want to use both lists (inputDSet, end of script) in order to create a GA Layer it doesn't work.

File "C:\Program Files\ArcGIS\Desktop10.0\ArcPy\arcpy\ga.py", line 921, in GACreateGeostatisticalLayer     raise e
ExecuteError: Failed to execute. Parameters are not valid.
ERROR 045001: Input dataset(s) error. Table of inputs is not complete.
Failed to execute (GACreateGeostatisticalLayer)

I don't know how correctly use lists? May I transform them or something else?

Here, is my code:

import sys, string, os, arcgisscripting
# Create the Geoprocessor object
gp = arcgisscripting.create()

def listFcsInGDB():
    # list all Feature Classes in a geodatabase 
    gp.workspace = "C:\\INTERP\\PJ_class.gdb"

    fcs = []
    for fds in gp.ListDatasets('','feature') + ['']:
        for fc in gp.ListFeatureClasses('','',fds):
            #yield os.path.join(fds, fc)
            fcs.append(os.path.join(fds, fc))
    return fcs

fcs = listFcsInGDB()
for fc in fcs:
    print fc

# List of value in fields
input_dataset =  fc   
Atts = 'PJ_RacPJ'    #field with rainfall values       

rows = gp.searchcursor(fc)
row = rows.next()

NewList = []
for fc in fcs:
    for row in gp.SearchCursor(fc):
        fcValue = row.getvalue(Atts)
        NewList.append(fcValue)
#print NewList

arcpy.CheckOutExtension("GeoStats")

# Load required toolboxes
arcpy.ImportToolbox("C:/Documents and Settings/ArcGIS/Toolbox.tbx")

# Local variables:
Cokriging_xml = "C:\\INTERP\\COK_.xml"
CK = "CK2"
ids1 = "C:\\INTERP\\data.gdb\\dem1000"


for fc in fcs:
    for fcValue in NewList:
        inputDset1 = fcValue          #Variables
        inputDset2 = ids1            #Covariable
        InputDset = "inputDset1;inputDset2"

# Create cokriging layer (code from model builder)

    tempEnvironment0 = gp.autoCommit
    arcpy.env.autoCommit = "1000"
    tempEnvironment1 = gp.spatialGrid1
    arcpy.env.spatialGrid1 = "0"

    arcpy.GACreateGeostatisticalLayer_ga(Cokriging_xml,InputDset,CK) ##Function

    arcpy.env.autoCommit = tempEnvironment0
    arcpy.env.spatialGrid1 = tempEnvironment1

I need help.

If someone has an idea? It would be great.

Thanks

Answer 1

If I'm interpreting the help on this correctly, in_datasets for GACreateGeostatisticalLayer should be a comma-separated list of space-separated input dataset/input field pairs.

[inputDset1 field1, inputDset2 field2]

You also have them in quotes - InputDset = "inputDset1;inputDset2" - which could be a problem since they are variables; when in the quotes they could be getting interpreted as strings.

I haven't tested any of this, just some thoughts of things to look at.