topic Re: NumPyArrayToFeatureClass (arcpy.da) in Python script in Python Questions

NumPyArrayToFeatureClass (arcpy.da) in Python script

AlastairOram1 — Fri, 29 Aug 2014 08:31:09 GMT

Hi

I'm using arcmap 10.2 and have been trying in vain to get the
numpyArrayToFeatureClass tool to work in a python script. It generates a shapefile
but the output I get only ever has an index column - no other columns are
output.

When I use the code example here ArcGIS Help (10.2, 10.2.1, and 10.2.2) I get exactly the same problem. Can you please
give me any indication where I am going wrong?

Thanks

Alastair

Re: NumPyArrayToFeatureClass (arcpy.da) in Python script

DanPatterson_Retired — Fri, 29 Aug 2014 12:01:46 GMT

Can you post your code

Re: NumPyArrayToFeatureClass (arcpy.da) in Python script

DanielGarcia — Fri, 29 Aug 2014 12:10:58 GMT

Hi Alastair,

This code works fine for me. Please, use it and tell us what you get.

import arcpy
import numpy

outFC = "E:/temp/testShape.shp"
array = numpy.array([(1, 'testText1', (471316.3835861763, 5000448.782036674)),
                     (2, 'testText2', (470402.49348005146, 5000049.216449278))],
                    numpy.dtype([('idfield',numpy.int32),('textField','<U10'),('XY', '<f8', 2)]))

arcpy.da.NumPyArrayToFeatureClass(array, outFC, ['XY'])

Regards

Dani

Re: NumPyArrayToFeatureClass (arcpy.da) in Python script

AlastairOram1 — Fri, 29 Aug 2014 12:58:24 GMT

Well a shape file gets created but only has an index column - no other fields - see below.

idfield

Thanks

Re: NumPyArrayToFeatureClass (arcpy.da) in Python script

DanPatterson_Retired — Fri, 29 Aug 2014 12:58:31 GMT

Dani...line 7 is truncated

Re: NumPyArrayToFeatureClass (arcpy.da) in Python script

AlastairOram1 — Fri, 29 Aug 2014 13:26:55 GMT

Note that I’m starting with a pandas data frame (df) and converting it using “to_records” (although I have tried several methods like 'df.as_matrix()';

turbinesNPArray = df.to_records()

dt = turbinesNPArray.dtype

#(Note this returns:

#[('index', '<i8'), ('TurbineID', '<f8'), ('X', '<f8'), ('Y', '<f8'), ('Z', '<f8'), ('HubH', '<f8'), ('RotorD', '<f8'), ('WGS84Lat', '<f8'), ('WGS84Long', '<f8'), ('OFFSETA', '<f8'), ('OFFSETB', '<f8')]

#type turbinesNPArray:<class 'numpy.core.records.recarray'>

outShapeFile = r"%s\%s.shp" %(rootPath + "Outputs", "outLayout_temp")

arcpy.da.NumPyArrayToFeatureClass(turbinesNPArray, outShapeFile, ("TurbineID","X","Y","Z","HubH","RotorD","WGS84Lat","WGS84Long","OFFSETA", "OFFSETB"), self.spRef)

Re: NumPyArrayToFeatureClass (arcpy.da) in Python script

AlastairOram1 — Fri, 29 Aug 2014 13:27:55 GMT

Not sure what you mean about line 7?

Re: NumPyArrayToFeatureClass (arcpy.da) in Python script

DanielGarcia — Mon, 01 Sep 2014 09:36:14 GMT

I attach you py file and the result shapefile.

These work for me in ArcGIS 10.2.2.3552 version and Python 2.7.

Re: NumPyArrayToFeatureClass (arcpy.da) in Python script

AlastairOram1 — Tue, 02 Sep 2014 11:20:05 GMT

I get;

Traceback (most recent call last):

File "<module1>", line 23, in <module>

RuntimeError: create table

Re: NumPyArrayToFeatureClass (arcpy.da) in Python script

JamesCrandall — Thu, 04 Sep 2014 11:40:58 GMT

I suspect that you are not correctly setting up the nparray:

turbinesNPArray = df.to_records()

Without knowing your fields info, I don't see you fully qualifying the fields in your conversion from Pandas DataFrame to NumPy Array. You can try the following but I have no real way to test:


turbinesNPArray = numpy.array(df.to_records(), numpy.dtype([('index', '<i8'), ('TurbineID', '<f8'), ('X', '<f8'), ('Y', '<f8'), ('Z', '<f8'), ('HubH', '<f8'), ('RotorD', '<f8'), ('WGS84Lat', '<f8'), ('WGS84Long', '<f8'), ('OFFSETA', '<f8'), ('OFFSETB', '<f8')]))