If you don't have any data already preserved in the output feature class, it might be faster to insert new rows instead of updating existing ones. For this, you would use da.InsertCursor.
If you still would like to update the rows, then do the following:
#23 data columns to be added
fieldlist = ["azimuth", "orient", "crit#f", "freq1", "amp1", "pwr%1", "freq2", "amp2", "pwr%2", "freq3", "amp3", "pwr%3", "freq4", "amp4", "pwr%4", "freq5", "amp5", "pwr%5", "Nspurs", "len_mean", "len_sd", "ht_mean", "ht_sd"]
for fieldname in fieldlist:
arcpy.AddField_management(gridcopy, fieldname, "FLOAT") #Add new fields
with arcpy.da.UpdateCursor(gridcopy, fieldlist) as cursor: #Add data
for row in cursor:
row = X
cursor.updateRow(row)
~50K rows is nothing, will be updated in no time at all. Since you are working with numpy arrays, it might also be worth to check arcpy.da methods for converting numpy arrays into feature classes such as NumPyArrayToFeatureClass which will convert a NumPy structured array to a point feature class. They are very fast, I've been using arcpy.da module functions for large numpy arrays (millions of records) and they are very fast.
Performance information update:
It will be faster to use the UpdateCursor than to create geodatabase tables first from numpy arrays as creating new objects in a geodatabase does have an overhead and using ArcGIS-based joins would be very slow. Look at the execution time results - updating actual rows takes just 2 seconds; the slowest part is importing arcpy and getting count of features:
from functools import wraps
import numpy as np
import time
fields_list = ['field{}'.format(i) for i in xrange(1,24)]
fc = r'C:\GIS\Temp\ArcGISHomeFolder\Default.gdb\points'
#----------------------------------------------------------------------
def report_time(func):
'''Decorator reporting the execution time'''
@wraps(func)
def wrapper(*args, **kwargs):
start = time.time()
result = func(*args, **kwargs)
end = time.time()
print(func.__name__, round(end-start, 3))
return result
return wrapper
@report_time
#----------------------------------------------------------------------
def setup():
"""set up sample data and add fields to feature class"""
import arcpy
total_features = int(arcpy.GetCount_management(fc).getOutput(0))
print total_features
#for f in fields_list:
#arcpy.AddField_management(in_table=fc, field_name=f, field_type='FLOAT')
data = np.random.uniform(low=0.5, high=13.3, size=(total_features, 23))
return data
@report_time
#----------------------------------------------------------------------
def update_cursor(data):
"""update feature class fields using da.UpdateCursor"""
with arcpy.da.UpdateCursor(fc, fields_list + ['OID']) as cursor: #Add data
for row in cursor:
row[0:23] = data[row[-1] - 1]
cursor.updateRow(row)
#----------------------------------------------------------------------
def main():
""""""
data = setup()
update_cursor(data)
if __name__ == '__main__':
main()
49696 ('setup', 20.678) ('update_cursor', 1.837)
