Tag Info

New answers tagged

0

If the Reference System for the shape is equal, you can use "intersect" tool (vector>geoprocessing tools>intersect). The vector result have, in the attribute table, the values of two vector that you have processed (wetland areas and boundary).


0

This question is a bit old - but thought I'd mention that there is a free script for erasing with an ArcView (standard) licence (from back in the 9 releases). Check it out for future reference: http://arcscripts.esri.com/details.asp?dbid=14625


0

Use the ST_Intersects() function then save your query as a new table by going to Options > Table From Query. SQL code would be something like this: SELECT a.* FROM polygon_table AS a, point_table AS b WHERE ST_Intersects(a.the_geom, b.the_geom)


2

I'm getting the same error in my Arc 10.2 installation. I've encountered issues in the past when using geometry objects with geoprocessing tools, as you're doing with the geom object from your cursor and the Clip tool. It's just ironic since this is almost exactly the example in the Esri help page. One solution is to avoid using geometry objects by exporting ...


3

I can't explain why Clip is failing when you pass in the geometry. It happens for me as well. I can offer this workaround, though. Persist that geometry out to in_memory (or to disk), then clip with it. with arcpy.da.SearchCursor(splitFC, [splitField,"SHAPE@"]) as cursor: for row in cursor: name = row[0] geom = row[1] ...


1

Try: polygone1 <- gBuffer(polygone1, byid=TRUE, width=0) polygone2 <- gBuffer(polygone2, byid=TRUE, width=0) clip2 <- gIntersection(Polygone1, Polygone2, byid=TRUE) It is ugly, but it usually solves this kind of problem. HTH


2

Rather than using the GNU Parallel method you could use the python mutliprocess module to create a pool of tasks and execute them. I don't have access to a QGIS setup to test it on but multiprocess was added in Python 2.6 so provided that you are using 2.6 or later it should be available. There are a lot of examples online on using this module.


0

I can't say I have a solution to the problem using subprocess as I had similar problems to what you describe before. In the end I decided to use os.system instead of subprocess (I am aware that subprocess is meant to replace os.system). So what I do is: cmd = "gdalwarp ..." %(...) os.system(cmd)


1

Here is the gnu parallel solution. With some care most emabrrassingly parallel linux based ogr or saga algorithms could be made to run with it inside your QGIS installation. Obviously this solution requires the installation of gnu parallel. To install gnu parallel in Ubuntu, for example, go to your terminal and type sudo apt-get -y install parallel NB: I ...


1

I am using 10.2 but had no issues clipping your TIN. Your polygon is not a 3D feature so you may have incorrectly set the Height Field parameter which by default is, in this case, set to iD - this will force to clip the TIN at that height which is 0. Set the Height Field to NONE and run it again: The resulting TIN will be clipped as follows:


8

If you change your program to read the file name from the command line and split up your input file in smaller chunks, you can do something like this using GNU Parallel: parallel my_processing.py {} /path/to/polygon_file.shp ::: input_files*.shp This will run 1 job per core. All new computers have multiple cores, but most programs are serial in nature ...


1

I figured it out, FYI, or rather I found the answer on ESRI help page: http://resources.arcgis.com/en/help/main/10.1/index.html#//01m60000002w000000 Clip the inner buffer in an editing session rather than with the Clip tool.



Top 50 recent answers are included