Latest Contributions by ArielleSimmons

I've been looking for a way to create a shapefile from a python list (in the case below, uniqueList). Basically, I have been using SearchCursor's to iterate through a shapefile, and I would like to export the result as a shapefile...I tried CopyFeatures...but it didn't even return a blank file (it returned nothing). import arcpy def find_overlaps(input_features, output_features):    uniqueList = []    for row in arcpy.da.SearchCursor(input_features, ('OBJECTID', 'SHAPE@', 'name')):       foundMatch = False       for row2 in uniqueList:           if row[1].equals(row2[1]):               foundMatch = True               break       if foundMatch == False           uniqueList.append(row)    return uniqueList    ## no work ##    arcpy.management.CopyFeatures(uniqueList, output_features) ... View more

Hi, I am trying to do the steps below (using Editor and ArcScan...w/ArcInfo 10 license). I am stuck on Step #2. I can't seem to find a way to get the ArcScan toolbar to recognize the Raster layer (which is in the same fgdb as the empty line feature which I created). The ArcScan toolbar just remains grayed out...I don't know what I'm missing... Any idea what might be up?? Here are the steps I've followed (from here: forums.arcgis.com/threads/11744-Polygon-to-centerline) Using ArcScan to create centerlines from polygons (roads, paths, water bodies): Step 1: Create a raster from your polygons. - Use �??features to raster�?� or �??polygon to raster�?� tool. Make sure that cell size for new raster is small enough that even narrowest sections of a polygon get represented. For sidewalks, cell size of 3ft worked, default cell size of ~10ft did not (missing parts of path, disconnected) - Imperative to Reclassify resulting raster into a �??bi-level image�?�, with all useful pixels (the ones you wish to convert to lines, typically non-0 values in your raster) reclassified as 1, and NoData reclassified as 0. After that, �??You can use either the ArcMap Unique Values or Classified rendering options to separate the raster into two colors�?�; raster must be symbolized in two unique colors. Step 2: Set up editing environment (ArcScan works inside Editing environment only) - Create a new empty line file in which newly created centerlines will be saved - Enable ArcScan extension under Customize-Extensions pulldown menu, if you haven�??t yet. - Start an Edit Session, designating the new empty line file/gdb as target; remember to create a feature template for new empty line fine - Add ArcScan toolbar (Customize �?? Toolbars) Step 3: Use ArcScan to create lines. - If everything is set up correctly, the rasterized polygons will show up as an option in ArcScan�??s pulldown menu for �??Raster:�?� target. - Click Vectorization pulldown menu, chose �??Generate Features�?� - Check to see that all of your original polygons generated a center line. Some of my shortest paths do not generate a centerline; turning Noise Level to 0 in Vectorization pulldown-Vectorization Settings helps some of these centerlines to come through; using �??Show Preview�?� in Vectorization pulldown helps you see whether changes you�??ve made to settings and options are having an affect (hit Apply for every change to show in preview) - For ESRI Help on ArcScan: http://help.arcgis.com/en/arcgisdesk...0001000000.htm ... View more

Hi all, I have been working on a data inventory script. Idea is to use various os and arcpy functions to build a simple csv (i. .txt file) catalog of my data. In my beta testing stages, this worked perfectly. However, when running it on my server, I am getting an ???999999: Error executing function??? when I am not even a 1/4th of the way through recording my server files.  Also, I noticed the .txt file I am outputting to is locked (which means my f.close() isn???t working right???but I don???t see how that could be???.looks right to me, and like I said, it works fine on my subset of data). Anyone know what could be going on here? Attached is the code.  Copy and paste as is to keep the indents. ############################################### ## Arielle Simmons                           ## ## Planner/GIS Specialist                    ## ## Pioneer Valley Planning Commission        ## ## Date created: September 22, 2012          ## ## Date modified:                            ## ############################################### ## This script records the complete pathways of ## specified GIS data for documentation and inventory ## purposes. Also prints the date of the shapefile ## (Warning: this date could be EITHER the date ## modified or the creation date of the file...it ## is only pointing to the date of the file with ## that has been identified by the search pattern). ## Results of this script are all printed to ## an excel readable file via csv. ## This script can be modified to search ## for any ESRI feature class type (i.e. shapefile, geodatabase, etc). ## Due to the volume of data present, it is recommended that this ## script only be run for one feature class at a time. To do this, ## modify the following: ## 1) the search pattern which is based ##     on file extensions (i.e. pattern = r"(file extension)") ## 2) the 'shp' variable name (rename for better identification of the ##     data type) ## 3) the directory where the feature classes you wish to search for ##     are contained (i.e. folder = r"(folder)" ## 4) BEFORE RUNNING THIS CODE, ##     In explorer, create and name a .txt file according to the ##     file extension you are configuring the code to search ##     on. MAKE SURE to set the f = open("filename.txt") path to it!!! ## Import modules and create geoprocessor object ## the csv module contains classes for working with csv files import arcpy import os import fnmatch import sys import string import os.path, time import datetime import csv import traceback from arcpy import env # Set the workspace. List all of the shapefiles # input workspace arcpy.env.workspace = r"P:\GIS\Python_test\Ari_Python_script" try: # returns a list of shapefiles in the specified folder # that match the specified wildcard search pattern        folder =r"P:\GIS\Raid"        pattern = "*.shp"        shpList = []        for path, dirs, files in os.walk(folder):               for filename in fnmatch.filter(files, pattern):                      shpList.append(os.path.join(path, filename)) # Warning: arcpy.AddMessage is only for importing this script into the # ArcGIS ModelBuilder. Uncomment out to have it run properly. # arcpy.AddMessage("Shapefiles returned:") # open an output csv file        f = open("P:\GIS\Python_test\Data_inventory\shapefiles_test.txt", 'w') # for loop for printing outputs        for shp in shpList: # Warning: arcpy.AddMessage is only for importing this script into the # ArcGIS ModelBuilder. Uncomment out to have it run properly. # arcpy.AddMessage(shp) # first create a datetime object (of when the file was modified OR created # could be EITHER) which captures the dates of the files that the # wildcard search is looking for...convert it to a string               date_object = datetime.datetime.fromtimestamp(os.path.getmtime(shp)) #               date = date_object.strftime('%m/%d/%Y') # this is the 'for' loop that generates... # ...a list of Fields, then it # prints the list of fields for each featureclass               allFields = ""               fieldList = arcpy.ListFields(shp)               for field in fieldList:                      allFields = allFields + field.name + ", " # concatenate the path, filename, date of when file was modified OR created # (COULD BE EITHER), and field names into the identified csv file                      f.write(shp + ", " + os.path.basename(shp) + ", " + date + ", " + allFields + "\n") # print path               print shp               # print file name               print os.path.basename(shp) + "," # print date of when file was modified and or created               print datetime.datetime.fromtimestamp(os.path.getmtime(shp))        f.close() except arcpy.ExecuteError: # get the geoprocessing error messages        print arcpy.GetMessages(2) except Exception as e: # print the error message        print e.message ... View more

Hi, I am in a roadblock with ModelBuilder (using ArcView 10, SP2). I have a rather simple bit of model where I am just adding a field to a road dataset, populating it, and then converting it (based on the populated field) to Raster.  Simple beans, right? I've gotten this process to work independently outside of ModelBuilder using the same tools I've put into the Model Builder GUI. However, it will NOT RUN in ModelBuilder  :/!!! Enclosing the piece of the Model. Notably it keeps shutting down at "Polyline to Raster" (I have also tried w/ Feature to Raster, same results). As far as datasets go, this seems like a ridiculously small dataset (only 9881 records). The error at crash (yes, ArcMap fully crashes routinely right at 'Polyline To Raster') is either a '99999' error or just a regular 'ArcMap has committed a serious error'. Any thoughts? ... View more

Hi, I have this following script that I am trying to test out using IDLE 2.6.5. Notably, I've had a very similar script run before using Arc 9.3... What I'd like to do: 1) Write a script that looks into a directory, creates a list of shapefiles, and from this list batch projects all the shapefiles in the directory into a specific gdb (i.e. always to be named 'input.gdb'...there is one already existing in each directory). Output feature classes will have a specific name. 2) I would like to be able to use the script in model builder and make the directory user defined (i.e. variable folder in the below script) will be sys.argv(1). A little uncertain how to do this..but theoretically makes sense, right? 3) I am fine hard-coding the search pattern with a wildcard..."*.shp"...mostly because I doubt some of my users will know what a wildcard is... SO far my script looks as thus: ##Convert and reporoject shapefiles from RR_Sierra_Nevada\Forests ## Import modules and create geoprocessor object import arcpy import os import fnmatch import sys import string # Set the workspace. List all of the shapefiles # input workspace arcpy.env.workspace = r"D:\RR_Sierra_Nevada" # returns a list of shapefiles in the specified folder # that match the specified wildcard search pattern folder =r"D:\RR_Sierra_Nevada\Forests\Inyo" pattern = "*.shp" shpList = [] for path, dirs, files in os.walk(folder):        for filename in fnmatch.filter(files, pattern):                shpList.append(os.path.join(path, filename)) arcpy.AddMessage("Shapefiles returned:") for shp in shpList:        arcpy.AddMessage(shp) print shp The error I keep getting is EITHER a blank return...or there is no such thing as 'import arcpy' (?). Any thoughts... ... View more

Hi all, I'm a py-noob....but not a total programmer-noob. Definitely very new to harnessing all the gp functions. I'm hoping someone can give me some good guidance. I have an insanely big folder full of (many) folders which contain (many) shapefiles and images. The projections of all the files are the same. My task is to: 1) Bulk Load the shapefiles into feature classes (according to their different geometry type- which are polygons, points, and lines). Can't do this using the Load function in Catalog (**unless I find a way to print the list of shapefiles I generated in Python to an xml document which I can load there**). 2) Compare the attributes from these different shapefiles. Create a list of where the attributes differ so that I can see to it to populate the geodb with only the COMMON and desired attributes.  In trying to DIY, I found this blog: http://www.aubreyrhea.net/gis/index.php/tag/geodatabase/ This blog (nearly) tells me how to do exactly what I need to do to get this job done...but doesn't add those critical details a noob needs to move past confounded to success (not blogger's fault: I just don't think noobs are the target audience). I did email blogger...but not that hopeful I'll hear back from her. Where I am: I have successfully set up the code to print a list of all the shapefiles in the folder. Now I'm trying to connect the dots. I would like to try her attribute tables comparison method....but it simply won't run.  I think I'm having problems with the "sys.argv[]" function (for example: in order to generate my shpList, I had to change the pattern = "*.shp"...not "sys.argv[1]" like the blog indicated). I'm also loss at what functions are used to check shapefile geometry (if it's point, polyline, or polygon)....and how to use the ArcToolbox's "To Geodatabase" function to move the items from my shpList into geodatabase feature classes. Enclosed is my code: # import modules and create the geoprocessor object import arcgisscripting, os, fnmatch, sys, string #Set the workspace. List all of the shapefiles # load toolbox # gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Conversion Tools.tbx") # returns a list of shapefiles in the specified folder # that match the specified wildcard search pattern folder = r"R:\Temp\arielle\CalEMA" pattern = "*.shp" shpList = [] for path, dirs, files in os.walk(folder): for filename in fnmatch.filter(files, pattern):   shpList.append(os.path.join(path, filename)) gp.AddMessage("Shapefiles returned:") for shp in shpList: gp.AddMessage(shp) for item in shpList:         print item # create a list of attribute field names targetLayer = sys.argv[2] targetList = gp.ListFields(targetLayer) targetNames = [] for field in targetList: targetNames.append(field.Name) # create a list of attribute field names for a list of shapefiles x = 0 fieldNames = [] for shps in shpList: fieldList = gp.ListFields(shpList ) for field in fieldList:   fieldNames.append(field.Name) x = x + 1 ... View more