In version 9.3 there is a toolbox called Create Features from Text File, which can do what you require. The help files are quite detailed and are found here.
The following code sample should be simple enough to adapt to reflect your dataset:
# Create geoprocessing dispatch object
import arcgisscripting
gp = arcgisscripting.create()
# Set up inputs to tool
inTxt = r"C:\temp\StreamPoints.txt"
inSep = "." #the symbol used to indicate a decimal point
output = r"C:\temp\Streams.shp"
# Run tool
gp.CreateFeaturesFromTextFile(inTxt, inSep, output, "#")
UPDATE:
If you require numerous attributes of a range of types to be written to a feature dataset from a text file or csv file, it my be easiest to use pyshp. I wrote up an answer to a similar question about pyshp a few days ago and will reproduce the relevant parts of this answer below.
The pyshp/xbase datatyping has always been tricky for me until I found this user guide for the xbase format and as a result of this question I have written a small note on my blog regarding the relevant pyshp datatypes, part of which I have pasted below:
- C is ASCII characters
- N is a double precision integer limited to around 18 characters in length
- D is for dates in the YYYYMMDD format, with no spaces or hyphens between the sections.
- F is for floating point numbers with the same length limits as N
- L is for logical data which is stored in the shapefile's attribute table as a short integer as a 1 (true) or a 0 (false). The values it can receive are 1, 0, y, n, Y, N, T, F or the python builtins True and False
There is a DateTime type as well which should allow you to store correctly formatted time data as well using the same syntax as the above datatypes.
Here is an example python script showing the use of the various pyshp datatypes being read in from lists, which could be populated with your data using the csv module:
import shapefile as shp
#Input parameters
out_file = 'GPS_Pts.shp'
#Set up blank lists for data
x,y,id_no=[1.0,2.4,3.9],[4.3,5.5,6.1],[55,16,10]
date,notes = ['20121012','20121012','20121011'],['This','is a','string']
logic = [True,False,True]
#Set up shapefile writer and create empty fields
w = shp.Writer(shp.POINT)
w.autoBalance = 1 #ensures gemoetry and attributes match
w.field('X','F',10,5)
w.field('Y','F',10,5) #float - needed for coordinates
w.field('Date','D') #date
w.field('Notes','C',50) #string
w.field('ID_N','N',12) #double
w.field('Logic','L',1)#boolean
#loop through the data and write the shapefile
for j,k in enumerate(x):
w.point(k,y[j]) #write the geometry
w.record(k,y[j],date[j], notes[j], id_no[j], logic[j]) #write the attributes
#Save shapefile
w.save(out_file)
Your questions are answered in Rtree's docs, which come with the source and are online at http://toblerity.github.com/rtree/tutorial.html.
You can make your shapefile indexing code much simpler:
for i, shape in enumerate(shapes):
idx.insert(i, shape.bbox)
Use Python's enumerate() function for this kind of thing whenever possible. It's tidy and foolproof.
Best Answer
From the pyshp documentation page:
Just substitute your list for the
parts=argument toWriter.poly, and set whatever fields you would like to associate with your shapefile.