Tag Info

New answers tagged

1

The single tool most suitable to your goal is Identity, which will cut up your A layer by your B layer and not leave any polygons based on the B layer that don't overlap A. However it requires an ArcInfo license. The next suggestion would be Union, which will return all polygons cut up by overlap. This will return not only the areas of common overlap ...


1

This can be done with the identity tool in one step, but identity needs an ArcInfo licence. Without this licence, you can use the union tool to create a new feature class. Then you delete the polygons where the old FID field corresponding with the polygon to split is equal to -1 (those are the polygons that only belong to the polygon that you use to split)


1

Can you do 2 operations, one a clip, the other an intersect and then append the 2 result sets together?


1

Using the Cut Polygons tool within an Edit Session allows you to split polygons by drawing a line, intersections in your case.


5

You can place all the feature classes within a single dataset then use model builder to clip all at once and output to a new dataset:


5

In you case (only one tool), you can simply right click on the clip tool in your toolbox, and launch it as a "batch". See the ArcGIS help for more details


6

You can use a simple block of code either directly from Python window in ArcGIS or as a Python script or as a custom script tool. Below is a sample code for a script tool. Just add your vector layers into a map document (.mxd) and specify the clip layer and the output geodatabase. import arcpy arcpy.env.overwriteOutput = True clip_layer = ...


1

1.Which line returns the error? Maybe you should check if the layer can be converted to a feature layer and clipped using if layer.isFeatureLayer: 2.Besides this, your script will overwrite clipped_fc at every iteration. Change the name to something like clipped_fc = env.scratchFolder + os.sep + layer.name + "_clipped.shp". Check that the layers' name is ...


0

You can use gdal_calc.py see: http://www.gdal.org/gdal_calc.html if your mask is a raster of same dimension gdal_calc.py -A input.tif -B mask.tif --A_band=3 --B_band=3 --outfile=result.tif --calc="A*B" That will multiply each cell of the raster with the corresponding cell of the mask. Where the mask cell is 0,0,0 the output cell will be 0,0,0 and where ...


3

If you have the purple polygon selected when you perform the clip it will only output the selected polygon.Deselect all the features you plan on clipping and try the clip again. That should clip everything to the larger polygon.


0

QGIS does not appear to have direct functionality to clip a vector layer by the extent of a raster layer. You can only clip by another vector layer. However, you can use the plug-in Image Boundary to create a polygon that matches the boundary of the raster layer. With that in hand, you can use Vector->GeoProcessing Tools->Clip to clip the vector you're ...


0

If you don't mind doing this manually: Start editing in ArcMap. Double click on the feature that needs holes to start editing the vertices. Find the continue feature tool that pops up (read http://resources.arcgis.com/en/help/main/10.1/index.html#//01m800000039000000). Right click and select finish part to tell ArcMap you don't want to modify the ...


3

Two ways spring to mind. The most efficient depends on how many inner polygons you have. Clip: Select one of your existing 'lake' polygons then choose clip from the editor menu. Use a buffer distance of 0 and discard the area that intersects. This will clip a hole in any overlapping polygons. You will have to do this for each inner polygon. Cut Polygon ...


1

Thats not a bug. You have a point layer, point dont have an area (only 1 dimension) so you can't clip it. For doing what you want, you have to intersect them or make a spatial query with a contain.


6

Seems to be a simple application of gDifference from the rgeos package: > require(rgeos) > ukhole = gDifference(uk, lnd) Warning message: In RGEOSBinTopoFunc(spgeom1, spgeom2, byid, id, "rgeos_difference") : spgeom1 and spgeom2 have different proj4 strings > plot(ukhole) The projection warning is because the LondonBoroughs shapefile doesn't ...


3

I've created a small function for this very purpose and it has been used by others with good reviews! gClip <- function(shp, bb){ if(class(bb) == "matrix") b_poly <- as(extent(as.vector(t(bb))), "SpatialPolygons") else b_poly <- as(extent(bb), "SpatialPolygons") gIntersection(shp, b_poly, byid = T) } This should solve your problem. Further ...



Top 50 recent answers are included