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; } } }

Export selected records to CSV using python (IDLE)

TheoFaull

Hello all!

I have a polygon shapefile. It contains 200 records and a (long) field with reference numbers (type codes) 0, 1, 2 or 3.

I have a point shapefile with 25,000 records. The points fall inside the various polygons throughout the map.

I want to-

Select all points that intersect the polygons that have a type code of 0, and export this selection of points to a new CSV file.

Then repeat this process for ref number = 1, 2 and 3...

4 seperate scripts (one for each ref number) or 1 script doing it all at once: either would be fine.

I want to do all this using the IDLE python GUI rather than going into ArcMap. So I want the script to point to shapefiles on my hard drive raher than layers in the Table Of Contents.

I'm still a noob at python so all I have on my mind so far is:

# Import the required python modules (I still don't quite know what this does)

??? import arcpy ???

??? from arcpy import env ???

# set the workspace environment

??? env.workspace = "C:\Data\" ???

# Select all polygons where type code = 0 and export it to a temporary feature class

arcpy.Select_analysis('C:\Data\Polygon.shp', r'\\somernt\curo\users\faullt\ArcGIS\Default.gdb\Type_code_0_polygon', '"Type_code" = 0')

# select all points that intersect the new 'type code = 0' feature class

arcpy.SelectLayerByLocation_management('C:\Data\Points.shp', 'INTERSECT', 'Scheme- unknown or no restriction', '#', 'NEW_SELECTION', 'NOT_INVERT')

# export the new selection of points to a CSV file in C:\Data

??? ???

# display a completion message

print "CSV created"

The actual lines of python I have typed above are taken from dragging the results in the geoprocessing results pane into an ArcMap python window...

Any help would be much appreciated!

Thanks

Najděte více příspěvků označených jako

Accepted answers

JamesCrandall

I'm typically using the numpy and pandas libraries and they have good export/saveas functionality. Just convert your Feature Class to a numpy array and then save it to csv.

<SPAN class="keyword token">import</SPAN> numpy
<SPAN class="keyword token">import</SPAN> pandas <SPAN class="keyword token">as</SPAN> pd
fc <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'H:\Documents\ArcGIS\Default.gdb\MyFeatureClass'</SPAN>
nparr <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>FeatureClassToNumPyArray<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="string token">'Field1'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Field2'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Field3'</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
pdarr <SPAN class="operator token">=</SPAN> pd<SPAN class="punctuation token">.</SPAN>DataFrame<SPAN class="punctuation token">(</SPAN>nparr<SPAN class="punctuation token">)</SPAN>

pdarr<SPAN class="punctuation token">.</SPAN>to_csv<SPAN class="punctuation token">(</SPAN>r<SPAN class="string token">'H:\nparr.csv'</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>

All comments

Anonymous User

James Crandall‌ Excellent advice to use pandas and numpy. If saved me a lot of headache trying to make arcpy.TableToTable_conversion work in a certain out field order.

By the way for others that don't know you can exclude the index added to the front of each row using index=False

<SPAN class="" style="color: #999999; border: 0px; font-weight: inherit;">pdarr.to_csv(fileName,header=False, index=False)
</SPAN>

JamesCrandall

I'm typically using the numpy and pandas libraries and they have good export/saveas functionality. Just convert your Feature Class to a numpy array and then save it to csv.

<SPAN class="keyword token">import</SPAN> numpy
<SPAN class="keyword token">import</SPAN> pandas <SPAN class="keyword token">as</SPAN> pd
fc <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'H:\Documents\ArcGIS\Default.gdb\MyFeatureClass'</SPAN>
nparr <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>FeatureClassToNumPyArray<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="string token">'Field1'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Field2'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'Field3'</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
pdarr <SPAN class="operator token">=</SPAN> pd<SPAN class="punctuation token">.</SPAN>DataFrame<SPAN class="punctuation token">(</SPAN>nparr<SPAN class="punctuation token">)</SPAN>

pdarr<SPAN class="punctuation token">.</SPAN>to_csv<SPAN class="punctuation token">(</SPAN>r<SPAN class="string token">'H:\nparr.csv'</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>

TheoFaull

I'm starting to believe I won't be able to do any of this without opening ArcMap, as the select by location tool can only deal with layer inputs from the Table of Contents, not shapefile inputs from a C drive location.

TheoFaull

No I don't want to join the datasets as I don't need any info from the polygon layer.

I just want a script that creates a CSV or Excel file of all points that fall within all the polygons where Type Code = 0 (and then 1, 2, 3)

NeilAyres

You say :

The points fall inside the various polygons throughout the map

And

Select all points that intersect the polygons that have a type code of 0, and export this selection of points to a new CSV

So you are not doing a join here? No info from the polygons is joined to the points?

So why not just create 4 layers in your map using definition queries, then use the export to excel tool?

If you know nothing about scripting or python, this will be the simpler way.

Připojte se k diskusi

Přihlášení členové mohou přispívat, sledovat aktualizace a další. Jste tu noví? Zaregistrujte si bezplatný účet.

Rychlé odkazy

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