Latest Contributions by Luke_Pinner

Post your code. From what little I can see in the screenshot, your import statement lacks commas between the modules you want to import. ... View more

I don't think there's any way of releasing the license without restarting the IDE's python process. That said... I use an IDE called PyScripter which has a remote debugger. This means it can be restarted (or crashed...) without restarting (or crashing...) the entire IDE.  To restart the remote python process and release a license after running an ArcPy script, right click the python interpreter section of the IDE and select Python Engine->Reinitialize Python Engine. See attached screenshot. Or hit Ctrl+F2. [ATTACH=CONFIG]23535[/ATTACH] ... View more

You can fill in the gaps with the nearest values using the Nibble (Spatial Analyst) tool. Step 1, create an input for Nibble. Nibble doesn�??t fill NoData, so the input raster needs to have (a) original values to use as the fill value, and b) a valid value (e.g. 0), in areas to be filled.  You could use the following raster calculator statement to generate an input: Con(IsNull("your_raster"),0,"your_raster") Step 2, create a Nibble mask. This identifies the areas to be filled (NoData). (You should just be able to use your original raster as the mask) Step 3, run the Nibble tool. Something like the following raster calculator statement: Nibble("nibble_input", "nibble_mask", "DATA_ONLY") Some examples from my own data (black = NoData) below. My steps were a little more complicated as I had to mask the white areas from being filled or being used to fill. ... View more

That's because you're passing just the filename and not the full path of the shapefile to the tool.  It works when the shapefile is in c:/test_define as you have set that as your workspace. To fix, change infc = name to infc = os.path.join(path,name) ... View more

For some reason its loosing all my formatting.. sorry. Read this: How to post Python code ... View more

Yes, this is how I do it. See http://forums.arcgis.com/threads/66230-Anyone-using-the-new-Python-Toolboxes-at-10.1?p=230255&viewfull=1#post230255 ... View more

import osgeo import gdal from osgeo import ogr This is unnecessary , as you're not using the GDAL/OGR python bindings. Note the correct way to import gdal is "from osgeo import gdal" import ogr2ogr import ogrinfo ogr2ogr and ogrinfo are commandline programs not python modules. This should fail with an ImportError exception. string=str("ogr2ogr -f"+ ' "KML"'+" "+ 'throwthisaway.kml'+" "+'test.dgn') Use a list instead, it's easier to read: cmd=["ogr2ogr", "-f", "KML", "throwthisaway.kml", "test.dgn"] You can use subprocess.Popen if you want to get the ogr2ogr error message: proc = subprocess.Popen(cmd, stdout=subprocess.PIPE,stderr=subprocess.PIPE) stdout,stderr=proc.communicate() exit_code=proc.wait() if exit_code: print stderr else: print stdout ... View more

Can you edit your post and enclose your script in CODE tags. See this post if you're not sure how. ... View more

For future reference, "doesn't work" is not useful and makes it difficult for people to answer your question. What was the problem exactly - i.e. did you get an error message, did your script crash, what line did it crash on, did python crash with some sort of windows error, did your computer reboot, did smoke come out of your computer vents, etc...? 😉 Luckily, I think I can see your problem. MX is a number, you can't add a number to a string, you'll get a TypeError. You can use the following instead: "CHARTSCALE = '" + str(MX) + "'" #or "CHARTSCALE = '%s'"%MX ... View more

Sorry, the masks are 1bit - either visible or not. ... View more

A couple of workarounds are: 1 - Manually set the 4th band as alpha in the layer properties->symbology 2 - Write an internal/external mask. ArcGIS understands these, though I have no idea if Flex can write these (I use gdal). ... View more

Yes. I'm not sure why you are concerned about the 3x3 window though. I interpreted your original question as wanting to calculate slope in specific direction, not downhill slope - the method I described does exactly that, it simply uses downhill slope and aspect to calculate slope for a specified direction.  Look at the hypothetical surface in my 2nd attachment, slope (D) i.e in an easterly direction, directly across slope, is 0 as expected. Unless I've misinterpreted your question and you want to calculate a single downhill slope value using only 4 adjacent (i.e. not diagonal) cells...? - aka "Rooks case": [ATTACH=CONFIG]22272[/ATTACH] Image source: http://www.lpc.uottawa.ca/publications/moransi/moran.htm#Defining%20Contiguity ... View more

Although there are a number of slope algorithms, only one is implemented in ArcGIS (How slope works). It uses the 8 neighbbours (i.e. 3x3 around the centre cell). You can't specify a different neighbourhood without writing your own slope algorithm using focal functions. ... View more

What data type are the bands?  If you have byte or unsigned integer data, you might be getting integer overflows as you're applying the float conversion after the operations. As an example - say you have unsigned 16bit integers, band 4 = 500 and band 43 = 750. Float(500 - 750) will return 65286 but Float(500) - 750 will return -250.0. Try: 2.5 * ((Float("%Band 4%") - "%Band 3%")) / (Float("%Band 4%") + 6 * "%Band 3%" - 7.5 *"%Band 1%" + 1)) ... View more