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 (arcpy.da) in Python script

AlastairOram1

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

Find more posts tagged with

Accepted answers

All comments

JamesCrandall

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')]))

DanielGarcia

I attach you py file and the result shapefile.

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

test.py.zip

testShape.zip

AlastairOram1

Not sure what you mean about line 7?

AlastairOram1

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)

DanPatterson_Retired

Dani...line 7 is truncated

AlastairOram1

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

idfield

Thanks

DanielGarcia

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

DanPatterson_Retired

Can you post your code

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.