19

A couple of general tricks I have found useful in the past in this situation: Run your Python script as stand-alone (e.g. from IDLE, PyWin, Eclipse or preferably CMD) to remove the overhead of ArcMap. Spawning subprocesses is an old trick to solving ArcGIS memory leaks even if you don't want to parallize a process. It works because the memory is released ...


16

You can use a conditional (Con) statement in the raster calculator or the con tool in ArcToolbox. The below statement is evaluated as [IF myrasters > 2000 THEN 1800 ELSE myraster] Con("myraster" > 2000, 1800, "myraster")


15

You can use a conditional statement. The issue with previous recommendations is that when you rasterize your polygons (which is necessary) the background, that does not contain polygons, will be NoData resulting in corresponding areas in the output also being NoData. You will need to set your analysis extent to your original raster and then set a background ...


11

There is no need to write any additional ArcPy for this. Simply use dynamic legends (which were added at ArcGIS 10.1) when you are authoring your maps. The Map Extent Options allow you to control how eligible legend items interact with the current map extent. If the Only show classes that are visible in the current map extent check box is checked, a ...


10

I have tried supervised classification in ArcGIS. Firstly I would say that it is not the best software for classification. As I did it, you can create training sites as points. Just create a shapefile (or geodatabase), add Integer field, click points over your image and assign classes as numbers. (I think you can also use polygon shapefile). For ...


9

Three options I can think of: Save your map document and open it when you subsequently want to work with the shapefile. This should be obvious so perhaps you have not explained why this is not working for you. Save the layer to a layer file and load that instead of the shapefile Load the shapefile into a geodatabase and set up a cartographic representation ...


8

I realise this question is tagged 9.3, but anyone looking for the same answer at 10.1 onwards is better off using arcpy.da.Walk. It is faster and more accurate than ListDatasets/FeatureClasses/Rasters/etc. import arcpy import os for root, dirs, datasets in arcpy.da.Walk('d:\scratch.gdb'): for ds in datasets: print os.path.join(root, ds) The ...


8

Please mention the sensor of Landsat 5, is it MSS or TM? Assuming it is Thematic Mapper data, you have visible red and shortwave infrared data. You can directly infer from the band reflectance values about where vegetation patches lie and hence moisture content. Band 3 (Red) can help you discriminate vegetation slopes and Band 5 (SWIR) can help you ...


8

The add-in model is a new feature in version 10.0. This means you will not be able to use an add-in at lower versions, which require you to implement your commands, extensions, toolbars etc. in the older COM-only model.


8

GeoChalkboard can give you a good overview of the changes! Here's the list, but check it out for more detail: Replace the import statement Replace all gp. with arcpy. Check your capitalization. The result tool is different. Replace any del gp statements. The overwrite tool is different. Import the env module from arcpy. List item


8

Your question is a mélange of misinformation and cross-threaded topics, which makes it exceedingly difficult to answer. First off, there is no comparison of the "SDE API" and "MSSQL Native geodata format". ArcSDE is many things, including a code library for accessing data (not spatial data, all types of data, including geometry) in multiple SQL (and no-SQL)...


7

Consider integrating DEMs into your research on soil moisture/exposure. I have used the following indices in the past for regression models (Davies et al. 2010): Site exposure index = slope∗cos(pi∗(aspect−180)/180) (Balice et al. 2000) Heat load index = 0.039 + [0.808 * cos(l) * cos(s)] – [0.196*sin(l)*sin(s)] – [0.482*cos(a)*sin(s)] (McCune and Keon 2002)...


7

After doing some research, I found a solution and it's called Assembly Binding Redirection Basically, for any .net .exe, you can add XML code in the application config file to redirect a referenced assembly from a version number to another. Esri implemented this in 10.0 and up versions in their EsriRegasm.exe utility. (located here : %commonprogramfiles(...


7

As @Baltok alluded to, you are trying to have the Selection on a particular layer be treated in the same way as other layers are, with respect to drawing order. I think that by default, and design, ArcGIS makes the Selection override the symbology of anything that is underneath it, as a means of making it easy to find the selection. Since you are selecting ...


7

Please have a look IFeatureClass.Search Method To get feature by mouse click: public List<IFeature> Search(IPoint mousePoint, double buffer, IFeatureClass featureClass) { var envelope = mousePoint.Envelope; envelope.Expand(buffer, buffer, false); var geodataset = (IGeoDataset)featureClass; string shapeFieldName = featureClass....


7

A solution to this is discussed on the ESRI forum here: Every planar surface has a simple formula, Jesse, in the form (something)*x + (something else)*y + (constant). In the Raster Calculator, x is $$XMap and y is $$YMap. The trick is figuring out the coefficients. Often you want to specify a plane in terms of its slope (s), aspect (a), ...


6

I had a similar problem with the Tabulate Intersection tool (using 10.4). I found the suggestion to place the resulting table into a geodatabase from here: https://geonet.esri.com/thread/13680 It fixed the problem!


6

Perhaps the features were converted to graphics at some point. In that case, they might appear exactly the same as the original feature class but would not turn off with the feature class, since they are graphics. This is just a guess but worth a look.


6

I don't think I can directly link, but what you want is American Community Survey Table S1810 with New Mexico County Geography. Start here: http://factfinder2.census.gov/faces/nav/jsf/pages/searchresults.xhtml?refresh=t Enter S1810 for your search topic. This should give you four tables named Disability Characteristics. From the options on the left, now ...


6

Not necessary a screen shot, however here is how you can export the active view (refer to code snippet below). Export active view


6

First answer - No, you can not have duplicate names in a geodatabase. Yes, this is the expected functionality since the Feature-Dataset is only a record reference in a database. Any database will balk at having tables with duplicate names. It sounds like you are trying to use the Feature-Dataset model for organizational purposes versus the true purpose of it....


6

You can do this through Spatial Analyst. Adding separate bands is very simple. When you are adding layer, go to your tiff and instead of adding, double-click on it. You will enter into its bands and now you can add these bands separately. So, from this you will enter into this Then you can use Raster calculator to do the sum. Open it from the Spatial ...


6

This question asks for two things: (1) identify the boundary cells and (2) "find" their elevations. Although "find" could be interpreted in many ways (do you want an output table of elevations? A DEM containing only those cells? etc.), doing it is straightforward once you have solved (1), especially if the boundary cells have been represented as an indicator ...


6

I've written a script for 10.1/10.2, however you might see if it is possible to rewrite it for 9.3. I've almost never used arcgisscripting, but I guess it would be hard to implement the same without arcpy in 10.0+. Bearing Distance To Line GP tool I've used is available only in 10.0+, too. If 9.3 is your only option, you could try to generate the lines ...


6

Thanks for your script Alex, very nice! I saw a small issue, since relies on several geoprocessing tools that run for each row: Create Table > Add Field > Cursor to Insert Rows > Bearing Distance To Line > Append > Delete Table > Delete Lines >>> Reiterate for next row. It took 10 secs per feature for me. Since, I needed to process around 18,000 features (e....


6

Your downloaded shapefile already had a coordinate system defined in its *.prj file: PROJCS["NAD_1983_StatePlane_Texas_South_Central_FIPS_4204_Feet",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],...


5

In ArcMap 10, use the MAKE QUERY TABLE. First, make sure the table and feature class are in the same geodatabase. The MAKE QUERY TABLE tool will create a feature layer so you'll have to use COPY to create a permanent feature class. You can then export this feature class to a shapefile.


5

The short answer is you can't produce a clipped sid image in ArcGIS 10, the format is proprietary. But you should be able to produce a clipped jpeg 2000 (jp2) which is similar. A typical sid compression is 20 to 1. In ArcGIS to get the equivalent compression on a jp2 you would use 5 in the Environment - Raster Storage form compression quality box. (That ...


5

Sounds like a permissions issue. Try working through these steps: http://blogs.esri.com/esri/arcgis/2008/09/05/digging-deeper-troubleshooting-geoprocessing-errors-when-using-arcsde


5

You can still do this, although Excel has made it more difficult by breaking the loading sequence for associated .dbf's. Export normally. Open Excel on its own. Go to File -> Open, select 'All Files', and find the .dbf.


Only top voted, non community-wiki answers of a minimum length are eligible