Latest Contributions by DanPatterson_Retired

I failed to mention that you should add a new text field and get a text representation of the elevation.  Floating point numbers aren't supported directly.  If that doesn't work, zip the shapefile and email it to me at Dan_Patterson AT carleton DOT ca and I will have a look ... View more

try http://resources.arcgis.com/gallery/file/geoprocessing/details?entryID=37AEB018-1422-2418-A036-CA6D9920F808 ... View more

stick with Python solutions, VBA is gone, perhaps VBScript will be next, Python is stable for the next three years (aka no major changes, just bug fixes) and it seems to be the preferred scripting/programming environment within Arcmap and open source GIS packages.  Pretty anything you see in VBA can be converted to Python easily (start by removing the Dim this, Dim that stuff), the only thing that isn't great is the code block section within the field calculator which eventually will allow one to load Python scripts like the old VB scripts....translate away ... View more

Did you check Ianko's website http://www.ian-ko.com/ to see if he has provided updates for version 10?  I did and there is stuff there.  You should bookmark the site and visit it regularly to keep current with the new tools provided for ArcMap ... View more

No need to debate the merits, what has happened is what I suspected, Python would go up, Geoprocessing would go down and the confused would increase the GP/Python threads in General...so all I want is for General to be terminated 🙂 ... View more

try >>> python_test = "complain "*500 >>> print python_test at the command line of your favourite IDE 🙂 ... View more

yet geoprocessing is suffering and the General tab is increasing with both geoprocessing and python stuff...kill general or institute the "Check/Look here first" sticky ... View more

does the script in arctoolbox recognize a straight import of the shapely module? If not, then it may be an installation issue of shapely.  This capability has been added to ArcMap 10 for future reference, through the arcpy module import arcpy square = [[0,0],[0,1],[1,1],[1,0]] small_square = [[0.1, 0.1],[0.1, 0.9], [0.9,0.9], [0.9,0.1]] pnt = arcpy.Point() array = arcpy.Array() for a_pnt in square:   pnt.X = a_pnt[0]; pnt.Y = a_pnt[1]   array.add(pnt) poly1 = arcpy.Polygon(array) array.removeAll() for a_pnt in small_square:   pnt.X = a_pnt[0]; pnt.Y = a_pnt[1]   array.add(pnt) poly2 = arcpy.Polygon(array) array.removeAll() print "poly 1 contains poly2 :", poly1.contains(poly2) print "poly 2 contains poly1 :", poly2.contains(poly1) ... View more

check the online help, specifically http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//002z0000000n000000.htm case-sensitive you can import arcpy and do a help(arcpy) or dir(arcpy) to get further information or import its submodules and get further info there. ... View more

From this thread, in version 10 http://forums.arcgis.com/threads/865-quot-built-in-quot-rasters-in-python-map-algebra try the example import numpy as np import arcpy #source http://forums.arcgis.com/threads/865-quot-built-in-quot-rasters-in-python-map-algebra #with some fiddling # Setup some rasters of rows and columns # (like old $$ROWMAP and $$COLMAP) arcpy.env.workspace = "c:/temp" arcpy.env.overwriteOutput = 1 nprows = np.indices((10,10))[0] npcols = np.indices((10,10))[1] # Convert the numpy arrays to ESRI rasters (ie Raster objects) # called 'rows' and 'cols' row_ras = arcpy.NumPyArrayToRaster(nprows) row_ras.save("rowraster") col_ras = arcpy.NumPyArrayToRaster(npcols) col_ras.save("colraster") print row_ras print col_ras ... View more

you are probably trying to cover too big an area, check the help topic, it is not designed to do broad regions covering many degrees of latitude http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Area_Solar_Radiation/009z000000t5000000/ ... View more

All the files should be projected. ... View more

Luke that appears to be the case, which appears to go against the documentation which suggests that X and Y would be assigned a value of none >>> from arcpy import Point >>> pnt = Point() >>> print str(pnt) 0 0 NaN NaN >>> |  Methods inherited from arcpy.arcobjects.mixins.PointMixin: |  |  __init__(self, X=None, Y=None, Z=None, M=None, ID=None) |      Point({X}, {Y}, {Z}, {M}, {ID}) which suggests that None is being translated to 0...perhaps ESRI could provide some rationale which I presume has something to do with None not being able to be represented as a numeric value for coordinates...go figure ... View more