nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

NumPyArrayToFeatureClass gives RuntimeError: create table

GrzegorzJeleniewicz

I've tried to change polygon feautre to point by using tools:

FeatureClassToNumPyArray
NumPyArrayToFeatureClass.
My code as below:
import arcpy
import numpy as np
in_file='C:/ios/OneDrive/siatki_obliczeniowe/siatki.gdb/siatka_005'
a = arcpy.da.FeatureClassToNumPyArray(in_file, ["*"], explode_to_points=False)
out_file='C:/ios/OneDrive/siatki_obliczeniowe/siatki.gdb/xxxx3'

arcpy.da.NumPyArrayToFeatureClass(a, out_file,["shape"])

I got an error:

Traceback (most recent call last):
File "C:/ios/OneDrive/python/feature_to_point.py", line 7, in <module>
arcpy.da.NumPyArrayToFeatureClass(a, out_file,["shape"])
RuntimeError: create table

I had changed out_file to shp file and it works but it is not solution for me. I would like to be able write to .gdb.
Any idea how to solve this problem?

Find more posts tagged with

Accepted answers

GrzegorzJeleniewicz

I got answer:

arcgis desktop - NumPyArrayToFeatureClass gives RuntimeError: create table - Geographic Information Systems Stack Exchan…

All comments

DanPatterson_Retired

for future reference

all_flds <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>f<SPAN class="punctuation token">.</SPAN>name <SPAN class="keyword token">for</SPAN> f <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListFields<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"*"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN>

editable_flds <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>f<SPAN class="punctuation token">.</SPAN>name <SPAN class="keyword token">for</SPAN> f <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListFields<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"*"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">if</SPAN> f<SPAN class="punctuation token">.</SPAN>editable<SPAN class="punctuation token">]</SPAN>

all_flds

<SPAN class="punctuation token">[</SPAN><SPAN class="string token">'OBJECTID'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Shape'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Grid_codes'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'F_'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Shape_Length'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Shape_Area'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Rand_0'</SPAN><SPAN class="punctuation token">,</SPAN>
 <SPAN class="string token">'Rand_1'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Shape_Length_1'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Shape_Area_1'</SPAN><SPAN class="punctuation token">]</SPAN>

editable_flds

<SPAN class="punctuation token">[</SPAN><SPAN class="string token">'Shape'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Grid_codes'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'F_'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Rand_0'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Rand_1'</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

OBJECTID, and Shape info can be accessed through their @ property, other variants of Shape aren't editable and removing them from the list of candidates is a safe bet... at least for a file gdb

GrzegorzJeleniewicz

I changed my code following your instructions:

import arcpy
import numpy as np
in_file='C:/ios/OneDrive/siatki_obliczeniowe/siatki.gdb/strefy_miejskie_siatka_01'
flds = ['OID@', 'SHAPE@X', 'SHAPE@Y']
other_flds = [i.name for i in arcpy.ListFields(in_file) if i.type not in ('OID', 'Geometry')]
flds.extend(other_flds)
a = arcpy.da.FeatureClassToNumPyArray(in_file, flds, explode_to_points=False)
out_file='C:/ios/OneDrive/siatki_obliczeniowe/siatki.gdb/axxxyz9'
print (a.dtype)
arcpy.da.NumPyArrayToFeatureClass(a, out_file,['SHAPE@X', 'SHAPE@Y'])

Unfortunately I got error again:

Traceback (most recent call last):
File "C:/ios/OneDrive/python/feature_to_point.py", line 10, in <module>
arcpy.da.NumPyArrayToFeatureClass(a, out_file,['SHAPE@X', 'SHAPE@Y'])
RuntimeError: create table

The most annoying fact is that if I change the 'out_file ' to shp file instead of feature class in geodatabase (.gdb) it works (with both code versions). For sure there is a problem with saving output to geodatabase (.gdb). I have no clue why.

DanPatterson_Retired

hmmm if you only want the centroid coordinate, you should avoid using the true shape, but use one of the alternates and specify your field names explictly

FeatureClassToNumPyArray—Data Access module | ArcGIS Desktop

you can easily parse the fields rather than using "*" which is ok in some situations, but the shape field should be parsed into x and y coordinates usually

flds = ['OID@', 'SHAPE@X', 'SHAPE@Y']

other_flds = [i.name for i in arcpy.ListFields(fc) if i.type not in ('OID', 'Geometry')]

flds.extend(other_flds)

flds

['OID@', 'SHAPE@X', 'SHAPE@Y', 'ID', 'Norm', 'Rank1', 'Text01', 'Text02', 'Xs', 'Ys',
 'Int_short', 'Float_fld', 'Grid_code', 'Rank_shp', 'PolyID', 'Grid_Count',
 'Diff_from_mean', 'x']<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

GrzegorzJeleniewicz

Thank you for your answer.

my geomtery field is ['Shape']. I've cheked it with a.dtype and got result:

[('OBJECTID', '<i4'), ('Shape', '<f8', (2,)), ('Shape_Length', '<f8'), ('Shape_Area', '<f8'), ('miasto', '<U50')]

I've changed small 's' for capital 'S', but it didn't help

Yes, I want the centroid of the original polygon, so explode_to_points=False is ok for me.

The output featureclass name is set to new each run.

DanPatterson_Retired

your geometry field... what was it named? and what was it in the array?

you indicate it was ['shape']? or ['Shape@'] or 'SHAPE@' or some other variant

use

a<SPAN class="punctuation token">.</SPAN>dtype‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN></SPAN>

and get the exact case and spelling of the geometry field.

Also you will need to set

<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>explode_to_points<SPAN class="operator token">=</SPAN><SPAN class="token boolean">True</SPAN><SPAN class="punctuation token">)</SPAN>  <SPAN class="comment token"># in order to get points, ‍</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN></SPAN>

you can't set explode_to_points=False, unless you want the centroid of the original polygon

Thirdly, if the output featureclass name exists, it will fail

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.