Latest Contributions by DuncanHornby

Ann, Looks like you are right those tools won't work with ArcSDE. I've never used ArcSDE but I know the problem you described occurs in personal and file geodatabases and the compact tool resolves that. According to this old thread it suggests the issue is with the underlying database. This page lists a whole bunch of commands one resetting minimum row ID, not sure how you would impliment those having never used ArcSDE? Duncan ... View more

Ann, Run the Compact GeoDatabase tool after you have deleted all your features. Duncan ... View more

Aulton, Firstly when you post code can you place it inside the tags, use the # button in the editor bar to do this. It makes your code easier to read. I think your problem is this line: PublicFunction GetMxApplicationFromArcMap(ByVal application As IApplication) As ESRI.ArcGIS.Carto.IMap Look at what your function is returning, an IMap object you want to return an IMXApplication object. Duncan ... View more

Just an idea (without seeing the whole code) but with the cursor have you set SubFields property and not included "RIP_CYCLE" by mistake? ... View more

Michael, I'm just throwing some ideas around as I'm not so good with Python. I've never used that .join method before so I had to look at the online help. The format they give is: string.join(list,separator)   you have f.write('\n'.join(output)) maybe the following will work instead? string.join(output,'\n') Your error log says it found 4650229 rows, thats quite a lot of data! Maybe there is some sort of limit on doing a write to a text file with the method you are using? You could try stepping through the list and writing out each item one at a time? If none of these work it maybe a dreaded geoprocessor memory leak? I find I cannot run python scripts more that a few hundred times when calling Spatial Analyst tools, it just bombs out. Duncan ... View more

Cécile, Have you tried to recompile your VBA project? Go into VBA editor and click on the menu option Debug > Compile Project. If it bombs out at the same line I recon you have some sort of namespace conflict. In this case prefix it with the object library as shown below MsgBox esriArcMap.Application.Templates.Item(esriArcMap.Application.Templates.Count - 1) Duncan ... View more

Graeme, Just thinking off the top of my head here (also assuming you are using ArcGIS 10). Why don't you create a model and then use a get count tool to feed into the calculate value tool to do an if/then statement the output being a boolean.  This could act as a precondition to the copy features tool? Duncan ... View more

Kimo, It looks like I am going to have to apologies here. What I said is correct for ArcGIS 9.3 but ESRI has disappointingly removed this functionality from ArcGIS 10. I was not aware that they had killed this off until I just tried it. It's discussed here and it seems to be something to do with a limited version of VBA. So unless you can figure out how to access arcobjects from python (here is a presentation that discusses this, you may have seen it already) you'll need to get to grips with ArcObjects in VB or c#. It's a real shame that ESRI have removed this ability. From my developers point of view about the only thing ESRI have done good with the release of version 10 is the ESRIAddIn, this makes deploying buttons and toolbars developed in .Net a whole lot easier. So sorry for getting your hopes up, thats "progress" for you... Duncan ... View more

Claire, Instead of writing a Python script, why not use model builder? You could have a feature iterator that spews out a polygon one for each row (or even grouped by a case field) and this simply feeds into a clip raster tool as the clipping polygon? Duncan ... View more

I just want to compliment Richards answer. When a python script is run from within ArcGIS as a scripting tool it's possible to call arcobjects within the calculate field tool. I've used this approach sometimes when I wanted to access the fine grade geometry tools that ArcObjects has. The QueryPointAndDistance method on IPolyline is certainly a very useful bit of wizardry I use it regularly! ... View more

Hi, I tried this code (note I do not set a workspace) and it copied the featuredataset and the featureclasses within it without error. import arcpy i = r'C:\Scratch\fGDB_Scratch.gdb\test' o = r'C:\Scratch\fGDB_AIS_Cleaned.gdb\test' arcpy.Copy_management(i,o,"FeatureDataset") Duncan ... View more

Just an idea but have you tried to link to the database using ArcCatalog? This will at least prove that you can connect to it (so not some weird permission issue) and you'll be able to see the properties of the connection so you can ensure you have a correct connection string? Well thats what I would do first... ... View more

You can do this in one of 2 ways, create a geoprocessor object and call the spatial analyst CON tool. If you want to stay within the arcobject environment then you need to be using the IConditionalOp Interface and do a conditional operation (if nodata set to zero else use existing cell value). Duncan ... View more

Mustafa, Are you using ESRI software? If so why not simply use the geoprocessor and call the Project tool? If you don't want or can't use ESRI software then you could try to use something like OGR in GDAL which is open source? Duncan ... View more