Latest Contributions by DanPatterson_Retired

Luke...a quick test... # Globals outputFeatureClass = r"c:\temp\donutdemo.shp" coordList = [[[[0,0],[0,10],[10,10],[10,0],[0,0]],               [[1,1],[2,1],[2,2],[1,2],[1,1]],               [[8,8],[9,8],[9,9],[8,9],[8,8]],               [[8,1],[9,1],[9,2],[8,2],[8,1]],               [[1,8],[2,8],[2,9],[1,9],[1,8]]],              [[[10,10], [11,11], [12,10], [10,10]]]] def main():   array = arcpy.Array()   point = arcpy.Point()   # Create a list to store the features   features = []   # Read the coordinates   for feature in coordList:     print "feature", feature     for part in feature:       for coordPair in part:         point.X = coordPair[0]         point.Y = coordPair[1]         array.add(point)       null_point = arcpy.Point()       array.add(null_point)       # Create the polygon object       polygon = arcpy.Polygon(array)     # Clear the array for the next feature     array.removeAll()     # Append to the feature list     features.append(polygon)   # Copy the features to an output feature class   arcpy.CopyFeatures_management(features, outputFeatureClass) if __name__ == '__main__':   import arcpy   arcpy.env.overwriteOutput = True   main() ... View more

Karen check all syntax on the online help http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html?TopicName=welcome or for quick reference, launch your ide and use the following >>> from arcpy import sa >>> dir(sa) it will print a list and you will see what I mean ... View more

begin by checking case which may be the issue http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//009z000000m8000000.htm isnull != IsNull ... View more

can't comment on the whole problem but this line tsFolder = "D:\Timeseries_All" is going to cause problems since Python expects paths to be written as tsFolder = "D:\\Timeseries_All"   or tsFolder = r"D:\Timeseries_All"   or tsFolder = "D:/Timeseries_All" print it out as well, since you may be missing the required \\ or / in the full path/filename ... View more

selections are only honored when a layer is the input, if you then need a summary operation on the same file, then do it on the featureclass on disk associated with the layer. ... View more

what is the aversion to using tools in ArcToolbox?  they add the extra dimension of being used in models etc and the toolboxes and scripts are easily distributed...the only difference is clicking a button versus clicking a tool in terms of the interface ... View more

Sorry Cedric, don't have a 9.3.x installation here, acrpy does open up a whole new world which didn't exist in previous versions. ... View more

Cedric I haven't had much time to explore all the intricacies of arcpy but you could begin by exploring sa >>> from arcpy import sa >>> help(sa.Extent) provides some clues, if you know the extent that you need, then you can set it ... View more

perhaps     print gp.GetMessages() should be unindented such as print gp.GetMessages() ... View more

Given Bill's interesting example, what is returned...and what should be considered a cluster given such a circumstance? ... View more

start here http://resources.arcgis.com/gallery/file/geoprocessing download every toolbox/document/presentation that is associated with python, open them up and learn  Also, the online help files are generally pretty good once you get a sense of what is going on. ... View more

careful Chris, we don't want JB to get the idea that if he sides with you he can add his name to the BDFL list under the ESRI banner ... View more

Are you trying to define the projection for a file that doesn't have a defined coordinate system or are you trying to project geographic data (aka decimal degree data) to some other coordinate system?  If it is the former, then you need to find out from the source what the projection was before using the Define Projection Tool...if it is the latter, then select from the many projected coordinate systems that can be used for Canada or Ontario depending upon your needs (Project tool in Arctoolbox) ... View more

Perhaps the -ve was left out of the Lancashire data putting it to the east of 0 Long.  spatial join should be by coordinates, so I would check file separately in its own dataframe to ensure that the unprojected coordinates line up where they should ... View more