Tag Info

New answers tagged

0

The order of parameters is incorrect. Have a look at the Project tool help page, your parameter order must be the same as this. Also what is "GCD_WGC_1984"? To my knowledge that's not a coordinate system, I'm guessing you meant GCS_WGS_1984? But that's academic as your parameters are out of order!


0

I've seen this issue in one of my applications that has considerable Geoprocessing requirements. In my case I saw very similar symptoms running it with 4Gb RAM on several machines. I then tested on several machines with 12Gb RAM and it ran cleanly on them. My advice is, try to test it on a machine with more RAM before spending more effort trying to debug. ...


0

I recently dealt with a similar problem; I wanted to pass an SDE layer as .near_features and it only worked when I passed the layer's path string. Talked to the Esri Support and they came back a few days later and said it was a bug. Bug number: #NIM102924 Synopsis: Generate Near table geoprocessing tool s "Near Features" parameter doesn't accept ...


0

The answer seems to be that the Overlay tool in the Linear Referencing toolbox does not like GUID data type fields as the Route ID field. The solution was to transfer the Route ID field which is a guid to a string field. All is working with the Route ID field as a string data type.


0

This should go as a comment. But is too long. Thanks much Gene! There are both LineString and MultiLineString features. So this code for line in shape(elem['geometry']) shows TypeError: 'LineString' object is not iterable when feature is MultiLineString. So the final code that works is as below. with fiona.open(source_path) as source: with ...


0

Is the DEM in geographic coordinates? If there is projection-on-the-fly going on based on geoprocessing settings, those settings are not seen if running the python process standalone. Note, you only need to set overwriteOutput once and it remains in effect, you can move it out the loop (though it's not a performance thing, just style!)


2

A MultiLineString is a list of lines: from shapely.geometry import MultiLineString, mapping, shape coords = [((0, 0), (1, 1)), ((-1, 0), (1, 0))] lines = MultiLineString(coords) print lines MULTILINESTRING ((0 0, 1 1), (-1 0, 1 0)) for line in lines: print line LINESTRING (0 0, 1 1) LINESTRING (-1 0, 1 0) # convert to GeoJSON format: print ...


0

There is an ArcgisToolbox built for Arcgis 10.1 freely available here Seems it is developed for the multispectral raster classification or multi band raster classification.. I've just tried for some single band raster, just to check and actually got an error related to file mosaic'ing.. and renameing file (248 line). If you just change the name if should ...


1

One challenge faced by address verification services is that some addresses are valid, yet not deliverable by USPS. I live in a small town in which USPS provides no home delivery; all residents must use PO boxes. However, both UPS and FedEx do provide delivery to our homes, so clearly our addresses mean something. ;-) Anyhow, the SmartyStreets ...


1

There's another option that isn't as elegant as Erase or Union and involves a few more steps, but it can be done at Basic and without third party tools (which can be an issue if you don't have administrative rights on the machine). Make a copy of your buffer layer, bufferdummy. If you have more than one layer, you can select them all and export/copy them ...


1

If you don't have an info level of license you can do it using union, which is limited to two feature classes in basic and standard, but that's enough. Union your two layers, to keep it simple join the attributes with FID only, and then delete the ones that have a FID value from the second feature class. There are a few ways to delete but if you wish to do ...


2

The ArcGIS Erase tool will work, but as user two seven three nine says, you have to have the required license level to use it. The Xtools Pro extension suite has an Erase tool built in that will do the same thing and is free. Some of the features of Xtools require that you buy a license after the 14-day evaluation period is over, but the Erase tool (and a ...


4

Erase is the tool for this in ArcGIS. http://resources.arcgis.com/en/help/main/10.1/index.html#//00080000000m000000 Kind of like an "opposite clip" You need a higher level license to do this.


0

Your getting the non-descriptive error "IDW failed" because of your error handler. Take the code out of the error handler, re-run, then report the traceback.It would be far more useful to everyone then what your error handler is reporting. I don't have the cred to add this as a comment but your code needs to get cleaned up too. There is some redundancy ...


1

I'm only addressing the "How does one take a text string, and pass it off as a text expression?" part here but try using Python string formatting as below: RC1 = str(BotVal) RC2 = str(MegaCount - 1) NODATA = "NODATA" RC3 = str(MegaCount) RC4 = str(TopVal) Val = "1" Qu = "" Equ = "{0} {1} {2}; {3} {4} {5} {6}".format(RC1,RC2,NODATA,RC3,RC4,Val,Qu) print Equ ...


0

I think you'll find your problem is: arcpy.gp.Reclassify_sa should be: arcpy.sa.Reclassify it goes like this: Raster = arcpy.sa.Reclassify(InRaster,"Value",Equ, "NODATA") Raster.save(SingleValueNamePath) The spatial analyst / raster objects have changed since 9.x. BTW, there's nothing wrong with your iteration but the pythonic way to do it is: ...


0

The progress bar only appears to be jumping up and down. What is actually happening is that the geoprocessing tool has an individual percentage/progress bar for each of its subset calculations. A subset calculation will complete at 100% and then a new calculation will begin at 0%. ArcGIS does not always inform you that this is happening, and when it ...


1

i solved the problem by converting the tow TIN to a raster before execute the cut fill 3D analyst tool :


2

This has nothing to do with your Python version or your os. The Python Networkx module (nx_spatial is very old, use directly NetworkX: read_shp) do not cut the lines at each intersection and simplify the line geometries into start and end coordinates. Generates a networkx.DiGraph from shapefiles. Point geometries are translated into nodes, lines into ...


2

There is a way to do this directly in ArcGIS and it's fairly simple. It's the same answer I posted to this question, with a slight change in process output since you have overlapping buffers. Give your buffers a unique name attribute (rather than just using oid, which you could do, but having a (short) name to follow will make it easier). Now you need your ...


0

Alright this is a great use of SQL Spatial. This answer is using MS SQL Spatial, but could just as easily be done with PostGIS and I suppose Oracle Spatial (never used it...) For this example, I used school points from DPS (sorry can't find them online right now...) and the street centerlines from City and County of Denver loaded into a SQL Server database ...


0

If you don't have an Advanced license: Union the old and the updated polygons. Select the features having FIDs from both input feature classes and delete them. Merge the updated feature class with the output of the previous step.


2

For mean values of elevation in building polygons... In QGIS: Activate Zonal Statistics plugin Run it - "Raster/Zonal Statistics", it will add new columns to polygon layer the one needed is "mean". Duplicate your elevation layer (raster/grid) - you will use it for final file Use "Raster / Conversion / Rasterize" select you polygon (building) layer, ...


1

You are right that dumping the multipolygons to polygons makes the index more efficient, for the reason you mention. So if the size of the dataset makes performance an issue you should: 1) dump the multipolygons into a new table. 2) create a spatial index on the polygons 3) get the intersecting part and restrict the query to polygons theat intersect. ...


1

This answer works only if you want constant value - in example 0. In QGIS: Clip elevation by polygons (buildings) to new raster (grid) - buildings_grd Raster/Extraction/Clipper (Use Clipping mode "Mask Layer" and your building shapefile) Run raster calculator and calculate Raster/Raster Calculator... Expression: ("elevation" = "buildings_grd")*0 + ...



Top 50 recent answers are included