¿Cómo obtener el valor de salida de Feature Compare (esquema) y añadirlo al csv existente?

Question

Estoy trabajando en un proyecto de integración de datos y estoy tratando de identificar las clases de características que tienen un esquema que varía de dos clases de características base (una para las características de puntos y otra para las de polilíneas).

He jugado con la herramienta Feature Compare para comparar esquemas, pero me resulta tedioso utilizar esta herramienta, incluso en modo batch, porque genera un único archivo de texto con unos pocos valores para cada comparación realizada .

Empecé a escribir un script en Python para recorrer una lista de características en una lista y comparar su esquema basado en el tipo de geometría. Sin embargo, el problema de integrar la salida en un solo archivo de texto persiste. ¿Hay alguna manera de obtener este valor y anexarlo a un archivo de texto existente?

El código (desordenado) hasta ahora es el siguiente:

# Import arcpy module
import arcpy

# Import environments; set workspace
from arcpy import env

# Set workspace
env.workspace = "Data_Comparison.gdb"

# List features to be compared
fclist = arcpy.ListFeatureClasses()
print fclist

# Loop through list to separate point and line feature classes
for fc in fclist:
    desc = arcpy.Describe(fc)
    name = desc.file
    geo = desc.shapeType

    print ("Feature name: " + name)
    print ("Feature geometry: " + geo)

    if geo == "Point":
        arcpy.FeatureCompare_management("swr_net_nodal", fc, "ENTID", "SCHEMA_ONLY", "", "", "", "", "", "", "CONTINUE_COMPARE", "fc + '_compare'")
        output = arcpy.GetMessages()
        print output

Answer 1

Una posible solución es crear un archivo de texto con antelación y leer el archivo compare_file línea por línea y añadirlo al archivo de texto.

    # import arcpy module
    import arcpy

    # import os module
    import os

    # set workspace
    workspace = r"some\path"
    arcpy.env.workspace = workspace

    # prepare compare file
    out_compare_file = r"some\path\to\file.txt"

    # list of features to compare
    fclist = arcpy.ListFeatureClasses()

    # loop through all feature class in workspace
    for fc in fclist:
        fc_desc = arcpy.Describe(fc)
        # getting name without extention
        fc_name = os.path.splitext(fc_desc.file)[0]
        fc_shape_type = fc_desc.shapeType

        if fc_shape_type == 'Point':
            # Set variables
            base_features = r"some\path\to\base_feature"
            test_features = fc
            sort_field = 'some_field'
            compare_type = "SCHEMA_ONLY"
            ignore_option = "#"
            xy_tolerance = "#"
            m_tolerance = "#"
            z_tolerance = "#"
            attribute_tolerance = "#"
            omit_field = "#"
            continue_compare = "CONTINUE_COMPARE"
            compare_file = workspace + "\\" + fc_name + "_compare.txt"
            try:
                compare_result = arcpy.FeatureCompare_management(
                    base_features,
                    test_features,
                    sort_field,
                    compare_type,
                    ignore_option,
                    xy_tolerance,
                    m_tolerance,
                    z_tolerance,
                    attribute_tolerance,
                    omit_field,
                    continue_compare,
                    compare_file
                )
                print compare_result.getOutpuput(1)

                # open prepared text file to append
                with open(out_compare_file, 'a') as outfile:
                    # open compare file to read
                    with open(compare_file) as infile:
                        for line in infile:
                            outfile.write(infile)
                        # do not forget to close file handlers
                        infile.close()
                    outfile.close()                           
            except:
                print arcpy.GetMessages()