Экспорт таблицы в XYZ ASCII файл через ArcPy?

23

Я ищу способ экспортировать таблицу ArcGIS (созданную с помощью инструмента Sample ) в текстовый файл через ArcPy.

Я могу сделать это в ArcGIS через контекстное меню, щелкнув правой кнопкой мыши по таблице, но не нашел способа написать это.

Toni
источник

Ответы:

31

Вы можете сделать это с помощью курсора, чтобы получить данные из таблицы и записать в текстовый файл с разделителями-запятыми.

РЕДАКТИРОВАТЬ: я добавляю более краткий блок кода для выполнения задачи с использованием csvмодуля Python

Новый ответ с использованием курсора arcpy.da:

import arcpy,csv

table =r'c:\path\to\table'
outfile = r'c:\path\to\output\ascii\text\file'

#--first lets make a list of all of the fields in the table
fields = arcpy.ListFields(table)
field_names = [field.name for field in fields]

with open(outfile,'wb') as f:
    dw = csv.DictWriter(f,field_names)
    #--write all field names to the output file
    dw.writeheader()

    #--now we make the search cursor that will iterate through the rows of the table
    with arcpy.da.SearchCursor(table,field_names) as cursor:
        for row in cursor:
            dw.writerow(dict(zip(field_names,row)))

Новый ответ, используя курсор в старом стиле:

import arcpy,csv

table =r'c:\path\to\table'
outfile = r'c:\path\to\output\ascii\text\file'      

#--first lets make a list of all of the fields in the table
fields = arcpy.ListFields(table)
field_names = [field.name for field in fields]

with open(outfile,'wb') as f:
    w = csv.writer(f)
    #--write all field names to the output file
    w.writerow(field_names)

    #--now we make the search cursor that will iterate through the rows of the table
    for row in arcpy.SearchCursor(table):
        field_vals = [row.getValue(field.name) for field in fields]
        w.writerow(field_vals)
    del row

Старый ответ:

import arcpy

table =r'c:\path\to\table'
outfile = r'c:\path\to\output\ascii\text\file'


#--first lets make a list of all of the fields in the table
fields = arcpy.ListFields(table)

i = 1
f = open(outfile,'w')
for field in fields:
    #--write all field names to the output file
    if i < len(fields):
        f.write('%s,' % field.name)
        i += 1
    else:
        f.write('%s\n' % field.name)

#--now we make the search cursor that will iterate through the rows of the table
rows = arcpy.SearchCursor(table)
for row in rows:
    i = 1
    for field in fields:
        if i < len(fields):
            f.write('%s,' % row.getValue(field.name))
            i += 1
        else:
            f.write('%s\n' % row.getValue(field.name))
del rows
f.close()
Джейсон
источник
Рад, что я мог помочь вам @Toni
Джейсон
1
@ Джейсон - Спасибо, это было очень полезно. Я новичок, поэтому у меня нет репутации, чтобы комментировать ваш принятый ответ. Я думаю, что есть небольшая ошибка в новом ответе, который использует курсор arcpy.da. with arcpy.da.SearchCursor(table) as cursor:должно бытьwith arcpy.da.SearchCursor(table, field_names) as cursor:
Хороший улов @TylerG, я отредактировал ответ, включив в него список полей, требуемых курсором доступа к данным. Спасибо.
Джейсон
8

Вы можете захотеть «Экспортировать атрибут объекта в ASCII», хитро названный arcpy.ExportXYv_stats

http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//005p0000003v000000

import arcpy

feature = "path to feature here"
# fieldnames must be explicitly provided. Note that you will get additional fields based on the feature type (e.g., "XCoord" and "YCoord" for point features)
fieldnames = [X.name for X in arcpy.ListFields(feature)]
# delimiter options "SPACE", "COMMA", or "SEMI-COLON"
# header options "ADD_FIELD_NAMES" or "NO_FIELD_NAMES"
arcpy.ExportXYv_stats(feature, fieldnames, "SPACE", "path to outfile", "ADD_FIELD_NAMES")
Майк Ричардс
источник
+1 за убийство! Это работает интерактивно, но не так хорошо в модели или скрипте, потому что имена полей должны быть указаны.
Мэтт Вилки
1

Вот кусок кода, который я использую. Это помогает мне генерировать все мои выходные файлы в файл .txt с диапазоном от 0,100. Надеюсь, это поможет

for x in xrange(0,100):
    if os.path.isfile(outfolder + "/" + "outputs" + str(x) +".shp" ):
       inFeatures = "selected_features" + str(x) +".shp"
       export_ASCII = "ASCII " + str(x) +".txt"
       arcpy.ExportXYv_stats(inFeatures, ["Cur1_pr2","Cur3_pl1","slp1"],"SPACE", export_ASCII,"ADD_FIELD_NAMES")
MAXg
источник