Latest Contributions by BobNoyes

So, I have created a script that adds 6 fields to a shapefile: Situs_city, Situs_St, Situs_Zip, Latitude, Longitude, and GIS_Acres. Further, the script is set to calculate the values for the Latitude, Longitude, and GIS_Acres fields. The script will work fine until after creating the GIS_acres field. Once this field is created the next step is to calculate and that is when I get an Error 999999: cannot acquire a lock. The weird aspect is once in a while the script will run all the way through just fine... anyway, I am stumped. I am pretty new to creating python scripts. If anyone has a suggestion that would be awesome. import arcpy import os # Step 1  -  Copy new [Taxlots_TEST_DeleteME.shp] from \\cove\Department Shares\Common\Assessor\ArcGIS to E:\STAGING (IKRIT) #set variables in_data = "S:\Common\Assessor\ArcGIS\Taxlots_TEST_DeleteME.shp" Shapefile = "E:\\STAGING\\Taxlots_TEST_DeleteME.shp" #out_data = "E:\STAGING\Taxlots_TEST_DeleteME.shp" arcpy.Copy_management(in_data, Shapefile) print "Step 1 complete" #Step 2: Add Situs fields # Create a new field - Situs_City (string, 25) arcpy.AddField_management(Shapefile, "Situs_City", "TEXT", "","","25","","NULLABLE","NON_REQUIRED","") print "Step 2a: Add Situs_City complete." # Create a new field - Situs_St (string, 2) arcpy.AddField_management(Shapefile, "Situs_St", "TEXT", "#", "#", "2", "#", "NULLABLE", "NON_REQUIRED", "#") print "Step 2b: Add Situs_St complete." # Create a new field - Situs_Zip (string, 10) arcpy.AddField_management(Shapefile, "Situs_Zip", "TEXT", "#", "#", "10", "#", "NULLABLE", "NON_REQUIRED", "#") print "Step 2c: Add Situs_Zip complete." print "Step 2: Create situs fields complete." # Step 3: Create and calculate Latitude and Longitude fields. latLonRef = "Coordinate Systems\Geographic Coordinate Systems\World\WGS 1984.prj" Taxlot_shp = "E:\STAGING\Taxlots_TEST_DeleteME.shp" featureClassesList = Taxlot_shp.split(";") field_Type = "DOUBLE" field_precision_1 = 12 field_scale_1 = 8 for featureClass in featureClassesList:  arcpy.AddMessage("Calculating XY coordinates for: " + featureClass)  arcpy.AddField_management(featureClass, "Latitude", field_Type, field_precision_1, field_scale_1)  arcpy.AddField_management(featureClass, "Longitude", field_Type, field_precision_1, field_scale_1)  rows = arcpy.UpdateCursor(featureClass, "", latLonRef) for row in rows:  feat = row.getValue("shape")  cent = feat.centroid                 # To get the polygon area: cent = feat.area  row.Latitude = cent.Y  row.Longitude = cent.X  rows.updateRow(row)   #arcpy.AddMessage(str(lat) + ", " + str(lon)) print "Step 3: Add Lat and Long complete" # Step 4: Create a new field - GIS_Acres (Double, 15, 3) Shapefile3 = "E:\STAGING\Taxlots_TEST_DeleteME.shp" #Set local variables field_Name = "GIS_Acres" field_Type = "DOUBLE" field_Precision = 15 #total number of digits stored field_Scale = 4 #number of decimals places arcpy.AddField_management(Shapefile3, field_Name, field_Type, field_Precision, field_Scale) arcpy.CalculateField_management(Shapefile3, field_Name, '!SHAPE.area@ACRES!', "PYTHON_9.3") print "Step 4: Calculate acres complete" print "Congratulations! You have completed adding fields to to the Taxlots feature class. " ... View more

I am building a wall map for our SAR team and the regional SAR group which includes 9 counties in southern Oregon and Northern California. They would like to have a UTM grid on it and while it is primarily in UTM T 10 the eastern portion is in UTM T 11. This creates issues with grid angles and at this point I haven’t found a projection that works with 2 different utm zones on the same map. Any ideas suggestions or ideas? What am I missing? ? ... View more

Here is the problem I am hoping to get assistance with. I have a large area of points (covers 9 counties) for which I have created a grid polygon that has 382 boxes . I would like to write a python script or create a model that will select each box of the grid and then using that selected box do a "selection by location" of the point feature class finding the point with the high value in a column in the attribute table. The highest value would get a 1 in a field called "grid_ht".  The code or model would work sequentially through each box of the grid and finding the highest valued point in that location. What I am trying to do is find and mark the highest point (elevation) within each part of the grid. It would need to be able ignore boxes that do not have any point features with their boundaries. Anyone have suggestions? So: -  Select Grid ID #1 - Select by location - Summit points - find highest point in that box and put a 1 in it's Grid_Ht column - select Grid ID #2 - Select by location - Summit points - find highest point in that box and put a 1 in it's Grid_Ht column - Select Grid ID #3 and so on. ... View more

I would like to add a feature in a script that I run monthly that will find a replace some text in the script. In other words when I run the script the first thing that happens is a popup opens which prompts for the text to find and the text to replace. When finished and the popup closed the scripts runs through completion. The text that needs to be found and replaced is a .dbf file that has a name change each month. So, 16AprRes.DBF becomes 16MayRes.DBF. In the attached script the goal is to add the find and replace popup that will cover steps 7 and 10. Anyone have suggestions? Thanks! ... View more