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

Creating Many Shapefiles/Features using Python/Model Builder from one Shapefile.

MarkFoster1

We have a .shp with thousands of tracts each with specific number (i.e. tract 1, 2, 3....etc). Each tract has multiple contract numbers. (i.e. tract 1 = contract A, B, C...etc.) The contracts are in one field based on the tract number. Many of the contracts share tracts, but we want to see what each contract holds.

IS there a way to take the contracts and create them as individual shapefiles from the Tract shapefile. For example, we want to take contract A and all of the tracts associated with it to a shapefile.

We are new to python and any help would be helpful.

We are using Arc 10 and have version 2.6 of python

Thanks

タグ付けされた投稿をもっと見る

Accepted answers

All comments

SenthilnathanThyagarajan

You have been of great help Wayne. And sure I shall mark this as answered 🙂

T__WayneWhitley

Not much change to the script, just have to set it to accept params from the toolbox script tool you create.

This pretty much explains it better than I can because there are some good pics, etc.:
http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//00150000000n000000.htm

Anyway, the basic steps are:
1- From Catalog (or Toolbox), you can right-click and choose to create a custom toolbox.
2- Right-click the toolbox you created and select 'add script'...navigate to your script and select it, name how it is referenced, etc...
3- Create the params you need, following the help ref above....you need the 'multivalue' type to enter all your yr vals by.
4- oh, and of course, modify the script to take the input you just defined...

You can do it this way, for example, just for your multivalue param (if it is the 1st one you defined):

yrs = arcpy.GetParameterAsText(0)
yrs = yrs.split(';')

Hope that helps - and please don't forget to mark this post as answered!

SenthilnathanThyagarajan

Yes, I would like to make it more interactive. So I guess in the parameters I use feature class and for output I define the workspace. But how much changes would have to be done to the code and how to set the sql query.

Thank You.

T__WayneWhitley

So it works now?

I think what you are asking is how to make this a script tool. A script tool accepts parameters you generally run from the Toolbox (or Catalog) and has the same gui interface appearance as other ArcGIS system tools. Is this what you want, so that you can manually enter the years into a text box, several of them as what is called a multi-value input? The script will operate the same way, just 'grabbing' the list input from what you input at your interface...instead of 'hard-coding' it in the script (as it is now).

SenthilnathanThyagarajan

Wow, you have been of such great help.It works. I know too many questions I have asked but I need to know is their a way if I can make the years as a parameter. So that this tool can become more generic in sense of input,output and the field(year) can be any field I want.

T__WayneWhitley

OK, test this:

import arcpy
fc = "D:\\2010.shp"
yrs = [1991,1992,1996,1998,1999,2000,2003,2004,2005,2006,2008,2010]
subtract = ''
for yr in yrs:
     where = '"Year" <= ' + str(yr) + subtract
     subtract = ' AND "Year" > ' + str(yr)
     filename = str(yr) + '.shp'
     arcpy.FeatureClassToFeatureClass_conversion(fc, "D:\\test output", filename, where)

I just read your last message...we are almost there, but I must have screwed up on the quotes.
Standy by please, this should be a minor fix...

SenthilnathanThyagarajan


import arcpy
fc = "D:\\2010.shp"
yrs = [1991,1992,1996,1998,1999,2000,2003,2004,2005,2006,2008,2010]
for yr in yrs:
      where = '"Year" <=' + str(yr) + "'"
      filename = str(yr) + '.shp'     
      arcpy.FeatureClassToFeatureClass_conversion(fc, "D:\\test output", filename, where)

del yr

I tried the code and it gives me output for just 1991.shp

T__WayneWhitley

Yes, I just corrected that...see the last post.

There is another problem, I forgot about the 'less than or equal to' so you're actually querying for ranges correct?
So that for the 3rd item in the list, 1996, that range would actually be 1992 on, or:

1992 < the range <= 1996

- give me a minute.

SenthilnathanThyagarajan

Is this how it should be

import arcpy
fc = "D:\\2010.shp"
yrs = [1991,1992,1996,1998,1999,2000,2003,2004,2005,2006,2008,2010]
for yr in yrs:
     where = '"Year" <=' + yr
     filename = str(yr) + '.shp'
arcpy.FeatureClassToFeatureClass_conversion(fc, "D:\\test output", filename, where)

Now the error is

An invalid SQL statement was used.
An invalid SQL statement was used.
Failed to execute (FeatureClassToFeatureClass).

T__WayneWhitley

So, make sure I'm understanding, you want individual yearly outputs, correct? If so you only need one list:

yrs = [1991,1992,1996,1998,1999,2000,2003,2004,2005,2006,2008,2010]

Then you can dynamically create both the SQL qry and filenames...as in:

for yr in yrs:
     where = '"Year" <=' + str(yr) + "'"
     filename = str(yr) + '.shp'

So, modifying your existing code to add a list variable and loop:

import arcpy
fc = "D:\\2010.shp"
yrs = [1991,1992,1996,1998,1999,2000,2003,2004,2005,2006,2008,2010]

for yr in yrs:
     where = '"Year" <=' + str(yr) + "'"
     filename = str(yr) + '.shp'     
     arcpy.FeatureClassToFeatureClass_conversion(fc, "D:\\test output", filename, where)

The only problem I see with this is I thought originally you wanted some yr ranges?

SenthilnathanThyagarajan

I have to my list of years ready for which I want to create the shapefiles starting from (1991,1992,1996,1998,1999,2000,2003,2004,22005,2006,2008,2010).

So simultaneously I need a shapefile for each of these years. The first one was pretty straightforward but now looping is not that I understand.

T__WayneWhitley

Great, the rest should be relatively easy...

So do you have just a few queries to run? If so, we can save the inline sub for later - what I mean is you can put your whole queries in a list:

myList = ['strQuery1', 'strQuery2', 'strQuery3']

Then loop on the list as Mathew did... if the code bails on a loop, then you now know it's likely it failed on a where clause...

Do something like:

myList = ['strQuery1', 'strQuery2', 'strQuery3']

for eachWhere in myList:
     # execute Feature Class To Feature class
     # (substitute 'where' variable with the loop variable, 'eachWhere' [no quotes])

I did forget something critical though - you need something to create dynamically changing output shapefile names so you don't end up overwriting your output. So what you could do is either in-line var substitution or if it is a standard query like the one you already have successfully executed, you could only make a list of the values, say 1990, 1991, etc., and name the shapefiles accordingly...

EDIT: Okay, the following is just a little bit clunky, but if you want to come up with your own string filenames and you just can't parse it very readily related to vals in your SQL strings, you can do this:

outNames = ['name1', 'name2', 'name3']
whereClauses = ['where1', 'where2', 'where3']
for i in range(len(outNames)):
     arcpy.FeatureClassToFeatureClass_conversion(fc, "D:\\test output", outNames, whereClauses)

SenthilnathanThyagarajan

I ran the code above and it creates the 1990.shp perfectly 🙂

How do i get the loop for all the years so that it creates shapefile for each of the years.

T__WayneWhitley

Try the code you posted again - I think your only problem there was your value was quoted - since it is a double val, remove the quotes:
where = '"Year"<= 1990'

...run the below again (I inserted the correction - also, make sure the 'test output' folder exists):

import arcpy
fc = "D:\\2010.shp"
where = '"Year" <= 1990'
arcpy.FeatureClassToFeatureClass_conversion(fc, "D:\\test output", "1990.shp", where)

SenthilnathanThyagarajan

This is the error when I run the code :

AttributeError: 'module' object has no attribute 'workspace'

yes the field Year used in Double


import arcpy

# lines preceded with the '#' symbol are comment lines (not code)
# your workspace where all input/output will go for now
arcpy.env.workspace = r'D:\\test output'

# your input shapefile
inputShp = '"D:\\2010.shp'  # What is your input shapefile?

# your output shapefile
outputShp = 'yr1991.shp'

# where clause -- substitute 'somefield' with your field (keeping the quotes)
# also, is your yr value numeric?
where = "'Year' = 1990"

# Form Feature Class to Feature Class command based on above params
# FeatureClassToFeatureClass_conversion (in_features, out_path, out_name, {where_clause}, {field_mapping}, {config_keyword})
arcpy.FeatureClassToFeatureClass_conversion(inputShp, arcpy.workspace, outputShp, where)

SenthilnathanThyagarajan


import arcpy
## below is the shapefile which I want to split it into
fc = "D:\\2010.shp" #

#the sql string which I want all the features to be selected for
where = '"Year"<= "1990"'
arcpy.FeatureClassToFeatureClass_conversion(fc, "D:\\test output", "1990.shp", where)

T__WayneWhitley

Bear in mind, that was pseudocode Mathew provided, just a general 'untested' guide... it is good though. Let's take this a step at a time. First, for readability of your code, enclose it with tags, at the end, and

 at the beginning.

Forget about the loop and inline substitution for a moment - the 1st priority is to get a single execution of Feature Class To Feature Class going with 1 properly formed query...so do something like this:


import arcpy

# lines preceded with the '#' symbol are comment lines (not code)
# your workspace where all input/output will go for now
arcpy.env.workspace = r'D:\\'

# your input shapefile
inputShp = 'yr2010.shp'  # What is your input shapefile?

# your output shapefile
outputShp = 'yr1991.shp'

# where clause -- substitute 'somefield' with your field (keeping the quotes)
# also, is your yr value numeric?
where = "'somefield' = 1991"

# Form Feature Class to Feature Class command based on above params
# FeatureClassToFeatureClass_conversion (in_features, out_path, out_name, {where_clause}, {field_mapping}, {config_keyword})
arcpy.FeatureClassToFeatureClass_conversion(inputShp, arcpy.workspace, outputShp, where)

SenthilnathanThyagarajan

I tried something like this , but this throws and exception error.

import arcpy
fc = "D:\\2010.shp"
Year = ["1991", "1992", "1997"]
for yr in Year:
    where = "contractfield like '%"+1+"%'"
    arcpy.FeatureClassToFeatureClass_conversion(fc, "D:\\test output", "somename"+yr+ as integer".shp", where)

T__WayneWhitley

You can still modify Mathew's code. You'd need to set up a SQL query for the year range and substitute that into say a Make Feature Layer, then feed that result into the feature export (think he use Feature Class To Feature Class but I think you can also use Copy Features). You may need to test the queries 1st to make sure they work properly. If it doesn't make the list of query parameters too long, you could conceivably place the entire query strings in there, then call each in turn in the loop.

EDIT: Actually, I see why Mathew suggested the use of Feature Class To Feature Class - you can enter the SQL query (where_clause) there and make your export in a single line of code - here's the syntax:

FeatureClassToFeatureClass_conversion (in_features, out_path, out_name, {where_clause}, {field_mapping}, {config_keyword})

Feature Class To Feature Class (Conversion)
Resource Center » Professional Library » Geoprocessing » Geoprocessing tool reference » Conversion toolbox » To Geodatabase toolset
http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//001200000020000000

SenthilnathanThyagarajan

Thank you for the response.I would like to use python so that I can loop the years and then accordingly create shapefiles such as 1990.shp ,1995.shp.

T__WayneWhitley

Either solution above, slightly adapted, would work for you. If you want to avoid modifying Python code, then you can 'preprocess' what you need manually before feeding into the ready-made tool Dan was talking about above.

For example, you said you had a single shapefile with a 'year' field which you want to use year ranges with to produce the 'split layers by attributes', so create an extra field in which to code your ranges and select by attributes and code each case you want segregated...once completed, you can run the tool execution to wham bam create the 'split layers', voila! (assuming you are already familiar with 'Select by Attributes' and the 'Field Calculator')

So to be clear, this would be split into 2 by the new field, 'newYrRange':

'oldYr', 'newYrRange'
1985, 'before 1990'
1980, 'before 1990'
1776, 'before 1990'
1991, '1990 to 1995'
1994, '1990 to 1995'

SenthilnathanThyagarajan

Any solutions for this yet. I have a similar question though. I have a shapefile with a field year in it. I want to create multiple shapefiles based on the years. For instance a shapefile for all building before 1990, then a shapefile for building before 1995. Any suggestions here ?

Thank you.

DanPatterson_Retired

Split layers by attributes
http://resources.arcgis.com/gallery/file/geoprocessing/details?entryID=37AEB018-1422-2418-A036-CA6D9920F808

MathewCoyle

Some sample data/screen shots would help with next steps, without knowing details, my best guess would be to do something like this.

Get a list of contracts, export to shapefile with a where clause of each contract possibility and export to shapefile.

Code example, untested.

fc = "C:/GIS/test.shp"
contracts = ["A", "B", "C"]
for cont in contracts:
    where = "contractfield like '%"+i+"%'"
    arcpy.FeatureClassToFeatureClass_conversion(templayer, exportdir, "somename"+cont+".shp", where)

会話に参加する

サインインしたメンバーは投稿、更新のフォローなどができます。初めてですか？無料アカウントを登録してください。

クイックリンク

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