Latest Contributions by RoyHewitt

‎07-15-2014

QGIS has allowed you to use shapefiles stored in .zip files in their Desktop GIS for quite some time. Seems like a feature that should have been added by now.

‎07-15-2014

QGIS has allowed you to use shapefiles stored in .zip files in their Desktop GIS for quite some time. Seems like a feature that should have been added by now.

‎10-18-2013

Relationship classes are great because they allow relationships between two datasets to persist between ArcMap sessions. When dealing with relationships in scripts generally I create two separate loops in order to retrieve related records. The first loop iterates a spatial file's primary key, while a second, nested loop iterates on the tabular file's foreign key. A list of related records can be created where these two keys match. An arcpy method that notes the relationship between the two datasets and allows the user to pull out related records more easily would be helpful. Perhaps the relationship class could be used as a parameter along with the two datasets since the relationship class already contains the logic for how the two datasets are related.

‎07-02-2012

I'm trying to do something similiar to this ESRI Blog Post to detail focus areas. I need to Force RGB in order to add Image Analysis functions to achieve the effect. Since I'm doing this on several focus areas I'd like to automate it. Is there any functionality in arcpy that allows you to FORCE RGB??

‎04-20-2012

I can't for the life of me find a fix to this. I need to change the extent of my view extent to a given parcel of land (which is a user input in a script tool). My current methodology is to create a feature layer, then select it, then zoomToSelectedFeatures(). The problem is, the feature layer isn't automatically added to the Map when that line of code is executed. I've tried many fixes, and been back and forth with an ESRI Rep (who apparently grew tired of my problem and has since stopped responding to my inqueries). For more background, including photos and suggested solutions, you can check this out: http://gis.stackexchange.com/questions/23348/add-feature-layer-to-display-with-in-a-script-tool-for-desktop-mapping-app I've been trying to find an answer to this for a while, but I suspect it's the fault of the Script tool Properties.

‎02-17-2012

I've discovered the problem! The default action for both the arcpy.MakeFeatureLayer_management() and arcpy.Select_analysis() tools when run in the python window add the new layer to the map by default. When I run the same script through ArcToolbox as a custom script tool, the layer isn't added to the map by default. Calls to the feature class in the TOC fail, therefore I can't make a selection on the fclass to change the data frame extent. This seems to be a problem that people have had in the past (http://forums.esri.com/Thread.asp?c=93&f=1729&t=128409) I added the feature layer as a user input parameter, but the feature layer doesn't appear in my TOC until the script has ended (despite using the arcpy.RefreshTOC() method).

‎02-14-2012

I rewrote the code to use the for loop. I got the same problem at first (all null values in the attribute table), I went into field calculator and closed it without entering any code and VOILA, the correct values were in the table. I'm not sure if they would have been had I gone into the field calculator after running the while loop or not, but success none the less. Is there any "refresh" method that would fix this?? Thanks again!

‎02-14-2012

I'm writing a script for field staff that clips all wetlands in a watershed, calculates area (hectares) and multiples by a sampling factor to see how many hectares need to be sample Here is a sample of the update cursor part of the code: rows = arcpy.UpdateCursor(newWetland) row = rows.next() while row: if row.getValue(wetlandTypeField) == 'Estuarine and Marine Wetland': totalHectares = row.getValue(areaField) sampleArea = (totalHectares * .25) # Sample 25% of this type row.setValue(sampleField, sampleArea) rows.updateRow(row) row = rows.next() elif row.getValue(wetlandTypeField) == 'Freshwater Emergent Wetland': totalHectares = row.getValue(areaField) sampleArea = (totalHectares * .25) # Sample 25% of this type row.setValue(sampleField, sampleArea) rows.updateRow(row) row = rows.next() There are several more elif blocks for each wetland type in the dataset. When I type in the python window: print sampleArea I get the correct area, however, it does not update the row. The help files indicate that the second parameter of the row.setValue() method is supposed to be an object. Could that be the problem? Any help would be appreciated. Roy

‎02-07-2012

I'd like to write a simple script that automates the arcpad check out / check in process. I've been unable to find any helpfiles for the geoprocessing tools (I know they're only associated with the ArcPad extension). I put the tool in model builder and exported as a python script to find that the tool is: arcpy.ArcPadCheckin_ArcPad() I don't know what the parameters for the tool are. I've looked through google, the ArcMap helpfiles and ArcPad help files and haven't found anything of use. Can someone please point me in the right direction -- anything describing the parameters would be much appreciated.

‎01-27-2012

I'm writing a script that automates the creation of a map for a watershed (user input). My workflow is below; I need to maintain the correct area for wetland polygons throughout the script as I will be using the area for analysis when the map is complete. User Input Watershed Code => Clip All wetlands file with current watershed => Dissolve clipped wetlands file by wetland type to get acreage of each wetland type within shed => generate report with report layout file (.rlf) of acreage by wetland type Additionally, I'd like to know if the "repair geometry" geoprocessing tool will update the acreage, or if I should use the area property of the shape attribute in field calculator to update the area (!shape.area@acres!). The acreages will be multiplied by a ratio developed in house to detmermine the linear distance that must be walked through each type of wetland to consider the watershed as a whole "survey completed."

‎01-20-2012

Updated Script and Error: Script runs correctly up to the point of IDW. Error: 010092: Invalid output extent. # Author: Roy Hewitt # US Fish and Wildlife Service # December 2011 # Import modules, setup overwrite in environments import arcpy from arcpy.sa import * arcpy.env.overwriteOutput = True # Create file paths barrier = "P:/Employee_GIS_Data/Hewitt_GIS_Data/Python/NutriaProj.mdb/WetlandBarrier/" rhaPoints = "P:/Employee_GIS_Data/Hewitt_GIS_Data/Python/NutriaProj.mdb/RHA_Waypoints" habitatData = "P:/Employee_GIS_Data/Hewitt_GIS_Data/Python/NutriaProj.mdb/habitatDetails" outPath = "P:/Employee_GIS_Data/Hewitt_GIS_Data/Python/" rasterPath = "P:/Employee_GIS_Data/Hewitt_GIS_Data/Python/Rasters/" rhaLayer = "P:/Employee_GIS_Data/Hewitt_GIS_Data/Python/RHA_Waypoints.lyr" # Create field variables rhaField = "IDENT" habitatField = "PointName" rhaUnit = "ModelUnit" barrierUnit = "ModelUnit" # Create List Variables fieldList = [] unitList = [1,2,3,4,5,6,7,8] # Create join between RHA points and Habitat data table try: print "Joining waypoints with habitat data..." arcpy.MakeFeatureLayer_management(rhaPoints, "rhaLyr") arcpy.AddJoin_management("rhaLyr", rhaField, habitatData, habitatField) arcpy.SaveToLayerFile_management("rhaLyr", rhaLayer) except: print arcpy.GetMessages(0) # Create list of vegetation attributes for loop, local variables print "Creating list of vegetation fields for IDW modeling..." for field in arcpy.ListFields(habitatData): fieldList.append(field.name) fList = fieldList[2:17] power = 2 # As distance increases, point has less impact on interpolation cellSize = 60 # Raster cell size # Check out Spatial Analyst extension try: if arcpy.CheckExtension("spatial")== "Available": arcpy.CheckOutExtension("spatial") print "Spatial license checked out." except: print "Spatial Analyst extension not available." print arcpy.GetMessages(2) try: for unit in unitList: # Create where clause for current study area. print "Creating where clause for study unit " + str(unit) + "..." whereRHA = '['+ rhaUnit + '] = ' + "'" + str(unitList[unit - 1]) + "'" whereBarrier = '['+ barrierUnit + '] = ' + "'" + str(unitList[unit - 1]) + "'" # Create feature layer for RHA waypoints and Barrier for current study unit arcpy.MakeFeatureLayer_management("rhaLyr", "currentRHA", whereRHA) arcpy.MakeFeatureLayer_management(barrier, "currentBarrier", whereBarrier) for feat in fList: try: print "Running IDW model for " + feat + "..." # IDW geostat analyst #arcpy.Idw_ga("currentRHA", feat, "", rasterOutpath + feat + "_" + str(unit),cellSize, power, "", "currentBarrier") # IDW spatial analyst outRaster = arcpy.sa.Idw("currentRHA", "habitatDetails." + feat, cellSize, power, "", "currentBarrier") outRaster.save(outRaster + feat + "_" + str(unit)) print "Successfully ran IDW for " + feat + "." except: print feat + " interpolation failed." print arcpy.GetMessages(2) # Delete feature layer arcpy.Delete_management("currentRHA") arcpy.Delete_management("currentBarrier") except: print arcpy.GetMessages(2) finally: arcpy.CheckInExtension("spatial") print "Spatial license checked back in." print "Finished running IDW module."

‎01-19-2012

I added the check out / in methods to my script, but I don't think that's the problem (I'm loading and running the script in the python window with the extension checked in ArcMap. I suspect the problem is with the field names. Since the fields that I'm using for the z-value are from a joined table, the IDW isn't finding them. For the Z-value field I've tried: 1. Looping through Python list of fields (ex. SaltmarshBulrushPercentage) 2. Adding name of joined table (ex. habitatDetails.SaltmarshBulrushPercentage) Both of these methods return the same error: ERROR 001000: Z value field: Field SaltmarshBulrushPercentage does not exist RESOLVED Error 001000 AddJoin does not work with a feature class, only a feature layer. I rewrote the top of the script to get the join to work, now the script recognizes the field name.

‎01-19-2012

Hello All, I'm working on a script that creates an interpolated raster surface of wetland vegetation. Data points were collected for coastal marshes in boat surveys. Each point contains fields for each of the 15 different types of vegetation documented (as percent cover). The rivers are disjoint, and one river's vegetation should not affect the others (because of variables including salinity, tides, etc.). Each set of points has been separated into "ModelUnits" of which there are 8. My script has a nested loop: First, it loops through the study area, then loops through each vegetation type to create an interpolated surface (IDW). Everything seems to work until the script reaches the IDW tool: There are several IDW tools in different toolboxes (3D, Geostatistical Analysis, and Spatial Analysis). I have chosen the spatial analysis version, but I'm confused as to syntax. The help documents have only: arcpy.IDW(parameters) as opposed to: arcpy.IDW_ga(parameters). Trying IDW in the ArcMap python window, there is no option for the spatial analyst version of IDW. I've attached a map image to give you an idea of what the files look like and my python script. Any and all help is appreciated! import arcpy arcpy.env.overwriteOutput = True # Create file paths barrier = "P:/Employee_GIS_Data/Hewitt_GIS_Data/Python/NutriaProj.mdb/WetlandBarrier/" rhaPoints = "P:/Employee_GIS_Data/Hewitt_GIS_Data/Python/NutriaProj.mdb/RHA_Waypoints" habitatData = "P:/Employee_GIS_Data/Hewitt_GIS_Data/Python/NutriaProj.mdb/habitatDetails" rasterOutput = "P:/Employee_GIS_Data/Hewitt_GIS_Data/Python/NutriaProj.mdb/" # Create field variables rhaField = "IDENT" habitatField = "PointName" rhaUnit = "ModelUnit" barrierUnit = "ModelUnit" # Create List Variables fieldList = [] unitList = [1,2,3,4,5,6,7,8] # Create join between RHA points and Habitat data table try: arcpy.AddJoin_management(rhaPoints, rhaField, habitatData, habitatField) except: print arcpy.GetMessages(2) # Create list of vegetation attributes for loop, local variables for field in arcpy.ListFields(habitatData): fieldList.append(field.name) fList = fieldList[2:17] power = 2 # As distance increases, point has less impact on interpolation cellSize = 60 # Raster cell size try: for unit in unitList: # Create where clause for current study area. whereRHA = '['+ rhaUnit + '] = ' + "'" + str(unitList[unit - 1]) + "'" whereBarrier = '['+ barrierUnit + '] = ' + "'" + str(unitList[unit - 1]) + "'" # Create feature layer for RHA waypoints and Barrier for current study unit arcpy.MakeFeatureLayer_management(rhaPoints, "currentRHA", whereRHA) arcpy.MakeFeatureLayer_management(barrier, "currentBarrier", whereBarrier) for feat in fList: try: # IDW geostat analyst #arcpy.Idw_ga("currentRHA", feat, "", rasterOutpath + feat + "_" + str(unit),cellSize, power, "", "currentBarrier") # IDW spatial analyst outRaster = arcpy.Idw_sa("currentRHA", feat, cellSize, 2, "", "currentBarrier") outRaster.save(rasterOutput + feat + "_" + str(unit)) print "Successfully ran IDW for " + feat + "." except: print arcpy.GetMessages(2) # Delete feature layers arcpy.Delete_management("currentRHA") arcpy.Delete_management("currentBarrier") except: print arcpy.GetMessages(2)

‎01-17-2012

Solved the problem - Selection was made on the feature layer, which was deleted. Selection should have been made on the layer that persists throughout the script. I've read through all of the forums relating to the zoomToSelectedFeature() method and have not had any success changing the dataframe extent in my python script. I'm creating a script tool for less ArcMap inclined employees to enter in a Maryland Tax Parcel ID Code and Landowner name as parameters. The tool clips the national wetlands inventory layer and zooms to the Tax parcel, finally a PDF is exported. Everything runs as expected, but the PDF created at the end of the script doesn't represent the data frame of the selected feature. # Import modules, setup overwrite in environments import arcpy, os arcpy.env.overwriteOutput = True # User input variables mapTitle = arcpy.GetParameterAsText(0) userParcel = arcpy.GetParameterAsText(1) ownerName = arcpy.GetParameterAsText(2) # Map variables mxd = arcpy.mapping.MapDocument("P:/Employee_GIS_Data/Hewitt_GIS_Data/Parcel Automapper/ParcelAutomapper_Debug.mxd") df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0] # File Paths parcels = ---Personal FILE PATH---"/Current_Parcel_Data" captures = ---Personal FILE PATH---"/Capture_Locations" wetlands = ---Personal FILE PATH---"/Delmarva_Wetlands_MDSP" # Output Paths [removed to shorten script for forum] # Local Variables parcelIdField = "ACCTID" wetlandType = "WETLAND_TY" # Create where clause based on user input whereClause = '"' + parcelIdField + '" = ' + "'" + userParcel + "'" # Feature Class Lists fcList = ['Current_Parcel_Data', 'Capture_Locations', 'Delmarva_Wetlands_MDSP'] try: # Select parcels based on user input, zoom to selected parcel(s) parcelLayer = arcpy.mapping.ListLayers(mxd, "Current_Parcel_Data", df)[1] arcpy.SelectLayerByAttribute_management(parcelLayer,"NEW_SELECTION", whereClause) df.zoomToSelectedFeature() arcpy.RefreshActiveView() except: print arcpy.GetMessages() try:# Clip wetlands layer to extent of user dictated parcel, calculate staticstics arcpy.Clip_analysis(wetlands, "parcelLayer", shapePath + userParcel + "_Wtld") except: print arcpy.GetMessages() try: # Edit text elements (arcpy.mapping module) mxd.title = mapTitle + ' Parcel ID: ' + userParcel # update map title # Export image as .png for sharing print "Exporting map to image file." arcpy.mapping.ExportToPNG(mxd, imagePath + ownerName + "_" + userParcel, "PAGE_LAYOUT") except: print arcpy.GetMessages() finally: #mxd.save() arcpy.RefreshActiveView() print "Finished Automapper." del mxd

‎08-26-2011

Can you run IDW using this data? What version of ArcGIS are you using? Input data format - .shp, file or P geodatabase or ? Output format - grid, fGDB raster or ? Just tried to run IDW -- returned the same error message 999999 ArcInfo Desktop 10 Input data format: multiple feature classes in a feature dataset within a file geodatabase. Output format: Raster (default, seems to be the only option). Thanks, Roy

Online Status	Offline
Date Last Visited	‎11-11-2020 02:23 AM

My Ideas

Latest Contributions by RoyHewitt

Re: introduce a zipped shapefile format (SHZ)

Re: introduce a zipped shapefile format (SHZ)

An Arcpy method for dealing with related records

Re: Export Raster Data

How to programatically change view extent based on user input from script tool)

Re: Script tool yeilds different (wrong) result compared to Script run in python wind

Re: Update Cursor not updating field with value

Update Cursor not updating field with value

Accessing ArcPad Geoprocessing tools with Python

Automated Map Creation & Export - Maintain correct geometry / area

Re: IDW in python script

Re: IDW in python script

IDW in python script

dataframe: zoom to selected feature

Re: Failure or Kriging Tool

An Arcpy method for dealing with related records