Tag Info

New answers tagged

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.


0

My suspicion is that Clip is bailing when the shapefile it tries to write hits the 2Gb size limit. It would be good to confirm, by using Copy As Python Snippet, that you are using precisely the same inputs and outputs when this works in ArcMap and appears not to from ArcPy.


1

I am probably late but it may be helpful to others. I was in need of the same functionality in 10.2 that we had in ArcView 3.2, in winch we can cut, union, merge, clip, intersect, and other overlaying operations directly inside the editing session. I have just found my way by following this procedure (providing you have ArcEditor, as I do and do not know if ...


4

Reorganise your shapefile so that one shapefile contains one feature (A,B,C in your case) only Then use a loop like for i in A B C; do gdalwarp -cutline $i.shp ... $i.tif done to create each output raster. Example of script: #!/bin/sh # "shp" - folder for shapefiles # "outputraster" - folder for output rasters cd /home/user/cliprasters/ # ...


3

The USGS DEM format is a default GDAL input and output format. See this link for more information. This means you can "clip" is using the -a_ullr command in gdal_translate or use gdal_warp using the -cl and -crop commands. Here is some code to get you moving gdalwarp -of DEM -cutline C:\temp\area_of_interest.shp -crop_to_cutline C:\temp\input.dem ...


3

It is an open format See Link but not one typically used as an Output. I believe GRID is a far more commonly used format for DEM data but that may just be a personal preference. I did find a reference for converting "TO" that format but it is an old solution (written in Avenue ... the old language used with ArcView 3.X). A link I found indicated that it may ...


0

A workaround would be to re-add the attributes after doing the clip, while converting from SpatialPolygons to SpatialPolygonsDataFrame. clip = gDifference(spdf1,spdf2, byid=TRUE) row.names(clip) = c("p1","p2") spdf3 <-SpatialPolygonsDataFrame(clip,data=as.data.frame(spdf1@data)) spdf3@data variable1 variable2 p1 232 235 p2 242 ...


2

Using QGIS, Vector -> Reserach Tools -> Polygon from layer extent gives you a polygon that you can use for clipping the other vector data to. There is a GDAL module called gdaltindex that does the same without using other software: gdaltindex N51E007.shp N51E007.hgt creates a bounding polygon for a sample SRTM data file. For a python solution to ...



Top 50 recent answers are included