Latest Contributions by ShaunWalbridge

Chuck, OGR will work fine with the File GDB, just build against the API and you'll have a FGDB driver from within GDAL/OGR and its SWIG Python bindings. If you're working with features, you can also use Sean Gilles' Pythonic wrapper for OGR, Fiona. And as Vince mentioned, you can use SWIG to create a wrapper to the API as it's just C++, there's also an idea for this if you think it should be something that Esri supports. That said, for many use-cases GDAL/OGR works well directly. cheers, Shaun ... View more

Hi Ryan, I'm looking into it, but it may take a while before I hear back anything definitive -- today's the end of the user conference, and many folks are traveling or otherwise unavailable at the moment. cheers, Shaun ... View more

I can confirm this is a bug at 10.1SP1. I've made a simple Python toolbox to recreate the issue. It's the UC next week, so many of us will be traveling, but if you don't hear back soon, remind me to submit a bug report. Here's the full Problem signature: Problem signature:   Problem Event Name: BEX   Application Name: ArcMap.exe   Application Version: 10.1.1.3143   Application Timestamp: 5058f74e   Fault Module Name: MSVCR90.dll   Fault Module Version: 9.0.30729.6161   Fault Module Timestamp: 4dace5b9   Exception Offset: 00074361   Exception Code: c0000417   Exception Data: 00000000   OS Version: 6.1.7601.2.1.0.256.4   Locale ID: 1033   Additional Information 1: 3694   Additional Information 2: 3694321e13ff624a62165c908c6077b3   Additional Information 3: 9413   Additional Information 4: 941366b9226eef57f2089e0bc55a2c5f cheers, Shaun ... View more

Just to add a tip for those who may be experiencing the same problem: If you're testing a Python Toolbox from within ArcMap, you need to start a new instance of the tool in order to overwrite any previous Map spatial reference settings that you have since changed. One way around this is to force the relevant code to reload each time the toolbox is opened. That's what I generally do when developing a toolbox, then remove them for a release. ... View more

Yes, I agree with Dan, that's a good place to start if you're just interested in dipping your toes into programming. Eric Pimpler's book, Programming ArcGIS 10.1 with Python is a good follow-up once you have some of the basics down, and if you do decide to learn more Python, consider reading Learn Python the Hard Way, which is free online with an optional for-pay video guide. cheers, Shaun ... View more

Austin, That sounds like quite an annoying bug, I imagine it took some time to track down the hacky solution you ended up with. I'd like to reproduce the bug so we can get an issue on this, as from my reading of the help, NumPyArrayToRaster should respect arcpy.env.outputCoordinateSystem. Could you just share a couple of details: What are the versions of OS and ArcGIS that you're using to produce this issue? Do I need to use a certain projection for the map, or anything that differs from the data will trigger this? cheers, Shaun ... View more

The page size property (pageSize) is available, and often you can do something with that and scaled offsets. From what I know, currently page units aren't exposed to ArcPy, though you can get to them via comtypes or ctypes if they're critical. cheers, Shaun ... View more

I don't have an exact answer, but my understanding is, this is one of those questions which is easy to ask, but in turn generates many related questions before you can figure out what specifically you'd like. See, for example, the answer on finding the center of geometry of an object, and computing an average polygon width, in addition to the thread you'd linked to. Can you provide more details on the kind of analysis you'd like to perform, and what goals you'd like to accomplish? That'll probably help as there are likely many subjective alternatives to compute these averages. cheers, Shaun ... View more

You can do this using arcpy.mapping, which exposes many of the symbolization options via Python. Check into RasterClassifiedSymbology, which should let you do what you need. First, save a .lyr file with the palette you'd like to use, then use arcpy.mapping.updateLayer to apply the symbology to your layer, and finally perform reclassification (e.g. lyr.symbology.reclassify())or other changes to make the color map match the input data. cheers, Shaun ... View more

One way of getting this information is by running patchfinder.exe (download link for 10.1 version), which will scan the registry for any and all applications and patches currently installed. It can be run from a command-line context, and have its output written to a specific location, by specifying the output file after the /s option:   PatchFinder.exe /s patches.log You could call that utility and parse the resulting text all in Python. Not a direct solution, but may help in a pinch. ... View more

Do you have the input raster data which the contours are derived from? If you do have them, you can create new contours at the desired intervals directly, and not have to do any 'tricks' to guess where they might be. If you don't have them, you could probably make a function which measured the distances between the 25/50/75m lines and made a guess at where the next intervals would lie, but it'd be just a simple picture and not based on the geomorphology of the landscape. cheers, Shaun ... View more

You could subclass arcpy.Array to include a __getitem__ method, and then handle slicing and indexing there, if you want to be Pythonic about accessing the elements. However, Mark's solution is probably the simplest answer. cheers, Shaun ... View more

Sam, Can you provide a small sample of data that contains the issue you're seeing? cheers, Shaun ... View more