Latest Contributions by curtvprice

Hi! Im having some 10-15 rasters consisting of different soils which I want to be only one raster file. Im having some problems finding out how to proceed as the mosaic operators doesn´t seem to work. The values are just 1 where there this soil type occur, eg. gleysol - 1. The rest is no data. The Cell Statistics Tool calculates stats cell by cell stats across a list of rasters. The MAXIMUM option for that tool sounds like what you're looking for. Be sure and check the box to ignore NODATA cells and the cells that have a NODATA in the stack will still come through to the output, as long as there is a data cell in the stack of inputs. ... View more

PLEASE HELP URGENT - Is there anyone at all who can help with this?  I wrote a script based on the Extract Value To Points tool but it doesn't work properly. You have incorrect syntax for ExtractValuesToPoints. You need to add the toolbox name after it if you prefix the tool name with arcpy. I'm pretty sure your script is failing (the tool is not running). You are getting a python error which is not getting printed by your except block. There are good examples in the online help on how to use python exceptions to capture both kinds of errors and report more helpful messages. # script to extract raster values to point shapefile # adapted by Scott MacDonald Kingston University k0848626 GSM655Programming Asssignment C # import system modules import arcpy from arcpy import env from arcpy.sa import * # get user supplied path, layers and fields path = arcpy.GetParameterAsText(0) arcpy.env.workspace = path # in Feature Layer is Confluence_Pts in_point_features = arcpy.GetParameterAsText(1) # in Raster Layer is DEM in_raster = arcpy.GetParameterAsText(2) # out Feature Class out_point_features = arcpy.GetParameterAsText(3) # error trapping measures try: # check out spatial analyst extension arcpy.CheckOutExtension("Spatial") # run the ExractValuesToPoints tool arcpy.ExtractValuesToPoints_sa(in_point_features, in_raster, out_point_features)         # because of your "from arcpy.sa import *" above you could perhaps do this too: # ExtractValuesToPoints(in_point_features, in_raster, out_point_features) # check in spatial analyst extension arcpy.CheckInExtension("Spatial") except Exception, msg: print arcpy.GetMessages() # ArcGIS tool messages         print "Python errors:\n",str(msg) # python errors ... View more

Batch is an older functionality that does not work with ModelBuilder. I recommend instead creating a new model with an iterator and embedding your model inside it. This has the added benefit that you can use ModelBuilder to name the outputs using model variables, instead of having to create and validate each output dataset name. ... View more

My advice is to do this process interactively and look at the results window. "Copy As Python Snippet" is your friend! ... View more

I think the root of your problems may be this line: dbo_GeocodeAddresses3 = "\\\\MyDrive\\ArcGIS\\Default.gdb\\dbo_GeocodeAddresses3.shp" You can't have a shapefile inside a geodatabase; this may be messing up your parameter validation and adding .shp to the output feature class name. Just a tip, try "raw" strings to make it easier to copy and paste paths into code. arcpy.env.workspace = r"\\MyDrive\ArcGIS\Default.gdb" ... View more

The geoprocessing team posted a blog post just on this topic! http://blogs.esri.com/esri/arcgis/2012/12/14/how-to-debug-python-toolboxes-in-3-easy-steps/ ... View more

I have a script which I've added to a toolbox in ArcGIS 10.0.  The script asks the user for a geodatabase (I'm creating and storing rasters from a NetCDF file) and either uses an existing geodatabase or creates a new one if it doesn't exist.  In adding the script to a toolbox the problem I'm having is what parameter data type to select and/or what parameter properties need to be set to allow me to put in a geodatabase that doesn't exist. The script runs if an existing geodatabase is selected but not if a new one is selected. The only way you can do this and have the validation not delete the output geodatbase before the script begins is to provide two parameters: 1) an existing folder (Input,Folder,Required) where the geodatabase exists and 2) a geodatabase name (Input,Folder,String) (to which you can add a ".gdb" to within your script [or validation code] if it isn't provided).  If the geodatabase does not exist, your script must create it using the CreateFileGDB tool. If you want your script to return the gdb as an output (say, for modelbuilder to chain into the next tool) you put it last as a parameter of type (Derived,Workspace) and return the path in your script using arcpy.SetParameterAsText() in your script. ... View more

I never got a direct answer to my question in the post you replied to, but maybe you're having issues I posted in the link I provided Linda, have you updated to SP4 or SP5? NIM069368 (reportedly fixed in 10.0 SP 4) When running the Spatial Analyst > Hydrology > Fill tool and specifying a Z limit, ArcGIS Version 10 SP2 fills less sinks than ArcGIS Version 9.3.1 SP2. ... View more

This has certainly got to be a bug that needs to be addressed.  The best way to make sure a bug fix happens is to submit a support incident. If multiple people report the problem, it's even more likely to be prioritized by the development teams. Although it's possible developers may read this forum, a support incident captures all the information the developers need for the fix. So please do report bugs like this, and if the bug is logged feel free to post the NIM# here in the forums to further grease the wheels. (this saves time for other users reporting the same bug). ... View more

I certainly agree this is a bug and encourage you to report to support so they can get it on the NIM train. I have used this functionality to create a event table from geographic coordinates stored as fields in a UTM shapefile; I don't think it's far-fetched at all to expect this to work! ... View more

This little function should do the trick for you. I did not get fancy with the xslt, just used the canned "exact copy of" stylesheet that's shipped with ArcGIS to dump the entire metadata record out of the feature class for you. If the metadata is FGDC and hasn't been upgraded to 10x, the description is in the FGDC tag, so I handled that. This function could be modified to get the title or other metadata content you're interested in, as long as you know its "tag path" in the XML. def GetDescription(item):     """Get Description text from item metadata"""     import os     import arcpy     from xml.etree.ElementTree import ElementTree as ET     # get description     AGSHOME = arcpy.GetInstallInfo("Desktop")["InstallDir"]     exactCopy = os.path.join(AGSHOME,"Metadata/Stylesheets/gpTools/exact copy of.xslt")     tmpXML = arcpy.CreateScratchName("",".xml","file",os.environ["TEMP"])     arcpy.XSLTransform_conversion(item, exactCopy, tmpXML)     tree = ET().parse(tmpXML)     os.remove(tmpXML)     node = tree.find("dataIdInfo/idPurp")     # if not found try to extract FGDC CSDGM purpose element     if node == None:         node = tree.find("idinfo/descript/purpose")        try:         return node.text     except:         return None Here's a more generalized version: def GetMetadataText(item,elementPath):     """Get metadata node text attribute from item metadata     example     GetMetadataText(r"C:\Users\cprice\Documents\boxdd.shp","dataIdInfo/idPurp")     """     import os     import arcpy     from xml.etree.ElementTree import ElementTree as ET     AGSHOME = arcpy.GetInstallInfo("Desktop")["InstallDir"]     exactCopy = os.path.join(AGSHOME,"Metadata/Stylesheets/gpTools/exact copy of.xslt")     tmpXML = arcpy.CreateScratchName("",".xml","file",os.environ["TEMP"])     arcpy.XSLTransform_conversion(item, exactCopy, tmpXML)     tree = ET().parse(tmpXML)     os.remove(tmpXML)     node = tree.find(elementPath)      if node == None:         raise Exception, "Element %s not found in metadata" % elementPath            try:         return node.text     except:         raise Exception, "No text content for element %s" % elementPath ... View more

Is it possible to share your export metadata python script?  I'm heading down this path and need something to export HTML.  Thanks, Mike This is an old post, but the Export Metadata model can be run as a tool from python. ... View more

  There must be something simple I'm missing to access the featureclass .description programmatically.   The feature class has no description property. This is a property of a layer (read and write) - but not all types of layers support it. '>arcpy.mapping: Layer I don't see a description property exposed on the arcpy layer object, although it is clearly saved in the lyr file. Are you saying when you programmatically add a .lyr file to the mxd with arcpy.mapping.AddLayer, the description is reset? If you are adding a feature class, there is no description on a FC object so the path must be being added when the layer is created in the arcmap document. ... View more

Try copying the .dbf (CopyRows) to a new table and dropping the shape field (DeleteField) before you feed it to MakeXYEventLayer. CopyRows will drop the shape field, as the output is a non-spatial attribute table. ... View more